1.一种用于确定污染边界的污染地块采样布点加密方法,其特征在于,其包含步骤如下:S1,获取初始样本点集合;初始样本点集合从初始样本数据库中获得,初始样本数据库至少包括样本编号、空间位置和污染物浓度值三个字段;S2,将待调查工业场地生成待调查工业场地网格图;将待调查工业场地划分为规则的网格,待调查工业场地网格图的分辨率为每一个网格为一个像素点;S3,生成预测结果图;对包含有 个初始样本点的初始样本点集合进行 次有放回抽样,每次抽取的样本量为 ,其中 、 、 取值都为正整数;采用IDW空间插值方法得到 幅预测结果图,再根据 幅预测结果图得到均值预测结果图;S4,根据预测结果图得到像素值标准差;根据 幅预测结果图,计算每个像素点的像素值标准差;S5,根据均值预测结果图、像素值标准差与阈值确定待定区;根据均值预测结果图、像素值标准差与阈值将整个待调查工业场地划分为3个区域,划分标准如下: (2);其中, 为待调查工业场地网格图中的第 个像素点; 为第 个像素点的分类,1为污染区,2为清洁区,3为待定区; 为预设的污染物浓度阈值; 为在显著性水平 下的标准正态分布表 值, 表示第 个像素点的像素值标准差; 是均值预测结果图中第 个像素点的像素值;S6,确定待定区内每个网格的第一权重;针对待定区内的每个网格,计算每个网格的第一权重;S7,对待定区增加样本点;根据待定区面积的大小,设待增加的样本点总数不超过 ,分 次增加样本点,即 ,第 次增加的样本点数为 ,进行第 次加密布点选择;S8,判断增加样本点的次数是否达到 次,如果未达到就返回S3,进行第 次的布点加密,否则就结束布点加密;S9,根据布点加密完成后的样本点集合,以污染物浓度值作为待调查工业场地网格图的像素值,得到待调查工业场地污染区与清洁区分界线。
2.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S1中初始样本数据库,具体为:样本编号为每个初始样本点的ID,具有唯一性;初始样本点的空间位置为初始样本点的经纬度或者与待调查工业场地的相对位置;污染物浓度值是根据初始样本点的空间位置经实地采样后测试获得。
3.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S3中, , 。
4.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S3中,采用IDW空间插值方法得到 幅预测结果图,再根据 幅预测结果图得到均值预测结果图,具体为:对每一次抽样,根据所抽取初始样本的空间位置,利用IDW空间插值方法将待调查工业场地网格图进行污染物浓度插值生成一幅预测结果图,预测结果图中每一个像素点的像素值为污染物浓度值;再得到每一个像素点的污染物浓度均值 ,作为像素点的像素值,得到一幅均值预测结果图。
5.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S4中,根据 幅预测结果图,计算每个像素点的像素值标准差具体为: (1);其中, 为待调查工业场地网格图中的第 个像素点, 表示第 个像素点的像素值标准差, 表示第 个像素点在第 幅预测结果图中的像素值, 是均值预测结果图中第 个像素点的像素值。
6.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S6确定待定区内每个网格的第一权重,具体为:针对待定区内的每个网格,计算每个网格的第一权重: (3);其中, 第 个像素点的第一权重; 是均值预测结果图中第 个像素点的像素值; 是积分变量; 为预设的污染物浓度阈值; 是第 个像素点的像素值标准差。
7.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S7中 ,分配比例为 , ,其中 。
8.根据权利要求1所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S7中进行第 次加密布点选择,具体过程如下:S71,已确定样本点集合 ,备选样本点集合 ,加密布点集合 ,加密布点集合 初始值为空;S72,计算 中每个备选样本点到 中样本点的最近距离 作为第二权重,将每个备选点的第一权重和第二权重相乘得到各个备选样本点的最终权重,并根据最终权重找出最大值所对应的在 中的备选样本点 ;将 作为已选点添加到集合 和集合 中,并从集合 中移除;S73,判断第 次增加的样本点数是否达到 ,如果未达到就返回S72,否则执行S74;S74,获取样本点空间位置和污染物浓度值;S75,加入初始样本点;将集合 中的样本点加入初始样本点集合,也作为初始样本点。
9.根据权利要求8所述的用于确定污染边界的污染地块采样布点加密方法,其特征在于:所述S74中获取样本点空间位置和污染物浓度值,具体方式为:根据集合 中的样本点在待调查工业场地网格图中的位置,确定样本点的空间位置,待调查工业场地进行实地采样和测试后获得污染物浓度值。