全国高分辨率土地利用数据服务 土地利用数据服务 土地覆盖数据服务 坡度数据服务 土壤侵蚀数据服务 全国各省市DEM数据服务 耕地资源空间分布数据服务 草地资源空间分布数据服务 林地资源空间分布数据服务 水域资源空间分布数据服务 建设用地空间分布数据服务 地形、地貌、土壤数据服务 分坡度耕地数据服务 全国大宗农作物种植范围空间分布数据服务
多种卫星遥感数据反演植被覆盖度数据服务 地表反照率数据服务 比辐射率数据服务 地表温度数据服务 地表蒸腾与蒸散数据服务 归一化植被指数数据服务 叶面积指数数据服务 净初级生产力数据服务 净生态系统生产力数据服务 生态系统总初级生产力数据服务 生态系统类型分布数据服务 土壤类型质地养分数据服务 生态系统空间分布数据服务 增强型植被指数数据服务
多年平均气温空间分布数据服务 多年平均降水量空间分布数据服务 湿润指数数据服务 大于0℃积温空间分布数据服务 光合有效辐射分量数据服务 显热/潜热信息数据服务 波文比信息数据服务 地表净辐射通量数据服务 光合有效辐射数据服务 温度带分区数据服务 山区小气候因子精细数据服务
全国夜间灯光指数数据服务 全国GDP公里格网数据服务 全国建筑物总面积公里格网数据服务 全国人口密度数据服务 全国县级医院分布数据服务 人口调查空间分布数据服务 收入统计空间分布数据服务 矿山面积统计及分布数据服务 载畜量及空间分布数据服务 农作物种植面积统计数据服务 农田分类面积统计数据服务 农作物长势遥感监测数据服务 医疗资源统计数据服务 教育资源统计数据服务 行政辖区信息数据服务
Landsat 8 高分二号 高分一号 SPOT-6卫星影像 法国Pleiades高分卫星 资源三号卫星 风云3号 中巴资源卫星 NOAA/AVHRR MODIS Landsat TM 环境小卫星 Landsat MSS 天绘一号卫星影像
空间分析中的掩膜提取和栅格处理中的裁剪。
在ArcGIS中,空间分析中的掩膜提取是通过ArcToolbox->Spatial Analyst Tools->Extraction->Extract byMask实现的。在Extract by Mask对话框中设置好输入输出参数以后,直接点击OK即可得到掩膜后的栅格数据。然而看似简单的过程实则包含了几个比较棘手的问题。
裁切前后两个栅格的像元并不一定是完全重叠的,有可能错位,也就是说对应像元的中心点不是一个点,如下图1所示。这时裁切后像元的值可能是裁切前相邻像元的值,所以对应值就可能不一致,尤其在高程变化比较大的地方表现得更为明显。
图1 栅格数据裁切像元不完全重叠
若想让裁切前后的每个像元都完全重叠,可在裁切界面的Environment Settings里设置snap raster成裁切前的数据即可。
用不规则矢量多边形裁切栅格数据时,Extract by Mask对话框中没有提供设置Nodata的选项,这导致有时会出现裁剪后的pixel depth由原来的8bit变成了16bit(如下图2所示)。
原因为Unsigned 8 bit栅格的值的范围是0~255,当用不规则polygon裁切时,为了保留mask的准确的形状,系统会在polygon外的部分生成nodata区域,而nodata在栅格里是以某个极大值存储的(一般是65535)。当裁切部分的值已经占满0~255的所有值时,nodata的值只能是255以外的某个值了,系统自动会将栅格改为Unsigned 16 bit,导致数据量倍增。
解决方法为在arctoolbox中使用栅格计算器工具,运用setnull命令将选定(个数最少的那类象元值)的值(如0值)设为nodata,然后再裁剪。
在栅格数据中,除了符合比特深度(Pixel Depth)的值之外,用于标记此处Nodata的值----就是空值。产生空值的原因有很多种,有可能是在数据采集的时候,传感器没有记录下某个位置的反射值,或者是由于噪声。然而最常见的则是对卫星影像做了地理配准后,影像的覆盖范围变为菱形/平行四边形所造成的边角是黑色的情况。
用Identify工具查看黑边所示的像元值,通常会发现用0或者极小值来标记Nodata,而在ArcGIS中0值是一个有效值而不是空值。出于这个原因,在渲染这类影像的时候,ArcGIS会将0值渲染成黑色,而不是作为空值处理成透明色。那么我们如何消除这个黑边的问题呢?
在ArcGIS中,我们可以在图层属性Symbology页下看到这样的一个设置:Display Background value…as…。通过输入0,0,0,以及将颜色设置成透明,就可以达到让黑边消失的渲染效果(如下图3)。
图3 在Symbology中设置黑边
以上黑边处理是对单一的栅格数据,在实际应用中有时会遇到Mosaic Dataset这种数据模型,ArcGIS中提供了一个解决Mosaic Dataset中黑边问题的工具Define Mosaic DatasetNoData。在工具中可以对每一个波段单独设置空值,完成后也可以得到理想的黑边变透明的效果。
然而以上处理并没有从根本上裁切掉那些黑边,要从根本上解决问题可以通过导出的方式将黑边设置为空。具体操作就是在Export Data的时候,在NoData的选项中输入0,0,0(如图4)。看似问题得到了解决,然而这种做法会将影像中那些本来就是0值的像元也抹掉,所以这个办法要慎用。
图4 在导出栅格数据时设置Nodata
那么目前有没有既有效又彻底的方法来解决不规则矢量多边形裁切栅格数据时出现的上述问题呢,答案依然是裁剪,接下来介绍栅格处理中的裁剪方法。
在ArcGIS中,空间分析中的掩膜提取是通过ArcToolbox->Data ManagementTools->Raster->Raster Processing->Clip实现的。在Clip对话框中输入待裁切的栅格数据以及矢量多边形数据等参数,直接点击OK即可完成裁剪。通常输入矢量多边形时,系统会自动获取该矢量多边形的最小外包矩形的范围作为实际裁剪的区域,而我们只需要勾选Use Input Features forCliping Geometry(Optional)选项(如图5所示)就可以保证裁剪出来的栅格区域就是我们真正想要的。同时还需要注意的是,Clip对话框中同样还提供了Nodata Value选项,我们可以将那些需要设置为空值的像元过滤掉。
图5 Clip中参数设置
通过对比以上两种矢量裁切方法,我们发现栅格处理中的裁剪方法既简单又有效。对于习惯用Extraction by Mask方法的用户来说,与其舍近求远不如直接用Data Management Tool中的Clip工具。