栅格面积统计方法

《栅格面积统计方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、栅格面积统计统计栅格影像的面积，要求栅格数据的坐标必须为平面直角坐标；如不是则需要进行投影转换。经纬度到平面直角坐标转换创建一个本地数据库FileGeodatabase，右键菜单New—FileGeoDataBase为创建的数据库输入名称，如图中的rasterdata，也可使用默认名称导入要处理的栅格影像。在上一步创建的rasterdata数据库节点上鼠标右键Import—RasterDatasets选择要导入的影像数据，然后点击“OK”按钮导入数据。设置栅格影像的空间投影。在导入的影像数据节点上鼠标右键，选择Properties在RaseterDatasetPropertie

2、s窗口中点击SpatialReference右侧的“Edit”按钮在弹出的SpatialReferenceProperties窗口中，点击“Select…”按钮；在BrowseforcoordinateSystem对话框中双击“GeographicCoordinatesystems”进入；双击“Asia”进入；选择栅格影像数据的坐标系，如“Xian1980.prj”，点击“Add”按钮完成添加在ArcMap下加载导入到rasterdata数据库中的栅格影像；设置投影。双击ArcToolbox—DataManagementTools—ProjectionsandTransform

3、ations—Raster—ProjectRasterInputRaster：点击右侧的下拉箭头，选择加载到ArcMap中的影像OutputCoordinateSystem：点击右侧的按钮；在弹出的SpatialReferenceProperties窗口中，点击“Select…”按钮；在BrowseforcoordinateSystem对话框中双击“ProjectedCoordinatesystems”进入；选择合适的椭球体，如“GaussKruger”高斯克吕格，双击进入；选择栅格影像的坐标系，如“xian1980”，双击进入；根据山歌影像数据经纬度范围等信息，选择合适的坐标

4、带。如选择“西安80,123度中央经线所在的3度分带”。点击“OK”按钮，完成栅格影像的投影转换。投影转换后的数据被重新加载到ArcMap中点击View菜单下的DataFrameProperties菜单项，在DataFrameProperties窗口中选择“CoordinateSystem”选项栏，选择之前投影转换时选中的坐标带，如下图所示。（predefined—ProjectedCoordinatesystems—GaussKruger—Xian1980—Xian19803DegreeGKCM123E）选择“General”选项卡，更改下面的Display项为“Meters

5、”，点击“确定”按钮。此时，地图的坐标系以更改为平面直角坐标系，见ArcMap窗口右下角状态栏的坐标信息便可知。统计栅格不同颜色像素的面积更改栅格影像专题图样式。在投影转换后的栅格影像上点击鼠标右键，选择Properties菜单；在弹出的LayerProperties对话框中，左侧列表选择”Classified”项，可在右侧的“Classes”项选择级数，点击“确定”按钮。选择SpatialAnalyst工具栏下拉菜单的“Reclassify…”项；在弹出的窗口中点击“OK”按钮，按当前的分类方法重分类。重分类后的栅格影像，如下图所示为重分类后的影像添加面积字段。在重分类后的影

6、像上点击鼠标右键，选择“OpenAttributeTable”属性表窗口中点击“Options”按钮，在弹出的菜单中选择“AddField…”输入面积字段的名称，如“area”，字段类型选择“Double”，点击“OK”按钮添加字段获得像素单元格大小。在重分类后的影像上点击鼠标右键，选择“properties…”,在弹出的layerproperties窗口中选择“Source”选项栏，CellSize项的值为单元格大小信息。如下图所示：2179.961151*2179.961151即为像素单元代表的实际面积。（由于之前设置的地图坐标单位为米，所以这里计算出来的单元格的面积单位为

7、平方米）计算每类颜色单元格面积。在重分类后影像属性表的代表面积的列（如图中的AREA）上点击鼠标右键，选择“FieldCalculator”；输入如下图所示的表达式，点击“OK”按钮计算面积。每种颜色像素代表的面积，见下图所示；（由于之前设置的地图坐标单位为米，所以这里计算出来的单元格的面积单位为平方米）