How can I get the sea ice concentration data corresponding to the latitude and longitude from the geotiff file?

1 view (last 30 days)
Hello,
I got daily sea ice concentration data from https://nsidc.org/data/G02135/versions/3(please see attachment), but the file is in geotiff format, and I don’t know about geotiff format. I don't know how to get sea ice concentration data corresponding to longitude and latitude from geotiff?
Hope to get your help! Thank you very much!
Bei.

Accepted Answer

KSSV
KSSV on 18 Oct 2021
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) ;
figure(1)
mapshow(A,R)
figure(2)
pcolor(X,Y,A);
shading interp ;
  6 Comments
Bei JIANG
Bei JIANG on 18 Oct 2021
Now I know the projection is “NSIDC Polar Stereographic Projection", and I also know this link provided transformations to geographic coordinates from map coordinates: polar stereographic coordinate transformation (map to lat/lon) - File Exchange - MATLAB Central (mathworks.cn)
However, I dont konw if the X,Y in this command (below) is the XY you provided me in [X,Y] = meshgrid(x,y)?[LAT,LON]=POLARSTEREO_INV(X,Y,EARTHRADIUS,ECCENTRICITY,LAT_TRUE,LON_POSY)
I look forward to hearing from you very much, it is very important to me. Thank you very much !
Bei JIANG
Bei JIANG on 19 Oct 2021
There are two lines wrong in your answer, corrected as follows :
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.XWorldLimits(2),n) ;
y = linspace(R.YWorldLimits(1),R.YWorldLimits(2),m) ;
[X,Y] = meshgrid(x,y) ;
figure(1)
mapshow(A,R)
figure(2)
pcolor(X,Y,A);
shading interp ;

Sign in to comment.

More Answers (0)

Products

Community Treasure Hunt

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

Start Hunting!