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

3 次查看（过去 30 天）

显示更早的评论

Bei JIANG 2021-10-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1565991-how-can-i-get-the-sea-ice-concentration-data-corresponding-to-the-latitude-and-longitude-from-the-ge

评论： Bei JIANG 2021-10-19

采纳的回答： KSSV

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）, 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.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

KSSV 2021-10-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1565991-how-can-i-get-the-sea-ice-concentration-data-corresponding-to-the-latitude-and-longitude-from-the-ge#answer_810646

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 个评论
显示 4更早的评论隐藏 4更早的评论

Bei JIANG 2021-10-18

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 2021-10-19

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 ;

请先登录，再进行评论。

类别

Mathematics and Optimization Mapping Toolbox Coordinate Reference Systems

Radar Mapping Toolbox Coordinate Reference Systems

在 Help Center 和 File Exchange 中查找有关 Coordinate Reference Systems 的更多信息

产品

MATLAB

Community Treasure Hunt

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

Start Hunting!

Translated by