CN111915669B - 一种基于总量控制的土地调查线状地物图斑化方法 - Google Patents
一种基于总量控制的土地调查线状地物图斑化方法 Download PDFInfo
- Publication number
- CN111915669B CN111915669B CN202010770183.0A CN202010770183A CN111915669B CN 111915669 B CN111915669 B CN 111915669B CN 202010770183 A CN202010770183 A CN 202010770183A CN 111915669 B CN111915669 B CN 111915669B
- Authority
- CN
- China
- Prior art keywords
- current
- area
- linear
- ground object
- linear ground
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 77
- 238000013507 mapping Methods 0.000 claims abstract description 19
- 239000002689 soil Substances 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000011835 investigation Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于总量控制的土地调查线状地物图斑化方法，该方法包括：获取目标图像中待处理线状地物的当前面积以及当前宽度，根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；本发明通过当前调整系数确定当前缓冲距离，并根据当前缓冲距离对所述待处理线状地物进行图斑化处理，从而能够使图斑化处理后的线状地物面积更接近原有线状地物总面积。
Description
技术领域
本发明涉及数据处理技术领域，尤其涉及一种基于总量控制的土地调查线状地物图斑化方法。
背景技术
目前，线状地物图斑化一般采用缓冲分析的方法进行，由于线状地物缓冲区分析后存在空间重叠的情况，为了保证数据的拓扑正确性，需要去掉重叠部分，从而导致线状地物的总面积与原来相比减少，影响地类面积汇总的正确性。
上述内容仅用于辅助理解本发明的技术方案，并不代表承认上述内容是现有技术。
发明内容
本发明的主要目的在于提供一种基于总量控制的土地调查线状地物图斑化方法，旨在解决如何使线状地物图斑化处理后的面积更接近原有线状地物总面积的技术问题。
为实现上述目的，本发明提供一种基于总量控制的土地调查线状地物图斑化方法，所述基于总量控制的土地调查线状地物图斑化方法包括以下步骤：
获取目标图像中待处理线状地物的当前面积以及当前宽度；
根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离；
根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果。
优选地，所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤之后，所述基于总量控制的土地调查线状地物图斑化方法还包括：
根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值；
判断所述当前面积差值是否大于预设阈值；
在所述当前面积差值小于或等于所述预设阈值时，保存所述目标处理结果。
优选地，所述判断所述当前面积差值是否大于预设阈值的步骤之后，所述基于总量控制的土地调查线状地物图斑化方法包括：
在所述当前面积差值大于所述预设阈值时，返回所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至所述当前面积差值小于或等于所述预设阈值。
优选地，所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，具体包括：
根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区；
获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物；
根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果。
优选地，所述根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区的步骤，具体包括：
获取所述待处理线状地物的当前中心线；
根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度；
根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区。
优选地，所述根据所述待处理线状地物以及所述重叠线状地物生成目标处理结果的步骤，具体包括：
查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果；
根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
优选地，所述获取目标图像中待处理线状地物的当前面积以及当前宽度的步骤之前，所述基于总量控制的土地调查线状地物图斑化方法还包括：
获取目标图像中待处理线状地物的初始宽度，并根据所述初始宽度生成初始缓冲距离；
根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果；
根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积；
根据所述初始线状地物面积以及预设控制面积计算初始面积差值，并判断所述初始面积差值是否大于预设阈值；
相应地，所述获取目标图像中待处理线状地物的当前面积以及当前宽度，包括：
在所述初始面积差值大于所述预设阈值时，获取目标图像中待处理线状地物的当前面积以及当前宽度。
此外，为实现上述目的，本发明还提出一种基于总量控制的土地调查线状地物图斑化装置，所述基于总量控制的土地调查线状地物图斑化装置包括：获取模块、确定模块和处理模块；
所述获取模块，用于获取目标图像中待处理线状地物的当前面积以及当前宽度；
所述确定模块，用于根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离；
所述处理模块，用于根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果。
本发明中，获取目标图像中待处理线状地物的当前面积以及当前宽度，根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；本发明通过当前调整系数确定当前缓冲距离，并根据当前缓冲距离对所述待处理线状地物进行图斑化处理，从而能够使图斑化处理后的线状地物面积更接近原有线状地物总面积。
附图说明
图1为本发明基于总量控制的土地调查线状地物图斑化方法第一实施例的流程示意图；
图2为本发明基于总量控制的土地调查线状地物图斑化方法第二实施例的流程示意图；
图3为本发明基于总量控制的土地调查线状地物图斑化方法第三实施例的流程示意图；
图4为本发明基于总量控制的土地调查线状地物图斑化装置第一实施例的结构框图。
本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参照图1，图1为本发明基于总量控制的土地调查线状地物图斑化方法第一实施例的流程示意图，提出本发明基于总量控制的土地调查线状地物图斑化方法第一实施例。
在第一实施例中，所述基于总量控制的土地调查线状地物图斑化方法包括以下步骤：
步骤S10：获取目标图像中待处理线状地物的当前面积以及当前宽度。
应当理解的是，本实施例的执行主体是所述基于总量控制的土地调查线状地物图斑化设备，其中，所述基于总量控制的土地调查线状地物图斑化设备可为个人电脑或服务器等电子设备。
需要说明的是，目标图像可以是第二次全国国土调查图像。由于，第二次全国国土调查图像的土地利用要素为地类图斑、线状地物和零星地物；而第三次全国国土调查图像的土地利用要素只有地类图斑图层，取消了线状地物和零星地物图层，线状地物和零星地物全部按照实际边界调查为图斑。因此，需要将第二次全国国土调查图像作为目标图像进行后续处理。
线状地物可以是目标图像中的河流、铁路、公路、管道用地、农村道路、林带、沟渠和田坎等；当前面积可以是待处理线状地物在当前时刻的面积数据；当前宽度可以是待处理线状地物在当前时刻的宽度数据。
步骤S20：根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离。
需要说明的是，当前缓冲距离可以是以待处理线状地物的当前中心线为扩展中心，对待处理线状地物的当前宽度的外扩距离。
可以理解的是，根据所述当前面积确定当前调整系数可以是根据当前面积通过预设调整系数计算公式计算当前调整系数，其中，所述预设调整系数计算公式如下所示：
式中，XS为当前调整系数，S为所述待处理线状地物的实际面积，S1为当前面积。
应当理解的是，根据所述当前调整系数以及所述当前宽度确定当前缓冲距离可以是根据当前调整系数以及当前宽度通过预设当前缓冲距离计算公式计算当前缓冲距离，其中，所述预设当前缓冲距离计算公式如下所示：
式中，d1为当前缓冲距离，KD1为当前宽度，XS为当前调整系数。
步骤S30：根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果。
需要说明的是，《第三次全国国土调查实施方案》中明确规定需要将第三次全国国土调查数据和第二次全国国土调查数据进行对比分析。而由于第二次全国国土调查数据利用数据采用点、线、面三个图层表达，且三个图层的逻辑关系比较复杂，增加了空间比对分析的难度。为了简化第二次全国国土调查数据和第三次全国国土调查数据进行叠加分析复杂度，提出将线状地物、零星地物进行空间模拟，变成近似的“图斑”，并和原有地类图斑进行融合，形成与三调类似的只采用面图层表达的第二次全国国土调查数据利用数据，即线状地物、零星地物图斑化。
应当理解的是，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果可以是根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区，获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物，查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果，根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
在第一实施例中，获取目标图像中待处理线状地物的当前面积以及当前宽度，根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；本实施例通过当前调整系数确定当前缓冲距离，并根据当前缓冲距离对所述待处理线状地物进行图斑化处理，从而能够使图斑化处理后的线状地物面积更接近原有线状地物总面积。
参照图2，图2为本发明基于总量控制的土地调查线状地物图斑化方法第二实施例的流程示意图，基于上述图1所示的第一实施例，提出本发明基于总量控制的土地调查线状地物图斑化方法的第二实施例。
在第二实施例中，所述步骤S30之后，还包括：
步骤S40：根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值。
可以理解的是，根据所述目标处理结果确定待计算线状地物面积可以是直接获取目标处理结果中待计算线状地物对应的图斑面积，并将图斑面积作为待计算线状地物面积。
步骤S50：判断所述当前面积差值是否大于预设阈值。
进一步地，所述步骤S50之后，还包括：
在所述当前面积差值大于所述预设阈值时，返回所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至所述当前面积差值小于或等于所述预设阈值。
应理解的是，在当前面积差大于预设阈值时，说明待计算线状地物面积与预设控制面积之间面积差值的绝对值大于预设阈值，也就意味着待计算线状地物面积与预设控制面积之间的当前面积差值过大，需要进一步精细化处理。因此，需要循环执行根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至当前面积差值小于或等于所述预设阈值。
步骤S60：在所述当前面积差值小于或等于所述预设阈值时，保存所述目标处理结果。
应理解的是，在当前面积差值小于或等于预设阈值时，说明待计算线状地物面积与预设控制面积之间面积差值的绝对值小于或等于预设阈值，也就意味着待计算线状地物面积与预设控制面积之间的当前面积差值满足条件，可以保存所述目标处理结果，并将所述目标处理结果发送至目标客户端。
在第二实施例，所述步骤S10之前，还包括：
步骤S01：获取目标图像中待处理线状地物的初始宽度，并根据所述初始宽度生成初始缓冲距离。
需要说明的是，初始宽度可以是待处理线状地物在第二次国土调查图像中的宽度。
应当理解的是，根据所述初始宽度生成初始缓冲距离可以是根据初始宽度通过预设缓冲距离公式计算初始缓冲距离，其中，预设初始缓冲距离公式如下所示：
式中，d为初始缓冲距离，KD为初始宽度，k为预设缓冲外扩参数。
预设缓冲外扩参数可以根据用户需求进行设置，本实施例中以1.007为例进行举例说明。
步骤S02：根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果。
可以理解的是，根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果可以是根据初始缓冲距离生成待处理线状地物的初始缓冲区，获取初始缓冲区的线状地物信息，并根据初始缓冲区的线状地物信息确定初始重叠线状地物，查找待处理线状地物以及初始重叠线状地物对应的地类优先级，获得初始查找结果，根据初始查找结果确定初始待删除区域，并从初始缓冲区中删除初始待删除区域，获得初始图斑化处理结果。
步骤S03：根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积。
可以理解的是，根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积可以是直接获取初始图斑化处理结果中待计算线状地物对应的初始图斑面积，并将初始图斑面积作为初始线状地物面积。
步骤S04：根据所述初始线状地物面积以及预设控制面积计算初始面积差值，并判断所述初始面积差值是否大于预设阈值。
需要说明的是，预设控制面积可以是所述待处理线状地物的实际面积；所述预设阈值可以是用户根据实际需求预先设置的参数值，也可以是根据试验获得的最佳参数值，本实施对此不加以限制。
可以理解的是，根据所述初始线状地物面积以及预设控制面积计算初始面积差值可以是计算初始线状地物面积与预设控制面积之间面积差值的绝对值，并将面积差值的绝对值作为初始面积差值。
应当理解的是，判断所述初始面积差值是否大于预设阈值，在所述初始面积差值小于或等于所述预设阈值时，说明初始线状地物面积与预设控制面积之间面积差值的绝对值小于或等于预设阈值，也就意味着初始线状地物面积与预设控制面积之间的初始面积差值满足条件，可以直接将初始图斑化处理结果作为目标处理结果。
相应地，所述步骤S10，包括：
步骤S10'：在所述初始面积差值大于所述预设阈值时，获取目标图像中待处理线状地物的当前面积以及当前宽度。
应当理解的是，在初始面积差值大于预设阈值时，说明初始线状地物面积与预设控制面积之间面积差值的绝对值大于预设阈值，也就意味着初始线状地物面积与预设控制面积之间的初始面积差值过大，需要进一步精细化处理。因此，需要返回到获取目标图像中待处理线状地物的当前面积以及当前宽度的步骤。
在第二实施例中，根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值，判断所述当前面积差值是否大于预设阈值，在所述当前面积差值小于或等于所述预设阈值时，保存所述目标处理结果；本实施例通过当前面积差值判断是否保存目标处理结果，从而能够将满足精确度条件的目标处理结果进行保存；
在第二实施例中，根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值，判断所述当前面积差值是否大于预设阈值，在所述当前面积差值大于所述预设阈值时，返回所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至所述当前面积差值小于或等于所述预设阈值；本实施例通过在当前面积差值大于预设阈值时，不断调整待处理线状地物的宽度，直至当前面积差值满足预设条件，从而能够得到图斑化处理后的线状地物面积与原有线状地物总面积一致的图斑化成果；
在第二实施例中，获取目标图像中待处理线状地物的初始宽度，并根据所述初始宽度生成初始缓冲距离，根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果，根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积，根据所述初始线状地物面积以及预设控制面积计算初始面积差值，并判断所述初始面积差值是否大于预设阈值，在所述初始面积差值大于所述预设阈值时，获取目标图像中待处理线状地物的当前面积以及当前宽度，从而能够快速确定当前调整系数。
参照图3，图3为本发明基于总量控制的土地调查线状地物图斑化方法第三实施例的流程示意图，基于上述图1所示的第一实施例，提出本发明基于总量控制的土地调查线状地物图斑化方法的第三实施例。
在第三实施例中，所述步骤S30，包括：
步骤S301：根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区。
应当理解的是，根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区可以是获取所述待处理线状地物的当前中心线，根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度，根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区。
进一步地，所述步骤S301，包括：
获取所述待处理线状地物的当前中心线；
根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度；
根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区。
需要说明的是，当前中心线可以是待处理线状地物的几何中心。
应当理解的是，根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度可以是将所述当前宽度扩展当前缓冲距离，获得目标宽度。
步骤S302：获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物。
需要说明的是，线状地物信息可以是线状地物的边界信息等。
应当理解的是，根据所述线状地物信息确定重叠线状地物可以是根据边界信息确定重叠线状地物。
步骤S303：根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果。
可以理解的是，根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果可以是查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果，根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
进一步地，所述步骤S303，包括：
查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果；
根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
可以理解的是，查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级可以是在预设映射关系表中分别查找待处理线状地物对应的第一优先级以及重叠线状地物对应的第二优先级，并将第一优先级以及第二优先级作为查找结果，所述预设映射关系表中包含线状地物与优先级之间的对应关系，所述对应关系可以是用户根据实际需求预先设置。
应当理解的是，根据所述查找结果确定待删除区域可以是比较第一优先级和第二优先级，将优先级低的线状地物对应的重叠区域作为待删除区域。
在第三实施例中，根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区，获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物，根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果；本实施例通过获取当前缓冲区的线状地物信息确定重叠线状地物，再根据待处理线状地物以及重叠线状地物确定目标处理结果，从而能够快速生成目标处理结果；
在第三实施例中，获取所述待处理线状地物的当前中心线，根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度，根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区；本实施例通过当前缓冲距离确定目标宽度，再根据当前中心线以及目标宽度生成所述待处理线状地物的当前缓冲区，从而能够准确生成待处理线状地物的当前缓冲区；
在第三实施例中，通过查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果，根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果，从而能够删除重叠面积，提高准确性。
此外，参照图4，本发明实施例还提出一种基于总量控制的土地调查线状地物图斑化装置，所述基于总量控制的土地调查线状地物图斑化装置包括：获取模块10、确定模块20和处理模块30；
所述获取模块10，用于获取目标图像中待处理线状地物的当前面积以及当前宽度。
需要说明的是，目标图像可以是第二次全国国土调查图像。由于，第二次全国国土调查图像的土地利用要素为地类图斑、线状地物和零星地物；而第三次全国国土调查图像的土地利用要素只有地类图斑图层，取消了线状地物和零星地物图层，线状地物和零星地物全部按照实际边界调查为图斑。因此，需要将第二次全国国土调查图像作为目标图像进行后续处理。
线状地物可以是目标图像中的河流、铁路、公路、管道用地、农村道路、林带、沟渠和田坎等；当前面积可以是待处理线状地物在当前时刻的面积数据；当前宽度可以是待处理线状地物在当前时刻的宽度数据。
所述确定模块20，用于根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离。
需要说明的是，当前缓冲距离可以是以待处理线状地物的当前中心线为扩展中心，对待处理线状地物的当前宽度的外扩距离。
可以理解的是，根据所述当前面积确定当前调整系数可以是根据当前面积通过预设调整系数计算公式计算当前调整系数，其中，所述预设调整系数计算公式如下所示：
式中，XS为当前调整系数，S为所述待处理线状地物的实际面积，S1为当前面积。
应当理解的是，根据所述当前调整系数以及所述当前宽度确定当前缓冲距离可以是根据当前调整系数以及当前宽度通过预设当前缓冲距离计算公式计算当前缓冲距离，其中，所述预设当前缓冲距离计算公式如下所示：
式中，d1为当前缓冲距离，KD1为当前宽度，XS为当前调整系数。
所述处理模块30，用于根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果。
需要说明的是，《第三次全国国土调查实施方案》中明确规定需要将第三次全国国土调查数据和第二次全国国土调查数据进行对比分析。而由于第二次全国国土调查数据利用数据采用点、线、面三个图层表达，且三个图层的逻辑关系比较复杂，增加了空间比对分析的难度。为了简化第二次全国国土调查数据和第三次全国国土调查数据进行叠加分析复杂度，提出将线状地物、零星地物进行空间模拟，变成近似的“图斑”，并和原有地类图斑进行融合，形成与三调类似的只采用面图层表达的第二次全国国土调查数据利用数据，即线状地物、零星地物图斑化。
应当理解的是，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果可以是根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区，获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物，查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果，根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
在本实施例中，获取目标图像中待处理线状地物的当前面积以及当前宽度，根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离，根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；本实施例通过当前调整系数确定当前缓冲距离，并根据当前缓冲距离对所述待处理线状地物进行图斑化处理，从而能够使图斑化处理后的线状地物面积更接近原有线状地物总面积。
在一实施例中，所述基于总量控制的土地调查线状地物图斑化装置还包括：判断模块：
所述判断模块，用于根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值，判断所述当前面积差值是否大于预设阈值，在所述当前面积差值小于或等于所述预设阈值时，保存所述目标处理结果；
在一实施例中，所述基于总量控制的土地调查线状地物图斑化装置还包括：循环模块：
所述循环模块，用于在所述当前面积差值大于所述预设阈值时，返回所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至所述当前面积差值小于或等于所述预设阈值；
在一实施例中，所述处理模块30，还用于根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区，获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物，根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果；
在一实施例中，所述处理模块30，还用于获取所述待处理线状地物的当前中心线，根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度，根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区；
在一实施例中，所述处理模块30，还用于查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果，根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果；
在一实施例中，所述基于总量控制的土地调查线状地物图斑化装置还包括：预处理模块；
所述预处理模块，用于获取目标图像中待处理线状地物的初始宽度，并根据所述初始宽度生成初始缓冲距离，根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果，根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积，根据所述初始线状地物面积以及预设控制面积计算初始面积差值，并判断所述初始面积差值是否大于预设阈值；
相应地，所述获取模块10，还用于在所述初始面积差值大于所述预设阈值时，获取目标图像中待处理线状地物的当前面积以及当前宽度。
本发明所述基于总量控制的土地调查线状地物图斑化装置的其他实施例或具体实现方式可参照上述各方法实施例，此处不再赘述。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。词语第一、第二、以及第三等的使用不表示任何顺序，可将这些词语解释为名称。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如只读存储器镜像(Read Only Memory image，ROM)/随机存取存储器(Random AccessMemory，RAM)、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，空调器，或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (8)
1.一种基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述基于总量控制的土地调查线状地物图斑化方法包括以下步骤：
获取目标图像中待处理线状地物的当前面积以及当前宽度；
根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离；
根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；
其中，所述根据所述当前面积确定当前调整系数包括：
根据当前面积通过预设调整系数计算公式计算当前调整系数，所述预设调整系数计算公式为：
式中，XS为当前调整系数，S为所述待处理线状地物的实际面积，S1为当前面积；
所述根据所述当前调整系数以及所述当前宽度确定当前缓冲距离，包括：
根据当前调整系数以及当前宽度通过预设当前缓冲距离计算公式计算当前缓冲距离，所述预设当前缓冲距离计算公式为：
式中，d1为当前缓冲距离，KD1为当前宽度，XS为当前调整系数。
2.如权利要求1所述的基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤之后，所述基于总量控制的土地调查线状地物图斑化方法还包括：
根据所述目标处理结果确定待计算线状地物面积，并根据所述待计算线状地物面积以及预设控制面积确定当前面积差值；
判断所述当前面积差值是否大于预设阈值；
在所述当前面积差值小于或等于所述预设阈值时，保存所述目标处理结果。
3.如权利要求2所述的基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述判断所述当前面积差值是否大于预设阈值的步骤之后，所述基于总量控制的土地调查线状地物图斑化方法包括：
在所述当前面积差值大于所述预设阈值时，返回所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，直至所述当前面积差值小于或等于所述预设阈值。
4.如权利要求1-3中任一项所述基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果的步骤，具体包括：
根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区；
获取所述当前缓冲区的线状地物信息，并根据所述线状地物信息确定重叠线状地物；
根据所述待处理线状地物以及所述重叠线状地物确定目标处理结果。
5.如权利要求4所述的基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述根据所述当前缓冲距离生成所述待处理线状地物的当前缓冲区的步骤，具体包括：
获取所述待处理线状地物的当前中心线；
根据所述当前缓冲距离对所述当前宽度进行调整，获得目标宽度；
根据所述当前中心线以及所述目标宽度生成所述待处理线状地物的当前缓冲区。
6.如权利要求4所述的基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述根据所述待处理线状地物以及所述重叠线状地物生成目标处理结果的步骤，具体包括：
查找所述待处理线状地物以及所述重叠线状地物对应的地类优先级，获得查找结果；
根据所述查找结果确定待删除区域，并从所述当前缓冲区删除所述待删除区域，获得目标处理结果。
7.如权利要求1所述的基于总量控制的土地调查线状地物图斑化方法，其特征在于，所述获取目标图像中待处理线状地物的当前面积以及当前宽度的步骤之前，所述基于总量控制的土地调查线状地物图斑化方法还包括：
获取目标图像中待处理线状地物的初始宽度，并根据所述初始宽度生成初始缓冲距离；
根据所述初始缓冲距离对所述待处理线状地物进行图斑化处理，获得初始图斑化处理结果；
根据所述初始图斑化处理结果确定所述待处理线状地物的初始线状地物面积；
根据所述初始线状地物面积以及预设控制面积计算初始面积差值，并判断所述初始面积差值是否大于预设阈值；
相应地，所述获取目标图像中待处理线状地物的当前面积以及当前宽度的步骤，具体包括：
在所述初始面积差值大于所述预设阈值时，获取目标图像中待处理线状地物的当前面积以及当前宽度。
8.一种基于总量控制的土地调查线状地物图斑化装置，其特征在于，所述基于总量控制的土地调查线状地物图斑化装置包括：获取模块、确定模块和处理模块；
所述获取模块，用于获取目标图像中待处理线状地物的当前面积以及当前宽度；
所述确定模块，用于根据所述当前面积确定当前调整系数，并根据所述当前调整系数以及所述当前宽度确定当前缓冲距离；
所述处理模块，用于根据所述当前缓冲距离对所述待处理线状地物进行图斑化处理，获得目标处理结果；
其中，所述确定模块，还用于根据当前面积通过预设调整系数计算公式计算当前调整系数，所述预设调整系数计算公式为：
式中，XS为当前调整系数，S为所述待处理线状地物的实际面积，S1为当前面积；
所述确定模块，还用于根据当前调整系数以及当前宽度通过预设当前缓冲距离计算公式计算当前缓冲距离，所述预设当前缓冲距离计算公式为：
式中，d1为当前缓冲距离，KD1为当前宽度，XS为当前调整系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010770183.0A CN111915669B (zh) | 2020-08-03 | 2020-08-03 | 一种基于总量控制的土地调查线状地物图斑化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010770183.0A CN111915669B (zh) | 2020-08-03 | 2020-08-03 | 一种基于总量控制的土地调查线状地物图斑化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111915669A CN111915669A (zh) | 2020-11-10 |
| CN111915669B true CN111915669B (zh) | 2024-04-05 |
Family
ID=73287107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010770183.0A Active CN111915669B (zh) | 2020-08-03 | 2020-08-03 | 一种基于总量控制的土地调查线状地物图斑化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111915669B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113127485B (zh) * | 2021-03-24 | 2023-02-03 | 中国电力工程顾问集团西南电力设计院有限公司 | 一种输电线路工程测量软件中线状地物的处理方法 |
| CN113688198B (zh) * | 2021-08-27 | 2024-03-05 | 湖南省国土资源规划院 | 城区范围划定中线状地物截断处理方法、存储器及装置 |
| CN117274357A (zh) * | 2023-09-28 | 2023-12-22 | 中国国土勘测规划院 | 线状地物的图斑化方法、装置、电子设备及可读存储介质 |
| CN117853764B (zh) * | 2023-12-25 | 2024-07-23 | 广东国地规划科技股份有限公司 | 一种基于区域生长约束的图斑筛选方法及系统 |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102298705A (zh) * | 2011-08-12 | 2011-12-28 | 北京师范大学 | 一种分析景观特征对遥感分类图斑精度影响的方法 |
| CN102902747A (zh) * | 2012-09-19 | 2013-01-30 | 浙江大学 | 基于空间聚合的土地利用数据更新方法 |
| CN106530170A (zh) * | 2016-09-21 | 2017-03-22 | 环境保护部南京环境科学研究所 | 一种区域生态保护红线的边界优化方法 |
| CN106846332A (zh) * | 2016-12-30 | 2017-06-13 | 中国人民解放军61540部队 | 一种基于dsm的遥感影像变化检测方法及装置 |
| CN108132974A (zh) * | 2017-12-07 | 2018-06-08 | 陈菡 | 用于将线状地物转为图斑的方法和装置 |
| WO2018107939A1 (zh) * | 2016-12-14 | 2018-06-21 | 国家海洋局第二海洋研究所 | 一种基于边缘完备度的图像分割最优识别方法 |
| CN109410227A (zh) * | 2018-05-15 | 2019-03-01 | 刘宇轩 | 一种基于gvf模型的土地利用图斑轮廓提取算法 |
| CN109658509A (zh) * | 2018-12-11 | 2019-04-19 | 中国测绘科学研究院 | 一种狭长图斑分块方法、装置、计算机设备和存储介质 |
| CN110162650A (zh) * | 2019-04-18 | 2019-08-23 | 中国测绘科学研究院 | 一种兼顾局部最优与整体面积平衡的小图斑融解方法 |
| CN110176018A (zh) * | 2019-04-18 | 2019-08-27 | 中国测绘科学研究院 | 一种保持结构化地物轮廓特征的图斑合并方法 |
| CN110688756A (zh) * | 2019-09-28 | 2020-01-14 | 滁州学院 | 一种采集和制图一体化的竣工图制作方法 |
| CN110955663A (zh) * | 2019-12-02 | 2020-04-03 | 重庆市勘测院 | 一种大尺度区域土地资源资产负债表编制方法 |
| CN111145157A (zh) * | 2019-12-27 | 2020-05-12 | 国交空间信息技术（北京）有限公司 | 基于高分辨率遥感影像的路网数据自动化质量检查方法 |
-
2020
- 2020-08-03 CN CN202010770183.0A patent/CN111915669B/zh active Active
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102298705A (zh) * | 2011-08-12 | 2011-12-28 | 北京师范大学 | 一种分析景观特征对遥感分类图斑精度影响的方法 |
| CN102902747A (zh) * | 2012-09-19 | 2013-01-30 | 浙江大学 | 基于空间聚合的土地利用数据更新方法 |
| CN106530170A (zh) * | 2016-09-21 | 2017-03-22 | 环境保护部南京环境科学研究所 | 一种区域生态保护红线的边界优化方法 |
| WO2018107939A1 (zh) * | 2016-12-14 | 2018-06-21 | 国家海洋局第二海洋研究所 | 一种基于边缘完备度的图像分割最优识别方法 |
| CN106846332A (zh) * | 2016-12-30 | 2017-06-13 | 中国人民解放军61540部队 | 一种基于dsm的遥感影像变化检测方法及装置 |
| CN108132974A (zh) * | 2017-12-07 | 2018-06-08 | 陈菡 | 用于将线状地物转为图斑的方法和装置 |
| CN109410227A (zh) * | 2018-05-15 | 2019-03-01 | 刘宇轩 | 一种基于gvf模型的土地利用图斑轮廓提取算法 |
| CN109658509A (zh) * | 2018-12-11 | 2019-04-19 | 中国测绘科学研究院 | 一种狭长图斑分块方法、装置、计算机设备和存储介质 |
| CN110162650A (zh) * | 2019-04-18 | 2019-08-23 | 中国测绘科学研究院 | 一种兼顾局部最优与整体面积平衡的小图斑融解方法 |
| CN110176018A (zh) * | 2019-04-18 | 2019-08-27 | 中国测绘科学研究院 | 一种保持结构化地物轮廓特征的图斑合并方法 |
| CN110688756A (zh) * | 2019-09-28 | 2020-01-14 | 滁州学院 | 一种采集和制图一体化的竣工图制作方法 |
| CN110955663A (zh) * | 2019-12-02 | 2020-04-03 | 重庆市勘测院 | 一种大尺度区域土地资源资产负债表编制方法 |
| CN111145157A (zh) * | 2019-12-27 | 2020-05-12 | 国交空间信息技术（北京）有限公司 | 基于高分辨率遥感影像的路网数据自动化质量检查方法 |
Non-Patent Citations (5)
| Title |
|---|
| Area discription method of land use features based on sets algebra;Gao Yunbing等;《2016 Fifth International Conference on Agro-Geoinformatics》;1-6 * |
| Detecting peatland drains with pbject based image analysis and Geoeye-1 imagery;J. Connolly等;《Carbon Balance and Management》;第07卷(第2017期);1-13 * |
| 乡级土地规划数据库质量控制关键技术探讨;马黎莉等;《测绘科学》;第37卷(第03期);191-193 * |
| 土地利用数据多尺度表达模型构建及实现;陈青;《中国优秀硕士学位论文全文数据库 (经济与管理科学辑)》(第2017(07)期);J149-21 * |
| 第三次全国国土调查线状地物图斑化关键技术研究——以和林格尔县为例;刘志涛;《中国优秀硕士学位论文全文数据库 (信息科技辑)》(第2019（08）期);I138-608 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111915669A (zh) | 2020-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111915669B (zh) | 一种基于总量控制的土地调查线状地物图斑化方法 | |
| CN110443786B (zh) | 激光雷达点云滤波方法、装置、计算机设备和存储介质 | |
| US20100290708A1 (en) | Image retrieval apparatus, control method for the same, and storage medium | |
| CN112432647B (zh) | 车厢的定位方法、装置、系统及计算机可读存储介质 | |
| CN111626348B (zh) | 自动泊车测试模型构建方法、设备、存储介质及装置 | |
| CN117195361B (zh) | 高边坡稳定性监测方法、装置、设备及可读存储介质 | |
| CN110097108B (zh) | 非机动车的识别方法、装置、设备及存储介质 | |
| CN111079730A (zh) | 一种在界面图中确定样图所在区域的方法和电子设备 | |
| CN117994250B (zh) | 晶圆缺陷检测方法、装置、设备及存储介质 | |
| CN113436223A (zh) | 点云数据的分割方法、装置、计算机设备和存储介质 | |
| CN116385745A (zh) | 图像识别方法、装置、电子设备和存储介质 | |
| CN113808142B (zh) | 一种地面标识的识别方法、装置、电子设备 | |
| CN113012128B (zh) | 缺陷表征方法和装置 | |
| CN111340788B (zh) | 硬件木马版图检测方法、装置、电子设备和可读存储介质 | |
| CN113449021B (zh) | 空间坐标转换方法、装置、计算机设备和可读存储介质 | |
| CN111626288B (zh) | 数据处理方法、装置、计算机设备和存储介质 | |
| CN111060127B (zh) | 车辆起点定位方法、装置、计算机设备和存储介质 | |
| US20210405197A1 (en) | GLOBAL LOCALIZATION APPARATUS AND METHOD IN DYNAMIC ENVIRONMENTS USING 3D LiDAR SCANNER | |
| CN114821173A (zh) | 图像分类方法、装置、设备及存储介质 | |
| KR102350111B1 (ko) | 사용자의 특정 입력 데이터로부터 추출되는 적어도 하나의 특정 식별어구 및 이에 대응되는 적어도 하나의 도킹결과데이터에 접근할 수 있는 적어도 하나의 위치지정자를 포함한 복합웹페이지를 발행하는 방법 및 이를 이용한 복합웹페이지발행서버 | |
| CN115830248B (zh) | 一种基于无人机环绕规划的超高精细三维模型建立方法 | |
| CN112492010B (zh) | 基于定位信息的消息推送系统、方法、装置及存储介质 | |
| CN116663492B (zh) | 交叉仓所覆盖的交叉项数量的获取方法、电子设备和介质 | |
| JP4950970B2 (ja) | 画像特徴抽出装置 | |
| US12067704B2 (en) | Defect characterization method and apparatus |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |