相同投影带有/无带号的栅格转换，“平移”与“重投影”是不是等效的？
手头在处理一批没有定义过投影带的影像数据（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工具在进行这样的操作的时候也会进行重采样，还是仅作平移和投影重设处理？两种方法哪一种比较稳妥？
我知道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工具在进行这样的操作的时候也会进行重采样，还是仅作平移和投影重设处理？两种方法哪一种比较稳妥？
2 个回复
慕晓燕 - 长期ArcGIS爱好者
赞同来自: yxiegis 、nfc 、杨乐
采用 Project Raster 方法也可以的，为了避免重采样，最好的做法是 执行工具前，设置 Output Cell size 、Snap Raster 或者 Registration Point 参数，来限定输出栅格的像元大小和范围。像元大小和范围的变化是导致重采样的根本原因。
我用相同的数据测试了一下，也是前者的效率更高。建议做个Model，调用 shift，Define Projection 工具批量实现。
徐珂 - ArcGIS for Desktop
赞同来自: yxiegis
您说的这种情况，如果不想改变像元值，那么可以使用您说的第一种方法。如果能接受像元值的变动（当然，变化不大），可以使用第二种方法。您这个只是换带，所以可以使用Shift。如果是基准面变了，那么只能使用Project Raster了。
PS：对于换带的这种情况，如果数据带投影，同事还提出了可以使用程序修改投影文件中坐标的想法，tif格式的是tfw。栅格数据集的坐标文件请参考：http://desktop.arcgis.com/zh-c ... s.htm
要回复问题请先登录或注册
发起人
相关问题
问题状态