[go: up one dir, main page]

CN118865265B - A method and system for quality control of side-scan flow radar data - Google Patents

A method and system for quality control of side-scan flow radar data Download PDF

Info

Publication number
CN118865265B
CN118865265B CN202411319934.1A CN202411319934A CN118865265B CN 118865265 B CN118865265 B CN 118865265B CN 202411319934 A CN202411319934 A CN 202411319934A CN 118865265 B CN118865265 B CN 118865265B
Authority
CN
China
Prior art keywords
data
area
vortex
monitoring area
point cloud
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
Application number
CN202411319934.1A
Other languages
Chinese (zh)
Other versions
CN118865265A (en
Inventor
黄豫源
孙益阔
徐争
黄磊
欧钰婷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Wangshe Technology Industry Co ltd
Original Assignee
Nanjing Wangshe Technology Industry Co ltd
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 Nanjing Wangshe Technology Industry Co ltd filed Critical Nanjing Wangshe Technology Industry Co ltd
Priority to CN202411319934.1A priority Critical patent/CN118865265B/en
Publication of CN118865265A publication Critical patent/CN118865265A/en
Application granted granted Critical
Publication of CN118865265B publication Critical patent/CN118865265B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01PMEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
    • G01P5/00Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering
    • 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/50Context or environment of the image
    • G06V20/52Surveillance or monitoring of activities, e.g. for recognising suspicious objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/64Three-dimensional objects
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computing Systems (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Databases & Information Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供一种侧扫测流雷达数据质量控制方法及系统,其中,所述方法包括以下步骤:对待监测水域的每一监测区域进行数据采集,得到当前点云数据、当前图像数据以及当前流速数据;基于当前点云数据、当前图像数据确定对应监测区域的水域的水域状态,并确定对应处于异常状态下的异常时间;进行与异常时间具有相同时长的数据采集,得补充点云数据、补充图像数据以及补充流速数据;响应于基于补充点云数据、补充图像数据确定对应监测区域处于正常状态,将补充流速数据对对应所述异常时间的异常流速部分进行替换;基于各当前流速数据确定与所述待监测水域对应的水域流速。本发明至少提高了对于水域流速的测量精准性。

The present invention provides a side-scanning current radar data quality control method and system, wherein the method comprises the following steps: data collection is performed on each monitoring area of the monitored water area to obtain current point cloud data, current image data and current flow rate data; based on the current point cloud data and current image data, the water state of the water area corresponding to the monitored area is determined, and the abnormal time corresponding to the abnormal state is determined; data collection with the same duration as the abnormal time is performed to obtain supplementary point cloud data, supplementary image data and supplementary flow rate data; in response to determining that the corresponding monitoring area is in a normal state based on the supplementary point cloud data and supplementary image data, the abnormal flow rate portion corresponding to the abnormal time is replaced by the supplementary flow rate data; based on each current flow rate data, the water flow rate corresponding to the monitored water area is determined. The present invention at least improves the measurement accuracy of the water flow rate.

Description

Data quality control method and system for side-sweep flow-measuring radar
Technical Field
The invention relates to a data processing technology, in particular to a data quality control method and system for a side-sweep flow radar.
Background
The water flow rate is generally referred to as the velocity of the water flow and can be expressed in different units, such as meters per second (m/s) or centimeters per second (cm/s). The flow rate may vary depending on the nature of the body of water flow and the environmental conditions, such as the flow rate in rivers, streams, oceans.
In the prior art, a professional scanning radar is often used for measuring the flow rate of the water area to finish the corresponding measuring process, but the water area is always in different water area states due to different natural factors or human factors, for example, when a large vortex area, a large ripple area and a large noise area exist in the water area, the measuring accuracy of the flow rate of the water area can be influenced, the accurate measurement of the flow rate of the water area can not be realized, and the corresponding measuring efficiency is reduced.
Disclosure of Invention
The present invention has been made in view of the above problems, and it is an object of the present invention to provide a method and system for controlling data quality of a side-scan flow radar that overcomes or at least partially solves the above problems.
According to one aspect of the present invention, there is provided a method for controlling the quality of data of a side-scan flow radar, comprising the steps of:
carrying out grating treatment on a water area to be monitored, and respectively determining the obtained grating point positions as monitoring areas;
data acquisition is carried out on each monitoring area for preset time, and current point cloud data, current image data and current flow rate data of the corresponding monitoring area are obtained;
determining a water area state of a water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determining an abnormal time corresponding to the abnormal state;
Performing data acquisition of the supplementary time with the same duration as the abnormal time to obtain supplementary point cloud data, supplementary image data and supplementary flow rate data corresponding to the monitoring area;
In response to determining that the water area corresponding to the monitoring area is in a normal state within the replenishment time based on the replenishment point cloud data and the replenishment image data, replacing an abnormal flow rate portion corresponding to the abnormal time in the current flow rate data with the replenishment flow rate data;
And determining the water area flow rate corresponding to the water area to be monitored based on the current flow rate data corresponding to each monitoring area.
Optionally, in the method according to the present invention, data acquisition is performed for a preset time for each monitoring area to obtain current point cloud data, current image data, and current flow rate data of the corresponding monitoring area, including:
Triggering a point cloud radar, an image radar and a flow rate radar to acquire data for maintaining preset time in each monitoring area based on preset acquisition frequency;
Determining acquisition times based on the preset acquisition frequency and the preset time, and acquiring secondary point cloud data, secondary image data and secondary flow rate data corresponding to each acquisition time respectively;
Grouping each secondary point cloud data, each secondary image data and each secondary flow rate data based on each time node corresponding to each acquisition time respectively to obtain each node data group of the secondary point cloud data, the secondary image data and the secondary flow rate data corresponding to the same time node respectively;
and determining each node data set as current point cloud data, current image data and current flow rate data corresponding to the monitoring area.
Optionally, in the method according to the present invention, determining a water area state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determining an abnormal time corresponding to the abnormal state includes:
Establishing a two-dimensional plane based on a monitoring area, and establishing point clouds of each area corresponding to each point cloud height on the two-dimensional plane based on secondary point cloud data in any node data set;
Acquiring the area height corresponding to the monitoring area, respectively comparing the cloud height of each point with the area height, and determining a first morphological attribute of the monitoring area under a time node corresponding to the node data set based on a comparison result;
retrieving secondary image data in the same node data set as the secondary point cloud data, and determining a second shape attribute of the monitoring area under a time node corresponding to the node data set based on the secondary image data;
Determining a water area state of the monitoring area under the time node based on the first morphological attribute and the second morphological attribute corresponding to the monitoring area;
And summing all the time nodes corresponding to the abnormal state to obtain the abnormal time corresponding to the abnormal state.
Optionally, in the method according to the present invention, determining a first morphological attribute of the monitoring area under a time node corresponding to the node data set based on a comparison result includes:
Determining all the regional point clouds with the corresponding point cloud heights smaller than the regional heights as vortex point cloud groups according to the comparison result, and determining all the regional point clouds with the corresponding point cloud heights larger than the regional heights as ripple point cloud groups according to the comparison result;
performing point cloud connection on each area point cloud in adjacent position relation in the vortex point cloud group in the two-dimensional plane to obtain each vortex area in the two-dimensional plane;
performing point cloud connection on each area point cloud in adjacent position relation in the corrugated point cloud group in the two-dimensional plane to obtain each corrugated area in the two-dimensional plane;
Carrying out merging processing based on surrounding relation on each vortex area based on a preset merging strategy;
After the merging processing, determining the first area size and the first area number corresponding to each vortex area, and the second area size and the second area number corresponding to each ripple area as first morphological attributes of the monitoring area under the time node corresponding to the node data set.
Optionally, in the method according to the present invention, the merging process based on the surrounding relation is performed on each vortex area based on a preset merging policy, including:
carrying out coordinated processing on the two-dimensional plane to obtain vortex coordinate points respectively positioned in each vortex area, and obtaining vortex coordinate groups respectively corresponding to the vortex areas;
And comparing the vortex coordinate sets in pairs, determining whether the vortex coordinate sets have a surrounding relation or not based on a comparison result, and merging the vortex areas corresponding to the vortex coordinate sets with the surrounding relation into one vortex area.
Optionally, in the method according to the present invention, determining a second shape attribute of the monitoring area under a time node corresponding to the node data set based on the secondary image data comprises:
Performing binarization processing on the secondary image data to obtain a binarized image, wherein the binarized image comprises all water area pixel points corresponding to first pixel values and all noise pixel points corresponding to second pixel values respectively;
Pixel point connection is carried out on each noise pixel point in adjacent position relation in the secondary image data, so that each noise region in the secondary image data is obtained;
And determining the size and the number of the third areas corresponding to the noise areas as second shape attributes of the monitoring areas under the time nodes corresponding to the node data sets.
Optionally, in the method according to the present invention, determining the water status of the monitoring area under the time node based on the first morphological attribute and the second morphological attribute corresponding to the monitoring area includes:
summing calculation is carried out based on the sizes of all the first areas to obtain a first size total value, and ratio calculation is carried out based on the first size total value and the monitoring size of the corresponding monitoring area to obtain a first size ratio;
Respectively carrying out normalization processing on the first size ratio and the first area number to obtain a first size coefficient and a first quantity coefficient, and carrying out summation calculation on the basis of the first size coefficient and the first quantity coefficient to obtain a vortex coefficient;
Summing calculation is carried out based on the sizes of the second areas to obtain a second size total value, and ratio calculation is carried out based on the second size total value and the monitoring sizes of the corresponding monitoring areas to obtain a second size ratio;
respectively carrying out normalization processing on the second size ratio and the second area number to obtain a second size coefficient and a second number coefficient, and carrying out summation calculation on the basis of the second size coefficient and the second number coefficient to obtain a ripple coefficient;
summing calculation is carried out based on the sizes of all the third areas to obtain a third size total value, and ratio calculation is carried out based on the third size total value and the monitoring size of the corresponding monitoring area to obtain a third size ratio;
Respectively carrying out normalization processing on the third size ratio and the third area number to obtain a third size coefficient and a third number coefficient, and carrying out summation calculation based on the third size coefficient and the third number coefficient to obtain a noise coefficient;
And determining the water area state of the monitoring area under the time node based on the vortex coefficient, the ripple coefficient and the noise coefficient.
Optionally, in the method according to the present invention, determining the water state of the monitoring area under the time node based on the swirl coefficient, the ripple coefficient and the noise coefficient includes:
invoking vortex weight, ripple weight and noise weight;
performing product calculation on the vortex coefficient and the vortex weight to obtain a vortex evaluation value;
performing product calculation on the ripple coefficient and the ripple weight to obtain a ripple evaluation value;
performing product calculation on the noise coefficient and the noise weight to obtain a noise evaluation value;
When the vortex evaluation value is larger than a first preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
when the ripple evaluation value is larger than a second preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
When the noise evaluation value is larger than a third preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
When the vortex evaluation value is smaller than or equal to a first preset value, the ripple evaluation value is smaller than or equal to a second preset value and the noise evaluation value is smaller than or equal to a third preset value, summing the vortex evaluation value, the ripple evaluation value and the noise evaluation value to obtain a comprehensive evaluation value;
And when the comprehensive evaluation value is larger than a fourth preset value, determining that the water area state of the monitoring area under the time node is an abnormal state.
Optionally, in the method according to the present invention, the acquiring management end obtains the actual determining data based on the state of any node corresponding to any time sent by any monitoring area, and compares the water area state corresponding to the monitoring area with the actual determining data of the state;
If the state actually determines that the data is in a normal state and the water area state is in an abnormal state, performing reduction training on the vortex weight, the ripple weight and the noise weight;
if the state actually determines that the data is in an abnormal state and the water area state is in a normal state, training the vortex weight, the ripple weight and the noise weight in an increasing way;
The trained vortex weight, ripple weight, and noise weight are obtained by the following formula:
Wherein, Weighting eddiesThe number of times of training is increased,Weighting eddiesIs used for the training of the constant value of (a),Weighting eddiesThe number of times of training is reduced and,In order for the trained vortex weights to be applied,Is the ripple weightThe number of times of training is increased,Is the ripple weightIs used for the training of the constant value of (a),Is weighted by ripple coefficientThe number of times of training is reduced and,For the trained ripple weight,Weighting noiseThe number of times of training is increased,Weighting noiseIs used for the training of the constant value of (a),Weighting noiseThe number of times of training is reduced and,Is the weight of the vortex after training.
According to yet another aspect of the present invention, there is provided a side sweep radar data quality control system comprising:
The grid processing module is configured to perform grid processing on a water area to be monitored, and each obtained grid point position is respectively determined to be each monitoring area;
the data acquisition module is configured to acquire data of each monitoring area for a preset time to obtain current point cloud data, current image data and current flow rate data of the corresponding monitoring area;
The state determining module is configured to determine a water area state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determine an abnormal time corresponding to the abnormal state;
The supplementary acquisition module is configured to acquire data of supplementary time with the same duration as the abnormal time, and obtain supplementary point cloud data, supplementary image data and supplementary flow rate data corresponding to the monitoring area;
A data replacement module configured to replace the supplementary flow rate data with an abnormal flow rate portion corresponding to the abnormal time in the current flow rate data in response to determining that the water area corresponding to the monitoring area is in a normal state for the supplementary time based on the supplementary point cloud data, the supplementary image data;
And the flow rate determining module is configured to determine the water area flow rate corresponding to the water area to be monitored based on each current flow rate data corresponding to each monitoring area.
According to the scheme of the invention, through carrying out grating treatment on the water area to be monitored, respectively carrying out corresponding data acquisition based on the obtained monitoring areas, carrying out statistical calculation on the current flow velocity data corresponding to each monitoring area, obtaining the water area flow velocity corresponding to the water area to be monitored, finely realizing flow velocity measurement of the water area to be monitored, improving corresponding measurement accuracy, and simultaneously, when carrying out flow velocity measurement on each monitoring area, carrying out corresponding determination on the water area state of each monitoring area based on vortex dimension, ripple dimension and noise dimension, and when the water area state of the corresponding monitoring area is in an abnormal state, carrying out corresponding supplementary acquisition on the monitoring area, thereby eliminating interference on flow velocity measurement due to natural factors or artificial factors and further improving corresponding measurement accuracy.
Drawings
FIG. 1 illustrates a flow chart of a method of side sweep radar data quality control according to one embodiment of the present invention;
Fig. 2 shows a schematic diagram of flow rate measurement of a monitored area by the scanning radar in the present embodiment;
fig. 3 shows a system block diagram of a side-scan flow radar data quality control system according to another embodiment of the present invention.
Detailed Description
Exemplary embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
The water flow rate is generally referred to as the velocity of the water flow and can be expressed in different units, such as meters per second (m/s) or centimeters per second (cm/s). The flow rate may vary depending on the nature of the body of water flow and the environmental conditions, such as the flow rate in rivers, streams, oceans.
In the prior art, a professional scanning radar is often used for measuring the flow rate of the water area to finish the corresponding measuring process, but the water area is always in different water area states due to different natural factors or human factors, for example, when a large vortex area, a large ripple area and a large noise area exist in the water area, the measuring accuracy of the flow rate of the water area can be influenced, the accurate measurement of the flow rate of the water area can not be realized, and the corresponding measuring efficiency is reduced.
The proposal of the invention is provided for solving the problems in the prior art. One embodiment of the present invention provides a method for controlling the quality of data of a side-scan flow radar, which can be performed in a computing device, wherein the computing device can be understood as a terminal having a data processing function, such as a mobile phone or a computer.
As shown in fig. 1, the objective of the present embodiment is to implement a method for controlling the quality of data of a side-scan radar, which starts in step S102, and includes the following steps:
And carrying out rasterization treatment on the water area to be monitored, and respectively determining each obtained grille area as each monitoring area.
For example, in this embodiment, the flow rate measurement of the water area is generally performed by using a corresponding scanning radar, as shown in fig. 2, the scanning radar forms a corresponding scanning surface when scanning, and the scanning radar is controlled to move the scanning surface to the water surface of the water area to be monitored, so as to perform the flow rate scanning for a preset time on the water surface area overlapped with the scanning surface, so as to determine the water area flow rate of the water surface area; in a real situation, since the water area to be monitored generally has a larger water area, that is, the water area may be far larger than the area corresponding to the scanning surface, the water area to be monitored is subjected to corresponding grating treatment to obtain a plurality of corresponding grating areas, and the plurality of grating areas are respectively determined as each monitoring area, so that in a subsequent step, corresponding flow velocity measurement can be performed based on each monitoring area, and then the water area flow velocity corresponding to the water area to be monitored is determined according to the measurement result of each monitoring area, thereby improving the corresponding measurement accuracy.
In step S104, the following are included:
And acquiring data of each monitoring area for a preset time to obtain current point cloud data, current image data and current flow rate data of the corresponding monitoring area.
It should be noted that, in the embodiment, when the scanning radar is used to measure the flow rate of the water area to be monitored, the bad water surface condition of the water area to be monitored may affect the corresponding measurement accuracy, for example, when the corresponding large ripple, large vortex and large floater exist on the water surface, the measurement accuracy of the water area to be monitored is correspondingly reduced, and in order to reduce the influence of the bad water surface conditions on the flow rate measurement, the bad water surface conditions need to be identified correspondingly; it is known that when no wave or vortex appears on the water surface in any monitoring area, the area height corresponding to the monitoring area should be kept uniform, i.e. in a horizontal state, and when the corresponding wave or vortex appears on the water surface, the area height corresponding to the monitoring area may change correspondingly, i.e. when the wave appears, the area height corresponding to the position of the wave increases, and when the vortex appears, the area height corresponding to the position of the vortex decreases, based on the characteristic, the current point position data can be correspondingly acquired while the flow rate of the monitoring area is measured to determine whether the corresponding wave or vortex exists, and for the floating object, the current image data can be correspondingly acquired to further determine whether the corresponding monitoring area has larger floating object according to the current image data.
Further, in this embodiment, the foregoing "performing data acquisition for a preset time for each monitoring area to obtain current point cloud data, current image data, and current flow rate data of the corresponding monitoring area" may further include the following steps:
Triggering a point cloud radar, an image radar and a flow rate radar to acquire data for maintaining preset time in each monitoring area based on preset acquisition frequency;
Determining acquisition times based on the preset acquisition frequency and the preset time, and acquiring secondary point cloud data, secondary image data and secondary flow rate data corresponding to each acquisition time respectively;
Grouping each secondary point cloud data, each secondary image data and each secondary flow rate data based on each time node corresponding to each acquisition time respectively to obtain each node data group of the secondary point cloud data, the secondary image data and the secondary flow rate data corresponding to the same time node respectively;
and determining each node data set as current point cloud data, current image data and current flow rate data corresponding to the monitoring area.
For example, in this embodiment, the current point cloud data may be obtained based on the corresponding point cloud radar, the current image data may be obtained based on the corresponding image radar, the current flow rate data may be obtained based on the corresponding flow rate radar, and since the accuracy of the corresponding flow rate measurement is guaranteed when the data is collected in the water area to be monitored for a corresponding preset time, and the water generally has the corresponding fluidity, for the same monitoring area, different water surface conditions may be corresponding at different time instants, for example, when the preset time is 30s, only the corresponding larger ripple appears in 15s, and the corresponding larger vortex and float appear in 16s, and only the smaller float appears in 17s, in this case, it may be seen that the accuracy of the flow rate measurement for 17s is higher, and in order to be able to divide the point cloud data, the image data and the flow rate data at the same time instant, so as to determine the corresponding water surface state at any area may be realized by:
firstly, triggering a point cloud radar, an image radar and a flow rate radar to respectively acquire data for maintaining preset time for each monitoring area based on the same preset acquisition frequency;
Then, the acquisition times corresponding to the preset time can be determined based on the preset acquisition frequency, for example, when the preset time is 30s and the preset acquisition frequency is 2 s/time, the corresponding acquisition times are 15 times, and after the corresponding acquisition times are determined, the point cloud radar, the image radar and the flow rate radar can be controlled to acquire data of the monitoring area in the corresponding acquisition times at the same time so as to obtain secondary point cloud data, secondary image data and secondary flow rate data corresponding to the same acquisition times;
Then, the secondary point cloud data, the secondary image data and the secondary flow rate data corresponding to the same acquisition times can be divided into the same node data set to obtain each node data set corresponding to the monitoring area;
and finally, determining each node data set as current point cloud data, current image data and current flow rate data obtained after data acquisition of the corresponding monitoring area in a preset time.
In step S106, the following are included:
And determining the water area state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determining the abnormal time corresponding to the abnormal state.
For example, in this embodiment, after the corresponding current point cloud data and the corresponding current image data are obtained, the water area state of the monitoring area within the preset time may be corresponding to the two data, and according to the above, when at least one of the larger ripples, the eddies or the floaters occurs in the water surface condition of any monitoring water area, the measurement accuracy of the flow rate measurement corresponding to the time node is lower, that is, the water area state corresponding to the time node may be determined as the abnormal state, and in order to improve the corresponding measurement accuracy, it is necessary to locate the time nodes and determine the abnormal time of the corresponding monitoring area in the abnormal state, so as to perform the supplementary monitoring on the monitoring area with the same duration as the abnormal time in the subsequent steps.
Further, in this embodiment, the above-mentioned "determining the water area state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determining the abnormal time corresponding to the abnormal state" may further include the following steps:
Establishing a two-dimensional plane based on a monitoring area, and establishing point clouds of each area corresponding to each point cloud height on the two-dimensional plane based on secondary point cloud data in any node data set;
Acquiring the area height corresponding to the monitoring area, respectively comparing the cloud height of each point with the area height, and determining a first morphological attribute of the monitoring area under a time node corresponding to the node data set based on a comparison result;
retrieving secondary image data in the same node data set as the secondary point cloud data, and determining a second shape attribute of the monitoring area under a time node corresponding to the node data set based on the secondary image data;
Determining a water area state of the monitoring area under the time node based on the first morphological attribute and the second morphological attribute corresponding to the monitoring area;
And summing all the time nodes corresponding to the abnormal state to obtain the abnormal time corresponding to the abnormal state.
For example, in the present embodiment, the acquisition of the abnormal time may be specifically realized by the following scheme:
Firstly, a corresponding two-dimensional plane can be established based on a horizontal plane of a monitoring area in a static state without external influence, and point clouds of each area corresponding to each point cloud height are established on the two-dimensional plane based on secondary point cloud data in a node data set; it should be noted that, because the point cloud radar may acquire height information of the corresponding point cloud, that is, acquire the height information based on the generated regional point cloud in the corresponding monitoring region, so as to obtain corresponding secondary point cloud data, and further generate each regional point cloud corresponding to the height of each point cloud on the two-dimensional plane according to the secondary point cloud data;
Then, in order to determine the waves and/or eddies located in the monitoring area based on each area point cloud, it is necessary to acquire an area height corresponding to the monitoring area, where the area height may be understood as a height of a horizontal plane of the monitoring area in a stationary state without external influence; because the point cloud heights at the areas corresponding to the corrugated areas are larger than the area heights and the point cloud heights at the areas corresponding to the vortexes are smaller than the area heights, based on the characteristics, the point cloud heights of the point clouds of each area in the same monitoring area can be respectively compared with the area heights, and the first morphological attribute of the monitoring area under the time node of the corresponding node data set can be determined based on the corresponding comparison result;
Then, in order to be able to determine the float located in the monitoring area, secondary image data located in the same node data set as the secondary point cloud data may be retrieved, and a second shape attribute of the monitoring area under the time node of the corresponding node data set may be determined based on the corresponding secondary image data;
then, after determining the corresponding first morphological attribute and second morphological attribute, determining the water area state of the monitoring area under the time node based on the first morphological attribute and the second morphological attribute;
And finally, summing the time nodes corresponding to the abnormal states, thereby obtaining the abnormal time corresponding to the abnormal states.
It should be noted that, in this embodiment, the first shape attribute is determined based on whether there is a corresponding ripple or vortex in the monitoring area, that is, by the point cloud height, and the second shape attribute is determined based on whether there is a corresponding float in the monitoring area, that is, by image recognition.
Further, in this embodiment, the above-mentioned "determining the first morphological attribute of the monitoring area under the time node corresponding to the node data set based on the comparison result" may further include the following steps:
Determining all the regional point clouds with the corresponding point cloud heights smaller than the regional heights as vortex point cloud groups according to the comparison result, and determining all the regional point clouds with the corresponding point cloud heights larger than the regional heights as ripple point cloud groups according to the comparison result;
performing point cloud connection on each area point cloud in adjacent position relation in the vortex point cloud group in the two-dimensional plane to obtain each vortex area in the two-dimensional plane;
performing point cloud connection on each area point cloud in adjacent position relation in the corrugated point cloud group in the two-dimensional plane to obtain each corrugated area in the two-dimensional plane;
Carrying out merging processing based on surrounding relation on each vortex area based on a preset merging strategy;
After the merging processing, determining the first area size and the first area number corresponding to each vortex area, and the second area size and the second area number corresponding to each ripple area as first morphological attributes of the monitoring area under the time node corresponding to the node data set.
For example, in this embodiment, in the monitoring area, all area point clouds with a corresponding point cloud height smaller than the area height may be determined as a vortex point cloud group, and all area point clouds with a corresponding point cloud height greater than the area height may be determined as a vortex point cloud group, by performing point cloud connection on each area point cloud located in the vortex point cloud group and in an adjacent position relationship, each vortex area located in the two-dimensional plane may be obtained, and similarly, each area point cloud located in the corrugated point cloud group and in an adjacent position relationship may be performed, each corrugated area located in the two-dimensional plane may be obtained, where, since the corrugated areas and the vortex areas may have an array state, that is, there may be a corresponding surrounding relationship between each corrugated area, so, in order to be able to determine the number of corrugated areas and vortex areas respectively, it is necessary to determine whether each corrugated area and each area may be combined respectively based on a preset combining strategy, after the corresponding combining process is completed, each corrugated area may be obtained, and the first size of each area and the corresponding number of vortex areas and the second size of each corrugated area may be further determined as the first node size and the second size of the corresponding corrugated area.
Here, the above-mentioned "merging processing based on the surrounding relation for each vortex region based on the preset merging policy" may further include the following steps:
carrying out coordinated processing on the two-dimensional plane to obtain vortex coordinate points respectively positioned in each vortex area, and obtaining vortex coordinate groups respectively corresponding to the vortex areas;
And comparing the vortex coordinate sets in pairs, determining whether the vortex coordinate sets have a surrounding relation or not based on a comparison result, and merging the vortex areas corresponding to the vortex coordinate sets with the surrounding relation into one vortex area.
For example, in the present embodiment, in order to determine whether or not there is a corresponding surrounding relationship between the vortex regions, this may be achieved by:
Firstly, carrying out coordinated processing on two-dimensional coordinates so as to obtain vortex coordinate points forming each vortex region, and respectively forming vortex coordinate groups corresponding to the vortex regions;
Then, based on the characteristics of the vortexes, the vortexes are generally spiral water vortexes, namely, the vortexes are radiated outwards in an array mode through the center points of the vortexes, so that the vortexes can be compared in pairs, whether the vortexes have an enclosing relation or not is determined based on the comparison result, when at least two vortexes have the enclosing relation, at least two vortexes which have the same vortexes can be determined, and at the moment, the vortexes with the enclosing relation can be combined into one vortexes, and therefore the combined screening of the vortexes is completed, and therefore the first area number and the first area size of the vortexes can be accurately obtained in the subsequent method steps.
Here, after the determination of the first morphological attribute of the monitoring area within the preset time is completed based on the corresponding current point cloud data, the second morphological attribute of the monitoring area within the preset time may be determined based on the corresponding current image data, that is, in this embodiment, "determining the second morphological attribute of the monitoring area under the time node corresponding to the node data set based on the secondary image data" may further include the following steps:
Performing binarization processing on the secondary image data to obtain a binarized image, wherein the binarized image comprises all water area pixel points corresponding to first pixel values and all noise pixel points corresponding to second pixel values respectively;
Pixel point connection is carried out on each noise pixel point in adjacent position relation in the secondary image data, so that each noise region in the secondary image data is obtained;
And determining the size and the number of the third areas corresponding to the noise areas as second shape attributes of the monitoring areas under the time nodes corresponding to the node data sets.
For example, in this embodiment, in order to determine the second shape attribute under the time node of the corresponding node data set, the floaters located in the monitored area may be identified based on the image identification method, where the area formed by the floaters may be referred to as a noise area, and the floaters may include leaves, branches, or other objects;
In order to quickly identify the floating objects floating on the monitoring area, binarization processing can be performed on the acquired secondary image data to obtain corresponding binarized images; because the pixel values of all the water area pixel points corresponding to the monitoring area are the same, other pixel points with different pixel values with the water area pixel points are noise pixel points corresponding to floaters, the corresponding water area pixel points can be respectively converted into a first pixel value and the corresponding noise pixel points can be respectively converted into a second pixel value through the binarization processing of the secondary image data, so that the identification of floaters is completed, compared with the existing mode of adopting an image identification model for identification, the method has lower data processing amount, and the corresponding identification precision can be ensured;
After the determination of each noise pixel point is completed, further connecting the noise pixel points in adjacent position relation to obtain each noise region in the secondary image data, and determining the size of each third region and the number of the third regions corresponding to each noise region as a second shape attribute of the monitoring region under the time node of the corresponding node data set.
For example, in this embodiment, after the determination of the first morphological attribute and the second morphological attribute is completed, the water area state of the corresponding time node may be determined based on the first morphological attribute and the second morphological attribute of the corresponding monitoring area, where specific method steps may include:
summing calculation is carried out based on the sizes of all the first areas to obtain a first size total value, and ratio calculation is carried out based on the first size total value and the monitoring size of the corresponding monitoring area to obtain a first size ratio;
Respectively carrying out normalization processing on the first size ratio and the first area number to obtain a first size coefficient and a first quantity coefficient, and carrying out summation calculation on the basis of the first size coefficient and the first quantity coefficient to obtain a vortex coefficient;
Summing calculation is carried out based on the sizes of the second areas to obtain a second size total value, and ratio calculation is carried out based on the second size total value and the monitoring sizes of the corresponding monitoring areas to obtain a second size ratio;
respectively carrying out normalization processing on the second size ratio and the second area number to obtain a second size coefficient and a second number coefficient, and carrying out summation calculation on the basis of the second size coefficient and the second number coefficient to obtain a ripple coefficient;
summing calculation is carried out based on the sizes of all the third areas to obtain a third size total value, and ratio calculation is carried out based on the third size total value and the monitoring size of the corresponding monitoring area to obtain a third size ratio;
Respectively carrying out normalization processing on the third size ratio and the third area number to obtain a third size coefficient and a third number coefficient, and carrying out summation calculation based on the third size coefficient and the third number coefficient to obtain a noise coefficient;
And determining the water area state of the monitoring area under the time node based on the vortex coefficient, the ripple coefficient and the noise coefficient.
For example, in the present embodiment, the corresponding swirl coefficient, ripple coefficient, and noise coefficient can be determined based on the obtained first morphology attribute, respectively, where the determination of the swirl coefficient includes the following processes:
Firstly, acquiring the monitoring size of a corresponding monitoring area, and simultaneously carrying out summation calculation on the sizes of the first areas corresponding to the vortex areas to obtain a corresponding first size total value;
Then, calculating the ratio of the first dimension total value to the monitoring dimension to obtain a corresponding first dimension ratio;
and finally, carrying out normalization processing on the first size proportion and the first area number to obtain a corresponding first size coefficient and a first number coefficient, and obtaining the vortex coefficient of the corresponding vortex area through summation calculation.
Further, the determination of the ripple factor includes the following process:
Firstly, acquiring the monitoring size of a corresponding monitoring area, and simultaneously carrying out summation calculation on the sizes of the second areas corresponding to the corrugated areas to obtain a corresponding total value of the second sizes;
Then, calculating the ratio of the second dimension total value to the monitoring dimension to obtain a corresponding second dimension ratio;
And finally, carrying out normalization processing on the second size proportion and the second area number to obtain a corresponding second size coefficient and a second number coefficient, and obtaining the ripple coefficient of the corresponding ripple area through summation calculation.
Further, the determination of the noise figure includes the following process:
Firstly, acquiring the monitoring size of a corresponding monitoring area, and simultaneously carrying out summation calculation on the sizes of all third areas corresponding to all noise areas to obtain a corresponding third size total value;
Then, calculating the ratio of the third dimension total value to the monitoring dimension to obtain a corresponding third dimension ratio;
And finally, carrying out normalization processing on the third size proportion and the third region number to obtain a corresponding third size coefficient and a third number coefficient, and obtaining the noise coefficient of the corresponding noise region through summation calculation.
For example, in this embodiment, since when a large ripple, a vortex or a floater appears in the monitoring area at a certain time node, the corresponding measurement accuracy is affected, in this case, the water area state of the monitoring area under the corresponding time node needs to be determined based on the obtained vortex coefficient, the obtained ripple coefficient and the obtained noise coefficient, and the specific method process may include the following steps:
invoking vortex weight, ripple weight and noise weight;
performing product calculation on the vortex coefficient and the vortex weight to obtain a vortex evaluation value;
performing product calculation on the ripple coefficient and the ripple weight to obtain a ripple evaluation value;
performing product calculation on the noise coefficient and the noise weight to obtain a noise evaluation value;
When the vortex evaluation value is larger than a first preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
when the ripple evaluation value is larger than a second preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
When the noise evaluation value is larger than a third preset value, determining that the water area state of the monitoring area under the time node is an abnormal state;
When the vortex evaluation value is smaller than or equal to a first preset value, the ripple evaluation value is smaller than or equal to a second preset value and the noise evaluation value is smaller than or equal to a third preset value, summing the vortex evaluation value, the ripple evaluation value and the noise evaluation value to obtain a comprehensive evaluation value;
And when the comprehensive evaluation value is larger than a fourth preset value, determining that the water area state of the monitoring area under the time node is an abnormal state.
For example, in this embodiment, by taking the corresponding vortex weight, ripple weight, and noise weight, and multiplying each weight by each corresponding coefficient to obtain the corresponding vortex evaluation value, ripple evaluation value, and noise evaluation value, it can be known through the above calculation process that the larger the vortex evaluation value, the more the corresponding vortex size and number of vortices, the larger the corresponding ripple size and number of ripples, the larger the noise evaluation value, the more the corresponding noise size and number of noises, based on which it can be determined that when any evaluation value is greater than the corresponding preset value, the water area state under the time node of the monitoring area is determined as an abnormal state, and when all evaluation values are less than or equal to the corresponding preset value, in order to further determine the water area state, it is necessary to determine the comparison result between the integrated evaluation value obtained by summing the vortex evaluation value, the ripple evaluation value, and the noise evaluation value and the corresponding preset value, and when the integrated evaluation value is greater than the corresponding preset value, the water area state can also be indicated as the abnormal state under the time node of the monitoring area.
In this embodiment, the weights and the preset values may be set by the corresponding management end, and the specific numerical values are not limited in this embodiment.
Further, in this embodiment, since the determination of the water area state is determined based on the corresponding numerical comparison result, there may be a certain deviation between the water area state determined by the numerical comparison result and the actual situation, and in order to improve the corresponding determination accuracy, the corresponding weights may be updated and trained in real time, so as to implement the process of updating and iterating the weights, and further, the corresponding determination accuracy may be continuously improved.
Here, the update iteration for each weight can be achieved by the following scheme steps:
Acquiring actual state determination data of a management end under any time node corresponding to any monitoring area, and comparing the water area state corresponding to the monitoring area with the actual state determination data;
If the state actually determines that the data is in a normal state and the water area state is in an abnormal state, performing reduction training on the vortex weight, the ripple weight and the noise weight;
if the state actually determines that the data is in an abnormal state and the water area state is in a normal state, training the vortex weight, the ripple weight and the noise weight in an increasing way;
The trained vortex weight, ripple weight, and noise weight are obtained by the following formula:
Wherein, Weighting eddiesThe number of times of training is increased,Weighting eddiesIs used for the training of the constant value of (a),Weighting eddiesThe number of times of training is reduced and,In order for the trained vortex weights to be applied,Is the ripple weightThe number of times of training is increased,Is the ripple weightIs used for the training of the constant value of (a),Is weighted by ripple coefficientThe number of times of training is reduced and,For the trained ripple weight,Weighting noiseThe number of times of training is increased,Weighting noiseIs used for the training of the constant value of (a),Weighting noiseThe number of times of training is reduced and,Is the weight of the vortex after training.
In step S108, the following are included:
And acquiring data of the supplementing time with the same duration as the abnormal time, and obtaining supplementing point cloud data, supplementing image data and supplementing flow rate data corresponding to the monitoring area.
For example, in this embodiment, by counting the time nodes in the abnormal state, it is possible to determine the abnormal time existing in the preset time, and since the abnormal time corresponds to each time node affecting the measurement accuracy, it is necessary to remove each corresponding time node and perform data acquisition of the complementary time having the same duration as the abnormal time, thereby obtaining the complementary point cloud data, the complementary image data, and the complementary flow rate data corresponding to the monitoring area.
In step S110, the following are included:
and replacing the abnormal flow rate part corresponding to the abnormal time in the current flow rate data in response to the fact that the water area corresponding to the monitoring area is determined to be in a normal state in the supplementing time based on the supplementing point cloud data and the supplementing image data.
For example, since the supplementary point cloud data and the supplementary image data are obtained based on the supplementation, the corresponding supplementary flow rate data needs to be ensured to be in a normal state so as to replace the abnormal flow rate part corresponding to the abnormal time in the current flow rate data, thereby improving the corresponding measurement accuracy, and if the corresponding supplementary point cloud data and the supplementary image data correspond to the abnormal state, the corresponding supplementation needs to be performed again so as to ensure that the current flow rate data after the supplementation has higher measurement accuracy.
In step S112, the following are included:
And determining the water area flow rate corresponding to the water area to be monitored based on the current flow rate data corresponding to each monitoring area.
For example, after the data supplementation for each monitoring area is completed, the water area flow rate corresponding to the water area to be monitored can be determined according to the current flow rate data corresponding to each monitoring area, wherein the water area flow rate obtaining mode can be obtained based on mean value calculation, that is, when the monitoring areas comprise three monitoring areas, the current flow rate data corresponding to the three monitoring areas can be subjected to mean value calculation, and the corresponding mean value result can be determined to be the water area flow rate.
According to the scheme of the embodiment, through carrying out grating treatment on the water area to be monitored, respectively carrying out corresponding data acquisition on the basis of the obtained monitoring areas, carrying out statistical calculation on the current flow velocity data corresponding to each monitoring area, the water area flow velocity corresponding to the water area to be monitored can be obtained, flow velocity measurement of the water area to be monitored is realized in a refined mode, corresponding measurement accuracy is improved, meanwhile, when the flow velocity measurement of each monitoring area is carried out, the water area state of each monitoring area is correspondingly determined on the basis of vortex dimension, ripple dimension and noise dimension, when the water area state of the corresponding monitoring area is in an abnormal state, corresponding supplementary acquisition is needed on the monitoring area, and therefore interference of natural factors or artificial factors on flow velocity measurement can be eliminated, and corresponding measurement accuracy is further improved.
Another embodiment of the present invention provides a system for controlling data quality of a side-scan flow radar, and fig. 3 is a corresponding system block diagram, where the system includes:
The grid processing module is configured to perform grid processing on a water area to be monitored, and each obtained grid point position is respectively determined to be each monitoring area;
the data acquisition module is configured to acquire data of each monitoring area for a preset time to obtain current point cloud data, current image data and current flow rate data of the corresponding monitoring area;
The state determining module is configured to determine a water area state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determine an abnormal time corresponding to the abnormal state;
The supplementary acquisition module is configured to acquire data of supplementary time with the same duration as the abnormal time, and obtain supplementary point cloud data, supplementary image data and supplementary flow rate data corresponding to the monitoring area;
A data replacement module configured to replace the supplementary flow rate data with an abnormal flow rate portion corresponding to the abnormal time in the current flow rate data in response to determining that the water area corresponding to the monitoring area is in a normal state for the supplementary time based on the supplementary point cloud data, the supplementary image data;
And the flow rate determining module is configured to determine the water area flow rate corresponding to the water area to be monitored based on each current flow rate data corresponding to each monitoring area.
In the description provided herein, algorithms and displays are not inherently related to any particular computer, virtual system, or other apparatus. Various general-purpose systems may also be used with examples of the invention. The required structure for a construction of such a system is apparent from the description above. In addition, the present invention is not directed to any particular programming language. It should be appreciated that the teachings of the present invention as described herein may be implemented in a variety of programming languages and that the foregoing description of specific languages is provided for disclosure of preferred embodiments of the present invention.
In the description provided herein, numerous specific details are set forth. However, it is understood that embodiments of the invention may be practiced without these specific details. In some instances, well-known methods, structures and techniques have not been shown in detail in order not to obscure an understanding of this description.
Similarly, it should be appreciated that in the foregoing description of exemplary embodiments of the invention, various features of the invention are sometimes grouped together in a single embodiment, figure, or description thereof for the purpose of streamlining the disclosure and aiding in the understanding of one or more of the various inventive aspects.
Those skilled in the art will appreciate that the modules or units or components of the devices in the examples disclosed herein may be arranged in a device as described in this embodiment, or alternatively may be located in one or more devices different from the devices in this example. The modules in the foregoing examples may be combined into one module or may be further divided into a plurality of sub-modules.
Those skilled in the art will appreciate that the modules in the apparatus of the embodiments may be adaptively changed and disposed in one or more apparatuses different from the embodiments. The modules or units or components of the embodiments may be combined into one module or unit or component and, furthermore, they may be divided into a plurality of sub-modules or sub-units or sub-components.
Furthermore, those skilled in the art will appreciate that while some embodiments described herein include some features but not others included in other embodiments, combinations of features of different embodiments are meant to be within the scope of the invention and form different embodiments.
Furthermore, some of the embodiments are described herein as methods or combinations of method elements that may be implemented by a processor of a computer system or by other means of performing the functions. Thus, a processor with the necessary instructions for implementing the described method or method element forms a means for implementing the method or method element. Furthermore, the elements of the apparatus embodiments described herein are examples of apparatus for performing the functions performed by the elements for the purpose of practicing the invention.
As used herein, unless otherwise specified the use of the ordinal terms "first," "second," "third," etc., to describe a general object merely denote different instances of like objects, and are not intended to imply that the objects so described must have a given order, either temporally, spatially, in ranking, or in any other manner.
While the invention has been described with respect to a limited number of embodiments, those skilled in the art, having benefit of the above description, will appreciate that other embodiments are contemplated within the scope of the invention as described herein. Furthermore, it should be noted that the language used in the specification has been principally selected for readability and instructional purposes, and may not have been selected to delineate or circumscribe the inventive subject matter.

Claims (10)

1.一种侧扫测流雷达数据质量控制方法,其特征在于,包括以下步骤:1. A method for controlling the quality of side-scan current radar data, comprising the following steps: 对待监测水域进行格栅化处理,并将得到的各格栅区域分别确定为各监测区域;The water area to be monitored is gridded, and each grid area obtained is determined as each monitoring area; 对每一监测区域进行预设时间的数据采集,得到对应监测区域的当前点云数据、当前图像数据以及当前流速数据;Collect data for each monitoring area at a preset time to obtain the current point cloud data, current image data and current flow velocity data of the corresponding monitoring area; 基于所述当前点云数据、所述当前图像数据确定对应所述监测区域的水域在所述预设时间内的水域状态,并确定对应处于异常状态下的异常时间;Determine the state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determine the abnormal time corresponding to the abnormal state; 进行与所述异常时间具有相同时长的补充时间的数据采集,得到对应所述监测区域的补充点云数据、补充图像数据以及补充流速数据;Performing data collection for a supplementary time having the same duration as the abnormal time to obtain supplementary point cloud data, supplementary image data and supplementary flow velocity data corresponding to the monitoring area; 响应于基于所述补充点云数据、所述补充图像数据确定对应所述监测区域的水域在所述补充时间内处于正常状态,将所述补充流速数据对位于所述当前流速数据中的对应所述异常时间的异常流速部分进行替换;In response to determining that the water area corresponding to the monitoring area is in a normal state during the supplementary time based on the supplementary point cloud data and the supplementary image data, replacing the abnormal flow velocity portion corresponding to the abnormal time in the current flow velocity data with the supplementary flow velocity data; 基于与各监测区域分别对应的各当前流速数据确定与所述待监测水域对应的水域流速。The water area flow velocity corresponding to the water area to be monitored is determined based on each current flow velocity data corresponding to each monitoring area. 2.根据权利要求1所述的侧扫测流雷达数据质量控制方法,其特征在于,2. The method for controlling the quality of side-scan flow radar data according to claim 1, characterized in that: 对每一监测区域进行预设时间的数据采集,得到对应监测区域的当前点云数据、当前图像数据以及当前流速数据,包括:Data is collected for each monitoring area at a preset time to obtain the current point cloud data, current image data and current flow rate data of the corresponding monitoring area, including: 触发点云雷达、图像雷达以及流速雷达分别基于预设采集频率对每一监测区域进行维持预设时间的数据采集;The point cloud radar, image radar and flow velocity radar are triggered to collect data for each monitoring area for a preset time based on a preset collection frequency; 基于所述预设采集频率以及所述预设时间确定采集次数,并获取分别对应每个采集次数的次级点云数据、次级图像数据以及次级流速数据;Determine the number of acquisitions based on the preset acquisition frequency and the preset time, and obtain secondary point cloud data, secondary image data, and secondary flow velocity data corresponding to each number of acquisitions; 基于与各采集次数分别对应的各时间节点对各次级点云数据、各次级图像数据以及各次级流速数据进行分组,得到分别对应同一时间节点的次级点云数据、次级图像数据以及次级流速数据的各节点数据组;Grouping each secondary point cloud data, each secondary image data and each secondary flow velocity data based on each time node corresponding to each acquisition number, to obtain each node data group of the secondary point cloud data, the secondary image data and the secondary flow velocity data corresponding to the same time node; 将各节点数据组确定为对应所述监测区域的当前点云数据、当前图像数据以及当前流速数据。Each node data group is determined to be the current point cloud data, current image data and current flow velocity data corresponding to the monitoring area. 3.根据权利要求2所述的侧扫测流雷达数据质量控制方法,其特征在于,3. The method for controlling the quality of side-scan flow radar data according to claim 2, characterized in that: 基于所述当前点云数据、所述当前图像数据确定对应所述监测区域的水域在所述预设时间内的水域状态,并确定对应处于异常状态下的异常时间,包括:Determining the water state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determining the abnormal time corresponding to the abnormal state, including: 基于监测区域建立二维平面,并基于位于任一节点数据组中的次级点云数据在所述二维平面上建立分别对应各点云高度的各区域点云;A two-dimensional plane is established based on the monitoring area, and point clouds of each region corresponding to the height of each point cloud are established on the two-dimensional plane based on the secondary point cloud data located in any node data group; 获取与所述监测区域对应的区域高度,并将各点云高度与所述区域高度分别进行高度比较,并基于比较结果确定所述监测区域在对应所述节点数据组的时间节点下的第一形态属性;Acquire a region height corresponding to the monitoring region, compare the height of each point cloud with the region height, and determine a first morphological attribute of the monitoring region at a time node corresponding to the node data group based on the comparison result; 将与所述次级点云数据位于同一节点数据组中的次级图像数据进行调取,并基于所述次级图像数据确定所述监测区域在对应所述节点数据组的时间节点下的第二形态属性;Retrieving secondary image data located in the same node data group as the secondary point cloud data, and determining a second morphological attribute of the monitoring area at a time node corresponding to the node data group based on the secondary image data; 基于对应所述监测区域的所述第一形态属性以及第二形态属性确定所述监测区域在所述时间节点下的水域状态;Determine the water state of the monitoring area at the time node based on the first morphological attribute and the second morphological attribute corresponding to the monitoring area; 将对应异常状态的所有时间节点进行求和计算,得到对应处于异常状态下的异常时间。All time nodes corresponding to the abnormal state are summed up to obtain the abnormal time corresponding to the abnormal state. 4.根据权利要求3所述的侧扫测流雷达数据质量控制方法,其特征在于,4. The method for controlling the quality of side-scan flow radar data according to claim 3, characterized in that: 基于比较结果确定所述监测区域在对应所述节点数据组的时间节点下的第一形态属性,包括:Determining a first morphological attribute of the monitoring area at a time node corresponding to the node data group based on the comparison result includes: 将比较结果为对应点云高度小于所述区域高度的所有区域点云确定为漩涡点云组、将比较结果为对应点云高度大于所述区域高度的所有区域点云确定为波纹点云组;Determine all regional point clouds whose corresponding point cloud heights are less than the regional height as a vortex point cloud group, and determine all regional point clouds whose corresponding point cloud heights are greater than the regional height as a ripple point cloud group; 在所述二维平面中将位于所述漩涡点云组中的处于相邻位置关系的各区域点云进行点云连接,得到位于所述二维平面中的各漩涡区域;Connecting point clouds of each region located in the vortex point cloud group and in adjacent position in the two-dimensional plane to obtain each vortex region located in the two-dimensional plane; 在所述二维平面中将位于所述波纹点云组中的处于相邻位置关系的各区域点云进行点云连接,得到位于所述二维平面中的各波纹区域;Connecting point clouds of the respective regions in the ripple point cloud group that are in an adjacent position relationship in the two-dimensional plane to obtain ripple regions in the two-dimensional plane; 基于预设合并策略对各漩涡区域进行基于包围关系的合并处理;Based on the preset merging strategy, each vortex area is merged based on the encirclement relationship; 在经过所述合并处理后,将各漩涡区域分别对应的各第一区域尺寸以及第一区域数量、与各波纹区域分别对应的各第二区域尺寸以及第二区域数量确定为所述监测区域在对应所述节点数据组的时间节点下的第一形态属性。After the merging process, the first area sizes and the number of first areas corresponding to each vortex area, and the second area sizes and the number of second areas corresponding to each ripple area are determined as the first morphological attributes of the monitoring area at the time node corresponding to the node data group. 5.根据权利要求4所述的侧扫测流雷达数据质量控制方法,其特征在于,5. The method for controlling the quality of side-scan flow radar data according to claim 4, characterized in that: 基于预设合并策略对各漩涡区域进行基于包围关系的合并处理,包括:Based on the preset merging strategy, each vortex area is merged based on the encirclement relationship, including: 对所述二维平面进行坐标化处理,获取分别位于每个漩涡区域的各漩涡坐标点,得到分别对应各漩涡区域的各漩涡坐标组;Performing coordinate processing on the two-dimensional plane, obtaining the vortex coordinate points located in each vortex area, and obtaining the vortex coordinate groups corresponding to the vortex areas; 对各漩涡坐标组进行两两比较,基于比较结果确定各漩涡坐标组之间是否具有包围关系,并将具有包围关系的各漩涡坐标组所分别对应的各漩涡区域合并为一个漩涡区域。The vortex coordinate groups are compared in pairs, and based on the comparison results, it is determined whether the vortex coordinate groups have an enclosing relationship, and the vortex regions corresponding to the vortex coordinate groups having an enclosing relationship are merged into one vortex region. 6.根据权利要求5所述的侧扫测流雷达数据质量控制方法,其特征在于,6. The method for controlling the quality of side-scan flow radar data according to claim 5, characterized in that: 基于所述次级图像数据确定所述监测区域在对应所述节点数据组的时间节点下的第二形态属性,包括:Determining a second morphological attribute of the monitoring area at a time node corresponding to the node data group based on the secondary image data includes: 对所述次级图像数据进行二值化处理,得到二值化图像,其中,所述二值化图像包括分别对应第一像素值的各水域像素点以及分别对应第二像素值的各噪声像素点;Binarizing the secondary image data to obtain a binary image, wherein the binary image includes water area pixel points corresponding to the first pixel value and noise pixel points corresponding to the second pixel value; 在所述次级图像数据中将处于相邻位置关系的各噪声像素点进行像素点连接,得到位于所述次级图像数据中的各噪声区域;Connecting the noise pixels in adjacent positions in the secondary image data to obtain noise regions in the secondary image data; 将各噪声区域对应的各第三区域尺寸以及第三区域数量确定为所述监测区域在对应所述节点数据组的时间节点下的第二形态属性。The sizes of the third regions and the number of the third regions corresponding to the noise regions are determined as the second morphological attributes of the monitoring region at the time node corresponding to the node data group. 7.根据权利要求6所述的侧扫测流雷达数据质量控制方法,其特征在于,7. The method for controlling the quality of side-scan flow radar data according to claim 6, characterized in that: 基于对应所述监测区域的所述第一形态属性以及第二形态属性确定所述监测区域在所述时间节点下的水域状态,包括:Determining the water state of the monitoring area at the time node based on the first morphological attribute and the second morphological attribute corresponding to the monitoring area includes: 基于各第一区域尺寸进行求和计算,得到第一尺寸总值,并基于所述第一尺寸总值与对应所述监测区域的监测尺寸进行比值计算,得到第一尺寸比值;Performing a sum calculation based on the sizes of the first regions to obtain a first size total value, and performing a ratio calculation based on the first size total value and a monitoring size corresponding to the monitoring region to obtain a first size ratio; 分别对所述第一尺寸比值、第一区域数量进行归一化处理,得到第一尺寸系数以及第一数量系数,并基于所述第一尺寸系数与所述第一数量系数进行求和计算,得到漩涡系数;Normalizing the first size ratio and the first area quantity respectively to obtain a first size coefficient and a first quantity coefficient, and performing a sum calculation based on the first size coefficient and the first quantity coefficient to obtain a vortex coefficient; 基于各第二区域尺寸进行求和计算,得到第二尺寸总值,并基于所述第二尺寸总值与对应所述监测区域的监测尺寸进行比值计算,得到第二尺寸比值;Performing a sum calculation based on the sizes of the second regions to obtain a second size total value, and performing a ratio calculation based on the second size total value and a monitoring size corresponding to the monitoring region to obtain a second size ratio; 分别对所述第二尺寸比值、第二区域数量进行归一化处理,得到第二尺寸系数以及第二数量系数,并基于所述第二尺寸系数与所述第二数量系数进行求和计算,得到波纹系数;Normalizing the second size ratio and the second area quantity respectively to obtain a second size coefficient and a second quantity coefficient, and performing a sum calculation based on the second size coefficient and the second quantity coefficient to obtain a ripple coefficient; 基于各第三区域尺寸进行求和计算,得到第三尺寸总值,并基于所述第三尺寸总值与对应所述监测区域的监测尺寸进行比值计算,得到第三尺寸比值;Performing a sum calculation based on the sizes of the third regions to obtain a total value of the third sizes, and performing a ratio calculation based on the total value of the third sizes and the monitoring size of the corresponding monitoring region to obtain a third size ratio; 分别对所述第三尺寸比值、第三区域数量进行归一化处理,得到第三尺寸系数以及第三数量系数,并基于所述第三尺寸系数与所述第三数量系数进行求和计算,得到噪声系数;Normalizing the third size ratio and the third area quantity respectively to obtain a third size coefficient and a third quantity coefficient, and performing a sum calculation based on the third size coefficient and the third quantity coefficient to obtain a noise coefficient; 基于所述漩涡系数、波纹系数以及噪声系数确定所述监测区域在所述时间节点下的水域状态。The water state of the monitoring area at the time node is determined based on the vortex coefficient, ripple coefficient and noise coefficient. 8.根据权利要求7所述的侧扫测流雷达数据质量控制方法,其特征在于,8. The method for controlling the quality of side-scan flow radar data according to claim 7, characterized in that: 基于所述漩涡系数、波纹系数以及噪声系数确定所述监测区域在所述时间节点下的水域状态,包括:Determining the water state of the monitoring area at the time node based on the vortex coefficient, the ripple coefficient, and the noise coefficient includes: 调取漩涡权重、波纹权重以及噪声权重;Call up the vortex weight, ripple weight and noise weight; 将所述漩涡系数与所述漩涡权重进行乘积计算,得到漩涡评价值;Calculate the product of the vortex coefficient and the vortex weight to obtain a vortex evaluation value; 将所述波纹系数与所述波纹权重进行乘积计算,得到波纹评价值;Calculate the product of the ripple coefficient and the ripple weight to obtain a ripple evaluation value; 将所述噪声系数与所述噪声权重进行乘积计算,得到噪声评价值;Calculate the product of the noise coefficient and the noise weight to obtain a noise evaluation value; 当所述漩涡评价值大于第一预设值时,确定所述监测区域在所述时间节点下的水域状态为异常状态;When the vortex evaluation value is greater than a first preset value, determining that the water state of the monitoring area at the time node is an abnormal state; 当所述波纹评价值大于第二预设值时,确定所述监测区域在所述时间节点下的水域状态为异常状态;When the ripple evaluation value is greater than a second preset value, determining that the water state of the monitoring area at the time node is an abnormal state; 当所述噪声评价值大于第三预设值时,确定所述监测区域在所述时间节点下的水域状态为异常状态;When the noise evaluation value is greater than a third preset value, determining that the water state of the monitoring area at the time node is an abnormal state; 当所述漩涡评价值小于等于第一预设值、所述波纹评价值小于等于第二预设值、所述噪声评价值小于等于第三预设值时,将所述漩涡评价值、波纹评价值以及噪声评价值进行求和计算,得到综合评价值;When the vortex evaluation value is less than or equal to a first preset value, the ripple evaluation value is less than or equal to a second preset value, and the noise evaluation value is less than or equal to a third preset value, the vortex evaluation value, the ripple evaluation value, and the noise evaluation value are summed to obtain a comprehensive evaluation value; 当所述综合评价值大于第四预设值时,确定所述监测区域在所述时间节点下的水域状态为异常状态。When the comprehensive evaluation value is greater than a fourth preset value, it is determined that the water state of the monitoring area at the time node is abnormal. 9.根据权利要求8所述的侧扫测流雷达数据质量控制方法,其特征在于,9. The method for controlling the quality of side-scan flow radar data according to claim 8, characterized in that: 获取管理端基于任一监测区域发送的对应任一时间节点下的状态实际确定数据,并将对应该监测区域的水域状态与所述状态实际确定数据进行比较;Obtaining actual state determination data corresponding to any time node sent by the management end based on any monitoring area, and comparing the water state corresponding to the monitoring area with the actual state determination data; 若状态实际确定数据为正常状态,水域状态为异常状态,则对所述漩涡权重、波纹权重以及噪声权重进行减小训练;If the state is actually determined to be a normal state and the water state is an abnormal state, the vortex weight, ripple weight and noise weight are reduced by training; 若状态实际确定数据为异常状态,水域状态为正常状态,则对所述漩涡权重、波纹权重以及噪声权重进行增大训练;If the state actually determines that the data is in an abnormal state and the water state is in a normal state, the vortex weight, ripple weight and noise weight are increased and trained; 通过以下公式得到训练后的漩涡权重、波纹权重以及噪声权重:The trained vortex weight, ripple weight, and noise weight are obtained through the following formula: 其中,为漩涡权重增大训练的次数,为漩涡权重的训练常数值,为漩涡权重减小训练的次数,为经过训练后的漩涡权重,为波纹权重增大训练的次数,为波纹权重的训练常数值,为波纹系数权重减小训练的次数,为经过训练后的波纹权重,为噪声权重增大训练的次数,为噪声权重的训练常数值,为噪声权重减小训练的次数,为经过训练后的漩涡权重。in, Vortex Weight Increase the number of training sessions. Vortex Weight The training constant value of Vortex Weight Reduce the number of training sessions. is the trained vortex weight, is the ripple weight Increase the number of training sessions. is the ripple weight The training constant value of is the ripple coefficient weight Reduce the number of training sessions. is the trained ripple weight, is the noise weight Increase the number of training sessions. is the noise weight The training constant value of is the noise weight Reduce the number of training sessions. is the trained vortex weight. 10.一种侧扫测流雷达数据质量控制系统,其特征在于,包括:10. A side-scan flow radar data quality control system, characterized by comprising: 格栅处理模块,被配置为对待监测水域进行格栅化处理,并将得到的各格栅点位分别确定为各监测区域;A grid processing module is configured to perform grid processing on the water area to be monitored, and determine each grid point obtained as each monitoring area; 数据采集模块,被配置为对每一监测区域进行预设时间的数据采集,得到对应监测区域的当前点云数据、当前图像数据以及当前流速数据;The data acquisition module is configured to collect data for each monitoring area for a preset time to obtain current point cloud data, current image data and current flow velocity data of the corresponding monitoring area; 状态确定模块,被配置为基于所述当前点云数据、所述当前图像数据确定对应所述监测区域的水域在所述预设时间内的水域状态,并确定对应处于异常状态下的异常时间;A state determination module is configured to determine the state of the water area corresponding to the monitoring area within the preset time based on the current point cloud data and the current image data, and determine the abnormal time corresponding to the abnormal state; 补充采集模块,被配置为进行与所述异常时间具有相同时长的补充时间的数据采集,得到对应所述监测区域的补充点云数据、补充图像数据以及补充流速数据;A supplementary acquisition module is configured to collect data for a supplementary time having the same duration as the abnormal time, and obtain supplementary point cloud data, supplementary image data and supplementary flow velocity data corresponding to the monitoring area; 数据替换模块,被配置为响应于基于所述补充点云数据、所述补充图像数据确定对应所述监测区域的水域在所述补充时间内处于正常状态,将所述补充流速数据对位于所述当前流速数据中的对应所述异常时间的异常流速部分进行替换;A data replacement module is configured to replace the abnormal flow velocity portion corresponding to the abnormal time in the current flow velocity data with the supplementary flow velocity data in response to determining that the water area corresponding to the monitoring area is in a normal state during the supplementary time based on the supplementary point cloud data and the supplementary image data; 流速确定模块,被配置为基于与各监测区域分别对应的各当前流速数据确定与所述待监测水域对应的水域流速。The flow rate determination module is configured to determine the water area flow rate corresponding to the water area to be monitored based on the current flow rate data corresponding to each monitoring area.
CN202411319934.1A 2024-09-23 2024-09-23 A method and system for quality control of side-scan flow radar data Active CN118865265B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202411319934.1A CN118865265B (en) 2024-09-23 2024-09-23 A method and system for quality control of side-scan flow radar data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202411319934.1A CN118865265B (en) 2024-09-23 2024-09-23 A method and system for quality control of side-scan flow radar data

Publications (2)

Publication Number Publication Date
CN118865265A CN118865265A (en) 2024-10-29
CN118865265B true CN118865265B (en) 2024-11-29

Family

ID=93157597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202411319934.1A Active CN118865265B (en) 2024-09-23 2024-09-23 A method and system for quality control of side-scan flow radar data

Country Status (1)

Country Link
CN (1) CN118865265B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114997250A (en) * 2022-08-01 2022-09-02 北京海兰信数据科技股份有限公司 Radar radial flow post-processing method and system
CN116125459A (en) * 2023-02-03 2023-05-16 水利部南京水利水文自动化研究所 An Effective Velocity Unit Determination Method for Side Scan Radar

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11506512B2 (en) * 2019-05-22 2022-11-22 TDK Japan Method and system using tightly coupled radar positioning to improve map performance
CN113064129B (en) * 2021-03-03 2023-06-16 湖北中南鹏力海洋探测系统工程有限公司 High-frequency ground wave radar ocean current synthesis method
CN118425555B (en) * 2024-07-03 2025-02-18 武汉大水云科技有限公司 Photoelectric fusion-based water flow velocity measurement method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114997250A (en) * 2022-08-01 2022-09-02 北京海兰信数据科技股份有限公司 Radar radial flow post-processing method and system
CN116125459A (en) * 2023-02-03 2023-05-16 水利部南京水利水文自动化研究所 An Effective Velocity Unit Determination Method for Side Scan Radar

Also Published As

Publication number Publication date
CN118865265A (en) 2024-10-29

Similar Documents

Publication Publication Date Title
CN110287932B (en) Road blocking information extraction method based on deep learning image semantic segmentation
Xiang et al. Superpixel generating algorithm based on pixel intensity and location similarity for SAR image classification
CN103268496B (en) SAR Image Target Recognition Method
CN109272016A (en) Target detection method, device, terminal equipment and computer readable storage medium
CN109344812A (en) An improved clustering-based method for denoising single-photon point cloud data
CN118097207B (en) Ocean engineering surface corrosion image identification method and related device
CN107590512A (en) The adaptive approach and system of parameter in a kind of template matches
CN117635615B (en) Defect detection method and system for realizing punching die based on deep learning
CN107301649A (en) A Coastline Detection Algorithm Based on Superpixel Merging in SAR Image
CN112508970B (en) Point cloud data segmentation method and device
CN117593657A (en) Method and system for processing refined weather forecast data and readable storage medium
CN114742849B (en) Leveling instrument distance measuring method based on image enhancement
CN119741360A (en) An AI calculation method and system for sand ship measurement
CN118865265B (en) A method and system for quality control of side-scan flow radar data
CN106600616A (en) Image background clutter measurement method and system
CN108509835A (en) PolSAR image terrain classification methods based on DFIC super-pixel
CN116740332B (en) Method for positioning center and measuring angle of space target component on satellite based on region detection
US20150308809A1 (en) Apparatus and method for calculating ground surface roughness
CN118364663B (en) Tidal flat terrain remote sensing inversion method and system based on tide level complementary cumulative distribution function
CN113159426A (en) Weather type similarity judgment method and device, electronic equipment and readable storage medium
CN118603050B (en) High-standard farmland monitoring method and system based on multi-source remote sensing technology
CN119178929B (en) SVG-based three-phase unbalanced current detection method and system
CN118965658B (en) Arrangement method and system for transmission line crossing measurement
CN118395190B (en) A training data generation method, device, terminal and storage medium
CN117593409B (en) Topographic map mapping method, device and equipment for reservoir

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant