0

我想使用 matlab 创建图,以径向方式表示质量的数值评估。

我发现的最佳方法似乎无法正常工作。一个运行以下代码:

theta = (0 : (360/11) : 360)*pi/180;
r = 0 : 2 : 20 ;
[TH,R] = meshgrid(theta,r);
[X,Y] = pol2cart(TH,R);
Z = meshgrid(Data);
surf(X,Y,Z);

Data是包含 11 个数字的数据向量,示例数据集如下:

Data = 0.884, 0.882, 0.879, 0.880, 0.8776, 0.871, 0.8587, 0.829, 0.811, 0.803, 0.780 

这里的输出surf是这样的:

不正确的冲浪图

我想制作这种类型图像的更精致版本:

大部分正确的输出

我使用以下代码生成的:

for theta = 0 : pi/100 : pi;  
    v = [InterpolatedImageHeight;LengthVector];
    x_center = InterpolatedImageHeight((HorizontalRes+1)/2);
    y_center = 0; %InterpolatedImageHeight((HorizontalRes+1)/2);
    center = repmat([x_center; y_center], 1, length(InterpolatedImageHeight));
    R = [cos(theta) -sin(theta); sin(theta) cos(theta)];
    vo = R*(v - center) + center;
    x_rotated = vo(1,:);
    y_rotated = vo(2,:);
    scatter(x_rotated,y_rotated,DotSize,InterpolatedData,'filled'); %x,y,area,color,properties
end

这个问题是它是一个散点图,我基本上使用它plot(r,Data),绘制许多副本,并增加点大小。图形本身有很多接缝,这会占用大量内存,并且在运行速度极快且占用内存最少的地方surf非常mesh耗时。

如何产生具有可变颜色输入的同心环?

4

2 回答 2

1

您的问题中有两个完全不同的情节。第一个将数据表示为从原点到圆外的射线。数据点逆时针放置。可以这样实现它的改进版本:

Data = [0.884, 0.882, 0.879, 0.880, 0.8776, 0.871,...
        0.8587, 0.829, 0.811, 0.803, 0.780]; 

theta = linspace(0,2*pi,length(Data));
r = linspace(0,20,length(Data));
[TH,R] = meshgrid(theta,r);

Z = meshgrid(Data);
[X,Y,Z] = pol2cart(TH,R,Z);

surf(X,Y,Z);
view(2);
shading interp

请注意,我曾经linspace生成thetar始终匹配Data. Z也是通过低谷pol2cart。然后您可以使用shading interp删除补丁之间的线条并插入颜色。view(2)您可以设置透视图,就像拥有 2d 绘图一样。

这是结果:

结果1


获得第二个示例中的结果相对容易。那里的数据点代表围绕原点的同心圆,并从原点向外放置。因此,只需Z使用以下行转置网格网格:

Z = meshgrid(Data)';

这就是结果:

结果2

于 2015-07-10T20:27:08.777 回答
0

基于此线程中 Darren Rowland 的代码,我提出了以下解决方案:

x = interp1(1:length(data),datax,(datax(1):datax(end)/f:datax(end)),'linear'); 
y = interp1(1:length(datay),datay,datay(1):datay(end)/f:datay(end),'spline');
theta = linspace(0,2*pi,n);

xr = x.'*cos(theta);
zr = x.'*sin(theta);
yr = repmat(y.',1,n);

figure;
surf(xy,yr,zr,zr*numcolors);

它优雅,运行迅速,并产生美丽的人物。这是带有一些额外图表元素的输出示例:

代码生成的样图:Canon EF 50mm f/1.4 Average MTF 30lp/mm map

于 2015-07-13T16:48:20.637 回答