CN1887232A - Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration - Google Patents
Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration Download PDFInfo
- Publication number
- CN1887232A CN1887232A CN 200610029686 CN200610029686A CN1887232A CN 1887232 A CN1887232 A CN 1887232A CN 200610029686 CN200610029686 CN 200610029686 CN 200610029686 A CN200610029686 A CN 200610029686A CN 1887232 A CN1887232 A CN 1887232A
- Authority
- CN
- China
- Prior art keywords
- characteristic point
- registration
- section
- cross
- virtual
- 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
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
The method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration belongs to the field of image treating technology. The method includes finding out the central path of the whole tubular organ and obtaining the vertical cross sections along the central path; finding out the registration source characteristic points in the cross sections and calculating the registration target characteristic points based on these registration source characteristic points; and final 3D non-rigid characteristic point registration based on these characteristic points to obtain virtual inflated expansion organ. The present invention can obtain virtual inflated expansion results to raise diagnosis accuracy.
Description
Technical field
What the present invention relates to is a kind of method of technical field of image processing, is specifically related to a kind of method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration.
Technical background
At present, utilize medical imaging datas such as CT, MRI, the inwall of body interior conduit internal organs such as colon, esophagus, trachea etc. is carried out three-dimensional visualization is rebuild and what carry out that clinical diagnosis mainly adopts is virtual endoscopic techniques.This technical modelling advance in the traditional optical endoscope tract, to observe and diagnose in the visual angle of endoscope optical.The advantage of this technology is that it is a kind of method of not having intrusive mood, to the patient bring painful less.
Find through literature search prior art, in " development of the virtual scope of 3 D human body and applied research " literary composition that the leaf Kien Giang is delivered, introduced the characteristics and main application of virtual endoscopic techniques on " medical equipment information " (2002 the 6th phase 48-50 pages or leaves).Pipeline internal organs in human body under in the raw, as colon, esophagus etc., a lot of places can be the flat state that contracts, and can influence the reconstruction of virtual endoscopic techniques and visual like this.In this literary composition, solution to this problem is when carrying out CT or MRI inspection colon to be carried out actual inflation.The deficiency of this method is: then powerless for some internal organs that can not inflate; The process of inflation can be brought discomfort to the patient.
Summary of the invention
The present invention is directed to existing virtual endoscopic techniques and the pipeline internal organs that are the flat state that contracts are carried out the deficiency of visual reconstruction aspect, a kind of method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration has been proposed, make before the patient who need is taked the virtual endoscope diagnosis carries out CT or MRI inspection, need not patient's internal organs are inflated, reduce patient's misery.Simultaneously, because of inflating the internal organs that make that effect of visualization is very poor, peep effect of visualization in also obtaining preferably for original.
The present invention realizes by following technical method, the present invention at first finds out the center path of whole pipe internal organs, obtain and vertical each cross section of center path along center path, in each cross section, find out registration source characteristic point then respectively, according to these characteristic points, calculate registration target characteristic point, more at last based on characteristic point, carry out three-dimensional non-rigid registration, thereby obtain the internal organs behind the virtual inflation.
The described center path of finding out the whole pipe internal organs, obtain and vertical each cross section of center path along center path, be meant: to cutting apart the image of the pipeline internal organs after the reconstruction, seek the method for shortest path after method by the topology refinement or the range conversion, find out the center path of whole internal organs.Along this paths, just get one every a small distance and follow this paths then in the vertical cross section of the tangential direction of this point.
The described registration source characteristic point of in each cross section, finding out respectively, be meant: at first in each cross section, find out the edge of internal organs, on the edge, be the center then with the center path point, get a point every an angle, these points that will take out are at last followed according to needs and are carried out smoothing processing with marginal point on every side respectively, can obtain registration source characteristic point.
Described according to these characteristic points, calculate registration target characteristic point, be meant: each the registration source characteristic point with in each cross section corresponds on the position behind the virtual inflation.These points after the correspondence promptly are registration target characteristic points.
Position behind the described virtual inflation is meant: on same cross section, with the point on the circumference of a circle of this registration source characteristic point correspondence.The center of circle of this circle is the intersection point that center path is followed this cross section, and radius is that the long sum of camber line is divided by 2 π between each registration source characteristic point, and the value that calculates with near cross section is made the smothing filtering gained again.
Described based on characteristic point, carry out three-dimensional non-rigid registration, be meant: the registration source characteristic point and the registration target characteristic point that utilize said method to try to achieve, method according to characteristic point non-rigid transformation registration, obtain that each puts the corresponding relation of target image in the source images, thereby calculate target image.
The method of described characteristic point non-rigid transformation registration, be meant: according to the corresponding relation between known characteristic point, methods such as application of thin dish batten, approximate thin dish batten, volume batten or elastomer batten obtain any some corresponding relation to target image in the image.
The present invention uses above method, has realized the flat pipeline internal organs that contract are carried out virtual inflation, makes these internal organs under the situation without true inflation, still can realize visual under the technology of virtual endoscope well.And, choose under the appropriate situation in each parameter, after process this method is carried out virtual inflation, still can keep the original tiny anatomical structure of internal organs inwall substantially, thereby guarantee the accuracy of medical diagnosis.
Advantage of the present invention is: (1) avoided clinical in real gas replenishment process, reduce patient's misery.(2) make the internal organs that originally can't truly inflate,, also can obtain the effect behind the inflation as pipeline internal organs such as blood vessel, fallopian tube.(3) can increase the visual field of method for visualizing such as follow-up virtual endoscope, improve diagnosis effect.
The specific embodiment
The inventive method concrete steps are as follows: (1) is read in primitive medicine CT or MRI data and is carried out three-dimensional reconstruction.(2) carry out data and cut apart, obtain interested pipeline internal organs part.(3) obtain center path.(4) obtain and vertical each cross section of center path.(5) find out registration source characteristic point.(6) try to achieve registration target characteristic point.(7) carry out the characteristic point non-rigid registration.(8) calculate target image behind the virtual inflation.Below in conjunction with a specific embodiment technical scheme of the present invention is described in further detail.
This embodiment adopts the CT slice of data of human body esophagus as sample, and the implementation procedure of whole invention is as follows:
1. read in human chest CT data, and these successive section sequences are reconstructed into a three-dimensional data.
2. volume data is cut apart, adopted the algorithm of region growing, obtain the three-dimensional binaryzation volume data of esophagus.
3. utilize the method for morphology refinement, obtain the skeleton of esophagus, and, make it enough level and smooth this skeleton application of low-pass filters.This promptly can be used as the center path of this esophagus.
4. from then on the summit in path begins, and gets a cross section every 2mm, and this cross section is perpendicular to the tangential direction of path at this point.
5. in each cross section, obtain each marginal point of esophagus.With the center path center then, from horizontal positive direction, every marginal point of 5 ° of taking-ups, as registration source characteristic point.
6. in each cross section, try to achieve the distance of each characteristic point, and utilize following formula to center path point:
Obtain j the expansion radius R in the cross section
jWherein, r
iBe the distance of i characteristic point to center path, n is the sum of characteristic point in this cross section.
Then, to R
jCarry out Filtering Processing, mode as shown in the formula:
Wherein, weight w
hBy a low pass filter decision, and satisfy
In the present embodiment, may be selected to be ramp filter, get m=9.
Successively in each cross section, be the center of circle again with the center path point, with
Value is justified for radius.From horizontal positive direction, on this circle, get characteristic point then every 5 °.These characteristic points are promptly as registration target characteristic point.
7. to known source and target characteristic point, try to achieve each to the conversion vector between the characteristic point.Utilize thin dish spline interpolation again, draw every bit in the original image to the conversion vector of target image.
8. utilize the conversion vector of having tried to achieve, original image is transformed in the object space.The method that cooperates neighbour's interpolation again calculates the 3-D view volume data after the virtual expansion.
Such embodiment can be realized making originally the esophagus of comparatively flat wrinkle is the tympanites state, and not lose the original tiny anatomical structure of esophagus, can conveniently carry out further visual and other processing operations.
Claims (7)
1, a kind of method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration, it is characterized in that: the center path of at first finding out the whole pipe internal organs, obtain and vertical each cross section of center path along center path, in each cross section, find out registration source characteristic point then respectively, again according to these characteristic points, calculate registration target characteristic point, carry out three-dimensional non-rigid registration based on characteristic point at last, thereby obtain the internal organs behind the virtual inflation.
2, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 1, it is characterized in that, the described center path of finding out the whole pipe internal organs, obtain and vertical each cross section of center path along center path, be meant: to cutting apart the image of the pipeline internal organs after the reconstruction, seek the method for shortest path after method by the topology refinement or the range conversion, find out the center path of whole internal organs, along this paths, just get one every a small distance and follow this paths then in the vertical cross section of the tangential direction of this point.
3, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 1, it is characterized in that, the described registration source characteristic point of in each cross section, finding out respectively, be meant: at first in each cross section, finding out the edge of internal organs, is the center with the center path point on the edge then, gets a point every an angle, these points that will take out are at last followed according to needs and are carried out smoothing processing with marginal point on every side respectively, obtain registration source characteristic point.
4, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 1, it is characterized in that, described according to these characteristic points, calculate registration target characteristic point, be meant: with each the registration source characteristic point in each cross section, correspond on the position behind the virtual inflation, these points after the correspondence are exactly registration target characteristic point.
5, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 4, it is characterized in that, position behind the described virtual inflation, be meant: on same cross section, with the point on the circumference of a circle of this registration source characteristic point correspondence, the center of circle of this circle is the intersection point that center path is followed this cross section, and radius is that the long sum of camber line is divided by 2 π between each registration source characteristic point, and the value that calculates with near cross section is made the smothing filtering gained again.
6, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 1, it is characterized in that, describedly carry out three-dimensional non-rigid registration based on characteristic point, be meant: the registration source characteristic point and the registration target characteristic point that utilize said method to try to achieve, method according to characteristic point non-rigid transformation registration, obtain that each puts the corresponding relation of target image in the source images, thereby calculate target image.
7, the method that realizes virtual inflated expansion of tubular organ with the characteristic point non-rigid registration according to claim 6, it is characterized in that, the method of described characteristic point non-rigid transformation registration, be meant: according to the corresponding relation between known characteristic point, application of thin dish batten, approximate thin dish batten, volume batten or elastomer spline method obtain any some corresponding relation to target image in the image.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006100296862A CN100434044C (en) | 2006-08-03 | 2006-08-03 | Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CNB2006100296862A CN100434044C (en) | 2006-08-03 | 2006-08-03 | Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1887232A true CN1887232A (en) | 2007-01-03 |
| CN100434044C CN100434044C (en) | 2008-11-19 |
Family
ID=37576511
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNB2006100296862A Expired - Fee Related CN100434044C (en) | 2006-08-03 | 2006-08-03 | Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN100434044C (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102034275A (en) * | 2010-12-10 | 2011-04-27 | 中国人民解放军国防科学技术大学 | Large-scale deformation-oriented non-rigid registration method |
| CN101661627B (en) * | 2008-08-29 | 2012-07-04 | 株式会社东芝 | Image processing apparatus and image processing method |
| CN102846331A (en) * | 2011-06-30 | 2013-01-02 | 上海西门子医疗器械有限公司 | X-ray computed tomography scanning system and method |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5611025A (en) * | 1994-11-23 | 1997-03-11 | General Electric Company | Virtual internal cavity inspection system |
| JP2002238887A (en) * | 2001-02-19 | 2002-08-27 | National Cancer Center-Japan | Virtual endoscope |
| US7417636B2 (en) * | 2003-05-08 | 2008-08-26 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for automatic setting of rendering parameter for virtual endoscopy |
| US7457444B2 (en) * | 2003-05-14 | 2008-11-25 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for fast automatic centerline extraction for virtual endoscopy |
| US20050152588A1 (en) * | 2003-10-28 | 2005-07-14 | University Of Chicago | Method for virtual endoscopic visualization of the colon by shape-scale signatures, centerlining, and computerized detection of masses |
-
2006
- 2006-08-03 CN CNB2006100296862A patent/CN100434044C/en not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101661627B (en) * | 2008-08-29 | 2012-07-04 | 株式会社东芝 | Image processing apparatus and image processing method |
| CN102034275A (en) * | 2010-12-10 | 2011-04-27 | 中国人民解放军国防科学技术大学 | Large-scale deformation-oriented non-rigid registration method |
| CN102034275B (en) * | 2010-12-10 | 2013-12-18 | 中国人民解放军国防科学技术大学 | Large-scale deformation-oriented non-rigid registration method |
| CN102846331A (en) * | 2011-06-30 | 2013-01-02 | 上海西门子医疗器械有限公司 | X-ray computed tomography scanning system and method |
| CN102846331B (en) * | 2011-06-30 | 2015-04-08 | 上海西门子医疗器械有限公司 | X-ray computed tomography scanning system and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN100434044C (en) | 2008-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110969622B (en) | Image processing method and system for assisting pneumonia diagnosis | |
| CN112150524B (en) | Two-dimensional and three-dimensional medical image registration method and system based on deep learning | |
| JP4899068B2 (en) | Medical image observation support device | |
| US8579800B2 (en) | Systematic chromoendoscopy and chromocolonoscopy as a novel systematic method to examine organs with endoscopic techniques | |
| CN113066026B (en) | Endoscope image smoke purification method based on deep neural network | |
| WO2017027638A1 (en) | 3d reconstruction and registration of endoscopic data | |
| CN113112609A (en) | Navigation method and system for lung biopsy bronchoscope | |
| WO2014081725A2 (en) | Electromagnetic sensor integration with ultrathin scanning fiber endoscope | |
| WO2007129616A1 (en) | Insertion assist system of endoscope and insertion assist method of endoscope | |
| WO2005041761A1 (en) | Insertion support system | |
| CN107204045A (en) | Virtual endoscope system based on CT images | |
| CN112700551A (en) | Virtual choledochoscope interventional operation planning method, device, equipment and storage medium | |
| CN100434044C (en) | Method of realizing virtual inflated expansion of tubular organ through non-rigid characteristic point registration | |
| CN116421313A (en) | Augmented reality fusion method in navigation of lung tumor resection operation under thoracoscope | |
| Gu et al. | Vision–kinematics interaction for robotic-assisted bronchoscopy navigation | |
| CN115830016A (en) | Medical image registration model training method and equipment | |
| Robb | Virtual endoscopy: evaluation using the visible human datasets and comparison with real endoscopy in patients | |
| CN100454340C (en) | Visual method for virtual incising tubular organ | |
| Khmelinskii et al. | Atlas-based organ & bone approximation for ex-vivo μMRI mouse data: a pilot study | |
| CN115953377A (en) | Digestive tract ultrasonic endoscope image fusion method and system | |
| CN116993805A (en) | Intraoperative residual organ volume estimation system oriented to operation planning assistance | |
| CN100418478C (en) | Virtual endoscope surface color mapping method based on blood flow imaging | |
| CN100561519C (en) | Visual method for virtual incising tubular organ based on traversal excision vector | |
| CN115919464A (en) | Tumor positioning method, system and device and tumor development prediction method | |
| Wang et al. | A monocular SLAM system based on SIFT features for gastroscope tracking |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| C17 | Cessation of patent right | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20081119 Termination date: 20110803 |