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

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) ; 

Answer 1

[Zhaoyi Zhang]

File Exchange - MATLAB Central (mathworks.cn)

[Zhaoyi Zhang]

不知道有没有帮助，我现在也在读取NSIDC的数据。