CN1875380A - Rendering for coronary visualization - Google Patents
Rendering for coronary visualization Download PDFInfo
- Publication number
- CN1875380A CN1875380A CNA2004800324109A CN200480032410A CN1875380A CN 1875380 A CN1875380 A CN 1875380A CN A2004800324109 A CNA2004800324109 A CN A2004800324109A CN 200480032410 A CN200480032410 A CN 200480032410A CN 1875380 A CN1875380 A CN 1875380A
- Authority
- CN
- China
- Prior art keywords
- organ
- mip
- ray
- outside surface
- rays
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
A system (100) and corresponding method (200) for organ rendering are provided, where the system (100) includes a processor (102), an imaging adapter (130) in signal communication with the processor for receiving organ scan data indicative of an organ, a segmentation unit (170) in signal communication with the processor for segmenting an outer surface of the organ, and a ray-casting unit (180) in signal communication with the processor for providing a maximum intensity projection (MIP) substantially normal to the segmented outer surface for each of a plurality of rays, and for forming a ray collection by casting each of the plurality of rays from a center of the organ through the respective MIP to the segmented outer surface; and where the method (200) includes segmenting (216) an outer surface of an organ, providing (218) a maximum intensity projection (MIP) substantially normal to the segmented outer surface for each of a plurality of rays, and forming a ray collection by casting (220) each of the plurality of rays from a center of the organ through the respective MIP to the segmented outer surface.
Description
The cross reference of related application
That the application requires to submit on November 3rd, 2003, title intactly is attached to it herein by reference for the rights and interests of the U.S. Provisional Application sequence number 60/516953 (attorney docket 2003P16622US) of " being used for crown visual high-quality normal direction MIP presents ".
Background
For example, the medical image scan-data obtains with sliced form in various types of imaging forms usually.These sections then are stacked to form three-dimensional (3D) volume.Hope is carried out visual to volume and is cut apart.
In the current method of medical image scanning, the researchist has developed the various cutting techniques that are used to isolate heart coronary artery.A large amount of patient that the research in this field is subjected to coronary artery disease promotes.Heart coronary artery is because its size and be difficult to usually cut apart near heart and blood bank surface.
Therefore, required is a kind ofly can be the crown visual system and method that presents.The disclosure is at these and other problem.
Summary of the invention
These and other defective of prior art and shortcoming solve by being used for the crown visual system and method that presents.
Being used for crown visual organ presents system and comprises: processor; The imaging adapter carries out the signal transmission with processor, is used to receive the organ scan-data of expression organ; Cutting unit carries out the signal transmission with processor, is used to cut apart the outside surface of organ; And ray projecting unit, carry out the signal transmission with processor, be used in a plurality of rays each that maximum intensity projection perpendicular to cutting apart of outside surface is provided in fact, and be used for cutting apart outside surface by each of a plurality of rays projected institute from the organ center via homolographic projection, form ray collection.
Being used for the correlation method that crown visual organ presents comprises: the outside surface of cutting apart organ; For in a plurality of rays each provides in fact maximum intensity projection perpendicular to cutting apart of outside surface; And cut apart outside surface and form ray collection by in a plurality of rays each projected institute from the center of organ via homolographic projection.
Read the explanation of following example embodiment in conjunction with the accompanying drawings, these and other aspect of the present disclosure, feature and advantage will become obvious.
Brief description
The disclosure is lectured according to following exemplary drawings and is used for the crown visual system and method that presents, and accompanying drawing comprises:
Fig. 1 explanation according to an illustrative embodiment of the present disclosure, be arranged for the synoptic diagram of the crown visual system that presents;
Fig. 2 explanation according to an illustrative embodiment of the present disclosure, be arranged for the process flow diagram of the crown visual method that presents;
Fig. 3 explanation is according to the figure of the computed segmentation of an illustrative embodiment of the present disclosure;
Fig. 4 explanation is according to the figure of the distance map of an illustrative embodiment of the present disclosure;
Fig. 5 explanation is according to the figure of an illustrative embodiment of the present disclosure, the maximum intensity projection (MIP) that shows on three-dimensional surface models such as (3D);
Fig. 6 explanation according to an illustrative embodiment of the present disclosure, be used to make up the synoptic diagram of ray collection;
Fig. 7 explanation according to an illustrative embodiment of the present disclosure, be used to analyze the figure of ray collection; And
Fig. 8 explanation according to an illustrative embodiment of the present disclosure, be used for figure that ray collection is filtered.
DETAILED DESCRIPTION OF THE PREFERRED
According to preferred embodiment of the present disclosure, this paper describes and is used for the crown visual system and method that presents.This method allows the user to carry out visual near heart coronaries the heart surface and vascular better.The very high-quality that this method is used for calculating heart coronaries presents.Achieve this end, ray projects the surface on the sheet (slab) of heart surface from the heart center.Then, calculate maximum intensity projection (MIP) for all rays.Gained presents and is called " normal direction MIP ", because all rays are orthogonal in fact or perpendicular to the general surface of heart.The disclosure is described the high-quality that makes up crown tree and is presented required step.
The researchist has developed various cutting techniques, is used for the isolation of heart coronary artery and visual.A large amount of patient that the research in this field is subjected to coronary artery disease promotes.Heart coronary artery is because its size and be difficult to visual usually near heart and blood bank surface.In addition, Pulmonic near having blocked the heart surface at coronary artery place.
Owing to cut apart the difficulty of existence, and since coronary artery near heart surface, the unfolded surface method can be used to solve visualization problem.This technology provides the visual improvement for heart coronary artery.Current disclosed method comprises the surface of " expansion " heart and creates this surperficial maximum intensity projection (MIP).The result is two dimension (2D) figure of heart surface, comprises peripheral vessels, and it is easier to read for the cardiologist.
As shown in Figure 1, the crown visual system that presents according to an illustrative embodiment of the present disclosure is generally represented by reference number 100.System 100 comprises at least one processor or the central processing unit (CPU) 102 that carries out the signal transmission with system bus 104.ROM (read-only memory) (ROM) 106, random-access memory (ram) 108, display adapter 110, I/O adapter 112, user interface adapter 114, communication adapter 128 and imaging adapter 130 also carry out the signal transmission with system bus 104.Display unit 116 carries out the signal transmission via display adapter 110 and system bus 104.Disk storage unit 118, for example magnetic or rom memory cell carry out the signal transmission via I/O adapter 112 and system bus 104.Mouse 120, keyboard 122 and eye tracking apparatus 124 carry out the signal transmission via user interface adapter 114 and system bus 104.MR imaging apparatus 132 is carried out the signal transmission via imaging adapter 130 and system bus 104.
Fig. 2 is generally represented by reference number 200 according to the crown visual process flow diagram that presents that is used for of an illustrative embodiment of the present disclosure.Flow process Figure 200 comprises the beginning frame 210 that passes control to functional block 212.Functional block 212 is initiated preliminary heart scanning session, and passes control to input frame 214.Input frame 214 receives preliminary heart scan data, and passes control to functional block 216.
Functional block 216 is cut apart the outside surface of heart, and passes control to functional block 218.Functional block 218 is carried out in fact maximum intensity projection (MIP) perpendicular to cutting apart of outside surface for a plurality of rays, and passes control to functional block 220.Functional block 220 is cut apart outside surface and is formed ray collection by each of ray projected institute from the center of organ by corresponding MIP, and passes control to functional block 222.Functional block 222 is launched the 3D model, and passes control to end block 224.
Refer now to Fig. 3, cut apart by reference number 300 and generally represent.Cut apart and cut apart 310 and cut apart 320 afterwards before comprising, the pattern cut algorithm that their employings have the single seed point in the heart obtains.
As shown in Figure 4, distance map is generally represented by reference number 400.Distance map 400 show ray (r) and with the inside and outside distance of center (c) (D1, D2).(x, y z) are distance from the voxel to the heart surface to each value on the 3D distance map 400.
Fig. 5 waits surface model generally to be represented by reference number 500.Here, the image that produces from ray projection directly texture on the face such as is shown as with 3D.
Refer now to Fig. 6, ray collection unit is generally represented by reference number 600.Ray collection unit 600 comprises ray processing unit 610, cuts apart/projecting unit 620 and dish I/O unit 630.Ray processing unit 610 comprises collection of rays part 640, collection of rays control section 642, maximum intensity projection (MIP) part 644, ray control section 646 and voxel part 648.Cut apart/projecting unit 620 comprises cutting unit 670 and ray projecting unit 680, and dish I/O unit 630 comprises loading unit 650 and preserve unit 660.Ray collection unit makes up ray collection by creating the data structure that comprises voxel for each ray.That data structure is used for again is fast mutual, as changing MIP thickness and being used for further vessel detection.In this exemplary internal data structure, current ray collection 640 for example can comprise 400 * 400 * 40 elements.
As shown in Figure 7, the ray collection analysis is generally represented by reference number 700.Analyzing 700 comprises simple venation pipe 710 and intersects 720.The distribution curve 712 of simple venation pipe 710 demonstrates single intensity peak 714, and 720 the distribution curve 722 of intersecting has then shown the dual intensity peak value 724 and 726 of first and second vasculars respectively.Ray collection by analyzing each ray distribution curve and search known pattern analysis.When can choose from launch texture by focus, the 3D position is known for each voxel on the distribution curve accurately.
Fig. 8, ray collection is filtered and is generally represented by reference number 800.Filter 800 and comprise radiation profile 810 and vascular tracing figure 820.By search the local maximum on each ray and analyze slope and intensity to obtain first group of point, come ray is filtered.The 3D point cloud chart of filtered points is provided.
In the operation of the method 200 of Fig. 2, example embodiment is used for following technology the isolation of heart surface.For computing method to MIP, the surface of heart is identified, and makes up the 3D distance map then.Those of ordinary skill in the related art are known, and 3D pattern cut algorithm can be used to the cutting cardiac data.The result of this isolation is the clear boundary between the remainder of definition cardiac muscle and volume.
Partitioning algorithm output 3D mask, it will be converted to the 3D distance map in next procedure.Notice that the quality of final image does not depend on result itself that cut apart on certain program, this makes this technology very healthy and strong.In addition, other method can be used for isolating heart, for example cutting apart based on model in step.
In cutting apart in 300 the operation of reference Fig. 3, the calculating of distance map can realize from the 3D mask.Distance map can be used to the normal vector that calculates the thickness of MIP and obtain lip-deep any point, as in the following Table 1 through pseudo-code shown in the table.
The initialization of // distance map
For each voxel
A for n=14 different directions (6 axial+8 diagonal angles)
Mask value on the Get direction n
If is on the edge
Distance=to the distance at edge;
The distance_map[voxel]=distance;
// the first traversal
For each voxel from start to end
Temp=distance_map[voxel];
for(n=1;n<14;n+=2)
Temp=min (temp, " incremental modulus "+distance_map[voxel+increment])
The distance_map[voxel]=temp;
}
// the second traversal (with reverse order)
For each voxel from finish to start
{
Temp=distance_map[voxel];
for(n=0;n<14;n+=2)
Temp=min (temp, " incremental modulus "+distance_map[voxel+increment])
distance_map[vlnd]=temp;
}
Table 1.
In the operation of the distance map 400 of reference Fig. 4, in case obtain distance map, then for each voxel projection radiation on surface.Each ray defines by vector and the predefine scope perpendicular to the surface.Voxel along the maximum intensity of ray shows on heart surface.Be generally each voxel on zero the heart surface for distance map value wherein, normal vector is by adopting central point (c) or calculating by the gradient of the distance map on the employing voxel coordinate.For each distance of increment of heart inside, ray coordinates is calculated as r and equals n1 and multiply by normal vector, and the intensity on each ray coordinates is stored.In addition, for each distance of increment of heart outside, ray coordinates is calculated as r and equals n2 and multiply by normal vector, and the intensity on each ray coordinates is stored.By reading the voxel of the maximum intensity between D1 and the D2, correlation line computation maximum intensity projection (MIP), and the MIP that shows current ray.Then, the gained image can be directly texture on the face such as be shown as with 3D.
Embodiment of the present disclosure provides crown isolation, wherein have the visual improvement quality of unfolded image and angiogram pattern, algorithm improvement speed, have blind spot that automatic crown tree is visually cut apart, vascular is checked improvement, reduce and the useful orientation that does not have distortion.
In the operation of example embodiment,, from launch view, obtain three-dimensional (3D) and isolate by cutting apart heart surface and ray being projected the surface from the center.Ray collection makes up by creating the structure that comprises voxel for each ray.That structure can be used for again is fast mutual, as changing MIP thickness and being used for further vessel detection.Ray collection by analyzing each ray distribution curve and search known pattern analysis, therein, the 3D position is known for each voxel on the distribution curve accurately.By obtain the local maximum on each ray and analyze slope and intensity so that first group of point to be provided, come ray is filtered.Concurrent through overregulating and handle, so as to approach better vascular center line, eliminate noise and increase coherence between the point on the same vessel.Adopt eigenvalue point of adjustment in 3D, so that help to determine the center of point, and a part of eliminating noise.
Focus becomes increases the density of the point of pipe along the pulse, so that make link easier.Adopt one group of heuristics that each point is realized the continuity filtrator, so that create other point on the same vessel.Heuristics can comprise that minor increment, the ridge between the existing point detects and/or the velocity of previous point.Click vessel detection is a little by the institute that has similar direction in the same vicinity that is linked at a seed point and is realized.
Those of ordinary skill in the related art will appreciate that alternative can detect crown two dimension (2D) center line, the separation artery launches the calcium in the view and/or separates the lung vascular with spot, the detection that vein, assessment launch view.
This demonstration methods is that the vascular on heart coronaries and the heart surface brings improved visualization technique.Feature is surface of " expansion " heart and the MIP that creates this unfolded surface.The gained 2D figure of heart surface comprises the high-contrast vascular.
Therefore, the exemplary teachings that is used to launch heart surface is finished by four steps: 1) outside surface of heart cuts apart; 2) the fixedly 3D model parameterization of heart; 3) but ray is projected the 3D mold surface of applying MIP filtrator from the center of heart; And 4) expansion of 3D model.
For example, pattern cut algorithm known in the art can be used for the outside surface of heart is cut apart.Create distance map from this result of cutting apart, so that the each point of assessment on volume is to the distance of the heart surface of cutting apart.Then, known 3D model can be fitted to heart, makes the surface fitting mold surface of heart.After this step, ray projects the 3D mold surface from the center of heart.Distribution curve is created, and simultaneously ray is by heart, and response filter is employed, so that detect the final position of vascular.If find the position, then algorithm is presented at the result on the 3D mold surface.Those of ordinary skill in the related art will appreciate that the expansion of 3D model is the problem of broad research, and can adopt some different algorithms.
A preferred embodiment adopts spheroid as the 3D model, and adopts maximum intensity projection (MIP) as profile curve filter.Though can adopt alternative 3D model and ray filter,, a part may influence result's quality unfriendly.
Therefore, preferred embodiment of the present disclosure is for crown isolation and visually provide powerful heart to launch instrument, makes the user can extract obvious characteristic and is subjected to region-of-interest.Preferred embodiment can be used as acquisition time modeling extremely useful in the clinical practice and automatic aftertreatment instrument.
May to be those of ordinary skill in the related art be easy to determine according to the theory of this paper these and other feature and advantage of the present disclosure.Everybody understands, and theory of the present disclosure can make up by various forms of hardware, software, firmware, application specific processor or its and realize.
Most preferably, theory of the present disclosure realizes with the array configuration of hardware and software.In addition, software is preferably realized with the tangible application program that is embodied on the program storage unit (PSU).Application program can upload on the machine that comprises any suitable architecture and by its execution.This machine is preferably realized on the computer platform with the hardware such as one or more central processing units (CPU), random-access memory (ram) and I/O (I/O) interface.Computer platform also can comprise operating system and micro-instruction code.Various process as herein described and function can be can be by ingredient or the ingredient of application program or their any combination of the micro-instruction code of CPU operation.In addition, other various peripheral cells, can be connected to computer platform as additional-data storage unit and print unit.
It is also understood that because a part of construction system assembly shown in the drawings and method preferably realize with software, so actual between system component or the process function frame connect may be according to the mode of disclosure programming and different.Provide the theory of this paper, those of ordinary skill in the related art can consider that of the present disclosure these reach similarly realization or configuration.
Though this paper describes illustrative embodiment with reference to the accompanying drawings, should be appreciated that the disclosure is not limited to those clear and definite embodiment, various changes and modification can be implemented by those of ordinary skill in the art, and do not deviate from the scope of the present disclosure or spirit.All these class changes and modification all are intended to be included within the scope of the present disclosure of setting forth as claims.
Claims (20)
1. organ rendering method comprises:
The outside surface of cutting apart organ;
For in a plurality of rays each provides in fact maximum intensity projection (MIP) perpendicular to cutting apart of outside surface; And
Cut apart outside surface by each center from organ in described a plurality of rays is projected institute by corresponding MIP, form ray collection.
2. the method for claim 1 is characterized in that, also comprises two dimension (2D) image that three-dimensional (3D) model expands into normal direction MIP is corresponding in fact with organ.
3. the method for claim 1 is characterized in that:
Described organ is a heart; And
The described MIP of normal direction in fact shows the feature that comprises coronary vasodilator.
4. the method for claim 1 is characterized in that, also comprises receiving preliminary organ scan-data.
5. method as claimed in claim 2 is characterized in that, described 3D model is fixed.
6. method as claimed in claim 2 is characterized in that, described 3D model is chosen from the group that comprises spheroid, cylinder and ellipsoidal 3D shape.
7. the method for claim 1 is characterized in that, projection radiation comprises maximum intensity projection's filter application at least one ray.
8. one kind is used for the device that organ presents, and comprising:
Partition member is used to cut apart the outside surface of organ;
Projection part is used in a plurality of rays each that maximum intensity projection (MIP) perpendicular to cutting apart of outside surface is provided in fact; And
Collecting part is used for projecting institute by each center from organ with described a plurality of rays by corresponding MIP and cuts apart outside surface, forms ray collection.
9. device as claimed in claim 8 is characterized in that, also comprises the sweep unit that is used to receive the organ scan-data.
10. device as claimed in claim 8 is characterized in that, also comprises the display unit that is used to show the organ scan-data.
11. one kind is used for the system that organ presents, comprises:
Processor;
The imaging adapter carries out the signal transmission with described processor, is used to receive the organ scan-data of expression organ;
Cutting unit carries out the signal transmission with described processor, is used to cut apart the outside surface of organ; And
The ray projecting unit, carry out the signal transmission with described processor, be used in a plurality of rays each that maximum intensity projection (MIP) perpendicular to cutting apart of outside surface is provided in fact, and be used for cutting apart outside surface and forming ray collection by each of described a plurality of rays projected institute from the center of organ by corresponding MIP.
12. system as claimed in claim 11 is characterized in that, described processor presents described organ scan-data as three-dimensional (3D) image relevant with being subjected to region-of-interest.
13. system as claimed in claim 12 is characterized in that, also comprises display adapter, it and described processor carry out the signal transmission, are used to show the 3D rendering that is presented.
14. system as claimed in claim 13 is characterized in that, also comprises the user interface adapter that is used to check quality of scanning.
15. system as claimed in claim 11 is characterized in that, described processor provides the distribution of each ray, and the identification known pattern.
16. a machine-readable program storage device visibly is embodied as the program of the executable instruction of machine, is used for the program step that organ presents so that carry out, described program step comprises:
The outside surface of cutting apart organ;
For in a plurality of rays each provides in fact maximum intensity projection (MIP) perpendicular to cutting apart of outside surface; And
Cut apart outside surface by each center from organ in described a plurality of rays is projected institute by corresponding MIP, form ray collection.
17. device as claimed in claim 16 is characterized in that, described program step also comprises two dimension (2D) image with the three-dimensional of organ (3D) model expands into normal direction MIP is corresponding in fact.
18. device as claimed in claim 16 is characterized in that:
Described organ is a heart; And
The described MIP of normal direction in fact shows the feature that comprises coronary vasodilator.
19. device as claimed in claim 16 is characterized in that, the described program step of projection radiation comprises maximum intensity projection's filter application at least one ray.
20. device as claimed in claim 16 is characterized in that, described program step also comprises:
The distribution of each ray is provided; And
In distributing, each discerns known pattern.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US51695303P | 2003-11-03 | 2003-11-03 | |
| US60/516,953 | 2003-11-03 | ||
| US10/978,654 | 2004-11-01 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN1875380A true CN1875380A (en) | 2006-12-06 |
Family
ID=37484942
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNA2004800324109A Pending CN1875380A (en) | 2003-11-03 | 2004-11-03 | Rendering for coronary visualization |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1875380A (en) |
-
2004
- 2004-11-03 CN CNA2004800324109A patent/CN1875380A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Leung et al. | Automated border detection in three-dimensional echocardiography: principles and promises | |
| US20070165917A1 (en) | Fully automatic vessel tree segmentation | |
| US8107703B2 (en) | Quantitative real-time 4D stress test analysis | |
| JP4918048B2 (en) | Image processing apparatus and method | |
| US7233329B2 (en) | Rendering for coronary visualization | |
| CN101596111B (en) | Automatic localization of the left ventricle in cardiac cine magnetic resonance imaging | |
| CN1640360A (en) | Automatic coronary isolation using a n-Mip (normal maximum intensity projection) ray casting technique | |
| US8948484B2 (en) | Method and system for automatic view planning for cardiac magnetic resonance imaging acquisition | |
| US7711160B2 (en) | Automatic optimal view determination for cardiac acquisitions | |
| RU2669680C2 (en) | View classification-based model initialisation | |
| US7542595B2 (en) | Virtual organ unfolding for visualization | |
| US8428316B2 (en) | Coronary reconstruction from rotational X-ray projection sequence | |
| CN1748646A (en) | Method and apparatus for medical intervention procedure planning and location and navigation of an intervention tool | |
| EP3432797B1 (en) | Systems and methods for characterizing a central axis of a bone from a 3d anatomical image | |
| CN108846830A (en) | The method, apparatus and storage medium be automatically positioned to lumbar vertebrae in CT | |
| CN1918601A (en) | Apparatus and method for registering images of a structured object | |
| CN1862596A (en) | System and method for fused PET-CT visualization for heart unfolding | |
| US7233330B2 (en) | Organ wall analysis with ray-casting | |
| CN110335235A (en) | Processing unit, processing system and the medium of cardiologic medical image | |
| CN1620991A (en) | Cardiac display methods and apparatus | |
| US20220222812A1 (en) | Device and method for pneumonia detection based on deep learning | |
| CN114445429A (en) | Whole-heart ct segmentation method and device based on multiple labels and multiple decoders | |
| CN111429457B (en) | Intelligent evaluation method, device, equipment and medium for brightness of local area of image | |
| CN1531904A (en) | Method of operating imaging medical detection device | |
| CN1875380A (en) | Rendering for coronary visualization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |