Change coordinate system of geotif from WGS 34 to EPSG:32628

Question

Timbro 2023-6-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1979829-change-coordinate-system-of-geotif-from-wgs-34-to-epsg-32628

评论： Timbro 2023-6-14

test.mat

在 MATLAB Online 中打开

Hello

I have a geotif with data projected in WGS 84 (standart) and I want to convert it in EPSG:32628

since .tif files cannot be uploaded here, I saved it in a .mat (using geotiffread and save). Thus the data are in the file test.mat (attached). I want to read this data and export it in a new geotif but with different coordinate system, but I can't :

load test.mat
coordRefSysCode = 32628
geotiffwrite('test',A,R,'CoordRefSysCode',coordRefSysCode)

I get this error :

Error using geotiffwrite

Function GEOTIFFWRITE expected input number 3, R, to be either a 3-by-2 referencing matrix or a scalar map

raster reference object. Instead its type was: map.rasterref.GeographicCellsReference.

How can I do it ?

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

Timbro 2023-6-9

I understand from the error message that R should be different, but then the question is how can I get R as a 3-by-2 referencing matrix or a scalar map raster reference object ?

请先登录，再进行评论。

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

Answer 1

Daniele Sportillo 2023-6-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1979829-change-coordinate-system-of-geotif-from-wgs-34-to-epsg-32628#answer_1255814

在 MATLAB Online 中打开

Hi Timbro,

you can try doing this:

load test.mat
coordRefSysCode = 32628;
M = makerefmat('RasterSize',R.RasterSize,'LatitudeLimits',R.LatitudeLimits,'LongitudeLimits',R.LongitudeLimits,'ColumnsStartFrom', R.ColumnsStartFrom);
geotiffwrite('test',A,M,'CoordRefSysCode',coordRefSysCode);