CN106803241A - The processing method and processing device of angiographic image - Google Patents
The processing method and processing device of angiographic image Download PDFInfo
- Publication number
- CN106803241A CN106803241A CN201710041232.5A CN201710041232A CN106803241A CN 106803241 A CN106803241 A CN 106803241A CN 201710041232 A CN201710041232 A CN 201710041232A CN 106803241 A CN106803241 A CN 106803241A
- Authority
- CN
- China
- Prior art keywords
- image
- angiographic
- frisket
- gray value
- pixel
- 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.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 230000002708 enhancing effect Effects 0.000 claims abstract description 29
- 230000000052 comparative effect Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 14
- 230000009467 reduction Effects 0.000 claims description 18
- 238000003708 edge detection Methods 0.000 claims description 10
- 239000008280 blood Substances 0.000 claims description 8
- 210000004369 blood Anatomy 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 18
- 210000004204 blood vessel Anatomy 0.000 description 11
- 238000003745 diagnosis Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 5
- 230000036541 health Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 210000000056 organ Anatomy 0.000 description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000001225 therapeutic effect Effects 0.000 description 2
- 238000002583 angiography Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 210000004209 hair Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a kind of processing method and processing device of angiographic image, the method includes：With the full picture as reference picture,Image registration is carried out to the frisket image,Obtain the frisket image after registration,The frisket image for obtaining and the full picture are carried out to subtract shadow treatment,Obtain the first angiographic image,And,Frisket image and the full picture after the registration is carried out to subtract shadow treatment,Obtain the second angiographic image,The gray value summation of the pixel in the gray value summation of the pixel in preset regions in first angiographic image and the preset regions in second angiographic image is compared,And first angiographic image or second angiographic image are chosen as target image according to comparative result,Image definition enhancing treatment is carried out to the target image,Obtain result images,And the result images are shown,Which enhance the identification of image,Improve the effect of the image for subtracting movie queen.
Description
Technical field
The invention belongs to image processing field, more particularly to a kind of processing method and processing device of angiographic image.
Background technology
Digital subtraction angiography is widely used diagnosis and treatment in vascular diseases, and its cardinal principle is to remove and blood
The background structure that pipe is overlapped.Because doctor needs that patient is diagnosed and treated by subtracting the image of movie queen, so subtracting movie queen
The definition of image be particularly important.The method of blood vessel substraction has various, is broadly divided into temporal subtraction, energy subtraction and mixes
Conjunction subtracts shadow.
In the prior art, generally using Digital Subtraction blood-vessel image processing system carries out blood vessel substraction to image, conventional
Subtract image method for temporal subtraction.Digital Subtraction blood-vessel image processing system (does not have the image that initial several width do not have blood vessel to manifest
Have injection contrast agent, referred to as frisket), and image that several follow-up blood vessels manifest (injection contrast agent, referred to as fill piece) subtracts each other,
The image of blood vessel is obtained comprising only, while manifesting blood vessel influence dynamic change situation over time in real time.But by
In external factor, such as noise and motion blur, the image for subtracting movie queen can be caused not clear enough, and then influence diagnosis.
The content of the invention
The present invention provides a kind of processing method and processing device of angiographic image, it is intended to solve what is caused due to external factor
The image for subtracting movie queen is not clear enough, and then influences the problem of diagnosis.
A kind of processing method of angiographic image that the present invention is provided, including：Frisket image and full picture are obtained,
And with the full picture as reference picture, image registration is carried out to the frisket image, obtain the frisket after registration
Image；The frisket image for obtaining and the full picture are carried out to subtract shadow treatment, the first angiographic image is obtained, with
And, the frisket image and the full picture after the registration are carried out subtracting shadow treatment, obtain the second angiogram
Picture；By the gray value summation of the pixel in preset regions in first angiographic image and second angiogram
The gray value summation of the pixel as described in preset regions is compared, and chooses first blood vessel according to comparative result
Contrastographic picture or second angiographic image are used as target image；The target image is carried out at image definition enhancing
Reason, obtains result images, and the result images are shown.
A kind of processing unit of angiographic image that the present invention is provided, including：Acquisition module, for obtaining frisket image
With full picture；Image registration module, for the full picture as reference picture, figure being carried out to the frisket image
As registration, the frisket image after registration is obtained；Photographic Subtraction module, for the frisket image for obtaining and described filling
The picture that is full of carries out subtracting shadow treatment, obtains the first angiographic image, and, to the frisket image after the registration and institute
Stating full picture carries out subtracting shadow treatment, obtains the second angiographic image；Computing module, for by first angiogram
In the gray value summation of the pixel in image in preset regions and preset regions described in second angiographic image
The gray value summation of pixel is compared, and chooses first angiographic image or second blood according to comparative result
Pipe contrastographic picture is used as target image；Image processing module, for carrying out image definition enhancing treatment to the target image,
Result images are obtained, and the result images are shown.
The processing method and processing device of the angiographic image that the present invention is provided, obtains frisket image and full picture, and
With the full picture as reference picture, image registration is carried out to the frisket image, the frisket image after registration is obtained, to obtaining
The frisket image and the full picture for taking carry out subtracting shadow treatment, obtain the first angiographic image, and, after the registration
The frisket image and the full picture subtract shadow treatment, the second angiographic image is obtained, by first angiogram
Pixel in the gray value summation of the pixel in image in preset regions and second angiographic image in the preset regions
The gray value summation of point is compared, and chooses first angiographic image or second angiogram according to comparative result
Picture shows to the result images as target image, compared to prior art, increases by image registration and definition
After the reason of strength, the identification of image is improve, improve the effect of the image for subtracting movie queen, and then effective auxiliary doctor to patient
Diagnosis, be beneficial to the health of patient.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
The accompanying drawing to be used needed for having technology description is briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention.
Fig. 1 is that the processing method of the angiographic image that first embodiment of the invention is provided realizes schematic flow sheet；
Fig. 2 is that the processing method of the angiographic image that second embodiment of the invention is provided realizes schematic flow sheet；
Fig. 3 is the schematic diagram of frisket image provided in an embodiment of the present invention；
Fig. 4 is the schematic diagram of full picture provided in an embodiment of the present invention；
Fig. 5 is the schematic diagram of the first angiographic image provided in an embodiment of the present invention；
Fig. 6 is the schematic diagram of the second angiographic image provided in an embodiment of the present invention；
Fig. 7 is the partial schematic diagram of the first angiographic image provided in an embodiment of the present invention；
Fig. 8 is that the independent part by enhanced first angiographic image in image border provided in an embodiment of the present invention is shown
It is intended to；
Fig. 9 is the partial schematic diagram after the independent noise reduction to the first angiographic image provided in an embodiment of the present invention；
Figure 10 is the partial schematic diagram after the independent image enhaucament to the first angiographic image provided in an embodiment of the present invention；
Figure 11 is the structural representation of the processing unit of the angiographic image that third embodiment of the invention is provided；
Figure 12 is the structural representation of the processing unit of the angiographic image that fourth embodiment of the invention is provided.
Specific embodiment
To enable that goal of the invention of the invention, feature, advantage are more obvious and understandable, below in conjunction with the present invention
Accompanying drawing in embodiment, is clearly and completely described to the technical scheme in the embodiment of the present invention, it is clear that described reality
It is only a part of embodiment of the invention to apply example, and not all embodiments.Based on the embodiment in the present invention, people in the art
The every other embodiment that member is obtained under the premise of creative work is not made, belongs to the scope of protection of the invention.
Refer to Fig. 1, Fig. 1 realizes that flow is shown for what first embodiment of the invention provided the processing method of angiographic image
It is intended to, can be applied on the imaging device of therapeutic treatment or diagnosis, such as C-arm X-ray machine, the angiographic image shown in Fig. 1
Processing method, mainly includes the following steps that：
S101, acquisition frisket image and full picture, and with the full picture as reference picture, to the frisket image
Image registration is carried out, the frisket image after registration is obtained.
Easily cause the deformation of organ due to breathing or movement, frisket image can be caused by image registration and filled
Characteristic point on picture reaches matching.
S102, the frisket image and the full picture to acquisition carry out subtracting shadow treatment, obtain the first angiogram
Picture, and, the frisket image and the full picture after the registration are carried out to subtract shadow treatment, obtain the second angiogram
Picture.
The step has been done and subtract twice shadow treatment, and one is that the frisket image and full picture before registration subtract at shadow
Reason, two is that the full picture before the frisket image after registration and registration is carried out to subtract shadow treatment.
S103, by the gray value summation of the pixel in preset regions in first angiographic image and second blood vessel
The gray value summation of the pixel in contrastographic picture in the preset regions is compared, and chooses first blood according to comparative result
Pipe contrastographic picture or second angiographic image are used as target image.
Because image border part can interfere to result of calculation, thus choose first angiographic image and this second
The rule of the preset regions is identical in angiographic image, i.e., the preset regions are to remove the region of image border part.
S104, image definition enhancing treatment is carried out to the target image, obtain result images, and the result images are entered
Row display.
There may be noise and/or fuzzy by the target image obtained after step S103, by image
Definition enhancing treatment can eliminate the noise and/or fuzzy situation of target image, improve the definition of the target image.
In the embodiment of the present invention, frisket image and full picture are obtained, and with the full picture as reference picture, it is right
The frisket image carries out image registration, obtains the frisket image after registration, to the frisket image and the full piece figure for obtaining
As carrying out subtracting shadow treatment, the first angiographic image is obtained, and, to the frisket image and the full picture after the registration
Carry out subtracting shadow treatment, the second angiographic image is obtained, by the pixel in preset regions in first angiographic image
The gray value summation of the pixel in gray value summation and the preset regions in second angiographic image is compared, and according to
First angiographic image or second angiographic image are chosen as target image according to comparative result, to the target image
Image definition enhancing treatment is carried out, result images are obtained, and the result images are shown, compared to prior art, warp
After crossing image registration and definition enhancing treatment, the identification of image is improve, improve the effect of the image for subtracting movie queen, entered
And diagnosis of the doctor to patient is effectively aided in, it is beneficial to the health of patient.
Refer to Fig. 2, the processing method of the angiographic image that Fig. 2 is provided for second embodiment of the invention realizes flow
Schematic diagram, can be applied on the imaging device of therapeutic treatment or diagnosis, such as C-arm X-ray machine, the angiographic image shown in Fig. 2
Processing method, mainly include the following steps that：
S201, acquisition frisket image and full picture, and with the full picture as reference picture, to the frisket image
Image registration is carried out, the frisket image after registration is obtained.
Easily cause the deformation of organ due to breathing or movement, frisket image can be caused by image registration and filled
Characteristic point on picture reaches matching.As shown in Figure 3 and Figure 4, the schematic diagram of frisket image, Fig. 4 before Fig. 3 is image registration
To fill the schematic diagram of picture before image registration.
S202, the frisket image and the full picture to acquisition carry out subtracting shadow treatment, obtain the first angiogram
Picture, and, the frisket image and the full picture after the registration are carried out to subtract shadow treatment, obtain the second angiogram
Picture.
The step has been done and subtract twice shadow treatment, and one is that the frisket image and full picture before registration subtract at shadow
Reason, two is that the full picture before the frisket image after registration and registration is carried out to subtract shadow treatment.As shown in Figure 5 and Figure 6,
Fig. 5 is the schematic diagram of the first angiographic image, and Fig. 6 is the schematic diagram of the second angiographic image.
S203, by the gray value summation of the pixel in preset regions in first angiographic image and second blood vessel
The gray value summation of the pixel in contrastographic picture in the preset regions is compared, and chooses first blood according to comparative result
Pipe contrastographic picture or second angiographic image are used as target image.
Because image border part can interfere to result of calculation, thus choose first angiographic image and this second
The rule of the preset regions is identical in angiographic image, i.e., the preset regions are to remove the region of image border part.
Further, by the gray value summation of the pixel in preset regions in first angiographic image and this second
The gray value summation of the pixel in angiographic image in the preset regions is compared, and according to comparative result choose this
One angiographic image or second angiographic image are specially as target image：
The gray value summation of the pixel in first angiographic image in the preset regions is calculated, and, this second
The gray value summation of the pixel in angiographic image in the preset regions；
It is corresponding with second angiographic image in the gray value summation of the corresponding pixel of the first angiographic image
Pixel gray value summation in choose the small angiographic image of numerical value as the target image.
If the gray value summation of the corresponding pixel of the first angiographic image is less than second angiographic image pair
The gray value summation of the pixel answered, then choose the first angiographic image as the target image；Conversely, then choosing the second blood
Pipe contrastographic picture is used as the target image.The numerical value of the gray value summation is smaller, represents through over-subtraction shadow, eliminates most non-
The image of angiosomes.
S204, the edge in the target image medium vessels region is detected according to preset edge detection operator, and according to
Preset weighted value carries out enhancing treatment to the edge of the angiosomes after detection, obtains the enhanced target image in edge.
The preset edge detection operator is not limited, and the preset edge detection operator can be Sobel operators, it is also possible to
Laplacian operators, can also be the edge detection operators such as Canny operators.The selection of weighted value can be according to real image
Enhancing effect chosen, the preset weighted value is bigger, and the enhancing effect to image border is stronger, and the preset weighted value is got over
Small, the enhancing effect to image border is weaker.As shown in Figure 7 and Figure 8, Fig. 7 is the partial schematic diagram of the first angiographic image,
Fig. 8 is partial schematic diagram individually by enhanced first angiographic image in image border.It should be noted that Fig. 7 and Fig. 8
It is in order to simply to illustrate that the difference of the same position image before and after the enhancing of image border.
S205, the target image enhanced to edge carry out noise remove and image enhaucament, obtain the result images, and
The result images are shown.
Further, the target image enhanced to edge carries out noise remove and image enhaucament, obtains the result figure
As being specially：
Noise is carried out according to integration noise reduction, frequency domain noise reduction or Recursive Filtering target image enhanced to the edge
Remove, obtain image to be reinforced；
Grey scale curve stretching is carried out to the image to be reinforced；
Using image to be reinforced after stretching as the result images.
As shown in Figure 9 and Figure 10, Fig. 9 is to the partial schematic diagram after the independent noise reduction of the first angiographic image.Figure 10 is
To the partial schematic diagram after the independent image enhaucament of the first angiographic image.The process of grey scale curve stretching is because curve is oblique
The maximum region of rate just falls the region that gray value is concentrated in the image to be reinforced, and being stretched by grey scale curve can be by the ash
The region that angle value is concentrated is stretched to bigger scope.
In actual applications, the mode of the removal noise is not only comprising integration noise reduction, frequency domain noise reduction or Recursive Filtering, ability
Field technique personnel are according to the mode for removing noise for integrating noise reduction, frequency domain noise reduction or Recursive Filtering described by the embodiment of the present invention
Associate the mode of other removal noises.
Alternatively, above step may operate in CUDA (Compute Unified Device Architecture)
Computing platform on, to ensure the speed of image procossing.
It should be noted that the result schematic diagram in order to more intuitively find out each stage image procossing, Fig. 8, Fig. 9 and Tu
10 is respectively individually to carry out a result schematic diagram for stage image procossing.
In the embodiment of the present invention, frisket image and full picture are obtained, and with the full picture as reference picture, it is right
The frisket image carries out image registration, obtains the frisket image after registration, to the frisket image and the full piece figure for obtaining
As carrying out subtracting shadow treatment, the first angiographic image is obtained, and, to the frisket image and the full picture after the registration
Carry out subtracting shadow treatment, the second angiographic image is obtained, by the pixel in preset regions in first angiographic image
The gray value summation of the pixel in gray value summation and the preset regions in second angiographic image is compared, and according to
First angiographic image or second angiographic image are chosen as target image according to comparative result, to the target image
Image definition enhancing treatment is carried out, result images are obtained, and the result images are shown, compared to prior art, warp
After crossing image registration and definition enhancing treatment, the identification of image is improve, improve the effect of the image for subtracting movie queen, entered
And diagnosis of the doctor to patient is effectively aided in, it is beneficial to the health of patient.
Figure 11 is referred to, Figure 11 is that the structure of the processing unit of the angiographic image that third embodiment of the invention is provided is shown
It is intended to, for convenience of description, illustrate only the part related to the embodiment of the present invention.The place of the angiographic image of Figure 11 examples
Reason device can be the executive agent of the processing method of the angiographic image that earlier figures 1 and embodiment illustrated in fig. 2 are provided.Figure 11
The processing unit of the angiographic image of example, mainly includes：Acquisition module 1101, image registration module 1102, Photographic Subtraction
Module 1103, computing module 1104 and image processing module 1105.Each functional module describes in detail as follows above：
Acquisition module 1101, for obtaining frisket image and full picture；
Image registration module 1102, for the full picture as reference picture, to the frisket image and the full piece
Image carries out image registration, obtains the frisket image after registration；
Photographic Subtraction module 1103, carries out subtracting shadow treatment for the frisket image and the full picture to obtaining, and obtains
To the first angiographic image, and, the frisket image and the full picture after the registration are carried out to subtract shadow treatment, obtain
Second angiographic image；
Computing module 1104, for by the gray value summation of the pixel in preset regions in first angiographic image
Gray value summation with the pixel in the preset regions in second angiographic image is compared, and according to comparative result
First angiographic image or second angiographic image are chosen as target image；
Image processing module 1105, for carrying out image definition enhancing treatment to the target image, obtains result images,
And the result images are shown.
Easily cause the deformation of organ due to breathing or movement, frisket image can be caused by image registration and filled
Characteristic point on picture reaches matching.Because image border part can interfere to result of calculation, therefore choose first blood
The rule of the preset regions is identical in pipe contrastographic picture and second angiographic image, i.e., the preset regions are removal figure
As the region of marginal portion.
The present embodiment details not to the greatest extent, refers to the description of foregoing embodiment illustrated in fig. 1, and here is omitted.
It should be noted that in the implementation method of the processing unit of the angiographic image of figure 11 above example, each function
The division of module is merely illustrative of, can as needed in practical application, such as the configuration requirement or software of corresponding hardware
Realization convenient consideration, and above-mentioned functions distribution is completed by different functional modules, will said apparatus internal structure
Different functional modules are divided into, to complete all or part of function described above.And, in practical application, this implementation
Corresponding functional module in example can be realized by corresponding hardware, it is also possible to which it is complete to perform corresponding software by corresponding hardware
Into.Each embodiment that this specification is provided can all apply foregoing description principle, below repeat no more.
In the embodiment of the present invention, acquisition module 1101 obtains frisket image and full picture, image registration module 1102
With the full picture as reference picture, image registration is carried out to the frisket image and the full picture, after obtaining registration
The frisket image, the frisket image and the full picture of the 1103 pairs of acquisitions of Photographic Subtraction module carry out subtracting shadow treatment, obtain
First angiographic image, and, to the frisket image and the full picture after the registration subtract shadow treatment, obtain the
Two angiographic images, computing module 1104 is total by the gray value of the pixel in preset regions in first angiographic image
And be compared with the gray value summation of pixel in the preset regions in second angiographic image, and according to comparing knot
Fruit chooses first angiographic image or second angiographic image as target image, and image processing module 1105 pairs should
Target image carries out image definition enhancing treatment, obtains result images, and the result images are shown, compared to existing
Technology, by after image registration and definition enhancing treatment, improve the identification of image, improves the image that subtracts movie queen
Effect, and then diagnosis of the effective auxiliary doctor to patient, are beneficial to the health of patient.
Figure 12 is referred to, the structure of the processing unit of the angiographic image that Figure 12 is provided for fourth embodiment of the invention is shown
It is intended to, for convenience of description, illustrate only the part related to the embodiment of the present invention.The place of the angiographic image of Figure 12 examples
Reason device can be the executive agent of the processing method of the angiographic image that earlier figures 1 and embodiment illustrated in fig. 2 are provided.Figure 12
The processing unit of the angiographic image of example, mainly includes：Acquisition module 1201, image registration module 1202, Photographic Subtraction
Module 1203, computing module 1204, image processing module 1205, wherein computing module 1204 include：The He of calculating sub module 12041
Choose submodule 12042.Each functional module describes in detail as follows above：
Acquisition module 1201, for obtaining frisket image and full picture.
Image registration module 1202, for the full picture as reference picture, carrying out image to the frisket image and matching somebody with somebody
Standard, obtains the frisket image after registration.
Easily cause the deformation of organ due to breathing or movement, frisket image can be caused by image registration and filled
Characteristic point on picture reaches matching.
Photographic Subtraction module 1203, carries out subtracting shadow treatment for the frisket image and the full picture to obtaining, and obtains
To the first angiographic image, and, the frisket image and the full picture after the registration are carried out to subtract shadow treatment, obtain
Second angiographic image.
The step has been done and subtract twice shadow treatment, and one is that the frisket image and full picture before registration subtract at shadow
Reason, two is that the full picture before the frisket image after registration and registration is carried out to subtract shadow treatment.
Computing module 1204, for by the gray value summation of the pixel in preset regions in first angiographic image
Gray value summation with the pixel in the preset regions in second angiographic image is compared, and according to comparative result
First angiographic image or second angiographic image are chosen as target image.
Because image border part can interfere to result of calculation, thus choose first angiographic image and this second
The rule of the preset regions is identical in angiographic image, i.e., the preset regions are to remove the region of image border part.
Further, computing module 1204 includes：Calculating sub module 12041 and selection submodule 12042；
Calculating sub module 12041, the ash for calculating the pixel in first angiographic image in the preset regions
Angle value summation, and, the gray value summation of the pixel in second angiographic image in the preset regions；
Choose submodule 12042, for the gray value summation in the corresponding pixel of the first angiographic image with should
The small angiographic image of numerical value is chosen in the gray value summation of the corresponding pixel of the second angiographic image as the target
Image.
If the gray value summation of the corresponding pixel of the first angiographic image is less than second angiographic image pair
The gray value summation of the pixel answered, then choose submodule 12042 and choose the first angiographic image as the target image；Instead
It, then choose submodule 12042 and choose the second angiographic image as the target image.The numerical value of the gray value summation is got over
It is small, represent through over-subtraction shadow, eliminate the image in most non-vascular region.
Image processing module 1205, for according to preset edge detection operator to the edge in the target image medium vessels region
Detected, and enhancing treatment is carried out to the edge of the angiosomes after detection according to preset weighted value, obtained edge enhancing
The target image afterwards.
The preset edge detection operator is not limited, and the preset edge detection operator can be Sobel operators, it is also possible to
Laplacian operators, can also be the edge detection operators such as Canny operators.The selection of weighted value can be according to real image
Enhancing effect chosen, the preset weighted value is bigger, and the enhancing effect to image border is stronger, and the preset weighted value is got over
Small, the enhancing effect to image border is weaker.
Image processing module 1205, being additionally operable to the target image enhanced to edge carries out noise remove and image increasing
By force, the result images are obtained, and the result images is shown.
Further, image processing module 1205, are additionally operable to perform following steps：
The image processing module 1205, is additionally operable to strengthen the edge according to integration noise reduction, frequency domain noise reduction or Recursive Filtering
The target image afterwards carries out noise remove, obtains image to be reinforced；
The image processing module 1205, is additionally operable to carry out grey scale curve stretching to the image to be reinforced；
The image processing module 1205, is additionally operable to image to be reinforced after stretching as the result images.
The process of grey scale curve stretching is because the maximum region of the slope of curve just falls the ash in the image to be reinforced
The region that angle value is concentrated, being stretched by grey scale curve can be stretched to bigger scope by the region that the gray value is concentrated.
In actual applications, the mode of the removal noise is not only comprising integration noise reduction, frequency domain noise reduction or Recursive Filtering, ability
Field technique personnel are according to the mode for removing noise for integrating noise reduction, frequency domain noise reduction or Recursive Filtering described by the embodiment of the present invention
Associate the mode of other removal noises.
The present embodiment details not to the greatest extent, refers to the description of earlier figures 1 and embodiment illustrated in fig. 2, and here is omitted.
In the embodiment of the present invention, acquisition module 1201 obtains frisket image and full picture, image registration module 1202
With the full picture as reference picture, image registration is carried out to the frisket image and the full picture, after obtaining registration
The frisket image, the frisket image and the full picture of the 1203 pairs of acquisitions of Photographic Subtraction module carry out subtracting shadow treatment, obtain
First angiographic image, and, to the frisket image and the full picture after the registration subtract shadow treatment, obtain the
Two angiographic images, computing module 1204 is total by the gray value of pixel in first angiographic image in preset regions
And with the preset regions in second angiographic image in the gray value summation of pixel be compared, and according to comparing knot
Fruit chooses first angiographic image or second angiographic image as target image, and image processing module 1205 pairs should
Target image carries out image definition enhancing treatment, obtains result images, and the result images are shown, compared to existing
Technology, by after image registration and definition enhancing treatment, improve the identification of image, improves the image that subtracts movie queen
Effect, and then diagnosis of the effective auxiliary doctor to patient, are beneficial to the health of patient.
In multiple embodiments provided herein, it should be understood that disclosed system, apparatus and method, can be with
Realize by another way.For example, device embodiment described above is only schematical, for example, the module
Divide, only a kind of division of logic function there can be other dividing mode when actually realizing, such as multiple module or components
Can combine or be desirably integrated into another system, or some features can be ignored, or do not perform.It is another, it is shown or
The coupling each other for discussing or direct-coupling or communication linkage can be the indirect couplings of device or module by some interfaces
Close or communication linkage, can be electrical, mechanical or other forms.
The module that is illustrated as separating component can be or may not be it is physically separate, it is aobvious as module
The part for showing can be or may not be physical module, you can with positioned at a place, or can also be distributed to multiple
On mixed-media network modules mixed-media.Some or all of module therein can be according to the actual needs selected to realize the mesh of this embodiment scheme
's.
In addition, during each functional module in each embodiment of the invention can be integrated in a processing module, it is also possible to
It is that modules are individually physically present, it is also possible to which two or more modules are integrated in a module.Above-mentioned integrated mould
Block can both be realized in the form of hardware, it would however also be possible to employ the form of software function module is realized.
If the integrated module is to realize in the form of software function module and as independent production marketing or use
When, can store in a computer read/write memory medium.Based on such understanding, technical scheme is substantially
The part for being contributed to prior art in other words or all or part of the technical scheme can be in the form of software products
Embody, the computer software product is stored in a storage medium, including some instructions are used to so that a computer
Equipment (can be personal computer, server, or network equipment etc.) performs the complete of each embodiment methods described of the invention
Portion or part steps.And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only
Memory), random access memory (RAM, Random Access Memory), magnetic disc or CD etc. are various can store journey
The medium of sequence code.
It should be noted that for foregoing each method embodiment, in order to simplicity is described, therefore it is all expressed as a series of
Combination of actions, but those skilled in the art should know, the present invention not by described by sequence of movement limited because
According to the present invention, some steps can sequentially or simultaneously be carried out using other.Secondly, those skilled in the art should also know
Know, embodiment described in this description belongs to preferred embodiment, and involved action and module might not all be this hairs
Necessary to bright.
In the above-described embodiments, the description to each embodiment all emphasizes particularly on different fields, and does not have the portion described in detail in certain embodiment
Point, may refer to the associated description of other embodiments.
It is more than the description of the processing method and processing device to angiographic image provided by the present invention, for this area
Technical staff, according to the thought of the embodiment of the present invention, will change in specific embodiments and applications, comprehensive
On, this specification content should not be construed as limiting the invention.
Claims (8)
1. a kind of processing method of angiographic image, it is characterised in that including：
Frisket image and full picture are obtained, and with the full picture as reference picture, the frisket image is carried out
Image registration, obtains the frisket image after registration；
The frisket image for obtaining and the full picture are carried out to subtract shadow treatment, the first angiographic image is obtained, with
And, the frisket image and the full picture after the registration are carried out subtracting shadow treatment, obtain the second angiogram
Picture；
By the gray value summation of the pixel in preset regions in first angiographic image and second angiogram
The gray value summation of the pixel described in image in preset regions is compared, and chooses first blood according to comparative result
Pipe contrastographic picture or second angiographic image are used as target image；
Image definition enhancing treatment is carried out to the target image, result images is obtained, and the result images are shown
Show.
2. method according to claim 1, it is characterised in that described by preset regions in first angiographic image
The gray value of the pixel in the gray value summation of interior pixel and preset regions described in second angiographic image
Summation is compared, and chooses first angiographic image or the second angiographic image conduct according to comparative result
Target image includes：
The gray value summation of the pixel described in first angiographic image in preset regions is calculated, and, described
The gray value summation of the pixel described in two angiographic images in preset regions；
It is corresponding with second angiographic image in the gray value summation of the corresponding pixel of first angiographic image
Pixel gray value summation in choose the small angiographic image of numerical value as the target image.
3. method according to claim 1 and 2, it is characterised in that described that image definition is carried out to the target image
Enhancing is processed, and obtaining result images includes：
The edge in the target image medium vessels region is detected according to preset edge detection operator, and according to preset weight
Value to detection after the edge of the angiosomes carry out enhancing treatment, obtain the enhanced target image in edge；
The target image enhanced to edge carries out noise remove and image enhaucament, obtains the result images.
4. method according to claim 3, it is characterised in that the target image enhanced to edge is made an uproar
Sound is removed and image enhaucament, and obtaining the result images includes：
Noise is carried out according to integration noise reduction, frequency domain noise reduction or the Recursive Filtering target image enhanced to the edge
Remove, obtain image to be reinforced；
Grey scale curve stretching is carried out to the image to be reinforced；
Image described to be reinforced after using stretching is used as the result images.
5. a kind of processing unit of angiographic image, it is characterised in that described device includes：
Acquisition module, for obtaining frisket image and full picture；
Image registration module, for the full picture as reference picture, image registration being carried out to the frisket image, obtains
The frisket image after to registration；
Photographic Subtraction module, shadow treatment subtract for the frisket image and the full picture to obtaining, and obtains the
One angiographic image, and, the frisket image and the full picture after the registration are carried out subtracting shadow treatment, obtain
To the second angiographic image；
Computing module, for by the gray value summation of the pixel in preset regions in first angiographic image with it is described
The gray value summation of the pixel described in the second angiographic image in preset regions is compared, and according to comparative result choosing
First angiographic image or second angiographic image are taken as target image；
Image processing module, for carrying out image definition enhancing treatment to the target image, obtains result images, and to institute
Result images are stated to be shown.
6. device according to claim 5, it is characterised in that the computing module includes：
Calculating sub module, the gray value for calculating the pixel described in first angiographic image in preset regions is total
With, and, the gray value summation of the pixel described in second angiographic image in preset regions；
Submodule is chosen, for the gray value summation in the corresponding pixel of first angiographic image and second blood
The small angiographic image of numerical value is chosen in the gray value summation of the corresponding pixel of pipe contrastographic picture as the target image.
7. the device according to claim 5 or 6, it is characterised in that
Described image processing module, is additionally operable to the edge to the target image medium vessels region according to preset edge detection operator
Detected, and according to preset weighted value to detection after the edge of the angiosomes carry out enhancing treatment, obtain edge increasing
The target image after strong；
Described image processing module, being additionally operable to the target image enhanced to edge carries out noise remove and image enhaucament,
Obtain the result images.
8. device according to claim 7, it is characterised in that
Described image processing module, is additionally operable to enhanced to the edge according to integration noise reduction, frequency domain noise reduction or Recursive Filtering
The target image carries out noise remove, obtains image to be reinforced；
Described image processing module, is additionally operable to carry out grey scale curve stretching to the image to be reinforced；
Described image processing module, image described to be reinforced after being additionally operable to stretching is used as the result images.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710041232.5A CN106803241A (en) | 2017-01-20 | 2017-01-20 | The processing method and processing device of angiographic image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710041232.5A CN106803241A (en) | 2017-01-20 | 2017-01-20 | The processing method and processing device of angiographic image |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106803241A true CN106803241A (en) | 2017-06-06 |
Family
ID=58986941
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710041232.5A Pending CN106803241A (en) | 2017-01-20 | 2017-01-20 | The processing method and processing device of angiographic image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106803241A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108198184A (en) * | 2018-01-09 | 2018-06-22 | 北京理工大学 | The method and system of contrastographic picture medium vessels segmentation |
| CN108814633A (en) * | 2018-03-05 | 2018-11-16 | 深圳市安健科技股份有限公司 | Digital subtraction angiography method, computer storage medium |
| CN111915538A (en) * | 2020-08-19 | 2020-11-10 | 南京普爱医疗设备股份有限公司 | Image enhancement method and system for digital blood vessel subtraction |
| CN114533096A (en) * | 2022-02-21 | 2022-05-27 | 郑州市中心医院 | Artifact removing method and artifact removing system in cerebrovascular angiography |
| WO2023030344A1 (en) * | 2021-08-30 | 2023-03-09 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical image processing |
| CN116342603A (en) * | 2023-05-30 | 2023-06-27 | 杭州脉流科技有限公司 | Method for obtaining arterial input function |
| CN117274502A (en) * | 2023-11-17 | 2023-12-22 | 北京唯迈医疗设备有限公司 | Image processing method and device for assisting interventional operation |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101051384A (en) * | 2007-05-14 | 2007-10-10 | 华中科技大学 | Self adaptive reinforced digital image reducing blood vessle image |
| CN101127117A (en) * | 2007-09-11 | 2008-02-20 | 华中科技大学 | Method for segmenting blood vessel data using serial DSA image |
| US20120163686A1 (en) * | 2010-09-23 | 2012-06-28 | Siemens Corporation | System and method for 2-d/3-d registration between 3-d volume and 2-d angiography |
| CN103544690A (en) * | 2012-07-10 | 2014-01-29 | 伊姆普斯封闭式股份有限公司 | Method for acquisition of subtraction angiograms |
| CN104318567A (en) * | 2014-10-24 | 2015-01-28 | 东北大学 | Method for segmenting kidney blood vessel atrioventricular on basis of medical images |
| CN104700397A (en) * | 2013-12-06 | 2015-06-10 | 株式会社东芝 | Medical image processing device used for segmentation of medical image structures |
| CN105005998A (en) * | 2015-08-05 | 2015-10-28 | 大连理工大学 | Cerebrovascular image segmentation method based on multi-angle serialized image space feature point set |
-
2017
- 2017-01-20 CN CN201710041232.5A patent/CN106803241A/en active Pending
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101051384A (en) * | 2007-05-14 | 2007-10-10 | 华中科技大学 | Self adaptive reinforced digital image reducing blood vessle image |
| CN101127117A (en) * | 2007-09-11 | 2008-02-20 | 华中科技大学 | Method for segmenting blood vessel data using serial DSA image |
| US20120163686A1 (en) * | 2010-09-23 | 2012-06-28 | Siemens Corporation | System and method for 2-d/3-d registration between 3-d volume and 2-d angiography |
| CN103544690A (en) * | 2012-07-10 | 2014-01-29 | 伊姆普斯封闭式股份有限公司 | Method for acquisition of subtraction angiograms |
| CN104700397A (en) * | 2013-12-06 | 2015-06-10 | 株式会社东芝 | Medical image processing device used for segmentation of medical image structures |
| CN104318567A (en) * | 2014-10-24 | 2015-01-28 | 东北大学 | Method for segmenting kidney blood vessel atrioventricular on basis of medical images |
| CN105005998A (en) * | 2015-08-05 | 2015-10-28 | 大连理工大学 | Cerebrovascular image segmentation method based on multi-angle serialized image space feature point set |
Non-Patent Citations (2)
| Title |
|---|
| 周正东 等: "脑部数字减影图象增强研究", 《中国图象图形学报》 * |
| 黄惠芬 等: "《数字图像司法取证技术》", 30 April 2015, 山东大学出版社 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108198184A (en) * | 2018-01-09 | 2018-06-22 | 北京理工大学 | The method and system of contrastographic picture medium vessels segmentation |
| CN108198184B (en) * | 2018-01-09 | 2020-05-05 | 北京理工大学 | Method and system for vessel segmentation in contrast images |
| CN108814633A (en) * | 2018-03-05 | 2018-11-16 | 深圳市安健科技股份有限公司 | Digital subtraction angiography method, computer storage medium |
| CN111915538A (en) * | 2020-08-19 | 2020-11-10 | 南京普爱医疗设备股份有限公司 | Image enhancement method and system for digital blood vessel subtraction |
| CN111915538B (en) * | 2020-08-19 | 2024-03-19 | 南京普爱医疗设备股份有限公司 | Image enhancement method and system for digital blood vessel subtraction |
| WO2023030344A1 (en) * | 2021-08-30 | 2023-03-09 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical image processing |
| CN114533096A (en) * | 2022-02-21 | 2022-05-27 | 郑州市中心医院 | Artifact removing method and artifact removing system in cerebrovascular angiography |
| CN116342603A (en) * | 2023-05-30 | 2023-06-27 | 杭州脉流科技有限公司 | Method for obtaining arterial input function |
| CN116342603B (en) * | 2023-05-30 | 2023-08-29 | 杭州脉流科技有限公司 | Method for obtaining arterial input function |
| CN117274502A (en) * | 2023-11-17 | 2023-12-22 | 北京唯迈医疗设备有限公司 | Image processing method and device for assisting interventional operation |
| CN117274502B (en) * | 2023-11-17 | 2024-03-01 | 北京唯迈医疗设备有限公司 | Image processing method and device for assisting interventional operation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106803241A (en) | The processing method and processing device of angiographic image | |
| Hatvani et al. | Deep learning-based super-resolution applied to dental computed tomography | |
| CN111540025B (en) | Predicting images for image processing | |
| Liu et al. | Motion artifacts reduction in brain MRI by means of a deep residual network with densely connected multi-resolution blocks (DRN-DCMB) | |
| CN107533755B (en) | Apparatus and method for improving medical image quality | |
| Xian et al. | Main coronary vessel segmentation using deep learning in smart medical | |
| US9275437B2 (en) | Method for efficient digital subtraction angiography | |
| Unberath et al. | Consistency‐based respiratory motion estimation in rotational angiography | |
| Fazlali et al. | Vessel region detection in coronary X-ray angiograms | |
| CN107680057A (en) | The method and device of ultrasonoscopy enhancing | |
| CN106530236A (en) | Medical image processing method and system | |
| Sadda et al. | Real-time medical video denoising with deep learning: application to angiography | |
| CN104463815B (en) | DSA image generating methods and system | |
| US12039718B2 (en) | Inference apparatus, medical apparatus, and program | |
| Liu et al. | SDCNet: Smoothed dense-convolution network for restoring low-dose cerebral CT perfusion | |
| CN112991220A (en) | Method for correcting image artifacts by convolutional neural network based on multiple constraints | |
| Wagner et al. | Noise reduction for curve‐linear structures in real time fluoroscopy applications using directional binary masks | |
| EP4070236A1 (en) | Medical image synthesis for motion correction using generative adversarial networks | |
| EP3658031B1 (en) | Motion compensated cardiac valve reconstruction | |
| Căleanu et al. | A study over the importance of arterial phase temporal parameters in focal liver lesions CEUS based diagnosis | |
| CN113538419B (en) | Image processing method and system | |
| Ahmad et al. | CPS-based fully automatic cardiac left ventricle and left atrium segmentation in 3D MRI | |
| Taleb et al. | A 3D space–time motion evaluation for image registration in digital subtraction angiography | |
| Solak et al. | Adrenal Tumor Segmentation on U-Net: A Study About Effect of Different Parameters in Deep Learning. | |
| CN113989171A (en) | Subtraction map generation method and device, storage medium and computer equipment |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170606 |