WO2024125434A1 - Regional-consistency-based building principal angle correction method - Google Patents

Regional-consistency-based building principal angle correction method Download PDF

Info

Publication number
WO2024125434A1
WO2024125434A1 PCT/CN2023/137754 CN2023137754W WO2024125434A1 WO 2024125434 A1 WO2024125434 A1 WO 2024125434A1 CN 2023137754 W CN2023137754 W CN 2023137754W WO 2024125434 A1 WO2024125434 A1 WO 2024125434A1
Authority
WO
WIPO (PCT)
Prior art keywords
angle
main
building
building structure
regional
Prior art date
Application number
PCT/CN2023/137754
Other languages
French (fr)
Chinese (zh)
Inventor
刘杰
陈洋涛
董铱斐
邹圣兵
Original Assignee
北京数慧时空信息技术有限公司
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 北京数慧时空信息技术有限公司 filed Critical 北京数慧时空信息技术有限公司
Publication of WO2024125434A1 publication Critical patent/WO2024125434A1/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/13Satellite images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/26Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/44Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/761Proximity, similarity or dissimilarity measures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/762Arrangements for image or video recognition or understanding using pattern recognition or machine learning using clustering, e.g. of similar faces in social networks
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/764Arrangements for image or video recognition or understanding using pattern recognition or machine learning using classification, e.g. of video objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/176Urban or other man-made structures

