根据湖底DEM计算湖泊蓄水量

1.湖泊面积的计算

S:湖泊面积

N：湖泊所在像元个数

R:影像空间分辨率（比如30m)

在ENVI里统计像元个数，是在阈值确定后得到二值图后，quickstates,或者做密度分割后，统计个数。

2.形成湖泊的shp文件

（1）得到湖泊的二值图在NEVI classic里加载出来，

（2）选择Vector--Raster to Vector，获得evf格式的矢量数据；

（3）在矢量文件窗口，选择file--Export Layer to Shapefile,d得到shp格式的矢量文件

3.由得到的湖泊面积的shp加载到Arcgis中，面转线得到湖的边界线（数据管理工具--要素--面转线），属性表面积降序排列后反向选择删除小的面积--做缓冲区，地理处理--缓冲区--左右都缓冲15m(FULL)

4.用缓冲的边界裁剪DEM：Spatial...--提取分析--按掩膜提取--提取完后打开属性，平均值得到湖面的平均水位

5.用整个湖面的shp去裁剪DEM，还是提取分析按掩膜提取，得到湖底的DEM，然后选择3D Analyst 工具--功能性表面--表面体积--输入表面是刚才用shp裁剪的湖底的DEM（不是原始的DEM），因为水量是从湖底往上的体积，所以选择ABOVE，平面高度是湖面的平均水位高程值，也就是4中的平均值，确定即可，得到的表最后一格就是水量。