Question

你好
我想以2D格式表示具有2个变量（纬度和经度）的数据。该值由颜色表示，将2个变量表示为2轴，我正在使用Contourf函数绘制数据。所有数据均来自xlsx文件，我将其放在矩阵中。

Locations = xlsread('Availability results.xlsx');
column_numberloc = 1; % Column in the locations file containing the number of the locations
column_latitude = 2; % Column in the locations file containing the latitude of the locations
column_longitude = 3; % Column in the locations file containing the longitude of the locations
column_availability = 4; % Column in the locations file containing the availability of the locations

min_latitude = min(Locations(:,column_latitude));
max_latitude = max(Locations(:,column_latitude));
min_longitude = min(Locations(:,column_longitude));
max_longitude = max(Locations(:,column_longitude));
max_availability = max(Locations(:,column_availability));
min_availability = min(Locations(:,column_availability));
longitude = Locations(:,column_longitude);
latitude = Locations(:,column_latitude);
Contour = zeros(23,17);

for numerofile=1:204
    [coord_x,coord_y] =transformation(Locations(numerofile,column_latitude),Locations(numerofile,column_longitude));
    Contour(coord_x,coord_y) = Locations(numerofile,column_availability);
end

for i=1:23
    for j=1:17
        if Contour(i,j) == 0
            Contour(i,j) = NaN;
        end
    end
end

cMap=jet(256);
figure(1);

x = linspace(min_longitude,max_longitude,17);
y = linspace(min_latitude,max_latitude,23);
newpoints = 100;

[xq,yq] = meshgrid(linspace(min(x),max(x),newpoints),linspace(min(y),max(y),newpoints ));
Contourq = interp2(x,y,Contour,xq,yq,'linear',max_availability);


[c,h]=contourf(xq,yq,Contourq,100);

%[c,h]=contourf(x,y,Contour,50);
set(h, 'edgecolor','none');
colormap(cMap);
cb=colorbar;
caxis([min_availability max_availability]);

transformation 函数允许我将所有数据放置在Contour矩阵中，因为它将经度和纬度与行和列相关联。

我为每个等于零的数据设置了一个NaN以便更好地查看我的数据，我得到了： interpolation_linear

哪个很好，但我希望此数据接近： Without interpolation

因此，我决定将线性插值更改为“最近”插值，我得到了： interpolation_nearest

我可以看到更多数据，但是轮廓图并不像线性插值法那么平滑。

我看过很多关于如何制作平滑轮廓图的帖子（这就是我发现函数“ interp2”的方式），但是我认为我的问题来自NaN数据，这使我无法在边缘获得平滑轮廓图在NaN值和其余值之间（如第一张图片），但具有足够的数据（如第三张图片）。

我的问题是：您是否知道如何获得平滑的边缘轮廓图，这要归功于最近的插值法，但具有与第一个图像一样好的视觉效果？

非常感谢您

Answer 1

由于要在正方形网格上进行插值，因此可以直接使用imagesc显示2D图像。好处是您可以访问图像对象的AlphaData属性，该属性可用作显示蒙版。

r=rand(50);            % random 50x50 array
r(11:20,11:20)=NaN;    % some hole filled with NaN
imagesc(r)             % show the image, with NaN considered as the lowest value in color scale
imagesc(r,'AlphaData',~isnan(r)) % show the image, with NaN values set as fully transparent

您还可以：

首先设置显示蒙版
用一些有意义的值（最近的非NaN值？）替换零或NaN
使用interp2进行插值，甚至可以使用'cubic'参数进行插值，以提高平滑度
由于AlphaData中设置了显示蒙版，因此仅显示图像中有意义的部分。

平滑边缘轮廓图

1 个答案: