plotting rainfall data, having latitude and longitude, on map

Question

Madan Kumar 2017-5-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/342579-plotting-rainfall-data-having-latitude-and-longitude-on-map

评论： Arnob Das 2021-11-27

Hi, I have problem with plotting rainfall data (with latitude and longitude) on map. I googled it, there were many answers but couldn't do it. My data has latitude(1x86400),longitude(1x86400) and corresponding rainfall (1x86400). I have to plot rainfall data on map corresponding to latitude and longitude. My code is similar to below one but not working though not showing any error. Need help.Thank you.

A=importdata('new.dat');
lat=A(:,1);
lon=A(:,2);
lat1=meshgrid(lat,1:25);
lon1=meshgrid(lon,1:25);
z=A(:,3);
z1=meshgrid(z,1:25);
worldmap('world')      % initializes map
contourm(lat1,lon1,z1)    % plots contours
c = load('coast.mat'); % loads coastlines
plotm(c.lat,c.long)    % plots coastlines

2 个评论
显示无隐藏无

the cyclist 2017-5-30

Can you be more specific about what is not working?

Is it possible for you to upload your file new.dat, or MAT file with the data, so that people can run your code? (I won't be able to, since I do not have the Mapping Toolbox.)

Madan Kumar 2017-5-31

new.txt

I have attached new.txt file since new.dat file was not supported. When I am running the program there is no error but the data are not plotted on map.

请先登录，再进行评论。

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

Answer 1

KSSV 2017-5-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/342579-plotting-rainfall-data-having-latitude-and-longitude-on-map#answer_269000

在 MATLAB Online 中打开

Your data is a scattered data.....you cannot use meshgrid as the way you have used. You need to interpolate the data into a grid and then plot. Check the below code:

A=importdata('new.dat');
lat=A(:,1);
lon=A(:,2);
z=A(:,3);
lon0 = min(lon) ; lon1 = max(lon) ;
lat0 = min(lat) ; lat1 = max(lat) ;
N = 100 ;
x = linspace(lon0,lon1,N) ;
y = linspace(lon1,lat1,N) ;
[X,Y] = meshgrid(x,y) ;
F = scatteredInterpolant(lon,lat,z) ;
Z = F(X,Y) ;
worldmap('world')      % initializes map
contourm(X,Y,Z)    % plots contours
c = load('coast.mat'); % loads coastlines
plotm(c.lat,c.long)    % plots coastlines

2 个评论
显示无隐藏无

Madan Kumar 2017-5-31

Thank you so much. It plots world map and data over it. But I need to plot only the region required for me and data over it. I tried using worldmap([-50 50],[160 -30]). Although it is showing the required map but not plotting data over it. Need help.

Arnob Das 2021-11-27