Change only coordinates of a georeferenced tif file from utm to degree

Question

Jonas Müller 2022-7-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1756410-change-only-coordinates-of-a-georeferenced-tif-file-from-utm-to-degree

评论： Jonas Müller 2022-7-9

I want to open a bathymetry as tif file and only changed the UtM coordinates into degree. The problem is that in info are important information for the tif and I don't know how to put this information in geotiffwrite. So the new tif file changed more than only the coordinate transformation. Is it possible to change only the coordinates?

cheers Jonas

Code:

infilename = 'helgoland_bathy.tif';
[A,R] = readgeoraster(infilename);
info = geotiffinfo(infilename);
geoTags = info.GeoTIFFTags.GeoKeyDirectoryTag;
x =R.XWorldLimits;
y=R.YWorldLimits;
utmzone=['32 U';'32 U'];
[LatN, LonN]=utm2deg(x,y,utmzone); 
R.XWorldLimits = LonN';
R.YWorldLimits =LatN';
geotiffwrite('testtiforiginal_grad.tif',A,R,'GeoKeyDirectoryTag',geoTags);

sfsdfds

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

请先登录，再进行评论。

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

Answer 1

KSSV 2022-7-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1756410-change-only-coordinates-of-a-georeferenced-tif-file-from-utm-to-degree#answer_1003595

You may use the following file exchange function to convert utm to degree.

https://in.mathworks.com/matlabcentral/fileexchange/45699-ll2utm-and-utm2ll

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

KSSV 2022-7-9

Refer this to write into tiff file.

https://in.mathworks.com/matlabcentral/answers/355554-problem-creating-a-geotiff-from-a-non-georeferenced-image

Jonas Müller 2022-7-9

在 MATLAB Online 中打开

ok but how can i write all information from the old tiff to the new one. Also the infromation in the info variable? Something like this:

infilename = 'helgoland_bathy.tif';
[A,R] = readgeoraster(infilename);
info = geotiffinfo(infilename);
geoTags = info.GeoTIFFTags.GeoKeyDirectoryTag;
x =R.XWorldLimits;
y=R.YWorldLimits;
utmzone=['32 U';'32 U'];
[LatN, LonN]=utm2deg(x,y,utmzone); 
R.XWorldLimits = LonN';
R.YWorldLimits =LatN';
geotiffwrite('testtiforiginal_grad.tif',A,R,info);