CN109949233A - Filter out the method, system, equipment and storage medium of scattered rays in radioscopic image - Google Patents
Filter out the method, system, equipment and storage medium of scattered rays in radioscopic image Download PDFInfo
- Publication number
- CN109949233A CN109949233A CN201910123586.3A CN201910123586A CN109949233A CN 109949233 A CN109949233 A CN 109949233A CN 201910123586 A CN201910123586 A CN 201910123586A CN 109949233 A CN109949233 A CN 109949233A
- Authority
- CN
- China
- Prior art keywords
- image
- scattered rays
- original image
- details
- rays
- 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
Abstract
This application provides a kind of method, system, equipment and storage medium for filtering out scattered rays in radioscopic image, include the following steps: to extract the scattered rays characteristics of image in original image；The scattered rays iconic element in the original image is removed according to the scattered rays characteristics of image；Original image after removal scattered rays iconic element is subjected to details reduction treatment, obtains details also original image；The gray value of the details also original image is adjusted to state up to standard, acquisition filters out the radioscopic image after scattered rays.The interference of scattered rays can be removed while retaining the original details in image to the maximum extent by carrying out details reduction treatment by the original image after removal scattered rays iconic element and removal scattered rays iconic element；By choosing specified extracting parameter k, b, c to obtain with the good direct projection line image of contrast；Noise reduction process is completely avoided in complete image processing process, greatly remains the information in radioscopic image.
Description
Technical field
This application involves field of computer technology, filters out the method for scattered rays in radioscopic image more particularly to a kind of, is
System, equipment and storage medium.
Background technique
X-ray has very strong penetrability, can make to be placed in the light reaching the film after object, according to Beer-Lambert law, X
With exponential form energy attenuation can occur for ray, this is the theoretical basis of x-ray imaging during incident medium.In reality
In the process, the microcosmic interaction of multiple complexity can occur for X-ray and medium.According to exercising result, it can substantially be classified as two aspects:
On the one hand lose ray energy, i.e. energy attenuation；On the other hand change directions of rays, i.e. generation scattered rays.For imaging
For, energy attenuation is conducive to improve picture quality (beam hardening effect etc. is put aside), on the contrary, scattered rays then will lead to figure
As quality degradation, image low contrast resolution and spatial resolution decline are shown as, to influence clinical diagnosis, in some instances it may even be possible to
It causes to fail to pinpoint a disease in diagnosis.
It is influenced to eliminate scattered rays, traditional method is placed between human body and detector using mechanical grid device, right
Line is directly scattered across the radioparent of human body to filter out.However, this method drawback is obvious, it is larger not require nothing more than exposure dose,
Grid shadow can also be brought.Another way can effectively overcome these problems, i.e., without using mechanical grid, and use pure image procossing
Means, i.e., digital grid are scattered line inhibition processing to original image.Currently, existing number grid technology is generally adopted
Enhance image detail with multiresolution gain process, it is multiple dimensioned therebetween under the gain of multiplying property can multiple amplification noise, and most
The amplification noise under each scale is superimposed again in whole image reconstruction, noise problem is further aggravated.On the other hand, too strong
Noise reduction process also inevitably brings potential information loss side effect.
Summary of the invention
In view of the above problems, it proposes the embodiment of the present invention and overcomes the above problem or at least partly in order to provide one kind
A kind of method, system, equipment and storage medium for filtering out scattered rays in radioscopic image to solve the above problems.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of method for filtering out scattered rays in radioscopic image, packets
Include following steps:
Extract the scattered rays characteristics of image in original image；
The scattered rays iconic element in above-mentioned original image is removed according to above-mentioned scattered rays characteristics of image；
Original image after removal scattered rays iconic element is subjected to details reduction treatment, obtains details also original image；
The gray value of above-mentioned details also original image is adjusted to state up to standard, acquisition filters out the radioscopic image after scattered rays.
Optionally, the step of extracting the scattered rays characteristics of image in original image, comprising:
Above-mentioned original image is decomposed into areas of dielectric and non-areas of dielectric；
Scattered rays characteristics of image in Extraction medium region and non-areas of dielectric respectively.
The step of optional scattered rays characteristics of image distinguished in Extraction medium region and non-areas of dielectric, comprising:
Obtain the difference of the minimum relative intensity of direct rays and direct rays highest and minimum relative intensity in above-mentioned areas of dielectric
Different parameter, and according to the difference of the minimum relative intensity of direct rays and direct rays highest and minimum relative intensity in above-mentioned areas of dielectric
Different parameter obtains the scattered rays characteristics of image in above-mentioned areas of dielectric；
Obtain in above-mentioned non-areas of dielectric the minimum relative intensity of direct rays and according to direct rays in above-mentioned non-areas of dielectric most
Low relative intensity obtains the scattered rays characteristics of image in above-mentioned non-areas of dielectric.
Optionally, the scattered rays iconic element in above-mentioned original image is calculated according to following equation group:
Ist=I-Isc
Wherein, I indicates original image, IstIndicate direct rays iconic element, IscIndicate scattered rays iconic element, ΩmIt indicates
Areas of dielectric in image,Indicate image in non-areas of dielectric, k indicate areas of dielectric in direct rays highest with it is minimum relatively strong
The difference parameter of degree, b indicate the minimum relative intensity of direct rays in areas of dielectric, and c indicates that direct rays are relatively strong in non-areas of dielectric
Degree.
Optionally, the original image after removal scattered rays iconic element is subjected to details reduction treatment, obtains details reduction
The step of image, comprising:
Original image after above-mentioned removal scattered rays iconic element is subjected to picture breakdown processing, obtains picture breakdown number
According to；
Image detail in above-mentioned picture breakdown data is subjected to gain amplification disposal, the picture breakdown number after obtaining gain
According to；
Picture breakdown data after above-mentioned gain are subjected to image reconstruction, obtain above-mentioned details also original image.
Optionally, it in the gray value for adjusting above-mentioned details also original image to state up to standard, obtains the X after filtering out scattered rays and penetrates
Before the step of line image, further includes:
Judge whether the gray value of above-mentioned details also original image respectively reaches storage request and display requires respectively；
If it is not, the gray value for then adjusting above-mentioned details also original image obtains to state up to standard and filters out the X-ray after scattered rays
Image.
Optionally, the gray value of above-mentioned details also original image is adjusted to state up to standard, and acquisition filters out the X-ray after scattered rays
The step of image, comprising:
The gray scale domain of above-mentioned details also original image is adjusted to state up to standard by LUT mapping processing.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of system for filtering out scattered rays in radioscopic image, packets
Include following specific module:
Extraction module, for extracting the scattered rays characteristics of image in original image；
Remove module, for according to above-mentioned scattered rays characteristics of image remove the scattering line image in above-mentioned original image at
Point；
Recovery module carries out details reduction treatment for that will remove the original image after scattered rays iconic element, obtains thin
Section goes back original image；
Module is adjusted, the gray value for adjusting above-mentioned details also original image obtains after filtering out scattered rays to state up to standard
Radioscopic image.
To solve the above-mentioned problems, the embodiment of the invention discloses a kind of computer equipment, including memory, processor with
And the computer program that can be run on a memory and on a processor is stored, the processor realizes this when executing described program
The method of scattered rays in radioscopic image is filtered out described in any one of inventive embodiments.
To solve the above-mentioned problems, it the embodiment of the invention discloses a kind of computer readable storage medium, is stored thereon with
Computer program is realized when the program is executed by processor and filters out radioscopic image described in any one of embodiment of the present invention
The method of middle scattered rays.
Compared with prior art, the application includes following advantages:
In the embodiment of the present invention, pass through the original image after removal scattered rays iconic element and removal scattered rays iconic element
The interference of scattered rays can be removed while retaining the original details in image to the maximum extent by carrying out details reduction treatment；It is logical
It crosses and chooses specified extracting parameter k, b, c to obtain with the good direct projection line image of contrast；By choosing specified high frequency
Signal gain parameter makes the detail recovery of image work well；It is completely avoided at noise reduction in complete image processing process
Reason, greatly remains the information in radioscopic image.
Detailed description of the invention
Fig. 1 is the step flow diagram for filtering out the method for scattered rays in radioscopic image of one embodiment of the invention；
Fig. 2 is the step flow diagram for filtering out the method for scattered rays in radioscopic image of one embodiment of the invention；
Fig. 3 is the die body original image schematic diagram of specific embodiments of the present invention；
Fig. 4 is the die body scattered rays characteristic image schematic diagram of specific embodiments of the present invention；
Fig. 5 is the die body original image schematic diagram after the removal scattered rays characteristics of image of specific embodiments of the present invention；
Fig. 6 is the die body details also original image schematic diagram of specific embodiments of the present invention；
Fig. 7 is the die body X-ray image schematic diagram of the method for specific embodiments of the present invention filtered out after scattered rays；
Fig. 8 is the modular structure schematic diagram for filtering out the device of scattered rays in radioscopic image of one embodiment of the invention；
Fig. 9 is a kind of structural schematic diagram of computer equipment of one embodiment of the invention.
1, extraction module；2, module is removed；3, recovery module；4, module is adjusted；12, computer equipment；14, outside is set
It is standby；16, processing unit；18, bus；20, network adapter；22, (I/O) interface；24, display；28, system storage；30,
Random access memory (RAM)；32, cache memory；34, storage system；40, program/utility；42, program mould
Block.
The embodiments will be further described with reference to the accompanying drawings for the realization, the function and the advantages of the object of the present invention.
Specific embodiment
In order to make the above objects, features, and advantages of the present application more apparent, with reference to the accompanying drawing and it is specific real
Applying mode, the present application will be further described in detail.
Obviously, described embodiment is only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
It in addition, the technical solution between each embodiment can be combined with each other, but must be with ordinary skill
Based on personnel can be realized, this technical side will be understood that when the combination of technical solution appearance is conflicting or cannot achieve
The combination of case is not present, also not the present invention claims protection scope within.
Finally, all the embodiments in this specification are described in a progressive manner, what each embodiment stressed
It is the difference from other embodiments, the same or similar parts between the embodiments can be referred to each other.
Referring to Fig.1, a kind of step process for the embodiment of the method 1 for filtering out scattered rays in radioscopic image of the application is shown
Figure, can specifically include following steps:
Scattered rays characteristics of image in S1, extraction original image；
S2, scattered rays iconic element in above-mentioned original image is removed according to above-mentioned scattered rays characteristics of image；
S3, the original image after removal scattered rays iconic element is subjected to details reduction treatment, obtains details also original image；
S4, the gray value for adjusting above-mentioned details also original image obtain to state up to standard and filter out the X ray picture after scattered rays
Picture.
As described in above-mentioned steps S1, the scattered rays characteristics of image in original image is extracted, above-mentioned original image generally arrives
Up to the synthesis quantum image I of detector, wherein above-mentioned comprehensive quantum image I generally comprises direct projection quantum characteristics of image Ist, scattering
Quantum characteristics of image IscWith quantum noise characteristics of image In, it should be noted that due to quantum noise characteristics of image InFor X-ray
Therefore the non-uniform Distribution of photon and power spectrum disperse do not consider individually in embodiments of the present invention, and be included into direct projection quantum together
Characteristics of image IstWith scattered quantum characteristics of image IscIn, therefore, the original image described in various embodiments of the present invention is (that is, comprehensive
Quantum image I) characteristics of image only includes direct projection quantum characteristics of image IstWith scattered quantum characteristics of image Isc。
It should be noted that said extracted process, which is generally required, is split as areas of dielectric and non-dielectric area for original image
Domain, then the scattered rays characteristics of image in two regions is extracted respectively.
As described in above-mentioned steps S2, the scattering line image in above-mentioned original image is removed according to above-mentioned scattered rays characteristics of image
Ingredient, in the present embodiment, above-mentioned minimizing technology are generally the pixel value for modifying corresponding region, and concrete modification is will be above-mentioned original
Image subtracts the scattered rays image pixel being made of above-mentioned scattered rays characteristics of image and obtains, for the figure after removal scattered rays ingredient
As pixel value range, do not adjust.
It include: by being deleted the pixel coincideing with the scattered rays characteristics of image or being adjusted to the pixel value of the pixel
0 or 255, to obtain the original image after removing scattered rays iconic element；By by with the scattered rays characteristics of image coincide
The pixel value of pixel adjusts the original obtained after removal scattered rays iconic element as the identical pixel value of near zone pixel value
Beginning image.
As described in above-mentioned steps S3, the original image after removal scattered rays iconic element is subjected to details reduction treatment, is obtained
Obtaining details also original image, it should be noted that above-mentioned details reduction process generally comprises three steps, specifically: firstly, being gone above-mentioned
Except the original image after scattered rays iconic element carries out picture breakdown processing；It then, will be specified in the image after resolution process
Parts of images details carries out gain amplification disposal；Finally, the image after it will carry out details gain amplification disposal carries out image weight
Structure processing, obtains above-mentioned details also original image.
As described in above-mentioned steps S4, the gray value of above-mentioned details also original image is adjusted to state up to standard, acquisition filters out scattering
Radioscopic image after line, it should be noted that the gray value of above-mentioned details also original image to state up to standard generally comprises two aspects
Judgement: whether reach the memory requirement of X-ray image and whether reach the display requirement of X-ray image, above-mentioned details also original image
The gray value of picture needs to reach memory requirement and display requires just to meet transition up to standard.
Referring to Fig. 2, a kind of step process for the embodiment of the method 2 for filtering out scattered rays in radioscopic image of the application is shown
Figure, can specifically include following steps:
In embodiments of the present invention, the step of extracting the scattered rays characteristics of image in original image, comprising:
S11, above-mentioned original image is decomposed into areas of dielectric and non-areas of dielectric；
S12, respectively the scattered rays characteristics of image in Extraction medium region and non-areas of dielectric.
As described in above-mentioned steps S11, above-mentioned original image is decomposed into areas of dielectric and non-areas of dielectric, wherein medium
Region is the tissue regions in above-mentioned original image, and non-areas of dielectric is the non-tissue regions in above-mentioned original image.
As described in above-mentioned steps S12, the scattered rays characteristics of image in Extraction medium region and non-areas of dielectric, needs respectively
Illustrate, obtained by following sub-step:
Sub-step S12-1, obtain in above-mentioned areas of dielectric the minimum relative intensity of direct rays and direct rays highest with it is minimum
The difference parameter of relative intensity, and according to the minimum relative intensity of direct rays in above-mentioned areas of dielectric and direct rays highest with it is minimum
The difference parameter of relative intensity obtains the scattered rays characteristics of image in above-mentioned areas of dielectric；
Sub-step S12-2, the minimum relative intensity of direct rays is obtained in above-mentioned non-areas of dielectric and according to above-mentioned non-dielectric area
The minimum relative intensity of direct rays obtains the scattered rays characteristics of image in above-mentioned non-areas of dielectric in domain.
Sub-step S12-3, in actual operation calculates above-mentioned original image areas of dielectric and non-Jie according to following equation group
Scattered rays iconic element in matter region:
Ist=I-Isc
Wherein, I indicates original image, IstIndicate direct rays iconic element, IscIndicate scattered rays iconic element, ΩmIt indicates
Areas of dielectric in image,Indicate image in non-areas of dielectric, k indicate areas of dielectric in direct rays highest with it is minimum relatively strong
The difference parameter of degree, b indicate the minimum relative intensity of direct rays in areas of dielectric, and c indicates that direct rays are relatively strong in non-areas of dielectric
Degree.
In embodiments of the present invention, the original image after removal scattered rays iconic element is subjected to details reduction treatment, obtained
Details also original image the step of, comprising:
S31, the original image after above-mentioned removal scattered rays iconic element is subjected to picture breakdown processing, obtains picture breakdown
Data；
S32, the image detail in above-mentioned picture breakdown data is subjected to gain amplification disposal, the image after obtaining gain point
Solve data；
S33, the picture breakdown data after above-mentioned gain are subjected to image reconstruction, obtain above-mentioned details also original image.
It should be noted that resolution ratio is higher, and corresponding frequency band is higher, such as in Laplacian pyramid: when Decomposition order is
At 8 layers, first layer is identical as original image resolution, frequency highest；8th layer of resolution ratio is minimum, and frequency is minimum；Therein
Two layers then successively reduce to layer 7 resolution ratio and frequency.
It should be noted that high band described in various embodiments of the present invention and low-frequency range are non-limiting for relative value
Value, that is, primary completely in image processing process, current highest frequency range and peak low band, the as secondary image procossing mistake
The high band and low-frequency range of journey.
Similarly, high-frequency signal and low frequency signal described in various embodiments of the present invention are also relative value and non-limiting values,
That is, high-frequency signal represents finer textual details picture signal in same tomographic image；Low frequency signal is indicated compared with thick edge
Picture signal, such as: bone edges.
As described in above-mentioned steps S31, the original image after above-mentioned removal scattered rays iconic element is carried out at picture breakdown
Reason obtains picture breakdown data, it should be noted that and above-mentioned picture breakdown processing is generally Laplacian pyramid decomposition, point
Solution mode generally comprises direct sub-sample resolution, Gaussian filtering is decomposed and mean filter decomposes etc., wherein direct sub-sample resolution
It is specific as follows:
(1) using original image as first layer input picture, the drop that resolution ratio halves directly is carried out to first layer input picture
Sampling obtains second layer input picture；
(2) it is identical as first layer input picture to rise sampling acquisition for the interpolation for carrying out double resolution to second layer input picture
The forecast image of resolution ratio, wherein interpolation method includes in arest neighbors interpolation, bilinear interpolation, bi-cubic interpolation and Gauss
It inserts, specific method is depending on required precision and practical problem；
(3) error image is acquired according to first layer input picture, forecast image, error image is Laplacian gold word
Tower first layer (base)；
Error image=first layer input picture-forecast image
(4) continue to repeat the above process with jth layer input picture (j=2,3 ..., n) for input picture, be finally
Column error image (is opened, j=2,3 ..., n) comprising j, i.e. each tomographic image of Laplacian pyramid, finally enables (n+1)th differential chart
As=the (n+1)th layer of input picture (top layer), to obtain n+1 layers of Laplacian pyramid.
Gaussian filtering is decomposed specific as follows:
Gaussian filtering decompose with direct sub-sample resolution only difference is that, obtain the side of approximate jth layer input picture j
Formula is different: direct sub-sample resolution is that the direct resolution ratio that carries out halves down-sampled, and Gaussian filtering is decomposed and then first carried out
Then Gaussian filtering again halves filtered image progress resolution ratio down-sampled.As remaining process the two, therefore not
It repeats again.
Mean filter decomposes specific as follows:
Unique difference is decomposed with Gaussian filtering to be to generate the difference of filter used in approximate jth layer input picture,
As remaining process the two, therefore repeat no more.
Specific isolation is depending on the required precision of image and problem weighting, as long as final reduction result can't produce
The isolation of raw difference substantially can be used as the isolation of any embodiment of that present invention, wherein pyramid decomposition
The number of plies is generally 4 layers or more, preferably 4~8 layers, certainly, the number of plies can be reduced under the premise of not influencing detail recovery effect to 4
Layer or less or increase the number of plies to 8 layers or more；
As described in above-mentioned steps S32, the image detail in above-mentioned picture breakdown data is subjected to gain amplification disposal, is obtained
Picture breakdown data after gain, wherein during carrying out to image detail gain, the image detail frequency of suffered gain
Section includes the full range segment data in above-mentioned picture breakdown data, the gain ranges of the data of each frequency range by the frequency range curvature control
Parameter processed determines, wherein is to be positively correlated each other between the curvature control parameter, frequency range and enhancing intensity, that is, frequency range gets over higher curvature
Control parameter is bigger, and enhancing intensity is bigger.
In practical applications, the image detail of specified portions in above-mentioned picture breakdown data is subjected to gain amplification disposal,
Specially Sigmoid will be successively carried out by the picture breakdown data that Laplacian pyramid decomposition obtains in above mentioned step S3 1
Curve transformation, specific formula is as follows:
σ (x)=1/ (1+e-x)
In formula:
(a, b) is Sigmoid curve symmetric centre coordinate；M is output codomain control parameter, for controlling Sigmoid song
Line output amplitude range；C is Sigmoid curvature of curve control parameter, for controlling signal amplification dynamics.
As described in above-mentioned steps S33, the picture breakdown data after above-mentioned gain are subjected to image reconstruction, obtain above-mentioned details
Also original image, it should be noted that above-mentioned image reconstruction process is generally Laplacian Pyramid Reconstruction, it should be noted that
Reconstruct mode is the reverse process process of isolation, that is, the reconstruct mode implemented in this step is in same step embodiment
The reverse process of isolation used by above mentioned step S3 1 in the process, it should be noted that alternative abovementioned steps if it exists
The reconstruct mode obtained after isolation used by S31 is reverse can also be as the reconstruct mode of this step.
In embodiments of the present invention, in the gray value for adjusting above-mentioned details also original image to state up to standard, acquisition filters out scattered
Before the step of radioscopic image after ray, further includes:
S41, judge whether the gray value of above-mentioned details also original image respectively reaches storage request and display requires respectively；
S42, if it is not, then adjusting the gray value of above-mentioned details also original image to state up to standard, acquisition filters out the X after scattered rays
Ray image.
S41, judge whether the gray value of above-mentioned details also original image respectively reaches storage request and display requires, and needs respectively
It is noted that only have in this judgment step when storage request and show that judging result is yes when requiring while reaching standard,
That is, another is not reached when requiring or when storage request and display are wanted when any one of storage request and display requirement reach requirement
It asks and is not up to when requiring, judging result is no；
S42, if it is not, then adjusting the gray value of above-mentioned details also original image to state up to standard, acquisition filters out the X after scattered rays
Ray image；S43, if so, exporting above-mentioned details also original image for the above-mentioned radioscopic image filtered out after scattered rays, wherein
When determine result be it is no, that is, it is general to pass through mapping processing progress gray scale when needing to adjust the gray value of above-mentioned details also original image
The adjustment of value generally comprises the screening to pixel before adjustment, to filter out the region for needing to carry out gray value adjustment.
In practical applications, general to be adjusted in the gray scale domain of above-mentioned details also original image to up to standard by LUT mapping processing
State.
Referring to Fig. 2-7, in this embodiment, to include the radioscopic image of several pieces of organic glass die bodys (as schemed
Shown in 3) for, wherein the thickness of the organic glass die body is incremented by successively with ROI1-ROI7, hereinafter referred to as die body X ray picture
Picture:
Inputted that (x, y is pixels in image using the die body radioscopic image as die body original image Img_org (x, y)
Point position coordinates), the specific steps are as follows:
It inputs die body original image Img_org (x, y), according to following equation, from die body original image Img_org (x, y)
Middle extraction scattered rays characteristics of image:
In formula, Img_sc (x, y) is die body scattered rays characteristics of image；K indicate areas of dielectric in direct rays highest with it is minimum
The difference parameter of relative intensity, b indicate the minimum relative intensity of direct rays in areas of dielectric, and c indicates direct rays in non-areas of dielectric
Relative intensity；T=(c-b)/k is the gray threshold for distinguishing medium and non-areas of dielectric；
Scattered rays ingredient is reduced from original image by following equation:
Img_st (x, y)=Img_org (x, y)-Img_sc (x, y)
In formula, Img_st (x, y) is the die body original image removed after scattered rays ingredient；
Laplacian pyramid decomposition is carried out to the die body original image Img_st (x, y) after removal scattered rays ingredient,
In this specific embodiment, is decomposed using Gaussian filtering and image is decomposed:
A, with Img_st (x, y) for first layer input picture G1, and to G1The filtering of 5 × 5 Gaussian kernels is carried out, is then divided
The down-sampled acquisition second layer input picture G that resolution halves2；
B, to G2The interpolation for carrying out double resolution rises sampling acquisition and G1The forecast image P of equal resolution1, interpolation side
Method includes but is not limited to: arest neighbors interpolation value method, bilinear interpolation value method, bi-cubic interpolation value method and Gaussian interpolation value method, root
Depending on required precision and practical problem, Gaussian interpolation value method is selected in this embodiment；
C, according to G1、P1Seek out error image L1, wherein L1=G1-P1, L1For Laplacian pyramid first layer (base
Layer)；
With Gj(j=2,3 ..., n) repeats the above steps a-c for input picture, to obtain serial error image Lj(j=2,
3 ..., n), i.e. each tomographic image of Laplacian pyramid, finally enable Ln+1=Gn+1(top layer), to obtain n+1 layers
Laplacian pyramid
Following Sigmoid curve transformation is executed to each layer of Laplacian pyramid, with also each level of detail of basic stitch:
In formula, Lin(x, y) is input Laplacian pyramidal layer image, Lout(x, y) is output Laplacian pyramid
Tomographic image；σ (x)=1/ (1+e-x)；Parameter is Sigmoid conversion curve symmetrical centre coordinate to (a, b)；M is output codomain control
Parameter processed controls Sigmoid conversion curve output amplitude range；C is Sigmoid curvature of curve control parameter, controls each level
Details amplifies dynamics.
Laplacian pyramid is reconstructed, to obtain the die body details also original image Img_r (x, y) after details reduction,
In, restructuring procedure is the inverse process of pyramid decomposition
The gray scale domain of Img_r (x, y) is adjusted to optimum value range by LUT mapping mode, and exports final result figure
Picture, that is, the die body X-ray image after filtering out scattered rays.
For system embodiments, since it is basically similar to the method embodiment, related so being described relatively simple
Place illustrates referring to the part of embodiment of the method.
Referring to Fig. 3, a kind of step process for filtering out the system embodiment of scattered rays in radioscopic image of the invention is shown
Figure, can specifically include following module:
Extraction module 1, for extracting the scattered rays characteristics of image in original image；
Remove module 2, for according to above-mentioned scattered rays characteristics of image remove the scattering line image in above-mentioned original image at
Point；
Recovery module 3 carries out details reduction treatment for that will remove the original image after scattered rays iconic element, obtains thin
Section goes back original image；
Module 4 is adjusted, the gray value for adjusting above-mentioned details also original image obtains after filtering out scattered rays to state up to standard
Radioscopic image.
Said extracted module 1 is generally used for extracting the scattered rays characteristics of image in original image, and above-mentioned original image is general
For the synthesis quantum image I for reaching detector, wherein above-mentioned comprehensive quantum image I generally comprises direct projection quantum characteristics of image Ist、
Scattered quantum characteristics of image IscWith quantum noise characteristics of image In, it should be noted that due to quantum noise characteristics of image InFor X
Therefore the non-uniform Distribution of ray photons and power spectrum disperse do not consider individually in embodiments of the present invention, and be included into direct projection together
Quantum characteristics of image IstWith scattered quantum characteristics of image IscIn, therefore, original image described in various embodiments of the present invention (that is,
Comprehensive quantum image I) characteristics of image only includes direct projection quantum characteristics of image IstWith scattered quantum characteristics of image Isc。
It should be noted that said extracted process, which is generally required, is split as areas of dielectric and non-dielectric area for original image
Domain, then the scattered rays characteristics of image in two regions is extracted respectively.
Above-mentioned removal module 2 is generally used for removing the scattering in above-mentioned original image according to above-mentioned scattered rays characteristics of image
Line image ingredient, in the present embodiment, above-mentioned removal generally comprise two methods and are eliminated, comprising: by that will dissipate with described
The identical pixel of ray image feature deletes or the pixel value of the pixel is adjusted to 0 or 255, removes scattering line image to obtain
Original image after ingredient；It is adjusted by the pixel value for the pixel that will be coincide with the scattered rays characteristics of image as near zone
The identical pixel value of pixel value is to obtain the original image after removal scattered rays iconic element.
Above-mentioned recovery module 3 is generally used for that the original image progress details also original place after scattered rays iconic element will be removed
Reason obtains details also original image, it should be noted that and above-mentioned details reduction process generally comprises three steps, specifically: firstly, will
Original image after above-mentioned removal scattered rays iconic element carries out picture breakdown processing；It then, will be in the image after resolution process
Specified portions image detail carry out gain amplification disposal；Finally, the image after it will carry out details gain amplification disposal carries out
Image reconstruction process obtains above-mentioned details also original image.
Above-mentioned adjustment module 4 is generally used for adjusting the gray value of above-mentioned details also original image to state up to standard, is filtered out
Radioscopic image after scattered rays, it should be noted that the gray value of above-mentioned details also original image to state up to standard generally comprises two
The judgement of aspect: whether reach the memory requirement of X-ray image and whether reach the display requirement of X-ray image, above-mentioned details is also
The gray value of original image needs to reach memory requirement and display requires just to meet transition up to standard.
Referring to Fig. 4, show it is of the invention it is a kind of realize described in filter out the computer of the method for scattered rays in radioscopic image
Equipment can specifically include as follows:
Above-mentioned computer equipment 12 is showed in the form of universal computing device, the component of computer equipment 12 may include but
Be not limited to: one or more processor or processing unit 16, system storage 28, connecting different system components (including is
Unite memory 28 and processing unit 16) bus 18.
Bus 18 indicates one of a few 18 structures of class bus or a variety of, including memory bus 18 or memory control
Device, peripheral bus 18, graphics acceleration port, processor or the office using 18 structure of any bus in a variety of 18 structures of bus
Domain bus 18.For example, these architectures include but is not limited to industry standard architecture (ISA) bus 18, microchannel
Architecture (MAC) bus 18, enhanced isa bus 18, audio-video frequency electronic standard association (VESA) local bus 18 and outer
Enclose component interconnection (PCI) bus 18.
Computer equipment 12 typically comprises a variety of computer system readable media.These media can be it is any can be by
The usable medium that computer equipment 12 accesses, including volatile and non-volatile media, moveable and immovable medium.
System storage 28 may include the computer system readable media of form of volatile memory, such as arbitrary access
Memory (RAM) 30 and/or cache memory 32.Computer equipment 12 may further include other movement/it is not removable
Dynamic, volatile/non-volatile computer decorum storage medium.Only as an example, storage system 34 can be used for read and write can not
Mobile, non-volatile magnetic media (commonly referred to as " hard disk drive ").Although not shown in fig 4, it can provide for can
The disc driver of mobile non-volatile magnetic disk (such as " floppy disk ") read-write, and to removable anonvolatile optical disk (such as CD~
ROM, DVD~ROM or other optical mediums) read-write CD drive.In these cases, each driver can pass through one
A or multiple data mediums interface is connected with bus 18.Memory may include at least one program product, the program product
With one group of (for example, at least one) program module 42, these program modules 42 are configured to perform the function of various embodiments of the present invention
Energy.
Program/utility 40 with one group of (at least one) program module 42, can store in memory, for example,
Such program module 42 includes --- but being not limited to --- operating system, one or more application program, other program moulds
It may include the realization of network environment in block 42 and program data, each of these examples or certain combination.Program mould
Block 42 usually executes function and/or method in embodiment described in the invention.
Computer equipment 12 can also with one or more external equipments 14 (such as keyboard, sensing equipment, display 24,
Camera etc.) communication, the equipment interacted with the computer equipment 12 can be also enabled a user to one or more to be communicated, and/
Or with enable the computer equipment 12 and one or more other calculate any equipment that equipment are communicated (such as network interface card,
Modem etc.) communication.This communication can be carried out by interface input/output (I/O) 22.Also, computer equipment
12 can also by network adapter 20 and one or more network (such as local area network (LAN)), wide area network (WAN) and/or
Public network (such as internet) communication.As shown, network adapter 20 passes through other of bus 18 and computer equipment 12
Module communication.It should be understood that although not shown in fig 4, other hardware and/or software mould can be used in conjunction with computer equipment 12
Block, including but not limited to: microcode, device driver, redundant processing unit 16, external disk drive array, RAID system, magnetic
Tape drive and data backup storage system 34 etc..
Processing unit 16 by the program that is stored in system storage 28 of operation, thereby executing various function application and
Data processing, such as realize the method that scattered rays in radioscopic image is filtered out provided by the embodiment of the present invention.
That is, above-mentioned processing unit 16 is realized when executing above procedure: extracting the scattered rays characteristics of image in original image；
The scattered rays iconic element in above-mentioned original image is removed according to above-mentioned scattered rays characteristics of image；Scattered rays iconic element will be removed
Original image afterwards carries out details reduction treatment, obtains details also original image；Adjust the gray value of above-mentioned details also original image extremely
State up to standard, acquisition filter out the radioscopic image after scattered rays.
In embodiments of the present invention, the present invention also provides a kind of computer readable storage medium, it is stored thereon with computer
Program, realized when which is executed by processor as the application all embodiments offers filter out scattered rays in radioscopic image
Method:
That is, realization when being executed by processor to program: extracting the scattered rays characteristics of image in original image；According to above-mentioned
Scattered rays characteristics of image removes the scattered rays iconic element in above-mentioned original image；It will be original after removal scattered rays iconic element
Image carries out details reduction treatment, obtains details also original image；The gray value of above-mentioned details also original image is adjusted to state up to standard,
Acquisition filters out the radioscopic image after scattered rays.
It can be using any combination of one or more computer-readable media.Computer-readable medium can be calculating
Machine gram signal media or computer readable storage medium.Computer readable storage medium for example can be --- but it is unlimited
In system, device or the device of --- electricity, magnetic, optical, electromagnetic, infrared ray or semiconductor, or any above combination.Computer
The more specific example (non exhaustive list) of readable storage medium storing program for executing includes: electrical connection with one or more conducting wires, portable
Formula computer disk, hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory
(EPOM or flash memory), optical fiber, portable compact disc read-only memory (CD~ROM), light storage device, magnetic memory device or
Above-mentioned any appropriate combination.In this document, computer readable storage medium can be it is any include or storage program
Tangible medium, the program can be commanded execution system, device or device use or in connection.
Computer-readable signal media may include in a base band or as carrier wave a part propagate data-signal,
Wherein carry computer-readable program code.The data-signal of this propagation can take various forms, including --- but
It is not limited to --- electromagnetic signal, optical signal or above-mentioned any appropriate combination.Computer-readable signal media can also be
Any computer-readable medium other than computer readable storage medium, change computer-readable medium can send, propagate or
Transmission is for by the use of instruction execution system, device or device or program in connection.
The computer for executing operation of the present invention can be write with one or more programming languages or combinations thereof
Program code, above procedure design language include object oriented program language --- such as Java, Smalltalk, C+
+, further include conventional procedural programming language --- such as " C " language or similar programming language.Program code
It can fully execute on the user computer, partly execute, held as an independent software package on the user computer
Part executes on the remote computer or holds on a remote computer or server completely on the user computer for row, part
Row.In situations involving remote computers, remote computer can pass through the network of any kind --- including local area network
(LAN) or wide area network (WAN) --- it is connected to subscriber computer, or, it may be connected to outer computer (such as using because of spy
Service provider is netted to connect by internet).
In any of the above-described inventive embodiments, after removal scattered rays iconic element and removal scattered rays iconic element
Original image, which carries out details reduction treatment, to remove scattered rays while retaining the original details in image to the maximum extent
Interference；By choosing specified extracting parameter k, b, c to obtain with the good direct projection line image of contrast；Referred to by choosing
Fixed high-frequency signal gain parameter makes the detail recovery of image work well；It is completely avoided in complete image processing process
Noise reduction process greatly remains the information in radioscopic image.
It a kind of filter out the method, system, equipment of scattered rays and storage in radioscopic image to provided herein above and is situated between
Matter is described in detail, and specific examples are used herein to illustrate the principle and implementation manner of the present application, above
The explanation of embodiment is merely used to help understand the present processes and its core concept；Meanwhile for the general skill of this field
Art personnel, according to the thought of the application, there will be changes in the specific implementation manner and application range, in conclusion this
Description should not be construed as the limitation to the application.
Claims (10)
1. a kind of method for filtering out scattered rays in radioscopic image, which comprises the steps of:
Extract the scattered rays characteristics of image in original image；
The scattered rays iconic element in the original image is removed according to the scattered rays characteristics of image；
Original image after removal scattered rays iconic element is subjected to details reduction treatment, obtains details also original image；
The gray value of the details also original image is adjusted to state up to standard, acquisition filters out the radioscopic image after scattered rays.
2. the method according to claim 1, wherein extracting the step of the scattered rays characteristics of image in original image
Suddenly, comprising:
The original image is decomposed into areas of dielectric and non-areas of dielectric；
Scattered rays characteristics of image in Extraction medium region and non-areas of dielectric respectively.
3. according to the method described in claim 2, it is characterized in that, the scattering in Extraction medium region and non-areas of dielectric respectively
The step of line image feature, comprising:
The difference of the minimum relative intensity of direct rays and direct rays highest and minimum relative intensity in the areas of dielectric is obtained to join
Number, and joined according to the difference of the minimum relative intensity of direct rays in the areas of dielectric and direct rays highest and minimum relative intensity
Number obtains the scattered rays characteristics of image in the areas of dielectric；
Obtain in the non-areas of dielectric the minimum relative intensity of direct rays and according to the minimum phase of direct rays in the non-areas of dielectric
Scattered rays characteristics of image in the non-areas of dielectric is obtained to intensity.
4. the method according to claim 1, wherein calculating dissipating in the original image according to following equation group
Ray image ingredient:
Ist=I-Isc
Wherein, I indicates original image, IstIndicate direct rays iconic element, IscIndicate scattered rays iconic element, ΩmIndicate image
Middle areas of dielectric,Indicate non-areas of dielectric in image, k indicates direct rays highest and minimum relative intensity in areas of dielectric
Difference parameter, b indicate the minimum relative intensity of direct rays in areas of dielectric, and c indicates direct rays relative intensity in non-areas of dielectric.
5. the method according to claim 1, wherein the original image after removal scattered rays iconic element is carried out
Details reduction treatment, obtain details also original image the step of, comprising:
Original image after the removal scattered rays iconic element is subjected to picture breakdown processing, obtains picture breakdown data；
Image detail in described image decomposition data is subjected to gain amplification disposal, the picture breakdown data after obtaining gain；
Picture breakdown data after the gain are subjected to image reconstruction, obtain the details also original image.
6. the method according to claim 1, wherein in the gray value for adjusting the details also original image to up to standard
State, before the step of obtaining the radioscopic image filtered out after scattered rays, further includes:
Judge whether the gray value of the details also original image respectively reaches storage request and display requires respectively；
If it is not, the gray value for then adjusting the details also original image obtains to state up to standard and filters out the X ray picture after scattered rays
Picture.
7. the method according to claim 1, wherein adjusting the gray value of the details also original image to shape up to standard
State, the step of obtaining the radioscopic image filtered out after scattered rays, comprising:
The gray scale domain of the details also original image is adjusted to state up to standard by LUT mapping processing.
8. a kind of system for filtering out scattered rays in radioscopic image, which is characterized in that including following specific module:
Extraction module, for extracting the scattered rays characteristics of image in original image；
Module is removed, for removing the scattered rays iconic element in the original image according to the scattered rays characteristics of image；
Recovery module carries out details reduction treatment for that will remove the original image after scattered rays iconic element, obtains details also
Original image；
Module is adjusted, the gray value for adjusting the details also original image obtains the X after filtering out scattered rays and penetrate to state up to standard
Line image.
9. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage
Computer program, which is characterized in that the processor is realized when executing described program such as any one of claim 1~7 institute
The method stated.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is by processor
The method as described in any one of claim 1~7 is realized when execution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910123586.3A CN109949233B (en) | 2019-02-18 | 2019-02-18 | Method, system, device and storage medium for filtering scattered rays in X-ray image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910123586.3A CN109949233B (en) | 2019-02-18 | 2019-02-18 | Method, system, device and storage medium for filtering scattered rays in X-ray image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109949233A true CN109949233A (en) | 2019-06-28 |
| CN109949233B CN109949233B (en) | 2022-12-13 |
Family
ID=67006765
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910123586.3A Active CN109949233B (en) | 2019-02-18 | 2019-02-18 | Method, system, device and storage medium for filtering scattered rays in X-ray image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109949233B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111462004A (en) * | 2020-03-30 | 2020-07-28 | 北京推想科技有限公司 | Image enhancement method and device, computer equipment and storage medium |
| CN113134171A (en) * | 2021-04-23 | 2021-07-20 | 上海联影医疗科技股份有限公司 | Positioning method of radiotherapy quality control measuring point and computer equipment |
Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5648997A (en) * | 1995-12-29 | 1997-07-15 | Advanced Optical Technologies, Inc. | Apparatus and method for removing scatter from an x-ray image |
| CN101030298A (en) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | Method for enhancing medical image with multi-scale self-adaptive contrast change |
| CN101109718A (en) * | 2006-11-14 | 2008-01-23 | 北京国药恒瑞美联信息技术有限公司 | Virtual grid imaging method and system used for eliminating influence of scattered radiation |
| WO2010058329A1 (en) * | 2008-11-21 | 2010-05-27 | Philips Intellectual Property & Standards Gmbh | System and method for x-ray scatter correction |
| CN101779962A (en) * | 2010-01-19 | 2010-07-21 | 西安华海医疗信息技术股份有限公司 | Method for enhancing medical X-ray image display effect |
| CN101939763A (en) * | 2007-12-10 | 2011-01-05 | 爱克发医疗保健公司 | Method of generating a multiscale contrast enhanced image |
| US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
| CN102737395A (en) * | 2011-04-15 | 2012-10-17 | 深圳迈瑞生物医疗电子股份有限公司 | Method and apparatus for image processing in medical X-ray system |
| CN102793551A (en) * | 2011-05-24 | 2012-11-28 | 柯尼卡美能达医疗印刷器材株式会社 | Chest diagnostic support information generation system |
| CN103079466A (en) * | 2010-08-27 | 2013-05-01 | 柯尼卡美能达医疗印刷器材株式会社 | Thoracic diagnosis assistance system and program |
| US20130223717A1 (en) * | 2010-06-08 | 2013-08-29 | Zakrytoe Akcionernoe Obshchestvo "Impul's" | Method for Correction of Digital Images |
| CN103578082A (en) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | Cone beam CT scatter correction method and system |
| CN104299187A (en) * | 2013-07-15 | 2015-01-21 | 上海联影医疗科技有限公司 | Image processing method and image processing device for DR (Digital Radiography) system |
| CN104915930A (en) * | 2014-03-14 | 2015-09-16 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for carrying out gray scale compensation and noise suppression on image |
| CN106033598A (en) * | 2015-03-18 | 2016-10-19 | 北京纳米维景科技有限公司 | Digital filter wire grid imaging method used for eliminating scattered radiation influence |
| US20170055933A1 (en) * | 2014-03-24 | 2017-03-02 | Fujifilm Corporation | Radiographic image processing device, method, and recording medium |
| CN106504222A (en) * | 2016-11-21 | 2017-03-15 | 河海大学常州校区 | A kind of Polarization Image Fusion system under water based on bionic visual mechanism |
| CN106570836A (en) * | 2016-11-03 | 2017-04-19 | 北京航星机器制造有限公司 | Radiation image local enhancement method |
| CN107811647A (en) * | 2017-11-24 | 2018-03-20 | 上海联影医疗科技有限公司 | CT equipment, the X-ray detection X method with reference to detection device and radiographic source |
| CN108227029A (en) * | 2017-12-26 | 2018-06-29 | 北京无线电计量测试研究所 | A kind of mm-wave imaging compensation method, equipment and readable computer storage medium |
| CN108447028A (en) * | 2018-02-12 | 2018-08-24 | 中国科学院西安光学精密机械研究所 | Underwater picture quality improving method based on Multiscale Fusion |
| CN108577872A (en) * | 2017-12-29 | 2018-09-28 | 上海联影医疗科技有限公司 | Medical image removes scattering method, system and storage medium |
| CN109345608A (en) * | 2018-10-22 | 2019-02-15 | 中国人民解放军战略支援部队信息工程大学 | A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal |
-
2019
- 2019-02-18 CN CN201910123586.3A patent/CN109949233B/en active Active
Patent Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5648997A (en) * | 1995-12-29 | 1997-07-15 | Advanced Optical Technologies, Inc. | Apparatus and method for removing scatter from an x-ray image |
| CN101109718A (en) * | 2006-11-14 | 2008-01-23 | 北京国药恒瑞美联信息技术有限公司 | Virtual grid imaging method and system used for eliminating influence of scattered radiation |
| CN101030298A (en) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | Method for enhancing medical image with multi-scale self-adaptive contrast change |
| CN101939763A (en) * | 2007-12-10 | 2011-01-05 | 爱克发医疗保健公司 | Method of generating a multiscale contrast enhanced image |
| WO2010058329A1 (en) * | 2008-11-21 | 2010-05-27 | Philips Intellectual Property & Standards Gmbh | System and method for x-ray scatter correction |
| CN101779962A (en) * | 2010-01-19 | 2010-07-21 | 西安华海医疗信息技术股份有限公司 | Method for enhancing medical X-ray image display effect |
| US20130223717A1 (en) * | 2010-06-08 | 2013-08-29 | Zakrytoe Akcionernoe Obshchestvo "Impul's" | Method for Correction of Digital Images |
| CN103079466A (en) * | 2010-08-27 | 2013-05-01 | 柯尼卡美能达医疗印刷器材株式会社 | Thoracic diagnosis assistance system and program |
| US20120207370A1 (en) * | 2010-12-20 | 2012-08-16 | Benjamin Pooya Fahimian | Systems and Methods for Simultaneous Acquisition of Scatter and Image Projection Data in Computed Tomography |
| CN102737395A (en) * | 2011-04-15 | 2012-10-17 | 深圳迈瑞生物医疗电子股份有限公司 | Method and apparatus for image processing in medical X-ray system |
| CN102793551A (en) * | 2011-05-24 | 2012-11-28 | 柯尼卡美能达医疗印刷器材株式会社 | Chest diagnostic support information generation system |
| CN103578082A (en) * | 2012-08-09 | 2014-02-12 | 江苏超惟科技发展有限公司 | Cone beam CT scatter correction method and system |
| CN104299187A (en) * | 2013-07-15 | 2015-01-21 | 上海联影医疗科技有限公司 | Image processing method and image processing device for DR (Digital Radiography) system |
| CN104915930A (en) * | 2014-03-14 | 2015-09-16 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for carrying out gray scale compensation and noise suppression on image |
| US20170055933A1 (en) * | 2014-03-24 | 2017-03-02 | Fujifilm Corporation | Radiographic image processing device, method, and recording medium |
| CN106033598A (en) * | 2015-03-18 | 2016-10-19 | 北京纳米维景科技有限公司 | Digital filter wire grid imaging method used for eliminating scattered radiation influence |
| CN106570836A (en) * | 2016-11-03 | 2017-04-19 | 北京航星机器制造有限公司 | Radiation image local enhancement method |
| CN106504222A (en) * | 2016-11-21 | 2017-03-15 | 河海大学常州校区 | A kind of Polarization Image Fusion system under water based on bionic visual mechanism |
| CN107811647A (en) * | 2017-11-24 | 2018-03-20 | 上海联影医疗科技有限公司 | CT equipment, the X-ray detection X method with reference to detection device and radiographic source |
| CN108227029A (en) * | 2017-12-26 | 2018-06-29 | 北京无线电计量测试研究所 | A kind of mm-wave imaging compensation method, equipment and readable computer storage medium |
| CN108577872A (en) * | 2017-12-29 | 2018-09-28 | 上海联影医疗科技有限公司 | Medical image removes scattering method, system and storage medium |
| CN108447028A (en) * | 2018-02-12 | 2018-08-24 | 中国科学院西安光学精密机械研究所 | Underwater picture quality improving method based on Multiscale Fusion |
| CN109345608A (en) * | 2018-10-22 | 2019-02-15 | 中国人民解放军战略支援部队信息工程大学 | A kind of pyramidal CT image method for reconstructing based on asymmetric scatter removal |
Non-Patent Citations (3)
| Title |
|---|
| KYUNG SANG KIM ET AL: "Ultra-fast hybrid CPU-GPU Monte Carlo simulation for scatter correction in 3D PETs", 《2011 IEEE NUCLEAR SCIENCE SYMPOSIUM CONFERENCE RECORD》 * |
| 刘建邦: "锥束CT图像散射伪影校正方法研究", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 * |
| 胡勇: "一种基于图像的非均匀单散射参与介质建模方法", 《中国科学:信息科学》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111462004A (en) * | 2020-03-30 | 2020-07-28 | 北京推想科技有限公司 | Image enhancement method and device, computer equipment and storage medium |
| CN111462004B (en) * | 2020-03-30 | 2023-03-21 | 推想医疗科技股份有限公司 | Image enhancement method and device, computer equipment and storage medium |
| CN113134171A (en) * | 2021-04-23 | 2021-07-20 | 上海联影医疗科技股份有限公司 | Positioning method of radiotherapy quality control measuring point and computer equipment |
| CN113134171B (en) * | 2021-04-23 | 2022-12-30 | 上海联影医疗科技股份有限公司 | Positioning method of radiotherapy quality control measuring point and computer equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109949233B (en) | 2022-12-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huang et al. | Metal artifact reduction on cervical CT images by deep residual learning | |
| Qian et al. | Tree structured wavelet transform segmentation of microcalcifications in digital mammography | |
| CN107481297B (en) | CT image reconstruction method based on convolutional neural network | |
| Chun et al. | MRI super‐resolution reconstruction for MRI‐guided adaptive radiotherapy using cascaded deep learning: In the presence of limited training data and unknown translation model | |
| Na et al. | Guided filter‐based images fusion algorithm for CT and MRI medical images | |
| Liu et al. | Deep learning with noise‐to‐noise training for denoising in SPECT myocardial perfusion imaging | |
| Zhu et al. | Lesion focused super-resolution | |
| CN104574316B (en) | A kind of digitized X-ray machine system | |
| Ma et al. | Sinogram denoising via attention residual dense convolutional neural network for low-dose computed tomography | |
| Bhairannawar | Efficient medical image enhancement technique using transform HSV space and adaptive histogram equalization | |
| CN109949233A (en) | Filter out the method, system, equipment and storage medium of scattered rays in radioscopic image | |
| Li et al. | Incorporation of residual attention modules into two neural networks for low‐dose CT denoising | |
| CN106530236A (en) | Medical image processing method and system | |
| Xia et al. | Endoscopic image enhancement with noise suppression | |
| CN110070510A (en) | A kind of CNN medical image denoising method for extracting feature based on VGG-19 | |
| Wang et al. | Adaptive decomposition method for multi‐modal medical image fusion | |
| Li et al. | An adaptive self-guided wavelet convolutional neural network with compound loss for low-dose CT denoising | |
| Zhang et al. | Medical image fusion based on quasi-cross bilateral filtering | |
| Zhu et al. | STEDNet: Swin transformer‐based encoder–decoder network for noise reduction in low‐dose CT | |
| Xu et al. | Patient-specific hyperparameter learning for optimization-based CT image reconstruction | |
| Li et al. | Unsupervised data fidelity enhancement network for spectral CT reconstruction | |
| Zhang et al. | Medical image fusion and noise suppression with fractional‐order total variation and multi‐scale decomposition | |
| Zhu et al. | MRI enhancement based on visual-attention by adaptive contrast adjustment and image fusion | |
| Zhou et al. | Image illumination adaptive correction algorithm based on a combined model of bottom-hat and improved gamma transformation | |
| Li et al. | Medical image fusion based on local Laplacian decomposition and iterative joint filter |
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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: No.103, baguang District Service Center, No.2 BaiShaWan Road, baguang community, Kuiyong street, Dapeng New District, Shenzhen, Guangdong 518000 Applicant after: Shenzhen Lanying Medical Technology Co.,Ltd. Address before: 518000 Shenzhen, Baoan District, Guangdong Shiyan street, Lang Xin neighborhood committee, brick factory village, Liyuan Industrial Zone, hung lung hi tech Industrial Park (blue Yun Industrial Park), 1 1 buildings. Applicant before: SHENZHEN LANYUN MEDICAL IMAGE CO.,LTD. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |