CN108961177B - OCT image scatter noise suppression method - Google Patents
OCT image scatter noise suppression method Download PDFInfo
- Publication number
- CN108961177B CN108961177B CN201810622714.4A CN201810622714A CN108961177B CN 108961177 B CN108961177 B CN 108961177B CN 201810622714 A CN201810622714 A CN 201810622714A CN 108961177 B CN108961177 B CN 108961177B
- Authority
- CN
- China
- Prior art keywords
- data
- column
- cos
- oct
- sliding window
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000001629 suppression Effects 0.000 title abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 16
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000009499 grossing Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 7
- 230000003287 optical effect Effects 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 12
- 230000009467 reduction Effects 0.000 abstract description 10
- 238000012935 Averaging Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012014 optical coherence tomography Methods 0.000 description 39
- 230000000694 effects Effects 0.000 description 6
- 239000013535 sea water Substances 0.000 description 6
- 239000012528 membrane Substances 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000013329 compounding Methods 0.000 description 2
- 210000003298 dental enamel Anatomy 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 208000010412 Glaucoma Diseases 0.000 description 1
- 208000017442 Retinal disease Diseases 0.000 description 1
- 208000010641 Tooth disease Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 210000004087 cornea Anatomy 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 210000000695 crystalline len Anatomy 0.000 description 1
- 238000005202 decontamination Methods 0.000 description 1
- 230000003588 decontaminative effect Effects 0.000 description 1
- 208000002925 dental caries Diseases 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 208000027993 eye symptom Diseases 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 208000038015 macular disease Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 210000000214 mouth Anatomy 0.000 description 1
- 238000010428 oil painting Methods 0.000 description 1
- 210000001328 optic nerve Anatomy 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 229910052573 porcelain Inorganic materials 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 210000001525 retina Anatomy 0.000 description 1
- 230000002207 retinal effect Effects 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/10—Image enhancement or restoration using non-spatial domain 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- 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/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a method for suppressing scattered point noise of an OCT image. The method comprises the steps of collecting an OCT image of a target object through an OCT system, performing curve smoothing on line data, performing wavelet transformation on the line data after the curve smoothing, constructing window functions with four different window widths, performing windowing processing on the line data after the wavelet transformation, performing weighted averaging on the line data after the four windows are added, performing wavelet inverse transformation again to obtain line data after noise suppression, and performing splicing processing on each line of data to further obtain the OCT image after the noise suppression. The invention furthest reserves the detail characteristics, improves the detail quality of OCT imaging, is beneficial to improving the noise reduction speed, and can be used for online noise reduction application.
Description
Technical Field
The invention belongs to the field of biological image processing, relates to an OCT image processing method, and particularly relates to a method for suppressing speckle noise of an OCT image.
Background
Optical Coherence Tomography (OCT) is a high resolution cross-sectional imaging technique that achieves imaging by collecting photons of the microstructural back-reflected light from the surface region of a sample and interfering with the reference path light. OCT devices were first applied to the in vivo analysis of biological tissues, mainly in the field of ophthalmology; the OCT can be used for high-resolution imaging of structures such as cornea, retina, crystalline lens and the like, photographing of macular diseases, measuring of optic nerve fiber thickness, monitoring and diagnosis of retinal diseases, measuring of retinal structures and the like, so that eye symptoms such as macular degeneration, glaucoma and the like can be diagnosed accurately at an early stage, and the eye is hardly damaged. Then OCT is applied to dentistry to obtain clear OCT image of separated enamel cementum, successfully finish the imaging of the tooth with carious lesion, and can quantitatively and qualitatively analyze the health condition of the oral cavity, estimate the thickness of the enamel, and have certain reports on the pre-diagnosis of tooth diseases. The OCT imaging system not only plays a great role in the medical field, but also has many successful applications in the industrial field. For example, in the field of production of filter membrane/nano-membrane materials, the quality of the membrane can be evaluated in the production process, the local appearance characteristic parameters of the membrane can be measured, and meanwhile, the decontamination and filtration capacity of the membrane can be objectively and accurately evaluated. In the field of printed circuit board manufacturing, particularly for a radio frequency circuit with high requirement on manufacturing precision, the appearance parameters of the radio frequency circuit can be rapidly measured, and the electrical performance and the wave performance can be evaluated. In the field of cultural relic protection and identification, the three-dimensional imaging of OCT is adopted, the organization structure evaluation can be carried out on the paint surfaces of oil paintings and porcelain, an effective observation tool is provided for cleaning and repairing the cultural relics, and meanwhile, partial jades and glassware can also be identified. With the development of the technology, the application field of OCT is being continuously expanded, and the industrial value thereof is being gradually increased.
However, due to the nature of interferometric imaging, speckle noise is inevitably introduced into the image, and this type of noise greatly reduces image contrast. Speckle suppression techniques typically use the acquisition of multiple images of the same object volume, relying on the irrelevancy of these several images to achieve noise suppression. These include spatial, phase, angle and frequency irrelevancy. These techniques, however, require multiple acquisitions at the same location in order to eliminate irrelevant scatter later using averaging or other algorithms.
Disclosure of Invention
In view of the problems in the background art, the present invention is directed to a method for suppressing speckle noise in an OCT image. The method combines frequency and multi-window processing and weighted averaging, can obviously reduce scattered point noise in the OCT image through the uncorrelated characteristic of the frequency, can consider a plurality of windows in the frequency window selection process, improves the frequency resolution ratio, simultaneously furthest retains the detail characteristics, finally improves the detail quality of the OCT imaging, and lays a technical foundation for improving the performance of the OCT technology in a plurality of detection applications by matching with other image processing methods of appearance detection such as imaging and the like.
The present invention proposes to suppress speckle noise using a frequency compounding technique, but this technique is divided into narrow bands by spectral bandwidth to produce uncorrelated speckle patterns, resulting in reduced image axial resolution. The method of the invention independently determines different frequency spectrums and time resolutions by using four orthogonal frequency spectrum windows, and retains the original information of the image to the maximum extent.
The method of the invention directly completes the non-correlated weighting of the image in the frequency domain by utilizing the time-frequency irrelevance under different window functions, has less quantity of images needing to be simultaneously acquired, is favorable for improving the noise reduction speed, and can be used for online imaging.
The technical scheme adopted by the invention comprises the following steps:
firstly, an OCT image of a target object is acquired through an OCT system, and the following modes are adopted for processing aiming at each column of data in the image:
1) performing curve smoothing on each line of data of the image;
2) performing wavelet transformation on a column of data after curve smoothing;
3) selecting four window widths and constructing respective window functions to construct four different window functions, marking as windows 1-4, and respectively windowing the column data after wavelet transformation in the step 2) by using the four different window functions to respectively obtain four windowed column data;
4) and taking weighted average of the four windowed column data, performing wavelet inverse transformation to obtain noise-suppressed column data, and performing splicing processing on each column data to obtain a noise-suppressed OCT image.
The line data is pixel points of one line of the image.
The target object is a biological tissue or a sample, and specifically can be a pearl sample, a human skin tissue, a fruit tissue and the like.
The step 1) is specifically as follows:
1.1) taking a sliding window with a fixed size and taking the maximum value of gray values of all pixel points in the sliding window as Zmax and the minimum value as Zmin, wherein the line data is X;
1.2) judging the gray value of a central pixel point in the sliding window, and if the gray value is in a range of [ Zmin x 1.05, Zmax x 0.95], reserving the gray value of the central pixel point; otherwise, expanding the sliding window and carrying out the next step;
1.3) taking the maximum value of the gray value of the pixel point in the expanded sliding window as Zmax and the minimum value as Zmin, and repeating the step 1.2) to judge whether the gray value of the central pixel point in the expanded sliding window is in the interval of [ Zmin x 1.05, Zmax x 0.95] or not again, and if the gray value is in the interval of [ Zmin x 1.05, Zmax x 0.95], keeping the gray value of the central pixel point;
otherwise, replacing the gray value of the central pixel point in the sliding window by using the median value Zmed of all the gray values of the pixel points in the expanded sliding window, expanding the sliding window again and repeating the steps until the gray value of the pixel point in the current sliding window is in the interval of [ Zmin x 1.05, Zmax x 0.95 ].
In specific implementation, the neighborhood of one pixel point is enlarged each time the sliding window is enlarged.
In the step 2), the data column is X, and the following formula discretization processing is adopted for wavelet transformation:
wherein, WTx(m, n) is the result after wavelet transform, xiIs the ith pixel point in the data column X, k isTotal number of pixels in a data column, tiThe wavelet time variable of the ith pixel point in the data column is represented, m and n respectively represent a time domain position parameter and a frequency domain component, and the time domain position and the frequency domain resolution are adjusted by adjusting the values of m and n; Ψ () represents a discrete wavelet basis function, typically an optional Harr wavelet, Daubechies wavelet, Symlet wavelet, or the like.
In step 3), four different window functions are constructed as follows:
wherein n is a frequency domain resolution parameter, k represents a light source wave number bandwidth of the OCT system, Win1 (n)1)~Win4(n2) Respectively representing four different window functions;
the windowing process with four different window functions is calculated as:
DW1(m,n1)=∫∫WTx(m,n1)·cos(k1·l)·cos(k2·l)·Win1(k1-n1)·Win2(k2-n1)·dk1dk2
DW2(m,n2)=∫∫WTx(m,n2)·cos(k1·l)·cos(k3·l)·Win1(k1-n2)·Win3(k3-n2)·dk1dk3
DW3(m,n3)=∫∫WTx(m,n3)·cos(k1·l)·cos(k3·l)·Win1(k1-n3)·Win4(k4-n3)·dk1dk4
DW4(m,n4)=∫∫WTx(m,n4)·cos(k3·l)·cos(k4·l)·Win1(k3-n4)·Win4(k4-n4)·dk3dk4
wherein, DW1(m,n1)~DW4(m,n4) The results after windowing of four different window functions are respectively shown, i is the optical path difference of a reference arm and a sample arm in the OCT system, and k 1-k 4 respectively show the frequency sliding components in the windowing process.
In the step 4), the weighted average of the four windowed column data is specifically:
DW=(8DW1+4DW2+2DW3+DW4)/15
DW=(8DW1(m,n1)+4DW2(m,n2)+2DW3(m,n3)+DW4(m,n4))/15
wherein DW represents the weighted average of the column data, DW1(m,n1)~DW4(m,n2) Respectively representing the results of windowing processing of four different window functions.
The invention has the beneficial effects that:
the invention uses a frequency complex windowing method to control the scatter noise in the image. The novel noise reduction scheme only needs to acquire images for 4 times, can work on the existing OCT small-sized instrument, and does not need to further increase the computing power of the system.
The sizes of the four spectral windows related in the invention balance the scattered point inhibition capability, the resolution loss and the calculation speed, and an optimal balance point is found. More importantly, the method can be applied in deployment at a video rate, and through the parallel computing capability of the GPU of the modern computer, the method has less insertion computation amount and less resources required by signal processing and image rendering.
In contrast, conventional speckle suppression requires averaging over about 100 images with uncorrelated speckle, with a large normalized standard deviation. The composite processing mode of the invention can obviously improve the signal-to-noise ratio of the images by averaging 4 images with angle intervals. Although the signal-to-noise ratio improvement is low compared to other computationally intensive methods, these hybrid weighting processes do not physically modify the OCT system or require multiple image acquisitions, and the number of scatter points applicable in real-time imaging systems and post-processing is rapidly reduced.
The invention adopts the serialized preprocessing, has less loss of the original information of the image and has certain robustness of the detection effect.
Frequency compounding is achieved by processing each of the images with a plurality of narrow spectral windows, each having independent speckle removal, which are then combined to minimize axial resolution loss, equivalent to a low resolution windowed image multiplied by a full bandwidth image.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is an OCT image of a seawater nucleated pearl input for an embodiment.
FIG. 3 shows the effect of wavelet transform of the intermediate data column of the OCT image of the seawater nucleated pearl sample in FIG. 1; (a) original data; (b) approximation coefficients after wavelet transform; (c) detail coefficient after wavelet transform;
fig. 4 is a final effect diagram of the seawater nucleated pearl sample OCT image of fig. 1 after noise reduction.
Detailed Description
The present invention will be described in further detail below with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The embodiment of the invention and the implementation process thereof are as follows:
1) using a common OCT system to acquire a batch of OCT images of pearl samples serving as target objects; FIG. 2 is an OCT image showing a seawater nucleated pearl.
2) Performing curve smoothing on each line of data of the image; the image resolution of the OCT system is 512 x 512, the length of an image data column is 512, and the bit width is 8 bits, so the range of the gray value is 0-255;
2) performing curve smoothing on each line of data of the image; the image resolution of the OCT system is 512 x 512, the length of an image data column is 512, and the bit width is 8 bits, so the range of the gray value is 0-255;
2.1) taking a sliding window with a fixed size, wherein the maximum value of gray values of all pixel points in the sliding window is Zmax, and the minimum value is Zmin;
2.2) judging the gray value of the central pixel point in the sliding window, and if the gray value is in the interval of [ Zmin x 1.05, Zmax x 0.95], reserving the gray value of the central pixel point; otherwise, expanding a sliding window of one pixel and carrying out the next step;
2.3) taking the maximum value of the gray value of the pixel point in the expanded sliding window as Zmax and the minimum value as Zmin, judging whether the gray value of the central pixel point in the expanded sliding window is in the interval of [ Zmin x 1.05, Zmax x 0.95], and if the gray value is in the interval of [ Zmin x 1.05, Zmax x 0.95], keeping the gray value of the central pixel point; otherwise, replacing the gray value of the central pixel point in the sliding window by the median value Zmed of all the gray values of the pixel points in the expanded sliding window, expanding the sliding window of one pixel again and repeating the steps until the gray value of the pixel point in the current sliding window is in the interval of [ zminx 1.05, Zmax x 0.95 ].
3) Performing wavelet transformation on each line of data after curve smoothing;
the wavelet transformation is carried out by adopting the following formula discretization treatment:
the wavelet basis function of an embodiment takes the sym4 wavelet.
Wherein, WTx(m, n) is the result after wavelet transform, xiIs the ith pixel in the data column X, k is the total number of pixels in the data column, tiRepresenting a column of dataThe wavelet time variable of the ith pixel point, m and n respectively represent a time domain position parameter and a frequency domain component,
FIG. 3 shows the results of wavelet transform of the intermediate data column of the OCT image of the seawater nucleated pearl sample of FIG. 1; the original data is shown in fig. 3(a), the approximation coefficients after wavelet transform are shown in fig. 3(b), and the detail coefficients after wavelet transform are shown in fig. 3 (c).
4) Selecting four window widths, constructing four different window functions, and recording the window functions as windows 1-4;
the wave number bandwidth kappa of the light source of the OCT system built in the implementation is 1/200nm-1Four different window functions are calculated and established as follows:
where n is the frequency domain resolution parameter, Win1 (n)1)~Win4(n2) Respectively representing four different window functions;
5) windowing 1-4 processing is carried out on the column data after the wavelet transformation in the step 3), and a window function is constructed in the step 4) to respectively obtain four windowed column data 1-4;
DW1(m,n1)=∫∫WTx(m,n1)·cos(k1·l)·cos(k2·l)·Win1(k1-n1)·Win2(k2-n1)·dk1dk2
DW2(m,n2)=∫∫WTx(m,n2)·cos(k1·l)·cos(k3·l)·Win1(k1-n2)·Win3(k3-n2)·dk1dk3
DW3(m,n3)=∫∫WTx(m,n3)·cos(k1·l)·cos(k3·l)·Win1(k1-n3)·Win4(k4-n3)·dk1dk4
DW4(m,n4)=∫∫WTx(m,n4)·cos(k3·l)·cos(k4·l)·Win1(k3-n4)·Win4(k4-n4)·dk3dk4
6) taking weighted average of the 1-4 windowed column data obtained in the step 5) to obtain windowed column data DW:
DW=(8DW1(m,n1)+4DW2(m,n2)+2DW3(m,n3)+DW4(m,n4))/15
7) and (3) performing wavelet inverse transformation on each column of data processed in the step 6), and finally combining each column of data to obtain the OCT image after noise suppression.
According to the steps, each column of data is processed to obtain an OCT image after noise suppression.
FIG. 4 shows the final effect of noise reduction of the OCT image of the seawater nucleated pearl sample of FIG. 1;
in the present embodiment, the standard deviation σ of the entire OCT image gradation is calculated; manually selecting a pearl target area of an OCT image of a pearl sample, and calculating the average value beta of the gray scale of the target area; the noise reduction effect of different methods is measured by the uniform noise intensity alpha-sigma/beta; in this embodiment, 40 pearl images are counted, and if the conventional noise reduction method of averaging multiple images is adopted, the average value of α is 0.6, whereas the average value of α is 0.48 by adopting the method proposed by the present invention. Experimental results show that the method has a good noise reduction effect on 40 OCT images. Compared with the existing reported frequency, amplitude and phase averaging scheme, the method obviously improves the noise reduction speed and the number of acquired images, and shows the advantages of the method.
In the embodiment of the present invention, it can be further understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiment may be implemented by instructing the relevant hardware through a program, where the program may be stored in a computer-readable storage medium, where the storage medium includes a ROM/RAM, a magnetic disk, an optical disk, and the like.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.
Claims (5)
1. A method for suppressing OCT image scatter noise is characterized by comprising the following steps: firstly, an OCT image of a target object is acquired through an OCT system, and the following modes are adopted for processing aiming at each column of data in the image:
1) performing curve smoothing on the column data;
2) performing wavelet transformation on a column of data after curve smoothing;
3) selecting four window widths and constructing respective window functions to construct four different window functions, and windowing the column data subjected to wavelet transform in the step 2) by using the four different window functions respectively to obtain four windowed column data respectively;
4) and taking weighted average of the four windowed column data, performing wavelet inverse transformation to obtain noise-suppressed column data, and performing splicing processing on each column data to obtain a noise-suppressed OCT image.
2. The method for suppressing speckle noise of an OCT image according to claim 1, wherein: the step 1) is specifically as follows:
1.1) taking a sliding window with a fixed size, wherein the maximum value of gray values of all pixel points in the sliding window is Zmax, and the minimum value is Zmin;
1.2) judging the gray value of a central pixel point in the sliding window, and if the gray value is in a range of [ Zmin x 1.05, Zmax x 0.95], reserving the gray value of the central pixel point; otherwise, expanding the sliding window and carrying out the next step;
1.3) taking the maximum value of the gray value of the pixel point in the expanded sliding window as Zmax and the minimum value as Zmin, judging whether the gray value of the central pixel point in the expanded sliding window is in the interval of [ Zmin x 1.05 and Zmax x 0.95], and if the gray value is in the interval of [ Zmin x 1.05 and Zmax x 0.95], keeping the gray value of the central pixel point;
otherwise, jumping to step 1.4);
1.4) replacing the gray value of the central pixel point in the sliding window by using the median value Zmed of all the gray values of the pixel points in the expanded sliding window, expanding the sliding window again and repeating the step 1.3) until the gray value of the pixel point in the current sliding window is in the interval of [ zminx 1.05, Zmax x 0.95 ].
3. The method for suppressing speckle noise of an OCT image according to claim 1, wherein: in the step 2), the following formula discretization processing is adopted for wavelet transformation:
wherein, WTx(m, n) is the result after wavelet transform, xiIs the ith pixel in the data column X, k is the total number of pixels in the data column, tiThe wavelet time variable of the ith pixel point in the data column is represented, and m and n respectively represent a time domain position parameter and a frequency domain component; Ψ () represents a discrete wavelet basis function.
4. The method for suppressing speckle noise of an OCT image according to claim 1, wherein: in step 3), four different window functions are constructed as follows:
wherein n is1~n4Is the frequency domain component of four window functions, k represents the wave number bandwidth of the light source of the OCT system, Win1 (n)1)~Win4(n2) Respectively representing four different window functions;
the windowing process with four different window functions is calculated as:
DW1(m,n1)=∫∫WTx(m,n1)·cos(k1·l)·cos(k2·l)·Win1(k1-n1)·Win2(k2-n1)·dk1dk2
DW2(m,n2)=∫∫WTx(m,n2)·cos(k1·l)·cos(k3·l)·Win1(k1-n2)·Win3(k3-n2)·dk1dk3
DW3(m,n3)=∫∫WTx(m,n3)·cos(k1·l)·cos(k3·l)·Win1(k1-n3)·Win4(k4-n3)·dk1dk4
DW4(m,n4)=∫∫WTx(m,n4)·cos(k3·l)·cos(k4·l)·Win1(k3-n4)·Win4(k4-n4)·dk3dk4
wherein, DW1(m,n1)~DW4(m,n4) The results after windowing of four different window functions are respectively shown, i is the optical path difference of a reference arm and a sample arm in the OCT system, and k 1-k 4 respectively show the frequency sliding components in the windowing process.
5. The method for suppressing speckle noise of an OCT image according to claim 1, wherein: in the step 4), the weighted average of the four windowed column data is specifically:
DW=(8DW1(m,n1)+4DW2(m,n2)+2DW3(m,n3)+DW4(m,n4))/15
wherein DW represents the weighted average of the column data, DW1(m,n1)~DW4(m,n2) Respectively representing the results of windowing processing of four different window functions.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810622714.4A CN108961177B (en) | 2018-06-15 | 2018-06-15 | OCT image scatter noise suppression method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810622714.4A CN108961177B (en) | 2018-06-15 | 2018-06-15 | OCT image scatter noise suppression method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108961177A CN108961177A (en) | 2018-12-07 |
CN108961177B true CN108961177B (en) | 2021-08-06 |
Family
ID=64489259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810622714.4A Active CN108961177B (en) | 2018-06-15 | 2018-06-15 | OCT image scatter noise suppression method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108961177B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113409302B (en) * | 2021-07-13 | 2023-07-07 | 浙江科技学院 | OCT image-based corn kernel early mildew identification method |
CN113643396B (en) * | 2021-08-05 | 2023-05-12 | 中国科学院苏州生物医学工程技术研究所 | Swept OCT image processing method and device |
CN116091758B (en) * | 2022-09-09 | 2023-08-25 | 锋睿领创(珠海)科技有限公司 | Image processing method, device, equipment and storage medium based on detail extraction |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073994A (en) * | 2010-12-31 | 2011-05-25 | 哈尔滨工业大学 | Ultrasonic medical image speckle noise inhibition method based on multi-scale anisotropic diffusion |
CN102800064A (en) * | 2012-07-12 | 2012-11-28 | 南京航空航天大学 | OCT (Optical Coherence Tomography) image speckle noise reducing algorithm based on adaptive bilateral filtering |
CN105748041A (en) * | 2016-02-15 | 2016-07-13 | 苏州大学 | Speckle noise suppression system and method in optical coherence tomography imaging |
CN107167805A (en) * | 2017-04-19 | 2017-09-15 | 西安电子科技大学 | Based on the common sparse ISAR high-resolution imaging method of multilayer |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9255914B2 (en) * | 2009-08-13 | 2016-02-09 | Kabushiki Kaisha Toshiba | Ultrasonic diagnosis apparatus and program |
KR102025756B1 (en) * | 2013-01-04 | 2019-09-27 | 삼성전자주식회사 | Method, Apparatus and system for reducing speckles on image |
-
2018
- 2018-06-15 CN CN201810622714.4A patent/CN108961177B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073994A (en) * | 2010-12-31 | 2011-05-25 | 哈尔滨工业大学 | Ultrasonic medical image speckle noise inhibition method based on multi-scale anisotropic diffusion |
CN102800064A (en) * | 2012-07-12 | 2012-11-28 | 南京航空航天大学 | OCT (Optical Coherence Tomography) image speckle noise reducing algorithm based on adaptive bilateral filtering |
CN105748041A (en) * | 2016-02-15 | 2016-07-13 | 苏州大学 | Speckle noise suppression system and method in optical coherence tomography imaging |
CN107167805A (en) * | 2017-04-19 | 2017-09-15 | 西安电子科技大学 | Based on the common sparse ISAR high-resolution imaging method of multilayer |
Non-Patent Citations (5)
Title |
---|
"A Novel Approach for Minimization of Speckle Noise in Ultrasound Imaging Using Redundant Discrete Wavelet Transform";Sanket, at el.;《IOSR Jouranl of Electrical and Electronics Engineering》;20141231;第57-63页 * |
"一种抑制声呐图像散斑噪声的形态学滤波器";郭海涛等;《仪器仪表学报》;20151231;第36卷(第03期);第654-660页 * |
"再现图像细节抑制散斑噪声技术研究";吴育民等;《影像科学与光化学》;20180430;第36卷(第2期);第187-197页 * |
"基于自适应高斯滤波的斑点噪声抑制研究";邓阳阳;《中国优秀硕士学位论文全文数据库 信息科技辑》;20170215(第2期);I138-3267 * |
"小波域高斯混合模型方差估计近红外降噪方法";周扬等;《光电工程》;20111231;第38卷(第8期);第96-100页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108961177A (en) | 2018-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111095349B (en) | Reducing noise in images | |
CN110390650B (en) | OCT image denoising method based on dense connection and generation countermeasure network | |
CN108961177B (en) | OCT image scatter noise suppression method | |
Soomro et al. | Role of image contrast enhancement technique for ophthalmologist as diagnostic tool for diabetic retinopathy | |
CN109961411B (en) | Non-downsampling shear wave transformation medical CT image denoising method | |
WO2009128213A1 (en) | Medical diagnostic device and method of improving image quality of medical diagnostic device | |
CN108272434B (en) | Method and device for processing fundus images | |
CN110101362B (en) | Method for removing image noise related to OCT and OCTA | |
Shinde et al. | Pixelate removal in an image fiber probe endoscope incorporating comb structure removal methods | |
CN108288058B (en) | Improved wavelet threshold knee joint swing signal denoising algorithm | |
Chithra et al. | Hybrid denoising technique for suppressing Gaussian noise in medical images | |
KR100760251B1 (en) | Ultrasound image processing system and method | |
CN109389567B (en) | Sparse filtering method for rapid optical imaging data | |
Saroj | An efficient hybrid approach for medical images enhancement | |
Ramasubramanian et al. | A comprehensive review on various preprocessing methods in detecting diabetic retinopathy | |
EP1309943A2 (en) | Image enhancement | |
JP6740177B2 (en) | Image processing apparatus, image processing method and program | |
CN109146812B (en) | Method for removing hexagonal noise from endoscope image based on frequency domain filtering | |
Sivakumar et al. | Computed radiography skull image enhancement using Wiener filter | |
Ahmed | Image enhancement and noise removal by using new spatial filters | |
Gupta et al. | An efficient technique for speckle noise reduction in ultrasound images | |
Biswas et al. | Modified Gaussian Fuzzy Membership Function for Mixed Noise Reduction from Ultrasound Images | |
Janušauskas et al. | The empirical mode decomposition and the discrete wavelet transform for detection of human cataract in ultrasound signals | |
Al-jaburi et al. | Medical ultrasound image quality enhancement and regions segmentation | |
Anantrasirichai et al. | Adaptive-weighted bilateral filtering for optical coherence tomography |
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 |