本次鄱阳湖高程—面积、高程—容积关系计算，采用的方法：一是利用DEM进行基本地形因子计算；二是MAPGIS线画图分图幅量算。最终成果采用DEM计算结果。MAPGIS分图幅、逐级量算各等高线所包围的面积，其成果用于与DEM计算结果的校验和合理性分析。

3.2.4.1 基于DEM的面积、容积计算（简称DEM计算）

DEM是在某一投影平面（如高斯投影平面）上规则格网点的平面坐标（X，Y）及高程（Z）的数据集。DEM作为地球空间框架数据的基本内容和其他各种地理信息的载体，是各种地学分析的基础数据，自然也是GIS的基本内容。本次鄱阳湖基础地理测量根据实测1∶10000数字地形图，建立鄱阳湖1∶10000DEM数据覆盖鄱阳湖测区范围，共计439幅。利用DEM计算面积、容积的图解模型见图3.2-1。

1.GRID数据整理

（1）通过ASCII to Raster命令将grd格式ArcInfo Grid交换格式文件转换为GRID格式的Raster数据集。

（2）通过Define Projection命令添加投影。

图3.2-1　DEM计算面积和容积的图解模型

（3）通过Mosaic命令将439幅数据合成一幅GRID格式的Raster数据集。

2.制作计算边界范围面

从DEM中提取21m等高线，与数字线划地图DLG的21m等高线、堤防等信息进行对比分析，局部参考数字正射影像图DOM资料，先形成各计算区域的边界线，然后通过ArcInfo Workstation将边界线转换为区域面图。

3.生成不同区域的独立DEM数据(www.daowen.com)

（1）边界范围面保存成shp文件，打开属性表，修改边界面字段“ID”的值为1，通过Polygon to Raster_conversion命令转换成栅格数据，指定栅格大小为1。

（2）通过栅格计算器Times_sa命令，执行边界范围栅格图层和用以剪切的整个鄱阳湖的栅格之间的相乘运算，通过这种方法生成各类型分区的栅格数据。

4.面积容积计算

（1）通过Reclassify_sa命令将栅格数据进行重分类，Classify分类标准为固定0.5m一级，从-20m到100m共制作分级数为241级。从重分类后的DEM上导出数据表dbf文件。

（2）通过Delphi编程实现对指定文件夹下的所有dbf文件的自动循环计算不同高程级的面积和容积，从而取得鄱阳湖区实测21m高程以下范围及各类分区高程—面积、高程—容积计算成果。

3.2.4.2 基于MapGIS面积量算（简称MapGIS量算）

Map GIS主要实现制图、空间分析、属性库管理等功能，分为输入、编辑、输出、空间分析、库管理、实用程序六大部分。

（1）采用成图软件MAPGIS的部分存文件的功能，将鄱阳湖区实测21m以下范围内4～21m等高线提取出来，并将有关的地性线（如堤坝、陡坎等）拷贝出来后叠加图幅裁截线生成Map GIS的*.wl线文件；在输入编辑中对线文件（*.wl）进行线转为可以生成面的弧段，进行拓扑重建生成面文件（*.wp）存盘。

（2）将面文件叠加到图形中，编辑区属性结构，在区属性结构中增加高程属性。

（3）利用自动区标注功能选择相应的属性字段进行区标注，按照地形图上的高程点、等高线、陡坎等高程属性给每个面赋相应值，如4m以下的面赋4m，4～5m面赋5m等。

（4）将赋好高程属性的面按属性→表格导出高程与面积值，然后统计每幅图各级等高线所包围的面积，最终累计得出鄱阳湖区实测21m以下范围高程面积成果。