CN108805815A - Blood vessel based on X-ray angiographic image stretches method for reconstructing - Google Patents
Blood vessel based on X-ray angiographic image stretches method for reconstructing Download PDFInfo
- Publication number
- CN108805815A CN108805815A CN201810612090.8A CN201810612090A CN108805815A CN 108805815 A CN108805815 A CN 108805815A CN 201810612090 A CN201810612090 A CN 201810612090A CN 108805815 A CN108805815 A CN 108805815A
- Authority
- CN
- China
- Prior art keywords
- point
- blood vessel
- image
- ray
- reconstructing
- 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
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000003384 imaging method Methods 0.000 description 4
- 230000003902 lesion Effects 0.000 description 4
- 238000002583 angiography Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002980 postoperative effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4023—Scaling of whole images or parts thereof, e.g. expanding or contracting based on decimating pixels or lines of pixels; based on inserting pixels or lines of pixels
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/04—Context-preserving transformations, e.g. by using an importance map
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a kind of blood vessels based on X-ray angiographic image to stretch method for reconstructing, including：Obtain X-ray angiogram data and vessel centerline data；Certain distance N is extended respectively along normal both forward and reverse directions to point each of in vessel centerline, both forward and reverse directions respectively obtain a point、；By four adjacent points、、、Be connected quadrangularly；Obtain image coordinate of each quadrangle in X-ray angiogram data；Corresponding pixel grey scale is obtained according to image coordinate, obtains gray areas；Obtained gray areas is merged along straight line, obtain blood vessel stretches reconstruction image.Blood vessel can be carried out along its center line stretching reconstruction, doctor is contributed to rapidly and accurately to analyze blood vessel qualitatively or quantitatively.
Description
Technical field
The present invention relates to technical field of medical image processing, more particularly to a kind of based on X-ray angiographic image
Blood vessel stretches method for reconstructing.
Background technology
Currently, most domestic hospital all be that X-ray single armed angiography system does x-ray imaging to patient, pass through rotation
Turn radiography arm and obtains a contrastographic picture for corresponding to different radiography angles.The mainstream diagnosis side of current X-ray angiography
Method is all doctor by experience or using in the mating image playout software of X-ray single armed angiography system, by simple
2 distance measurement tools are in art or the postoperative blood vessel in contrastographic picture carries out quantitative or qualitative analysis.
But it has the following defects：
1, the personal experience of doctor and level will appear differentiation, and contrastographic picture is different from caused by other image documentation equipments
It is necessary to have the long-term experiences for carrying out contrastographic picture processing for image；
2, at present in image processing software, entrained survey tool is all more simple, can meet measurement demand, but not
Can meet the needs of of accurately measuring, and there is not in the contrastographic picture blood vessel for measuring complicated bend for simple survey tool
Hold the error ignored；
3, in the analytic process man-machine interactively uncertainty；
4, due to the characteristic of X-ray angiographic image, in specific time, contrast agent flows through blood vessel, enables the clinician to see
Blood vessel dynamic tube chamber situation is examined, but just due to dynamic observation, the observation difficulty being likely to result in time increases, and then loses
Leak parts of lesions situation on blood vessel；
5, it is subjective judgement that all judging results to blood vessel, which are all people, do not follow systematically, the data of uniform rules
It generates.
Invention content
In order to solve the above technical problems, the purpose of the present invention is：It provides a kind of based on X-ray angiographic image
Blood vessel stretches method for reconstructing, can carry out along its center line blood vessel to stretch reconstruction, contribute to doctor rapidly and accurately qualitative or
Quantitatively analyze blood vessel.
The technical scheme is that：
A kind of blood vessel based on X-ray angiographic image stretches method for reconstructing, includes the following steps：
S01：Obtain X-ray angiogram data and vessel centerline data；
S02：Certain distance N, both forward and reverse directions are extended respectively along normal both forward and reverse directions to point each of in vessel centerline
Respectively obtain a point An、Bn；
S03：By four adjacent point An、Bn、An+1、Bn+1Be connected quadrangularly；
S04：Obtain image coordinate of each quadrangle in X-ray angiogram data；
S05：Corresponding pixel grey scale is obtained according to image coordinate, obtains gray areas；
S06：Obtained gray areas is merged along straight line, obtain blood vessel stretches reconstruction image.
In preferred technical solution, in the step S02,Wherein,For method
Line direction vector, PnCentered on point on line.
In preferred technical solution, the step S04 further includes, using image interpolation method in X-ray angiographic image
There is no the point inserted with clear gray value in clear gray value point in data.
In preferred technical solution, described image interpolation method is closest first method, including：In the neighbouring of gray scale point to be asked
In four known gray scale points, the gray value of the distance nearest point of gray scale point to be asked is given to gray scale point to be asked.
In preferred technical solution, described image coordinate PtCoordinate acquired by following formula：
P0For origin, d is contrastographic picture pel spacing, PwFor world coordinates.
Compared with prior art, it is an advantage of the invention that：
Blood vessel can stretch reconstruction by the method for the present invention along its center line, and blood vessel stretches reconstruction can be more intuitive
Blood vessel is presented to doctor, is conducive to doctor's multiple field, the multi-angle observation of realization to blood vessel when carrying out clinical diagnosis.Simultaneously
After being straightened due to blood vessel, hemadostewnosis situation and patch situation are observed what is be more prone to, and doctor can more quickly
Blood vessel is accurately analyzed qualitatively or quantitatively, finds lesion locations, weighs lesion situation.Therefore, it is X-ray that blood vessel, which stretches reconstruction,
The important component of Angiographic findings, to patient's lesion vessels analysis important in inhibiting.
Description of the drawings
The invention will be further described with reference to the accompanying drawings and embodiments：
Fig. 1 is flow chart of the method for the present invention；
Fig. 2 is X-ray angiographic image schematic diagram；
Fig. 3 is the blood vessel straightening effect figure in Fig. 2.
Specific implementation mode
In order to make the objectives, technical solutions and advantages of the present invention clearer, With reference to embodiment and join
According to attached drawing, the present invention is described in more detail.It should be understood that these descriptions are merely illustrative, and it is not intended to limit this hair
Bright range.In addition, in the following description, descriptions of well-known structures and technologies are omitted, to avoid this is unnecessarily obscured
The concept of invention.
As shown in Figure 1, the blood vessel based on X-ray angiographic image of the present invention stretches method for reconstructing, including following step
Suddenly.
Step S1：Obtain image data and vessel centerline data；
The acquisition of image data has three approach：First way is used and is directly connected to x-ray imaging machine, is passed through
Wired mode carries out the transmission of image data；The second way, which is used, to be stored and is transmitted with heart internal medicine in hospital medical image
System is attached, and image data is obtained from the system by wired mode；The third mode is used to be situated between with storing
Matter, in the control machine of x-ray imaging machine or the server of the medical image storage and Transmission system that are connected with x-ray imaging machine
On, carry out the copy of image data.
The acquisition modes of vessel centerline data, the grid side for being mainly based upon unit circle and rolling tracking used herein
Method.Rasterizing is carried out using by the region of required calculating, then calculates the central point of each grid to the most narrow spacing on both sides boundary
From, then set two minimum range absolute value of the difference permissible range, if result in permissible range, can be explained the point
For the central point of this grid region.Repeat above step and obtain two central points A1, A2, be with radius | A1A2 | draw circle
Shape should can circumferentially find two points, and the minimum range on the two point to both sides boundaries is equal, wherein the point not overlapped with A1
A3 is another central point, repeats above step, until finding all central points for meeting length, is connected and can be obtained center
Line.
Step S2：One is respectively obtained respectively along normal both forward and reverse directions extended distance N to point each of in vessel centerline
Point An(n=1,2n), Bn(n=1,2n)；Specifically, it is assumed that line direction vector isSo on center line
Any point PnThe coordinate for the corresponding points point that the online direction (n=0,1 ..., n-1) extends can acquire as follows：
Step S3：Four adjacent point An(n=1,2n), Bn(n=1,2n), An+1(n=1,
2···n)、Bn+1(n=1,2n) is connected quadrangularly；Since this paper is for the display requirement of the center line of calculating
It is smooth as possible, therefore the very little that the range accuracy of adjacent center point will be prescribed, and because obtain the positive and negative side of central point normal
It is fixed N to the distance of extension, therefore can is almost rectangle by the way that four adjacent points are connected the quadrangle to be formed, repeats
Several times, can be considered can obtain several with successional rectangle.
Step S4：Due to the quadrangle of drafting resolution ratio and practical X-ray angiogram data resolution ratio not
Equally, X-ray angiogram data contains not specific enough the point of some gray values, in order to make these point can on quadrangle
Point it is corresponding as possible and possess specific gray value, this purpose is reached using image interpolation method.The effect of image interpolation method
It is the point by X-ray angiogram data, being inserted into more clear gray scales, is such as but not limited to using closest
First method gives the gray scale of the distance nearest point of gray scale point to be asked in the known gray scale point of neighbouring four of gray scale point to be asked
Gray scale point to be asked, the quadrangle resolution ratio that the resolution ratio of all X-ray angiographic images can in this way drawn with previous step
The gray scale of close and contained point is clear.Then the world coordinates each put in quadrangle is corresponded into X-ray angiographic image number
Image coordinate in；Assuming that origin is P0, contrastographic picture pel spacing is d, world coordinates Pw, then image coordinate Pt
Coordinate can acquire as follows：
Step S5：According to image coordinate, searching and the most similar pixel of the image coordinate obtain corresponding pixel ash
Degree, obtains gray areas Sn(n=1,2n)；
Step S6：It will obtain gray areas SnMerge along straight line (central axes), obtain blood vessel stretches reconstruction image.
It should be understood that the above-mentioned specific implementation mode of the present invention is used only for exemplary illustration or explains the present invention's
Principle, but not to limit the present invention.Therefore, that is done without departing from the spirit and scope of the present invention is any
Modification, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.In addition, appended claims purport of the present invention
Covering the whole variations fallen into attached claim scope and boundary or this range and the equivalent form on boundary and is repairing
Change example.
Claims (5)
1. a kind of blood vessel based on X-ray angiographic image stretches method for reconstructing, which is characterized in that include the following steps：
S01：Obtain X-ray angiogram data and vessel centerline data；
S02：Certain distance N is extended respectively along normal both forward and reverse directions to point each of in vessel centerline, both forward and reverse directions are each got
To a point An、Bn；
S03：By four adjacent point An、Bn、An+1、Bn+1Be connected quadrangularly；
S04：Obtain image coordinate of each quadrangle in X-ray angiogram data；
S05：Corresponding pixel grey scale is obtained according to image coordinate, obtains gray areas；
S06：Obtained gray areas is merged along straight line, obtain blood vessel stretches reconstruction image.
2. the blood vessel according to claim 1 based on X-ray angiographic image stretches method for reconstructing, which is characterized in that
In the step S02,Wherein,For normal direction vector, PnCentered on line
Point.
3. the blood vessel according to claim 1 based on X-ray angiographic image stretches method for reconstructing, which is characterized in that
The step S04 further includes not having to insert in clear gray value point in X-ray angiogram data using image interpolation method
Enter to have the point of clear gray value.
4. the blood vessel according to claim 3 based on X-ray angiographic image stretches method for reconstructing, which is characterized in that
Described image interpolation method is closest first method, including：In the known gray scale point of neighbouring four of gray scale point to be asked, distance is waited for
The gray value for the point for asking gray scale point nearest gives gray scale point to be asked.
5. the blood vessel according to claim 1 based on X-ray angiographic image stretches method for reconstructing, which is characterized in that
Described image coordinate PtCoordinate acquired by following formula：
P0For origin, d is contrastographic picture pel spacing, PwFor world coordinates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810612090.8A CN108805815B (en) | 2018-06-14 | 2018-06-14 | Blood vessel straightening reconstruction method based on X-ray angiography image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810612090.8A CN108805815B (en) | 2018-06-14 | 2018-06-14 | Blood vessel straightening reconstruction method based on X-ray angiography image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108805815A true CN108805815A (en) | 2018-11-13 |
| CN108805815B CN108805815B (en) | 2023-02-17 |
Family
ID=64085882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810612090.8A Active CN108805815B (en) | 2018-06-14 | 2018-06-14 | Blood vessel straightening reconstruction method based on X-ray angiography image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108805815B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110033442A (en) * | 2019-04-01 | 2019-07-19 | 数坤（北京）网络科技有限公司 | A kind of angiosteosis method for detecting area and system based on analysis line drawing |
| CN112116615A (en) * | 2019-11-19 | 2020-12-22 | 苏州润迈德医疗科技有限公司 | Method and device for obtaining blood vessel contour line according to blood vessel central line |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060079746A1 (en) * | 2004-10-11 | 2006-04-13 | Perret Florence M | Apparatus and method for analysis of tissue classes along tubular structures |
| US20060269108A1 (en) * | 2005-02-07 | 2006-11-30 | Viswanathan Raju R | Registration of three dimensional image data to 2D-image-derived data |
| US20110096996A1 (en) * | 2009-10-28 | 2011-04-28 | Siemens Corporation | Method for processing vascular structure images |
| US20130064343A1 (en) * | 2011-09-13 | 2013-03-14 | Boudewijn J.A. Verstraelen | Method and Apparatus for Determining Three-Dimensional Reconstruction of an Object |
| KR101472558B1 (en) * | 2013-10-04 | 2014-12-16 | 원광대학교산학협력단 | The system and method for automatic segmentation of lung, bronchus, pulmonary vessels images from thorax ct images |
| CN105913479A (en) * | 2016-04-05 | 2016-08-31 | 苏州润心医疗科技有限公司 | Vascular curved surface reconstruction method based on heart CT image |
| US20170018116A1 (en) * | 2015-07-14 | 2017-01-19 | Siemens Medical Solutions Usa, Inc. | 3-d vessel tree surface reconstruction method |
| CN107451406A (en) * | 2017-07-28 | 2017-12-08 | 海纳医信（北京）软件科技有限责任公司 | Vessels analysis method, apparatus, storage medium and processor |
-
2018
- 2018-06-14 CN CN201810612090.8A patent/CN108805815B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060079746A1 (en) * | 2004-10-11 | 2006-04-13 | Perret Florence M | Apparatus and method for analysis of tissue classes along tubular structures |
| US20060269108A1 (en) * | 2005-02-07 | 2006-11-30 | Viswanathan Raju R | Registration of three dimensional image data to 2D-image-derived data |
| US20110096996A1 (en) * | 2009-10-28 | 2011-04-28 | Siemens Corporation | Method for processing vascular structure images |
| US20130064343A1 (en) * | 2011-09-13 | 2013-03-14 | Boudewijn J.A. Verstraelen | Method and Apparatus for Determining Three-Dimensional Reconstruction of an Object |
| KR101472558B1 (en) * | 2013-10-04 | 2014-12-16 | 원광대학교산학협력단 | The system and method for automatic segmentation of lung, bronchus, pulmonary vessels images from thorax ct images |
| US20170018116A1 (en) * | 2015-07-14 | 2017-01-19 | Siemens Medical Solutions Usa, Inc. | 3-d vessel tree surface reconstruction method |
| CN105913479A (en) * | 2016-04-05 | 2016-08-31 | 苏州润心医疗科技有限公司 | Vascular curved surface reconstruction method based on heart CT image |
| CN107451406A (en) * | 2017-07-28 | 2017-12-08 | 海纳医信（北京）软件科技有限责任公司 | Vessels analysis method, apparatus, storage medium and processor |
Non-Patent Citations (3)
| Title |
|---|
| PAPITHA JAYABAL，ETC.: "Rician Noise Removal and Straightening of Blood Vessel in MR Images", 《JOURNAL OF BIOMEDICAL ENGINEERING AND MEDICAL IMAGING》 * |
| 肖若秀等: "应用于造影图像的血管中心线全自动提取方法", 《计算机辅助设计与图形学学报》 * |
| 郝爽: "三维血管拉直算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110033442A (en) * | 2019-04-01 | 2019-07-19 | 数坤（北京）网络科技有限公司 | A kind of angiosteosis method for detecting area and system based on analysis line drawing |
| CN112116615A (en) * | 2019-11-19 | 2020-12-22 | 苏州润迈德医疗科技有限公司 | Method and device for obtaining blood vessel contour line according to blood vessel central line |
| WO2021097817A1 (en) * | 2019-11-19 | 2021-05-27 | 苏州润迈德医疗科技有限公司 | Method and apparatus for acquiring contour lines of blood vessel according to center line of blood vessel |
| CN112116615B (en) * | 2019-11-19 | 2023-12-05 | 苏州润迈德医疗科技有限公司 | Method and device for acquiring blood vessel contour line according to blood vessel center line |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108805815B (en) | 2023-02-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2550542C2 (en) | Method and device for shaping computer tomographic images using geometries with offset detector | |
| US7546154B2 (en) | Method and apparatus for automatic detection of anomalies in vessel structures | |
| US11607189B2 (en) | Method and apparatus for processing blood vessel image, storage medium, and imaging device | |
| JP2004243117A (en) | Method for obtaining physical parameters of physiological structure | |
| Prause et al. | Binary reconstruction of the heart chambers from biplane angiographic image sequences | |
| Almazroa et al. | Agreement among ophthalmologists in marking the optic disc and optic cup in fundus images | |
| JP5952564B2 (en) | Image processing apparatus and image processing method | |
| JP2002177254A (en) | Method and device for restructuring dynamic three- dimensional image of subject covered with contrast medium | |
| JP3927488B2 (en) | Image processing method, image processing apparatus, and program | |
| CN108320314A (en) | A kind of image processing method and device based on the cross-section images of CT, Bone images display system | |
| CN108805815A (en) | Blood vessel based on X-ray angiographic image stretches method for reconstructing | |
| US7706589B2 (en) | Analysis of a multi-dimensional structure | |
| JP4701063B2 (en) | Tomography system | |
| JP2007159934A (en) | Comparative diagnostic reading supporting apparatus | |
| US10080492B2 (en) | Method and system for motion artefacts removal in optical coherence tomography | |
| JPH07239933A (en) | Image processing unit | |
| JP3630812B2 (en) | X-ray diagnostic imaging equipment | |
| JP2022547463A (en) | Confidence Map for Limited Angle Artifact Mitigation Based on Neural Networks in Cone-Beam CT | |
| US20220202376A1 (en) | Medical imaging apparatus including biological signal processing system, medical imaging system, and biological signal processing method | |
| US7116808B2 (en) | Method for producing an image sequence from volume datasets | |
| EP3781014B1 (en) | Post-processing method to improve lso-based tracking in oct | |
| WO2014096835A1 (en) | Improvements in and relating to ophthalmoscopes | |
| US11145094B2 (en) | Image reconstruction apparatus and image reconstruction method | |
| JP6103829B2 (en) | Image diagnosis support apparatus, image diagnosis support method, image diagnosis support system, and program | |
| CN113538298B (en) | Focal high-definition image reconstruction method, device, equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |