在拼接2010年的N50_30和N51_30的太湖区域landcover拼接时(数据来源于全球地表覆盖-数据浏览和成果查询),envi拼接不行,改成arcGIS拼接(data manage tools-栅格-镶嵌至新栅格),发现各种拼接都有问题,具体问题和原因如下:
(1)中间有原图黑边:因为原图附近有黑边,选了First或last拼接法时,都会留下黑边;
(2)中间有款黑边:用了blend/mean拼接;
(3)中间重叠区域地物分类异常:用了sum拼接,相当于把重叠区域的value相加了,比如value30的地物叠了两次变成value60的地物了;
(4)没有输出TIFF:“具有扩展名的栅格数据集名称”这里没有给文件名加.tif
(5)中间有条虚线:用max和min拼接时都是这样,最后发现是原图左边数据有虚线状的地物是错分的;
解决方案:最后尝试多种方案,发现既然右图没有问题,那就用右图叠加左图,用右图和First拼接法就行,那么只剩下去黑边了,去黑边的方法就用“Spatial Analyst Tools-重分类-重分类”的方法,在重分类时直接不填写黑边对应的value值,然后勾选上“将缺失值更改为NoData (可选)”,黑边就会自动被去除,然后直接First拼接就行;
(原图层还有一些小问题需要处理:1.左右两边的原图同样地物对应的value不一致,需要先重分类到同一value;2.重分类后的地物颜色就变了,可以先看原图的rgb值然后一个个的在图层属性里改,改完后导出为色彩映射表,后面就用这个映射表导入就可以改图层颜色了,但要注意图层种类的一致性;3.原图层没有landuse这个字段,要加上;4.原图层没有属性表,在属性里选择”符号系统-唯一值“时会提示创建)
具体步骤:1、建立原图属性表:在属性里选择”符号系统-唯一值“时会提示创建
2、重分类左右原图”N50_30“和”N51_30“,使其同地物的value值一样,然后添加字段landuse并填写字段内容(耕地,水体...)
3、去除原图黑边:重分类时删除背景黑边对应的那一行,不给它赋值value,然后勾选上“将缺失值更改为NoData (可选)
(选唯一;然后新值删除,什么也不填;填新值;最后勾选上(Nodata))
4、处理后在属性里导入色彩映射表,更改回颜色
5、两张图都以上处理后,开始镶嵌(注意的是处理时“新值”那里两个图要同地物用同一个value,比如都是水体,那么新值都取60)
(重采样前后对比图:不过图是2000年的,但是流程是差不多的,后面一张图的黑边被去掉了)
6、镶嵌到新栅格:
Data management tools-栅格-栅格数据集-镶嵌到新栅格
7、裁剪感兴趣区域:掩膜裁剪
Spatial analyse tools-提取分析-按掩膜提取
部分参考:土地利用分类日志1——分幅土地利用图拼接融合step by step教程_wenhao048的博客-CSDN博客_土地利用数据拼接