相同投影带有/无带号的栅格转换,“平移”与“重投影”是不是等效的?

手头在处理一批没有定义过投影带的影像数据(ERDAS IMG),因为没定义过投影,所以是坐标系是undefined的。联系过数据发布方,影像对应的平面投影应该是CGCS2000_3_Degree_GK_CM_111E,也就是没带号的。
我知道ArcMap有动态投影的功能,但是为了数据在我们其他的软件中能够和已有数据(有带号)正常叠加,打算全部转换为带号37的投影坐标,也就是CGCS2000_3_Degree_GK_Zone_37。
目前看来有两个办法:
1、因为原数据没定义过投影,可以直接用Shift工具将影像的X坐标平移+37000000个单位,然后将投影定义为CGCS2000_3_Degree_GK_Zone_37。
2、先用Define Projection将影像定义为CGCS2000_3_Degree_GK_CM_111E,然后使用Project Raster将影像重新投影到CGCS2000_3_Degree_GK_Zone_37。
 
问题:Project Raster有一个重采样的参数,但是这样一个操作应该仅仅平移,不涉及重采样的问题。Shift可以确保不会执行任何重采样或扭曲操作,不知道Project Raster工具在进行这样的操作的时候也会进行重采样,还是仅作平移和投影重设处理?两种方法哪一种比较稳妥?
已邀请:

慕晓燕 - 长期ArcGIS爱好者

赞同来自: yxiegis nfc 杨乐

看了你的题目,这两种投影变换的方法都可以使用的。 如果为了避免栅格数据重采样,建议使用 Shift 方法,然后重新指定坐标系。因为这个坐标系变换的需求本质就是坐标轴的伪东平移。
 
 
采用 Project Raster 方法也可以的,为了避免重采样,最好的做法是 执行工具前,设置 Output Cell size 、Snap Raster 或者 Registration Point  参数,来限定输出栅格的像元大小和范围。像元大小和范围的变化是导致重采样的根本原因。
 
我用相同的数据测试了一下,也是前者的效率更高。建议做个Model,调用 shift,Define Projection 工具批量实现。

徐珂 - ArcGIS for Desktop

赞同来自: yxiegis

Project Raster是会进行重采样。如果是离散的数据,如土地利用数据,使用最邻近法,是不会更改像元值的。对于连续的数据,如影像、DEM数据,重采样后会更改像元值。
 
您说的这种情况,如果不想改变像元值,那么可以使用您说的第一种方法。如果能接受像元值的变动(当然,变化不大),可以使用第二种方法。您这个只是换带,所以可以使用Shift。如果是基准面变了,那么只能使用Project Raster了。
 
PS:对于换带的这种情况,如果数据带投影,同事还提出了可以使用程序修改投影文件中坐标的想法,tif格式的是tfw。栅格数据集的坐标文件请参考:http://desktop.arcgis.com/zh-c ... s.htm
tfw.jpg

 

要回复问题请先登录注册