Changing the orientation Arctic Projection maps

5 次查看(过去 30 天)
[LON,LAT] = meshgrid(lon1,lat1);
% hold on
fs=20;
load coastlines
figure('Color','w');
axesm('stereo','MapLatLimit',[65 90],'FontSize',fs,'FontWeight','bold');
axis off
framem off
gridm on
mlabel on
plabel on;
setm(gca,'MLabelParallel',0);
setm(gca,'MLabelLocation',30);
setm(gca,'PLabelLocation',15);
geoshow(LAT,LON,D,'DisplayType','texturemap');
hold on;
geoshow(coastlat,coastlon,'Color','k','LineWidth',2);
hold on;
cmap = jet(256);
cmap(1,:)=[0.9 0.9 0.9];
colormap(cmap)
what modification should I do in my code to make the projection similar to the map I have attached? The position of Greenland should be same in both the plots

回答(1 个)

Chad Greene
Chad Greene 2022-11-8
I typically bypass the Matlab Mapping Toolbox for stereographic projections, because the functions are kind of clunky and slow.
Instead, I use my Arctic Mapping Tools, which puts Greenland in the center (central meridian 45 degrees west) by default. The added advantage of the Arctic Mapping Tools is it plots everything in standard polar stereographic meters, meaning you can easily pair plotpsn(lat,lon) with pcolor(x,y,Z) of your favorite dataset (like BedMachine, etc) where x and y are standard meters.

Community Treasure Hunt

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

Start Hunting!

Translated by