栅格面积统计方法.docx

栅格面积统计 统计栅格影像的面积，要求栅格数据的坐标必须为平面直角坐标；如不是则需要进行投影转换。 经纬度到平面直角坐标转换 创建一个本地数据库File Geodatabase，右键菜单New—File GeoDataBase 为创建的数据库输入名称，如图中的rasterdata，也可使用默认名称 导入要处理的栅格影像。在上一步创建的rasterdata数据库节点上鼠标右键Import—RasterDatasets 选择要导入的影像数据，然后点击“OK”按钮导入数据。 设置栅格影像的空间投影。在导入的影像数据节点上鼠标右键，选择Properties 在Raseter Dataset Properties窗口中点击Spatial Reference 右侧的“Edit”按钮 在弹出的Spatial Reference Properties窗口中，点击“Select…”按钮； 在Browse for coordinate System对话框中双击“Geographic Coordinate systems”进入； 双击“Asia”进入； 选择栅格影像数据的坐标系，如“Xian 1980.prj”，点击“Add”按钮完成添加 在ArcMap下加载导入到rasterdata数据库中的栅格影像； 设置投影。双击ArcToolbox—DataManagement Tools—Projections and Transformations—Raster—Project Raster Input Raster：点击右侧的下拉箭头，选择加载到ArcMap中的影像 Output Coordinate System：点击右侧的按钮； 在弹出的Spatial Reference Properties窗口中，点击“Select…”按钮； 在Browse for coordinate System对话框中双击“Projected Coordinate systems”进入； 选择合适的椭球体，如“Gauss Kruger”高斯克吕格，双击进入； 选择栅格影像的坐标系，如“xian 1980”，双击进入； 根据山歌影像数据经纬度范围等信息，选择合适的坐标带。如选择“西安80,123度中央经线所在的3度分带”。 点击“OK”按钮，完成栅格影像的投影转换。 投影转换后的数据被重新加载到ArcMap中 点击View菜单下的Data Frame Properties菜单项， 在Data Frame Properties窗口中选择 “Coordinate System”选项栏， 选择之前投影转换时选中的坐标带，如下图所示。（predefined—Projected Coordinate systems—Gauss Kruger—Xian 1980—Xian 1980 3Degree GK CM 123E） 选择“General”选项卡，更改下面的Display项为“Meters”，点击“确定”按钮。此时，地图的坐标系以更改为平面直角坐标系，见ArcMap窗口右下角状态栏的坐标信息便可知。 统计栅格不同颜色像素的面积 更改栅格影像专题图样式。在投影转换后的栅格影像上点击鼠标右键，选择Properties菜单； 在弹出的Layer Properties对话框中，左侧列表选择”Classified”项，可在右侧的“Classes”项选择级数，点击“确定”按钮。 选择Spatial Analyst工具栏下拉菜单的“Reclassify…”项； 在弹出的窗口中点击“OK”按钮，按当前的分类方法重分类。 重分类后的栅格影像，如下图所示 为重分类后的影像添加面积字段。在重分类后的影像上点击鼠标右键，选择“Open Attribute Table” 属性表窗口中点击“Options”按钮，在弹出的菜单中选择“Add Field…” 输入面积字段的名称，如“area”，字段类型选择“Double”，点击“OK”按钮添加字段 获得像素单元格大小。在重分类后的影像上点击鼠标右键，选择“properties…”, 在弹出的layer properties窗口中选择“Source”选项栏，CellSize项的值为单元格大小信息。如下图所示：2179.961151*2179.961151即为像素单元代表的实际面积。（由于之前设置的地图坐标单位为米，所以这里计算出来的单元格的面积单位为平方米） 计算每类颜色单元格面积。在重分类后影像属性表的代表面积的列（如图中的AREA）上点击鼠标右键，选择“Field Calculator”； 输入如下图所示的表达式，点击“OK”按钮计算面积。 每种颜色像素代表的面积，见下图所示； （由于之前设置的地图坐标单位为米，所以这里计算出来的单元格的面积单位为平方