Definitions

  • the invention relates to the field of remote sensing image processing, and in particular to a method for correcting a main angle of a building based on regional consistency.
  • Deep learning is a new stage in the development of machine learning in artificial intelligence, which effectively solves the problems of characterizing the features of complex objects and analyzing the association of complex scenes.
  • the deep learning method for extracting buildings from high-resolution remote sensing images can automatically extract the characteristic information of buildings and achieve high-precision and high-efficiency building extraction.
  • buildings due to the complexity of remote sensing images, buildings are affected by noise, occlusion, shadows, and low contrast.
  • the current automatic building extraction method has low reliability and obvious defects. For example, the spectral characteristics of buildings and roads are similar and cannot be effectively distinguished, and it is impossible to directly obtain high-precision extraction results through fully automatic methods.
  • Another existing method is to combine computer automatic extraction technology with manual interaction, that is, interactive object extraction, and the extraction accuracy is guaranteed to a certain extent.
  • this method requires the initial position of the building to be given manually, relies on edge information to extract buildings, has complex interactions, requires high accuracy for the manually given position, and is difficult to implement.
  • the main angle of the building structure obtained by calculation can help improve the accuracy of the building extraction by the deep learning method.
  • the existing main angle calculation method still has two difficult problems to overcome. First, due to the diverse shapes of the buildings themselves, coupled with the errors caused by automated extraction, it is difficult to obtain the main angle of the building; second, different main angle calculation methods have certain scopes of application and limitations, especially for buildings with poor extraction effects, the results obtained by different calculation methods are quite different, making it difficult to achieve good post-processing results.
  • the present invention proposes a method for correcting the main angle of a building based on regional consistency.
  • the method corrects the main angle of the irregular building according to the consistency of the orientation of the buildings in the area where the irregular building is located, thereby improving the accuracy of the main angle calculation in post-processing, thereby greatly improving the post-processing optimization effect of building extraction.
  • a method for correcting a main angle of a building based on regional consistency comprising the following steps:
  • S1 obtains a target remote sensing image, and obtains candidate main angles of buildings and structures in the target remote sensing image by a preset method, wherein the candidate main angles include a first candidate value and a second candidate value;
  • S2 calculates the difference between the first candidate value and the second candidate value, and according to the relationship between the difference and the first preset threshold, divides the building structure into a standard building structure and a building structure to be corrected, and simultaneously obtains the true main angle of the standard building structure;
  • S3 obtains the regional main angle of the building structure to be corrected through a regional statistical strategy, and the regional statistical strategy includes external expansion main angle statistics and line segment histogram statistics;
  • S4 calculates the difference between the candidate main angle and the regional main angle, and selects the real main angle of the building structure to be corrected from the candidate main angles according to the relationship between the difference and the second preset threshold.
  • step S1 comprises:
  • S11 obtains the angles of the two major axis center lines of the circumscribed rectangle of the minimum area of the building structure as the first candidate value of the building structure;
  • S12 obtains the vector contour line segment of the building structure by a vector contour extraction method, and obtains the angle of the vector contour line segment;
  • S13 groups vector contour line segments with similar angles into one group, selects a group of vector contour line segments with the longest total length as a reference vector contour line segment group, and obtains a second candidate value of the building structure by a vertical equivalent method based on the angle of the reference vector contour line segment;
  • the vector contour segments of similar angles include those with angle differences between [0°, 2°] ⁇ [88°, 90°] Vector contour segments of .
  • the vector contour extraction method comprises:
  • the vertical equivalent method includes:
  • the auxiliary vector angle is vertically transformed toward the main vector angle to obtain a vertical auxiliary vector angle, and the average of the main vector angle and the vertical auxiliary vector angle is used as the second candidate value.
  • step S2 includes:
  • the building structure is used as a standard building structure, and the second candidate value is used as the true main angle of the standard building structure;
  • the building structure is used as the building structure to be corrected.
  • step S3 includes:
  • the minimum circumscribed rectangle of the building structure to be modified is expanded by a preset distance to obtain an expanded rectangle
  • the standard housing structures within the outer expansion rectangle and intersecting with the outer expansion rectangle are used as reference housing structures;
  • the average of the real main angles of the reference building structures is used as the regional main angle of the building structure to be corrected.
  • step S3 includes:
  • the houses and structures in the target remote sensing image are partitioned by hierarchical clustering method based on distance matrix;
  • the angle histogram of each partition is constructed by using the angle of the vector contour line segment of the building structure in each partition;
  • the highest frequency pair in the angle histogram of the partition where the building to be corrected is located The corresponding angle is used as the main angle of the area of the building structure to be modified.
  • step S4 includes:
  • a candidate main angle close to the regional main angle is selected as the true main angle of the building structure to be corrected
  • the second candidate value is selected as the true main angle of the building structure to be corrected.
  • the present invention proposes a new method for correcting the main angle of a building based on regional consistency.
  • the regional main angle of the building structure to be corrected is calculated by a regional statistical strategy, and the optimal value of the candidate main angle is selected based on the regional main angle to obtain the true main angle.
  • the present invention solves the problem of difficulty in calculating the main angle of irregular buildings, thereby greatly improving the post-processing optimization effect of building extraction in remote sensing images.
  • FIG1 is a flow chart of a method for correcting a main angle of a building based on regional consistency according to the present invention
  • FIG2 is a schematic diagram of the original distribution of buildings provided in Example 1 of the present invention.
  • FIG3 is a schematic diagram of a correction result of a main angle of a building in Embodiment 1 of the present invention.
  • FIG. 4 is a schematic diagram of the original distribution of buildings provided in the third embodiment of the present invention.
  • FIG. 1 is a flow chart of an embodiment of a method for correcting a main angle of a building based on regional consistency according to the present invention.
  • the method comprises the following steps:
  • S1 obtains a candidate main angle of a building structure by a preset method, wherein the candidate main angle includes a first candidate value and a second candidate value;
  • S2 classifies the housing structures into standard housing structures and housing structures to be corrected according to the difference between the first candidate value and the second candidate value;
  • S3 obtaining the real main angle of the standard building structure within the preset range, and calculating the regional main angle of the building structure to be corrected according to the real main angle of the standard building structure within the preset range;
  • S4 selects the value of the main angle of the region closest to the building structure to be corrected from the candidate main angles of the building structure to be corrected, and obtains the real main angle of the building structure to be corrected.
  • the technical idea of the present invention is as follows: 1.
  • the current building extraction method based on remote sensing images can only extract relatively regular buildings with obvious features, and its versatility is poor. When the buildings are dense or irregular buildings are encountered, the extraction effect is often not ideal, and appropriate post-processing steps are required to optimize the accuracy of the extracted buildings. Among them, obtaining the main angle of the building is of great help to post-processing optimization.
  • the method of obtaining the main angle of the building usually includes: 1) obtaining the minimum circumscribed rectangle of the building, and the direction of the two major axis center lines of the minimum circumscribed rectangle is used as the first candidate value of the main angle of the building; 2) obtaining the contour line segment of the building, calculating the total length of all similar angle contour line segments, and the direction with the longest total length is used as the second candidate value of the main angle of the building.
  • different methods for calculating the main angle may also obtain different results. How to further determine the most realistic main angle is emphasized in the subsequent steps of the present invention; 2.
  • the present invention introduces the concept of regional main angle to help determine the real main angle of the building.
  • the regional main angle is the average value of the main angle of a group of buildings with similar main angles within a certain area where the building is located.
  • the regional main angle can reflect the regional consistency of the building area and provide guidance for the main angle of the building to be corrected that is consistent with the regional consistency; 3. Select the candidate main angle of the building to be corrected that is closest to the regional main angle as the true main angle. This method is guided by the regional consistency of the building area, can determine the true main angle of the irregular building to be corrected, and further improve the effect of building extraction post-processing, thereby achieving high-precision building extraction.
  • This embodiment provides a method for obtaining the regional main angle of the building structure to be corrected based on the statistics of the external expansion main angle. At the same time, this embodiment provides a method for obtaining the real main angle of the building structure to be corrected when the target difference is less than the second preset threshold:
  • S1 obtains a target remote sensing image, and obtains candidate main angles of buildings and structures in the target remote sensing image by a preset method, wherein the candidate main angles include a first candidate value and a second candidate value;
  • S11 obtains the angles of the two major axis center lines of the circumscribed rectangle of the minimum area of the building structure as the first candidate value of the building structure;
  • S12 obtains the vector contour line segment of the building structure by a vector contour extraction method, and obtains the angle of the vector contour line segment;
  • S13 groups vector contour segments with similar angles into a group and selects the segment with the longest total length.
  • a set of vector contour line segments is used as a reference vector contour line segment set, and based on the angles of the reference vector contour line segments, a second candidate value of the building structure is obtained by a vertical equivalent method;
  • the vector contour line segments of similar angles are vector contour line segments whose angle difference is [0°, 2°] ⁇ [88°, 90°].
  • the vector contour extraction method includes:
  • This embodiment uses mask-RCNN to perform semantic segmentation on the target remote sensing image.
  • Mask-RCNN follows the idea of Faster RCNN, and uses the architecture of ResNet-FPN for feature extraction.
  • a Mask prediction branch is added, which can complete high-quality semantic segmentation tasks. It has a good extraction effect on ordinary buildings, and the effect is significantly reduced only when facing complex irregular buildings.
  • This embodiment uses ArcGIS to make samples by manual sample outlining. During the training process, pixel blocks cropped to 512 ⁇ 512 are used as samples, and about 7,000 instances are used as training data sets and about 3,000 instances are used as verification data sets.
  • Computer configuration CPU: Xeon(R)CPU E5-2620 v4@2.10GHz ⁇ 32; Graphics card: Quadro M4000; Memory: 128G; Operating system: Ubuntu 16.04. The training process took more than 40 hours and 90,000 iterations.
  • the segmented data is binarized, and 0.5 is used as the binarization threshold of the building target probability map to obtain a segmented binary image;
  • the segmentation binary image of the building structure is used to highlight the outline of the target building and compress the overall image data volume
  • the segmented binary image is post-processed by mathematical morphology method to obtain the post-processed image:
  • the mathematical morphology method consists of a set of morphological algebraic operators, and the basic operations include: dilation, erosion, opening and closing;
  • the segmented binary image is denoised by mathematical morphology methods to remove useless information such as classification noise and small non-building structures.
  • Image shape and structure analysis and processing can be performed based on mathematical morphology methods, including image segmentation, feature extraction, edge detection, image filtering, image enhancement and restoration.
  • the boundary of the post-processed image is traced by Moore's field boundary tracing method to generate the contour line segments of the building structure;
  • vector extraction is performed on the post-processed image to obtain the vector contour line segments of the building structures.
  • the vertical equivalent method includes:
  • the total length of the line segments of the main vector contour line segment group is greater than that of the auxiliary vector contour line segment group.
  • the angle difference between the main vector contour line segment group and the auxiliary vector contour line segment group is [0°, 2°], and the angle difference between the groups is [88°, 90°], so the main vector angle and the auxiliary vector angle are close to a vertical relationship;
  • xi is the main vector angle
  • yi is the auxiliary vector angle
  • f(xi) is the vertical auxiliary vector angle
  • the angle difference between the vertical auxiliary vector angle and the main vector angle is [0°, 2°];
  • the average of the main vector angle and the vertical auxiliary vector angle is used as the second candidate value b2.
  • S2 calculates the difference between the first candidate value and the second candidate value, and divides the building structures into standard building structures and building structures to be corrected according to the relationship between the difference and the first preset threshold.
  • the building structure is used as a standard building structure, and the second candidate value is used as the true main angle of the standard building structure;
  • the building structure is used as the building structure to be corrected.
  • the preset threshold is set to 3°. If the difference between the first candidate value and the second candidate value is less than 3°, it is considered that the main direction of the building structure does not need to be corrected, and the building structure is taken as a standard building structure, and the second candidate value is taken as the true main angle of the standard building structure.
  • S3 obtains the real main angle of the standard building structure within the preset range, and obtains the regional main angle of the building structure to be corrected according to the shape information of the standard building structure within the preset range.
  • the regional statistical strategy is external expansion main angle statistics
  • step S3 includes:
  • A is the building to be corrected, and other houses within a certain distance d are reference buildings.
  • the houses involved in the main angle calculation of the area include Includes all houses within a distance d from A and those intersecting the range line.
  • the real main angle of the reference building structure is clustered by the nearest neighbor clustering method:
  • the nearest neighbor clustering method is classified according to the Euclidean distance of the true main angle, and the Euclidean distance is expressed as:
  • dist(a i ,a j ) is the Euclidean distance from the true main angle a i to the true main angle a j ;
  • a reference group with the largest number of reference buildings and structures is selected, and the mean of the true main angles of the reference buildings and structures in the group is calculated to obtain the regional main angle of the building and structure to be corrected;
  • the reference group with the largest number of reference buildings and structures is A1 .
  • the mean of the main angles of the buildings and structures in this group is calculated to obtain the main angle of the area of the buildings and structures to be corrected.
  • S4 selects the value of the main angle of the region closest to the building structure to be corrected from the candidate main angles of the building structure to be corrected, and obtains the real main angle of the building structure to be corrected.
  • step S4 includes:
  • the target difference is smaller than the second preset threshold, and the candidate main angle corresponding to the target difference is selected as the true main angle of the building structure to be corrected.
  • the acquisition of the main angle of the area in this embodiment is shown in Figure 2.
  • the target difference is less than the second preset threshold, which is visually manifested as the building structure to be corrected.
  • the candidate main angle of the object is consistent with its regional main angle, that is, it is possible to obtain a better main angle correction result by using the regional main angle through the method of the present invention.
  • the value of the regional main angle closest to the building structure to be corrected is selected from the two candidate main angles of the building structure to be corrected to obtain the real main angle of the building structure to be corrected.
  • the result of the main angle correction in the embodiment is shown in FIG. 3 .
  • This embodiment provides a method for obtaining the main angle of the area of the building structure to be corrected based on line segment histogram statistics.
  • the regional statistical strategy is segment histogram statistics
  • step S3 includes:
  • the hierarchical clustering method based on the distance matrix is used to partition the standard housing structures within the preset range:
  • the values of the elements in the distance matrix are obtained by weighting the spatial distance matrix and the angular distance matrix between buildings.
  • d i,j is the element value corresponding to the i-th building and the j-th building in the distance matrix
  • ⁇ i,j is the spatial distance between the i-th building and the j-th building
  • ⁇ i,j is the angular distance between the i-th building and the j-th building
  • w 1 is the weight value of the spatial distance
  • w 2 is the weight value of the angular distance
  • w 1 +w 2 1.
  • the shortest distance between the corner points of the boundary contours of the two buildings is taken as the spatial distance between the two buildings.
  • the formula is as follows:
  • x p,i , y p,i represent the X and Y coordinates of the pth corner point in the i-th building contour curve
  • x q,j , y q,j represent the X and Y coordinates of the qth corner point in the j-th building contour curve
  • Ni represents the number of line segments of the i-th building contour curve
  • N j represents the number of line segments of the j-th building contour curve
  • the difference between the mean angles of the boundary contours of the two buildings is taken as the angular distance between the two buildings, and the formula is as follows:
  • ⁇ k ,i represents the angle of the kth line segment in the outer contour curve of the i-th building
  • ⁇ l ,j represents the angle of the lth line segment in the outer contour curve of the j-th building
  • S i represents the number of line segments of the i-th building
  • S j represents the number of line segments of the contour curve of the j-th building.
  • House structures are composed of multiple vector contour segments
  • the angle histograms of each partition are constructed using the angles of the vector contour segments of the reference buildings and structures in each partition.
  • the angle corresponding to the highest frequency in the histogram of each partition is taken as the regional main angle of the corresponding building.
  • the target difference is greater than the second preset threshold:
  • S3 obtains the regional main angle of the building structure to be corrected by expanding the main angle statistics.
  • A is the house to be corrected, and the real main angles of other houses within a certain distance d are a 1 , a 2 , ..., a 8 .
  • the houses participating in the regional main angle calculation include all houses within the distance d of A and those intersecting with the range line.
  • the main angles of the houses and structures within the range of the house and structure to be corrected are counted, and the houses and structures with similar main angles are grouped together to obtain 4 groups of results:
  • a 1 ⁇ a 1 , a 2 , a 3 , a 4 , a 7 ⁇ ,
  • a 2 ⁇ a 5 ⁇ ,
  • a 3 ⁇ a 6 ⁇ ,
  • a 4 ⁇ a 8 ⁇ , and the group of houses and structures with the largest number is selected as A1.
  • S4 calculates the difference between the candidate main angle and the regional main angle, and selects the real main angle of the building structure to be corrected from the candidate main angles according to the relationship between the difference and the second preset threshold.
  • the target difference is greater than a second preset threshold value, and the second candidate value is selected as the true main angle of the building structure to be corrected.
  • the main angle candidate values b1 and b2 of house A are both significantly different from the regional main angle a, so the regional main angle of the house structure to be corrected cannot be used as a guiding parameter for correcting the main angle of the house structure to be corrected. Therefore, the second candidate value b2 of the main angle of the house structure to be corrected is finally used as the final true main angle.
  • the present invention proposes a new method for correcting the main angle of a building based on regional consistency.
  • the main angle of the area where the building to be corrected is located is calculated by measuring the main angle of standard buildings within a certain range around the building to be corrected.
  • the regional main angle is used as a reference to select the optimal value among the candidate main angles to obtain the true main angle.
  • the present invention solves the problem of difficulty in calculating the main angle of irregular buildings, thereby greatly improving the post-processing optimization effect of building extraction in remote sensing images.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Astronomy & Astrophysics (AREA)
  • Image Analysis (AREA)

Abstract

The present invention relates to the field of remote sensing image processing. Disclosed is a regional-consistency-based building principal angle correction method. The method comprises acquiring a regional principal angle and two candidate principal angles of a building structure, and selecting a real principal angle from the candidate principal angles by means of the relationships between the regional principal angle and the candidate principal angles, so as to correct the principal angle of the building. The present method as a post-processing method for building extraction corrects the principal angle of an irregular building structure according to the alignment consistency of building structures in the range of a region where the irregular building structure is located, thus improving the accuracy of principal angle calculation in post-processing, and thereby greatly improving the post-processing optimization effect of building structure extraction.

Description

