CN102592274B - Processing procedure for the radiation image of stenosis detection - Google Patents
Processing procedure for the radiation image of stenosis detection Download PDFInfo
- Publication number
- CN102592274B CN102592274B CN201110387801.4A CN201110387801A CN102592274B CN 102592274 B CN102592274 B CN 102592274B CN 201110387801 A CN201110387801 A CN 201110387801A CN 102592274 B CN102592274 B CN 102592274B
- Authority
- CN
- China
- Prior art keywords
- image
- tremulous pulse
- interested
- line
- interest
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The present invention is provided to the processing procedure of the radiation image of stenosis detection.The present invention relates to the processing procedure of the radiation image of the region of interest of patient, these radiation images are 2D projection pictures, the most elongated utensil has been previously inserted in tremulous pulse, and this process comprises the following steps: to obtain (S1, S2) by the first image (I of this region of interest1) and the second image (I2) the least one set image that constitutes, this first image (I1) it is the image of this region of interest being wherein previously implanted into reference substance, this second image is not have reference substance or have the image of this region of interest of minimum number reference substance, and often group is corresponding to given angle；This process comprises the following steps: to split (S3) this first image (I for often group1) to detect multiple tremulous pulsies of this region of interest；Segmentation (S4) this second image (I2) to detect and to isolate this utensil；(S5) a plurality of line (C is limited in the first segmentation imagei), every line (Ci) limit tremulous pulse；From the second segmentation image and the line of restriction determine (S6) corresponding to wherein already inserted into the tremulous pulse interested of tremulous pulse of this utensil, the line (C7) of this tremulous pulse interested is closest to this utensil；(S7) is applied to detect the damage of this tremulous pulse interested in this tremulous pulse interested the quantitative analysis algorithm of coronary artery injury.
Description
Technical field
The present invention relates to medical imaging field and relate more specifically to be applied to blood vessel radiation Jie
Enter the radiation field of imaging field.
Background technology
Blood vessel intervention radiation is included in imaging and controls the lower process carried out and specifically treat the heart
Myocardial ischemia.
Myocardial ischemia be a kind of developed country's impact more than the disease of 1/3rd populations and
It embodies via narrow (that is, arteriarctia).
There is some treatments for treating this disease, it is included in coronary artery and is inserted into propping up
Frame.
To this, cardiologist utilizes intervention imaging, and it provides coronarius possible to him
Damage and quantify and characterize and the narrowest characterizing to select to be introduced propping up
The size of frame.
Cardiologist generally uses the area image including tremulous pulse to be treated, thereon his edge
The some labellings of tremulous pulse manual placement to be treated to carry out narrow analysis, i.e. determine narrow position
In place and determine the size for support necessary to narrow treatment.
This detection and quantization need radiologist and medical imaging apparatus mutual and detection is wanted
The previous step of the tremulous pulse for the treatment of.
Here inferior position detection and quantify time-consuming and be to be difficult to carry out for doctor.
Summary of the invention
The present invention eliminates these inferior positions.
According to first aspect, the present invention relates to the radiation image of region of interest for processing patient
Process, these radiation images be 2D projection in image, the most elongated utensil is the most prior
Inserting in tremulous pulse, this process comprises the following steps:
Obtain the first image by region of interest and the least one set image of the second image construction, should
First image is the image of the region of interest being wherein previously implanted into reference substance, and this second image is not
Having reference substance or have the image of region of interest of few quantity reference substance, often group is corresponding to given
Angle；This process also comprises the following steps: for often organizing
Split this first image to detect multiple tremulous pulsies of region of interest；Split this second image with
Detect and isolate this utensil；Limiting a plurality of line in the first segmentation image, every line limits dynamic
Arteries and veins；Determine corresponding to wherein moving already inserted into utensil from the second segmentation image and the line limited
The tremulous pulse interested of arteries and veins, the line of this tremulous pulse interested is closest to this utensil；By coronary artery injury
Quantitative analysis algorithm be applied to this tremulous pulse interested to detect the damage of this tremulous pulse interested.
Owing to detecting tremulous pulse interested, it is possible to quantification of stenosis is without along tremulous pulse artificial limit
Put and thus the most mutual with medical image system.
This provides quick and easy in performing blood vessel injury detection.
Other aspects of this process are as follows:
-by using the point of the line placement along tremulous pulse interested to carry out answering of quantitative analysis algorithm
With；
-determine that the step of tremulous pulse interested includes: utilizing distance criterion, it combines the Europe of tremulous pulse line
Must be apart from (Euclidian distance) and the device all along the curve limiting elongated utensil in several
The assessment of the orientation difference of tool and this line；
The line of-restriction tremulous pulse is the centrage of the axis of symmetry limiting tremulous pulse respectively；
-obtaining at least two group images, this process includes: complete the tremulous pulse interested for often organizing
Really timing, carries out 3D reconstruction step to obtain the 3D rendering of tremulous pulse interested；
-on 3D tremulous pulse interested, use quantitative analysis algorithm；
-obtain a series of groups of image, wherein at the tremulous pulse each interested for often organizing detection
On complete the identification of tremulous pulse interested to obtain the two-dimensional movement field of the line of tremulous pulse interested；
-obtain the groups of image of multiple series, wherein use quantitative analysis to calculate on image often organizing
Method obtains the moving field of projection in the image of the line of tremulous pulse interested；
The radiation image of-synthesis is previously acquired and is stored in the memorizer of medical image system
Image in unit.
According to second aspect, the present invention relates to medical image system, it includes for performing basis
The parts of the process of a first aspect of the present invention.
According to the third aspect, the present invention relates to computer program, it includes for performing according to this
The machine instruction of the process of the first aspect of invention.
Accompanying drawing explanation
Other characteristics of the present invention and advantage will from purely illustrative and nonrestrictive and
The following explanation that must be considered in conjunction with the accompanying appears, wherein:
Fig. 1 diagram is according to the medical image system 5 of the present invention；
Fig. 2 schematically illustrates the step of the processing procedure of the radiation image according to the present invention；
The first figure that Fig. 3 diagram is obtained by the processing procedure of the radiation image according to the present invention
Picture；
The second figure that Fig. 4 diagram is obtained by the processing procedure of the radiation image according to the present invention
Picture；
Fig. 5 is schematically illustrated according to acquisition during the processing procedure of the radiation image of the present invention
Tremulous pulse interested.
Detailed description of the invention
Medical image system
Fig. 1 schematically illustrates the medical image system 100 for gathering radiation image.
This medical image system 100 includes the support 1 for housing patient 10 to be checked, sets
Count into the source 2 of emitting x-ray 3, be placed on this source 2 offside and be arranged to check for by this
The detector 4 of X-ray, control unit 6, memory element 7 and the display unit that source 2 is launched
8。
Source 2 and the detector 4 of X-ray are attached by C-arm 5.Such arm 5 more generally claims
For arch.This arm 5 can orient according to three degree of freedom.
Detector 4 can be semiconductor image sensor, and it includes the transistor of such as non-crystalline silicon
Cesium iodide phosphorus (scintillator) on/photodiode array.Other detectors being suitable for are: CCD
Sensor, a kind of Direct Digital detector that directly X-ray is converted into digital signal.At figure
In 1, the detector 4 of diagram is plane and limits plane picture surface, and other geometric forms
Shape can be suitable.
Control unit 6 is connected to encircle 5 by wired or wireless connection.Control unit 6 is passed through
Some parameters such as the angle location of the fixing radiation dose such as launched by x-ray source and arm 5
Control to gather.The position of control unit 6 control arm 5, i.e. source 2 is relative to detector 4
Position.
Control unit 6 can include reading device (not shown), such as disk player or
CD-ROM reader, DVD-ROM reader or connectivity port, connectivity port is used for reading
The such as instruction media (not shown) such as disk, CD-ROM, DVD-ROM or usb key
Processing procedure instruction or more generally by any removable storage medium or even via
Network connects the instruction of the processing procedure reading instruction media (not shown).
Memory element 7 is connected to control unit 6 for recording parameter and the image of collection.Can
Memory element 7 is made to be positioned at inside control unit 6 or outside it.
Memory element 7 can be by hard-drive or SSD or any other removable and rewritable storage
Parts (usb key, memory card etc.) are formed.Memory element 7 can be control unit 6
ROM/RAM memorizer, usb key, memory card, the memorizer of central server.
Display unit 8 is connected to control unit 6 for showing the image of collection and/or about collection
Control the information of parameter.
Display unit 8 can be such as computer screen, monitor, flat screen, plasma
The display device of any other type of screen or known type.
Such display unit 8 allows the collection that doctor controls radiation image.
Medical image system 100 is coupled in processing system 200.This processing system 200 includes meter
Calculate unit 9 and memory element 10.
Processing system 200 receives the memory element gathering and being stored in medical image system 100
Image in 4, processing system 200 uses this image to implement the process (seeing below) of some quantity.
The data of the memory element 4 of medical image system 100 are sent to the meter of processing system 200
Calculate unit 9 can via interiorly or exteriorly information network or with such as disk, CD-ROM,
Any applicable physical store such as DVD-ROM, outside hard-drive, usb key, SD card
Device medium is carried out.
Computing unit 9 be such as computer or multiple computer, processor or multiple processor,
Microcontroller or multiple microcontroller, microcomputer or multiple microcomputer, able to programme from
Motivation or multiple automat, special IC or multiple such circuit, other able to programme
Circuit or include other devices of the computers such as such as work station.
As version, computer 9 can include reading device (not shown), such as magnetic
Dish reader, CD-ROM or DVD-ROM reader or connectivity port, connectivity port is used
In reading such as, the instruction media such as disk, CD-ROM, DVD-ROM or usb key is not (
Illustrate) processing procedure instruction or more generally by any removable storage medium or even
The instruction of the processing procedure reading instruction media (not shown) is connected via network.
And, processing system includes the storage list of the data generated for storage by computing unit 9
Unit 11.
Computing unit 9 may be coupled to display unit 8 (as in FIG) or additionally connects
To another display unit (not shown).
The processing procedure of radiation image
The embodiment of the processing procedure of the radiation image according to the present invention is described below.Fig. 2 shows
Meaning ground illustrates the step of this process.
Oneth I1With 5 the 2nd I2Acquisition S1, S2 of image.
The processing procedure of radiation image utilizes two radiation image I of the region of interest of patient1、I2,
Utensil introduces in this patient the most in advance, is elongated utensil below.
The utensil of the type is such as seal wire, conduit or its combination, and preferably can
It is easily introduced into the seal wire of support.
First image I1Corresponding to wherein having been injected into the radiation image of the region of interest of reference substance.
Such reference substance is such as iodine, and region of interest is typically intended to patient crown for the treatment of
Tremulous pulse district.
Second image I2Corresponding to without reference substance or to have the sense of few quantity reference substance emerging
The radiation image in interest district.
Few quantity reference substance mean to show utensil but this utensil not by tremulous pulse cover right
Quantity according to product.
First and second image I1、I2Can be by adopting in the whole processing procedure of radiation image
Acquisition component and obtain, or can obtain from the memory element of medical image system.
It is noted that the first and second image I1、I2Come from by above-described medical imaging system
For given angle, (that is, relative to being orthogonal to supporting station, (patient is placed on it to x-ray source to system
On) orientation in direction) collection made.
First and second image I1、I2It it is 2D projection picture.
In remaining describes, one group of image definition is to be by first gathered for given angle
Pairing with the second image construction.
Further contemplate a series of image in groups, i.e. for given angle, gather in the different time
Some groups of image I1、I2。
First image I1Segmentation S3
This step includes: from the first image I1And the reference substance owing to injecting extracts and isolation is suffered from
The tremulous pulse of person.
More accurately, extract and isolate the tremulous pulse network of region of interest.
Fig. 3 schematically illustrates the first segmentation image I1' example, it includes patient dynamic of display
Network of vein network 101, include the tremulous pulse 103 of the constriction 102 corresponding to damage to be detected.
Second image I2Segmentation S4
This step S4 includes: extract from the region of interest not having reference substance of display and isolation is inserted
Enter the utensil of the tremulous pulse interested of patient.
Fig. 4 schematically illustrates the second segmentation image I2' example, it include and first segmentation image
I1' the identical network 101 by dotted line diagram, and utensil 201, seal wire inserts tremulous pulse here
In the tremulous pulse of network 101.It practice, when completing this segmentation step S4, at image I2In '
Only show utensil, and network 101 represents here by instruction.
It is noted that segmentation image is those skilled in that art's widely-known technique and will not
Carry out any described in further detail.
The restriction S5 of the line of tremulous pulse
This step S5 includes: limit line C for each tremulous pulsei, such as, limit the symmetry of tremulous pulse
The centrage of axis.
Limit centrage CiThe technology described in documents below is such as used to realize:
Karl Krissian, Gregoire Malandain, Nicholas Ayache,
R é gis Vaillant and Yves Trousset: " Model-Based Detection of
Tubular Structures in 3D Images”Computer Vision and Image
Understanding, Vol.80, num.2, p.130-171,2000.
Fig. 5 pictorial images I3', for this image I3' in tremulous pulse network 101, limited 9
Bar centrage Ci(i=1 to 9).
Detection S5 including the tremulous pulse of utensil
This step S5 includes: detect tremulous pulse interested (wherein already inserted into the tremulous pulse of utensil).
The detection of tremulous pulse interested is from the second image I2' and from image I3' implement.
Such detection uses splits image I second2The utensil 201 of isolation and detection in ' and limit
Fixed centrage CiBetween distance criterion implement.
It is that tremulous pulse interested corresponds to the distance for tremulous pulse between utensil and centrage C
Little tremulous pulse.
This can be Euclidean distance and utensil and this centrage of the centrage C combining tremulous pulse
The distance criterion of the assessment of the orientation difference of C.
The possible narrow determination S6 of tremulous pulse interested
Once identify the tremulous pulse interested wherein already inserted into utensil, quantitative analysis can be carried out.
In order to accomplish this point, coronary artery injury (" quantitative coronary analysis " (QCA))
The process of quantitative analysis be applied to tremulous pulse interested.
Such analysis process describes the most in the following documents: Paul L. Van Herck et al.
" Quantitative Coronary Arteriography on Digital Flat-Panel System "
Catheterization and Cardiovascular Interventions, p.192-200,2004.
QCA analyzes based on being efficiently used for clinically detecting the algorithm of profile and can determine
Inaccessible percentage ratio, natural arterial and represent the diameter of narrow tremulous pulse and the obturation of size.
Specifically, such Algorithm Analysis is all along the tremulous pulse of centrage C and the most true
The apparent diameter of the tremulous pulse centered at each point of line.The change of this diameter and specifically reducing
It it is the indicant of the pathology existence such as such as arterial injury.Quantify to include by relative to normal diameter
Area percentage determines this reduction.
Obtain (S8, S8 ') 3D rendering
This process above-described can apply to the some 2D projection taked for different angles
Image.
Use known type and can apply to line 3D reconstruct can produce 3D line, its describe
The geometry of tremulous pulse.More specifically, be obtained in that the 3D view of tremulous pulse interested.
This can apply to a series of images to obtain a series of 3D rendering and therefore obtain under
Row key element:
-analyze the projection moved by the tremulous pulse developing description of the line identified in this image series；
The projection in the continuous print image of-quantitative analysis tremulous pulse；
-analyze tremulous pulse in three-dimensional from the process of the some image series according to different angles collection
Movement in space；
Three peacekeeping time quantitative analyses of-tremulous pulse interested.
Computer program
The processing procedure of radiation image can use computer program, and (it includes for performing this mistake
The machine instruction of journey) form be advantageously carried out.
List of parts
| Fig. 1 | |||
| 1 | Supporting station | 2 | Source |
| 3 | X-ray beam | 4 | Detector |
| 5 | Arm | 6 | Control unit |
| 7 | Memory element | 8 | Display unit |
| 9 | Computing unit | 10 | Memory element |
| 11 | Memory element | 100 | Medical image system |
| 200 | Processing system | ||
| Fig. 2 | |||
| S1 | Obtain the step of image | S2 | Obtain the step of image |
| S3 | The step of segmentation image | S4 | The step of segmentation image |
| S5 | Detecting step | S6 | Determine step |
| S7 | Applying step | ||
| Fig. 3 | |||
| 101 | Network | 102 | Constriction |
| 103 | Tremulous pulse | ||
| Fig. 4 | |||
| 201 | Utensil | 101 | Network |
| Fig. 5 | |||
| C1-C9 | Line | 101 | Network |
| 102 | Constriction |
Claims (9)
1. a medical imaging apparatus, for carrying out the radiation image of the region of interest of patient
Processing, described radiation image is 2D projection picture, and the most elongated utensil has been previously inserted
In tremulous pulse, medical imaging apparatus includes:
-be used for obtaining (S1, S2) by the first image (I of described region of interest1) and second
Image (I2) device of least one set image that constitutes, described first image (I1) it is wherein
Having been injected into the image of the described region of interest of reference substance in advance, described second image is not have
Reference substance or there is the image of described region of interest of few quantity reference substance, often group corresponding to
Given angle；
-be used for splitting (S3) described first image (I1) many to detect described region of interest
The device of individual tremulous pulse；
-be used for splitting (S4) described second image (I2) to detect and to isolate described utensil
Device；
-for limiting (S5) a plurality of line (C in the first segmentation imagei) device, every
Line (Ci) limit tremulous pulse；
-for determining that (S6) is corresponding to the most from the line of the second segmentation image and described restriction
Through inserting the device of the tremulous pulse interested of the tremulous pulse of described utensil, the line of described tremulous pulse interested
(C7) closest to described utensil；
-for applying (S7) in described interested the quantitative analysis algorithm of coronary artery injury
Tremulous pulse detects the device of the damage of described tremulous pulse interested.
2. device as claimed in claim 1, also includes: for described quantitative analysis being calculated
The device that the application of method is carried out by using the point of the line placement along described tremulous pulse interested.
3. device as claimed in claim 1, also includes: for performing the dress of distance criterion
Putting, the described device for performing distance criterion is configured to comment along the curve of described elongated utensil
Estimate the orientation difference of described elongated utensil and described line (Ci).
4. device as claimed in claim 1, the line wherein limiting described tremulous pulse is to limit respectively
Determine the centrage of the axis of symmetry of tremulous pulse.
5. device as claimed in claim 1, wherein obtains at least two group image (I1、
I2), described device also includes: for having worked as the tremulous pulse described interested for often organizing really
Regularly, 3D reconstruction step is carried out to obtain the device of the 3D rendering of described tremulous pulse interested.
6. device as claimed in claim 1, also includes: for described quantitative analysis being calculated
The device that method uses on 3D tremulous pulse interested.
7. device as claimed in claim 1, also includes: be used for obtaining a series of groups of
Image, and wherein using described sense emerging on the tremulous pulse each interested of detection for often organizing
The identification of interest tremulous pulse obtains the device of the two-dimensional movement field of the line of described tremulous pulse interested.
8. device as claimed in claim 1, also includes: for obtaining the one-tenth of multiple series
The image of group and wherein use described quantitative analysis algorithm to obtain in institute on image often organizing
The device of the moving field of projection in the image of the line stating tremulous pulse interested.
9. device as claimed in claim 1, the described radiation image wherein obtained is previous
Gather and be stored in the image in the memory cell of described medical imaging apparatus.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP10306249.3A EP2453408B1 (en) | 2010-11-12 | 2010-11-12 | Method for processing radiographic images for stenosis detection |
| EP10306249.3 | 2010-11-12 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102592274A CN102592274A (en) | 2012-07-18 |
| CN102592274B true CN102592274B (en) | 2016-12-14 |
Family
ID=
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101422352A (en) * | 2008-12-10 | 2009-05-06 | 华北电力大学（保定） | Interactive coronary artery virtual angioscope implementation method |
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101422352A (en) * | 2008-12-10 | 2009-05-06 | 华北电力大学（保定） | Interactive coronary artery virtual angioscope implementation method |
Non-Patent Citations (2)
| Title |
|---|
| 《A comparison of line enhancement techniques-Applications to guide-wire detection and respiratory motion tracking》;Vincent Bismuth et al;《Proceddings of SPIE》;20090208;第7259卷;72591M-1-72591M-9 * |
| 《Robust guidewire segmentation through boosting,clustering and linear programming》;Nicolas Honnorat et al;《Biomedical imaging:from nano to macro,2010 IEEE international symposium on》;20100414;924-927 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8880148B2 (en) | Treatment process of radiological images for detection of stenosis | |
| CN105120738B (en) | Narrow treatment is planned | |
| CN109035355B (en) | System and method for PET image reconstruction | |
| CN104008540B (en) | Method and apparatus for executing the registration of medical image | |
| CN103889328B (en) | Perfusion imaging | |
| CN102622743B (en) | For comparing the method and apparatus of 3D and 2D view data | |
| CN107787203A (en) | Image registration | |
| US8831323B2 (en) | Method and apparatus for measuring activity of a tracer | |
| US10910101B2 (en) | Image diagnosis support apparatus, image diagnosis support method, and image diagnosis support program | |
| CN111260703B (en) | Method, system, medium and storage medium for obtaining spinal straightening image set | |
| CN115209808A (en) | Learning model creation method, image generation method, and image processing device | |
| WO2016037969A1 (en) | Medical imaging apparatus | |
| CN109313803B (en) | Method and apparatus for mapping at least part of a structure in an image of at least part of a body of a subject | |
| US11213220B2 (en) | Method for determining in vivo tissue biomarker characteristics using multiparameter MRI matrix creation and big data analytics | |
| CN106663321A (en) | Imaging data statistical testing including a stereotactical normalization with a personalized template image | |
| US10402991B2 (en) | Device and method for registration of two images | |
| JP2013514117A (en) | Associating acquired images with objects | |
| US20200261032A1 (en) | Automatic identification and segmentation of target regions in pet imaging using dynamic protocol and modeling | |
| CN102592274B (en) | Processing procedure for the radiation image of stenosis detection | |
| US20150025359A1 (en) | Method for evaluation and comparison of a chronological sequence of combined medical imaging examinations and also a medical imaging system which is designed for executing the inventive method | |
| CN103295246B (en) | Processing to interventional radiology image is analyzed by ECG | |
| KR20220121217A (en) | Device and method for diagnosing brain hemorrhage based on deep learning | |
| JP5700712B2 (en) | Patrac analysis of 3D dynamic myocardial nuclear medicine image data | |
| JP6253474B2 (en) | Information processing apparatus, computer program, and recording medium | |
| Noz et al. | Can the specificity of MRI breast imaging be improved by fusing 3D MRI volume data sets with FDG PET? |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |