栅格面积统计方法.docx
《栅格面积统计方法.docx》由会员分享,可在线阅读,更多相关《栅格面积统计方法.docx(21页珍藏版)》请在冰豆网上搜索。
栅格面积统计方法
栅格面积统计
统计栅格影像的面积,要求栅格数据的坐标必须为平面直角坐标;如不是则需要进行投影转换。
经纬度到平面直角坐标转换
创建一个本地数据库FileGeodatabase,右键菜单New—FileGeoDataBase
为创建的数据库输入名称,如图中的rasterdata,也可使用默认名称
导入要处理的栅格影像。
在上一步创建的rasterdata数据库节点上鼠标右键Import—RasterDatasets
选择要导入的影像数据,然后点击“OK”按钮导入数据。
设置栅格影像的空间投影。
在导入的影像数据节点上鼠标右键,选择Properties
在RaseterDatasetProperties窗口中点击SpatialReference右侧的“Edit”按钮
在弹出的SpatialReferenceProperties窗口中,点击“Select…”按钮;
在BrowseforcoordinateSystem对话框中双击“GeographicCoordinatesystems”进入;
双击“Asia”进入;
选择栅格影像数据的坐标系,如“Xian1980.prj”,点击“Add”按钮完成添加
在ArcMap下加载导入到rasterdata数据库中的栅格影像;
设置投影。
双击ArcToolbox—DataManagementTools—ProjectionsandTransformations—Raster—ProjectRaster
InputRaster:
点击右侧的下拉箭头,选择加载到ArcMap中的影像
OutputCoordinateSystem:
点击右侧的按钮;
在弹出的SpatialReferenceProperties窗口中,点击“Select…”按钮;
在BrowseforcoordinateSystem对话框中双击“ProjectedCoordinatesystems”进入;
选择合适的椭球体,如“GaussKruger”高斯克吕格,双击进入;
选择栅格影像的坐标系,如“xian1980”,双击进入;
根据山歌影像数据经纬度范围等信息,选择合适的坐标带。
如选择“西安80,123度中央经线所在的3度分带”。
点击“OK”按钮,完成栅格影像的投影转换。
投影转换后的数据被重新加载到ArcMap中
点击View菜单下的DataFrameProperties菜单项,
在DataFrameProperties窗口中选择“CoordinateSystem”选项栏,
选择之前投影转换时选中的坐标带,如下图所示。
(predefined—ProjectedCoordinatesystems—GaussKruger—Xian1980—Xian19803DegreeGKCM123E)
选择“General”选项卡,更改下面的Display项为“Meters”,点击“确定”按钮。
此时,地图的坐标系以更改为平面直角坐标系,见ArcMap窗口右下角状态栏的坐标信息便可知。
统计栅格不同颜色像素的面积
更改栅格影像专题图样式。
在投影转换后的栅格影像上点击鼠标右键,选择Properties菜单;
在弹出的LayerProperties对话框中,左侧列表选择”Classified”项,可在右侧的“Classes”项选择级数,点击“确定”按钮。
选择SpatialAnalyst工具栏下拉菜单的“Reclassify…”项;
在弹出的窗口中点击“OK”按钮,按当前的分类方法重分类。
重分类后的栅格影像,如下图所示
为重分类后的影像添加面积字段。
在重分类后的影像上点击鼠标右键,选择“OpenAttributeTable”
属性表窗口中点击“Options”按钮,在弹出的菜单中选择“AddField…”
输入面积字段的名称,如“area”,字段类型选择“Double”,点击“OK”按钮添加字段
获得像素单元格大小。
在重分类后的影像上点击鼠标右键,选择“properties…”,
在弹出的layerproperties窗口中选择“Source”选项栏,CellSize项的值为单元格大小信息。
如下图所示:
2179.961151*2179.961151即为像素单元代表的实际面积。
(由于之前设置的地图坐标单位为米,所以这里计算出来的单元格的面积单位为平方米)
计算每类颜色单元格面积。
在重分类后影像属性表的代表面积的列(如图中的AREA)上点击鼠标右键,选择“FieldCalculator”;
输入如下图所示的表达式,点击“OK”按钮计算面积。
每种颜色像素代表的面积,见下图所示;
(由于之前设置的地图坐标单位为米,所以这里计算出来的单元格的面积单位为平方米)