CN1862596A - System and method for fused PET-CT visualization for heart unfolding - Google Patents
System and method for fused PET-CT visualization for heart unfolding Download PDFInfo
- Publication number
- CN1862596A CN1862596A CN 200610084045 CN200610084045A CN1862596A CN 1862596 A CN1862596 A CN 1862596A CN 200610084045 CN200610084045 CN 200610084045 CN 200610084045 A CN200610084045 A CN 200610084045A CN 1862596 A CN1862596 A CN 1862596A
- Authority
- CN
- China
- Prior art keywords
- image
- heart
- projection
- pet
- ellipsoid
- 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
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Nuclear Medicine (AREA)
Abstract
A method for fusing two digitized image datasets for 2-dimensional visualization of the heart, includes providing a first digitized image and a second digitized image, each image acquired from a different imaging modality and each comprising a plurality of intensities corresponding to a domain of points on a 3-dimensional grid, and wherein the first and second images have been registered, segmenting the heart in one of said digitized images, forming a 2-dimensional maximum intensity projection of a surface of said heart from each of said digitized images, and fusing said images by overlaying one 2-dimensional projection of said surface over the other 2-dimensional projection of said surface.
Description
Technical field
The present invention relates to visual (visualization) of the image data amount in imaging of medical and computer vision application.
Background technology
The name that the application requires people such as Moreau-Gobard to submit on April 19th, 2005 is called the U.S. Provisional Application No.60/672 of " FusedPET-CT Visualization for Heart Unfolding (visual at the PET-CT that is merged that heart launches) ", 679 right of priority, the content of this application is hereby incorporated by.
The appearance that combines the multi-modal imaging scanner of positron emission computerized tomography (PET) and computed tomography (CT) causes the new example of image demonstration and performance, and this example produces new challenge at workstation interpretation software, image-guidance with in communicating by letter.The essence of multi-modal imaging is to superpose from the ability of the image-forming information of different mode in compulsory mode visually.This is by utilizing coloud coding (color-encoding) technology that function realized that with the synthetic multidimensional view of anatomical data set this coloud coding technology provides visual prompting on the space distribution of view data.Multi-mode is visualized as not only needs structured data but also need the diagnosis of performance data that extra dimension is provided.
Need special reconstruction software and image display program to come (native) image of nature is readjusted plane into the reformatting of quadrature or inclination from different spatial resolutions, in these planes, to carrying out coloud coding from the data of PET image and these data being overlapped on the corresponding anatomy CT image.The color addition technology allows the user visually to identify the zone of high tracer activity (tracer activity) and determines the anatomical structure on basis.Because the multidimensional property of data, the multidimensional airmanship that visual needs are mutual, this airmanship allow the observer on three direction in space mobile visualization planes.And CT provides anatomic information (such as location coronarius), and the PET scanning of heart can be used to determine that blood flows to cardiac muscle and the symptom that helps to estimate coronary artery disease.Present challenge is to use the CT and the PET volume that are merged to create " fusion " curved surface MPR image and more senior 2D projection, such as the U.S. Patent Application Publication No. of submitting on August 3rd, 2004 at RomainMoreau-Gobard is that 2005/0058328 name is called the disclosed technology of co-pending application of " Heart unfolding for coronary visualization (being used for crown visual heart launches) ", and the content of this application is hereby incorporated by.
Visual people's the reciprocation that needs usually of the 3D of cardiac image is prepared data and is removed the structure (such as lung, liver, rib etc.) that does not belong to heart.For the work efficiency that improves the cardiologist with improve diagnostic accuracy, the unfolded surface at the CT pattern has been described particularly in above-mentioned patented claim, so that improve visual.
Summary of the invention
Exemplary embodiment of the present invention as the described herein generally includes the method and system that is used at the heart expansion PET being become with the CT image co-registration image of single reformatting.This technology adopts volumetric layer to be used for the surface of visual heart and provides a kind of and the interactional new method of multi-mode 3D volume.According to embodiments of the invention, cut apart and " expansions " this surface after, flat (2D) image that synthesizes of the heart surface that obtains from the CT volume.The PET volume is repeated same operation, and by becoming the color M PR image of single reformatting to obtain final image with the PET information fusion CT.Final overlapping image helps the doctor to launch expression in conjunction with heart to watch the diagnostic message that is obtained from PET scanning, and this diagnostic message is relevant with anatomic information from CT scan.
According to aspects of the present invention, providing a kind of is used at the two-dimensional visualization of heart and merges the method for two digitizing image data sets, this method comprises provides first digitized image and second digitized image, each image of collection and each image comprise a plurality of intensity corresponding to the territory of the point on the 3D grid from different imaging patterns, and wherein first and second images are registered, cut apart heart in the digitized image in described digitized image, form the two-dimentional maximum intensity projection on the surface of described heart according to each described digitized image, and merge described image in another two-dimensional projection by described surface that a two-dimensional projection on described surface is added to.
According to another aspect of the present invention, from positron emission computerized tomography (PET) scanning, gather an image in the described image, and from computed tomography (CT) scanning, gather another image in described two images.
According to another aspect of the present invention, in the CT image, cut apart heart.
According to another aspect of the present invention, cardiac segmentation identifies the surface of heart.
According to another aspect of the present invention, described method comprise that formation will be from described heart each put the distance of heart surface and the distance map of described spot correlation connection.
According to another aspect of the present invention, described method comprises and uses described distance map by triaxial ellipsoid heart surface to carry out modeling.
According to another aspect of the present invention, triaxial ellipsoid is the twin shaft ellipsoid, and wherein the center is calculated as the center of gravity G that is had a few that is positioned on distance 0 the distance map contour surface, and major axis Z is a vector
Wherein M makes length
Point on the maximum described contour surface, minor axis X is a vector
Wherein N makes length
Point on the minimum described contour surface, and the 3rd Y and described Z and X quadrature.
According to another aspect of the present invention, forming described two-dimentional maximum intensity projection comprises ray is projected the described heart so that form described two-dimensional projection from described ellipsoid, (u, intensity v) is each intensity I along ray to some I (X)=I in the wherein said 2D projection
nNormalized and.
Wherein said 2D point X (u, v) calculate according to following formula:
Wherein ∈ [0,2 π] and θ ∈ [0, π] are respectively lip-deep each some P (, longitude θ) and latitudes of described ellipsoid.
According to another aspect of the present invention, described method comprises uses horizontal Mercator (Mercator) projection, is used for forming according to described ellipsoid the described two-dimentional maximum intensity projection of described heart surface.
According to another aspect of the present invention, the described two-dimentional maximum intensity projection that utilizes described PET image to form described heart surface comprises and utilizes described registration to calculate respective point on the described heart surface in the described PET image from described CT image, in described PET image, select together with inner and at heart slightly slightly at lip-deep a plurality of points of the point of heart outside, and to the intensity level summation of the ray that passes described selected a plurality of points along projection and ask average.
According to another aspect of the present invention, provide a kind of program storage device that can read by computing machine, this program storage device visibly comprises the executable programmed instruction of computing machine, so that carry out the described method step be used for merging at the two-dimensional visualization of heart two digitizing image data sets.
Description of drawings
Fig. 1 according to one exemplary embodiment of the present invention described point from 3d space to the 2D space the exemplary map of point.
Fig. 2 is used for launching and being the process flow diagram of method of the image of single reformatting with PET and CT image co-registration at heart according to an embodiment of the invention.
Fig. 3 has described example view with the PET volume of CT volume registration according to embodiments of the invention.
Fig. 4 has described the synoptic diagram of the PET scanning of diseased heart according to embodiments of the invention.
Fig. 5 has described the example of the n-MIP view that is unfolded of heart surface according to embodiments of the invention.
Fig. 6 is used to realize to be used for launch and be the block diagram of exemplary computer system of method of the image of single reformatting with PET and CT image co-registration at heart according to an embodiment of the invention.
Embodiment
Exemplary embodiment of the present invention as the described herein generally includes the visual system and method for the PET-CT that is merged that is used at the heart expansion.But the form that the present invention can multiple replacement is implemented should not be constructed to be limited to the embodiments set forth herein.
In view of the above, although the present invention is subject to the influence of the various modifications and the form of replacement, its specific embodiment is shown in the accompanying drawings by example, and will be described in detail here.But, should be understood that intention does not limit the invention to particular forms disclosed, but opposite, the present invention will be contained all modifications, equivalents and the replaceable form that falls in the spirit and scope of the present invention.Run through the description of the drawings, same numbering is represented same element.
Should be understood that,, should not limit these elements by these terms though use various elements such as first, second term description such as grade at this.These terms only are used to make an element and another element to differentiate.For example, first element can be called as second element, and similarly, second element can be called as first element, does not depart from the scope of the present invention simultaneously.As used in this, term " and/or " comprise the one or more listed projects that are associated arbitrarily and all combinations.
Should be understood that when an element being called " being connected " or " being coupled " to another element, this element can directly be connected or be coupled to other element maybe may exist insertion element.On the contrary, when an element being called " directly being connected " or " directly being coupled ", there is not insertion element to another element.Other speech that is used for describing the relation between the element should be explained (that is, " between " and " directly between ", " vicinity " and " directly contiguous " etc.) in the same way.
At this employed term only is not mean that it is limitation of the present invention in order to describe certain embodiments.As used herein, " " of singulative, " one " and " this " mean and also comprise plural form, unless context is clearly pointed out in addition.What it should be further apparent that is, when when this uses, term " comprises ", " having comprised ", " comprising " and/or the existence of stipulating " having comprised " described feature, integral body, step, operation, element and/or parts, does not exist or increases one or more further features, integral body, step, operation, element, parts and/or its group but do not get rid of.
As used in this, term " image " refers to the multidimensional data by discontinuous pictorial element (for example, the voxel of the pixel of 2D image and 3D rendering) formation.This image can for example be by computed tomography, magnetic resonance imaging, ultrasound wave or other well known to a person skilled in the art the person under inspection's that medical imaging system is collected medical image arbitrarily.This image also can provide from the non-medical environment, this non-medical environment such as remote sensing system, electron microscopy etc.Though image can be considered to from R
3To the function of R, but method of the present invention is not limited to these images, and can be used to the image of any dimension, for example 2D picture or 3D volume.For two dimension or 3-D view, the territory of image is typically two dimension or three-dimensional rectangle array, and wherein each pixel or voxel can come addressing with reference to the mutually orthogonal axle of a group 2 or 3.Term " numeral " and " digitized " refer to that (if suitably) is by digital collection system or image numeral or digitized format or volume by gathering from the conversion of analog image as used in this.
Along with adopting multi-mode PET/CT imaging technique more and more in routine clinical work, new merged visual on the surface of launching at heart had any different for accurate diagnosis of cardiovascular diseases.The 2D view that is unfolded of heart just demonstrates all coronary arteries and heart external perihaemal canal quickly subtly.The clinician is concerned about this view especially, because they can study the structure and the internal anatomy of blood vessel, and does not need and reciprocations such as 3D volume, transport function, shear plane.Use this view, the clinician can be apace and is located exactly and quantize vascular pathological.
Coronary artery is responsible for supplying with oxygenated blood (oxygenated blood) to cardiac muscle.If a disease has changed normal function coronarius, then the possibility of result is the dead of some tissue and is fatal.On the other hand, some vascular pathologicals can fully block coronary artery and to heart body without any influence.In order to recognize that those need the pathology of paying close attention to immediately, the clinician need use the PET imaging and need connect PET result and CT imaging.PET shows the chemical functionality or the death of organ and tissue.Along with PET and CT are superimposed upon on the view of heart expansion, the clinician can measure the seriousness of crown pathology immediately, and according to the CT image specific pathology and its influence to heart body is connected, as determined according to the PET image.
The process flow diagram that is used for merging for the surface of launching heart the illustrative methods of PET and CT data set has been shown among Fig. 2, and this method comprises the following steps.In step 21, gather PET and CT cardiac image data amount.According to embodiments of the invention, use the PET-CT scanner to come images acquired, and registration information also can be comprised with view data.In step 22, one of split image volume, and identify heart surface.Utilize registration transformation that be mapped on other image cutting apart of an image.According to another embodiment of the invention, if do not provide registration information to image, then heart surface and other image volume registration.
In step 23, come heart surface is carried out modeling by the 3D ellipsoid.In step 24, projection radiation from the 3D ellipsoid is so that calculate the 2D view of heart surface.In step 25, these rays are mapped to 2D maximum intensity projection (MIP), so that form the view that is unfolded.Be noted that step 23,24 and 25 can or carry out at the PET data set or at the CT data set.Can preserve the 3D coordinate of each ray get off and utilize two registration transformations between the set to be applied to other image.At last, in step 26, the view stack of 2D PET and 2D CT expansion, so that form the 2D view that is merged.Though can use interchangeable 3D model and expansion technique, influence result's quality on the contrary.
Should be noted that, though be the fusion that comes descriptor data set according to PET and CT data set, but among other embodiment of the present invention, be applicable to the combination of the image data set that from any imaging pattern combination (such as SPECT-CT and even MR-CT), collects in this several different methods.
Next step is to cut apart heart volume.According to embodiments of the invention, since specific, be more preferably on the CT volume and work, and the registration correspondence between the voxel in use CT volume and the PET volume, cut apart to carry out corresponding PET data set.Main cause is the low spatial resolution of PET data, and this has reduced the precision of cutting apart, because the PET volume collected on a plurality of cardiac cycles, this causes " spot (blob) " volume of bluring.The pattern cut algorithm for example can be used to cut apart the outside surface of CT heart as known in the art.Thisly cut apart the 3D mask that generation wherein is labeled as pixel object and background.
In order to merge PET data set and CT data set, the coordinate system of these two data sets needs coupling, and common coordinate system is gathered because these two kinds of data sets are not total use.Therefore need to carry out the 3D conversion, so that these two data sets of registration.Have a lot of explanations known in the art to be used for the method for the 3D conversion of two images of registration, and any suitable mapping algorithm can be used to registration PET and CT data set.Because PET-CT volume registration Algorithm is more and more effectively and in the widespread use of clinical boundary, thus the heart expansion technique that is used for CT now and PET use simultaneously.Explanatory view with the PET volume of CT volume registration has been shown among Fig. 3.The PET volume usually comprises registration data, because same equipment is normally used for not only producing the CT volume but also produce the PET volume, does not need the step of registration of adding in this case.But the better CT scanner that is used for the heart coronaries diagnosing image does not have the PET ability usually, and therefore needs registration.There are a lot of possible method for registering, are used to aim at two 3D " spot ".According to embodiments of the invention, by the PET volume being got threshold value and this volume being aimed at the registration of carrying out heart surface with the former CT heart surface of cutting apart.
According to this result of cutting apart, produced distance map, so that assessment is put to the distance of the heart surface of being cut apart from the volume each.Distance map from the surface of 3D mask can provide valuable constraint in the process of restoring external perihaemal canal.For example, because external perihaemal canal and heart surface keeping parallelism, so the distance of the two in the distance map (that is the distance on the dirty surface of these blood vessel centroids) variation reposefully.Distance map is the volume that the distance between each point and the object of reference (in this case, being the surface of heart) is associated with this point basically.The voxel of being cut apart carries out modeling by simple geometric body (such as ellipsoid) more to heart, so that easily use drawing practice that its surface is flattened, as described in greater detail.
Distance map is used to by ellipsoid or dual-axis rotation ellipsoid heart wall be carried out modeling.Though can consider more accurately model, such as triaxial ellipsoid, major defect is not to the simple mathematical solution of 3D-2D projection.
Twin shaft liptical projection is known drawing technique, but just is used to the imaging of medical of heart recently.If having reduced, the dual-axis rotation ellipsoid uses the contingent distortion of ellipsoid of revolution.For semiaxis length is the twin shaft ellipsoid of a and b, and the surface equation formula is:
Perhaps, with the form of parameter: x=a cos () sin (θ), y=a sin () sin (θ), z=b cos (θ), wherein ∈ [0,2 π] and θ ∈ [0, π].With employed the same in the common drawing, and θ are called as longitude and latitude respectively.The ellipsoid center is calculated as the center of gravity G that is had a few that is positioned on the distance map contour surface 0 (cardiac contour).Major axis Z is a vector
Wherein M makes length
Point on the maximum contour surface.Equally, minor axis X is a vector
Wherein N makes length
Point on the minimum contour surface.Derivation shaft Y has straight orthogonal basis, and B=(G, X, Y, Z).
Should be noted in the discussion above that the quality of modeling and do not rely on and cut apart quality (this makes that the selection of method as described herein and above-mentioned dividing method is irrelevant) and stable for noise.In addition, different with other organ, the natural shape of heart with ellipsoid.
Then, ray is from ellipsoid surface projection, calculating two dimension (" the 2D ") view of heart surface, and the intensity of collecting the voxel that ray crosses over.Lip-deep every the ray of ellipsoid is corresponding to the pixel on the 2D view.Therefore, draw out the 2D view by utilizing maximum intensity projection (MIP) to reproduce every ray.Collect voxel intensities in predetermined direction inside and outside heart wall and the scope.By following simple conversion, the node that calculates on the ellipsoid that is used for projection radiation distributes.With reference now to Fig. 1,, at each the some P in (on the ellipsoid) 3d space (, θ), according to following formula calculate (on the plane) 2D point X (u, v) 110:
(u v) is each intensity I of ray to each pixel I (X)=I on the 2D view
nMark one change and
Should be noted in the discussion above that if ray is directly throwed from distance map itself then ray may evenly can not utilize inadequately.In fact, because heart surface is irregular, so this ray can be omitted a lot of blood vessels.Another kind of solution can be based on the horizontal Mercator projection of twin shaft ellipsoid.The shortcoming of this method is that this volume must be the ellipsoid of rotation, this means that minor semi-axis is positioned on the pole-face, and major semi-axis is positioned on the equatorial plane.
These utilize ellipsoid that heart surface is carried out modeling, do not need the PET data set is carried out from this surface projection radiation and step that ray is mapped on the 2D MIP.The 3D coordinate that projects all rays on the CT volume is preserved.Utilize CT (X, Y, Z) and PET (the corresponding formula of the registration between Z) can use the CT ray coordinates to obtain intensity on the PET volume for X, Y, and needn't carry out identical complex operations set (cut apart, expansion, ray projection).
Usually, the PET volume has by low-down intensity level formed " hole ", and ill tissue is arranged in this " hole ".Described the synoptic diagram of the PET scanning of ill heart in Fig. 4, the low intensity value 41 of its cardiac 40 is centered on by the intensity values 42 of the heart tissue of health.According to embodiments of the invention, with the lip-deep intensity level of the corresponding 3D PET spot in CT surface in company with being considered in the value outside heart within the heart and slightly slightly.
Owing to kept about being projected onto CT volume { X
Ct, Y
Ct, Z
Ct, intensity (Intensity)
CtOn the information of ray, so can calculate registration { X
Pet, Y
Pet, Z
Pet, so that the intensity on the acquisition PET volume.Then intensity is asked on average to produce single value.Should be noted in the discussion above that PET is not an imaging pattern very accurately, therefore have thick registration and thick intensity level is an acceptable at the surface.Usually, the position of ill tissue is quite obvious on the PET image volume.
At last, in this stage, utilize the specified correspondence of registration of two 3D rendering volumes to superpose from the 2D expansion view of PET data set and CT data set.
Figure 5 illustrates the example of the n-MIP that is unfolded (that is, the standard maximum intensity projection) view of heart, this view produces by merging PET and CT view.In the flat surface of 2D this up-to-date generation, heart, can be clear that part and the unsound part (52,53) of supplying with myocardium oxygenated blood (51).
Should be understood that the present invention can various forms of hardware, software, firmware, dedicated process or its make up and realize.In one embodiment, the present invention can be implemented as the application program that visibly is included on the computer-readable program storage device by software.This application program can be uploaded and be carried out by the machine that comprises any appropriate configuration.
Fig. 6 is used to according to an embodiment of the invention carry out that to be used for launching PET and CT image co-registration at heart be the block diagram of exemplary computer system of method of the image of single reformatting.Referring now to Fig. 6, be used to realize that computer system 61 of the present invention especially comprises CPU (central processing unit) (CPU) 62, storer 63 and I/O (I/O) interface 64.Computer system 61 is coupled to display 65 and various input equipment 66 (such as mouse and keyboard) by I/O interface 64 usually.Auxiliary circuit can comprise the circuit such as cache memory, power supply, clock circuit and communication bus.Storer 63 can comprise random-access memory (ram), ROM (read-only memory) (ROM), disc driver, tape drive etc. or its combination.The present invention can be implemented as routine 67, and this routine 67 is stored in the storer 63 and by CPU62 to be carried out, so that handle the signal from signal source 68.Equally, computer system 61 is general-purpose computing systems, and when carrying out routine 67 of the present invention, this system becomes dedicated computer system.
Computer system 61 also comprises operating system and micro-instruction code.The part (or its combination) of the part of the micro-instruction code that various process described herein and function can or be carried out by operating system or the application program carried out by operating system.In addition, various other peripherals also can be connected to such as the additional data storage device and the computer platform of printing device.
Should further be understood that because some that described are in the accompanying drawings formed system units and method step can be realized by software, so the actual connection between the system unit (or process steps) can be depending on the present invention in its mode of programming difference.Be given in the instruction of the present invention that this provides, those of ordinary skill in the related art can design these and similar embodiment or configuration of the present invention.
Although with reference to preferred embodiment the present invention is described in detail, those skilled in the art will appreciate that, can carry out various modifications and replacement to it, and not break away from the spirit and scope of the present invention of being set forth as in the appended claim.
Claims (24)
1. one kind at the two-dimensional visualization of heart and merge the method for two digitizing image data sets, and it may further comprise the steps:
First digitized image and second digitized image are provided, each image be from different imaging patterns, gather and each image comprise a plurality of intensity corresponding to the territory of the point on the 3D grid, and wherein first and second images are registered;
Cut apart heart in the image in described digitized image;
Form the two-dimentional maximum intensity projection on the surface of described heart according to each described digitized image; With
Merge described image in another two-dimensional projection by described surface that a two-dimensional projection on described surface is added to.
2. the method for claim 1, wherein from positron emission computerized tomography (PET) scanning, gather an image in the described image, and from computed tomography (CT) scanning, gather another image in described two images.
3. method as claimed in claim 2 wherein, is cut apart heart in the CT image.
4. the method for claim 1, wherein described cardiac segmentation identifies the surface of described heart.
5. method as claimed in claim 4, each that comprises further that formation will be from described heart put the distance of heart surface and the distance map of described spot correlation connection.
6. method as claimed in claim 5 further comprises and uses described distance map by triaxial ellipsoid heart surface to carry out modeling.
7. method as claimed in claim 6, wherein, described triaxial ellipsoid is the twin shaft ellipsoid, and wherein, the center is calculated as the center of gravity G that is had a few that is positioned on distance 0 the distance map contour surface, and major axis Z is a vector
, wherein M makes length
Point on the maximum described contour surface, minor axis X is a vector
, wherein N makes length
Point on the minimum described contour surface, and the 3rd Y and described Z and X quadrature.
8. method as claimed in claim 6, wherein, forming described two-dimentional maximum intensity projection comprises ray is projected from described ellipsoid surface the described heart so that form described two-dimensional projection, wherein, (u, intensity v) is each intensity I along ray to some I (X)=I in the described 2D projection
nNormalized and
Wherein said 2D point X (u, v) calculate according to following formula:
Wherein ∈ [0,2 π] and θ ∈ [0, π] are respectively lip-deep each some P (, longitude θ) and latitudes of described ellipsoid.
9. method as claimed in claim 6 comprises and uses horizontal Mercator projection, is used for forming according to described ellipsoid the described two-dimentional maximum intensity projection of described heart surface.
10. method as claimed in claim 3, wherein, the described two-dimentional maximum intensity projection that utilizes described PET image to form described heart surface comprises and utilizes described registration to calculate respective point on the described heart surface in the described PET image from described CT image, in described PET image,, and the intensity level that passes the ray of described selected a plurality of points along projection sued for peace and ask average together with inner and select lip-deep a plurality of point at the point of heart outside slightly at heart slightly.
11. one kind at the two-dimensional visualization of heart and merge the method for two digitizing image data sets, it may further comprise the steps:
Digitizing computed tomography (CT) image and digitizing positron emission computerized tomography (PET) image are provided, and each image comprises a plurality of intensity corresponding to the territory of the point on the 3D grid;
In described CT image, cut apart heart surface;
Utilize triaxial ellipsoid that heart surface is carried out modeling and form the two-dimensional projection of described heart surface according to each described digitized image; With
Two-dimensional projection from described CT image is added to from the two-dimensional projection of described PET image.
12. method as claimed in claim 11, further comprise with the described heart surface of cutting apart from CT image registration to the PET image, wherein obtain corresponding formula from the CT image to the PET image.
13. method as claimed in claim 12, wherein, the two-dimensional projection that forms described heart surface according to described CT image comprises ray is projected the described heart so that form the maximum intensity projection of described heart surface from described ellipsoid that the coordinate of wherein said ray is saved.
14. method as claimed in claim 13, wherein, the two-dimensional projection that forms described heart surface according to described PET image comprises from described CT image and to described PET image described corresponding formula is applied to the described ray coordinates of preserving.
15. a computer-readable program storage device visibly comprises the executable instruction repertorie of computing machine, so that carry out the described method step be used for merging at the two-dimensional visualization of heart two digitizing image data sets, this method may further comprise the steps:
First digitized image and second digitized image are provided, each image be from different imaging patterns, gather and each image comprise a plurality of intensity corresponding to the territory of the point on the 3D grid, and wherein first and second images are registered;
Cut apart heart in the digitized image in described digitized image;
Form the two-dimentional maximum intensity projection on the surface of described heart according to each described digitized image; And
Merge described image in another two-dimensional projection by described surface that a two-dimensional projection on described surface is added to.
16. computer-readable program storage device as claimed in claim 16, wherein, an image from the described image of positron emission computerized tomography (PET) scanning collection, and another image from described two images of computed tomography (CT) scanning collection.
17. computer-readable program storage device as claimed in claim 16 wherein, is cut apart heart in the CT image.
18. computer-readable program storage device as claimed in claim 15, wherein, cardiac segmentation identifies the surface of heart.
19. computer-readable program storage device as claimed in claim 18, described method comprise further that formation will be from described heart each put the distance of heart surface and the distance map of described spot correlation connection.
20. further comprising, computer-readable program storage device as claimed in claim 19, described method use described distance map to come heart surface to be carried out modeling by triaxial ellipsoid.
21. computer-readable program storage device as claimed in claim 20, wherein, described triaxial ellipsoid is the twin shaft ellipsoid, and wherein, the center is calculated as the center of gravity G that is had a few that is positioned on distance 0 the distance map contour surface, and major axis Z is a vector
, wherein M makes length
Point on the maximum described contour surface, minor axis X is a vector
, wherein N makes length
Point on the minimum described contour surface, and the 3rd Y and described Z and X quadrature.
22. computer-readable program storage device as claimed in claim 20, wherein, forming described two-dimentional maximum intensity projection comprises ray is projected the described heart from described ellipsoid surface so that form described two-dimensional projection, wherein, (u, intensity v) is each intensity I along ray to some I (X)=I in the described 2D projection
nNormalized and
Wherein said 2D point X (u, v) calculate according to following formula:
Wherein ∈ [0,2 π] and θ ∈ [0, π] are respectively lip-deep each some P (, longitude θ) and the latitudes at described ellipsoid.
23. computer-readable program storage device as claimed in claim 20, described method further comprise the described two-dimentional maximum intensity projection that uses horizontal Mercator projection to be used for forming according to described ellipsoid described heart surface.
24. computer-readable program storage device as claimed in claim 17, wherein, the described two-dimentional maximum intensity projection that utilizes described PET image to form described heart surface comprises and utilizes described registration to calculate respective point on the described heart surface in the described PET image from described CT image, in described PET image,, and the intensity level that passes the ray of described selected a plurality of points along projection is sued for peace and asked average together with inner and select lip-deep a plurality of point at the point of heart outside slightly at heart slightly.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US67267905P | 2005-04-19 | 2005-04-19 | |
| US60/672679 | 2005-04-19 | ||
| US11/401640 | 2006-04-11 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN1862596A true CN1862596A (en) | 2006-11-15 |
Family
ID=37390018
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 200610084045 Pending CN1862596A (en) | 2005-04-19 | 2006-04-19 | System and method for fused PET-CT visualization for heart unfolding |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1862596A (en) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101856237A (en) * | 2009-02-10 | 2010-10-13 | 株式会社东芝 | Medical diagnostic imaging apparatus and method for correcting image |
| CN102395999A (en) * | 2009-04-15 | 2012-03-28 | 皇家飞利浦电子股份有限公司 | Quantification of medical image data |
| CN102397104A (en) * | 2010-07-29 | 2012-04-04 | 西门子公司 | Method for visualizing an atrium of the heart in a patient |
| CN102663824A (en) * | 2010-11-29 | 2012-09-12 | 西门子公司 | Method and system for heart isolation in cardiac computed tomography volumes |
| CN101903909B (en) * | 2007-12-18 | 2013-05-29 | 皇家飞利浦电子股份有限公司 | System for multimodality fusion of imaging data based on statistical models of anatomy |
| CN103345746A (en) * | 2013-06-25 | 2013-10-09 | 上海交通大学 | Medical image diagnostic method based on CT-PET |
| CN103687541A (en) * | 2011-03-02 | 2014-03-26 | 皇家飞利浦有限公司 | Visualization for navigation guidance |
| CN103826539A (en) * | 2012-09-20 | 2014-05-28 | 株式会社东芝 | Image processing system, X-ray diagnostic device, and image processing method |
| RU2524302C2 (en) * | 2009-02-17 | 2014-07-27 | Конинклейке Филипс Электроникс Н.В. | Extension on basis of model of vision field in radionuclide visualisation |
| CN105190692A (en) * | 2013-05-09 | 2015-12-23 | 皇家飞利浦有限公司 | Method and system for mesh segmentation and mesh registration |
| CN101536913B (en) * | 2008-03-21 | 2016-12-14 | 通用电气公司 | The method and apparatus of correcting multi-modality imaging data |
| CN107610095A (en) * | 2017-08-04 | 2018-01-19 | 南京邮电大学 | Heart CT coronary artery full-automatic partition methods based on image co-registration |
| CN111882513A (en) * | 2012-02-10 | 2020-11-03 | 皇家飞利浦有限公司 | Clinically driven image fusion |
-
2006
- 2006-04-19 CN CN 200610084045 patent/CN1862596A/en active Pending
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101903909B (en) * | 2007-12-18 | 2013-05-29 | 皇家飞利浦电子股份有限公司 | System for multimodality fusion of imaging data based on statistical models of anatomy |
| CN101536913B (en) * | 2008-03-21 | 2016-12-14 | 通用电气公司 | The method and apparatus of correcting multi-modality imaging data |
| CN101856237A (en) * | 2009-02-10 | 2010-10-13 | 株式会社东芝 | Medical diagnostic imaging apparatus and method for correcting image |
| RU2524302C2 (en) * | 2009-02-17 | 2014-07-27 | Конинклейке Филипс Электроникс Н.В. | Extension on basis of model of vision field in radionuclide visualisation |
| CN102395999A (en) * | 2009-04-15 | 2012-03-28 | 皇家飞利浦电子股份有限公司 | Quantification of medical image data |
| CN102397104A (en) * | 2010-07-29 | 2012-04-04 | 西门子公司 | Method for visualizing an atrium of the heart in a patient |
| CN102663824B (en) * | 2010-11-29 | 2016-06-29 | 西门子公司 | The method and system that heart in cardiac computer fault imaging volume separates |
| CN102663824A (en) * | 2010-11-29 | 2012-09-12 | 西门子公司 | Method and system for heart isolation in cardiac computed tomography volumes |
| CN103687541A (en) * | 2011-03-02 | 2014-03-26 | 皇家飞利浦有限公司 | Visualization for navigation guidance |
| CN111882513A (en) * | 2012-02-10 | 2020-11-03 | 皇家飞利浦有限公司 | Clinically driven image fusion |
| CN103826539B (en) * | 2012-09-20 | 2016-11-02 | 东芝医疗系统株式会社 | Image processing system, radiographic apparatus and image processing method |
| CN103826539A (en) * | 2012-09-20 | 2014-05-28 | 株式会社东芝 | Image processing system, X-ray diagnostic device, and image processing method |
| CN105190692B (en) * | 2013-05-09 | 2019-04-16 | 皇家飞利浦有限公司 | Method and system for mesh segmentation and mesh registration |
| CN105190692A (en) * | 2013-05-09 | 2015-12-23 | 皇家飞利浦有限公司 | Method and system for mesh segmentation and mesh registration |
| CN103345746A (en) * | 2013-06-25 | 2013-10-09 | 上海交通大学 | Medical image diagnostic method based on CT-PET |
| CN103345746B (en) * | 2013-06-25 | 2016-12-28 | 上海交通大学 | The method that CT picture reconstruction is gone out 3-D graphic |
| CN107610095A (en) * | 2017-08-04 | 2018-01-19 | 南京邮电大学 | Heart CT coronary artery full-automatic partition methods based on image co-registration |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1862596A (en) | System and method for fused PET-CT visualization for heart unfolding | |
| US7813535B2 (en) | System and method for fused PET-CT visualization for heart unfolding | |
| JP5632913B2 (en) | Multi-modality chest imaging system, method and program | |
| CN101796544B (en) | Visualization method and system of voxel data | |
| Kurgan et al. | Knowledge discovery approach to automated cardiac SPECT diagnosis | |
| US9020235B2 (en) | Systems and methods for viewing and analyzing anatomical structures | |
| US8107703B2 (en) | Quantitative real-time 4D stress test analysis | |
| CN104346821B (en) | Automatic planning for medical imaging | |
| US9741131B2 (en) | Anatomy aware articulated registration for image segmentation | |
| US7492968B2 (en) | System and method for segmenting a structure of interest using an interpolation of a separating surface in an area of attachment to a structure having similar properties | |
| JP5295562B2 (en) | Flexible 3D rotational angiography-computed tomography fusion method | |
| CN102231963A (en) | Reparametrized bull's eye plots | |
| CN104885126B (en) | The Computer assisted identification of tissue of interest | |
| AU2012350363A1 (en) | Method and apparatus for the assessment of medical images | |
| CN103093424A (en) | Method and apparatus for generating enhanced image from medical imaging data | |
| CN1920882A (en) | System and method for salient region feature based 3d multi modality registration of medical images | |
| CN107194909A (en) | Medical image-processing apparatus and medical imaging processing routine | |
| CN110807770A (en) | Medical image processing, recognizing and displaying method and storage medium | |
| EP4322836A1 (en) | Systems and methods for reconstruction of 3d images from ultrasound and camera images | |
| US20100303314A1 (en) | Systems and methods for detecting and visualizing correspondence corridors on two-dimensional and volumetric medical images | |
| Cucchiara et al. | An image analysis approach for automatically re-orienteering CT images for dental implants | |
| CN114943714A (en) | Medical image processing system, medical image processing apparatus, electronic device, and storage medium | |
| US9558568B2 (en) | Visualization method for a human skeleton from a medical scan | |
| US20030154201A1 (en) | Data storage format for topography data | |
| JP4807819B2 (en) | Image processing device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |