geoshow cuts off grid data

2014 2 14
1 个回答
回答已采纳
更新时间：2018 10 18
14 次查看（30 天）

0 个投票

Hello,
I have a global grid of data (nLat,nLon = 24,36) which I'm trying to map with the geoshow() function. However, it's cutting off several columns on the left and right side of the map. I've tried different ways of mapping as described in the help, but the same thing keeps happening.
Below a minimal example. What am I doing wrong? Should I use a different function?
Thanks & cheers,
lon = -5:10:355;
lat = -90:7.5:90;
[latM lonM] = meshgrat(lat,lon);
data = rand([24,36]);
figure
ax = axesm('MapProjection', 'eqdcylin','grid','on');
mlabel;plabel;
geoshow(latM,lonM,data,'DisplayType','texturemap')

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

 采纳的回答

Kelly Kearney
Kelly Kearney 2014-2-14

0 个投票

Unfortunately, all the mapping toolbox plotting functions clip data to the map frame. If part of a line segment, patch object, etc. overhangs the frame, then that graphics object is not plotted at all. Sometimes not noticeable, and other times, as in your case, very noticeable.
The only way around it (other than preprocessing your data, via maptrimp/l/s and upsampling of raster data) that I know of is to fiddle with the exact limits of the frame to match it up to the boundaries of your data.
It might help to rearrange your input data to the same order as it will be displayed. For example, you get less trimming if you sent your map limits to [0 360] (like your data) than [-180 180]:
lon = -5:10:355;
lat = -90:7.5:90;
[latM lonM] = meshgrat(lat,lon);
data = rand([24,36]);
subplot(2,1,1);
ax = axesm('MapProjection', 'eqdcylin','grid','on', ...
'frame', 'on');
mlabel;plabel;
geoshow(latM,lonM,data,'DisplayType','texturemap');
tightmap;
subplot(2,1,2);
ax = axesm('MapProjection', 'eqdcylin','grid','on', ...
'frame', 'on', 'maplonlimit', [0 360]);
mlabel;plabel;
geoshow(latM,lonM,data,'DisplayType','texturemap');
tightmap

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

Maarten
Maarten 2014-2-16
Thanks, Kelly, that was helpful. I solved the problem by shifting the map limits by 1/2 of a longitude step so that they coincide with the edges of the of columns.
It turns out I have the same problem for the latitudes. Here this solution does not work because you cannot specificy map limits outside the [-90 90] range. So I simply omitted the top and bottom rows (which do not contain interesting information anyway) and shrunk the map with 1/2 a latitude step at the top and bottom.

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Create Plots on Maps 的更多信息

产品

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by