栅格面积统计方法.doc

栅格面积统计统计栅格影像的面积，要求栅格数据的坐标必须为平面直角坐标；如不是则需要进行投影转换。经纬度到平面直角坐标转换创建一个本地数据库File Geodatabase，右键菜单NewFile GeoDataBase为创建的数据库输入名称，如图中的rasterdata，也可使用默认名称导入要处理的栅格影像。在上一步创建的rasterdata数据库节点上鼠标右键ImportRasterDatasets选择要导入的影像数据，然后点击“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数据库中的栅格影像；设置投影。双击ArcToolboxDataManagement ToolsProjections and TransformationsRasterProject RasterInput 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”选项栏，选择之前投影转换时选中的坐标带，如下图所示。（predefinedProjected Coordinate systemsGauss KrugerXian 1980Xian 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）上点击鼠标右键，选择“Fie