使用ginput在MATLAB中绘制数字化

Question

我正在尝试使用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');

Answer 1

您可以使用

以双对数刻度绘图

loglog(x,y);

help loglog或文档提供了其他信息。

对于单个对数刻度使用

semilogx(x,y);
semilogy(x,y);

Answer 2

以下是一个概念证明，可以帮助您走上正轨。我用原始代码替换了我认为是“良好做法”的东西。

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文件中获取此图像，它将显示为矢量图像（您可以通过尽可能多地放大来测试，而不会使图表变为像素化;此似乎< / em>从.png看起来的方式来看你的情况，你最好将它保存为矢量图像，然后你有两个选择：
     • 将图像导出到位图，大大提高分辨率，然后再次尝试数字化过程。
     • 将矢量图像保存为.svg，然后使用您喜欢的文本编辑器打开文件并获取点的确切坐标。
  2. 如果源图像是位图（与矢量图形相对），则可以"trace the bitmap"，然后将其转换为矢量，然后#GOTO步骤1.
• 此解决方案（目前）不支持调整图形的大小。
• Position设置中显示的幻数是下图中解释的缩放因子（size(im)也是[609 785 3]）。从技术上讲，这些可以使用“原始图像处理”找到，但在这种情况下，我只是明确地对它们进行了硬编码。