CN102142133A - Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing - Google Patents
Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing Download PDFInfo
- Publication number
- CN102142133A CN102142133A CN2011100982726A CN201110098272A CN102142133A CN 102142133 A CN102142133 A CN 102142133A CN 2011100982726 A CN2011100982726 A CN 2011100982726A CN 201110098272 A CN201110098272 A CN 201110098272A CN 102142133 A CN102142133 A CN 102142133A
- Authority
- CN
- China
- Prior art keywords
- coefficient
- matrix
- image
- sequence
- sampling
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
The invention discloses a mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing, and the method is mainly used for solving the defect of insufficient enhancement effect on the low-contrast medical image in the existing method. The image enhancement method is characterized by comprising the following steps: introducing non-subsampled Directionlet transform and compressive sensing into image enhancement, namely, firstly performing non-subsampled Directionlet transform on the image, and centralizing energy with a high-frequency coefficient by utilizing a compressive sensing technology; then enhancing the concentrated high-frequency coefficient by utilizing a linear enhancement algorithm; and finally reconstructing the enhanced frequency domain representation coefficient through non-subsampled Directionlet transform so as to obtain the enhanced mammary image. By utilizing the mammary X-ray image enhancement method, influence of pathological change region background can be better inhibited, features of a pathological change region in the low-contrast image are obviously enhanced, and the information quantity and readability of the image are improved, thus the method can be used in aided medical radiodiagnosis.
Description
Technical field:
The invention belongs to image processing field, relate to the image enchancing method of non-lower sampling Directionlet conversion and compressed sensing, this method can strengthen the mammary X-ray image of low contrast effectively, and the doctor of auxiliary radiation section carries out medical diagnosis.
Background technology:
Breast cancer is comparatively one of common malignancy disease of women, and women's life and health in serious threat.Development along with the modern medicine imaging technique, various new medical imaging technologies have been widely used in each links such as medical diagnosis, the preceding plan of art, treatment, monitoring after operation, these imaging techniques can obtain patient's various data comprehensively and exactly, for diagnosis, treatment, operation and postoperative evaluation provide information more accurately, wherein the soft X line of molybdenum target is higher with its spatial resolution, comparatively responsive to lump and calcification, and advantage such as equipment needed thereby is simple, cheap, become present early-stage breast cancer diagnosis means the most reliable and commonly used.Yet misdiagnosis rate and rate of missed diagnosis when using the soft X-ray film of breast molybdenum target to diagnose are still higher, and this mainly is because inferior picture quality, the optimum performance of malignant change and observer's visual fatigue or carelessness.Wherein inferior picture quality is embodied in the following aspects: the area-of-interest contrast is relatively poor, suspicious lesions is regional and its surrounding tissue between intensity difference very faint, the lesion region shape is changeable and not of uniform size and obscurity boundary etc.Along with computing machine and development of technologies thereof, utilize computer technology that the mammary X-ray image is strengthened and can effectively address this problem, the help doctor understands image better and judges.
In order to highlight the feature of lesion region, improve the visual effect of mammary X-ray image, the most frequently used method is carried out enhancement process to image exactly.Traditional image enhancement processing technology has obtained to a certain extent and has strengthened effect preferably, but for the mammary X-ray image that has than low contrast, their enhancing effect can not be satisfactory.So far, be directed to the mammary X-ray image, proposed multiple Enhancement Method.
1. anti-sharpening masking method
Anti-sharpening masking method is one of Enhancement Method commonly used in the Flame Image Process, and it is meant on the basis of original image, adds a certain proportion of image radio-frequency component, in the hope of reaching the effect that edge and detailed information strengthen.The advantage of anti-sharpening masking method is to give prominence to edge of image and detailed information preferably, thereby reaches the enhancing effect to image.But this anti-sharpening masking method is owing to the contrast information that does not take into full account image, thereby under the lower situation of contrast, it is unsatisfactory that it strengthens effect.
2. adaptive histogram equalization method
The adaptive histogram equalization method is the effective image enchancing methods of a kind of classics, it adopts sliding window technique, the window area that comprises processed point is carried out histogram equalization, be about to the region histogram distribution and change into even histogram distribution, and give pixel pending in window assignment again according to the mapping relations of histogram and gray scale on this basis.The advantage of adaptive histogram equalization method is the dynamic range that can be good at regulating image, strengthens the details of image simultaneously, but this method has also been amplified noise when improving contrast, so its enhancing effect is still waiting to improve.
3. wavelet transformation adaptive gain disposal route
Wavelet transformation adaptive gain disposal route is a kind of image enchancing method comparatively commonly used, it at first carries out wavelet transformation to image, then according to the adaptive enhancement process of carrying out of the distribution situation of image transformation coefficient, again the conversion coefficient after handling is carried out inverse transformation, the image after being enhanced at last.The advantage of wavelet transformation adaptive gain disposal route is to improve the contrast of image preferably, and noise is had certain robustness, but when the gray difference of the lesion region of image and background area was very little, its enhancing effect still had much room for improvement.
Though three kinds of above-mentioned methods have obtained to a certain extent and have strengthened effect preferably, but because the mammary X-ray image has that contrast is low, noise is than characteristics such as the gray difference of horn of plenty and lesion region and background area are little, so said method is not very good to the enhancing effect of mammary X-ray image.
Summary of the invention
The objective of the invention is to overcome the deficiency of above-mentioned prior art, a kind of mammary X-ray image enchancing method based on non-lower sampling Directionlet conversion and compressed sensing has been proposed, background with effective inhibition image, highlight the lesion region feature, improve the contrast of image, make the enhancing effect of mammary X-ray image more obvious.
Realize that the object of the invention technical thought is: by removing the down-sampling operation in the Directionlet conversion, realized the Directionlet conversion of non-lower sampling, and it is combined with compression sensing method, strengthened clear effect to the mammary X-ray image.Its concrete scheme comprises the steps:
(1) with input picture I
InTransform to the Directionlet transform domain, promptly utilize non-lower sampling Directionlet transfer pair input picture I
InCarry out sub-band division, obtain frequency domain representation coefficient D, this frequency domain representation coefficient D comprises high fdrequency component
And low frequency component
(2) generate white Gaussian noise observing matrix Φ at random.
(3) use the white Gaussian noise observing matrix Φ of generation to the high fdrequency component among the frequency domain representation coefficient D
Observe, obtain observed reading X.
(4) adopt the OMP algorithm, X recovers to observed reading, the high fdrequency component after being restored
(5) to the high fdrequency component after recovering
Carry out linear enhancement process, the high fdrequency component after being enhanced
(6) to the low frequency component among the frequency domain representation coefficient D
And the high fdrequency component after strengthening
Carry out non-lower sampling Directionlet inverse transformation, the image I after finally being enhanced
Out
The present invention has the following advantages:
(1) the present invention has been owing to adopted the Directionlet transform method, thereby can catch the direction detailed information of image effectively, highlighted the minutia of lesion region in the mammary X-ray image significantly, improves the sharpness of image.
(2) the present invention concentrates image energy effectively, thereby can carry out enhancement process to high-frequency information in the concentrated area owing to adopt high fdrequency component after the method for compressed sensing is restored, improves the contrast of image significantly.
(3) the present invention is simple in structure, and calculation cost is low, can carry out enhancement process to image simply and effectively.
Description of drawings:
Fig. 1 is a process flow diagram of the present invention.
Fig. 2 is the schematic diagram that the present invention uses sampling matrix that image is sampled.
Fig. 3 is that the present invention uses the figure as a result after sampling matrix is sampled to image.
Fig. 4 is that the present invention uses sampling matrix frequency coefficient to be carried out the sub-process figure of matrix interpolation and stack.
Fig. 5 uses the present invention and existing method to mammary X-ray treatment of picture effect contrast figure.
Specific embodiments:
With reference to Fig. 1, the present invention includes non-lower sampling Directionlet conversion, compressed sensing, enhancement process and non-lower sampling Directionlet inverse transformation.Concrete steps are as follows:
Step 1: generate sampling matrix.
1.1) be starting point with the true origin, produce the integer vectors d of two linear independences at random
1And d
2, respectively as image changing direction and formation direction, wherein d
1=[a
1, b
1], d
2=[a
2, b
2], a
1Be integer vectors d
1The terminal point horizontal ordinate, b
1Be integer vectors d
1The terminal point ordinate, a
2Be integer vectors d
2The terminal point horizontal ordinate, b
2Be integer vectors d
2The terminal point ordinate;
1.2) with integer vectors d
1And d
2Constitute sampling matrix M
Λ:
In the formula, Z is the integer complete or collected works.
Step 2: according to sampling matrix M
ΛWith input picture I
InBe divided into | det (M
Λ) | individual mutual incoherent subgraph sequence F, each subgraph among the F is F
k, k=0 ..., | det (M
Λ) |-1,
Wherein, | det (M
Λ) | be sampling matrix M
ΛDeterminant, F
kBe arbitrary subgraph, be expressed as: F
k(n)=I
In(M
Λ(n-S
k))
In the formula, n=(n
1n
2) be the position coordinates of pixel in image, F
k(n) the position coordinate is (n in k subgraph of expression
1n
2) pixel, S
kRepresent k subgraph F
kDisplacement vector, I
In(M
Λ(n-S
k)) be k subgraph, it is to utilize displacement vector S
kWith input picture I
InPixel be shifted, re-use sampling matrix M
ΛIt is carried out obtaining after the matrix sampling computing, and sampling principle is chosen integer vectors d respectively as shown in Figure 2 among Fig. 2
1=[1,1] and d
1=[1,1] is as changing direction and the formation direction, and the sampling matrix of generation is:
Utilize displacement vector S
kTo input picture I
InPosition coordinates be shifted after, re-use the sampling matrix M of generation
ΛWith the displacement after the picture position coordinate multiply each other, the subgraph sequence after obtaining sampling as shown in Figure 3, sampling matrix M
ΛInput picture is divided into two subgraphs, and wherein stain is represented subgraph F
1, white point is represented subgraph F
2
Step 3: antithetical phrase graphic sequence F carries out three layers of redundant wavelet decomposition, obtains the high-frequency sub-band coefficient sequence of subgraph sequence F under each yardstick: W
j=(H
j, V
j, D
j), j=1,2,3 and low frequency sub-band coefficient sequence A
3, high-frequency sub-band coefficient sequence W
jIn each sub-band coefficients be:
K=0 ..., | det (M
Λ) |-1, low frequency sub-band coefficient sequence A
3In each sub-band coefficients be:
K=0 ..., | det) M
Λ) |-1;
Wherein,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle horizontal direction detailed information,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle vertical direction detailed information,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle diagonal detailed information,
Be meant expression subgraph F after carrying out the 3rd layer of redundant wavelet decomposition
kThe frequency domain sub-band coefficients of medium and low frequency profile information, frequency domain sub-band coefficients under the j yardstick
Be expressed as successively:
In the formula,
Be the low pass resolution filter coefficient of " 97 " small echo,
Be the high pass resolution filter coefficient of " 97 " small echo,
With
Be respectively right
With
Carry out the filter coefficient that j interpolation obtains,
Represent subgraph to be transformed.
Step 4: according to sampling matrix M
ΛRespectively to high-frequency sub-band coefficient sequence W
j=(H
j, V
j, D
j) and low frequency sub-band coefficient sequence A
3By matrix interpolation and stack, obtain input picture I
InCarry out the frequency domain representation coefficient D=(DH after the non-lower sampling Directionlet conversion
j, DV
j, DD
j, DA
3).
Referring to Fig. 4, being implemented as follows of this step:
4.1) the matrix interpolation computing
Respectively to high-frequency sub-band coefficient sequence W
j=(H
j, V
j, D
j) and low frequency sub-band coefficient sequence A
3Carrying out sampling matrix is M
ΛThe matrix interpolation computing, obtain the high-frequency sub-band coefficient sequence after the interpolation: W
J '=(H
J ', V
J ', D
J ') and low frequency sub-band coefficient sequence A
3 ', be expressed as:
In the formula, M
ΛBe sampling matrix, n=(n
1n
2) position coordinates of remarked pixel point in image,
Expression is carried out the new position coordinates that obtains after the matrix interpolation computing to the position coordinates of pixel, and Λ represents the subscript collection of frequency domain sub-band coefficients;
4.2) stack of sub-band coefficients sequence
Respectively to the high-frequency sub-band coefficient sequence W behind the interpolation arithmetic
J '=(H
J ', V
J ', D
J ') and low frequency sub-band coefficient sequence A
3 'Superpose, obtain input picture I
InCarry out the frequency coefficient D=(DH after the non-lower sampling Directionlet conversion
j, DV
j, DD
j, DA
3),
Wherein: DH
jBe meant the change direction high fdrequency component of detailed information of presentation video, DV
jBe meant the high fdrequency component of presentation video formation direction detailed information, DD
jBe meant that presentation video is changed direction and the high fdrequency component of formation direction detailed information, DA
3Be meant the low frequency component of presentation video low frequency profile information, they are expressed as respectively:
In the formula, n=(n
1n
2) position coordinates of remarked pixel point in image, S
kRepresent k displacement vector, (n+S
k) expression new position coordinates that position coordinates is shifted and obtains, DH
j, DV
j, DD
jForm the high fdrequency component among the frequency domain representation coefficient D
DA
3Form the low frequency component among the frequency domain representation coefficient D
Step 5: generate a gaussian random matrix Φ at random as observing matrix, the size of Φ is M * N, wherein M is an observation frequency, N is the line number that needs the frequency domain representation coefficient of observation, the present invention is on the basis of carrying out repeatedly experiment test and checking, the span that draws observation frequency M is between 250~500, and the observation frequency of this example is taken as 400.
Step 6: use observing matrix Φ that the high fdrequency component among the frequency domain representation coefficient D is observed, promptly with observing matrix Φ that generates at random and the high fdrequency component among the frequency domain representation coefficient D
Multiply each other, obtain observed reading
Wherein:
Be to high fdrequency component
In the change direction coefficient DH of detailed information of presentation video
jThe result who observes is expressed as:
Be to high fdrequency component
The coefficient DV of middle presentation video formation direction detailed information
jThe result who observes is expressed as:
Be to high fdrequency component
In presentation video change direction and the coefficient DD of formation direction detailed information
jThe result who observes is expressed as:
Step 7: the observing matrix Φ that utilizes step 2 to generate, use orthogonal matching pursuit OMP method that observed reading X is recovered, the high frequency coefficient after the recovery is
7.1) iteration count t is initialized as 1, i.e. t=1;
7.2) the calculating observation matrix
In the column vector of each row
With observed reading
In i part coefficient
P row column vector v
pBetween inner product, obtain inner product sequence I, each inner product is I among the I
jBe expressed as:
7.3) calculate the inner product of numerical value maximum among the inner product sequence I:
7.5) the maximum inner product I of usefulness
MaxPosition number λ
tConstitute sequence number collection: Λ
t=(Λ
T-1, λ
t), in the formula, sequence number collection Λ
0Be empty set:
The expression empty set;
7.6) with λ among the observing matrix Φ
tThe column vector of row
Constitute augmented matrix:
In the formula, augmented matrix Φ
0Be empty set:
7.7) with column vector
From observing matrix Φ, remove, be about to λ among the observing matrix Φ
tThe column vector of row
Be changed to empty set:
7.8) utilize augmented matrix Φ
tTo column vector v
pEstimate, obtain new estimated value x
t:
In the formula,<Φ
t, v
pBe meant augmented matrix Φ
tWith column vector v
pInner product,<Φ
t, Φ
tBe meant augmented matrix Φ
tWith the inner product of itself;
7.9) according to new estimated value x
tWith position number λ
tCalculate the column vector after reconstruct recovers:
7.10) utilize new estimated value x
t, augmented matrix Φ
tWith column vector v
pCalculate residual values r
t: r
t=v
p-Φ
tx
t
7.11) increase iteration count t:t=t+1, use residual values r
tReplacement step 7.2) and 7.8) in column vector v
p, repeating step 7.2)~7.11) equal 1/4th of observation frequency M up to iteration count t, promptly
Till, the column vector that obtain this moment
P=1,2 ..., n has constituted through the high frequency coefficient after the reconstruct recovery
Be expressed as
N is a coefficient
Columns, the high frequency coefficient after reconstruct recovers
Constituted the frequency domain high fdrequency component:
This frequency domain high fdrequency component is the main high fdrequency component feature of representative image more concentratedly, and then improves the enhancing effect of subsequent treatment to image.
Step 8: to the high frequency coefficient after the reconstruct recovery
Multiply by amplification factor μ, the high frequency coefficient after being enhanced is:
Be expressed as
The present invention is on the basis of carrying out repeatedly experiment test and checking, and the span that draws reinforcing coefficient μ is between 3~6, and the observation frequency of this example is taken as 4.
Step 9: use sampling matrix M
ΛRespectively to the high fdrequency component after strengthening
With the low frequency component among the frequency domain representation coefficient D
Carry out the matrix sampling, obtain the frequency coefficient sequence
Wherein:
Be to use sampling matrix M
ΛTo high fdrequency component
Carry out the result of matrix sampling,
Be to use sampling matrix M
ΛTo low frequency component DA
3Carry out the result of matrix sampling, be expressed as
In the formula, n=(n
1n
2) be the position coordinates of pixel in image, S
kRepresent k displacement vector, M
Λ(n-S
k) represent by displacement vector S
kPosition coordinates is shifted, and uses sampling matrix the position coordinates after being shifted to be carried out the new position coordinates that obtains after the matrix sampling.
Step 10: with the frequency coefficient sequence
Carry out three layers of redundant wavelet reconstruction respectively, obtain the coefficient sequence Z behind the redundant wavelet reconstruction of j layer
j, coefficient sequence Z
jIn each subsystem number
Be expressed as:
In the formula,
With
Reconstruct low-pass filter coefficients and the reconstruct Hi-pass filter coefficient of representing " 97 " small echo respectively,
With
Be respectively right
With
Carry out the filter coefficient that j interpolation obtains; Obtain the image sequence after the conversion: Z this moment
0=Z
j(j=0).
Step 11: according to sampling matrix M
ΛRespectively to the image sequence Z after the conversion
0In each subgraph
By matrix interpolation and stack, obtain through the enhancing image I after non-lower sampling Directionlet conversion and the compressed sensing enhancing
Out
11.1) matrix interpolation
To image sequence Z
0In each subgraph
Carrying out sampling matrix is M
ΛThe matrix interpolation computing, obtain the image sequence Z after the interpolation, each subgraph Z among the image sequence Z
kBe expressed as:
In the formula, M
ΛBe sampling matrix, n=(n
1n
2) position coordinates of remarked pixel point in image,
The new position coordinates that obtains after the matrix interpolation is carried out in expression to position coordinates;
11.2) the image sequence stack
Image sequence Z after the matrix interpolation computing is superposeed, obtain through the enhancing image I after non-lower sampling Directionlet conversion and the compressed sensing enhancing
Out, be expressed as:
In the formula, n=(n
1n
2) position coordinates of remarked pixel point in image, S
kBe meant k displacement vector, Z
k(n+S
k) represent image sequence Z
kIn k subgraph carry out the new subgraph that obtains after the coordinate displacement.
Advantage of the present invention can further specify by following emulation experiment:
1. simulated conditions
Test pattern used in the present invention comes from the mammary X-ray image in the MIAS database.
2. emulation content
2.1) the present invention uses the sharpening mask method respectively, the adaptive histogram equalization method, these four kinds of methods of wavelet transformation adaptive gain facture and the inventive method have carried out strengthening test to four groups of mammary X-ray images, obtain four groups and strengthen the back image, as shown in Figure 5, wherein Fig. 5 (a) from top to bottom is four groups of original images successively, Fig. 5 (b) is to use existing sharpening mask method to the enhancing of Fig. 5 (a) figure as a result from top to bottom successively, Fig. 5 (c) is to use existing adaptive histogram equalization method to the enhancing of Fig. 5 (a) figure as a result from top to bottom successively, Fig. 5 (d) is to use existing wavelet transformation adaptive gain disposal route to the enhancing of Fig. 5 (a) figure as a result from top to bottom successively, and Fig. 5 (e) is to use the present invention to the enhancing of Fig. 5 (a) figure as a result from top to bottom successively.As seen Fig. 5 (e) is compared with Fig. 5 (b), Fig. 5 (c), Fig. 5 (d) respectively, the present invention can be effectively strengthens the lesion region of galactophore image, suppressed to belong in the image normal structure of background preferably for the influence that strengthens lesion region, make the complex background zone become smoothly, thereby further increased information content of image and readability.
2.2) the present invention is with detection target and background contrast ratio TB based on variance
cValue to using four groups of mammary X-ray images after these four kinds of methods of sharpening mask method, adaptive histogram equalization method, wavelet transformation adaptive gain facture and the inventive method strengthen to test, obtains four groups of contrast ratio TB respectively as judging basis
cValue, as shown in table 1, wherein first row are four groups of original sequence, TB after secondary series is to use the sharpening mask method to four picture group image intensifyings
cValue, TB after the 3rd row are to use the adaptive histogram equalization method to four picture group image intensifyings
cValue, TB after the 4th row are to use wavelet transformation adaptive gain facture to four picture group image intensifyings
cValue, TB after the 5th row are to use the inventive method to four picture group image intensifyings
cValue, wherein, contrast ratio TB
cValue representation is:
TB
c=δ
μ/σ
In the formula, δ
μMeasured the difference between the average gray ratio that detects target and background in original image and the enhancing back image, δ
μBe expressed as:
Wherein
With
Be meant the average that detects target T and background B,
With
Be meant the average that strengthens the back image.σ has measured to strengthen and has detected target in the image with respect to the reduction degree that detects target gray level divergence in the original image,
In the formula
With
Be respectively original image and strengthen the variance that detects target in the image of back, therefore detect target and background contrast ratio TB
cValue is big more to show that then the enhancing effect of image is good more.
Table 1 strengthens evaluation of result TB
c
By table 1 as seen, the TB of the present invention after four groups of images are tested
cValue is obviously greater than other three kinds of existing methods, and therefore enhancing effect of the present invention is better than these three kinds of existing methods of sharpening mask method, adaptive histogram equalization method and wavelet transformation adaptive gain disposal route on objective metric.
To sum up, the present invention has improved the contrast that detects target and background, has strengthened the detailed information of galactophore image effectively.
Claims (6)
1. the mammary X-ray image enchancing method based on non-lower sampling Directionlet conversion and compressed sensing comprises the steps:
(1) with input picture I
InTransform to the Directionlet transform domain, promptly utilize non-lower sampling Directionlet transfer pair input picture I
InCarry out sub-band division, obtain frequency domain representation coefficient D, this frequency domain representation coefficient D comprises high fdrequency component
And low frequency component
(2) generate white Gaussian noise observing matrix Φ at random.
(3) use the white Gaussian noise observing matrix Φ of generation to the high fdrequency component among the frequency domain representation coefficient D
Observe, obtain observed reading X.
(4) adopt the OMP algorithm, X recovers to observed reading, the high fdrequency component after being restored
(5) to the high fdrequency component after recovering
Carry out linear enhancement process, the high fdrequency component after being enhanced
2. according to claims 1 described mammary X-ray image enchancing method, the described non-lower sampling Directionlet transfer pair input picture I that utilizes of step (1) wherein
InCarry out sub-band division, carry out as follows:
(1a) produce the integer vectors d of two linear independences at random
1And d
2, respectively as image changing direction and formation direction, wherein d
1=[a
1, b
1], d
2=[a
2, b
2], a
1Be integer vectors d
1The terminal point horizontal ordinate, b
1Be integer vectors d
1The terminal point ordinate, a
2Be integer vectors d
2The terminal point horizontal ordinate, b
2Be integer vectors d
2The terminal point ordinate;
(1b) with integer vectors d
1And d
2Constitute sampling matrix M
Λ:
In the formula, Z is the integer complete or collected works;
(1c) according to sampling matrix M
ΛWith input picture I
InBe divided into | det (M
Λ) | individual mutual incoherent subgraph sequence F, each subgraph among the F is F
k, k=0 ..., | det (M
Λ) |-1,
Wherein, | det (M
Λ) | be sampling matrix M
ΛDeterminant, F
kBe arbitrary subgraph, be expressed as:
F
k(n)=I
in(M
Λ(n-S
k))
In the formula, n=(n
1n
2) be the position coordinates of pixel in image, F
k(n) the position coordinate is (n in k subgraph of expression
1n
2) pixel, S
kRepresent k subgraph F
kDisplacement vector, I
In(M
Λ(n-S
k)) be k subgraph, it is to utilize displacement vector S
kWith input picture I
InPixel be shifted, re-use sampling matrix M
ΛIt is carried out obtaining after the matrix sampling computing;
(1d) antithetical phrase graphic sequence F carries out three layers of redundant wavelet decomposition, obtains the high-frequency sub-band coefficient sequence of subgraph sequence F under each yardstick: W
j=(H
j, V
j, D
j), j=1,2,3 and low frequency sub-band coefficient sequence A
3, high-frequency sub-band coefficient sequence W
jIn each sub-band coefficients be:
K=0 ..., | det (M
Λ) |-1, low frequency sub-band coefficient sequence A
3In each sub-band coefficients be:
K=0 ..., | det (M
Λ) |-1;
Wherein,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle horizontal direction detailed information,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle vertical direction detailed information,
Be meant expression subgraph F after carrying out the redundant wavelet decomposition of j layer
kThe frequency domain sub-band coefficients of middle diagonal detailed information,
Be meant expression subgraph F after carrying out the 3rd layer of redundant wavelet decomposition
kThe frequency domain sub-band coefficients of medium and low frequency profile information, frequency domain sub-band coefficients under the j yardstick
Be expressed as successively:
In the formula,
Be the low pass resolution filter coefficient of " 97 " small echo,
Be the high pass resolution filter coefficient of " 97 " small echo,
With
Be respectively right
With
Carry out the filter coefficient that j interpolation obtains,
Represent subgraph to be transformed;
(1e) to high-frequency sub-band coefficient sequence W
j=(H
j, V
j, D
j) and low frequency sub-band coefficient sequence A
3Carrying out sampling matrix is M
ΛThe matrix interpolation computing, obtain the high-frequency sub-band coefficient sequence after the interpolation: W
J '=(H
J ', V
J ', D
J ') and low frequency sub-band coefficient sequence A
3 ', be expressed as:
In the formula, M
ΛBe sampling matrix, n=(n
1n
2) position coordinates of remarked pixel point in image,
Expression is carried out the new position coordinates that obtains after the matrix interpolation computing to the position coordinates of pixel, and Λ represents frequency domain sub-band coefficients W
kThe subscript collection;
(1f) respectively to the high-frequency sub-band coefficient sequence W behind the interpolation arithmetic
J '=(H
J ', V
J ', D
J ') and low frequency sub-band coefficient sequence A
3 'Superpose, obtain input picture I
InCarry out the frequency domain representation coefficient D=(DH after the non-lower sampling Directionlet conversion
j, DV
j, DD
j, DA
3),
Wherein: DH
jBe meant the change direction high fdrequency component of detailed information of presentation video, DV
jBe meant the high fdrequency component of presentation video formation direction detailed information, DD
jBe meant that presentation video is changed direction and the high fdrequency component of formation direction detailed information, DA
3Be meant the low frequency component of presentation video low frequency profile information, they are expressed as respectively:
In the formula, n=(n
1n
2) position coordinates of remarked pixel point in image, S
kRepresent k displacement vector, (n+S
k) expression new position coordinates that position coordinates is shifted and obtains, DH
j, DV
j, DD
jForm the high fdrequency component among the frequency domain representation coefficient D
DA
3Form the low frequency component among the frequency domain representation coefficient D
3. according to claims 1 described mammary X-ray figure Enhancement Method, wherein the white Gaussian noise observing matrix Φ of the described use generation of step (3) is to the high fdrequency component among the frequency domain representation coefficient D
Observe, be meant observing matrix Φ that use generates at random and the high fdrequency component among the frequency domain representation coefficient D
Multiply each other, obtain observed reading
Wherein:
Be to high fdrequency component
In the change direction coefficient DH of detailed information of presentation video
jThe result who observes is expressed as:
Be to high fdrequency component
The coefficient DV of middle presentation video formation direction detailed information
jThe result who observes is expressed as:
4. according to claims 1 described mammary X-ray figure Enhancement Method, the described employing of step (4) OMP algorithm wherein, X recovers to observed reading, carries out as follows:
(4a) iteration count t is initialized as 1, i.e. t=1;
(4b) calculating observation matrix
In the column vector of each row
With observed reading
In i part coefficient
P row column vector v
pBetween inner product, obtain inner product sequence I, each inner product is I among the I
jBe expressed as:
In the formula, n is a coefficient
Columns, N is the columns of observing matrix Φ;
(4e) with maximum inner product I
MaxPosition number λ
tConstitute sequence number collection: Λ
t=(Λ
T-1, λ
t), in the formula, sequence number collection Λ
0Be empty set:
The expression empty set;
(4f) with λ among the observing matrix Φ
tThe column vector of row
Constitute augmented matrix:
In the formula, augmented matrix Φ
0Be empty set:
(4g) with column vector
From observing matrix Φ, remove, be about to λ among the observing matrix Φ
tThe column vector of row
Be changed to empty set:
(4h) utilize augmented matrix Φ
tTo column vector v
pEstimate, obtain new estimated value x
t:
In the formula,<Φ
t, v
pBe meant augmented matrix Φ
tWith column vector v
pInner product,<Φ
t, Φ
tBe meant augmented matrix Φ
tWith the inner product of itself;
(4i) according to new estimated value x
tWith position number λ
tCalculate the column vector after reconstruct recovers:
(4j) utilize new estimated value x
t, augmented matrix Φ
tWith column vector v
pCalculate residual values r
t: r
t=v
p-Φ
tx
t
(4k) increase iteration count t:t=t+1, use residual values r
tReplacement step (4b) and (4h) in column vector v
p, repeating step (4b)~(4k) equals 1/4th of observation frequency M up to iteration count t, promptly
Till, the column vector that obtain this moment
P=1,2 ..., n has constituted through the high frequency coefficient after the reconstruct recovery
Be expressed as
N is a coefficient
Columns, the high frequency coefficient after reconstruct recovers
Constituted the frequency domain high fdrequency component:
This frequency domain high fdrequency component is the main high fdrequency component feature of representative image more concentratedly, and then improves the enhancing effect of subsequent treatment to image.
5. according to claims 1 described mammary X-ray figure Enhancement Method, wherein step (5) is described to the high fdrequency component after recovering
Carry out linear enhancement process, be meant the frequency domain high frequency coefficient after the reconstruct recovery
Multiply by amplification factor μ, the frequency domain high frequency coefficient after being enhanced is:
Be expressed as
The present invention is on the basis of carrying out repeatedly experiment test and checking, and the span that draws reinforcing coefficient μ is between 3~6.
6. according to claims 1 described mammary X-ray image enchancing method, wherein step (6) is described to the low frequency component among the frequency domain representation coefficient D
And the high fdrequency component after strengthening
Carry out non-lower sampling Directionlet inverse transformation, carry out as follows:
(6a) use sampling matrix M
ΛRespectively to the high fdrequency component after strengthening
With the low frequency component among the frequency domain representation coefficient D
Carry out the matrix sampling, obtain the frequency coefficient sequence
Wherein:
Be to use sampling matrix M
ΛTo high fdrequency component
Carry out the result of matrix sampling,
Be to use sampling matrix M
ΛTo low frequency component DA
3Carry out the result of matrix sampling, be expressed as
In the formula, n=(n
1n
2) be the position coordinates of pixel in image, S
kRepresent k displacement vector, M
Λ(n-S
k) represent by displacement vector S
kPosition coordinates is shifted, and uses sampling matrix the position coordinates after being shifted to be carried out the new position coordinates that obtains after the matrix sampling;
(6b) with the frequency coefficient sequence
Carry out three layers of redundant wavelet reconstruction respectively, obtain the coefficient sequence Z behind the redundant wavelet reconstruction of j layer
j, coefficient sequence Z
jIn each subsystem number
Be expressed as:
In the formula,
With
Reconstruct low-pass filter coefficients and the reconstruct Hi-pass filter coefficient of representing " 97 " small echo respectively,
With
Be respectively right
With
Carry out the filter coefficient that j interpolation obtains; Obtain the image sequence after the conversion: Z this moment
0=Z
j(j=0).
(6c) to the image sequence Z after the conversion
0In each subgraph
Carrying out sampling matrix is M
ΛThe matrix interpolation computing, obtain the image sequence Z after the interpolation, each subgraph Z among the image sequence Z
kBe expressed as:
In the formula, M
ΛBe sampling matrix, n=(n
1n
2) position coordinates of remarked pixel point in image,
The new position coordinates that obtains after the matrix interpolation is carried out in expression to position coordinates;
(6d) the image sequence Z after the matrix interpolation computing is superposeed, obtain through the enhancing image I after non-lower sampling Directionlet conversion and the compressed sensing enhancing
Out, be expressed as:
In the formula, n=(n
1n
2) position coordinates of remarked pixel point in image, S
kBe meant k displacement vector, Zk (n+S
k) represent image sequence Z
kIn k subgraph carry out the new subgraph that obtains after the coordinate displacement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110098272 CN102142133B (en) | 2011-04-19 | 2011-04-19 | Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110098272 CN102142133B (en) | 2011-04-19 | 2011-04-19 | Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102142133A true CN102142133A (en) | 2011-08-03 |
CN102142133B CN102142133B (en) | 2013-01-23 |
Family
ID=44409622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110098272 Expired - Fee Related CN102142133B (en) | 2011-04-19 | 2011-04-19 | Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102142133B (en) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663709A (en) * | 2012-04-28 | 2012-09-12 | 大连海事大学 | Vessel enhancement method for X-ray angiogram |
CN102903099A (en) * | 2012-09-05 | 2013-01-30 | 西安电子科技大学 | Color image edge detection method based on directionlet conversion |
CN103310429A (en) * | 2013-03-06 | 2013-09-18 | 西安电子科技大学 | Image enhancement method based on hidden Markov tree (HMT) model in directionlet domain |
CN103473797A (en) * | 2013-09-16 | 2013-12-25 | 电子科技大学 | Image space-domain reducible reconstructing method based on compressed sensing sample data correction |
CN104463783A (en) * | 2014-10-31 | 2015-03-25 | 杭州美诺瓦医疗科技有限公司 | Processing and displaying method for controlling image through local multi-parameter and multi-picture shortcut key |
CN104574361A (en) * | 2014-11-27 | 2015-04-29 | 沈阳东软医疗***有限公司 | Image processing method and device for mammary peripheral tissue equalization |
CN105894547A (en) * | 2016-05-06 | 2016-08-24 | 南昌航空大学 | Image processing method based on group-wave transformation compressed sensing |
CN107730476A (en) * | 2017-09-27 | 2018-02-23 | 温州大学 | A kind of image enchancing method and device based on compressed sensing |
CN107862666A (en) * | 2017-11-22 | 2018-03-30 | 新疆大学 | Mixing Enhancement Methods about Satellite Images based on NSST domains |
CN111507912A (en) * | 2020-04-08 | 2020-08-07 | 深圳市安健科技股份有限公司 | Mammary gland image enhancement method and device and computer readable storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101639537A (en) * | 2009-09-04 | 2010-02-03 | 西安电子科技大学 | SAR image noise suppression method based on direction wave domain mixture Gaussian model |
CN101719268A (en) * | 2009-12-04 | 2010-06-02 | 西安电子科技大学 | Generalized Gaussian model graph denoising method based on improved Directionlet region |
EP2232446A1 (en) * | 2007-12-20 | 2010-09-29 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction |
-
2011
- 2011-04-19 CN CN 201110098272 patent/CN102142133B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2232446A1 (en) * | 2007-12-20 | 2010-09-29 | Wisconsin Alumni Research Foundation | Method for prior image constrained image reconstruction |
CN101639537A (en) * | 2009-09-04 | 2010-02-03 | 西安电子科技大学 | SAR image noise suppression method based on direction wave domain mixture Gaussian model |
CN101719268A (en) * | 2009-12-04 | 2010-06-02 | 西安电子科技大学 | Generalized Gaussian model graph denoising method based on improved Directionlet region |
Non-Patent Citations (2)
Title |
---|
张新生,高新波,王颖,张士杰: "乳腺X线图像的增强与噪声抑制研究", 《红外与毫米波学报》, vol. 29, no. 1, 28 February 2010 (2010-02-28) * |
戴琼海: "压缩感知研究", 《计算机学报》, vol. 34, no. 3, 31 March 2011 (2011-03-31) * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102663709B (en) * | 2012-04-28 | 2014-09-03 | 大连海事大学 | Vessel enhancement method for X-ray angiogram |
CN102663709A (en) * | 2012-04-28 | 2012-09-12 | 大连海事大学 | Vessel enhancement method for X-ray angiogram |
CN102903099A (en) * | 2012-09-05 | 2013-01-30 | 西安电子科技大学 | Color image edge detection method based on directionlet conversion |
CN103310429A (en) * | 2013-03-06 | 2013-09-18 | 西安电子科技大学 | Image enhancement method based on hidden Markov tree (HMT) model in directionlet domain |
CN103310429B (en) * | 2013-03-06 | 2015-05-27 | 西安电子科技大学 | Image enhancement method based on hidden Markov tree (HMT) model in directionlet domain |
CN103473797B (en) * | 2013-09-16 | 2016-04-20 | 电子科技大学 | Spatial domain based on compressed sensing sampling data correction can downscaled images reconstructing method |
CN103473797A (en) * | 2013-09-16 | 2013-12-25 | 电子科技大学 | Image space-domain reducible reconstructing method based on compressed sensing sample data correction |
CN104463783A (en) * | 2014-10-31 | 2015-03-25 | 杭州美诺瓦医疗科技有限公司 | Processing and displaying method for controlling image through local multi-parameter and multi-picture shortcut key |
CN104574361A (en) * | 2014-11-27 | 2015-04-29 | 沈阳东软医疗***有限公司 | Image processing method and device for mammary peripheral tissue equalization |
CN104574361B (en) * | 2014-11-27 | 2017-12-29 | 沈阳东软医疗***有限公司 | A kind of mammary gland peripheral tissues balanced image processing method and device |
CN105894547A (en) * | 2016-05-06 | 2016-08-24 | 南昌航空大学 | Image processing method based on group-wave transformation compressed sensing |
CN107730476A (en) * | 2017-09-27 | 2018-02-23 | 温州大学 | A kind of image enchancing method and device based on compressed sensing |
CN107862666A (en) * | 2017-11-22 | 2018-03-30 | 新疆大学 | Mixing Enhancement Methods about Satellite Images based on NSST domains |
CN111507912A (en) * | 2020-04-08 | 2020-08-07 | 深圳市安健科技股份有限公司 | Mammary gland image enhancement method and device and computer readable storage medium |
CN111507912B (en) * | 2020-04-08 | 2023-03-24 | 深圳市安健科技股份有限公司 | Mammary gland image enhancement method and device and computer readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN102142133B (en) | 2013-01-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102142133A (en) | Mammary X-ray image enhancement method based on non-subsampled Directionlet transform and compressive sensing | |
Ayres et al. | Characterization of architectural distortion in mammograms | |
Sahiner et al. | Computer‐aided detection of clustered microcalcifications in digital breast tomosynthesis: a 3D approach | |
CN104616255B (en) | Self-adapting enhancement method based on mammography X | |
ITBO20000166A1 (en) | METHOD, AND RELATED EQUIPMENT, FOR THE AUTOMATIC DETECTION OF MICROCALCIFICATIONS IN DIGITAL SIGNALS OF BREAST FABRIC. | |
Nusantara et al. | Classification of digital mammogram based on nearest-neighbor method for breast cancer detection | |
CN104182755A (en) | Mammary gland molybdenum target X-ray image block feature extraction method based on tower-shaped principal component analysis (PCA) | |
Jones et al. | Improving mammography lesion classification by optimal fusion of handcrafted and deep transfer learning features | |
CN111784701A (en) | Ultrasonic image segmentation method and system combining boundary feature enhancement and multi-scale information | |
Beheshti et al. | Classification of abnormalities in mammograms by new asymmetric fractal features | |
Ning et al. | Multi-scale gradational-order fusion framework for breast lesions classification using ultrasound images | |
CN102682439B (en) | Medical image fusion method based on multidirectional empirical mode decomposition | |
Ganvir et al. | Filtering method for pre-processing mammogram images for breast cancer detection | |
Azevedo-Marques et al. | Classification of vertebral compression fractures in magnetic resonance images using spectral and fractal analysis | |
Fan et al. | Magnetic-resonance image segmentation based on improved variable weight multi-resolution Markov random field in undecimated complex wavelet domain | |
Zhang et al. | A novel automatic tumor detection for breast cancer ultrasound images | |
Das et al. | Multimodal classification on PET/CT image fusion for lung cancer: a comprehensive survey | |
Al-Absi et al. | Computer aided diagnosis system based on machine learning techniques for lung cancer | |
CN104331864B (en) | Based on the processing of the breast image of non-down sampling contourlet and the significant model of vision | |
Kim et al. | Breast mass detection using slice conspicuity in 3D reconstructed digital breast volumes | |
Oliver et al. | Automatic diagnosis of masses by using level set segmentation and shape description | |
CN109377461A (en) | A kind of breast X-ray image self-adapting enhancement method based on NSCT | |
Al-Azzawi et al. | An efficient medical image fusion method using contourlet transform based on PCM | |
Makandar et al. | Mammography image enhancement using linear, nonlinear and wavelet filters with histogram equalization | |
dos Santos Romualdo et al. | Mammography images restoration by quantum noise reduction and inverse MTF filtering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130123 Termination date: 20180419 |
|
CF01 | Termination of patent right due to non-payment of annual fee |