CN113592877B - Method and device for identifying red line exceeding of pumped storage power station - Google Patents
Method and device for identifying red line exceeding of pumped storage power station Download PDFInfo
- Publication number
- CN113592877B CN113592877B CN202110318387.5A CN202110318387A CN113592877B CN 113592877 B CN113592877 B CN 113592877B CN 202110318387 A CN202110318387 A CN 202110318387A CN 113592877 B CN113592877 B CN 113592877B
- Authority
- CN
- China
- Prior art keywords
- boundary
- red line
- construction
- area
- boundary point
- 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 47
- 238000010276 construction Methods 0.000 claims abstract description 162
- 239000013598 vector Substances 0.000 claims description 80
- 230000008859 change Effects 0.000 claims description 54
- 230000011218 segmentation Effects 0.000 claims description 43
- 238000005516 engineering process Methods 0.000 claims description 19
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 5
- 230000000007 visual effect Effects 0.000 claims description 5
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 238000007499 fusion processing Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 abstract description 10
- 230000008569 process Effects 0.000 abstract description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 239000002689 soil Substances 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/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration by the use of histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Abstract
The application provides a method and a device for identifying red line exceeding of a pumped storage power station. The method comprises the following steps: acquiring a pre-construction remote sensing image and an in-construction remote sensing image corresponding to a pumped storage power station; determining a construction disturbance area of the pumped storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image; determining all boundary points exceeding the red line area in the construction disturbance area based on a pre-stored red line area corresponding to the pumped storage power station, and integrating all boundary points exceeding the red line area into a plane to obtain a red line exceeding area; the red line area is used for indicating a standard construction disturbance area of the pumped storage power station. The scheme ensures that the construction disturbance area and the red line exceeding area in the construction process of the pumped storage power station are extracted more accurately, and compared with manual measurement, the efficiency is higher, the cost is lower, and the labor cost and the measurement time are greatly saved.
Description
Technical Field
The application relates to the technical field of power stations, in particular to a pumped storage power station red line exceeding identification method and device.
Background
With the economic growth of China, the demand for electric power is higher and higher. In order to save resources, a pumped storage power station which uses water resources as power generation power is generated. The construction disturbance range of the construction project of the pumped storage power station is large, and the construction time is long. In recent years, as the construction progress is accelerated, the construction range of the pumped storage power station is enlarged, but the construction occupation range of the pumped storage power station has clear requirements in a water and soil conservation scheme report. Therefore, how to timely acquire the construction disturbance range of the pumped storage power station and accurately judge whether the construction disturbance range exceeds the red line range is an urgent problem to be solved.
The traditional red line exceeding area checking mode mainly relies on manual field measurement to obtain a construction disturbance range firstly, and then the construction disturbance range and the red line range are subjected to comparison analysis. The manual field measurement requires at least two people to hold the instrument for coordinate determination. Firstly, selecting a certain point as a measuring station, and setting three-dimensional coordinates of the measuring station; and then setting parameters such as rear view point coordinates, prism constants and the like, calibrating the target prism, and pressing a coordinate measurement key so that the total station starts to acquire the three-dimensional coordinates of the measuring point. Because the pumped storage power station occupies a large area, the construction disturbance range is measured manually, and the problems of more investment personnel, low working efficiency, high instrument cost, inaccurate data acquisition and the like exist, so that the engineering requirements cannot be met.
Disclosure of Invention
Aiming at the problems in the prior art, the embodiment of the application provides a red line exceeding identification method and device for a pumped storage power station, which solve the technical problems of low efficiency and low accuracy of a method for manually checking a red line exceeding area.
In one aspect, an embodiment of the present application provides a pumped storage power station red line exceeding identification method, where the method includes: acquiring a pre-construction remote sensing image and an in-construction remote sensing image corresponding to a pumped storage power station; determining a construction disturbance area of the pumped storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image; determining all boundary points exceeding the red line area in the construction disturbance area based on a pre-stored red line area corresponding to the pumped storage power station, and integrating all boundary points exceeding the red line area into a plane to obtain a red line exceeding area; the red line area is used for indicating a standard construction disturbance area of the pumped storage power station.
According to the method and the device, remote sensing images of the pumped storage power station under different construction states are obtained through a satellite remote sensing technology, construction disturbance areas in the remote sensing images are obtained through a GIS technology, finally, the construction disturbance areas are compared with specified red line areas through a preset algorithm, the portions of the construction disturbance areas exceeding the red line areas are calculated, detection of the red line exceeding areas can be automatically completed, site measurement is not needed, the red line exceeding areas can be rapidly checked, and the accuracy of determining the red line exceeding areas is improved on the premise of ensuring efficiency.
In one embodiment, the determining the construction disturbance area of the pumped-storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image specifically includes: performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image; acquiring a pixel value corresponding to the differential image, and determining a pixel value distribution histogram corresponding to the differential image based on the pixel value; determining a first segmentation threshold corresponding to the differential image through the pixel value distribution histogram; dividing the differential image into a first change region and a first non-change region based on the first segmentation threshold, calculating an average pixel value of the first change region and calculating an average pixel value of the first non-change region; determining a second segmentation threshold corresponding to the differential image according to the average pixel value of the first change region and the average pixel value of the first non-change region; dividing the differential image into a second change region and a second non-change region based on the second segmentation threshold; and determining a construction disturbance area of the pumped storage power station in the second change area through a Geographic Information System (GIS) technology and a man-machine interaction visual interpretation method.
According to the embodiment of the application, the differential images of the remote sensing image before construction and the remote sensing image in construction are obtained through an image differential technology, then the segmentation threshold value is determined according to the distribution condition of the pixel values in the check image, and the change area and the non-change area of the pumped storage power station are distinguished. In order to improve the dividing accuracy, the embodiment of the application determines a new dividing threshold according to the average pixel value of the changed area after the first division and the average pixel value of the unchanged area, and finally obtains the changed area of the pumped storage power station. The accuracy of the extraction of the change region is improved by selecting multiple thresholds. The differential image and threshold segmentation method can accurately extract the change region, and has high extraction accuracy and high speed.
In one embodiment, the dividing the differential image into the first variable region and the first non-variable region based on the first segmentation threshold specifically includes: determining a region with pixel values larger than the first segmentation threshold value in the differential image as the first change region; determining a region with the pixel value smaller than or equal to the first segmentation threshold value in the differential image as the first unchanged region; based on the second segmentation threshold, dividing the differential image into a second variation region and a second non-variation region, specifically including: determining a region in the differential image, the pixel value of which is greater than the second segmentation threshold value, as the second change region; and determining the region with the pixel value smaller than or equal to the second segmentation threshold value in the differential image as the second unchanged region.
In one embodiment, the performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image specifically includes: according toObtaining the tif of the remote sensing image before construction 1 Is the kth wave band of the image tif of the remote sensing image in construction 2 A difference image between the kth bands; wherein (1)>For the pre-construction remote sensing image tif 1 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band; />For the remote sensing image tif in the construction 2 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band; />Representative remote sensing image tif before construction 1 Remote sensing image tif in construction 2 And the absolute value of the difference value of the pixel values corresponding to the ith row and the jth column of pixel points in the kth wave band.
In one embodiment, the determining, based on the pre-stored red line area corresponding to the pumped storage power station, all the boundary points exceeding the red line area in the construction disturbance area, and integrating the boundary points exceeding the red line area into a plane, to obtain a red line exceeding area specifically includes: selecting a first boundary point on the construction disturbance area boundary line, and enabling the first boundary point, a second boundary point and a third boundary point on the red line area boundary line to respectively form a first vector and a second vector; wherein the second boundary point and the third boundary point are any group of adjacent boundary points on the boundary line of the red line area; calculating an included angle between the first vector and the second vector according to the coordinate values of the first boundary point, the second boundary point and the third boundary point; and traversing all adjacent boundary points on the boundary line of the red line area in a preset direction by taking the second boundary point as a starting point and the third boundary point as an ending point so as to obtain the sum of vector included angles formed between the first boundary point and all the adjacent boundary points on the boundary line of the red line area.
In one embodiment, the calculating the included angle between the first vector and the second vector according to the coordinate values of the first boundary point, the second boundary point and the third boundary point specifically includes: constructing a two-dimensional plane coordinate system by taking a first boundary line of the differential image as a horizontal axis, a second boundary line of the differential image as a vertical axis and an intersection point of the first boundary line and the second boundary line as an origin; in the two-dimensional plane coordinate system, coordinate values respectively corresponding to the first boundary point, the second boundary point and the third boundary point are determined; and determining the coordinates of the first vector and the coordinates of the second vector based on the coordinate values respectively corresponding to the first boundary point, the second boundary point and the third boundary point, so as to determine the included angle between the first vector and the second vector.
In one embodiment, after the obtaining the sum of vector angles formed between the first boundary point and all adjacent boundary points on the red line region boundary line, the method further includes: determining the position of the first boundary point relative to the red line region according to the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line region, wherein the method specifically comprises the following steps: under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 360 degrees, determining that the first boundary point is positioned in the red line area; and determining that the first boundary point is positioned outside the red line area under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 0 degree.
In one embodiment, after determining the position of the first boundary point relative to the red line region according to the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the red line region boundary line, the method further includes: traversing each boundary point on the construction disturbance range boundary line in a preset direction by taking the first boundary point as a starting point so as to determine all boundary points on the construction disturbance range boundary line, which are positioned outside the red line area; and converting all boundary points positioned outside the red line area into planes through a Geographic Information System (GIS) technology so as to obtain the red line exceeding area.
According to the embodiment of the application, the construction disturbance area and the red line area are mutually overlapped and are subjected to space analysis, and whether the boundary point on the construction disturbance area is positioned inside the red line area or outside the red line area is determined according to the sum of vector included angles formed by the boundary points on the construction disturbance area and all groups of adjacent boundary points on the red line area, so that an accurate red line exceeding area is obtained. The red line exceeding area has high extraction speed and high extraction accuracy, and does not need to be manually participated.
In one embodiment, after the acquiring the pre-construction and in-construction remote sensing images of the pumped-storage power plant, the method further comprises: and preprocessing the pre-construction remote sensing image and the in-construction remote sensing image at least by one or more of the following: radiation correction processing, atmospheric correction processing, geometric correction processing, and image fusion processing.
On the other hand, the embodiment of the application provides a pumped storage power station red line exceeds standard recognition device, the device includes: the acquisition unit is used for acquiring a pre-construction remote sensing image and an in-construction remote sensing image of the pumped storage power station; the determining unit is used for determining a construction disturbance area of the pumped storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image; the determining unit is further used for determining all boundary points exceeding the red line area in the construction disturbance area based on the pre-stored red line area corresponding to the pumped storage power station, and integrating all boundary points exceeding the red line area into a plane to obtain a red line exceeding area.
According to the method and the device, the differential technology is utilized to obtain the differential image of the remote sensing image before construction and the remote sensing image in construction of the pumped storage power station, then the differential image is subjected to threshold segmentation twice to obtain an accurate change area, and the construction disturbance area is extracted by utilizing the GIS technology, so that the extraction of the construction disturbance area is more accurate, and compared with manual measurement, the efficiency is higher, and the cost is lower. And finally, determining a red line exceeding area of the construction disturbance area exceeding the red line area through accurate calculation, so that the red line exceeding area in the construction process of the pumped storage power station is more convenient and accurate to extract, manual on-site measurement is not needed, the labor cost and the measurement time are greatly saved, and meanwhile, the accuracy of the determined red line exceeding area is ensured.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this application, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute an undue limitation to the application. In the drawings:
fig. 1 is a flowchart of a method for identifying a red line exceeding of a pumped storage power station according to an embodiment of the present application;
FIG. 2 is a flow chart of a second method for identifying the red line exceeding of the pumped storage power station according to the embodiment of the present application;
FIG. 3 is a flowchart of a third method for identifying a red line exceeding of a pumped storage power station according to an embodiment of the present application;
FIG. 4 is a histogram of pixel value distribution according to an embodiment of the present disclosure;
fig. 5 is a schematic structural diagram of a pumped storage power station red line exceeding identification device according to an embodiment of the present application.
Detailed Description
For the purposes, technical solutions and advantages of the present application, the technical solutions of the present application will be clearly and completely described below with reference to specific embodiments of the present application and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
The embodiment of the application provides a method and a device for identifying the red line exceeding of a pumped storage power station, and the technical scheme provided by the embodiment of the application is described in detail through the attached drawings.
Fig. 1 is a flowchart of a method for identifying red line exceeding of a pumped storage power station according to an embodiment of the present application, as shown in fig. 1:
s101, acquiring a pre-construction remote sensing image and an in-construction remote sensing image corresponding to the pumped storage power station.
Specifically, before the construction of the pumped storage power station, acquiring a remote sensing image tif before the construction of the pumped storage power station through a network 1 The method comprises the steps of carrying out a first treatment on the surface of the In the construction process of the pumped storage power station, remote sensing image tif in construction of the pumped storage power station is acquired through a network 2 . The technology of acquiring the remote sensing image through the network is a mature prior art, and this will not be described in detail in the embodiments of the present application.
Further, the acquired pre-construction remote sensing image and the pre-construction remote sensing image are subjected to preprocessing such as radiation correction, atmospheric correction, geometric correction, image fusion and the like, so that errors generated by the fact that the remote sensing image is subjected to sensor, landform, atmosphere and other factor images are eliminated, and the recognition and calculation accuracy of subsequent steps is improved.
S102, determining a construction disturbance area of the pumped storage power station based on the acquired pre-construction remote sensing image and the in-construction remote sensing image.
Specifically, as shown in fig. 2, S102 may specifically include the following steps:
s201, performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image.
Specifically, one remote sensing image comprises a plurality of wave bands, and the remote sensing image tif before construction 1 And remote sensing image tif in construction 2 And respectively carrying out image difference calculation on each wave band to obtain difference images of each wave band, and finally synthesizing the difference images of each wave band into a true color difference image.
The image difference calculation process is as follows:
according to the formulaObtaining the tif of the remote sensing image before construction 1 Is the kth wave band of the image tif of the remote sensing image in construction 2 A difference image between the kth bands;
wherein,for remote sensing image tif before construction 1 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band; />Is a remote sensing image tif in construction 2 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band; />Representative remote sensing image tif before construction 1 Remote sensing image tif in construction 2 And the absolute value of the difference value of the pixel values corresponding to the ith row and the jth column of pixel points in the kth wave band.
Calculating the tif of the remote sensing image before construction according to the formula 1 And remote sensing image tif in construction 2 Absolute values of differences of all corresponding pixel values in the kth wave band to obtain a remote sensing image tif before construction 1 And remote sensing image tif in construction 2 Differential imaging between the kth bands. And finally, synthesizing the obtained differential images of all the wave bands into a complete differential image.
Each pixel value in the differential image represents a pixel change value between the pre-construction remote sensing image and the in-construction remote sensing image, so that the change between the corresponding pixels of the pre-construction remote sensing image and the in-construction remote sensing image can be distinguished through the size of the pixel value in the differential image. Through the pixel value distribution condition in the differential image, the change of which pixels is larger and the change of which pixels is smaller can be clearly seen, so that the change area and the non-change area of the two images can be more conveniently analyzed.
S202, after the complete differential image is obtained, obtaining a pixel value corresponding to the differential image, and determining a pixel value distribution histogram corresponding to the differential image based on the pixel value.
Specifically, the pixel values in the obtained differential image are arranged into a pixel value distribution histogram, and the distribution situation of the pixel values in the differential image can be clearly seen through the pixel value distribution histogram.
In one embodiment, fig. 4 is an example of a pixel value distribution histogram, and as shown in fig. 4, the abscissa of the pixel value distribution histogram is a range bin of pixel values, and the ordinate is the number of pixels occurring in a certain range bin.
S203, determining a first segmentation threshold corresponding to the differential image through the pixel value distribution histogram, and dividing the differential image into a first change area and a first non-change area based on the first segmentation threshold.
Specifically, based on the obtained pixel value distribution histogram, a range section to which most of the pixel values belong is determined, and the right boundary value of the range section is set as the first division threshold. Dividing a region with pixel values larger than a first segmentation threshold value in the differential image into a first change region; and dividing the region with the pixel value smaller than or equal to the first segmentation threshold value in the differential image into a first unchanged region.
In one embodiment, as shown in fig. 4, most pixel values in the differential image are between [0-100], so that the right boundary value 100 of the interval can be set as the first division threshold, the pixels with pixel values greater than 100 are divided into the first variation regions, and the pixels with pixel values less than or equal to 100 are divided into the first non-variation regions.
S204, after the differential image is divided into a first change area and a first non-change area, calculating an average pixel value of the first change area and an average pixel value of the first non-change area, and determining a second segmentation threshold corresponding to the differential image according to the average pixel value of the first change area and the average pixel value of the first non-change area.
Specifically, the average pixel values P of the first variation regions obtained in S203 are calculated respectively 1 And an average pixel value P of the first non-changing region 2 . Then calculate P 1 And P 2 As the second segmentation threshold.
In one example, if the average pixel value of the first variation region obtained for the first time is 80 and the average pixel value of the first non-variation region is 180, the second segmentation threshold is (80+180)/2=130.
S205, dividing the differential image into a second change area and a second non-change area based on a second segmentation threshold, and determining a construction disturbance area of the pumped storage power station in the second change area through a GIS technology, a man-machine interaction technology and a visual interpretation method.
Specifically, dividing a region in the differential image, wherein the pixel value of the region is larger than a second segmentation threshold value, into a second change region; and dividing the region with the pixel value smaller than or equal to the second segmentation threshold value in the differential image into a second unchanged region.
Further, a construction disturbance area vector is extracted from the second change area through a GIS technology and a man-machine interaction visual interpretation method.
According to the embodiment of the application, the segmentation threshold is selected for multiple times, rather than only the segmentation threshold determined for the first time is adopted to segment the differential image, and therefore a more accurate segmentation result can be obtained.
And S103, after the construction disturbance area of the pumped storage power station is determined, determining a red line exceeding area of the construction disturbance area exceeding the red line area based on the red line area corresponding to the pre-stored pumped storage power station.
The construction disturbance area vector is superimposed with a red line area vector corresponding to a pre-stored pumped storage power station, space analysis is carried out on the two superimposed vectors, and an area, beyond the red line area, of the construction disturbance area, namely, a red line exceeding area is determined.
As shown in fig. 3, S103 may specifically include the following steps:
s301, selecting a first boundary point on a boundary line of a construction disturbance area, and enabling the first boundary point, a second boundary point and a third boundary point on a boundary line of a red line area to respectively form a first vector and a second vector.
Specifically, the boundary line of the construction disturbance area and the boundary line of the red line area are regarded as a closed line composed of innumerable continuous geometric points, such as a= { a 0 ，A 1 ，A 2 …，A n ，A 0 }. And selecting any point on the boundary line of the construction disturbance area as a first boundary point, and selecting any group of adjacent points on the boundary line of the red line area as a second boundary point and a third boundary point. The first boundary point forms a first vector and a second vector with the second boundary point and the third boundary point, respectively.
In one embodiment, specifically, any point B (X B ，Y B ) As the first boundary point, a second boundary point A arbitrarily adjacent to the red line region boundary line 0 (X A0 ，Y A0 ) Third boundary point A 1 (X A1 ，Y A1 ) Constitute a first vectorSecond vector->
S302, calculating an included angle between the first vector and the second vector according to coordinate values of the first boundary point, the second boundary point and the third boundary point.
Specifically, a two-dimensional planar coordinate system is constructed with the first boundary line of the differential image as the horizontal axis, the second boundary line of the differential image as the vertical axis, and the intersection point of the first boundary line and the second boundary line as the origin. For example, a two-dimensional planar coordinate system is constructed with the upper boundary line of the differential image as the horizontal axis, the left boundary line as the vertical axis, and the intersection point of the upper boundary line and the left boundary line as the origin.
After the two-dimensional plane coordinate system is constructed, coordinate values corresponding to the first boundary point, the second boundary point and the third boundary point are determined in the two-dimensional plane coordinate system, further, the coordinates of the first vector and the coordinates of the second vector are determined, and an included angle between the first vector and the second vector is calculated according to the coordinates of the first vector and the coordinates of the second vector.
In one embodiment, the first boundary point B (X B ，Y B ) Second boundary point A 0 (X A0 ，Y A0 ) Third boundary point A 1 (X A1 ，Y A1 ) Is a vector point multiplication formulaCalculating the first vector +.>And a second vector->Included angle delta between i ：
Since the construction disturbance area and the red line area are vectors, it is necessary to define a positive direction, for example, a clockwise direction or a counterclockwise direction. In the case that the first vector is located in the positive direction of the second vector, the included angle between the vectors is a positive number; in the case where the first vector is located in the negative direction of the second vector, the angle between the vectors is negative.
S303, traversing all adjacent boundary points on the boundary line of the red line area in a preset direction by taking the second boundary point as a starting point and the third boundary point as an ending point, so as to obtain the sum of vector included angles formed between the first boundary point and all the adjacent boundary points on the boundary line of the red line area.
Specifically, on the boundary line of the red line region, each group of adjacent boundary points is traversed according to a preset positive direction by taking the second boundary point as a starting point until the third boundary point is traversed. Each traversal calculates the vector angle formed between the first boundary point and two adjacent boundary points traversed once. And after the traversal is finished, summing vector included angles calculated in each traversal to obtain the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area.
In one embodiment, the boundary line A= { A is in the red line region 0 ，A 1 ，A 2 …，A n ，A 0 A on } 0 The point is a second boundary point, A is taken according to the positive direction 0 Adjacent point A of the points 1 Is the third boundary point. Calculating first boundary points B and A 0 、A 1 First vector formed between pointsAnd a second vector->Included angle delta between 1 . Then in A 1 The point is a second boundary point, A is taken according to the positive direction 1 Adjacent point A of the points 2 For the third boundary point, calculate the first boundary points B and A 1 、A 2 First vector between dots +.>And a second vector->Included angle delta between 2 . And so on until A is taken n As the second boundary point, A 0 For the third boundary point, calculate the first boundary points B and A n 、A 0 First vector between dots +.>And a second step ofVector->Included angle delta between n+1 The traversal is then stopped.
Finally according to the summation formulaThe sum delta of the angles of the vectors formed by the first boundary point B and each pair of adjacent points on the boundary of the red line area is obtained.
S304, determining the position of the first boundary point relative to the red line region according to the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line region.
Specifically, under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 360 degrees, the first boundary point is indicated to be positioned in the red line area; when the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 0 degrees, the first boundary point is positioned outside the red line area.
In one embodiment, when δ=360°, it is stated that the first boundary point B is located inside the red line region. When δ=0°, it is explained that the first boundary point B is located outside the red line region.
S305, traversing each boundary point on the construction disturbance range boundary line in a preset direction by taking the first boundary point as a starting point, so as to determine all boundary points on the construction disturbance range boundary line, which are positioned outside the red line area.
Specifically, using a first boundary point B as a starting point, traversing each boundary point on a construction disturbance range boundary line in a specified positive direction, and calculating the sum of the included angles of vectors formed by the boundary point and each group of adjacent points on a red line area boundary line respectively in each traversal; then, whether the boundary point is located inside or outside the red line region is determined by the method described in S304. And finally, determining and marking the set of all points on the construction disturbance range boundary line, which are positioned outside the red line area.
S306, converting the set of all boundary points positioned outside the red line area into a surface through a Geographic Information System (GIS) technology to obtain the red line exceeding area.
And (3) converting the set of all the points which are marked in the S305 and are positioned outside the red line area into planar vectors by a GIS technology, thereby obtaining the red line exceeding area. And the obtained red line exceeding area is stored and archived for convenient subsequent reading and analysis. It should be noted that, the existing GIS technology can efficiently convert the point set into one plane, which is not described in the embodiments of the present application.
The size of the red line exceeding area obtained through analysis can clearly judge whether the disturbance range in the construction process of the pumped storage power station exceeds the standard or not, and data support is provided for guaranteeing the construction safety of the pumped storage power station.
In addition, the embodiment of the application further provides a pumped storage power station red line exceeding identification device 500, as shown in fig. 5, the device includes:
an acquiring unit 501, configured to acquire a pre-construction remote sensing image and an in-construction remote sensing image of a pumped storage power station;
the determining unit 502 is configured to determine a construction disturbance area of the pumped storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image;
the determining unit 502 is further configured to determine, based on a red line area corresponding to the pre-stored pumped storage power station, a red line exceeding area where the construction disturbance area exceeds the red line area.
According to the method and the device for identifying the red line exceeding of the pumped storage power station, on one hand, the differential remote sensing images before and during construction are obtained through the image differential technology, and then the construction disturbance area is accurately extracted from the differential remote sensing images by utilizing the multiple threshold segmentation and the GIS technology. Compared with the manual measurement of the construction disturbance area, the machine operation speed is faster and the efficiency is higher; multiple thresholding results in higher accuracy in extracting the construction disturbance zone. After the construction disturbance area is extracted, the embodiment of the application adopts a vector included angle calculation method to skillfully determine whether the point on the construction disturbance area exceeds the red line area, and the calculation process is realized through an algorithm, so that the method is more accurate and simpler than manual comparison, the time for extracting the red line exceeding area is also greatly shortened, and the accuracy of the extracted red line exceeding area is ensured.
The foregoing is merely a preferred embodiment of the present application and is not intended to limit the present application. It is noted that various modifications and variations can be made by persons skilled in the art. Any modification, equivalent replacement, improvement, etc. made without departing from the principles of the present application should be included in the scope of the claims of the present application.
Claims (6)
1. The method for identifying the red line exceeding of the pumped storage power station is characterized by comprising the following steps of:
acquiring a pre-construction remote sensing image and an in-construction remote sensing image corresponding to a pumped storage power station;
based on the pre-construction remote sensing image and the in-construction remote sensing image, determining a construction disturbance area of the pumped storage power station specifically comprises the following steps:
performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image;
acquiring a pixel value distribution histogram corresponding to the differential image;
determining a first segmentation threshold corresponding to the differential image through the pixel value distribution histogram;
based on the first segmentation threshold, the differential image is divided into a first change region and a first non-change region: determining a region with pixel values larger than the first segmentation threshold value in the differential image as the first change region; determining a region with the pixel value smaller than or equal to the first segmentation threshold value in the differential image as the first unchanged region;
calculating an average pixel value of the first variation region, and calculating an average pixel value of the first non-variation region;
determining a second segmentation threshold corresponding to the differential image according to the average pixel value of the first change region and the average pixel value of the first non-change region;
based on the second segmentation threshold, dividing the differential image into a second change region and a second non-change region: determining a region in the differential image, the pixel value of which is greater than the second segmentation threshold value, as the second change region; determining that a region in the differential image, in which the pixel value is smaller than or equal to the second segmentation threshold value, is the second unchanged region;
determining a construction disturbance area of the pumped storage power station in the second change area through a Geographic Information System (GIS) technology and a man-machine interaction visual interpretation method;
determining all boundary points exceeding the red line area in the construction disturbance area based on a pre-stored red line area corresponding to the pumped storage power station, and integrating all boundary points exceeding the red line area into a plane to obtain a red line exceeding area; the red line area is used for indicating a standard construction disturbance area of the pumped storage power station;
the method specifically comprises the following steps: selecting a first boundary point on the construction disturbance area boundary line, and enabling the first boundary point, a second boundary point and a third boundary point on the red line area boundary line to respectively form a first vector and a second vector; wherein the second boundary point and the third boundary point are any group of adjacent boundary points on the boundary line of the red line area;
calculating an included angle between the first vector and the second vector according to the coordinate values of the first boundary point, the second boundary point and the third boundary point;
traversing all adjacent boundary points on the red line area boundary line in a preset direction by taking the second boundary point as a starting point and the third boundary point as an ending point so as to obtain the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the red line area boundary line;
determining the position of the first boundary point relative to the red line region according to the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line region, wherein the method specifically comprises the following steps:
under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 360 degrees, determining that the first boundary point is positioned in the red line area;
and determining that the first boundary point is positioned outside the red line area under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 0 degree.
2. The method for identifying the red line exceeding of the pumped storage power station according to claim 1, wherein the performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image specifically comprises:
according toObtaining the tif of the remote sensing image before construction 1 Is the kth wave band of the image tif of the remote sensing image in construction 2 A difference image between the kth bands;
wherein,for the pre-construction remote sensing image tif 1 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band; />For the remote sensing image tif in the construction 2 The pixel value corresponding to the ith row and the jth column of pixel points in the kth wave band;representative remote sensing image tif before construction 1 Remote sensing image tif in construction 2 And the absolute value of the difference value of the pixel values corresponding to the ith row and the jth column of pixel points in the kth wave band.
3. The method for identifying the red line exceeding of the pumped storage power station according to claim 1, wherein the calculating the included angle between the first vector and the second vector according to the coordinate values of the first boundary point, the second boundary point and the third boundary point specifically comprises:
constructing a two-dimensional plane coordinate system by taking a first boundary line of the differential image as a horizontal axis, a second boundary line of the differential image as a vertical axis and an intersection point of the first boundary line and the second boundary line as an origin;
in the two-dimensional plane coordinate system, coordinate values respectively corresponding to the first boundary point, the second boundary point and the third boundary point are determined;
and determining the coordinates of the first vector and the coordinates of the second vector based on the coordinate values respectively corresponding to the first boundary point, the second boundary point and the third boundary point, so as to determine the included angle between the first vector and the second vector.
4. The method according to claim 1, wherein after determining the position of the first boundary point relative to the red line region according to the sum of vector angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line region, the method further comprises:
traversing each boundary point on the construction disturbance range boundary line in a preset direction by taking the first boundary point as a starting point so as to determine all boundary points on the construction disturbance range boundary line, which are positioned outside the red line area;
and converting all boundary points positioned outside the red line area into planes through a Geographic Information System (GIS) technology so as to obtain the red line exceeding area.
5. The method for identifying red line exceeding of a pumped-hydro power station as defined in claim 1, wherein after the acquiring the pre-construction and in-construction remote sensing images of the pumped-hydro power station, the method further comprises:
and preprocessing the pre-construction remote sensing image and the in-construction remote sensing image at least by one or more of the following: radiation correction processing, atmospheric correction processing, geometric correction processing, and image fusion processing.
6. The utility model provides a pumped storage power station red line exceeds standard recognition device which characterized in that includes:
the acquisition unit is used for acquiring a pre-construction remote sensing image and an in-construction remote sensing image of the pumped storage power station;
the determining unit is used for determining a construction disturbance area of the pumped storage power station based on the pre-construction remote sensing image and the in-construction remote sensing image, and specifically comprises the following steps:
performing image difference calculation on the pre-construction remote sensing image and the in-construction remote sensing image to obtain a difference image;
acquiring a pixel value distribution histogram corresponding to the differential image;
determining a first segmentation threshold corresponding to the differential image through the pixel value distribution histogram;
based on the first segmentation threshold, the differential image is divided into a first change region and a first non-change region: determining a region with pixel values larger than the first segmentation threshold value in the differential image as the first change region; determining a region with the pixel value smaller than or equal to the first segmentation threshold value in the differential image as the first unchanged region;
calculating an average pixel value of the first variation region, and calculating an average pixel value of the first non-variation region;
determining a second segmentation threshold corresponding to the differential image according to the average pixel value of the first change region and the average pixel value of the first non-change region;
based on the second segmentation threshold, dividing the differential image into a second change region and a second non-change region: determining a region in the differential image, the pixel value of which is greater than the second segmentation threshold value, as the second change region; determining that a region in the differential image, in which the pixel value is smaller than or equal to the second segmentation threshold value, is the second unchanged region;
determining a construction disturbance area of the pumped storage power station in the second change area through a Geographic Information System (GIS) technology and a man-machine interaction visual interpretation method;
the determining unit is further used for determining all boundary points exceeding the red line area in the construction disturbance area based on a pre-stored red line area corresponding to the pumped storage power station, and integrating all boundary points exceeding the red line area into a plane to obtain a red line exceeding area; the method specifically comprises the following steps: selecting a first boundary point on the construction disturbance area boundary line, and enabling the first boundary point, a second boundary point and a third boundary point on the red line area boundary line to respectively form a first vector and a second vector; wherein the second boundary point and the third boundary point are any group of adjacent boundary points on the boundary line of the red line area;
calculating an included angle between the first vector and the second vector according to the coordinate values of the first boundary point, the second boundary point and the third boundary point;
traversing all adjacent boundary points on the red line area boundary line in a preset direction by taking the second boundary point as a starting point and the third boundary point as an ending point so as to obtain the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the red line area boundary line;
determining the position of the first boundary point relative to the red line region according to the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line region, wherein the method specifically comprises the following steps:
under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 360 degrees, determining that the first boundary point is positioned in the red line area;
and determining that the first boundary point is positioned outside the red line area under the condition that the sum of vector included angles formed between the first boundary point and all adjacent boundary points on the boundary line of the red line area is 0 degree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110318387.5A CN113592877B (en) | 2021-03-25 | 2021-03-25 | Method and device for identifying red line exceeding of pumped storage power station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110318387.5A CN113592877B (en) | 2021-03-25 | 2021-03-25 | Method and device for identifying red line exceeding of pumped storage power station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113592877A CN113592877A (en) | 2021-11-02 |
| CN113592877B true CN113592877B (en) | 2024-04-12 |
Family
ID=78238141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110318387.5A Active CN113592877B (en) | 2021-03-25 | 2021-03-25 | Method and device for identifying red line exceeding of pumped storage power station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113592877B (en) |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976437A (en) * | 2010-09-29 | 2011-02-16 | 中国资源卫星应用中心 | High-resolution remote sensing image variation detection method based on self-adaptive threshold division |
| CN103473764A (en) * | 2013-09-03 | 2013-12-25 | 中国科学院深圳先进技术研究院 | Remote-sensing image target change detection method |
| KR101380329B1 (en) * | 2013-02-08 | 2014-04-02 | (주)나노디지텍 | Method for detecting change of image |
| KR101490769B1 (en) * | 2013-11-15 | 2015-02-09 | 동명대학교산학협력단 | Method and apparatus of detecting fire using brightness and area variation |
| CN207149032U (en) * | 2017-03-15 | 2018-03-27 | 国家电网公司 | A kind of construction disturbance area remote supervision system based on satellite remote sensing images |
| CN108805920A (en) * | 2018-06-01 | 2018-11-13 | 环境保护部卫星环境应用中心 | The recognition methods in soil pollution risk region and device on unused land |
| CN108830828A (en) * | 2018-04-28 | 2018-11-16 | 新疆大学 | A kind of method for detecting change of remote sensing image and device |
| CN110411370A (en) * | 2019-07-26 | 2019-11-05 | 北京住总集团有限责任公司 | A kind of the risk of tunnel construction managing and control system based on space time parameter |
| CN110432832A (en) * | 2019-07-03 | 2019-11-12 | 平安科技（深圳）有限公司 | Method of adjustment, device and the robot of robot motion track |
| CN110675350A (en) * | 2019-10-22 | 2020-01-10 | 普联技术有限公司 | Cloud deck camera view field coordinate mapping method and device, storage medium and cloud deck camera |
| CN111626131A (en) * | 2020-04-27 | 2020-09-04 | 国家电网有限公司 | Electric power construction area information extraction and change detection method based on unmanned aerial vehicle remote sensing |
| CN111898494A (en) * | 2020-07-16 | 2020-11-06 | 大同煤矿集团有限责任公司 | Mining disturbed block boundary identification method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8160296B2 (en) * | 2005-04-15 | 2012-04-17 | Mississippi State University Research And Technology Corporation | Change analyst |
| CN105426818B (en) * | 2015-10-30 | 2019-07-02 | 小米科技有限责任公司 | Method for extracting region and device |
-
2021
- 2021-03-25 CN CN202110318387.5A patent/CN113592877B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101976437A (en) * | 2010-09-29 | 2011-02-16 | 中国资源卫星应用中心 | High-resolution remote sensing image variation detection method based on self-adaptive threshold division |
| KR101380329B1 (en) * | 2013-02-08 | 2014-04-02 | (주)나노디지텍 | Method for detecting change of image |
| CN103473764A (en) * | 2013-09-03 | 2013-12-25 | 中国科学院深圳先进技术研究院 | Remote-sensing image target change detection method |
| KR101490769B1 (en) * | 2013-11-15 | 2015-02-09 | 동명대학교산학협력단 | Method and apparatus of detecting fire using brightness and area variation |
| CN207149032U (en) * | 2017-03-15 | 2018-03-27 | 国家电网公司 | A kind of construction disturbance area remote supervision system based on satellite remote sensing images |
| CN108830828A (en) * | 2018-04-28 | 2018-11-16 | 新疆大学 | A kind of method for detecting change of remote sensing image and device |
| CN108805920A (en) * | 2018-06-01 | 2018-11-13 | 环境保护部卫星环境应用中心 | The recognition methods in soil pollution risk region and device on unused land |
| CN110432832A (en) * | 2019-07-03 | 2019-11-12 | 平安科技（深圳）有限公司 | Method of adjustment, device and the robot of robot motion track |
| CN110411370A (en) * | 2019-07-26 | 2019-11-05 | 北京住总集团有限责任公司 | A kind of the risk of tunnel construction managing and control system based on space time parameter |
| CN110675350A (en) * | 2019-10-22 | 2020-01-10 | 普联技术有限公司 | Cloud deck camera view field coordinate mapping method and device, storage medium and cloud deck camera |
| CN111626131A (en) * | 2020-04-27 | 2020-09-04 | 国家电网有限公司 | Electric power construction area information extraction and change detection method based on unmanned aerial vehicle remote sensing |
| CN111898494A (en) * | 2020-07-16 | 2020-11-06 | 大同煤矿集团有限责任公司 | Mining disturbed block boundary identification method |
Non-Patent Citations (4)
| Title |
|---|
| 《中国水土保持》;郑文革;北京市生产建设项目水土保持遥感监管创新与实践;第2.2.5部分 * |
| 《绿色科技》;张琳琳;无人机遥感在生产建设项目水土保持监测数据获取中的应用;第3.1.1部分 * |
| 张峰水库输水工程建设对沿线植被指数的影响;范彩娟;关志刚;;山西水土保持科技(第02期);全文 * |
| 无人机低空遥感技术在生产建设项目水土保持监测中的应用――以某燃机电厂为例;刘灿;杨浩翔;王鹏;张洋;;亚热带水土保持(第01期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113592877A (en) | 2021-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11551341B2 (en) | Method and device for automatically drawing structural cracks and precisely measuring widths thereof | |
| CN107092877B (en) | Remote sensing image roof contour extraction method based on building base vector | |
| CN109886939B (en) | Bridge crack detection method based on tensor voting | |
| CN107248159A (en) | A kind of metal works defect inspection method based on binocular vision | |
| CN112906694A (en) | Reading correction system and method for inclined pointer instrument image of transformer substation | |
| CN100468457C (en) | Method for matching depth image | |
| CN111336990A (en) | Tunnel section convergence rapid analysis method and device | |
| CN110619258B (en) | Road track checking method based on high-resolution remote sensing image | |
| CN110008207B (en) | Airborne L iDAR point cloud data vulnerability rapid detection method based on density histogram | |
| CN103295232A (en) | SAR (specific absorption rate) image registration method based on straight lines and area | |
| CN112489099A (en) | Point cloud registration method and device, storage medium and electronic equipment | |
| CN115993096A (en) | High-rise building deformation measuring method | |
| CN111354083A (en) | Progressive building extraction method based on original laser point cloud | |
| CN113487722A (en) | Automatic concrete member detection method based on three-dimensional laser scanning method | |
| Li et al. | A deep learning-based indoor acceptance system for assessment on flatness and verticality quality of concrete surfaces | |
| CN114913134A (en) | Tunnel shotcrete roughness identification method, terminal device and storage medium | |
| CN113902792A (en) | Building height detection method and system based on improved RetinaNet network and electronic equipment | |
| CN104596484A (en) | Method of measuring drift ice density in ice flood season of Yellow River | |
| CN113420670A (en) | Environment-friendly supervision method for changing power transmission and transformation line migration based on high-resolution remote sensing | |
| CN113673011A (en) | Method for intelligently identifying tunnel invasion boundary in operation period based on point cloud data | |
| CN113592877B (en) | Method and device for identifying red line exceeding of pumped storage power station | |
| CN113280764A (en) | Power transmission and transformation project disturbance range quantitative monitoring method and system based on multi-satellite cooperation technology | |
| CN104392209A (en) | Evaluation model for image complexity of target and background | |
| CN111336991A (en) | Tunnel ellipticity analysis method, device and system based on laser point cloud | |
| CN112284287B (en) | Stereoscopic vision three-dimensional displacement measurement method based on structural surface gray scale characteristics |
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 |