how to show my output data in map.

5 次查看（过去 30 天）

devendra 2014-6-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133481-how-to-show-my-output-data-in-map

评论： Srishti Gaur 2018-5-3

sir i want to plot my own output data in map.please help me. how will do this.sir i have 8 years annual time series of rainfall frequency and sir my data size is (20x22x8)(lon,lat,time).

thank you in advance.please sir help me.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

José-Luis 2014-6-10

Duplicate question

How to ask a question

请先登录，再进行评论。

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

采纳的回答

Kelly Kearney 2014-6-12

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133481-how-to-show-my-output-data-in-map#answer_140609

编辑：Kelly Kearney 2014-6-12

在 MATLAB Online 中打开

I reaally shouldn't be giving you the answer, since you've shown absolutely no effort to create this plot yourself. But I'm bored at work waiting for some simulations to finish, so it's your lucky day.

First of all, I'm assuming that you already have a shapefile holding the India boundary, and that you have a raster grid of x (lon), y (lat), and z coordinates. Also, I'm assuming that your data grid isn't already trimmed to the India borders.

So, the first quick-and-dirty way to do this is to simply mask the data points outside the border:

file = 'TM_WORLD_BORDERS-0.3.shp';
S = shaperead(file, 'usegeocoords', true, 'selector', ...
    {@(x) strcmp(x,'India'), 'NAME'});
x = linspace(min(S.Lon), max(S.Lon), 100);
y = linspace(min(S.Lat), max(S.Lat), 100);
[x,y] = meshgrid(x,y);
z = rand(100);
isin = inpolygon(x,y,S.Lon,S.Lat);
z2 = z;
z2(~isin) = NaN;
figure('color','w');
worldmap('India');
pcolorm(y,x,z2);
plotm(S.Lat, S.Lon, 'k');

You'll notice that things get a little jagged around the edges, though. You can get a cleaner map by plotting all the data, then overlaying it with a polygon that masks out all but your region of interest. (Or in this case, polygons, since neither geoshow nor patchm tend to play well with polygons with holes):

lnlim = [65 100];
ltlim = [5 40];
lt = linspace(ltlim(1), ltlim(2), 3);
ln = linspace(lnlim(1), lnlim(2), 3);
for ii = 1:2
    ltbox{ii} = lt([1 2 2 1 1]'+(ii-1));
    lnbox{ii} = ln([1 1 2 2 1]'+(ii-1));
end
[lnmask, ltmask] = deal(cell(2));
for ii = 1:2
    for jj = 1:2
        [lnmask{ii,jj}, ltmask{ii,jj}] = polybool('-', ...
        lnbox{ii}, ltbox{jj}, S.Lon, S.Lat);
    end
end
ltboxall = ltlim([1 2 2 1 1]);
lnboxall = lnlim([1 1 2 2 1]);
[lnmaskall, ltmaskall] = polybool('-', lnboxall, ltboxall, S.Lon, S.Lat);
figure('color','w');
worldmap('India');
pcolorm(y,x,z);
for ii = 1:4
    patchm(ltmask{ii},  lnmask{ii}, 'w', 'edgecolor', 'none');
end
plotm(ltmaskall, lnmaskall, 'k');

9 个评论
显示 7更早的评论隐藏 7更早的评论

devendra 2014-6-14

Kelly Kearney ma'am thank you so much for your support and you are very nice person. ma'am now i am stopping my conversation and i will try this my self but this is very difficult for me. and sorry for calling you sir.i did not have any idea that you are female.again sorry ma'am

Srishti Gaur 2018-5-3

Dear Kelly Kearney, How to make multiple sub-plots for the same? I tried but it not able to do. Thank you

请先登录，再进行评论。

更多回答（3 个）

David Sanchez 2014-6-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133481-how-to-show-my-output-data-in-map#answer_140174

在 MATLAB Online 中打开

you could take a look at the demos within Matlab documentation. You will find many working examples on data representation, like this one:

z=peaks(25);
surf(z);
colormap(jet);

Adapt your data to any of those examples.

5 个评论
显示 3更早的评论隐藏 3更早的评论

José-Luis 2014-6-10

There is no way you can get such a high-resolution map with the data you have.

devendra 2014-6-10

编辑：devendra 2014-6-11

sir there is not necessary that high-resolution map but problem is that how to do this thing.

请先登录，再进行评论。

David Sanchez 2014-6-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133481-how-to-show-my-output-data-in-map#answer_140407

在 MATLAB Online 中打开

You will have to present your data by year:

XX YY]= meshgrid(X1,Y1);
z=ltm(:,:,1);
surf(XX,YY,z')
xlabel('XX')
ylabel('YY')
view([0 90])

11 个评论
显示 9更早的评论隐藏 9更早的评论

José-Luis 2014-6-12

编辑：José-Luis 2014-6-12

So please post some code showing what you have tried so far.

devendra 2014-6-12

geoshow('map.shp', 'FaceColor', [0.9 0.9 0.9])

this code is polting india map but i have no idea about next step. i also read about contourfm function but how will use this for my output data i have no idea.

请先登录，再进行评论。

peterhack 2016-11-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/133481-how-to-show-my-output-data-in-map#answer_242925

在 MATLAB Online 中打开

As I have a question on the above example I will utilize this thread, hope that is alright. I am using the presented code to mask the data points outside the border. Works fine so far. However, I want to use the DisplayType surface. Thus I am using geoshow to plot the data. Unfortunately it does not fill the whole map within the borders as a buffer is left. Any hint?

file = 'TM_WORLD_BORDERS-0.3.shp';
S = shaperead(file, 'usegeocoords', true, 'selector', ...
    {@(x) strcmp(x,'India'), 'NAME'});
x = linspace(min(S.Lon), max(S.Lon), 100);
y = linspace(min(S.Lat), max(S.Lat), 100);
[x,y] = meshgrid(x,y);
z = rand(100);
isin = inpolygon(x,y,S.Lon,S.Lat);
z2 = z;
z2(~isin) = NaN;
lnlim = [65 100];
ltlim = [5 40];
lt = linspace(ltlim(1), ltlim(2), 3);
ln = linspace(lnlim(1), lnlim(2), 3);
for ii = 1:2
    ltbox{ii} = lt([1 2 2 1 1]'+(ii-1));
    lnbox{ii} = ln([1 1 2 2 1]'+(ii-1));
end
[lnmask, ltmask] = deal(cell(2));
for ii = 1:2
    for jj = 1:2
        [lnmask{ii,jj}, ltmask{ii,jj}] = polybool('-', ...
        lnbox{ii}, ltbox{jj}, S.Lon, S.Lat);
    end
end
ltboxall = ltlim([1 2 2 1 1]);
lnboxall = lnlim([1 1 2 2 1]);
[lnmaskall, ltmaskall] = polybool('-', lnboxall, ltboxall, S.Lon, S.Lat);
figure('color','w');
worldmap('India');
geoshow(y, x, z2, 'DisplayType','surface')
contourcmap('jet',min(z):1:max(z),'colorbar','on','location','horizontal')
for ii = 1:4
    patchm(ltmask{ii},  lnmask{ii}, 'w', 'edgecolor', 'none');
end
plotm(ltmaskall, lnmaskall, 'k');