Plot grided latitude and longitude as pixels on map

Question

BN 2020-2-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/505780-plot-grided-latitude-and-longitude-as-pixels-on-map

编辑： Giuseppe Inghilterra 2020-2-16

I want to plot gridded latitude and longitude as pixels on the map, Latitude is 1 x 32 and longitude is 1 x 40 I don't know how to this I tried geoshow but it said latitude and longitude must have the same size.

Error using checklatlon (line 25)
Function GEOSHOW expected its first and second input arguments,
LAT and LON, to match in size.
Error in geovecshow (line 37)
checklatlon(lat, lon, 'GEOSHOW', 'LAT', 'LON', 1, 2);
Error in geoshow (line 278)
h = showFcn(varargin{:}); 

I attache latitude and longitude. latitude and longitude are taken from the NetCDF file for gridded precipitation data, but I just want to plot its pixels and I don't need precipitation data.

I want to plot them on this map:

Code to generate map:

% Produce the map
borders('Iran Islamic Republic of') % borders is a function from Matlab file exchange by Chad A Greene (https://www.mathworks.com/matlabcentral/fileexchange/50390-borders)
grid on
axis equal                       % I added this, important to maintain aspect ratio (@AD)
xlabel('Longitude')    
ylabel('Latitude')   
xtickangle(90)

I want to plot the pixels like this: (ignore the color bar and colors I just want pixels)

Thank you

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Giuseppe Inghilterra 2020-2-16

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/505780-plot-grided-latitude-and-longitude-as-pixels-on-map#answer_415814

在 MATLAB Online 中打开

I obtain this result:

by running the following code:

load lati.mat
load loni.mat
[LONG, LAT] = meshgrid(loni,lati);
Z = rand(length(lati),length(loni));
figure
surface(LONG, LAT, Z)
colorbar
xlabel('Longitude')
ylabel('Latitude')
xticks(loni)
yticks(lati)
xtickangle(90)
grid on
ax = gca;
ax.XAxis.FontSize = 7;
ax.YAxis.FontSize = 7;

See https://www.mathworks.com/help/matlab/ref/meshgrid.html for meshgrid function.

In my example I generate ramdomly Z matrix, instead you should have a matrix with data in function of latitude and longitude coordinates. From your plot I see that some "pixels" are blank. This means that you can fill Z matrix with NaN values, in this way you obtain easily blank pixels.

2 个评论
显示无隐藏无

BN 2020-2-16

在 MATLAB Online 中打开

Dear Giuseppe Inghilterra First of all I want to thank you.

So according to your answer, it is impossible to have just grids without any data. So this is a good I dead to define random data at least I can see pixels.

here is the figure using your method:

I achive it using this code:

borders('Iran Islamic Republic of') % from FEX
%grid on
axis equal                       % I added this, important to maintain aspect ratio
xlabel('Longitude')    
ylabel('Latitude')   
xtickangle(90)
hold on
%% plot gridded on map
[LONG, LAT] = meshgrid(loni,lati);
Z = rand(length(lati),length(loni));
figure
surface(LONG, LAT, Z)
colorbar
xlabel('Longitude')
ylabel('Latitude')
xticks(loni)
yticks(lati)
xtickangle(90)
grid on
ax = gca;
ax.XAxis.FontSize = 7;
ax.YAxis.FontSize = 7;

One main challenge is do you know how I can bring the boundary of the country on? or delete outside the boundary pixels?

Also, it would be awesome if it is possible to grids haven't any colors just be white. thank you I waiting for your reply.

Giuseppe Inghilterra 2020-2-16

编辑：Giuseppe Inghilterra 2020-2-16

在 MATLAB Online 中打开

See the following code:

I define a circumference in x,y coordinates;
then I evaluate if LONG,LAT coordinates are inside the defined circumference through "inpolygon" matlab function and I set to zero these entries of Z matrix;
then I take the current colormap cmap and redefine the first line with [1 1 1]. In this way I assign color white to "zero" value and then I update colormap;

load lati.mat
load loni.mat
t = 0:pi/10:2*pi;                   %(1)
x = 55 + 5*cos(t);
y = 30 + 5*sin(t);
[LONG, LAT] = meshgrid(loni,lati);
Z = rand(length(lati),length(loni));
Z(inpolygon(LONG,LAT,x,y)) = 0;     %(2)
figure
surface(LONG, LAT, Z)
cmap = colormap;                    %(3)
cmap(1,:) = [1 1 1];
colormap(cmap)
colorbar
xlabel('Longitude')
ylabel('Latitude')
xticks(loni)
yticks(lati)
xtickangle(90)
grid on
ax = gca;
ax.XAxis.FontSize = 7;
ax.YAxis.FontSize = 7;