I need to extract latitude and longitude from geotiff file
17 次查看(过去 30 天)
显示 更早的评论
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 个评论
回答(1 个)
Sachin Lodhi
2024-5-2
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!
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!