基于区域一致性的建筑物主角度修正方法Building main angle correction method based on regional consistency 技术领域Technical Field
本发明涉及遥感影像处理领域,具体涉及一种基于区域一致性的建筑物主角度修正方法。The invention relates to the field of remote sensing image processing, and in particular to a method for correcting a main angle of a building based on regional consistency.
背景技术Background technique
目前,随着智慧城市大数据时代的到来、以及卫星遥感技术的不断发展,对地观测、卫星遥感、生态评估、国土监管逐渐向宏观、动态、精细化方向发展,因而,大范围高效率的卫星遥感影像解译具有现实需求。高分辨率遥感影像中的建筑物提取对违建监测、城区自动提取、地图更新、城市变化监测、城市规划、三维建模、数字化城市建立等应用具有重要意义。高分辨率遥感影像在提高地物光谱特征,突出地物结构、纹理和细节等信息的同时,也因为卫星观测角度问题造成地物遮挡,尺度的增大带来了严重的异物同谱问题,同时增加了影像的噪声,因而限制了建筑物提取的精度,使得目视解译方法仍然是最普遍使用的判读方法,目视解译方法虽然精度相对有保证,但是其效率低、费时费力的缺点严重制约了高分辨率遥感影像的大规模应用,造成影像数据的极大浪费。At present, with the advent of the era of smart city big data and the continuous development of satellite remote sensing technology, earth observation, satellite remote sensing, ecological assessment, and land supervision are gradually developing in the direction of macro, dynamic, and refined. Therefore, there is a practical need for large-scale and efficient satellite remote sensing image interpretation. Building extraction in high-resolution remote sensing images is of great significance for applications such as illegal construction monitoring, automatic extraction of urban areas, map updates, urban change monitoring, urban planning, three-dimensional modeling, and digital city establishment. While high-resolution remote sensing images improve the spectral characteristics of objects and highlight information such as the structure, texture, and details of objects, they also cause object occlusion due to satellite observation angle problems. The increase in scale brings serious problems of heterogeneous objects with the same spectrum, and increases the noise of the image, thus limiting the accuracy of building extraction. Visual interpretation methods are still the most commonly used interpretation method. Although the accuracy of visual interpretation methods is relatively guaranteed, their low efficiency and time-consuming and labor-intensive shortcomings seriously restrict the large-scale application of high-resolution remote sensing images, resulting in a huge waste of image data.
深度学习是人工智能中机器学习发展的新阶段,有效的解决了对复杂对象特征的刻画和复杂场景的关联分析等问题。面向高分辨率遥感影像建筑物提取的深度学习方法,可以自动提取建筑物的特征信息,实现高精度高效率的建筑物提取。然而,由于遥感影像的复杂性,建筑物受噪声、遮挡、阴影、低对比度的影响,目前的建筑物自动提取方法结果可靠性不高,存在明显的缺陷,如建筑物与道路的光谱特征相似因而无法被有效区分开来,无法通过全自动的方法直接获取精度较高的提取结果。现有的另一种方式是将计算机自动提取技术与人工交互相结合,即交互式地物提取,提取精度有了一定的保障。但是,这种方法需要人工给定建筑物的初始位置,依赖边缘信息提取建筑物,交互复杂,对人工给定位置的精度要求高,实施困难。Deep learning is a new stage in the development of machine learning in artificial intelligence, which effectively solves the problems of characterizing the features of complex objects and analyzing the association of complex scenes. The deep learning method for extracting buildings from high-resolution remote sensing images can automatically extract the characteristic information of buildings and achieve high-precision and high-efficiency building extraction. However, due to the complexity of remote sensing images, buildings are affected by noise, occlusion, shadows, and low contrast. The current automatic building extraction method has low reliability and obvious defects. For example, the spectral characteristics of buildings and roads are similar and cannot be effectively distinguished, and it is impossible to directly obtain high-precision extraction results through fully automatic methods. Another existing method is to combine computer automatic extraction technology with manual interaction, that is, interactive object extraction, and the extraction accuracy is guaranteed to a certain extent. However, this method requires the initial position of the building to be given manually, relies on edge information to extract buildings, has complex interactions, requires high accuracy for the manually given position, and is difficult to implement.
为进一步提高建筑物自动化提取的精度,并将提取结果应用于工程实践中,研究人员在建筑物提取的后处理方向展开了各种研究。基于建筑物的几何形状、空间分布等规则特征,对建筑物 提取结果进行图形化的修正处理,优化后的成果可以直接用于具体工程项目。In order to further improve the accuracy of building automatic extraction and apply the extraction results to engineering practice, researchers have carried out various studies in the post-processing direction of building extraction. The extracted results are graphically modified and the optimized results can be directly used in specific engineering projects.
在后处理方法中,通过计算获得的房屋建构筑物的主角度能够帮助提升深度学习方法提取建筑物的精度。然而现有的主角度计算方法仍然存在两个难以克服的问题,其一是由于建筑物本身形态多样,再加上自动化提取时带来的误差,给建筑物主角度的求取带来一定难度;其二是不同的主角度求取方法都有一定的适用范围和局限性,特别是对于提取效果不好建筑物,不同求取方法求取结果差异较大,因而难以取得较好的后处理效果。In the post-processing method, the main angle of the building structure obtained by calculation can help improve the accuracy of the building extraction by the deep learning method. However, the existing main angle calculation method still has two difficult problems to overcome. First, due to the diverse shapes of the buildings themselves, coupled with the errors caused by automated extraction, it is difficult to obtain the main angle of the building; second, different main angle calculation methods have certain scopes of application and limitations, especially for buildings with poor extraction effects, the results obtained by different calculation methods are quite different, making it difficult to achieve good post-processing results.
发明内容Summary of the invention
本发明提出一种基于区域一致性的建筑物主角度修正方法。该方法根据不规则房屋建构筑物所在区域范围内房屋建构筑物的走向一致性对该不规则房屋建构筑物的主角度进行修正,提高后处理中主角度计算的准确率,从而能够大大提升房屋建构筑物提取的后处理优化效果。The present invention proposes a method for correcting the main angle of a building based on regional consistency. The method corrects the main angle of the irregular building according to the consistency of the orientation of the buildings in the area where the irregular building is located, thereby improving the accuracy of the main angle calculation in post-processing, thereby greatly improving the post-processing optimization effect of building extraction.
实现上述技术目的,本发明的技术方案如下:To achieve the above technical objectives, the technical solution of the present invention is as follows:
一种基于区域一致性的建筑物主角度修正方法,该方法包括以下步骤:A method for correcting a main angle of a building based on regional consistency, the method comprising the following steps:
S1获取目标遥感影像,通过预设方法获取目标遥感影像的房屋建构筑物的候选主角度,所述候选主角度包括第一候选值和第二候选值;S1 obtains a target remote sensing image, and obtains candidate main angles of buildings and structures in the target remote sensing image by a preset method, wherein the candidate main angles include a first candidate value and a second candidate value;
S2计算第一候选值与第二候选值的差值,根据该差值和第一预设阈值的关系,将房屋建构筑物分为标准房屋建构筑物和待修正房屋建构筑物,同时获取标准房屋建构筑物的真实主角度;S2 calculates the difference between the first candidate value and the second candidate value, and according to the relationship between the difference and the first preset threshold, divides the building structure into a standard building structure and a building structure to be corrected, and simultaneously obtains the true main angle of the standard building structure;
S3通过区域统计策略得到待修正房屋建构筑物的区域主角度,区域统计策略包括外扩主角度统计和线段直方图统计;S3 obtains the regional main angle of the building structure to be corrected through a regional statistical strategy, and the regional statistical strategy includes external expansion main angle statistics and line segment histogram statistics;
S4计算候选主角度与区域主角度的差值,根据该差值与第二预设阈值的关系,从候选主角度中筛选出待修正房屋建构筑物的真实主角度。S4 calculates the difference between the candidate main angle and the regional main angle, and selects the real main angle of the building structure to be corrected from the candidate main angles according to the relationship between the difference and the second preset threshold.
进一步地,步骤S1包括:Further, step S1 comprises:
S11获取房屋建构筑物的最小面积的外接矩形的两长轴中心线的角度,作为房屋建构筑物的第一候选值;S11 obtains the angles of the two major axis center lines of the circumscribed rectangle of the minimum area of the building structure as the first candidate value of the building structure;
S12通过矢量轮廓提取法获取房屋建构筑物的矢量轮廓线段,并获取矢量轮廓线段的角度;S12 obtains the vector contour line segment of the building structure by a vector contour extraction method, and obtains the angle of the vector contour line segment;
S13将相似角度的矢量轮廓线段归为一组,筛选出总长度最长的一组矢量轮廓线段作为参考矢量轮廓线段组,基于参考矢量轮廓线段的角度,通过垂直等效法得到房屋建构筑物的第二候选值;S13 groups vector contour line segments with similar angles into one group, selects a group of vector contour line segments with the longest total length as a reference vector contour line segment group, and obtains a second candidate value of the building structure by a vertical equivalent method based on the angle of the reference vector contour line segment;
所述相似角度的矢量轮廓线段包括角度差在[0°,2°]∪[88°,90°] 的矢量轮廓线段。The vector contour segments of similar angles include those with angle differences between [0°, 2°]∪[88°, 90°] Vector contour segments of .
进一步地,所述矢量轮廓提取法包括:Furthermore, the vector contour extraction method comprises:
通过对遥感图像进行语义分割、二值化处理和数学形态学处理,得到后处理图像;By performing semantic segmentation, binarization and mathematical morphology processing on remote sensing images, post-processing images are obtained;
通过对后处理图像进行边界追踪和矢量提取,得到房屋建构筑物的矢量轮廓线段。By performing boundary tracing and vector extraction on the post-processed image, the vector contour segments of the building structure are obtained.
进一步地,所述垂直等效法包括:Furthermore, the vertical equivalent method includes:
基于参考矢量轮廓线段的角度对参考矢量轮廓线段组进行聚类,得到主矢量轮廓线段组和辅矢量轮廓线段组,所述主矢量轮廓线段组的线段总长度大于所述辅矢量轮廓线段组的线段总长度,分别对主矢量轮廓线段组的角度和辅矢量轮廓线段组的角度求均值,得到主矢量角度和辅矢量角度;Clustering the reference vector contour line segment group based on the angle of the reference vector contour line segment to obtain a main vector contour line segment group and an auxiliary vector contour line segment group, wherein the total length of the main vector contour line segment group is greater than the total length of the auxiliary vector contour line segment group, and averaging the angles of the main vector contour line segment group and the auxiliary vector contour line segment group to obtain a main vector angle and an auxiliary vector angle;
将辅矢量角度向主矢量角度的方向进行垂直变换,得到垂直辅矢量角度,将主矢量角度和垂直辅矢量角度的均值作为第二候选值。The auxiliary vector angle is vertically transformed toward the main vector angle to obtain a vertical auxiliary vector angle, and the average of the main vector angle and the vertical auxiliary vector angle is used as the second candidate value.
进一步地,步骤S2包括:Further, step S2 includes:
若第一候选值和第二候选值的差值小于第一预设阈值,将该房屋建构筑物作为标准房屋建构筑物,并将第二候选值作为标准房屋建构筑物的真实主角度;If the difference between the first candidate value and the second candidate value is less than a first preset threshold, the building structure is used as a standard building structure, and the second candidate value is used as the true main angle of the standard building structure;
若第一候选值和第二候选值的差值大于第一预设阈值,将该房屋建构筑物作为待修正房屋建构筑物。If the difference between the first candidate value and the second candidate value is greater than the first preset threshold, the building structure is used as the building structure to be corrected.
进一步地,所述区域统计策略为外扩主角度统计,步骤S3包括:Furthermore, the regional statistics strategy is outward expansion main angle statistics, and step S3 includes:
对待修正房屋建构筑物的最小外接矩形外扩预设距离,得到外扩矩形;The minimum circumscribed rectangle of the building structure to be modified is expanded by a preset distance to obtain an expanded rectangle;
将外扩矩形内以及与外扩矩形相交的标准房屋建构筑物作为参考房屋建构筑物;The standard housing structures within the outer expansion rectangle and intersecting with the outer expansion rectangle are used as reference housing structures;
将参考房屋建构筑物的真实主角度的均值,作为待修正房屋建构筑物的区域主角度。The average of the real main angles of the reference building structures is used as the regional main angle of the building structure to be corrected.
进一步地,所述区域统计策略为线段直方图统计,步骤S3包括:Furthermore, the regional statistical strategy is line segment histogram statistics, and step S3 includes:
通过基于距离矩阵的层次聚类方法对目标遥感影像中房屋建构筑物进行分区;The houses and structures in the target remote sensing image are partitioned by hierarchical clustering method based on distance matrix;
获取各分区中房屋建构筑物的矢量轮廓线段的角度,所述房屋建构筑物由多条矢量轮廓线段组成;Obtaining the angle of the vector contour line segment of the building structure in each partition, wherein the building structure is composed of a plurality of vector contour line segments;
利用各分区中房屋建构筑物的矢量轮廓线段的角度分别构造各分区的角度直方图;The angle histogram of each partition is constructed by using the angle of the vector contour line segment of the building structure in each partition;
将待修正房屋建构筑物所在分区的角度直方图中最高频率对 应的角度作为待修正房屋建构筑物的区域主角度。The highest frequency pair in the angle histogram of the partition where the building to be corrected is located The corresponding angle is used as the main angle of the area of the building structure to be modified.
进一步地,步骤S4包括:Further, step S4 includes:
分别计算第一候选值与区域主角度的差值,以及第二候选值与区域主角度的差值,选择较小的差值作为目标差值;Calculate the difference between the first candidate value and the main angle of the region, and the difference between the second candidate value and the main angle of the region, and select the smaller difference as the target difference;
若目标差值小于第二预设阈值,选择与区域主角度接近的候选主角度作为待修正房屋建构筑物的真实主角度;If the target difference is less than a second preset threshold, a candidate main angle close to the regional main angle is selected as the true main angle of the building structure to be corrected;
若目标差值大于第二预设阈值,选择第二候选值作为待修正房屋建构筑物的真实主角度。If the target difference is greater than a second preset threshold, the second candidate value is selected as the true main angle of the building structure to be corrected.
本发明的有益效果为:本发明提出了一种全新的基于区域一致性的建筑物主角度修正方法。通过区域统计策略计算得到待修正房屋建构筑物的区域主角度,以区域主角度作为基准,选取候选主角度的最优值,得到真实主角度。本发明作为建筑物提取的后处理方法,解决了不规则建筑物主角度计算困难的问题,从而能够大大提升遥感影像中建筑物提取的后处理优化效果。The beneficial effects of the present invention are as follows: the present invention proposes a new method for correcting the main angle of a building based on regional consistency. The regional main angle of the building structure to be corrected is calculated by a regional statistical strategy, and the optimal value of the candidate main angle is selected based on the regional main angle to obtain the true main angle. As a post-processing method for building extraction, the present invention solves the problem of difficulty in calculating the main angle of irregular buildings, thereby greatly improving the post-processing optimization effect of building extraction in remote sensing images.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.
图1为本发明的一种基于区域一致性的建筑物主角度修正方法的流程图;FIG1 is a flow chart of a method for correcting a main angle of a building based on regional consistency according to the present invention;
图2为本发明的实施例一提供的建筑物原始分布示意图;FIG2 is a schematic diagram of the original distribution of buildings provided in Example 1 of the present invention;
图3为本发明的实施例一中建筑物主角度修正结果示意图;FIG3 is a schematic diagram of a correction result of a main angle of a building in Embodiment 1 of the present invention;
图4为本发明的实施例三提供的建筑物原始分布示意图。FIG. 4 is a schematic diagram of the original distribution of buildings provided in the third embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员所获得的所有其他实施例,都属于本发明保护的范围。The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field belong to the scope of protection of the present invention.
请参照图1,图1是本发明所述的一种基于区域一致性的建筑物主角度修正方法实施例的流程示意图,该方法包括以下步骤:Please refer to FIG. 1 , which is a flow chart of an embodiment of a method for correcting a main angle of a building based on regional consistency according to the present invention. The method comprises the following steps:
S1通过预设方法获取房屋建构筑物的候选主角度,所述候选主角度包括第一候选值和第二候选值;S1 obtains a candidate main angle of a building structure by a preset method, wherein the candidate main angle includes a first candidate value and a second candidate value;
S2根据第一候选值与第二候选值的差值,将房屋建构筑物分为标准房屋建构筑物和待修正房屋建构筑物; S2 classifies the housing structures into standard housing structures and housing structures to be corrected according to the difference between the first candidate value and the second candidate value;
S3获取预设范围的标准房屋建构筑物的真实主角度,并根据预设范围的标准房屋建构筑物的真实主角度计算得到待修正房屋建构筑物的区域主角度;S3: obtaining the real main angle of the standard building structure within the preset range, and calculating the regional main angle of the building structure to be corrected according to the real main angle of the standard building structure within the preset range;
S4在待修正房屋建构筑物的候选主角度中选取最接近待修正房屋建构筑物的区域主角度的数值,得到待修正房屋建构筑物的真实主角度。S4 selects the value of the main angle of the region closest to the building structure to be corrected from the candidate main angles of the building structure to be corrected, and obtains the real main angle of the building structure to be corrected.
本发明的技术思路为:1.目前的基于遥感影像的建筑物提取方法仅能提取较为规则且特征明显的建筑物,通用性较差,当建筑物较为密集或是遇到不规则建筑物的情况下,提取效果往往不理想,需要合适的后处理步骤来优化提取的建筑物的精度。其中,获取建筑物的主角度对后处理优化有较大帮助。通常获取建筑物主角度的方法包括:1)获取建构筑物的最小外接矩形,最小外接矩形两长轴中心线的方向作为建构筑物主角度的第一候选值;2)获取建构筑物的轮廓线段,计算所有相似角度轮廓线段的总长度,总长度最长的方向作为建构筑物主角度的第二候选值。对于不规则建筑物,不同的计算主角度的方法获得的结果也可能不同,如何进一步确定最真实的主角度是本发明后续步骤所着重介绍的;2.本发明引入区域主角度的概念来帮助确定建筑物的真实主角度。区域主角度为建筑物所在一定区域范围内拥有相似主角度的一组房屋建构筑物的主角度的平均值。区域主角度可以反映建筑区的区域一致性,并对与区域一致性保持一致的待修正建筑物的主角度提供指导作用;3.选取最接近区域主角度的待修正建筑物的候选主角度,作为真实主角度。该方法以建筑区的区域一致性为指导,能够确定待修正不规则建筑物的真实主角度,并进一步提高建筑物提取后处理的效果,从而实现高精度的建筑物提取。The technical idea of the present invention is as follows: 1. The current building extraction method based on remote sensing images can only extract relatively regular buildings with obvious features, and its versatility is poor. When the buildings are dense or irregular buildings are encountered, the extraction effect is often not ideal, and appropriate post-processing steps are required to optimize the accuracy of the extracted buildings. Among them, obtaining the main angle of the building is of great help to post-processing optimization. The method of obtaining the main angle of the building usually includes: 1) obtaining the minimum circumscribed rectangle of the building, and the direction of the two major axis center lines of the minimum circumscribed rectangle is used as the first candidate value of the main angle of the building; 2) obtaining the contour line segment of the building, calculating the total length of all similar angle contour line segments, and the direction with the longest total length is used as the second candidate value of the main angle of the building. For irregular buildings, different methods for calculating the main angle may also obtain different results. How to further determine the most realistic main angle is emphasized in the subsequent steps of the present invention; 2. The present invention introduces the concept of regional main angle to help determine the real main angle of the building. The regional main angle is the average value of the main angle of a group of buildings with similar main angles within a certain area where the building is located. The regional main angle can reflect the regional consistency of the building area and provide guidance for the main angle of the building to be corrected that is consistent with the regional consistency; 3. Select the candidate main angle of the building to be corrected that is closest to the regional main angle as the true main angle. This method is guided by the regional consistency of the building area, can determine the true main angle of the irregular building to be corrected, and further improve the effect of building extraction post-processing, thereby achieving high-precision building extraction.
实施例一Embodiment 1
本实施例提供一种基于外扩主角度统计获取待修正房屋建构筑物的区域主角度的方法,同时,本实施例提供在目标差值小于第二预设阈值的情况下,得到待修正房屋建构筑物的真实主角度的方法:This embodiment provides a method for obtaining the regional main angle of the building structure to be corrected based on the statistics of the external expansion main angle. At the same time, this embodiment provides a method for obtaining the real main angle of the building structure to be corrected when the target difference is less than the second preset threshold:
S1获取目标遥感影像,通过预设方法获取目标遥感影像的房屋建构筑物的候选主角度,所述候选主角度包括第一候选值和第二候选值;S1 obtains a target remote sensing image, and obtains candidate main angles of buildings and structures in the target remote sensing image by a preset method, wherein the candidate main angles include a first candidate value and a second candidate value;
S11获取房屋建构筑物的最小面积的外接矩形的两长轴中心线的角度,作为房屋建构筑物的第一候选值;S11 obtains the angles of the two major axis center lines of the circumscribed rectangle of the minimum area of the building structure as the first candidate value of the building structure;
S12通过矢量轮廓提取法获取房屋建构筑物的矢量轮廓线段,并获取矢量轮廓线段的角度;S12 obtains the vector contour line segment of the building structure by a vector contour extraction method, and obtains the angle of the vector contour line segment;
S13将相似角度的矢量轮廓线段归为一组,筛选出总长度最长 的一组矢量轮廓线段作为参考矢量轮廓线段组,基于参考矢量轮廓线段的角度,通过垂直等效法得到房屋建构筑物的第二候选值;S13 groups vector contour segments with similar angles into a group and selects the segment with the longest total length. A set of vector contour line segments is used as a reference vector contour line segment set, and based on the angles of the reference vector contour line segments, a second candidate value of the building structure is obtained by a vertical equivalent method;
所述相似角度的矢量轮廓线段为角度差在[0°,2°]∪[88°,90°]的矢量轮廓线段。The vector contour line segments of similar angles are vector contour line segments whose angle difference is [0°, 2°]∪[88°, 90°].
对步骤S12进行补充解释,矢量轮廓提取法包括:To provide additional explanation for step S12, the vector contour extraction method includes:
通过深度学习方法对目标遥感图像进行语义分割,得到分割数据;Perform semantic segmentation on the target remote sensing image through deep learning methods to obtain segmentation data;
本实施例使用mask-RCNN对目标遥感图像进行语义分割。mask-RCNN沿用了Faster RCNN的思想,特征提取采用ResNet-FPN的架构,另外多加了一个Mask预测分支,能够完成较高质量的语义分割任务,对于普通建筑物的提取效果较好,仅在面对复杂的不规则建筑物时效果明显下降。本实施例使用ArcGIS通过手动样本勾画的方式进行样本制作。训练过程中使用裁剪为512×512的像素块作为样本,使用了约7000个实例作为训练数据集,约3000个实例作为验证数据集。计算机配置:CPU:Xeon(R)CPU E5-2620 v4@2.10GHz×32;显卡:Quadro M4000;内存:128G;操作系统:Ubuntu 16.04。训练的过程总花费40多个小时,迭代9万次;This embodiment uses mask-RCNN to perform semantic segmentation on the target remote sensing image. Mask-RCNN follows the idea of Faster RCNN, and uses the architecture of ResNet-FPN for feature extraction. In addition, a Mask prediction branch is added, which can complete high-quality semantic segmentation tasks. It has a good extraction effect on ordinary buildings, and the effect is significantly reduced only when facing complex irregular buildings. This embodiment uses ArcGIS to make samples by manual sample outlining. During the training process, pixel blocks cropped to 512×512 are used as samples, and about 7,000 instances are used as training data sets and about 3,000 instances are used as verification data sets. Computer configuration: CPU: Xeon(R)CPU E5-2620 v4@2.10GHz×32; Graphics card: Quadro M4000; Memory: 128G; Operating system: Ubuntu 16.04. The training process took more than 40 hours and 90,000 iterations.
对分割数据进行二值化处理,以0.5作为建筑物目标概率图的二值化阈值,得到分割二值化图像;The segmented data is binarized, and 0.5 is used as the binarization threshold of the building target probability map to obtain a segmented binary image;
房屋建构筑物的分割二值化图像用来凸显目标建筑物的轮廓并压缩整体图像数据量;The segmentation binary image of the building structure is used to highlight the outline of the target building and compress the overall image data volume;
通过数学形态学方法对分割二值化图像进行后处理,得到后处理图像:The segmented binary image is post-processed by mathematical morphology method to obtain the post-processed image:
其中,数学形态学方法由一组形态学的代数运算子组成,基本运算包括:膨胀、腐蚀、开启和闭合;Among them, the mathematical morphology method consists of a set of morphological algebraic operators, and the basic operations include: dilation, erosion, opening and closing;
通过数学形态学方法对分割二值化图像进行去噪处理,去除分类噪声和小型非建筑结构等无用的信息,基于数学形态学方法可以进行图像形状和结构的分析及处理,包括图像分割、特征抽取、边缘检测、图像滤波、图像增强和恢复。The segmented binary image is denoised by mathematical morphology methods to remove useless information such as classification noise and small non-building structures. Image shape and structure analysis and processing can be performed based on mathematical morphology methods, including image segmentation, feature extraction, edge detection, image filtering, image enhancement and restoration.
通过摩尔领域边界追踪法对后处理图像追踪边界生成房屋建构筑物的轮廓线段;The boundary of the post-processed image is traced by Moore's field boundary tracing method to generate the contour line segments of the building structure;
基于房屋建构筑物的轮廓线段对后处理图像进行矢量提取,得到房屋建构筑物的矢量轮廓线段。Based on the contour line segments of the building structures, vector extraction is performed on the post-processed image to obtain the vector contour line segments of the building structures.
对步骤S12进行补充解释,垂直等效法包括:To provide additional explanation for step S12, the vertical equivalent method includes:
基于参考矢量轮廓线段的角度对参考矢量轮廓线段组进行聚类,得到主矢量轮廓线段组和辅矢量轮廓线段组;Clustering the reference vector contour line segment group based on the angle of the reference vector contour line segment to obtain a main vector contour line segment group and an auxiliary vector contour line segment group;
所述主矢量轮廓线段组的线段总长度大于所述辅矢量轮廓线 段组的线段总长度;The total length of the line segments of the main vector contour line segment group is greater than that of the auxiliary vector contour line segment group. The total length of the line segments of the segment group;
分别对主矢量轮廓线段组和辅矢量轮廓线段组求角度均值,得到主矢量角度和辅矢量角度;Calculate the angle mean of the main vector contour line segment group and the auxiliary vector contour line segment group respectively to obtain the main vector angle and the auxiliary vector angle;
所述主矢量轮廓线段组和辅矢量轮廓线段组的组内线段角度差在[0°,2°],组间线段角度差在[88°,90°],故主矢量角度和辅矢量角度接近于垂直的关系;The angle difference between the main vector contour line segment group and the auxiliary vector contour line segment group is [0°, 2°], and the angle difference between the groups is [88°, 90°], so the main vector angle and the auxiliary vector angle are close to a vertical relationship;
将辅矢量角度向主矢量角度的方向进行垂直变换,得到垂直辅矢量角度;Transform the auxiliary vector angle vertically toward the direction of the main vector angle to obtain the vertical auxiliary vector angle;
垂直变换表示如下:
The vertical transformation is expressed as follows:
其中,xi为主矢量角度,yi为辅矢量角度,f(xi)为垂直辅矢量角度;Among them, xi is the main vector angle, yi is the auxiliary vector angle, and f(xi) is the vertical auxiliary vector angle;
经过垂直变化操作后,得到的垂直辅矢量角度与主矢量角度的角度差在[0°,2°];After the vertical change operation, the angle difference between the vertical auxiliary vector angle and the main vector angle is [0°, 2°];
将主矢量角度和垂直辅矢量角度的均值作为第二候选值b2。The average of the main vector angle and the vertical auxiliary vector angle is used as the second candidate value b2.
S2计算第一候选值与第二候选值的差值,根据该差值和第一预设阈值的关系,将房屋建构筑物分为标准房屋建构筑物和待修正房屋建构筑物。S2 calculates the difference between the first candidate value and the second candidate value, and divides the building structures into standard building structures and building structures to be corrected according to the relationship between the difference and the first preset threshold.
若第一候选值和第二候选值的差值小于第一预设阈值,将该房屋建构筑物作为标准房屋建构筑物,并将第二候选值作为标准房屋建构筑物的真实主角度;If the difference between the first candidate value and the second candidate value is less than a first preset threshold, the building structure is used as a standard building structure, and the second candidate value is used as the true main angle of the standard building structure;
若第一候选值和第二候选值的差值大于第一预设阈值,将该房屋建构筑物作为待修正房屋建构筑物。If the difference between the first candidate value and the second candidate value is greater than the first preset threshold, the building structure is used as the building structure to be corrected.
对步骤S2进行补充解释,本具体实施例中,预设阈值设置为3°,若第一候选值和第二候选值的差值小于3°,则认为该房屋建构筑物的主方向不需要进行修正,将该房屋建构筑物作为标准房屋建构筑物,将其第二候选值作为标准房屋建构筑物的真实主角度。To provide additional explanation for step S2, in this specific embodiment, the preset threshold is set to 3°. If the difference between the first candidate value and the second candidate value is less than 3°, it is considered that the main direction of the building structure does not need to be corrected, and the building structure is taken as a standard building structure, and the second candidate value is taken as the true main angle of the standard building structure.
S3获取预设范围的标准房屋建构筑物的真实主角度,并根据预设范围的标准房屋建构筑物的形状信息得到待修正房屋建构筑物的区域主角度。S3 obtains the real main angle of the standard building structure within the preset range, and obtains the regional main angle of the building structure to be corrected according to the shape information of the standard building structure within the preset range.
在实施例一中,所述区域统计策略为外扩主角度统计,步骤S3包括:In the first embodiment, the regional statistical strategy is external expansion main angle statistics, and step S3 includes:
获取参考房屋建构筑物的真实主角度a2、a3、a4、a5、a6、a7Obtaining the real main angles a 2 , a 3 , a 4 , a 5 , a 6 , and a 7 of the reference building structure;
请参阅图2,A为待修正房屋建构筑物,周边一定距离d范围内其他房屋为参考房屋建构筑物,参与区域主角度计算的房屋包 括A距离d范围内以及与范围线相交的所有房屋。Please refer to Figure 2, A is the building to be corrected, and other houses within a certain distance d are reference buildings. The houses involved in the main angle calculation of the area include Includes all houses within a distance d from A and those intersecting the range line.
基于真实主角度通过近邻聚类法对参考房屋建构筑物的真实主角度进行聚类:Based on the real main angle, the real main angle of the reference building structure is clustered by the nearest neighbor clustering method:
(1)从以上6个参考房屋建构筑物的真实主角度中任取一个真实主角度作为第一个聚类中心,如令z1=a2,其中z1为参考组A1的聚类中心;(1) Select any one of the true main angles of the above six reference building structures as the first cluster center, such as z 1 = a 2 , where z 1 is the cluster center of reference group A 1 ;
(2)当前聚类中心为z1,计算a3到聚类中心z1的Euclidean距离dist(a3,z1);(2) The current cluster center is z 1 , and the Euclidean distance from a 3 to the cluster center z 1 is calculated as dist(a 3 ,z 1 );
若dist(a3,z1)∈[0,3]或∪[87,93],则a2∈A1If dist(a3,z1)∈[0,3] or ∪[87,93], then a 2 ∈A 1 ;
(3)当前聚类中心为z1,计算a5到聚类中心z1的Euclidean距离dist(a4,z1);(3) The current cluster center is z 1 , and the Euclidean distance from a 5 to the cluster center z 1 is calculated as dist(a 4 ,z 1 );
若dist(a4,z1)∈(3,87),则将a4定义为第二个聚类中心z2,z2=a4,其中z2为参考组A2的聚类中心;If dist(a 4 ,z 1 )∈(3,87), then a 4 is defined as the second cluster center z 2 , z 2 =a 4 , where z 2 is the cluster center of the reference group A 2 ;
(1)以此类推,直至完成以上8个参考房屋建构筑物的真实主角度的聚类,聚类结果为A1={a2、a3、a5、a6、a7}},A2={a4}。(1) This process is deduced in this way until the clustering of the true main angles of the above eight reference building structures is completed. The clustering results are A1 = {a 2 , a 3 , a 5 , a 6 , a 7 }}, A 2 = {a 4 }.
其中近邻聚类法根据真实主角度的Euclidean距离进行分类,Euclidean距离表示为:
The nearest neighbor clustering method is classified according to the Euclidean distance of the true main angle, and the Euclidean distance is expressed as:
其中,dist(ai,aj)为真实主角度ai到真实主角度aj的Euclidean距离;Wherein, dist(a i ,a j ) is the Euclidean distance from the true main angle a i to the true main angle a j ;
选取组内参考房屋建构筑物的数量最多的一组参考组,计算该组内参考房屋建构筑物的真实主角度的均值,得到待修正房屋建构筑物的区域主角度;A reference group with the largest number of reference buildings and structures is selected, and the mean of the true main angles of the reference buildings and structures in the group is calculated to obtain the regional main angle of the building and structure to be corrected;
其中组内参考房屋建构筑物的数量最多的一组参考组即为A1,计算该组房屋建构筑物主角度的均值,得到待修正房屋建构筑物的区域主角度 The reference group with the largest number of reference buildings and structures is A1 . The mean of the main angles of the buildings and structures in this group is calculated to obtain the main angle of the area of the buildings and structures to be corrected.
S4在待修正房屋建构筑物的候选主角度中选取最接近待修正房屋建构筑物的区域主角度的数值,得到待修正房屋建构筑物的真实主角度。S4 selects the value of the main angle of the region closest to the building structure to be corrected from the candidate main angles of the building structure to be corrected, and obtains the real main angle of the building structure to be corrected.
在实施例一中,目标差值小于第二预设阈值,步骤S4包括:In the first embodiment, the target difference is less than the second preset threshold, and step S4 includes:
分别计算第一候选值与区域主角度的差值,以及第二候选值与区域主角度的差值,选择较小的差值作为目标差值;Calculate the difference between the first candidate value and the main angle of the region, and the difference between the second candidate value and the main angle of the region, and select the smaller difference as the target difference;
目标差值小于第二预设阈值,选择目标差值对应的候选主角度作为待修正房屋建构筑物的真实主角度。The target difference is smaller than the second preset threshold, and the candidate main angle corresponding to the target difference is selected as the true main angle of the building structure to be corrected.
本实施例中区域主角度的获取情况如图2所示。该实施例中目标差值小于第二预设阈值,在视觉上表现为待修正房屋建构筑 物的候选主角度与其区域主角度保持一致,即因而能够通过本发明的方法借助区域主角度获得较好的主角度修正结果。The acquisition of the main angle of the area in this embodiment is shown in Figure 2. In this embodiment, the target difference is less than the second preset threshold, which is visually manifested as the building structure to be corrected. The candidate main angle of the object is consistent with its regional main angle, that is, it is possible to obtain a better main angle correction result by using the regional main angle through the method of the present invention.
接着在待修正房屋建构筑物主角度的2个候选主角度中选取最接近待修正房屋建构筑物的区域主角度的数值,得到待修正房屋建构筑物的真实主角度。Then, the value of the regional main angle closest to the building structure to be corrected is selected from the two candidate main angles of the building structure to be corrected to obtain the real main angle of the building structure to be corrected.
在实施例中主角度修正的结果如图3所示。The result of the main angle correction in the embodiment is shown in FIG. 3 .
实施例二Embodiment 2
本实施例提供一种基于线段直方图统计的获取待修正房屋建构筑物的区域主角度的方法。This embodiment provides a method for obtaining the main angle of the area of the building structure to be corrected based on line segment histogram statistics.
在实施例二中,所述区域统计策略为线段直方图统计,步骤S3包括:In the second embodiment, the regional statistical strategy is segment histogram statistics, and step S3 includes:
通过基于距离矩阵的层次聚类方法对预设范围的标准房屋建构筑物进行分区:The hierarchical clustering method based on the distance matrix is used to partition the standard housing structures within the preset range:
所述距离矩阵中元素的值由建筑物之间的空间距离矩阵、角距离矩阵进行加权获得,公式如下:
di,j=W1Δρi,j+w2Δθi,j
The values of the elements in the distance matrix are obtained by weighting the spatial distance matrix and the angular distance matrix between buildings. The formula is as follows:
d i,j =W 1 Δρ i,j +w 2 Δθ i,j
其中,di,j为在距离矩阵中与第i个建筑物和第j个建筑物对应的元素值,Δρi,j为第i个建筑物与第j个建筑物之间的空间距离,Δθi,j为第i个建筑物与第j个建筑物之间的角距离,w1为空间距离的权重值,w2为角距离的权重值,w1+w2=1。Among them, d i,j is the element value corresponding to the i-th building and the j-th building in the distance matrix, Δρ i,j is the spatial distance between the i-th building and the j-th building, Δθ i,j is the angular distance between the i-th building and the j-th building, w 1 is the weight value of the spatial distance, w 2 is the weight value of the angular distance, and w 1 +w 2 =1.
其中,将两建筑物边界轮廓之间角点的最短距离作为两建筑物之间的空间距离,公式如下:
Among them, the shortest distance between the corner points of the boundary contours of the two buildings is taken as the spatial distance between the two buildings. The formula is as follows:
其中,xp,i、yp,i代表第i个建筑物轮廓曲线中第p个角点的X、Y坐标;xq,j、yq,j代表着第j个建筑物轮廓曲线中第q个角点的X、Y坐标,Ni代表第i个建筑物轮廓曲线的线段数量,Nj代表第j个建筑物轮廓曲线的线段数量;Where x p,i , y p,i represent the X and Y coordinates of the pth corner point in the i-th building contour curve; x q,j , y q,j represent the X and Y coordinates of the qth corner point in the j-th building contour curve; Ni represents the number of line segments of the i-th building contour curve, and N j represents the number of line segments of the j-th building contour curve;
其中,将两建筑物边界轮廓之间的角度均值之差作为两建筑物的角距离,公式如下:
Among them, the difference between the mean angles of the boundary contours of the two buildings is taken as the angular distance between the two buildings, and the formula is as follows:
式中,θk,i表示第i个建筑的外轮廓曲线中第k条线段的角度;θl,j表示第j个建筑的外轮廓曲线中第l条线段的角度;Si表示第i个建筑物的线段数量,Sj表示第j个建筑物轮廓曲线的线段数量。Wherein, θk ,i represents the angle of the kth line segment in the outer contour curve of the i-th building; θl ,j represents the angle of the lth line segment in the outer contour curve of the j-th building; S i represents the number of line segments of the i-th building, and S j represents the number of line segments of the contour curve of the j-th building.
获取参考房屋建构筑物的矢量轮廓线段的角度,所述标准房 屋建构筑物由多条矢量轮廓线段组成;Get the angle of the vector contour line segment of the reference house structure. House structures are composed of multiple vector contour segments;
利用各分区中参考房屋建构筑物的矢量轮廓线段的角度分别构造各分区的角度直方图,各分区直方图中最高频率对应的角度作为相应建筑物的区域主角度。The angle histograms of each partition are constructed using the angles of the vector contour segments of the reference buildings and structures in each partition. The angle corresponding to the highest frequency in the histogram of each partition is taken as the regional main angle of the corresponding building.
实施例三:Embodiment three:
本实施例提供在目标差值大于第二预设阈值的情况下,得到待修正房屋建构筑物的真实主角度的方法:This embodiment provides a method for obtaining the true main angle of the building structure to be corrected when the target difference is greater than the second preset threshold:
请参阅图4,在实施例三中,目标差值大于第二预设阈值:Please refer to FIG. 4 , in the third embodiment, the target difference is greater than the second preset threshold:
S3通过外扩主角度统计得到待修正房屋建构筑物的区域主角度。S3 obtains the regional main angle of the building structure to be corrected by expanding the main angle statistics.
A为待修正房屋,周边一定距离d范围内其他房屋的真实主角度为a1、a2、…、a8。参与区域主角度计算的房屋包括A距离d范围内以及与范围线相交的所有房屋。统计待修正房屋建构筑物所在范围内的房屋建构筑物的主角度,将拥有相似主角度的房屋建构筑物归为一组,得到4组结果:A1={a1、a2、a3、a4、a7},A2={a5},A3={a6},A4={a8},选取数量最多的一组房屋建构筑物即为A1。A is the house to be corrected, and the real main angles of other houses within a certain distance d are a 1 , a 2 , …, a 8 . The houses participating in the regional main angle calculation include all houses within the distance d of A and those intersecting with the range line. The main angles of the houses and structures within the range of the house and structure to be corrected are counted, and the houses and structures with similar main angles are grouped together to obtain 4 groups of results: A 1 = {a 1 , a 2 , a 3 , a 4 , a 7 }, A 2 = {a 5 }, A 3 = {a 6 }, A 4 = {a 8 }, and the group of houses and structures with the largest number is selected as A1.
计算该组房屋建构筑物主角度的平均值,得到待修正房屋建构筑物的区域主角度 Calculate the average of the main angles of the group of buildings and structures to obtain the main angle of the area of the building and structure to be corrected
S4计算候选主角度与区域主角度的差值,根据该差值与第二预设阈值的关系,从候选主角度中筛选出待修正房屋建构筑物的真实主角度。S4 calculates the difference between the candidate main angle and the regional main angle, and selects the real main angle of the building structure to be corrected from the candidate main angles according to the relationship between the difference and the second preset threshold.
分别计算第一候选值与区域主角度的差值,以及第二候选值与区域主角度的差值,选择较小的差值作为目标差值;Calculate the difference between the first candidate value and the main angle of the region, and the difference between the second candidate value and the main angle of the region, and select the smaller difference as the target difference;
目标差值大于第二预设阈值,选择第二候选值作为待修正房屋建构筑物的真实主角度。The target difference is greater than a second preset threshold value, and the second candidate value is selected as the true main angle of the building structure to be corrected.
从图4中可以看出,房屋A的主角度候选值b1和b2都与区域主角度a相差较大,因此该待修正房屋建构筑物的区域主角度无法作为待修正房屋建构筑物的主角度修正的指导参数。故最终使用待修正房屋建构筑物主角度的第二候选值b2作为最终的真实主角度。As can be seen from Figure 4, the main angle candidate values b1 and b2 of house A are both significantly different from the regional main angle a, so the regional main angle of the house structure to be corrected cannot be used as a guiding parameter for correcting the main angle of the house structure to be corrected. Therefore, the second candidate value b2 of the main angle of the house structure to be corrected is finally used as the final true main angle.
本发明的有益效果为:本发明提出了一种全新的基于区域一致性的建筑物主角度修正方法。通过待修正房屋建构筑物周围一定范围内的标准房屋建构筑物的主角度,计算得到待修正房屋建构筑物所处区域的区域主角度,以区域主角度作为基准,选取候选主角度中的最优值,得到真实主角度。本发明作为建筑物提取的后处理方法,解决了不规则建筑物主角度计算困难的问题,从而能够大大提升遥感影像中建筑物提取的后处理优化效果。 The beneficial effects of the present invention are as follows: the present invention proposes a new method for correcting the main angle of a building based on regional consistency. The main angle of the area where the building to be corrected is located is calculated by measuring the main angle of standard buildings within a certain range around the building to be corrected. The regional main angle is used as a reference to select the optimal value among the candidate main angles to obtain the true main angle. As a post-processing method for building extraction, the present invention solves the problem of difficulty in calculating the main angle of irregular buildings, thereby greatly improving the post-processing optimization effect of building extraction in remote sensing images.
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者装置不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者装置所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者装置中还存在另外的相同要素。It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or device including the element.
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。 The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

  1. 基于区域一致性的建筑物主角度修正方法,其特征在于,包括以下步骤:The method for correcting the main angle of a building based on regional consistency is characterized by comprising the following steps:
    S1获取目标遥感影像,通过预设方法获取目标遥感影像的房屋建构筑物的候选主角度,所述候选主角度包括第一候选值和第二候选值;S1 obtains a target remote sensing image, and obtains candidate main angles of buildings and structures in the target remote sensing image by a preset method, wherein the candidate main angles include a first candidate value and a second candidate value;
    S2计算第一候选值与第二候选值的差值,根据该差值和第一预设阈值的关系,将房屋建构筑物分为标准房屋建构筑物和待修正房屋建构筑物,同时获取标准房屋建构筑物的真实主角度;S2 calculates the difference between the first candidate value and the second candidate value, and according to the relationship between the difference and the first preset threshold, divides the building structure into a standard building structure and a building structure to be corrected, and simultaneously obtains the true main angle of the standard building structure;
    S3通过区域统计策略得到待修正房屋建构筑物的区域主角度,区域统计策略包括外扩主角度统计和线段直方图统计;S3 obtains the regional main angle of the building structure to be corrected through a regional statistical strategy, and the regional statistical strategy includes external expansion main angle statistics and line segment histogram statistics;
    S4计算候选主角度与区域主角度的差值,根据该差值与第二预设阈值的关系,从候选主角度中筛选出待修正房屋建构筑物的真实主角度。S4 calculates the difference between the candidate main angle and the regional main angle, and selects the real main angle of the building structure to be corrected from the candidate main angles according to the relationship between the difference and the second preset threshold.
  2. 根据权利要求1所述的基于区域一致性的建筑物主角度修正方法,其特征在于,步骤S1包括:The method for correcting the main angle of a building based on regional consistency according to claim 1, characterized in that step S1 comprises:
    S11获取房屋建构筑物的最小面积的外接矩形的两长轴中心线的角度,作为房屋建构筑物的第一候选值;S11 obtains the angles of the two major axis center lines of the circumscribed rectangle of the minimum area of the building structure as the first candidate value of the building structure;
    S12通过矢量轮廓提取法获取房屋建构筑物的矢量轮廓线段,并获取矢量轮廓线段的角度;S12 obtains the vector contour line segment of the building structure by a vector contour extraction method, and obtains the angle of the vector contour line segment;
    S13将相似角度的矢量轮廓线段归为一组,筛选出总长度最长的一组矢量轮廓线段作为参考矢量轮廓线段组,基于参考矢量轮廓线段的角度,通过垂直等效法得到房屋建构筑物的第二候选值;S13 groups vector contour line segments with similar angles into one group, selects a group of vector contour line segments with the longest total length as a reference vector contour line segment group, and obtains a second candidate value of the building structure by a vertical equivalent method based on the angle of the reference vector contour line segment;
    所述相似角度的矢量轮廓线段包括角度差在[0°,2°]∪[88°,90°]的矢量轮廓线段。The vector contour line segments with similar angles include vector contour line segments with an angle difference of [0°, 2°]∪[88°, 90°].
  3. 根据权利要求2所述的基于区域一致性的建筑物主角度修正方法,其特征在于,所述矢量轮廓提取法包括:The method for correcting the main angle of a building based on regional consistency according to claim 2, characterized in that the vector contour extraction method comprises:
    通过对遥感图像进行语义分割、二值化处理和数学形态学处理,得到后处理图像;By performing semantic segmentation, binarization and mathematical morphology processing on remote sensing images, post-processing images are obtained;
    通过对后处理图像进行边界追踪和矢量提取,得到房屋建构筑物的矢量轮廓线段。By performing boundary tracing and vector extraction on the post-processed image, the vector contour segments of the building structure are obtained.
  4. 根据权利要求2所述的基于区域一致性的建筑物主角度修正方法,其特征在于,所述垂直等效法包括:The method for correcting the main angle of a building based on regional consistency according to claim 2 is characterized in that the vertical equivalence method comprises:
    基于参考矢量轮廓线段的角度对参考矢量轮廓线段组进行聚类,得到主矢量轮廓线段组和辅矢量轮廓线段组,所述主矢量轮廓线段组的线段总长度大于所述辅矢量轮廓线段组的线段总长度,分别对主矢量轮廓线段组的角度和辅矢量轮廓线段组的角度 求均值,得到主矢量角度和辅矢量角度;The reference vector contour line segment group is clustered based on the angle of the reference vector contour line segment to obtain a main vector contour line segment group and an auxiliary vector contour line segment group, wherein the total length of the line segments of the main vector contour line segment group is greater than the total length of the line segments of the auxiliary vector contour line segment group, and the angles of the main vector contour line segment group and the angles of the auxiliary vector contour line segment group are respectively Find the mean and get the main vector angle and auxiliary vector angle;
    将辅矢量角度向主矢量角度的方向进行垂直变换,得到垂直辅矢量角度,将主矢量角度和垂直辅矢量角度的均值作为第二候选值。The auxiliary vector angle is vertically transformed toward the main vector angle to obtain a vertical auxiliary vector angle, and the average of the main vector angle and the vertical auxiliary vector angle is used as the second candidate value.
  5. 根据权利要求1所述的基于区域一致性的建筑物主角度修正方法,其特征在于,步骤S2包括:The method for correcting the main angle of a building based on regional consistency according to claim 1, characterized in that step S2 comprises:
    若第一候选值和第二候选值的差值小于第一预设阈值,将该房屋建构筑物作为标准房屋建构筑物,并将第二候选值作为标准房屋建构筑物的真实主角度;If the difference between the first candidate value and the second candidate value is less than a first preset threshold, the building structure is used as a standard building structure, and the second candidate value is used as the true main angle of the standard building structure;
    若第一候选值和第二候选值的差值大于第一预设阈值,将该房屋建构筑物作为待修正房屋建构筑物。If the difference between the first candidate value and the second candidate value is greater than the first preset threshold, the building structure is used as the building structure to be corrected.
  6. 根据权利要求1所述的基于区域一致性的建筑物主角度修正方法,其特征在于,所述区域统计策略为外扩主角度统计,步骤S3包括:The method for correcting the main angle of a building based on regional consistency according to claim 1 is characterized in that the regional statistical strategy is external expansion main angle statistics, and step S3 comprises:
    对待修正房屋建构筑物的最小外接矩形外扩预设距离,得到外扩矩形;The minimum circumscribed rectangle of the building structure to be modified is expanded by a preset distance to obtain an expanded rectangle;
    将外扩矩形内以及与外扩矩形相交的标准房屋建构筑物作为参考房屋建构筑物;The standard housing structures within the outer expansion rectangle and intersecting with the outer expansion rectangle are used as reference housing structures;
    将参考房屋建构筑物的真实主角度的均值,作为待修正房屋建构筑物的区域主角度。The average of the real main angles of the reference building structures is used as the regional main angle of the building structure to be corrected.
  7. 根据权利要求1所述的基于区域一致性的建筑物主角度修正方法,其特征在于,所述区域统计策略为线段直方图统计,步骤S3包括:The method for correcting the main angle of a building based on regional consistency according to claim 1 is characterized in that the regional statistical strategy is line segment histogram statistics, and step S3 comprises:
    通过基于距离矩阵的层次聚类方法对目标遥感影像中房屋建构筑物进行分区;The houses and structures in the target remote sensing image are partitioned by hierarchical clustering method based on distance matrix;
    获取各分区中房屋建构筑物的矢量轮廓线段的角度,所述房屋建构筑物由多条矢量轮廓线段组成;Obtaining the angle of the vector contour line segment of the building structure in each partition, wherein the building structure is composed of a plurality of vector contour line segments;
    利用各分区中房屋建构筑物的矢量轮廓线段的角度分别构造各分区的角度直方图;The angle histogram of each partition is constructed by using the angle of the vector contour line segment of the building structure in each partition;
    将待修正房屋建构筑物所在分区的角度直方图中最高频率对应的角度作为待修正房屋建构筑物的区域主角度。The angle corresponding to the highest frequency in the angle histogram of the partition where the building or structure to be corrected is located is taken as the regional main angle of the building or structure to be corrected.
  8. 根据权利要求1所述的基于区域一致性的建筑物主角度修正方法,其特征在于,步骤S4包括:The method for correcting the main angle of a building based on regional consistency according to claim 1, characterized in that step S4 comprises:
    分别计算第一候选值与区域主角度的差值,以及第二候选值与区域主角度的差值,选择较小的差值作为目标差值;Calculate the difference between the first candidate value and the main angle of the region, and the difference between the second candidate value and the main angle of the region, and select the smaller difference as the target difference;
    若目标差值小于第二预设阈值,选择与区域主角度接近的候选主角度作为待修正房屋建构筑物的真实主角度;If the target difference is less than a second preset threshold, a candidate main angle close to the regional main angle is selected as the true main angle of the building structure to be corrected;
    若目标差值大于第二预设阈值,选择第二候选值作为待修正 房屋建构筑物的真实主角度。 If the target difference is greater than the second preset threshold, the second candidate value is selected as the value to be corrected. The true main angle of the building structure.
