在matlab循环中使用pcolor

Question

我有一些坐标为lat，lon和z的数据点（数百万个数据点）。每个坐标点代表网格的中心，大小为0.0005。    我想用颜色代表值Z来绘制它们。但是，我的数据不是网格数据，它们是散布数据点。我知道另一个函数＆＃34; geoshow＆＃34;也可以制作类似的地图，但是对于网格数据。使用matlab将散布数据转换为网格数据时内存不足。所以我用pcolor实现了这一点。以下代码：

clc
clear all
data = load('data.txt');% This data.txt file is large, about 100M and each line with lat, lon and value Z 
X = data(:,1); % Lat
Y = data(:,2); %Lon
Z = data(:,3); % Value Z
grid_dx=0.0005;
grid_dy=0.0005;
xl=X-grid_dx/2;
xr=X+grid_dx/2;
yu=Y+grid_dy/2;
yd=Y-grid_dy/2;
for i1=1:length(x)
xx=[xl(i1),xr(i1);xl(i1),xr(i1)];
yy=[yd(i1),yd(i1);yu(i1),yu(i1)];
zz=repmat(z(i1),2,2);
pcolor(xx,yy,zz)
colormap('jet')
end

对于小数据集，它运行良好。但是，它会花费太多时间来处理更大的数据集。     是否可以在没有网格化的情况下更快地绘制散点图数据？任何功能pcolor，geoshow或其他东西都可以。  谢谢你的帮助！

Answer 1

您可以使用scatter函数在不执行插值的情况下使用基于z的颜色绘制(x,y)点。

lat = sort(10*rand(100,1));  % sorting is optional, just to give a 
lon = sort(10*rand(100,1));  % gradient in the coloring
z = sort(rand(100,1));

scatter(lat, lon, [] , z, 'fill');

此代码生成