CN108805815B - Blood vessel straightening reconstruction method based on X-ray angiography image - Google Patents
Blood vessel straightening reconstruction method based on X-ray angiography image Download PDFInfo
- Publication number
- CN108805815B CN108805815B CN201810612090.8A CN201810612090A CN108805815B CN 108805815 B CN108805815 B CN 108805815B CN 201810612090 A CN201810612090 A CN 201810612090A CN 108805815 B CN108805815 B CN 108805815B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- image
- point
- ray angiography
- gray scale
- 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
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 39
- 238000002583 angiography Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000005096 rolling process Methods 0.000 claims description 2
- 238000002601 radiography Methods 0.000 description 7
- 238000005259 measurement Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 230000036285 pathological change Effects 0.000 description 3
- 231100000915 pathological change Toxicity 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction 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
- 208000037804 stenosis Diseases 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
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 blood vessel straightening and reconstructing method based on an X-ray angiography image, which comprises the following steps: obtaining X-ray angiography image data and blood vessel center line data; for each point on the center line of the blood vessel, respectively extending a certain distance N along the positive and negative directions of the normal line, and respectively obtaining a point in the positive and negative directions、(ii) a Four adjacent points、、、Connected to form a quadrangle; acquiring image coordinates of each quadrangle in X-ray angiography image data; acquiring corresponding pixel gray scale according to the image coordinates to obtain a gray scale area; and combining the obtained gray regions along a straight line to obtain a straightened reconstructed image of the blood vessel. The blood vessel can be straightened and reconstructed along the central line of the blood vessel, and a doctor can be helped to analyze the blood vessel qualitatively or quantitatively quickly and accurately.
Description
Technical Field
The invention relates to the technical field of medical image processing, in particular to a blood vessel straightening reconstruction method based on an X-ray angiography image.
Background
At present, most hospitals in China use an X-ray single-arm radiography system to carry out X-ray radiography on a patient, and radiography images corresponding to different radiography angles are obtained by rotating an radiography arm. The current mainstream diagnostic method of X-ray angiography is that doctors perform quantitative or qualitative analysis on blood vessels in an angiography image during or after operation through a simple two-point distance measuring tool through self experience or in image playing software matched with an X-ray single-arm angiography system.
However, there are the following drawbacks:
1. the personal experience and level of doctors can be differentiated, and the contrast images are different from images generated by other imaging equipment and need to have the experience of processing the contrast images for a long time;
2. at present, in image processing software, measurement tools carried by the image processing software are simple and easy, can meet measurement requirements, but cannot meet the requirement of accurate measurement, and the simple measurement tools have non-negligible errors in measurement of complex curved contrast image vessels;
3. uncertainty of manual interaction in the analysis process;
4. due to the characteristics of the X-ray angiography image, the contrast agent flows through the blood vessel in a specific time, so that a doctor can observe the dynamic lumen condition of the blood vessel, but due to dynamic observation, the difficulty of observation in time is possibly increased, and further partial pathological changes on the blood vessel are missed;
5. all the judgment results of the blood vessels are artificially subjective judgment, and data conforming to a systematic and uniform rule is not generated.
Disclosure of Invention
In order to solve the technical problems, the invention aims to: the blood vessel straightening reconstruction method based on the X-ray angiography image can straighten and reconstruct the blood vessel along the central line of the blood vessel and is beneficial to a doctor to analyze the blood vessel qualitatively or quantitatively quickly and accurately.
The technical scheme of the invention is as follows:
a blood vessel straightening reconstruction method based on an X-ray angiography image comprises the following steps:
s01: obtaining X-ray angiography image data and blood vessel center line data;
s02: each point on the center line of the blood vessel extends a certain distance N along the positive and negative directions of the normal line respectively to obtain a point A n 、B n ;
S03: four adjacent points A n 、B n 、A n+1 、B n+1 Connected to form a quadrangle;
s04: acquiring image coordinates of each quadrangle in X-ray angiography image data;
s05: acquiring corresponding pixel gray scale according to the image coordinates to obtain a gray scale area;
s06: and combining the obtained gray regions along a straight line to obtain a straightened reconstructed image of the blood vessel.
In a preferred embodiment, in step S02,wherein,is a normal direction vector, P n As a point on the centerline.
In a preferred embodiment, the step S04 further includes inserting a point with a clear gray value into a point without a clear gray value in the X-ray angiography image data by using an image interpolation method.
In a preferred technical solution, the image interpolation method is a nearest neighbor method, and includes: and in four adjacent known gray scale points of the gray scale point to be solved, assigning the gray scale value of the point closest to the gray scale point to be solved.
Preferred technical solutionIn said image coordinate P t The coordinates of (a) are obtained by the following formula:
P 0 is the origin coordinate, d is the contrast image pixel pitch, P w Is a world coordinate.
Compared with the prior art, the invention has the advantages that:
the method can straighten and reconstruct the blood vessel along the central line of the blood vessel, can more intuitively present the blood vessel to a doctor, and is favorable for the doctor to realize multi-layer and multi-angle observation of the blood vessel during clinical diagnosis. Meanwhile, after the blood vessel is straightened, the stenosis condition and the plaque condition of the blood vessel can be observed more easily, and a doctor can analyze the blood vessel qualitatively or quantitatively more quickly and accurately to find the position of a pathological change and measure the pathological change condition. Therefore, the blood vessel straightening reconstruction is an important component of X-ray angiography analysis and has important significance for analyzing diseased blood vessels of patients.
Drawings
The invention is further described below with reference to the following figures and examples:
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a schematic representation of an X-ray angiographic image;
fig. 3 is a graph showing the effect of straightening the blood vessel in fig. 2.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the accompanying drawings in combination with the embodiments. It is to be understood that these descriptions are only illustrative and are not intended to limit the scope of the present invention. Moreover, in the following description, descriptions of well-known structures and techniques are omitted so as to not unnecessarily obscure the concepts of the present invention.
As shown in fig. 1, the method for reconstructing blood vessel by straightening based on X-ray angiography image of the present invention includes the following steps.
Step S1: acquiring contrast image data and vessel centerline data;
acquisition of contrast image data has three approaches: the first mode adopts direct connection with an X-ray radiography machine and carries out transmission of radiography image data in a wired mode; the second mode is connected with a medical image storage and transmission system of the cardiology department of the hospital, and obtains contrast image data from the system in a wired mode; the third method is to copy the contrast image data with a storage medium on a controller of the X-ray contrast machine or a server of a medical image storage and transmission system connected to the X-ray contrast machine.
The acquisition mode of the blood vessel central line data is mainly a grid method based on unit circle rolling tracking. The method comprises the steps of rasterizing a region to be calculated, calculating the minimum distance from the center point of each grid to two side boundaries, setting an allowable range of the absolute value of the difference between the two minimum distances, and if the result is in the allowable range, indicating that the point is the center point of the grid region. Repeating the steps to obtain two central points A1 and A2, drawing a circle by taking the radius as | A1A2|, finding two points on the circumference, wherein the minimum distances from the two points to the boundaries of the two sides are equal, wherein a point A3 which is not coincident with the point A1 is the other central point, repeating the steps until all the central points meeting the length are found, and connecting to obtain the central line.
Step S2: for each point on the center line of the blood vessel, respectively extending the distance N along the positive and negative directions of the normal line to obtain a point A n (n=1,2···n)、B n (n =1,2 · n); specifically, assume that the line direction vector isThen any point P on the centerline n (n =0,1, \8230;, n-1) the coordinates of the corresponding point points extending in the line direction can be found by the following formula:
and step S3: four adjacent points A n (n=1,2···n)、B n (n=1,2···n)、A n+1 (n=1, 2···n)、B n+1 (n =1,2 · n) connected in a quadrilateral; since the display of the calculated center line is required to be as smooth as possible, the precision of the distance between the adjacent center points is specified to be small, and since the distance between the normal to the center point and the forward and backward directions is fixed N, the quadrangle formed by connecting the adjacent four points can be nearly rectangular, and it can be considered that several rectangles having continuity can be obtained by repeating several times.
And step S4: because the resolution of the rendered quadrangle is different from the resolution of the actual X-ray angiography image data, the X-ray angiography image data contains some points with unclear gray values, and in order to enable the points to correspond to the points on the quadrangle as much as possible and to have definite gray values, an image interpolation method is adopted to achieve the purpose. The image interpolation method is used for inserting more points with definite gray scales into the X-ray angiography image data, for example, but not limited to, the nearest neighbor method is adopted, and in four known gray scale points adjacent to the gray scale point to be solved, the gray scale of the point closest to the gray scale point to be solved is assigned to the gray scale point to be solved, so that the resolution of all the X-ray angiography images is close to the resolution of the quadrangle drawn in the previous step, and the gray scale of the included points is definite. Then, the world coordinates of each point in the quadrangle correspond to the image coordinates in the X-ray angiography image data; assuming origin coordinates P 0 The pixel pitch of the contrast image is d, and the world coordinate is P w Then image coordinate P t The coordinates of (c) can be obtained by the following formula:
step S5: according to the image coordinates, searching a pixel point closest to the image coordinates, acquiring the corresponding pixel gray level, and obtaining a gray level area S n (n=1,2···n);
Step S6: will be provided withObtaining a gray scale region S n And merging along a straight line (central axis) to obtain a straightened reconstructed image of the blood vessel.
It is to be understood that the above-described embodiments of the present invention are merely illustrative of or explaining the principles of the invention and are not to be construed as limiting the invention. Therefore, any modification, equivalent replacement, improvement and the like made without departing from the spirit and scope of the present invention should be included in the protection scope of the present invention. Further, it is intended that the appended claims cover all such variations and modifications as fall within the scope and boundaries of the appended claims or the equivalents of such scope and boundaries.
Claims (5)
1. A blood vessel straightening reconstruction method based on an X-ray angiography image is characterized by comprising the following steps:
s01: acquiring X-ray angiography image data and blood vessel center line data based on a unit circle rolling tracking grid method;
s02: for each point on the center line of the blood vessel, respectively extending a certain distance N along the positive and negative directions of the normal line, and respectively obtaining a point A in the positive and negative directions n 、B n ;
S03: four adjacent points A n 、B n 、A n+1 、B n+1 Connected to form a quadrangle;
s04: acquiring image coordinates of each quadrangle in X-ray angiography image data;
s05: acquiring corresponding pixel gray scale according to the image coordinates to obtain a gray scale area;
s06: and combining the obtained gray areas along the central axis to obtain a straightened reconstructed image of the blood vessel.
3. The method for reconstructing blood vessel by straightening based on X-ray angiography image according to claim 1, wherein the step S04 further comprises inserting points with definite gray values into the points without definite gray values in the X-ray angiography image data by image interpolation.
4. The method of claim 3, wherein the image interpolation is nearest neighbor, comprising: and in four adjacent known gray scale points of the gray scale point to be solved, the gray scale value of the point closest to the gray scale point to be solved is assigned to the gray scale point to be solved.
5. Method for X-ray angiography image based vessel straightening reconstruction according to claim 1, characterized in that the image coordinates P t The coordinates of (c) are obtained by the following formula:
P 0 is the origin coordinate, d is the contrast image pixel pitch, P w Is a world coordinate.
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 CN108805815A (en) | 2018-11-13 |
CN108805815B true 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) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110033442B (en) * | 2019-04-01 | 2020-12-25 | 数坤(北京)网络科技有限公司 | Vascular calcification area detection method and system based on analysis line extraction |
CN112116615B (en) * | 2019-11-19 | 2023-12-05 | 苏州润迈德医疗科技有限公司 | Method and device for acquiring blood vessel contour line according to blood vessel center line |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
Family Cites Families (6)
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 |
US7756308B2 (en) * | 2005-02-07 | 2010-07-13 | Stereotaxis, Inc. | Registration of three dimensional image data to 2D-image-derived data |
US8934690B2 (en) * | 2009-10-28 | 2015-01-13 | Siemens Aktiengesellschaft | Method for processing vascular structure images |
EP2570079B1 (en) * | 2011-09-13 | 2017-06-14 | Pie Medical Imaging BV | Method and apparatus for determining optimal 3D reconstruction of an object |
US9934566B2 (en) * | 2015-07-14 | 2018-04-03 | Siemens Healthcare Gmbh | 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 (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
Non-Patent Citations (1)
Title |
---|
三维血管拉直算法研究;郝爽;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130315;I138-1151 * |
Also Published As
Publication number | Publication date |
---|---|
CN108805815A (en) | 2018-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
NL1032927C2 (en) | Method and system for automatically determining zones in a scanned object. | |
JP6334141B2 (en) | Method and apparatus for navigating a CT scan by a marker | |
JP5129480B2 (en) | System for performing three-dimensional reconstruction of tubular organ and method for operating blood vessel imaging device | |
EP2570079B1 (en) | Method and apparatus for determining optimal 3D reconstruction of an object | |
US10052032B2 (en) | Stenosis therapy planning | |
JP2004243117A (en) | Method for obtaining physical parameters of physiological structure | |
Galassi et al. | 3D reconstruction of coronary arteries from 2D angiographic projections using non-uniform rational basis splines (NURBS) for accurate modelling of coronary stenoses | |
US20180204339A1 (en) | Blood vessel image processing apparatus, blood vessel image processing program, and blood vessel image processing method | |
CN112419484B (en) | Three-dimensional vascular synthesis method, system, coronary artery analysis system and storage medium | |
US10555712B2 (en) | Segmenting an angiography using an existing three-dimensional reconstruction | |
Hoffmann et al. | A system for determination of 3D vessel tree centerlines from biplane images | |
US10769763B2 (en) | Method for reconstructing a reconstruction data set of a vessel segment | |
CN111476791B (en) | Image processing method, image processing apparatus, and non-transitory computer readable medium | |
JP2004320771A (en) | Method for performing digital subtraction angiography | |
Sieren et al. | Automated segmentation and quantification of the healthy and diseased aorta in CT angiographies using a dedicated deep learning approach | |
CN108805815B (en) | Blood vessel straightening reconstruction method based on X-ray angiography image | |
JP4444100B2 (en) | Multidimensional structure analysis method | |
JP2002345796A (en) | Method and device for determining magnification of radiographic image | |
KR20100025431A (en) | Method for detection of hepatic tumors using registration of multi-phase liver ct images | |
US7457658B2 (en) | Algorithm for accurate three-dimensional reconstruction of non-linear implanted medical devices in VIVO | |
CN107773261A (en) | X-ray shooting with overlapping plan information | |
US7640136B2 (en) | System and method for determination of object location for calibration using patient data | |
CN117084706B (en) | Three-source CT-based locating plate scanning and 3D contour estimation method thereof | |
US20240127450A1 (en) | Medical image processing apparatus and non-transitory computer readable medium | |
Hsieh | Local motion detection, characterization, and quantification for x-ray CT |
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 |