PCT/CN2023/137754 2022-12-15 2023-12-11 Regional-consistency-based building principal angle correction method WO2024125434A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202211609894.5A CN115620169B (en) 2022-12-15 2022-12-15 Building main angle correction method based on regional consistency
CN202211609894.5 2022-12-15

Publications (1)

Publication Number Publication Date
WO2024125434A1 true WO2024125434A1 (en) 2024-06-20

Family

ID=84879859

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2023/137754 WO2024125434A1 (en) 2022-12-15 2023-12-11 Regional-consistency-based building principal angle correction method

Country Status (2)

Country Link
CN (1) CN115620169B (en)
WO (1) WO2024125434A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115620169B (en) * 2022-12-15 2023-04-07 北京数慧时空信息技术有限公司 Building main angle correction method based on regional consistency

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120239332A1 (en) * 2010-09-17 2012-09-20 Seer Technology, Inc. Building perpendicularity testing and adjustment
CN102938066A (en) * 2012-12-07 2013-02-20 南京大学 Method for reconstructing outer outline polygon of building based on multivariate data
US20130329940A1 (en) * 2012-06-07 2013-12-12 Nec Corporation Image processing apparatus, control method of image processing apparatus, and storage medium
CN114898119A (en) * 2022-07-08 2022-08-12 浙江大华技术股份有限公司 Building outline drawing method, device, equipment and medium
CN115272882A (en) * 2022-08-03 2022-11-01 山东省国土测绘院 Discrete building detection method and system based on remote sensing image
CN115620169A (en) * 2022-12-15 2023-01-17 北京数慧时空信息技术有限公司 Building main angle correction method based on regional consistency

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU2016204648A1 (en) * 2009-05-15 2016-07-21 Eagle View Technologies, Inc. Pitch determination systems and methods for aerial roof estimation
CN103065295B (en) * 2012-12-04 2016-01-20 南京大学 A kind of aviation based on buildings angle point self-correction and ground lidar data high-precision automatic method for registering
CN109583284B (en) * 2017-09-29 2023-09-12 中国科学院空天信息创新研究院 Urban high-rise building height extraction method and device based on high-resolution SAR image
CN108106594B (en) * 2017-12-06 2020-05-15 上海建工五建集团有限公司 Real-time measuring method for deformation of building
CN111652892A (en) * 2020-05-02 2020-09-11 王磊 Remote sensing image building vector extraction and optimization method based on deep learning
CN113139453B (en) * 2021-04-19 2023-04-07 中国地质大学(武汉) Orthoimage high-rise building base vector extraction method based on deep learning
CN114494905A (en) * 2022-01-26 2022-05-13 中科星图股份有限公司 Building identification and modeling method and device based on satellite remote sensing image
CN114241338A (en) * 2022-02-15 2022-03-25 中航建筑工程有限公司 Building measuring method, device, equipment and storage medium based on image recognition

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120239332A1 (en) * 2010-09-17 2012-09-20 Seer Technology, Inc. Building perpendicularity testing and adjustment
US20130329940A1 (en) * 2012-06-07 2013-12-12 Nec Corporation Image processing apparatus, control method of image processing apparatus, and storage medium
CN102938066A (en) * 2012-12-07 2013-02-20 南京大学 Method for reconstructing outer outline polygon of building based on multivariate data
CN114898119A (en) * 2022-07-08 2022-08-12 浙江大华技术股份有限公司 Building outline drawing method, device, equipment and medium
CN115272882A (en) * 2022-08-03 2022-11-01 山东省国土测绘院 Discrete building detection method and system based on remote sensing image
CN115620169A (en) * 2022-12-15 2023-01-17 北京数慧时空信息技术有限公司 Building main angle correction method based on regional consistency

Also Published As

Publication number Publication date
CN115620169A (en) 2023-01-17
CN115620169B (en) 2023-04-07

Similar Documents

Publication Publication Date Title
WO2024077812A1 (en) Single building three-dimensional reconstruction method based on point cloud semantic segmentation and structure fitting
CN109711416B (en) Target identification method and device, computer equipment and storage medium
WO2024125434A1 (en) Regional-consistency-based building principal angle correction method
CN113689445B (en) High-resolution remote sensing building extraction method combining semantic segmentation and edge detection
CN110619258A (en) Road track checking method based on high-resolution remote sensing image
CN115018249B (en) Subway station construction quality evaluation method based on laser scanning technology
WO2022110862A1 (en) Method and apparatus for constructing road direction arrow, electronic device, and storage medium
CN110175574A (en) A kind of Road network extraction method and device
CN116258722B (en) Intelligent bridge building detection method based on image processing
Zheng et al. Building recognition of UAV remote sensing images by deep learning
CN115457277A (en) Intelligent pavement disease identification and detection method and system
CN116152494A (en) Building foot point identification segmentation method based on two-stage 3D point cloud semantic segmentation
CN115222883A (en) Electric power tower reconstruction method based on foundation LiDAR point cloud
CN114549956A (en) Deep learning assisted inclined model building facade target recognition method
CN113409332B (en) Building plane segmentation method based on three-dimensional point cloud
CN112561989B (en) Recognition method for hoisting object in construction scene
CN113902792A (en) Building height detection method and system based on improved RetinaNet network and electronic equipment
CN116363319B (en) Modeling method, modeling device, equipment and medium for building roof
CN110889418A (en) Gas contour identification method
CN115661398A (en) Building extraction method, device and equipment for live-action three-dimensional model
CN114511582A (en) Automatic ancient city battlement extraction method
CN114638970A (en) Electric power facility standardization judgment system based on machine learning and edge detection
Song et al. A region-based approach to building detection in densely build-up high resolution satellite image
Zhang et al. Building extraction from high-resolution remote sensing images based on GrabCut with automatic selection of foreground and background samples
CN116452604B (en) Complex substation scene segmentation method, device and storage medium