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