MATLAB: How to plot polar stereographic projection

Question

Hello, I have concentration, latitude and longitude data as [m x n] matrices from NSIDC (<http://nsidc.org/data/docs/daac/nsidc0079_bootstrap_seaice.gd.html)>. The data are gridded over polar stereographic grid. Can someone help me to plot sea ice concentration over the Arctic?

Thanks very much, BD

Answer 1

I did this for the Antarctic in a function called seaice.

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.