CN103109308B - The quantization of the tube chamber characteristic of tubular structure - Google Patents
The quantization of the tube chamber characteristic of tubular structure Download PDFInfo
- Publication number
- CN103109308B CN103109308B CN201180045151.3A CN201180045151A CN103109308B CN 103109308 B CN103109308 B CN 103109308B CN 201180045151 A CN201180045151 A CN 201180045151A CN 103109308 B CN103109308 B CN 103109308B
- Authority
- CN
- China
- Prior art keywords
- tubular structure
- voxel
- interested
- distance map
- diameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 230000000875 corresponding Effects 0.000 claims abstract description 8
- 210000004204 Blood Vessels Anatomy 0.000 claims description 26
- 206010002329 Aneurysm Diseases 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 abstract description 15
- 238000011002 quantification Methods 0.000 abstract description 6
- 238000000034 method Methods 0.000 description 11
- 239000002872 contrast media Substances 0.000 description 10
- 238000007689 inspection Methods 0.000 description 9
- 200000000009 stenosis Diseases 0.000 description 6
- 230000036262 stenosis Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000903 blocking Effects 0.000 description 3
- 238000002595 magnetic resonance imaging Methods 0.000 description 3
- 210000004369 Blood Anatomy 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006011 modification reaction Methods 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 206010003210 Arteriosclerosis Diseases 0.000 description 1
- 210000001072 Colon Anatomy 0.000 description 1
- 241000565118 Cordylophora caspia Species 0.000 description 1
- 206010061529 Polyp Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 201000001320 atherosclerosis Diseases 0.000 description 1
- 230000003143 atherosclerotic Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- OYPRJOBELJOOCE-UHFFFAOYSA-N calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002452 interceptive Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000000414 obstructive Effects 0.000 description 1
- 210000000056 organs Anatomy 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000001932 seasonal Effects 0.000 description 1
- 210000001519 tissues Anatomy 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Abstract
A kind of method includes producing signal based on distance map, the property quantification of tubular structure tube chamber interested on the preset range of the tubular structure interested represented in described signal designation imaging in vivo data.A kind of system includes quantifying parts (216), and it produces signal based on distance map, the property quantification of tubular structure tube chamber interested on the preset range of the tubular structure interested represented in described signal designation imaging in vivo data.A kind of method includes the voxel identifying the 3D segmentation imaging data representing tubular structure interested；Determine the distance between identified voxel and the nearest voxel corresponding to other structures in multiple 2D section of described segmentation imaging data；Utilize and represent described voxel to described apart from proportional intensity level；And the signal of instruction 3D distance map is produced based on intensity level.
Description
Technical field
Hereafter relate in general to quantify the characteristic of the tube chamber of tubular structure, here in connection with being applied particularly to meter
Calculation tomographic (CT) and quantifying and narrow being been described by of human or animal's person under inspection's blood vessel.But,
Herein below can also be revised to adapt to other application and/or other image modes, such as nuclear magnetic resonance
(MRI), ultrasonic (US) imaging and/or other image modes.
Background technology
Accurately measuring for determining fitting of the patient suffering from atherosclerotic stenosis disease of angiostenosis degree
It is critically important when clinic gets involved.Use medical imaging acquisition data, can obtain narrow from data
The tolerance of narrow degree.The suitably example of medical imaging flow process is three-dimensional (3D) angiogram, its
There is in blood vessel and surrounding tissue different image pixel intensities, it is possible to use blood pool contrast agent or blood flow effect
Should realize.Adopted by computerized axial tomography (CTA), nuclear magnetic resonance (MRA) and doppler ultrasound
Collection high resolution 3 d angiogram.
Regrettably, in order to determine such tolerance, clinicist must manually pass through 3D volume
Roll a large amount of (such as more than 100) two dimension (2D) image or section (the most axial, crown,
Sagittal, inclination etc.) narrow to position.Once finding, it is narrow that clinicist just rolls through covering
Cut into slices with located slices, block for quantification of stenosis, the most visually seem that there is maximum obstruction
Section.But, for clinicist, this time-intensive of mass data collection is manually commented
Estimate and may be difficult to bear.Depend on observer additionally, measure and be vulnerable to personal error impact.Such as,
The section and the cutting angle that select can strengthen or weaken narrow visual appearance.
Once positioning and identify such image, it is narrow to find that clinicist just measures blood vessel diameter
The minimum diameter at place and blood vessel diameter.Based on software or other survey tool can be passed through, utilize multiple
Miscellaneous model do so, to calculate blood vessel diameter in orthogonal cross-sections.From diameter measurement, clinical doctor
Teacher can calculate various tolerance, and such as NASCET(North America symptom carotid endarterectomy is tested)
Standard, it is according to following percentage ratio quantification of stenosis: (1-(minimum blood vessel diameter/blood vessel diameter)) ×
100%.This provides the intuition numeral between 0% and 100%, and wherein 0% represents do not have stenosis,
100% represents total blockage.
Summary of the invention
The each side of the application solves the problems referred to above and other problems.
According to an aspect, a kind of method includes producing signal based on distance map, and described signal designation exists
The tube chamber of tubular structure interested on the predetermined of the tubular structure interested represented in body imaging data
Property quantification.
According on the other hand, a kind of system includes quantifying parts, and it produces signal based on voxel distance map,
In described signal designation imaging in vivo data represent tubular structure interested predetermined on interested
The quantization of the tube chamber characteristic of tubular structure.
According on the other hand, a kind of method includes: identify the three-dimensional segmentation representing tubular structure interested
The voxel of imaging data；Identified voxel is determined in multiple two dimension slicings of described segmentation imaging data
And corresponding to the distance between the nearest voxel of other structures；Utilize to described apart from proportional intensity
Value represents described voxel；And the signal of three-dimensional distance figure is shown based on intensity level generation compass.
Accompanying drawing explanation
The present invention can be arranged by various parts and parts, and various step and procedure and have
Body.Accompanying drawing, just to illustrating preferred embodiments, is not necessarily to be construed as limiting the present invention.
Fig. 1 shows the imaging system combining assessment processor.
Fig. 2 shows that example assesses processor.
Fig. 3 shows that the example of the tubular structure diameter of the two dimension slicing function as 3-D data set is bent
Line chart.
Fig. 4 shows the method for the tube chamber characteristic for quantifying tubular structure interested.
Fig. 5 shows the non-limiting method producing distance map.
Fig. 6 shows for quantifying the non-limiting method of characteristic on identified volume based on distance map.
Detailed description of the invention
Fig. 1 shows imaging system 100, such as computerized axial tomography (CT) scanner unit.At other
In embodiment, imaging system can additionally or alternatively include different image modes, such as magnetic resonance
Imaging (MRI), ultrasonic (US) imaging, combination type modality scans machine and/or other scanner units.Join
Put the imaging system 100 of diagram to perform contrast-agent-free scanning and scanning based on contrast agent, such as blood
The scanning of pipe visualization and other scannings based on contrast agent.
Imaging system 100 includes stationary gantry 102 and rotatable gantry 104, by stationary gantry
102 are pivotably supported rotatable gantry 104.Rotatable gantry 104 around inspection area 106 about
The longitudinal axis or z-axis rotate.Bearing 108, such as bed, support person under inspection in inspection area 106.?
Seat 108 can be used for before the scan, period and/or afterwards relative to x, y and/or z-axis with various sides
Formula location person under inspection.
Radiation source 110, such as x-ray tube are supported and along with rotatable gantry by rotatable gantry 104
104 rotate around inspection area 106 together, and launch radiation, radiate traverses 106.One-dimensional
Or two-dimensional radiation sensitive detector array 112 crosses over inspection area 106 and is positioned at the relative of radiation source 110
Side also detects the radiation of traverses 106.Detector array 112 produces instruction detected radiation
Signal or data for projection.
Reconstructor 114 reconstructed projection data also produces three-dimensional (3D) the body figure indicating inspection area 106
As data (voxel).General-purpose computing system serves as operator's console 116, including such as display
Output device and the such as input equipment such as keyboard, mouse.Software resident on control station 116 allows behaviour
The operation of work person's control system 100, such as, it is allowed to operator's selection technique (such as contrast agent enhanced),
Driver sweep etc..
Syringe 118 is configured to injection or administration of contrast agents material, such as, for contrast agent enhanced
Irnaging procedures.The syringe 118 of diagram is controlled the control of platform 116.In another case,
Syringe 118 is controlled by user's Artificial Control or by another device.In still another case, injection is omitted
Device 118, artificial administration of contrast agents.
Assessment processor 122 is configured to assess view data.As described in more detail below, in one
In the case of, assessment includes produce the tubular structure interested represented in instruction view data one or many
The signal of the quantization of individual characteristic.As non-limiting example, in one case, tubulose interested
Structure includes blood vessel, this diameter that can include quantifying the tube chamber of blood vessel, and this can be so that such as passing through
Accurately and/or blocking or blocking or the tube chamber of broadening blood vessel of blood vessel can be quantified by playback system.Again, this
Individual example is nonrestrictive, have also contemplated that other tubular structures here and/or can quantify other spies
Property.
It would be recognized that assessment processor 122 can be a part and/or other calculating of control station 116
System, the most one or more computers.Suitable calculating system includes one or more processor,
It performs one or more computer-readable instructions, at computer-readable recording medium, such as local,
Portable or remote physical memorizer encodes or realizes instructions.Additionally or alternatively, Ke Yiyou
The carrier of such as signal or ripple carries one or more instructions.One or more processors are performing one
Or multiple seasonal one or more processors of instruction perform function and/or other merits of assessment processor 122
Energy.
Fig. 2 provides the non-limiting example of assessment processor 122.
In the illustrated embodiment, assessment processor 122 from reconstructor 114, data base and/or other
Device obtains, receives and/or retrieve 3D rendering data.Assessment processor 122 also from control station 116,
User and/or other resources obtain, receives and/or retrieve various processing parameter, such as voxel segmentation threshold,
Voxel segmentation identification information, assessment beginning and/or end position etc..The output bag of assessment processor 122
Include the signal of instruction assessment.
Dispenser 200 passes through automatization, semi-automation and/or manual type by the voxel in view data
It is divided into tubular structure interested and other structures.Such as, in one case, dispenser 200
Use threshold division method.Utilize a kind of thresholding method, it is possible to use predetermined voxel intensities threshold value
Being tubular structure interested or other structures by the voxel classification of volumetric image data, such as, mode is:
Intensity is identified as tubular structure interested (or other structures) intensity is low higher than the voxel of threshold value
Voxel in threshold value is identified as other structures (tubular structure interested).Can additionally use one or
Other threshold values multiple.Use two threshold values to allow intensity voxel classification between two threshold values to be
Tubular structure interested, is other structures by every other voxel classification.
Utilizing another kind of threshold technique, user provides the input of instruction structures of interest and indicates other
The input of structure, dispenser 200 determines suitable threshold value.Such as, user can use mouse etc. also
The part clicking on structures of interest represents the voxel of structures of interest to identify and clicks on other structures
A part represents the voxel of structures of interest to identify.The respective intensities value of identified voxel can be asked
Meansigma methods or be otherwise used for determining threshold value.Other dividing methods can also be used.Such as, exist
In the case of another kind, such as, by various partition tools, user can be with artificial cognition structures of interest
With other structures, these partition tools allow user by free-hand sketching the contours instrument, can redefine size
Prespecified geometric etc. identify ad hoc structure.
The data of segmentation are formatted by data formatter 202 enters one for assessment processor 122
Step processes.In one case, this can include that the view data resampling to segmentation is to change structure
Become the slice thickness of the individual section of 3D volume.In still another case, this can include segmentation
View data resampling, to produce the section on other directions, such as sagittal, crown, tilt and/
Or other sections.In the data that assessment processor 122 is split based on the section format analysis processing carrying out self-scanning
Time, it is convenient to omit data formatter 202.The concrete formatting used can be that acquiescence, user refer to
Fixed, person under inspection is special and/or extended formatting.
Distance map maker 204 produces for the tubular structure interested of one or more sections of 3D volume
The signal of raw instruction distance map.In one non-limiting case, suitable algorithm includes (such as base
In segmentation) identify the voxel areas representing tubular structure interested, determine and represent the one of structures of interest
Individual or multiple voxel to the beeline of the nearest voxel representing different structure, and by with short distance
From proportional intensity level distribution voxel.Profit in this way, represents the voxel of tubular structure interested
With represent different structure nearest voxel distance the most remote, relative intensity value will be the highest (or the lowest).Distance
Diagram generator 204 can use one or more algorithm with generation figure, such as, set from nomography storage
One or more drafting algorithms 206 of standby 208.
Process volume evaluator 210 to identify and produce in instruction formatting to be processed segmentation view data
The volume of tubular structure interested.In one non-limiting case, volume evaluator 210 base is processed
Volume in one or more input mark structures of interest of instruction volume.Such as, in a kind of situation
Under, user can use mouse etc. to split the avatars of image (such as via by prison from formatting
The graphic user interface that visual organ, display etc. present) identify structures of interest beginning and end.With
Family can also identify the path between beginning and end and/or the one or more point-of-interests along path.
In another case, user can identify a point of structures of interest, it is possible to use about
The preset range (such as, ± 5mm) of point is to determine volume.Process volume evaluator 210 can adopt
With one or more algorithm identification volumes, such as from one or many of volume algorithm storage device 214
Individual volume selection algorithm 212.In alternative embodiments, use process volume evaluator 210 produce away from
Before figure, identify and produce the signal of the volume indicating tubular structure interested, next use distance
Diagram generator 204 produces distance map for the one or more sections covering identified volume.
Quantify parts 216 to cut into slices from the 3D segmentation number covering volume based on the corresponding 2D from distance map
According to the various characteristics of one or more slice quantization identified volume of tubular structure interested.Quantization unit
Part 216 can use one or more such information of algorithm identification, such as, store from quantization algorithm
One or more quantization algorithms 218 of equipment 220.
In one non-limiting embodiment, the pre-of such as elastic ball is moved by tubular structure interested
Determine geometry, at the one or more points along tubular structure interested, determine the Local Metric of ball
Scope, such as radius or diameter, wherein by tubular structure interested around while ball at a point
Place grows the biggest.In one case, moving ball, thus along subrange moving ball
Center.In another case, along subrange otherwise moving ball.Can be determined by
The 2D section maximum of (wherein cutting into slices by this point) value of distance map or be determined by distance map
The maximum (wherein section is orthogonal with the longitudinal axis of tubular structure) of 2D section value determines the radius of ball.
Scope can be defined as the tubular structure interested face in the section orthogonal with path direction at this point
Long-pending.
The letter that processing component 222 can be produced by reason quantization unit part 216 at various modes and be exported
Number.In the illustrated embodiment, processing component 222 includes measuring determiner 224, and tolerance determiner is joined
It is set to determine one or more state (such as healthy) tolerance based on from the signal quantifying parts 216.
The suitably example of tolerance includes the tolerance indicating percent blockage.Can based on equation 1 determine one this
The tolerance of sample:
Equation 1:
Wherein null value (0) represents substantially do not have obturation ..., the value (50) of 50 represents that half is inaccessible ...,
The value (100) of 100 represents substantially all or entirely shuts.
In equation 1, minimum diameter can represent 1) practical minimum diameter of diameter value, 2) minimum
The meansigma methods of diameter range, 3) selected along the user in path or the diameter of pre-position, or 4) its
His diameter.Average diameter can represent 1) meansigma methods of all diameters, 2) think and represent tubular structure
Diameter meansigma methods within the diameter range of normal condition, 3) selected along the user in path or preset range
On average diameter, 4) maximum gauge in diameter, or 4) other diameters.
In modification, can with maximum gauge replace the minimum diameter in equation 1 to determine, such as with
The characteristic that the broadening of tubular structure interested is corresponding.One or more method identification can also be used not bery
It is suitable for the diameter value in equation 1.Such as, it is connected with other structures at tubular structure interested
Region, the diameter of tubular structure can increase.This region is used when determining the parameter of equation 1
Value may cause the average diameter of mistake.Identify such value allow to ignore, abandon, correct and/or with
Other modes process these values.
The processing component 222 of diagram also includes mapping means 226, and mapping means produces and indicates identified body
The correspondence of the diameter value of long-pending tubular structure interested and the view data of the identified volume of covering cuts into slices it
Between mapping.In one case, as illustrated signal can be rendered as curve 300.
In figure 3, such as, y-axis 302 represents that diameter value, x-axis 304 represent slice number, wherein cut into slices
From the beginning of numeral one (1), it is incremented by successively, until the last section covered in the section of identified volume.
Section index in y-axis can be for the section orthogonal with tubular structure direction forms direction.
In figure 3, region 306 represents the diameter value around the non-congested areas of tubular structure interested,
Region 308 represents the diameter value around tubular structure congested areas interested.Other districts of curve 300
Territory represents the diameter value of transitional region, is not likely to be and is suitable for being included to determine tubulose interested
The diameter of structure or the minimum diameter of tubular structure interested.Diameter value can be obtained from curve 300,
These values can be used to determine percentage blocks or other information, and/or identify that the correspondence of view data is cut
Sheet.Signal acquisition can also be utilized such and/or other information and do not show curve 300.
Fig. 4 shows the method for the characteristic for quantifying tubular structure interested.
402, it is thus achieved that include the body image of structures of interest and other structures.
404, volumetric image data is divided into structures of interest and other structures.
406, as described herein, three-dimensional distance figure is produced for volumetric image data.
408, identify tubular structure volume interested to be quantified from the volumetric image data of segmentation.
410, characteristic based on the tubular structure interested on distance map the identified volume of quantization.
Fig. 5 shows the non-limiting method producing distance map, such as, for the action 406 of Fig. 4.
502, the voxel of the two dimension slicing from segmentation view data is identified as tubular structure interested
Or other structures.
If this voxel is identified as tubular structure interested, then 504, determine from this voxel to table
Show the distance of the nearest voxel of other structures.
506, produce and store to this apart from the intensity level of proportional voxel.
508, it is determined whether there is another voxel to assess.
If it were to be so, then palikinesia 502.
If not, then 510, produce the signal of instruction intensity level and the 2D as distance map cuts
Sheet exports.
If being identified as other structures at 502 voxels, then execution action 508.
512, it is determined whether there is another section to assess.
If it were to be so, so for this palikinesia 502 of cutting into slices.
If not, then as action 514, combined distance figure is cut into slices, is consequently formed 3D distance map.
Fig. 6 shows the non-limiting method of the tube chamber characteristic for quantifying tubular structure interested, example
As, as the action 410 of Fig. 4.
602, it is thus achieved that the 2D section of the distance map corresponding with the point along tubular structure interested.
604, in 2D cuts into slices, identify tubular structure interested.
606, geometry is made to grow in the tubular structure interested identified, until this shape is covered
Lid section is connected in 602 all voxels of the point selected.
608, identify and store the maximum of voxel intensity value of this shape 2D section.As it was previously stated,
Intensity level is corresponding to the distance from the voxel of the structures of interest nearest voxel to other structures.
610, for other palikinesia 602-608 one or more along tubular structure interested.
612, produce based on maximum intensity value and indicate the letter of the quantization of tubular structure interested on these aspects
Number.Such as, this maximum and minimum of intensity level maximum that can include identifying intensity level maximum
Value, and determine the ratio of maxima and minima, or deduct this ratio with one.
It would be recognized that the sequence of action is not to limit.So, in other embodiments, action
Sequence can be different.Furthermore, it is possible to omit one or more action and/or one or many can be increased
Other actions individual.
Above-mentioned action can be implemented by computer-readable instruction, when being performed by computer processor,
Computer-readable instruction performs action described here by processor.In this case, instruction storage
In a computer-readable storage medium, such as associate with correlation computer and/or otherwise can be by
In the memorizer that correlation computer accesses.
Foregoing can be applied to study tubular structure, such as blood vessel, colon etc., including, but do not limit
In quantification of stenosis and produce narrow index on CT, MRI or ultrasonic three-dimensional data set, quantify aneurysm,
Quantify polyp and/or other possible other structures blocking tubular structure.
As non-limiting example, it is possible to use system 100 is to one or more blood vessels (the most left neck
Tremulous pulse) carry out the research of contrast agent enhanced, it is used for characterizing narrow.The view data of self-scanning in the future
Three-D volumes is divided into blood vessel and non-vascular.In one embodiment, can be by internal blood vessel point
Hit once, click on once in outside, and use the meansigma methods of two intensity of these points to carry out thresholding
Realize this purpose.
Computed range figure, it comprises the instruction voxel voxel intensities to the distance of nearest non-vascular voxel.
Distance map is 3-D data set, has same scale with raw data set.User provides interested blood vessel
Beginning or end at least one.Likewise it is possible to existed by the three-D volumes of divided data
Click in blood vessel starting point and on blood vessel terminal and realize this purpose.
Adjust the distance figure resampling to obtain the equidistant section orthogonal with the line between beginning and end.?
The voxel inside interested blood vessel with distance map maximum intensity is found in each section.This voxel exists
The center of maximum fitting sphere in this section.This diameter can be found in any 2D cuts into slices.Can
To present the process from origin-to-destination or subrange to user, such as, as shown in Figure 4.
Narrow diameter can be found from the minima of this curve.Can be by complete curve or narrow
Calculate maximum in neighbouring subset and find distal diameter index.Clinicist is allowed (such as to radiate above
Section doctor and/or scientist) accurately and reproducibly to calculate by interactive selection blood vessel beginning and end
Stenosis.This acceptable results including providing thorny position, such as blood vessel joint, or touch not
The blood vessel that can be separated by segmentation.
In a specific embodiment, it is possible to use this method determines narrow diameter, for contrast agent
The dual energy CT strengthened.This may be particularly useful, because double-energy scanning allows separating blood vessel
With calcium speckle, they all seem the brightest in contrast agent enhanced CT.Another kind of situation about using is to dynamic
Research incubation period of pulse atherosclerosis.
The present invention is described by reference to preferred embodiment.Other people are reading and are understanding discussed in detail above
Afterwards can be it is contemplated that modifications and changes.All such amendments and change should be invention is constructed as including
More, as long as they are within the scope of claims or its equivalents thereto.
Claims (12)
1. for the method quantifying the characteristic of the tube chamber of tubular structure interested, including:
Generating 3D voxel distance map, wherein, described 3D voxel distance map includes intensity level, described by force
Angle value instruction is from the nearest voxel of the voxel to other structures of expression representing described tubular structure interested
Distance；
Signal is produced, the institute represented in described signal designation imaging in vivo data based on 3D voxel distance map
State the described spy of the described tube chamber of described tubular structure interested on the predetermined of tubular structure interested
The quantization of property,
Wherein, described characteristic determine by elastic ball is moved through described tubular structure, edge
Radius at the one or more points of described tubular structure or diameter, and at the one or more point
Place determines described radius or diameter, grows as far as possible at the one or more every described ball in place put
Greatly,
Wherein, the voxel of the 2D section being determined by the described 3D voxel distance map corresponding with every is strong
The maximum of angle value and make the voxel with the maximum of described voxel intensity value as described elasticity
Ball center in described 2D cuts into slices determines described radius or diameter, and wherein, described 2D section is
By what described 3D voxel distance map resampling was obtained.
Method the most according to claim 1, also includes:
Described body imaging data is divided into one or more regions of tubular structure and the one of other structures
Individual or multiple regions；
One or more voxels for one or more regions of described tubular structure determine and other knots
The distance of the nearest voxel of structure；
One or more intensity level, wherein, the intensity of voxel is determined for the one or more voxel
Value is proportional to the distance that should determine that mutually；And
Described 3D voxel distance map is produced based on described intensity level.
Method the most according to claim 2, also includes:
The subregion in one or more regions of described tubular structure is identified as tubular structure interested；
Value based on described 3D voxel distance map determines the Local Metric model of described tubular structure interested
Enclose；And
The quantization of the characteristic of described tubular structure interested is determined according to described Local Metric scope.
Method the most according to claim 3, also includes:
Determine the beginning and end on described tubular structure interested；
Determine from described starting point to the path of described terminal；
Determine the multiple points along described path；And
Determine the described Local Metric scope of described tubular structure interested at described point.
5., according to the method according to any one of claim 3 to 4, also include:
Described Local Metric scope for any is defined as described tubular structure interested with this point
Locate the area in the section that the direction in described path is orthogonal.
6. according to the method according to any one of claim 3 to 4, wherein, the quantization of described characteristic
Based on minimum subrange and maximum subrange.
7. according to the method according to any one of claim 2 to 4, wherein, instruction user select
The described tubular structure interested of input identification of beginning and end.
8. according to the method according to any one of claim 1 to 4, wherein, tubular structure interested
For blood vessel, quantified characteristic is diameter or the radius of described blood vessel.
Method the most according to claim 8, wherein, described diameter or radius indicate the narrow of blood vessel
Narrow or aneurysm.
Method the most according to claim 9, wherein, by the value of the multiple points along described path
It is mapped to the correspondence position of imaging data.
11. 1 kinds of systems being used for quantifying the characteristic of the tube chamber of tubular structure interested, including:
Distance map maker (204), it produces 3D voxel distance map, wherein, described 3D voxel away from
Including intensity level from figure, the instruction of described intensity level is from representing that the voxel of described tubular structure interested is to table
Show the distance of the nearest voxel of other structures；
Quantifying parts (216), it produces signal, described signal designation based on described 3D voxel distance map
Described tubulose interested on the predetermined of the tubular structure described interested represented in imaging in vivo data
The quantization of the described characteristic of the described tube chamber of structure,
Wherein, described characteristic determine by elastic ball is moved through described tubular structure, edge
Radius at the one or more points of described tubular structure or diameter, and at the one or more point
Place determines described radius or diameter, grows as far as possible at the one or more every described ball in place put
Greatly, wherein, the voxel of the 2D section being determined by the described 3D voxel distance map corresponding with every is strong
The maximum of angle value and make the voxel with the maximum of described voxel intensity value as described elasticity
Ball center in described 2D cuts into slices determines described radius or diameter, and wherein, described 2D section is
By what described 3D voxel distance map resampling was obtained.
12. systems according to claim 11, also include:
Dispenser (200), described body imaging data is divided into tubular structure interested and other knots by it
Structure.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US38434610P | 2010-09-20 | 2010-09-20 | |
| US61/384,346 | 2010-09-20 | ||
| PCT/IB2011/054034 WO2012038863A1 (en) | 2010-09-20 | 2011-09-15 | Quantification of a characteristic of a lumen of a tubular structure |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103109308A CN103109308A (en) | 2013-05-15 |
| CN103109308B true CN103109308B (en) | 2016-11-30 |
Family
ID=
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101203180A (en) * | 2005-06-21 | 2008-06-18 | 皇家飞利浦电子股份有限公司 | Method and device for imaging a blood vessel |
| CN101493940A (en) * | 2009-02-27 | 2009-07-29 | 东南大学 | Vascular parameter measurement method in medicine image based on elastic pellet dynamic balance |
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101203180A (en) * | 2005-06-21 | 2008-06-18 | 皇家飞利浦电子股份有限公司 | Method and device for imaging a blood vessel |
| CN101493940A (en) * | 2009-02-27 | 2009-07-29 | 东南大学 | Vascular parameter measurement method in medicine image based on elastic pellet dynamic balance |
Non-Patent Citations (2)
| Title |
|---|
| Vessels as 4-D Curves: Global Minimal 4-D Paths to Extract 3-D Tubular Surfaces and Centerlines;Hua Li等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20070904;第26卷(第9期);第1213-1220页 * |
| 凸阵探头对颈动脉狭窄病变的检测价值;华扬等;《中国超声医学杂志》;20040116;第20卷(第1期);第17-19页 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102269467B1 (en) | Measurement point determination in medical diagnostic imaging | |
| US6545678B1 (en) | Methods, systems, and computer program products for generating tissue surfaces from volumetric data thereof using boundary traces | |
| CN106456125B (en) | System for linking features in medical images to anatomical model and method of operation thereof | |
| US8858436B2 (en) | Systems and methods to identify interventional instruments | |
| CN100528078C (en) | Method and system for airway measurement | |
| CN101002205B (en) | Data set visualization | |
| EP2377095B1 (en) | Vessel analysis | |
| CN101410060A (en) | Determining tissue surrounding an object being inserted into a patient | |
| JP6496403B2 (en) | A processing system configured to cooperate with an optical shape sensitive interventional device | |
| JP5295562B2 (en) | Flexible 3D rotational angiography-computed tomography fusion method | |
| CN105392428A (en) | System and method for mapping ultrasound shear wave elastography measurements | |
| JP2006255412A (en) | Method and system for monitoring tumor burden | |
| CN103300891B (en) | Combined display method and system for quantitative elasticity information and structure information of tissues | |
| CN101681514A (en) | Inspection of tubular-shaped structures | |
| CN101902964A (en) | Hardware tumor phantom for improved computer-aided diagnosis | |
| US20060104495A1 (en) | Method and system for local visualization for tubular structures | |
| Huang et al. | 2.5-D extended field-of-view ultrasound | |
| CN105723423B (en) | Volumetric image data visualization | |
| US20240193764A1 (en) | Systems and methods for reconstruction of 3d images from ultrasound and camera images | |
| CN102118999B (en) | Determining foreshortening optimal view maps taking into account the shape of the vessel cross-section in cardiovascular X-ray systems | |
| EP2619729B1 (en) | Quantification of a characteristic of a lumen of a tubular structure | |
| CN106133797A (en) | There is the medical viewing system checking that plane determines | |
| US20160180520A1 (en) | Quantitative method for 3-d joint characterization | |
| US20120078101A1 (en) | Ultrasound system for displaying slice of object and method thereof | |
| US8724878B2 (en) | Ultrasound image segmentation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |