我理解的错误消息,R应该是不同的,但问题是我怎样才能R 3×2引用矩阵或标量光栅引用对象地图吗?
改变坐标系统的geotif WGS 34 EPSG: 32628
45视图(30天)
显示旧的评论
你好
我有一个geotif数据预计在WGS 84 (standart),我想把它EPSG: 32628
因为.tif文件不能上传在这里,在一个.mat(使用我救了它
geotiffread并保存)。因此,数据文件中的测试。垫(附加)。我想读数据和导出一个新的geotif但不同的坐标系统,但我不能:
负载test.mat
coordRefSysCode = 32628
geotiffwrite (“测试”R,“CoordRefSysCode”coordRefSysCode)
我得到这个错误:
错误使用geotiffwrite
函数GEOTIFFWRITE预期输入数字3,R, 3×2引用矩阵或标量映射
光栅的参考对象。相反,它的类型是:map.rasterref.GeographicCellsReference。
我能怎么做?
答案(1)
Daniele Sportillo
2023年6月14日17:22
你好哈马,
你可以试试这样做:
负载test.mat
coordRefSysCode = 32628;
M = makerefmat (“RasterSize”R.RasterSize,“LatitudeLimits”R.LatitudeLimits,“LongitudeLimits”R.LongitudeLimits,“ColumnsStartFrom”,R.ColumnsStartFrom);
geotiffwrite (“测试”米,“CoordRefSysCode”,coordRefSysCode);
您可以验证矩阵M产生相同的GeographicCellsReference对象与下面的代码:
R2 = refmatToGeoRasterReference (M, R.RasterSize);
正常工作吗?