I need to generate a mapping toolbox projection structure for the Scientific Committee for Antarctic Research (SCAR) standard projection. It is a polar stereographic projection with a standard parallel at 71S. None of the azimuthal projections accept standard parallels, so how do I calculate the scale factor I need to make -71 the latitude of true scale?
MATLAB: Need mapping toolbox mstruct for SCAR stereographic projection
latlonmapping tollboxMapping Toolboxpolar stereographicprojectionsscar
Related Solutions
When you combine plotm (or other *m.m functions in Mapping Toolbox) and geoshow, you need to define axis in GEOSHOW function. Also, set DisplayType option as 'point' in GEOSHOW (default value is 'line').
ax = worldmap('angola');load coastlinessetm(ax, 'Origin', [0 20 0])plotm(coastlat,coastlon);lat = -8;lon = 13;geoshow(ax, lat, lon, 'DisplayType', 'point', 'Marker', 'o')
You'll want to download psn25lats_v3.dat and pss25lats_v3.dat, which contain lat,lon grid information. The seaice function imports lat,lon grids for the southern hemisphere like this:
fid = fopen('pss25lons_v3.dat');lon = fread(fid,[316 332],'long','ieee-le')/100000;fclose(fid);fid = fopen('pss25lats_v3.dat');lat = fread(fid,[316 332],'long','ieee-le')/100000;fclose(fid);You can then use Andrew Bliss' polarstereo_fwd to convert lat,lon to x,y and plot with pcolor or imagesc.
Hope this helps.
Best Answer