CN1739453A - System and method for 3D visualization of lung perfusion or density and statistical analysis thereof - Google Patents
System and method for 3D visualization of lung perfusion or density and statistical analysis thereof Download PDFInfo
- Publication number
- CN1739453A CN1739453A CNA2005100911547A CN200510091154A CN1739453A CN 1739453 A CN1739453 A CN 1739453A CN A2005100911547 A CNA2005100911547 A CN A2005100911547A CN 200510091154 A CN200510091154 A CN 200510091154A CN 1739453 A CN1739453 A CN 1739453A
- Authority
- CN
- China
- Prior art keywords
- lung
- perfusion
- parenchyma
- cut apart
- view data
- 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 and method for 3D visualization of a pair of lungs are provided. The method comprises: segmenting (220)image data of the pair of lungs and lung parenchyma; generating a 3D map(230) as a function of the segmented image data; and rendering the 3D map (240)as a color-coded semi-transparent 3D volume, wherein an opaque region highlights an area of interest.
Description
The application requires in the U.S. Provisional Application No.60/600 of submission on August 10th, 2004,257 rights and interests, and the copy of this application is introduced into for referencial use at this.
Technical field
The present invention relates to medical image analysis, relate in particular to system and method that the three-dimensional (3D) that is used for lung perfusion or density shows and the lung perfusion of analyzing the patient or Density Distribution with diagnostic method.
Background technology
When a blood clot from venous thrombosis (DVT) comes off and moves to the tremulous pulse in the lung and when stopping up tremulous pulse, infringement lung there and heart being produced pressure, pulmonary infarction takes place.The dangerous potentially life of this short-term complication and in patient with acute DVT phenomenon incidence rate be approximately 10.Because most of thromboembolisms do not have symptom when taking place, thus its may in addition than common recognize more common.
Most of people can recover fully from DVT and pulmonary infarction.Yet, big pulmonary infarction almost can choked flow to all blood flows of pulmonary and cause sudden death.In addition, pulmonary infarction can produce serious pressure to heart.After ischemic heart desease and apoplexy, pulmonary infarction is that heart disease causes the 3rd dead main cause.Yet it may be the most general preventible cause of death in hospital.
If the character of known pulmonary infarction, in time diagnosis is crucial.Yet the blood flow how the assessment embolus influences in the lung also is important.Recently, growing for the research interest of the automated process that detects the lung embolus according to high resolution computer tomography angiography (CTA).In addition, to being used to show and assessing that the degree that the perfusion that caused by pulmonary infarction lacks and the method for position have produced interest.This technology is utilized multi-lamina computer fault imaging (CT) machine, and this multi-lamina computer tomograph produces cut apart with identification or the inferior embolus of cutting apart of two dimension (2D) lamella more than 600 or 600 for each patient routinely.Yet this meeting is consuming time and itself can not cause the demonstration immediately of lung perfusion.
Yet research worker has proposed a kind of intravenous contrast agent experimental technique of the 2D demonstration of the interior lung perfusion of parenchyma afterwards of using recently.This method produces the dabbling 2D of the coloud coding parenchyma that covers on the original CT image and shows.Although this information has been passed on some Useful Informations, be difficult to follow the tracks of unusual dabbling zone.And if there is the lung embolus in contiguous such zone, following the tracks of it so is very difficult with respect to this regional position.
Therefore, perfusion that reduce or unusual zone or abnormal area can be discerned lung embolus or other unusual technology in the highlight parenchyma thereby a kind of 3D of observation of needs schemes, and need a kind of technology of such data to diagnose that be used to analyze.
Summary of the invention
The system and method that shows by the 3D that is provided for lung perfusion or density and by being provided for analyzing lung perfusion among the patient or Density Distribution with diagnostic method, aforementioned and other problem that the present invention has overcome in known teachings to be run into, wherein said 3D shows and allows the doctor to check whole lung volume in single view.
In one embodiment of the invention, being used for the method that the 3D of a pair of lung shows comprises: this view data to lung and lung parenchyma is cut apart; Produce 3D figure as the function of the view data after cutting apart; With reproduce described 3D figure as the translucent 3D volume of coloud coding, wherein the zone of opacity highlight area-of-interest.
Described 3D figure is in perfusion figure or the density map and can calculates by a quilt in adaptive smooth or the texture filtering.The zone that described area-of-interest indication has unusual perfusion or density.It is described that to have unusual dabbling zone be because embolus.
Described method further comprises: the rectangular histogram that produces the 3D figure that is reproduced; Determine that described rectangular histogram represents that the existence of embolus is the positive or feminine gender; With embolus is classified as a kind of in acute or chronic.
In another embodiment of the present invention, being used for the method that the 3D of lung perfusion shows comprises: the view data to lung parenchyma is cut apart; The perfusion figure of the view data after generation is cut apart; With reproduce described perfusion figure as the translucent 3D volume of coloud coding, wherein opaque perfusion zone is visible.
The step of the perfusion figure of the view data after generation is cut apart comprises: lung parenchyma is cut apart; Carry out local level and smooth; Local neighborhood average density with definite lung parenchyma.The step that lung parenchyma is cut apart comprises: cut apart this volume to lung from the thorax volume; Air flue and blood vessel in the lung volume that identification is cut apart; With cover figure (mask) by what remove from the lung volume of being cut apart that described air flue and blood vessel structure produce described lung parenchyma.
Carrying out local level and smooth step comprises: move described view data; Cover the view data that moves in the figure shelter with described parenchyma to obtain the parenchyma image; Described parenchyma is covered figure and image to be carried out Gauss and smoothly covers figure and image to obtain level and smooth parenchyma; Cover figure with described parenchyma and shelter described level and smooth parenchyma image; The balanced level and smooth parenchyma image of being sheltered; With move the isostatic image of institute to produce described perfusion figure.
A kind of in visible perfusion zone indication lung embolus or the dispersivity pneumonopathy.Described view data is by using one in CT, spiral CT or the MR imaging technique to obtain.This method further comprises regulates cromogram to observe interested high-density region.Described visible perfusion zone be low perfusion, HT or pour into unusually in a kind of.
In yet another embodiment of the present invention, being used for the system that the 3D of a pair of lung shows comprises: be used for stored program memory device; With the processor of described memory device communication, described processor can with described program operate with: this view data to lung and lung parenchyma is cut apart; Generation is as the 3D figure of the function of the view data of being cut apart; With reproduce described 3D figure as the translucent 3D volume of coloud coding, wherein the zone of opacity highlight area-of-interest.
Described 3D figure is in perfusion figure or the density map and can calculates by a quilt in adaptive smooth or the texture filtering.The zone that described area-of-interest indication has unusual perfusion or density.It is described that to have unusual dabbling zone be because embolus.
Described processor can be further operated with the rectangular histogram that produces the 3D figure that is reproduced with program code; Determine that described rectangular histogram represents that the existence of embolus is the positive or feminine gender; With embolus is classified as a kind of in acute or chronic.
In another embodiment of the present invention, being used for the system that the 3D of lung perfusion shows comprises: be used for stored program memory device; With the processor of described memory device communication, described processor can with described program move with: the view data to lung parenchyma is cut apart; Produce the perfusion figure of the view data of being cut apart; With reproduce described perfusion figure as the translucent 3D volume of coloud coding, wherein opaque perfusion zone is visible.
When producing the perfusion figure of the view data cut apart described processor can be further with described program code move with: lung parenchyma is cut apart; Carry out local level and smooth; Local neighborhood average density with definite lung parenchyma.When lung parenchyma is cut apart described processor can be further with described program code move with: cut apart this volume from the thorax volume to lung; Air flue and blood vessel in the lung volume that identification is cut apart; With by remove the figure that covers that described air flue and blood vessel structure produce described lung parenchyma from the lung volume of being cut apart.
When carry out local when level and smooth described processor can be further with described program code move with: move described view data; Cover the view data that moves in the figure shelter with described parenchyma to obtain the parenchyma image; Described parenchyma is covered figure and image to be carried out Gauss and smoothly covers figure and image to obtain level and smooth parenchyma; Cover figure with described parenchyma and shelter described level and smooth parenchyma image; The balanced level and smooth parenchyma image of being sheltered; With move the isostatic image of institute to produce described perfusion figure.
A kind of in described visible perfusion zone indication lung embolus or the dispersivity pneumonopathy.Described view data is by using one in CT, spiral CT or the MR imaging device to obtain.Thereby described processor can further move to regulate cromogram with described program code and observe interested high-density region.Described visible perfusion zone be low perfusion, HT or pour into unusually in a kind of.
Aforesaid feature to belong to representational embodiment and to propose in order helping to understand the present invention.Should be appreciated that them and do not mean that to be counted as, perhaps to the restriction of the equivalent of claim to restriction of the present invention defined by the claims.Therefore, should not to be looked at as when determining equivalent be conclusive in the summary of feature.Other feature of the present invention will become apparent in the following description and from accompanying drawing and claim.
Description of drawings
Fig. 1 according to an exemplary embodiments of the present invention, be used for the block chart of the system that the 3D of lung perfusion or density shows;
Fig. 2 is the flow chart of explanation according to the method for a 3D demonstration exemplary embodiments of the present invention, that be used for lung perfusion or density;
Fig. 3 be explanation according to an exemplary embodiments of the present invention, perfusion is reproduced and a pair of image of CT data lamella and this CT data lamella of not reproducing;
Fig. 4 is explanation according to the perfusion figure of the lamella of a Fig. 3 exemplary embodiments of the present invention, that reproduce with 3D and not have a pair of image of perfusion figure of lamella of Fig. 3 of reproduction;
Fig. 5 a pair of image that to be explanation show according to the 3D of the positive patient's of a patient exemplary embodiments of the present invention, that pulmonary infarction is negative and pulmonary infarction perfusion figure;
Fig. 6 be explanation according to an exemplary embodiments of the present invention, be used for the flow chart of level and smooth lung parenchyma with the method that produces perfusion figure;
Fig. 7 is the flow chart of explanation according to a method exemplary embodiments of the present invention, that be used to analyze lung perfusion or density;
Fig. 8 is that explanation is according to a chart exemplary embodiments of the present invention, that be used for the cromogram of lung perfusion demonstration;
Fig. 9 be explanation according to an exemplary embodiments of the present invention, be used for lung perfusion and show and be used for that lung perfusion shows and a pair of chart of the opacity figure of edema;
Figure 10 is the histogrammic a pair of chart of explanation according to the positive patient's of a patient exemplary embodiments of the present invention, that pulmonary infarction is negative and pulmonary infarction perfusion figure;
Figure 11 be explanation according to an exemplary embodiments of the present invention, have the patient of chronic pulmonary thromboembolism and have the patient's of hydrops lung the histogrammic a pair of chart of perfusion figure;
Figure 12 is explanation according to a scatter diagram exemplary embodiments of the present invention, relevant with many patients' the variance of perfusion figure that is used for the patient is divided into some disease type, the peaked Hounsfeld of rectangular histogram (HU) value.
The specific embodiment
Fig. 1 according to an exemplary embodiments of the present invention, be used for the block chart of the system 100 that the 3D of lung perfusion or density shows.As shown in Figure 1, wherein system 100 comprises scanning device 105, personal computer (PC) 110 and for example is connected operator's control station 115 on the Ethernet 120.Scanning device 105 can be magnetic resonance (MR) imaging device, CT imaging device, spiral CT equipment, positron emission computerized tomography (PET) equipment, 2D or 3D fluorescence imaging equipment, 2D, 3D or the four-dimension (4D) supersonic imaging apparatus, or X-ray equipment.Scanning device 105 also can be the hybrid image forming apparatus that can realize CT, MR, PET or other imaging technique.
PC 110 can be work station, portable or laptop computer, PDA(Personal Digital Assistant) etc., and it comprises CPU (CPU) 125 and the memorizer 130 that is connected on input equipment 150 and the outut device 155.CPU 125 comprises display module 145, and this display module comprises the method that one or more 3D that are used for lung perfusion or density show.
Memorizer 130 comprises random-access memory (ram) 135 and read only memory (ROM) 140.Memorizer 130 also can comprise data base, disc driver, tape drive etc. or its combination.RAM 135 play the employed data of program run duration that are stored among the CPU 125 data storage effect and be used as the working area.ROM 140 plays the effect that is stored in the program storage of the program of operation among the CPU 125.Input equipment 150 is by keyboard, compositions such as mouse, and outut device 155 is made up of liquid crystal display (LCD), cathode ray tube (CRT) display or printer.
From the operation of operator's control station 115 control system 100, operator's control station comprises for example controller 165 of keyboard and the display 160 of for example CRT monitor.Operator's control station 115 is communicated by letter with scanning device 105 with PC 110, so that the 2D view data of being collected by scanning device 105 can be translated into the 3D data and be checked on display 160 by PC 110.Be to be understood that, PC 110 can be provided under the situation that does not have operator's control station 115 and for example to use input 150 and output 155 equipment to operate and the information that is provided by scanning device 105 is provided, with some task of implementing to be carried out by controller 165 and display 160.
Operator's control station 115 comprises that further the Digital Image Data that can handle the image data set (or its part) that is obtained is to produce and to show any suitable image reproducing system/tool/application of 2D and/or 3D rendering on display 160.More specifically, image reproducing system can provide the 2D/3D reproduction of medical image and the application program that shows, this application program moves on universal or special computer workstation.And image reproducing system can make the user navigate on 3D rendering or a plurality of 2D image slice.PC 110 can comprise that also the Digital Image Data that is used to handle the image data set that is obtained is to produce and to show the image reproducing system/tool/application of 2D and/or 3D rendering.
As shown in Figure 1, display module 145 also can be used to receive and to handle digital medical image data by PC 110, as mentioned above, the form of this digital medical image data can be a raw image data, 2D data reconstruction (for example, axial slice), or 3D data reconstruction, for example volumetric image data or many planar reformat, or any combination of such form.Data processed result can output to image reproducing system operator's control station 115 through network 120 from PC 110, with according to the cutting apart of described data processed result, for example organ or anatomical structure, 2D and/or 3D that colour or variable density or the like produce view data reproduce.
Fig. 2 is the flow chart that shows the operation of the method that shows according to 3 a D exemplary embodiments of the present invention, that be used for lung perfusion or density.As shown in Figure 2, the CT view data obtains (210) from a pair of lung of patient inside.This realizes that by using scanning device 105 scanning device is a CT scanner in this example, and at operator's control station 115 these scanning devices 105 of operation, with the scan patients lung, produces a series of 2D image slice relevant with lung thus.The 2D image slice of lung is combined then to form 3D rendering.Except lung, should be understood that the CT view data can or comprise any one acquisition other body part of blood vessel from lower limb, arm, brain.In addition, can use the data of other type according to the present invention, for example the MR view data.
After the CT view data was acquired, this was cut apart (220) to the parenchymatous view data in the lung.The fragmentation technique of the image of can use traditional cutting techniques, the zones of different of image being cut apart such as the cutting techniques of edge that detects the object that will cut apart or profile or around distinguishing with identification is cut apart parenchymatous view data.In case view data is cut apart, produce the perfusion figure (230) of the view data of being cut apart.
Yet should be understood that except produce perfusion figure in this step, the function of the view data in the view data of being cut apart can be obtained.This function can be any function that is used for exporting based on the local message in the zone of being cut apart digital value.What for example, rich such as adding (Gabor) filter was such can be used to produce the figure of highlight special texture for special texture provides the filter of high response value.In addition, the output that this figure can be by for example making up perfusion figure and texture filter produces with the combination of this function that produces new drawing.
Perfusion figure produces by the adaptive smooth that uses averaging operator to carry out the image of being cut apart.Pour into figure then as the translucent 3D volume of coloud coding reproduced (240).An example like this is shown in the image (a) of Fig. 3.As shown in Figure 3, image (a) is the original lamella with CT data that the perfusion of covering shows.Shown in image (a), spherical 310a represents the position of lung embolus, the low perfusion zone in the blue opaque speckle 320 expression parenchymas, the average perfusion of green translucent speckle 330 expressions zone, and red transparent speckle 340 expression HT zones.In other words, speckle 320 expressions lack the zone of blood flow, and speckle 330 expressions have health or normal dabbling zone, and speckle 340 expressions have the density of increase or the zone of unusual HT.
Shown in Fig. 3 was further, image (b) was the original lamella with CT data of perfusion reproduction.In other words, image (b) is the original lamella that has CT data in the volume range of lung embolus position 310b of doctor's labelling.The lung embolus position 310b of institute's labelling also can obtain by lung embolus detection algorithm.Fig. 3 that is discussed below and the view of Fig. 4 can be by the doctor by simple switching the between initial data and coloud coding data obtains from work station or operator's control station 115.
The image of Fig. 4 (a) has shown the 2D view of the perfusion figure of the 3D perfusion figure that has covering.In addition, the image of Fig. 4 (b) is the view that 3D reproduces pent image (a), thereby realizes the 2D side view of perfusion figure.Can find that from the review of the image of Fig. 3 and 4 no matter from which visual angle, the zone with minimum dabbling lung is least transparent and the most visible.
Therefore, the doctor can identify apace and for example have low dabbling suspicious region (for example, the zone that influenced by the lung embolus) as the image (a) of Fig. 5 with in the lung (b).As shown in Figure 5, image (a) is that the 3D of perfusion figure of patient's CTA shows, as indicated by the even dabbling speckle of expression 510, so the result of this patient's pulmonary infarction is negative.Image (b) is that the 3D of perfusion figure of patient's CTA shows, as by the transparent speckle 530 of the low dabbling opaque speckle 520 of expression, expression HT and indicated by the lung embolus 540 of doctor institute labelling, the result of this patient's pulmonary infarction is positive.
To discuss perfusion figure now in more detail calculates.More specifically, calculate perfusion figure by the local neighborhood average density of determining lung parenchyma.This is by at first cutting apart lung parenchyma, carrying out the part then and smoothly realize.
Cutting apart by three steps of lung parenchyma formed: cut apart lung volume from the thorax volume; Identification air flue and blood vessel; And by removing described air flue and blood vessel structure produces the described parenchymatous figure of covering from described cutting apart.Lung volume is cut apart and can be realized that by using high threshold values region growing this high threshold values region growing starts from a seed points in the trachea, carries out form closure (morphological closing) afterwards.Air and blood volume elements in the lung volume that can use basic threshold to discern then to be cut apart.For example, can use be used for air-the last threshold values of 990HU and be used for blood-the following threshold values of 330HU.Should be understood that, can utilize other being used to discern the dividing method of air flue and vascular tree since after will typically follow level and smooth step, therefore accurately cut apart fully and there is no need.
In case it is available that parenchyma is covered figure, just partly carry out adaptive smooth to dropping on this original CT volume of covering among the figure.Should be understood that, can use any level and smooth/averaging operator, as long as the air that deducts/blood volume elements does not have on average to arrive in the level and smooth parenchyma.The 3D gaussian kernel is used as smoothing operator in the example that is described below.
It is shown in Figure 6 and be described below with the exemplary steps that produces perfusion figure to be used for level and smooth lung parenchyma.At first, with the original CT volume 1024HU that moves up, thereby-it is zero (605) that 1024HU becomes.Then, cover figure (610) with parenchyma and shelter CT volume (615) after moving with the parenchyma image after obtaining to move, wherein drop on All Ranges outside the described figure of covering in this image and be set as zero, all other zones keep the value (620) after the moving of they simultaneously.Then the parenchyma image after moving and parenchyma being covered figure carries out Gauss and smoothly covers figure (630) to produce level and smooth parenchyma image (625) and level and smooth parenchyma.
Level and smooth parenchyma is covered figure and is comprised the relative fractional value of binary value of covering figure with initial parenchyma.Described fractional value defines the volume elements outer with covering figure and compares from covering the percentage ratio of the volume elements acquisition particular voxel smooth value in the figure.The respective element value of covering among the figure divided by level and smooth parenchyma by each nonzero element comes balanced (635) level and smooth parenchyma to cover figure then, only comprises from the image of the voxel values that has those acquisitions of covering figure with generation.The image that produces is covered figure with initial parenchyma then and is sheltered, and moves down 1024HU, so that element value is once more corresponding to HU scale (640).
Except being discerned, the doctor has the low dabbling suspicious region, the existence that the rectangular histogram of described perfusion figure or a plurality of perfusion figure can be used to definite patient's lung embolus is the positive or feminine gender, and helps classification to exist the lung embolus and the patient of which kind of type whether to have other situation, for example edema or pneumonia.Fig. 7 is that explanation is according to the flow chart that is used to analyze the method for a lung perfusion exemplary embodiments of the present invention, a kind of like this.
As shown in Figure 7, after producing perfusion figure (710), produce one or more patients' multiple rectangular histogram (720) according to above-mentioned technology.Before producing and analyzing described rectangular histogram, cromogram is used to the value of the perfusion figure that is reproduced painted.Typical cromogram is shown in Figure 8.According to colored (cold-to-hot color) scheme of heat that is as cold as shown in Figure 8, normal perfusion represents that with green the perfusion that reduces represents that with blueness HT is represented with redness.The central value of-800HU is shown and is selected in this example makes compromise and utilizes cromogram in the available patient between the comparability between the patient.
Should be understood that, cromogram can be by the interactively translation and be selected to for example have the 100HU width or such as the center for the window setting of the variation of-900HU between low perfusion zone, to distinguish, perhaps the center for the window setting of the variation of-600HU between the HT zone, to distinguish.In addition, center and window value can be carried out the adaptability setting based on the value that derives from histogram analysis.
The setting that use is limited by cromogram can produce the piecewise linearity opacity figure of perfusion figure.Typical piecewise linearity opacity figure is shown in the chart (a) of Fig. 9.Shown in chart (a), the shown dabbling zone that reduces that has has higher opacity than parenchymatous other parts basically.In other words, from approximately-1000HU is very opaque 3D perfusion reproduced image to the zone of a little higher than-900HU.Should be understood that this drawing can be adjusted by doctor's interactively.
Scheme as an alternative can produce highlight and have the opacity figure that the dabbling zone that reduces and highlight have unusual high density or dabbling parenchyma zone.The example of such opacity figure is shown in the chart (b) of Fig. 9.Because it is opaque having high dabbling zone and they are usually represented such as the problem area that has pressure in the lung, therefore can produce this figure.Should be understood that also because volume reproduces and to be used to reproduce perfusion figure, so opacity and value of color can be adjusted immediately.In addition, any manually or automatic technique can be used to determine these values.
After producing opacity figure, can produce and analyze the rectangular histogram of perfusion figure then.When analyzing rectangular histogram, the negative patient of lung embolus trends towards having such parenchyma Luminance Distribution, this parenchyma Luminance Distribution major part is symmetric and has low broadening (spread), and therefore can be characterized by their two and third moment (moment).In addition, the negative histogrammic maximum of lung embolus typically is positioned at limited HU scope.The patient's that the lung embolus is negative example histogram is shown in the chart (a) of Figure 10.
Further analyzing rectangular histogram and/or when wherein extracting statistics, can distinguish acute and the chronic pulmonary embolus.Because the latter's perfusion distributed pole is asymmetric and owing to the characteristic pattern of mosaic sign (mosaic attenuation) has big broadening, so this is possible.For this point is described, have acute pulmonary embolism, chronic pulmonary thromboembolism and have edema or the patient's of the pulmonary infarction of pneumonia example histogram respectively at the chart (a) of the chart (b) of Figure 10 and Figure 11 with (b).
Use above-mentioned data, the positive and negative pulmonary infarction case (730) of can being distinguished from each other and can discerning such as edema, acute and situation (740) that the chronic pulmonary thromboembolism is such.Can how to discern such as the example of edema, the acute and situation that the chronic pulmonary thromboembolism is such scatter diagram and illustrate by Figure 12.This scatter diagram is by making relevant generation of variance of the peaked HU value of rectangular histogram and each patient's perfusion figure.Especially, Figure 12 has shown 19 patients' that are considered to have pulmonary infarction scatter diagram, and these patients estimate by CT after having used the IV contrast agent.The CT data set is evaluated by experienced breast radiologist.
As shown in figure 12,17 patients are diagnosed as has a plurality of lung emboluses, and two patients are lung embolus feminine gender.Patient with a plurality of lung emboluses has 2-28 embolus and each patient on average has 13 emboluses.As shown in the figure, for example by the caption among Figure 12, three patients are diagnosed as has the chronic pulmonary thromboembolism; Yet, only have two to have the mosaic sign that causes because of the long-term perfusion shortage relevant among these patients with the chronic pulmonary thromboembolism.As the caption among Figure 12 further indication, the patient can be divided into four groups: pulmonary infarction (PE) feminine gender; The PE positive is acute; The PE positive chronic; Positive and have other pneumonopathy with PE.Figure 10 and 11 rectangular histogram be corresponding to these four groups, and to its analysis therefore the patient can be classified.
According to an exemplary embodiments of the present invention, parenchymatous perfusion figure is reproduced as the translucent 3D volume of coloud coding, wherein has the dabbling zone that reduces unusually and is reproduced relatively opaquely.This is convenient to not rely on observation place and direction and quick identification goes out to have the suspicious region of unusual lung perfusion.In addition, by extracting statistical nature from perfusion figure, these features can be used to distinguish patient with pulmonary infarction and the patient who does not have pulmonary infarction.For example, by discerning the mosaic sign that is caused by the chronic pulmonary thromboembolism usually, the chronic pulmonary thromboembolism can be distinguished mutually with acute pulmonary embolism.And because these features can be used to distinguish the unusual of lung perfusion or density, so they can be used to discern the patient and whether have other situation, for example the dispersivity pneumonopathy of pneumonia or other type.
Should be understood that further that the present invention can make up with various forms of hardware, software, firmware, application specific processor or its and realize.In one embodiment, the present invention can be used as the software that is comprised in the application program in the program storage device (for example, magnetic floppy disk, RAM, CD ROM, DVD, ROM, and flash memory) really and realizes.This application program can upload to the machine that comprises any suitable construction and by its execution.
Should be understood that further that because some the system's building blocks and the method step that are described in the drawings can realize with software, so the actual connection between the system unit (or method step) can be depended on mode that the present invention is programmed and different.If the instruction of the present invention that provides here is provided, those of ordinary skill in the art can imagine these and similarly embodiment or structure of the present invention.
Should be appreciated that also top description only represents schematic embodiment.For convenience of the reader, top description concentrates on the representative example of possible embodiments, and this example has illustrated principle of the present invention.Described description is not to attempt not have to enumerate all possible flexible program with omitting.Alternate embodiment is not to propose at special door section of the present invention, and perhaps other replacement scheme of also not describing can be used for a part, and what be not considered to those alternate embodiments negates.Under the situation of the spirit and scope of the present invention that do not come off, can realize other application and embodiment.
Owing to can produce above-mentioned multiple permutation and combination and the above-mentioned enforcement that comprises the non-creativeness replacement scheme, therefore mean that the present invention is not limited to specifically described embodiment, but the present invention will be defined according to claim subsequently.Be understandable that in many all below the literal scope of claim among those embodiment that do not describe, and other is of equal value.
Claims (29)
1. be used for the method for the 3D demonstration of a pair of lung, comprise:
This view data to lung and lung parenchyma is cut apart;
Generation is as the 3D figure of the function of the view data of being cut apart; With
Reproduce described 3D figure, wherein zone of opacity highlight area-of-interest as the translucent 3D volume of coloud coding.
2. method according to claim 1, wherein said 3D figure is a kind of in perfusion figure or the density map.
3. method according to claim 1, wherein said 3D figure is by a generation in adaptive smooth or the texture filtering.
4. method according to claim 1, the zone that wherein said area-of-interest indication has unusual perfusion or density.
5. method according to claim 4, wherein said to have unusual dabbling zone be because embolus.
6. method according to claim 1 further comprises:
Produce the rectangular histogram of the 3D figure that is reproduced; With
Determine that described rectangular histogram represents that the existence of embolus is the positive or feminine gender.
7. method according to claim 6 further comprises:
Embolus is classified as a kind of in acute or chronic.
8. be used for the method for the 3D demonstration of lung perfusion, comprise:
View data to lung parenchyma is cut apart;
Produce the perfusion figure of the view data of being cut apart; With
Reproduce described perfusion figure as the translucent 3D volume of coloud coding, wherein opaque perfusion zone is visible.
9. method according to claim 8, the step that wherein produces the perfusion figure of the view data of being cut apart comprises:
Lung parenchyma is cut apart;
Carry out local level and smooth; With
Determine the local neighborhood average density of lung parenchyma.
10. method according to claim 9 wherein comprises the step that lung parenchyma is cut apart:
Cut apart this volume from the thorax volume to lung;
Air flue and blood vessel in the lung volume that identification is cut apart; With
By from the lung volume of being cut apart, removing the figure that covers that described air flue and blood vessel structure produce described lung parenchyma.
11. method according to claim 10 is wherein carried out local level and smooth step and is comprised:
Move described view data;
Cover figure with described parenchyma and shelter view data after moving to obtain the parenchyma image;
Described parenchyma is covered figure and image to be carried out Gauss and smoothly covers figure and image to obtain level and smooth parenchyma;
Cover figure with described parenchyma and shelter described level and smooth parenchyma image;
The balanced level and smooth parenchyma image of being sheltered; With
Move the isostatic image of institute to produce described perfusion figure.
12. method according to claim 8, a kind of in wherein visible perfusion zone indication lung embolus or the dispersivity pneumonopathy.
13. method according to claim 8, wherein said view data is by using one in CT, spiral CT or the MR imaging technique to obtain.
14. method according to claim 8 further comprises:
Regulate cromogram to observe interested high-density region.
15. method according to claim 8, wherein visible perfusion zone be low perfusion, HT or pour into unusually in a kind of.
16. be used for the system that the 3D of a pair of lung shows, comprise:
Be used for stored program memory device;
With the processor of described memory device communication, described processor can with described program operate with:
This view data to lung and lung parenchyma is cut apart;
Generation is as the 3D figure of the function of the view data of being cut apart; With
Reproduce described 3D figure, wherein zone of opacity highlight area-of-interest as the translucent 3D volume of coloud coding.
17. system according to claim 16, wherein said 3D figure are in perfusion figure or the density map.
18. system according to claim 16, the zone that wherein said area-of-interest indication has unusual perfusion or density.
19. system according to claim 18, wherein said to have unusual dabbling zone be because embolus.
20. system according to claim 16, wherein said processor can be further with described program code operate with:
Produce the rectangular histogram of the 3D figure that is reproduced; With
Determine that described rectangular histogram represents that the existence of embolus is the positive or feminine gender.
21. system according to claim 20, wherein said processor can be further with described program code operate with:
Embolus is classified as a kind of in acute or chronic.
22. be used for the system that the 3D of lung perfusion shows, comprise:
Be used for stored program memory device;
With the processor of described memory device communication, described processor can with described program operate with:
View data to lung parenchyma is cut apart;
Produce the perfusion figure of the view data of being cut apart; With
Reproduce described perfusion figure as the translucent 3D volume of coloud coding, wherein opaque perfusion zone is visible.
23. system according to claim 22, wherein when producing the perfusion figure of the view data cut apart described processor can be further with described program code operate with:
Lung parenchyma is cut apart;
Carry out local level and smooth; With
Determine the local neighborhood average density of lung parenchyma.
24. system according to claim 23, wherein when lung parenchyma is cut apart described processor can be further with described program code operate with:
Cut apart this volume from the thorax volume to lung;
Air flue and blood vessel in the lung volume that identification is cut apart; With
By from the lung volume of being cut apart, removing the figure that covers that described air flue and blood vessel structure produce described lung parenchyma.
25. system according to claim 24, wherein when carry out local when level and smooth described processor can be further with described program code operate with:
Move described view data;
Cover figure with described parenchyma and shelter view data after moving to obtain the parenchyma image;
Described parenchyma is covered figure and image to be carried out Gauss and smoothly covers figure and image to obtain level and smooth parenchyma;
Cover figure with described parenchyma and shelter described level and smooth parenchyma image;
The balanced level and smooth parenchyma image of being sheltered; With
Move the isostatic image of institute to produce described perfusion figure.
26. system according to claim 22, a kind of in wherein visible perfusion zone indication lung embolus or the dispersivity pneumonopathy.
27. system according to claim 22, wherein said view data is by using one in CT, spiral CT or the MR imaging device to obtain.
28. system according to claim 22, wherein said processor can be further with described program code operation with:
Thereby regulate cromogram and observe interested high-density region.
29. system according to claim 22, wherein said visible perfusion zone be low perfusion, HT or pour into unusually in a kind of.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US60025704P | 2004-08-10 | 2004-08-10 | |
| US60/600257 | 2004-08-10 | ||
| US11/182,913 | 2005-07-15 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN1739453A true CN1739453A (en) | 2006-03-01 |
Family
ID=36092085
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNA2005100911547A Pending CN1739453A (en) | 2004-08-10 | 2005-08-10 | System and method for 3D visualization of lung perfusion or density and statistical analysis thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1739453A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101541242B (en) * | 2006-11-30 | 2012-06-13 | 皇家飞利浦电子股份有限公司 | Visualizing a vascular structure |
| CN111358484A (en) * | 2020-03-23 | 2020-07-03 | 广州医科大学附属第一医院（广州呼吸中心） | Nuclear medicine lung perfusion imaging quantitative analysis method, analysis equipment and storage medium |
| CN113469934A (en) * | 2020-03-31 | 2021-10-01 | 西门子医疗有限公司 | Assessment of abnormal regions associated with disease from chest CT images |
-
2005
- 2005-08-10 CN CNA2005100911547A patent/CN1739453A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101541242B (en) * | 2006-11-30 | 2012-06-13 | 皇家飞利浦电子股份有限公司 | Visualizing a vascular structure |
| CN111358484A (en) * | 2020-03-23 | 2020-07-03 | 广州医科大学附属第一医院（广州呼吸中心） | Nuclear medicine lung perfusion imaging quantitative analysis method, analysis equipment and storage medium |
| CN111358484B (en) * | 2020-03-23 | 2021-12-24 | 广州医科大学附属第一医院（广州呼吸中心） | Nuclear medicine lung perfusion imaging quantitative analysis method, analysis equipment and storage medium |
| CN113469934A (en) * | 2020-03-31 | 2021-10-01 | 西门子医疗有限公司 | Assessment of abnormal regions associated with disease from chest CT images |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7715608B2 (en) | System and method for 3D visualization of lung perfusion or density and statistical analysis thereof | |
| CN1127700C (en) | Data visualization enhancement through removal of dominating structures | |
| Eltonsy et al. | A concentric morphology model for the detection of masses in mammography | |
| Betal et al. | Segmentation and numerical analysis of microcalcifications on mammograms using mathematical morphology. | |
| JP5814504B2 (en) | Medical image automatic segmentation system, apparatus and processor using statistical model | |
| CN101036165A (en) | System and method for tree-model visualization for pulmonary embolism detection | |
| CN105559813B (en) | Medical diagnostic imaging apparatus and medical image-processing apparatus | |
| DE102005036412B4 (en) | Improved GGN segmentation in lung recordings for accuracy and consistency | |
| US20150356730A1 (en) | Quantitative predictors of tumor severity | |
| JP2004105731A (en) | Processing of computer aided medical image | |
| CN105719324A (en) | Image processing apparatus, and image processing method | |
| CN1818972A (en) | System and method for splicing medical image datasets | |
| DE102004055767A1 (en) | Method and apparatus for processing image data to assist in the detection of a disease | |
| JP2008515466A (en) | Method and system for identifying an image representation of a class of objects | |
| JP2011514822A (en) | Method and system for segmenting CT scan data | |
| CN1656499A (en) | Pulmonary nodule detection using cartwheel projection analysis | |
| JP2021073550A (en) | Method of assisting viewing of images and device using same | |
| Mabrouk et al. | Computer aided detection of large lung nodules using chest computer tomography images | |
| Heuberger et al. | Lung CT segmentation for image retrieval using the Insight Toolkit (ITK) | |
| CN1739453A (en) | System and method for 3D visualization of lung perfusion or density and statistical analysis thereof | |
| Yeshua et al. | Deep learning for detection and 3D segmentation of maxillofacial bone lesions in cone beam CT | |
| WO2005029410A1 (en) | Method and system for ground glass nodule (ggn) segmentation with shape analysis | |
| CA2531871C (en) | System and method for detecting a protrusion in a medical image | |
| Xu et al. | Difficulty-aware bi-network with spatial attention constrained graph for axillary lymph node segmentation | |
| Bellotti et al. | The MAGIC-5 project: Medical applications on a grid infrastructure connection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20061201 Address after: American Pennsylvania Applicant after: American Siemens Medical Solutions Inc. Address before: new jersey Applicant before: Siemens Corporate Research, Inc. |
|
| 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 |