Patent Application
20210192268
The present invention belongs to the field of three-dimensional (3D) sonar point cloud image processing, which is in particular related to a distance statistics based method for 3D sonar point cloud image enhancement.
The phase-calculated array 3D sonar system is a new type of submarine 3D acoustic imaging system, which has high application value in such fields as real-time monitoring of port and protection of frogmen (diver) engaged in underwater operations.
However, due to influence from such factors as equipment precision, experiences of operators, underwater environment, changes to surface properties of objects tested as well as data integration and registration, it is inevitable that some noise points may occur to 3D sonar point cloud data as obtained. In addition to noise points produced by such measurement deviations, there might be some discrete points—outliers farther away from principal point clouds, namely point clouds of objects tested in point cloud data during practical application due to influence from such factors as external interferences and barrier. As the first step of pre-processment, filtering during point cloud processing normally has high impact on follow-up processing pipeline. The only way for better follow-up processing, such as registration, feature extraction, curved planar reformation and visualization, is to customize noise points and outliers in reference to follow-up processing during filtering.
In view of aforesaid conditions, it is urgent to propose a reliable and efficient method for 3D sonar point cloud image enhancement with high engineering value.
In view of foregoing factors, the present invention provides a distance statistics based method for 3D sonar point cloud image enhancement; such method features in easy operation, high efficiency and convenience, which can effectively remove outliers to minimize noise, and enhance point cloud image.
On one hand, the present invention proposes a distance statistics based method for 3D sonar point cloud image enhancement, comprising the following steps:
(1) obtaining sonar data, and converting 3D sonar range image information corresponding to sonar data per frame into point cloud data in global coordinate;
(2) using a kd-tree to search the point cloud data, and calculating Euclidean distance Lij between each point cloud data Pi and the nearest K point cloud data; wherein, value range of i and j is 1≤i≤N and 1≤j≤K respectively; N refers to the total quantity of point cloud data;
(3) excluding point cloud data corresponding to mean value of Lij which do not match the certain Gaussian distribution, and complete enhancement of 3D sonar point cloud image.
Specific procedures of the Step (2) are stated as follows:
(2-1) establishing a kd-tree for N points cloud data, and using the such kd-tree to search each point Pi in the cloud data;
(2-2) For each point Pi, using K-NN to search its K nearest point cloud data, and calculating the Euclidean distance Lij between point cloud data Pi and the nearest point cloud data.
Specific procedures of the Step (3) are stated as follows:
(3-1) calculating mean value Li of K Euclidean distance Lij for point Pi;
(3-2) calculating mean value μ and standard deviation σ for N elements in Li;
(3-3) estimating mean value of μ and standard deviation of σ for all Li; selecting point cloud data whose value of corresponding Li element is outside of a--b as outlier; removing the outlier to complete enhancement of 3D sonar point cloud image; wherein, a=μ−α×σ and b=μ+α×σ; α is a real number, referring as expansion coefficient.
On the other hand, the present invention proposes a distance statistics based method for 3D sonar point cloud image enhancement, comprising the following steps:
(1′) obtaining sonar data, and convert 3D sonar range image information corresponding to sonar data per frame into point cloud data in global coordinate;
(2′) using a kd-tree to search the point cloud data, and calculating Euclidean distance Lij between point Pi and all other point cloud data within its neighborhood in distance r; wherein, value range of i and j is 1≤i≤N and 1≤j≤Mi respectively; N refers to the total quantity of point cloud data; Mi refers to the quantity of point cloud data within neighborhood in distance r of point cloud data Pi;
(3′) excluding point cloud data corresponding to mean value of Lij which do not match the certain Gaussian distribution, and complete enhancement of 3D sonar point cloud image.
Specific procedures of the Step (2′) are stated as follows:
(2-1′) establishing a kd-tree for N point cloud data, and use the kd-tree to search each point Pi in the cloud data;
(2-2′) For each point Pi, searching all point cloud data within neighborhood in distance r, and calculating the Euclidean distance Lij between point cloud data Pi and all point cloud data within its neighborhood in distance r.
Specific procedures for the Step (3′) are stated as follows:
(3-1′) calculating mean value Li′ of Mi Euclidean distance Lij for point cloud data Pi;
(3-2′) calculating mean value μ′ and standard deviation σ′ for N elements in Li;
(3-3′) for all Li′, calculating means value of μ′ and standard deviation of σ′ for Gaussian distribution statistics; selecting point cloud data whose value of corresponding Li element is outside of a′--b′ as outlier; remove the outlier to complete enhancement of 3D sonar point cloud image; wherein, a′=μ′−α×σ′ and b′=μ′+α×σ′; α is a real number, referring as expansion coefficient.
As compared with prior arts, the present invention is provided with the following beneficial technical results:
(1) It is applicable to make use of this method for selective removal of most of outliers prior to processing of a mass of point cloud data as collected; this can alleviate influence from system bias and ambient noise, reduce work load for follow-up processing, and improve data efficiency, which is favorable for post processing, such as image reconstruction and image enhancement.
(2) The present invention is available for handy setting of parameters to remove outliers at different degrees according to individual difference between systems and environments. It requires no program for re-modification of system, which is convenient and quick owing to its high practicability and flexibility.
(3) The enhancement method according to present invention has high efficiency, which can apply 3D sonar point cloud image enhancement through special algorithm. It features in short program running time and quick data processing, which can satisfy requirement for real-time performance.
To ensure comprehensive description of the present invention, technical solutions of the present invention are described in details as follows in combination with drawings and preferred embodiments.
S01, obtaining sonar data, and converting 3D sonar range image information corresponding to sonar data per frame into point cloud data in global coordinate;
S02, establishing a kd-tree for N point cloud data, and use the kd-tree to search each point Pi in point cloud data; wherein, value range of i is 1≤i≤N.
S03, for point Pi, using K-NN to search its K nearest point cloud data of point cloud data Pi, and calculating the Euclidean distance Lij between point Pi and the K nearest point cloud data; wherein, value range of j is 1≤j≤K.
S04, calculating mean value Li of K Euclidean distance Lij for point Pi;
S05, calculating mean value μ and standard deviation σ for N elements in Li;
S06, Estimating mean value of μ and standard deviation of σ for all Li; selecting point cloud data whose value of corresponding Li element is outside of a--b as outlier; remove the outlier to complete enhancement of 3D sonar point cloud image; wherein, a=μ−α×σ and b=μ+α×σ; α is a real number, referring as expansion coefficient.
S01′, obtaining sonar data, and converting 3D sonar range image information corresponding to sonar data per frame into point cloud datain global coordinate;
S02′, establishing a kd-tree for N point cloud data, and using the kd-tree to search each each point Pi in point cloud data; value range of i is 1≤i≤N.
S03′, for point Pi, searching all point cloud data within neighborhood in distance r, and calculating the Euclidean distance Lij between point cloud data Pi and all point cloud data within its neighborhood in distance r through calculation; value range of j is 1≤j≤Mi; Mi refers to the quantity of point cloud data within neighborhood in distance r of point cloud data Pi.
S04′, calculating mean value Li′ of Mi Euclidean distance Lij for point Pi;
S05′, calculating mean value μ′ and standard deviation σ′ for N elements in Li;
S06′, for all Li′, calculating mean value of μ′ and standard deviation of σ′ for Gaussian distribution statistics; selecting point cloud data whose value of corresponding Li element is outside of a′--b′ as outlier; removing the outlier to complete enhancement of 3D sonar point cloud image; wherein, a′=μ′−α×σ′ and b′=μ′+α×σ′; α is a real number, referring as expansion coefficient.
The process for establishment of aforesaid kd-tree is stated as follows:
(1) calculating variance of all 3D sonar point cloud data on three dimensions of overall coordinate (x, y and z);
(2) taking dimension as corresponding to the maximum value of three dimension variances as the value of split field;
(3) ordering all point cloud data by its value of split dimension, and select the central data point as the split node;
(4) taking value of the split dimension of split node as hyper-plane of vertical coordinate axis, and dividing the whole space into two parts;
(5) repeating aforesaid steps for data points in left and right subspaces to construct left and right sub-trees until the sub-sample-set divided is void.
Making use of the first method to remove noise of point cloud data, and proceeding with setting of α=1 and K=30. As indicated by data subjecting to numerous practical tests, when aforesaid parameters are used, around 1% points are removed in view of noise.
technical solutions and beneficial results of the present invention are detailed described by the detail method described above. It should be understood that what described above are only optimal embodiments of the present invention, which are not intended to restrict the present invention. Any modification, supplement and equivalent substitution which made according to principles of the present invention will fall into the protection scope of the present invention.
| Number | Date | Country | Kind |
|---|---|---|---|
| 201710135604.0 | Mar 2017 | CN | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/CN2017/115158 | 12/8/2017 | WO | 00 |
