I need to extract latitude and longitude from geotiff file

Question

Aswathi VK 2024-5-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2114176-i-need-to-extract-latitude-and-longitude-from-geotiff-file

评论： Aswathi VK 2024-5-3

I am working with SMAP soil moisture data of downscaled product of 1km. The image is in geotiff format, i couldnt extract latitude and longitude value for the data, while visualising the image the x and y axis appears to be row and coloumn number. Please help me for the same

matlab version: MATLAB R2019a

geotiff file: NSIDC-0779_EASE2_G1km_SMAP_SM_DS_20160101.tif

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

请先登录，再进行评论。

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

Answer 1

Sachin Lodhi 2024-5-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2114176-i-need-to-extract-latitude-and-longitude-from-geotiff-file#answer_1451246

在 MATLAB Online 中打开

Hello Aswathi,

I understand that you want to extract latitude and longitude from geotiff file in Matlab R2019a. For achieving this, you can use the following workaround :

[A,R] = readgeoraster('inputFile.tif');
[x,y] = pixcenters(R, size(A), 'makegrid'); 
info = geotiffinfo('inputFile.tif');
[lat,lon] = projinv(info,x,y);
figure
geoshow(lat,lon,A);
xlabel('Longitude (degrees)')
ylabel('Latitude (degrees)')

Moreover, to fit the figure's axes to exactly match the size of the geotiff file, use the following code in addition to the above snippet :

xlimits = R.XWorldLimits;
ylimits = R.YWorldLimits;
[latlim,lonlim] = projinv(info,xlimits,ylimits); % if using 2021a or later, use 'R.ProjectedCRS' instead of 'info'
xlim(lonlim);
ylim(latlim);

I hope this helps!

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Aswathi VK 2024-5-3

Thank you very much for the quick reply. I tried using that and unfortunately it does not work. This is proving to be quite a difficult problem for me.

请先登录，再进行评论。