基于多源遥感影像数据的国土开发密度三维综合评价方法与流程

技术特征：

1.一种基于多源遥感影像数据的国土开发密度三维综合评价方法，其特征在于，包括步骤：

对待评价区域的遥感图像进行预处理，提取待评价区域的建设用地空间分布信息，根据待评价区域的建设用地空间分布信息确定待评价区域的开发强度和开发紧凑度；

利用夜间灯光影像数据，模拟建设用地之上的经济密度空间分布，用以表示立体开发程度；

从开发强度、开发紧凑度、立体开发程度三个维度建立国土开发密度三维综合评价指标体系；

采用空间直角坐标系对评价指标进行三维综合表达。

2.根据权利要求1所述的基于多源遥感影像数据的国土开发密度三维综合评价方法，其特征在于，

所述利用夜间灯光影像数据模拟建设用地之上的经济密度空间分布，步骤包括：

获取待评价区域的夜间灯光影像，并对夜间灯光影像进行修正；

利用待评价区域建设用地边界切割相应区域经过修正的夜间灯光影像，获得待评价区域建设用地范围内的夜间灯光数据；

计算待评价区域内部各个次级空间单元建设用地范围内的夜间灯光辐射总量，获取待评价区域建设用地夜间灯光辐射总量的空间分布信息：

设RAD_i为第i个次级空间单元建设用地范围内的夜间灯光辐射总量，S_i为第i个次级空间单元的建设用地面积，BRI_i是第i个次级空间单元夜间灯光数据的稳定灯光值部分派生出的稳定灯光强度数据；

RAD_i＝BRI_i×S_i

模拟建设用地之上的非农GDP空间分布：利用夜间灯光辐射总量与非农GDP的空间映射关系，将待评价区域非农GDP统计值展布到各个次级空间单元上；

设GDP_fi’为第i个次级空间单元的非农GDP估算值，GDP_f^*为整个区域的非农GDP统计值，RAD_all为整个区域建设用地范围内的灯光辐射总量；

${{\mathrm{GDP}}_{f_i}}^{\′}={{\mathrm{GDP}}_f}^{*}\×\frac{{\mathrm{RAD}}_i}{{\mathrm{RAD}}_{all}}$

测算建设用地之上的经济密度空间分布；

设ED_i为第i个空间单元建设用地之上的经济密度。

${\mathrm{ED}}_i={{\mathrm{GDP}}_{f_i}}^{\′}/S_i$

3.根据权利要求2所述的基于多源遥感影像数据的国土开发密度三维综合评价方法，其特征在于，

所述开发强度，用以下方法计算：设I_i为第i个空间单元的开发强度，A_i为第i个空间单元的建设用地面积，S_i代表第i个空间单元的总面积；

I_i＝A_i/S_i

所述开发紧凑度，用Richardson和Gibbs的经典面积周长法计算，设C_i为第i个空间单元的开发紧凑度，P_i为第i个空间单元的建设用地总周长；

$C_i=2\sqrt{\mathrm{\π}A}/P_i$

所述立体开发程度，用建设用地之上的经济密度的模拟结果ED_i表示第i个空间单元的立体开发程度。

4.根据权利要求3所述的基于多源遥感影像数据的国土开发密度三维综合评价方法，其特征在于，

所述遥感图像来自于Landsat TM数据，对所述遥感图像的预处理方法包括：

对遥感图像进行几何纠正，使得所有的遥感影像具有相同的投影信息；

采用人工目视解译的方式获取待评价区域各个时相的建设用地信息，

所述夜间灯光影像来自于DMSP/OLS的数据，对夜间灯光影像的修正过程包括：

从DMSP/OLS的数据中提取待评价区域内的夜间灯光数据；

对夜间灯光影像进行融合去噪和年际序列修正。