How to convert EPSG:3413 WGS 84 / NSIDC Sea Ice Polar Stereographic coordinate system to EPSG:4326 WGS 84 coordinate system?

Question

Bei JIANG 2021-10-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1566456-how-to-convert-epsg-3413-wgs-84-nsidc-sea-ice-polar-stereographic-coordinate-system-to-epsg-4326-w

评论： Zhaoyi Zhang 2021-12-7

N_20180101_concentration_v3.0.zip

Hello，

I got daily sea ice concentration data from https://nsidc.org/data/G02135/versions/3（please see attachment).I got the XY data using the following code. However, I want to convert them into lat (range 90~-90) and lon (range 180~-180).

The projection of N_20180101_concentration_v3.0.tif is Table 3 in this link:https://nsidc.org/data/polar-stereo/ps_grids.html. I can use epsg.io to achieve, but it can only be converted one by one, which is very troublesome. How to use matlab code to achieve the conversion?

Looking forward to your reply, thank you very much.！

Bei

filename = 'N_20180101_concentration_v3.0.tif' ;
[A,R] = readgeoraster(filename) ;  % If not try geotiffread
[m,n] = size(A) ; 
x = linspace(R.XWorldLimits(1),R.XIntrinsicLimits(2),n) ; 
y = linspace(R.YWorldLimits(1),R.YIntrinsicLimits(2),m) ; 
[X,Y] = meshgrid(x,y) ; 