CN113222833A - Side-scan sonar image processing method and device - Google Patents
Side-scan sonar image processing method and device Download PDFInfo
- Publication number
- CN113222833A CN113222833A CN202110401562.7A CN202110401562A CN113222833A CN 113222833 A CN113222833 A CN 113222833A CN 202110401562 A CN202110401562 A CN 202110401562A CN 113222833 A CN113222833 A CN 113222833A
- Authority
- CN
- China
- Prior art keywords
- image
- data
- scan sonar
- denoising
- enhanced
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 13
- 238000001914 filtration Methods 0.000 claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 44
- 230000002146 bilateral effect Effects 0.000 claims abstract description 27
- 238000000034 method Methods 0.000 claims abstract description 22
- 230000002708 enhancing effect Effects 0.000 claims abstract description 7
- 230000004927 fusion Effects 0.000 claims abstract description 6
- 230000006835 compression Effects 0.000 claims description 5
- 238000007906 compression Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000003139 buffering effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20028—Bilateral filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a side-scan sonar image processing method and a side-scan sonar image processing device, wherein the method comprises the following steps: 1) setting a sliding window M for the acquired side scan sonar image to perform sliding caching; 2) denoising the cached image; 3) enhancing the image in a layering way; taking the denoised image as an input signal of bilateral filtering layering, calculating high-frequency image data of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement; 4) and carrying out gray level equalization processing on the enhanced image. The method effectively enhances the details of the image through processing, and ensures the real-time performance of the processing.
Description
Technical Field
The invention relates to an image processing technology, in particular to a side scan sonar image processing method and device.
Background
The side scan sonar is an effective detection device of an underwater detection sonar, and mainly comprises a system device which transmits sound wave signals through an acoustic transducer and collects scattered echo signals to be further processed to form underwater imaging. Therefore, an effective side-scan sonar image real-time processing method is researched, and effective assistance is realized for accelerating side-scan sonar image processing and underwater acoustic image processing engineering application
Disclosure of Invention
The invention aims to solve the technical problem of providing a side-scan sonar image processing method and device aiming at the defects in the prior art.
The technical scheme adopted by the invention for solving the technical problems is as follows: a side scan sonar image processing method comprises the following steps:
1) caching the acquired side-scan sonar image;
performing sliding buffering on an input image according to line frames: the sliding window should be set to an odd value M not smaller than the wavelet coefficient length according to the wavelet coefficient length, and should not be set too large to ensure the processing speed.
2) Denoising the cached image;
3) enhancing the image in a layering way;
taking the denoised image as an input signal of bilateral filtering layering, calculating high-frequency image information of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement, wherein the method specifically comprises the following steps:
3.1) taking the length M of the sliding window in the step 1) as the size of the bilateral filtering blocks, and only processing ceil (M/2) th line data to obtain high-frequency image information of the single-line frame data; wherein the ceil () function is an upgoing integer rounding function;
3.2) weighting and fusing the high-frequency image and the original image data to realize image enhancement;
4) and (5) gray level equalization processing.
According to the scheme, the denoising processing in the step 2) comprises the steps of performing wavelet transformation on the stored image, denoising the image by using soft threshold filtering, and reconstructing the image by using wavelet inverse transformation.
According to the scheme, db 3-order wavelets are adopted in wavelet transformation in the step 2).
According to the scheme, the weighting coefficients adopted by the bilateral filtering in the step 3.1) are as follows:
W(x,y)=s(x,y)*g(x,y)
wherein s (x, y) is a spatial weighting coefficient, g (x, y) is a luminance difference weighting coefficient, W (x, y) is a bilateral filtering weighting coefficient, (x, y) is an image position coordinatec,yc) As coordinates of the image center position, σdIs the spatial Gaussian standard deviation, σrIs the luminance gaussian standard deviation.
According to the scheme, the step 4) is as follows:
4.1) determining a proper sliding cache window according to the gray level distribution characteristics of the image, taking 20ping data with length as cache according to the empirical value, and storing the enhanced image data;
4.2) the collected image gray distribution based on the Gaussian model, the mean value and the variance of the enhanced image data are counted, and an effective value delta t is calculated by using a Gaussian confidence a value;
and 4.3) carrying out linear compression processing on the enhanced data according to the effective value delta t to form image data in an 8-bit standard image format.
According to the method, the side scan sonar image processing device comprises the following steps:
the cache module is used for performing sliding cache on the acquired side-scan sonar images according to line frames:
the denoising module is used for denoising the cached image;
the image layered enhancement module is used for enhancing the image of the denoised image; the method comprises the following steps: taking the denoised image as an input signal of bilateral filtering layering, calculating a high-frequency image and a low-frequency image of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement, wherein the specific steps are as follows:
1) taking the length M of a sliding window in a cache module as the size of bilateral filtering blocks, and only processing ceil (M/2) th line data to obtain high-frequency image and low-frequency image information of the single-line frame data;
2) weighting and fusing the high-frequency image and the original data image to realize image enhancement;
and the balancing module is used for carrying out gray level balancing processing on the enhanced image.
According to the scheme, the denoising processing in the denoising module comprises the steps of performing wavelet transformation on the stored image, denoising the image by using soft threshold filtering, and reconstructing the image by using wavelet inverse transformation.
According to the scheme, the wavelet transformation in the denoising module adopts db3 order wavelet.
According to the scheme, the weighting coefficients adopted by bilateral filtering in the image layering enhancement module are as follows:
W(x,y)=s(x,y)*g(x,y)
wherein s (x, y) is a spatial weighting coefficient, g (x, y) is a luminance difference weighting coefficient, W (x, y) is a bilateral filtering weighting coefficient, (x, y) is an image position coordinatec,yc) As coordinates of the image center position, σdIs the spatial Gaussian standard deviation, σrIs the luminance gaussian standard deviation.
According to the scheme, the equalization module performs gray level equalization processing on the enhanced image, and the method specifically comprises the following steps:
1) determining a proper sliding cache window according to the gray level distribution characteristics of the image, taking 20ping data with length as cache according to the empirical value, and storing the enhanced image data;
2) the method comprises the steps of collecting image gray distribution based on a Gaussian model, counting the mean value and variance of enhanced image data, and calculating an effective value delta t by using a Gaussian confidence a value;
3) and according to the effective value delta t, carrying out linear compression processing on the enhanced data to form image data in an 8-bit standard image format.
The invention has the following beneficial effects:
1. according to the invention, through a sliding cache window with reasonable design, data are rapidly and effectively extracted for two-dimensional image denoising, single-line frame data are enhanced by using two-dimensional bilateral filtering, and a bilateral filtering function is used as a filtering function for image layering, so that the problem of image edge blurring after Gaussian filtering is avoided, a detail layer obtained after filtering is relatively fine, and bilateral filtering can accelerate real-time processing of a side-scanning image, thereby effectively ensuring the real-time property of data processing;
2. the invention designs a reasonable sliding cache window, realizes the smooth processing of image Gaussian model parameters by using the Kalman filtering technology, and realizes the balanced processing of dynamic side scan sonar
3. The invention utilizes the wavelet de-noising image to carry out image enhancement, so that the high-frequency component in the image layering is relatively smooth, the details of the image are effectively enhanced, and excessive noise is not increased.
4. The invention utilizes the image gray scale statistical model to realize that sonar image data with large dynamic range is linearly compressed in a single-line frame mode and converted into standard image data suitable for human eye visual observation and display.
5. The method has high stability and is suitable for side-scan image processing under different underwater acoustic environments.
Drawings
The invention will be further described with reference to the accompanying drawings and examples, in which:
FIG. 1 is a flow chart of a method of an embodiment of the present invention;
FIG. 2 is a diagram illustrating a wavelet function and a scale function according to an embodiment of the present invention;
FIG. 3 is a diagram illustrating a reconstructed wavelet function and a scale function according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of a wavelet three-scale decomposition according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a wavelet single-scale decomposition according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of wavelet inverse transformation reconstruction according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of an image processing process according to an embodiment of the invention;
fig. 8 is a graph of image processing contrast effect according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in fig. 1, a side scan sonar image processing method includes the following steps:
1) caching the acquired side-scan sonar image;
performing sliding buffering on an input image according to line frames: according to the wavelet coefficient length, the sliding window is set to be an odd value M which is not less than the wavelet coefficient length in the denoising treatment, and in addition, the sliding window is not suitable to be set to be too large in order to ensure the treatment speed.
2) Denoising the cached image;
the image is subjected to denoising processing, a large amount of scattered point noise is removed, the detail information of the image is completely reserved, and noise possibly introduced by image enhancement is reduced.
The denoising processing comprises the steps of performing wavelet transformation on a stored image, denoising the image by using soft threshold filtering, and reconstructing the image by using wavelet inverse transformation, wherein the specific steps are as follows:
2.1) wavelet transformation
In the embodiment, the wavelet function adopted by wavelet transform is DaubechAn ies3 order wavelet (referred to as Db3), wherein the Db3 wavelet function and the scale function are shown in fig. 2, and a single-scale filter bank shown in fig. 5 can generate P-scale transformation in a "iterative" manner, where the similar scale J is J-1, J-2, …, and J-P, in the method, P is 3, and the original image is used as the original imageInput of, its row andsum wavelet function ahψThe (-n) convolution is performed and downsampled to obtain two sub-images. Then, the same filtering process is performed on the two sub-image information in the column direction and down-sampled. Four sub-images including three high-frequency component images and one low-frequency component image are obtained through single-scale transformation, the low-frequency component image is subjected to single-scale decomposition again, the steps are repeated, and three-scale decomposition is carried out in total as shown in fig. 4.
2.2) Soft threshold Filtering
Soft threshold filtering processing is carried out on all the obtained high-frequency coefficients, the adopted threshold algorithm is a Visushrink threshold, and the threshold isSigma is the standard deviation of the noise, N is the number of image pixels, and according to an empirical formula, the standard deviation sigma value of the noise is the product of the median value of the input image and 1/0.6745;
2.3) image reconstruction
Performing wavelet inverse transformation on the filtered high-frequency coefficient and the filtered low-frequency coefficient, as shown in fig. 6, wherein a wavelet function and a scale function of the inverse transformation are the order inversions of the wavelet function and the scale function in the wavelet transformation, as shown in fig. 3;
image reconstruction, four images are upsampled and convolved with two one-dimensional filters in each iteration and decomposed into the inverse transform with the wavelet. Adding the results to obtain an approximation with a dimension of j +1, and performing iterative processing until the original image is restoredAnd completing image reconstruction.
3) Enhancing the image in a layering way;
taking the denoised image as an input signal of bilateral filtering layering, calculating high-frequency image information of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement, wherein a processing result is shown in fig. 7 and specifically comprises the following steps:
3.1) taking the length M of the sliding window in the step 1) as the size of the bilateral filtering block, and only processing ceil (M/2) th line data to obtain high-frequency image and low-frequency image information of the single-line frame data;
3.2) weighting and fusing the high-frequency image and the original data image to realize image enhancement;
according to the method, the bilateral filtering function is used as the filtering function of image layering, so that the problem of image edge blurring after Gaussian filtering is avoided, the detail layer obtained after filtering is very fine, and on the other hand, bilateral filtering is easy to realize and can accelerate real-time processing of the side-scan image.
Bilateral filtering benefits from but is different from gaussian filtering, bilateral filtering considers both spatial factors and luminance difference factors, and the weighting formula is shown as follows,
W(x,y)=s(x,y)*g(x,y) (3)
the formula (1) is a spatial weighting coefficient, the formula (2) is a luminance difference weighting coefficient, the formula (3) is a bilateral filtering weighting coefficient, wherein (x, y) is an image position coordinate, (x)c,yc) As the image center position, σdIs the spatial Gaussian standard deviation, σrIs the luminance gaussian standard deviation.
4) And (5) gray level equalization processing.
The method adopts image gray information based on Gaussian distribution to calculate the effective value of the image, utilizes the effective value to carry out linear gray level adjustment with limited amplitude on data in a large dynamic range, converts the data into a standard image suitable for human vision and display, and utilizes Kalman filtering to carry out filtering processing on statistical information so as to realize real-time equalization of the rolling side-scan image.
The method comprises the following specific steps:
4.1) determining a proper sliding cache window according to the gray level distribution characteristics of the image, taking 20ping data with length as cache according to the empirical value, and storing the enhanced image data.
4.2) the collected image gray distribution based on the Gaussian model, the mean value and the variance of the enhanced image data are counted, and the effective value delta t is calculated by using the Gaussian confidence a value.
Due to the fact that the underwater environment is complex, the dynamic range of sonar data acquisition is large, some special bright spots or dark spots can appear in an image, the spots cannot reflect real image information, an effective value delta t of the image needs to be calculated, the delta t serves as the maximum value of the image, and image data are mapped linearly. The calculation of the effective value is performed according to the gray level histogram, because the original image data is in a 16-bit or 32-bit data format, the effective value is not suitable for direct statistical histogram calculation, in this embodiment, the gaussian distribution N (μ, δ) is converted into the standard positive distribution N (0,1) as formula (4), and the value table can be queried to know α according to the effective interval probability value by using the standard positive distribution value table.
According to the experimental effect, pixels with the gray information within 99.4% are within the visual range in a default mode, namely the pixels are within the visual range when the confidence interval is greater than 99.4%. The value α is found to be about 2.5 when the standard positive space book index is used to look up a value α of 99.4%.
The value of X is calculated according to equation (6),
X=αδ+μ (6)
the value of X is also the value of the effective value Δ t.
And 4.3) carrying out linear compression processing on the enhanced data according to the effective value delta t to form image data in an 8-bit standard image format.
The image gain value B is obtained by using delta t, as shown in formula (7),
output data value YiAs shown in formula (8),
Yi=B×xi (8)
in the formula, xiIs the input pixel gray value.
And (3) setting a reasonable sliding cache window, taking the enhanced image as the input of gray statistics, calculating the mean value and the variance of the enhanced image, calculating an effective value X by using an equation (6), calculating a gain value B according to the effective value, and correcting and smoothing the B by using Kalman filtering. Finally, the output data is calculated using equation (8), and the data larger than 255 is set to 255. Finally, the method has the advantages that the finished result is shown in fig. 8, the dynamic range of the original image is large, the image is not beneficial to visual observation of human eyes, after mapping after manual adjustment, the image has high noise and poor contrast, the image processed by the method is remarkably improved, manual adjustment is not needed, and the image data can be automatically processed in real time.
It will be understood that modifications and variations can be made by persons skilled in the art in light of the above teachings and all such modifications and variations are intended to be included within the scope of the invention as defined in the appended claims.
Claims (10)
1. A side scan sonar image processing method is characterized by comprising the following steps:
1) setting a sliding window M for the acquired side scan sonar image to perform sliding caching;
2) denoising the cached image;
3) enhancing the image in a layering way;
taking the denoised image as an input signal of bilateral filtering layering, calculating high-frequency image data of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement, wherein the method specifically comprises the following steps:
3.1) selecting ceil (M/2) th row of data by taking the sliding window length M in the step 1) as the size of the bilateral filtering block to obtain high-frequency image information of the single-line frame data;
3.2) weighting and fusing the high-frequency image and the original image data to realize image enhancement;
4) and carrying out gray level equalization processing on the enhanced image.
2. The side-scan sonar image processing method according to claim 1, wherein the denoising in step 2) includes performing wavelet transform on the stored image, denoising the image by using soft threshold filtering, and reconstructing the image by using inverse wavelet transform.
3. The side-scan sonar image processing method according to claim 2, wherein the sliding window M in step 1) is set to an odd value not less than the wavelet coefficient length according to the wavelet coefficient length.
4. The side-scan sonar image processing method according to claim 1, wherein weighting coefficients used for the bilateral filtering in step 3.1) are as follows:
W(x,y)=s(x,y)*g(x,y)
wherein s (x, y) is a spatial weighting coefficient, g (x, y) is a luminance difference weighting coefficient, W (x, y) is a bilateral filtering weighting coefficient, (x, y) is an image position coordinatec,yc) As coordinates of the image center position, σdIs the spatial Gaussian standard deviation, σrIs the luminance gaussian standard deviation.
5. The side-scan sonar image processing method according to claim 1, wherein the step 4) is specifically as follows:
4.1) determining a corresponding sliding cache window according to the gray level distribution characteristics of the image, and storing the enhanced image data;
4.2) the collected image gray distribution based on the Gaussian model, the mean value and the variance of the enhanced image data are counted, and an effective value delta t is calculated by using a Gaussian confidence a value;
and 4.3) carrying out linear compression processing on the enhanced data according to the effective value delta t to form image data in an 8-bit standard image format.
6. A side scan sonar image processing apparatus, comprising:
the cache module is used for performing sliding cache on the acquired side-scan sonar images according to line frames:
the denoising module is used for denoising the cached image;
the image layered enhancement module is used for enhancing the image of the denoised image; the method comprises the following steps: taking the denoised image as an input signal of bilateral filtering layering, calculating high-frequency image data of single-line frame data based on cache block filtering data, and performing weighted fusion on the high-frequency image and an original data image to realize image enhancement, wherein the method specifically comprises the following steps:
1) taking the length M of a sliding window in a cache module as the size of bilateral filtering blocks, and only selecting ceil (M/2) th line data to obtain high-frequency image information of the single-line frame data;
2) weighting and fusing the high-frequency image and the original data image to realize image enhancement;
and the balancing module is used for carrying out gray level balancing processing on the enhanced image.
7. The side-scan sonar image processing apparatus according to claim 6, wherein the denoising in the denoising module includes performing wavelet transform on the stored image, denoising the image by soft threshold filtering, and reconstructing the image by inverse wavelet transform.
8. The side-scan sonar image processing apparatus according to claim 6, wherein the wavelet transform in the denoising module employs a db3 order wavelet.
9. The side-scan sonar image processing apparatus according to claim 6, wherein weighting coefficients used for bilateral filtering in the image hierarchy enhancing module are as follows:
W(x,y)=s(x,y)*g(x,y)
wherein s (x, y) is a spatial weighting coefficient, g (x, y) is a luminance difference weighting coefficient, W (x, y) is a bilateral filtering weighting coefficient, (x, y) is an image position coordinatec,yc) As coordinates of the image center position, σdIs the spatial Gaussian standard deviation, σrIs the luminance gaussian standard deviation.
10. The side-scan sonar image processing apparatus according to claim 6, wherein the equalization module performs a gray scale equalization process on the enhanced image, specifically as follows:
1) determining a corresponding sliding cache window according to the gray level distribution characteristics of the image, and storing the enhanced image data;
2) the method comprises the steps of collecting image gray distribution based on a Gaussian model, counting the mean value and variance of enhanced image data, and calculating an effective value delta t by using a Gaussian confidence a value;
3) and according to the effective value delta t, carrying out linear compression processing on the enhanced data to form image data in an 8-bit standard image format.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110401562.7A CN113222833B (en) | 2021-04-14 | 2021-04-14 | Side-scan sonar image processing method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110401562.7A CN113222833B (en) | 2021-04-14 | 2021-04-14 | Side-scan sonar image processing method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113222833A true CN113222833A (en) | 2021-08-06 |
CN113222833B CN113222833B (en) | 2022-06-10 |
Family
ID=77087222
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110401562.7A Active CN113222833B (en) | 2021-04-14 | 2021-04-14 | Side-scan sonar image processing method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113222833B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113658081A (en) * | 2021-08-24 | 2021-11-16 | 湖南君瀚信息技术有限公司 | Low-delay image defogging system based on FPGA |
CN116070094A (en) * | 2023-03-14 | 2023-05-05 | 青岛科技大学 | Underwater sound signal processing method based on adaptive wavelet threshold function |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080087777A (en) * | 2008-09-11 | 2008-10-01 | 포항공과대학교 산학협력단 | Conflict evaluated maximum approximated likelihood approach for building maps with sonar sensors of robots |
CN101776452A (en) * | 2009-08-28 | 2010-07-14 | 浙江工业大学 | Active three-dimensional omnidirectional vision-based river width measuring device |
CN102663697A (en) * | 2012-04-01 | 2012-09-12 | 大连海事大学 | Enhancement method of underwater color video image |
CN103325098A (en) * | 2013-07-02 | 2013-09-25 | 南京理工大学 | High dynamic infrared image enhancement method based on multi-scale processing |
CN104240203A (en) * | 2014-09-09 | 2014-12-24 | 浙江工业大学 | Medical ultrasound image denoising method based on wavelet transform and quick bilateral filtering |
CN106127711A (en) * | 2016-06-23 | 2016-11-16 | 浙江工业大学之江学院 | Shearlet conversion and quick two-sided filter image de-noising method |
CN106971379A (en) * | 2017-03-02 | 2017-07-21 | 天津大学 | A kind of underwater picture Enhancement Method merged based on stratified calculation |
CN109738903A (en) * | 2018-12-05 | 2019-05-10 | 华南理工大学 | A kind of side scan sonar real-time two-dimensional imaging method and system |
CN111696061A (en) * | 2020-06-12 | 2020-09-22 | 荆门汇易佳信息科技有限公司 | Noise image reconstruction method with fusion of space domain and transform domain |
-
2021
- 2021-04-14 CN CN202110401562.7A patent/CN113222833B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080087777A (en) * | 2008-09-11 | 2008-10-01 | 포항공과대학교 산학협력단 | Conflict evaluated maximum approximated likelihood approach for building maps with sonar sensors of robots |
CN101776452A (en) * | 2009-08-28 | 2010-07-14 | 浙江工业大学 | Active three-dimensional omnidirectional vision-based river width measuring device |
CN102663697A (en) * | 2012-04-01 | 2012-09-12 | 大连海事大学 | Enhancement method of underwater color video image |
CN103325098A (en) * | 2013-07-02 | 2013-09-25 | 南京理工大学 | High dynamic infrared image enhancement method based on multi-scale processing |
CN104240203A (en) * | 2014-09-09 | 2014-12-24 | 浙江工业大学 | Medical ultrasound image denoising method based on wavelet transform and quick bilateral filtering |
CN106127711A (en) * | 2016-06-23 | 2016-11-16 | 浙江工业大学之江学院 | Shearlet conversion and quick two-sided filter image de-noising method |
CN106971379A (en) * | 2017-03-02 | 2017-07-21 | 天津大学 | A kind of underwater picture Enhancement Method merged based on stratified calculation |
CN109738903A (en) * | 2018-12-05 | 2019-05-10 | 华南理工大学 | A kind of side scan sonar real-time two-dimensional imaging method and system |
CN111696061A (en) * | 2020-06-12 | 2020-09-22 | 荆门汇易佳信息科技有限公司 | Noise image reconstruction method with fusion of space domain and transform domain |
Non-Patent Citations (12)
Title |
---|
CHAO ZUO 等: "Display and detail enhancement for high-dynamic-range infrared images", 《OPTICAL ENGINEERING》 * |
CHAO ZUO 等: "Display and detail enhancement for high-dynamic-range infrared images", 《OPTICAL ENGINEERING》, vol. 50, no. 12, 31 December 2011 (2011-12-31), pages 1 - 9 * |
CHRISTOPHE SINTES 等: "Interferometric side scan sonar signal denoised by wavelets", 《PROCEEDINGS OF SPIE》 * |
CHRISTOPHE SINTES 等: "Interferometric side scan sonar signal denoised by wavelets", 《PROCEEDINGS OF SPIE》, 31 December 2003 (2003-12-31), pages 45 - 55 * |
FRANCESCO BRANCHITTA 等: "New technique for the visualization of high dynamic range infrared images", 《OPTICAL ENGINEERING》 * |
FRANCESCO BRANCHITTA 等: "New technique for the visualization of high dynamic range infrared images", 《OPTICAL ENGINEERING》, vol. 48, no. 9, 9 March 2009 (2009-03-09), pages 1 - 9 * |
李斌: "基于双边滤波器的红外图像细节增强算法研究", 《万方》 * |
李斌: "基于双边滤波器的红外图像细节增强算法研究", 《万方》, 30 July 2015 (2015-07-30), pages 391 * |
焦俊俊: "侧扫声呐图像预处理与拼接技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
焦俊俊: "侧扫声呐图像预处理与拼接技术研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, no. 06, 15 June 2020 (2020-06-15), pages 036 - 187 * |
王晓 等: "一种侧扫声呐图像去噪及均衡化综合处理方法", 《测绘通报》 * |
王晓 等: "一种侧扫声呐图像去噪及均衡化综合处理方法", 《测绘通报》, no. 8, 31 December 2018 (2018-12-31), pages 51 - 55 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113658081A (en) * | 2021-08-24 | 2021-11-16 | 湖南君瀚信息技术有限公司 | Low-delay image defogging system based on FPGA |
CN116070094A (en) * | 2023-03-14 | 2023-05-05 | 青岛科技大学 | Underwater sound signal processing method based on adaptive wavelet threshold function |
Also Published As
Publication number | Publication date |
---|---|
CN113222833B (en) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111583123A (en) | Wavelet transform-based image enhancement algorithm for fusing high-frequency and low-frequency information | |
CN101727658B (en) | Image processing method and device | |
CN101482617B (en) | Synthetic aperture radar image denoising method based on non-down sampling profile wave | |
CN105844601A (en) | Mine image enhancement method based on bilateral filtering and multi-scale Retinex algorithm | |
CN113222833B (en) | Side-scan sonar image processing method and device | |
CN110533614B (en) | Underwater image enhancement method combining frequency domain and airspace | |
CN103295204B (en) | A kind of image self-adapting enhancement method based on non-down sampling contourlet transform | |
CN101847257A (en) | Image denoising method based on non-local means and multi-level directional images | |
CN107403416A (en) | Improvement filtering and the medical ultrasound image denoising method of threshold function table based on NSCT | |
CN109472756A (en) | Image de-noising method based on shearing wave conversion and with directionality local Wiener filtering | |
CN112862716A (en) | Sonar image target enhancement method based on dual-tree complex wavelet and fuzzy theory | |
CN109961408A (en) | The photon counting Image denoising algorithm filtered based on NSCT and Block- matching | |
CN112750090A (en) | Underwater image denoising method and system for improving wavelet threshold | |
CN117252773A (en) | Image enhancement method and system based on self-adaptive color correction and guided filtering | |
CN113850744B (en) | Image enhancement method based on self-adaptive Retinex and wavelet fusion | |
CN111192204A (en) | Image enhancement method, system and computer readable storage medium | |
CN109872289A (en) | Image enchancing method based on improved non-downsampling Contourlet conversion | |
CN117011192A (en) | NSCT domain sonar image denoising method and device based on neutral set and bilateral filtering | |
CN105096274B (en) | Infrared image noise-reduction method based on non-down sampling contourlet domain mixing statistical model | |
CN111652810A (en) | Image denoising method based on wavelet domain singular value differential model | |
CN116452439A (en) | Noise reduction method and device for laser radar point cloud intensity image | |
Bhonsle et al. | Medical Image De-Noising Using Combined Bayes Shrink and Total Variation Techniques | |
Tun et al. | Joint Training of Noisy Image Patch and Impulse Response of Low-Pass Filter in CNN for Image Denoising | |
Tayade et al. | Medical image denoising and enhancement using DTCWT and Wiener filter | |
Arivazhagan et al. | A new hybrid image restoration method based on fusion of spatial and transform domain methods |
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 |