基于循环相关的覆冰微地形识别方法及系统与流程

本发明涉及电网防护
技术领域：
，尤其涉及一种基于循环相关的覆冰微地形识别方法及系统。
背景技术：
：随着全球气候变化，近年来冬季冰灾频频发生。特别是微地形区域，覆冰十分严重。由于覆冰仅发生在几百米范围内，与大尺度覆冰迥然不同，且线路覆冰增长快，极难预测。为了更好提升预测微地形覆冰，首先需要定位覆冰微地形。目前，常用的微地形识别方法主要依据人工经验，根据现场实地的勘察，确定地形类型。这种方式需要消耗大量的人力物力，工作强度大，对全国区域进行这种方式识别，十分困难。同时，识别过程中人为主观性强，导致地形分类判识结果存在很大的不确定性。因此需要建立一种客观定量的判识方法，能够较好的实现全国范围的微地形判识计算。技术实现要素：本发明提供了一种基于循环相关的覆冰微地形识别方法及系统，用以解决现有用的微地形人工识别方法，消耗人力物力多，工作强度大，对大范围区域识别时十分困难，并且人为主观性强，难以统一标准的技术问题。为解决上述技术问题，本发明提出的技术方案为：一种基于循环相关的覆冰微地形识别方法，包括以下步骤：收集历史上同一类类型的微地形坐标，并组成该类微地形的地形数据库；在地形数据库中，通过网格划分，计算每个网格点的平均海拔，根据平均海拔，计算得到基准微地形特征参数；计算目标点的地形特征参数，并与基准微地形特征进行相关分析；根据相关分析确定目标区域是否属于微地形。优选地，根据平均海拔，计算得到基准微地形特征参数，包括：对微地形点为中心网格及其周围8个网格点组成的9个网格点按照计算的平均海拔从大到小排序，得到中心点海拔大小所处的排序位置ln；对该类微地形的地形数据库中，算得中心点海拔的排序位置集合lo，排序位置范围为[lmin，lmax]；将平均海拔，两两匹配一组{p,q}，固定p的顺序，对q进行循环排序，分别计算每一组的相关，得到8个相关系数；选取每一组{p,q}中最大的值组成集合，从而得到所有最大相关的最小值t。优选地，计算目标点的地形特征参数，并与基准微地形特征进行相关分析，包括：选取目标点t，划分网格点，对以目标点为中心网格点及其周围8个网格点组成的9个网格点分别计算平均海拔，按照计算的平均海拔对各网格点从大到小排序，得到目标点的排序位置tl；计算以目标点为中心网格点及其周围8个网格点组成的目标区域与地形数据库中每一个微地形的8个循环相关系数，并记录其中最大相关为tmax。优选地，根据相关分析确定目标区域是否属于微地形，包括：ⅰ、当tmax≥t，进行中心点判断：当tl∈[lmin，lmax]，则目标点t属于该类微地形；当则目标点t不属于该类微地形；ⅱ、当tmax＜t，则目标点t不属于该类微地形。本发明还提供一种计算机系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一方法的步骤。本发明具有以下有益效果：本发明的基于循环相关的覆冰微地形识别方法，从已有的微地形数据库中，通过网格划分，计算得到微地形特征，同时计算目标区域的地形特征参数，并与微地形特征进行相关分析，从而确定目标区域是否属于微地形。本发明计算流程简便，通过数值化微地形点，能够定量计算微地形的地形特征。本发明通用性好，能够用于不同地区的不同类型微地形特征判识。除了上面所描述的目的、特征和优点之外，本发明还有其它的目的、特征和优点。下面将参照附图，对本发明作进一步详细的说明。附图说明构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：图1是本发明优选实施例的基于循环相关的覆冰微地形识别方法的流程示意图；图2是本发明优选实施例1的特征参数示意图；图3是本发明优选实施例2的地形数据库的3组网格平均海拔高度数据图；图4是本发明优选实施例2的目标点t及周围网格点的海拔高度数据示意图；图5是本发明优选实施例3的目标点t及周围网格点的海拔高度数据示意图；图6是本发明优选实施例4的目标点t及周围网格点的海拔高度数据示意图。具体实施方式以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。参见图1，本发明的基于循环相关的覆冰微地形识别方法，包括以下步骤：收集历史上同一类类型的微地形坐标，并组成该类微地形的地形数据库；在地形数据库中，通过网格划分，计算每个网格点的平均海拔，根据平均海拔，计算得到基准微地形特征参数；计算目标点的地形特征参数，并与基准微地形特征进行相关分析；根据相关分析确定目标区域是否属于微地形。通过上述步骤，提供了一种客观定量的判识微地形的方法，能够较好的实现全国范围的微地形判识计算。实际实施时，以上的方法还能进行以下的扩充或应用，以下实施例中的技术特征都能相互组合，实施例仅作为示例，不作为对技术特征的正常组合限制。实施例1本实施例的基于循环相关的覆冰微地形识别方法，包括以下步骤：s1：历史典型微地形收集：收集历史上同一类类型的微地形坐标，并组成该类微地形的地形库记为o，其中微地形库o中包含m个相似微地形。s2：地形网格划分：采用分辨率为s的数字高程数据(dem)，以步骤s1建立的微地形库o中的微地形点为中心，进行网格划分。得到以微地形点为中心网格及其周围8个网格点，对于微地形库o而言，共有m组网格。s3：基准微地形特征参数计算：根据步骤s1建立的地形库o和步骤s2划分的网格，计算其中每个网格点的平均海拔。参见图2，对于微地形中心点海拔记录为an0，周围8个格点分别记录为an1，an2，an3，an4，an5，an6，an7，an8，其中n＝1…m。对微地形点为中心网格及其周围8个网格点组成的9个网格点按照计算的平均海拔从大到小排序，得到中心点海拔大小所处的排序位置ln。对微地形集合o，则有中心点海拔的排序位置集合lo，排序位置范围为[lmin，lmax]。s4：目标地形参数计算：选取目标点t，采用步骤s2和步骤s3的方法，以目标点为中心向周围划分网格点并进行地形参数计算：对目标点为中心网格点及其周围8个网格点组成的9个网格点分别计算平均海拔，其中目标点海拔记录为t0，周围8个点的海拔分别记录为t1，t2，t3，t4，t5，t6，t7，t8。采用步骤s3的方法对网格点按照海拔从大到小进行排序，得到目标点(中心点)的排序位置tl。s5：循环相关计算：s5.1：特征地形内循环相关根据步骤s3计算的微地形平均海拔，两两匹配一组{p,q}，固定p的顺序，对q进行循环排序，参见表1，然后分别计算每一组的相关，得到8个相关系数。表1实施例1中的循环排序表pn1pn2pn3pn4pn5pn6pn7pn8qn1qn2qn3qn4qn5qn6qn7qn8qn2qn3qn4qn5qn6qn7qn8qn1qn3qn4qn5qn6qn7qn8qn1qn2qn4qn5qn6qn7qn8qn1qn2qn3qn5qn6qn7qn8qn1qn2qn3qn4qn6qn7qn8qn1qn2qn3qn4qn5qn7qn8qn1qn2qn3qn4qn5qn6qn8qn1qn2qn3qn4qn5qn6qn7s5.2：特征地形与目标地形循环相关：与步骤s5.1类似，计算以目标点为中心网格点及其周围8个网格点组成的目标区域与库o中每一个微地形点的8个循环相关系数，共计算m组，并记录其中最大相关为tmax。s6：相关结果筛选：根据步骤s5.1计算的结果，选取每一组{p,q}中最大的值组成集合，从而得到所有最大相关的最小值t。s7：地形相似度判识：根据步骤s5.2计算的最大相关tmax与步骤s6计算的所有最大相关的最小值t，以及步骤s3得到的排序位置范围为[lmin，lmax]进行判断。ⅰ、当tmax≥t，进行中心点判断：当tl∈[lmin，lmax]，则目标区域t属于该类微地形；当则目标区域t不属于该类微地形。ⅱ、当tmax＜t，则目标区域t不属于该类微地形。实施例2本实施例以某地区的高山微地形为例，说明本发明的基于循环相关的覆冰微地形识别方法，具体步骤如下：s1：历史典型微地形收集：收集某地区高山覆冰微地形坐标，并组成该类微地形的地形库记为o，其中微地形库o中包含3个相似微地形。s2：地形网格划分：采用分辨率为500m的数字高程数据(dem)，以步骤s1建立的微地形库o中的微地形点为中心，进行网格划分。得到以微地形点为中心网格及其周围8个网格点，对于微地形库o而言，共有3组网格。s3：特征地形参数计算根据步骤s1建立的地形库o和步骤s2划分的网格，计算其中每一个网格点的平均海拔。参见图3，对于微地形中心点海拔分别记录为a10，a20，a30，对应周围8个格点分别记录为{a11，a12，a13，a14，a15，a16，a17，a18}，{a21，a22，a23，a24，a25，a26，a27，a28}，{a31，a32，a33，a34，a35，a36，a37，a38}。对9个格点按照从大到小排序，得到中心点海拔大小所处的排序位置{1，1，2}，则有中心点海拔位置集合lo，排序位置范围为[1，2]。s4：目标地形参数计算：选取目标点t，采用步骤s2和步骤s3的方法，以t为中心向周围划分网格点并进行地形参数计算。参见图4，其中目标点海拔记录为t0，周围8个点的海拔分别记录为t1，t2，t3，t4，t5，t6，t7，t8。采用步骤s3的方法对目标点按照海拔从大到小进行排序，得到目标中心点的排序位置为tl＝1。s5：循环相关计算：s5.1：特征地形内循环相关：根据步骤s3计算的微地形平均海拔，两两匹配一组{p,q}，共匹配3组，固定一组中p的顺序，对q进行循环排序，参见表2，然后分别计算每一组相关，得到3组，每组8个相关系数。表2实施例2～4的循环排序表第一组0.650.43-0.10-0.51-0.74-0.560.220.61第二组0.07-0.13-0.32-0.34-0.400.290.320.51第三组-0.040.05-0.38-0.470.23-0.240.120.73s5.2特征地形与目标地形循环相关：与步骤s5.1类似，计算以目标点为中心网格点及其周围8个网格点组成的目标区域与库o中每一个微地形点的8个循环相关系数，共计算3组，并记录3组中最大相关系数为tmax＝0.96。s6：相关结果筛选：根据步骤s5.1计算的结果，选取每一组{p,q}中最大的值组成集合为[0.65,0.51,0.73]，从而得到最大相关的最小值为0.51。s7：地形相似度判识：根据步骤s5.2计算的最大相关tmax与步骤s6计算的所有最大相关的最小值0.51，以及步骤s3得到的排序位置范围为[1，2]，进行判断。tmax≥0.51，进行中心点判断：tl＝1，因此tl∈[1，2]，则t属于该类微地形。实施例3本实施例的基于循环相关的覆冰微地形识别方法，具体步骤如下：s1～s3：重复实施例2中的步骤s1～s3；s4：目标地形参数计算：选取目标点t，采用步骤s2和步骤s3的方法，以t为中心向周围划分网格点并进行地形参数计算。参见图5，其中目标点海拔记录为t0，周围8个点的海拔分别记录为t1，t2，t3，t4，t5，t6，t7，t8。采用步骤s3的方法对目标点按照海拔从大到小进行排序，得到目标点(中心网格点)的排序位置为tl＝5。s5～s6：执行实施例2中的步骤s5～s6；s7：地形相似度判识：根据步骤s5.2计算的最大相关tmax与步骤s6计算的所有最大相关的最小值0.51，以及步骤s3得到的排序位置范围为[1，2]，进行判断。tmax≥0.51，进行中心点判断tl＝5，因此则t不属于该类微地形。实施例4本实施例的基于循环相关的覆冰微地形识别方法，具体步骤如下：s1～s3：重复实施例2中的步骤s1～s3；s4：目标地形参数计算：选取目标点t，采用步骤s2和步骤s3的方法，以t为中心向周围划分网格点并进行地形参数计算。参见图6，其中目标点海拔记录为t0，周围8个点的海拔分别记录为t1，t2，t3，t4，t5，t6，t7，t8。采用步骤s3的方法对目标点按照海拔从大到小进行排序，得到目标点的排序位置为tl＝4。s5：循环相关计算：s5.1：执行实施例2中的步骤s5.1；s5.2：特征地形内循环相关：与步骤s5.1类似，计算以目标点为中心网格点及其周围8个网格点组成的目标区域与库o中每一个微地形点的8个循环相关系数，共计算3组，并记录3组中最大相关系数为tmax＝0.41。s6：执行实施例2中的步骤s6；s7：地形相似度判识：根据步骤s5.2计算的最大相关tmax与步骤s6计算的所有最大相关的最小值0.51，进行判断。tmax＜0.51，则t不属于该类微地形。实施例5本实施例提供一种计算机系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一实施例的步骤。综上所述，本发明的基于循环相关的覆冰微地形识别方法，从已有的微地形数据库中，通过网格划分，计算得到微地形特征，同时计算目标区域的地形特征参数，并与微地形特征进行相关分析，从而确定目标区域是否属于微地形。本发明计算流程简便，通过数值化微地形点，能够定量计算微地形的地形特征。本发明通用性好，能够用于不同地区的不同类型微地形特征判识。以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页12