改变geotiff地理空间数据文件。

10视图(30天)
你好,
我的目标是改变geotiff的地理空间数据文件。
换句话说,我想搬迁到另一个coodrinates形象。
更准确地说,我们有一个试车跑道,我想带一块一个真实世界的道路。
所以当我走鹃导入图像,它我需要它的地方。
我是怎么做的:
1。我读到单独的图像和数据
data = geotiffinfo(文件名);
形象= imread(文件名);
2。我定义了一个矩阵,我需要申请的形象
矩阵= (0.990950218730870 - 0.136163909329197 0
-0.136163909329197 - 0.990950218730870 0
741682.197901774 - -43966.3723596139 1);
3所示。我创建了一个 tform 对象
tform = affine2d(矩阵);
4所示。我创建了一个对象,将空间数据和像素。
xWorldLimits =[数据。BoundingBox(1,1)数据。BoundingBox (2, 1)];
yWorldLimits =[数据。BoundingBox(1、2)数据。BoundingBox (2, 2)];
ref_data = imref2d(大小(图片)、xWorldLimits yWorldLimits);
5。我应用矩阵的图像和地理空间数据
[upd_image, upd_data] = imwarp(图像、ref_data tform);
6。引用对象需要创建一个新的tiff图像。
Rmap = maprefcells (upd_data。XWorldLimits upd_data。YWorldLimits、大小(upd_image),
“ColumnsStartFrom”,“北”,“RowsStartFrom”,“西方”);
7所示。最后创建一个文件。
new_filename = [“ra_transformed”“.tif”];
geotiffwrite (new_filename upd_image Rmap,“CoordRefSysCode”,25832);
最后,它没有工作。
图片在另一个方向旋转,移动不精确。
这是接近但不完全。
矩阵是正确的,我已经检查过两次。
有什么我错过了吗?
或者分割的整体方法,改变并将它是错误的吗?

答案(1)

近日,Abdrakhimov
近日,Abdrakhimov 2022年2月9日
我已经解决了这个问题。
问题是在不同坐标方向UTM和Matlab。
把角(旋转)的另一个迹象。

s manbetx 845


释放

R2021b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!