0

我正在尝试使用 MATLAB 将此图像数字化：

对数比例图

我有以下脚本：

%// Get data from plot
clear all; close all;
%// Input
fname = 'Fig15a.PNG'; 
 xvec = [1e3:1:1e8]; 
 yvec = [1e-4:1:1e-1]; 
   xt = [1e3 1e4 1e5 1e6 1e7 1e8]; 
   yt = [1e-4 1e-3 1e-2 1e-1];

%// Read and plot the image
im = imread(fname);
figure(1), clf
im = im(end:-1:1,:,:);
image(xvec,yvec,im)
axis xy; 
grid on;

%// Set ticks
set(gca,'xtick',xt,'ytick',yt); %// Match tick marks

%// Collect data
[x,y] = ginput; %// Click on points, and then hit ENTER to finish

%// Plot collected data
hold on; plot(x,y,'r-o'); hold off;

%// Then save data as:
save Fig15a.mat x y

该脚本工作正常 线刻度输出

有没有办法可以将x和y轴更改为对数刻度？我尝试在不同的地方添加以下代码，但没有运气：

%// Set Log scale on x and y axes
set(gca,'XScale','log','YScale','log');
4

2 回答 2

0

下面是一个概念证明，可以让您走上正轨。我已将您原始代码中的内容替换为我认为的“良好做法”。

function q36470836
%% // Definitions:
FIG_NUM = 36470836;
%% // Inputs:
fname = 'http://i.stack.imgur.com/2as4t.png'; 
xt = logspace(3,8,6);
yt = logspace(-4,-1,4);
%% // Init
figure(FIG_NUM); clf
% Read and plot the image
im = imread(fname);
hIMG = imshow(im); axis image;

%// Set ticks
hDigitizer = axes('Color','none',...
                  'XLim',[xt(1) xt(end)],'YLim',[yt(1) yt(end)],...
                  'XScale','log','YScale','log',...
                  'Position',hIMG.Parent.Position .* [1 1 696/785 (609-64+1)/609]);

uistack(hDigitizer,'top'); %// May be required in some cases
grid on; hold on; grid minor;

%// Collect data:
[x,y] = ginput; %// Click on points, and then hit ENTER to finish

%// Plot collected data:
scatter(x,y,'o','MarkerEdgeColor','r');

%// Save data:
save Fig15a.mat x y

这是它的外观示例：

输出

几点注意事项：

  • xt,yt可以使用logspace.
  • 很难（可能不可能）将数字化网格与图像正确对齐，这将不可避免地导致数据错误。尽管这在以下情况下会有所帮助（您需要矢量图形编辑器，例如免费软件InkScape）：
    1. 如果您有任何机会从 PDF 文件中获得此图像，该图像显示为矢量图像（您可以通过尽可能多地放大来测试此图像，而不会使图表像素化；这似乎是您的情况外观），您最好将其.png保存为矢量图像，然后您有两个选择：
      • 将图像导出为分辨率大大提高的位图，然后再次尝试数字化过程。
      • 保存矢量图像，.svg然后使用您喜欢的文本编辑器打开文件并获取点的确切坐标。
    2. 如果源图像是位图（与矢量图相反），您可以“跟踪位图”，从而将其转换为矢量，然后 #GOTO 步骤 1。
  • 此解决方案（当前）不支持调整图形大小。
  • 设置中出现的幻数Position是下图中解释的比例因子（也是size(im)）[609 785 3]。从技术上讲，这些可以使用“原始图像处理”找到，但在这种情况下，我只是明确地对它们进行硬编码。 幻数的解释
于 2016-04-07T11:01:46.113 回答
0

您可以使用双对数刻度绘制

loglog(x,y);

help loglog或文档提供更多信息。

对于单个对数刻度使用

semilogx(x,y);
semilogy(x,y);
于 2016-04-07T09:36:31.950 回答