CN102246206A - Vessel analysis - Google Patents

Vessel analysis Download PDF

Info

Publication number
CN102246206A
CN102246206A CN2009801493677A CN200980149367A CN102246206A CN 102246206 A CN102246206 A CN 102246206A CN 2009801493677 A CN2009801493677 A CN 2009801493677A CN 200980149367 A CN200980149367 A CN 200980149367A CN 102246206 A CN102246206 A CN 102246206A
Authority
CN
China
Prior art keywords
tubular structure
blood vessel
module
pointed position
branch
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2009801493677A
Other languages
Chinese (zh)
Other versions
CN102246206B (en
Inventor
J·J·松内曼斯
R·J·E·哈比斯
J·奥利万贝斯科斯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN102246206A publication Critical patent/CN102246206A/en
Application granted granted Critical
Publication of CN102246206B publication Critical patent/CN102246206B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46—Arrangements for interfacing with the operator or the patient
    • A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/08—Volume rendering
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00—Manipulating 3D models or images for computer graphics
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/24—Indexing scheme for image data processing or generation, in general involving graphical user interfaces [GUIs]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/20—Special algorithmic details
    • G06T2207/20092—Interactive image processing based on input by user
    • G06T2207/20101—Interactive definition of point of interest, landmark or seed
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • G06T2207/30101—Blood vessel; Artery; Vein; Vascular
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/004—Annotating, labelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/012—Dimensioning, tolerancing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Graphics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Software Systems (AREA)
  • Pathology (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Human Computer Interaction (AREA)
  • Quality & Reliability (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Computer Hardware Design (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

A system for performing vessel analysis uses display means (1) for displaying a three-dimensional image representing at least a tubular structure. Indicating means (2) are used for enabling a user to indicate a position on a vessel of the tubular structure, for obtaining an indicated position. Identifying means (3) are used for identifying a portion of the tubular structure situated around the indicated position, including any bifurcations, and extending up to a predetermined distance measure from the indicated position, for obtaining an identified portion. The display means (1) is also used for displaying a graphical annotation in the displayed three-dimensional image, indicative of the identified portion of the tubular structure.

Description

Vessels analysis
Technical field
The present invention relates in tubular construction the analysis of one or more blood vessel in tubular structure, the especially medical image.The invention still further relates to user interactions with the vessels analysis instrument.
Background technology
Medical images data sets comprises the medical image of vascular structure.These medical images can be two dimensions (such as in vascular x ray), or three-dimensional (such as in CT angiogram or MR angiogram).3-D view can be a volumetric image, and it means with gray-scale value represents vascular structure.Existing blood vessel tracer tools can be based on the each several part of these gray-scale value identification vascular structures.
For example, US 2008/0101667A1 discloses a kind of being used for and has presented the information relevant with blood vessel to assess the device of described blood vessel to the user.This device comprises storer and processor, and described storer can be stored in the imaging data that limits in the three-dimensional, and described processor operationally is coupled to described storer to receive imaging data.Processor is configured to present to the user two dimension slicing of the three-dimensional imaging data of blood vessel in first user interface; Receive blood vessel from the user and select, wherein, the user by with the mutual selection blood vessel of first user interface; The blood vessel path that identification is associated with the blood vessel selection that receives from three-dimensional imaging data; And present intensity along the selected blood vessel in the blood vessel path discerned for selected blood vessel is analyzed to the user.
Summary of the invention
Have that a kind of what be used to carry out vessels analysis will be comparatively favourable through improved system.In order to address this problem better, in a first aspect of the present invention, a kind of system has been proposed, this system comprises
Display module, it is used to show the 3-D view of representing tubular structure at least;
Indicating module, it is used to make the user can point out a position on blood vessel in tubular construction, with the position that obtains to point out; And
Identification module, it is used to discern the part that is positioned at pointed position tubular structure on every side, and described part comprises branch arbitrarily, and extends predetermined distance from pointed position, with the part that obtains to be discerned.
The user only needs to point out a single point on vascular structure.After this, system is can predict and consistent mode is discerned the part of tubular structure.Because the part of being discerned comprises complete tubular structure, comprises branch arbitrarily, so there is not the doubt that will discern which bar blood vessel.This has increased the confidence of user when Location.Same, because be on predetermined distance, to discern blood vessel, not it may be noted that any end points of vessel segment to be analyzed.This system provides the effective tool that is used for the part of vascular structure is carried out vessels analysis.
Display module can be arranged to be used for the 3-D view display graphics note shown, the part of being discerned of its indication tubular structure.This provides feedback to institute's identification division for the user, makes the user can examine the part of being discerned.For example, can point out another point after the user on vascular structure, identification module can be discerned the correspondence position of tubular structure then.Discern the part of tubular structure in reproducible mode, described mode makes shown graphical annotations also can reproduce.Such note also can be stored in patient's archives.
Identification module can comprise the measurement module that is used for measuring along tubular structure preset distance.This is the suitable manner in conjunction with the shape measure preset distance of blood vessel.The part of being discerned can comprise the blood vessel that extends to the endpoint location of blood vessel from pointed position, and the length of vessel from pointed position to endpoint location is substantially equal to predetermined distance.Alternatively, can use Euclid (Euclidean) distance.
Measurement module can be arranged to be used to measure the length of the blood vessel of following the trail of from pointed position along tubular structure.The consistent length that keeps the blood vessel of following the trail of is to make blood vessel follow the trail of consistent effective means.
Measurement module can be arranged to be used for determining some branches of the experience from pointed position along tubular structure.This is to make blood vessel follow the trail of consistent another kind of mode.
The reorganization recombination module can be arranged to be used to provide many curved surface reorganization reorganization (multi-curved reformat), its follow from pointed position to tubular structure branch and away from the curvature of at least a portion of the blood vessel of the branch of tubular structure, and follow curvature away from least a portion of two blood vessels of branch.This visualization model provides the useful view of branch.
Indicating module can be arranged to be used to respond the mouse moving event.The mouse moving event makes and can discern very effectively some part of vascular structure by mouse pointer being moved to top, interested blood vessel position.Whenever detecting the mouse moving event, just can activate identification module and/or display module so that the position of the identification through upgrading to be provided.
Described system can comprise estimation module, and it is used for estimating local trace parameters in pointed position according to volumetric image after the Location the user.This allows to utilize local trace parameters initialization blood vessel tracing algorithm.This initialization step can use the blood vessel tracer tools with large-scale blood vessel.
Estimation module can be arranged to be used for the weighted mean of the key element of computation structure tensor, calculates described weighted mean around pointed position above the zone of intended size.Such weighted mean of the key element of structure tensor can be used as the key element through average structure tensor.Can be used to determine orientation through average structure tensor at pointed position blood vessel.
A kind of medical imaging workstation can comprise said system.
A kind of method of carrying out vessels analysis can comprise:
Show the 3-D view of representing tubular structure at least;
Make the user can on blood vessel in tubular construction, point out a position, to obtain pointed position; And
Identification is positioned at the part of the tubular structure around the pointed position, and described part comprises branch arbitrarily, and extends predetermined apart from size, with the part that obtains to be discerned from pointed position.
A kind of computer program can comprise being used to make processor system to carry out the instruction of said method.
It should be recognized by those skilled in the art that can be to think that useful any-mode makes up two or more in the foregoing description, enforcement and/or the each side of the present invention.
Those skilled in the art can realize the modification and the modification of image capture device, workstation, system and/or the computer program corresponding with the described modification of system and modification based on this explanation.
Those skilled in the art will be appreciated that, described method can be applied to the multidimensional image data, for example be applied to 2 dimensions (2-D), 3 dimensions (3-D) or 4 dimension (4-D) images, described image can be gathered by various collection apparatuses, such as, but be not limited to standard X-ray imaging, computer tomography (CT), magnetic resonance imaging (MRI), ultrasonic (US), PET (positron emission tomography) (PET), single photon emission computed tomography (SPECT) and nuclear medicine (NM).
Description of drawings
To further illustrate and describe these and other aspects of the present invention with reference to the accompanying drawings.
Fig. 1 shows the block scheme of the system that is used to carry out vessels analysis; And
Fig. 2 shows the block scheme of the method for carrying out vessels analysis.
Embodiment
In some vascular applications, one of target is to use curved surface or as seen the reorganization repopulate views that stretches makes blood vessel.Similarly, the local vascular parameter such as area and radius is measured in some positions that for example can be in view data, thus the degree of quantification of stenosis or aneurysmal size.These visualization techniques can use by angiocentric path.The automated path tracer tools usually only is confined to several limited blood vessels.Semi-automatic path tracing instrument may be followed the trail of or two points to point out the single-point on the blood vessel from user's input.When only using single input point, the length of the uncontrollable blood vessel of following the trail of the path or being followed the trail of of user.The starting point that the use user provides and the instrument of terminal point have following shortcoming: promptly they need more user's input.Similarly, if the user is placed on distance position far away excessively or be placed in the different blood vessels each other with point, then or route result be wrong, perhaps computing time is long.Some blood vessel tracer toolses need at first be cut apart blood vessel structure before the center line of following the trail of by blood vessel structure.But this energy efficiency is not high.
Can use provides the tracer tools of distinct interaction possibility to improve interactivity.Use visualization technique arbitrarily, the user can point to optional position on the vascular structure (for example by rolling mouse arrow above image), and this algorithm can be automatically and alternatively follow the trail of local tree.This tracking can be carried out the length of vessel that reaches maximum.Exchange method described herein makes it possible to limit again the example of vascular research.This algorithm can be that the user provides the preview view of selection particular blood vessel with the effect of appearance in the mode of " will being what situation if I Click here ".Exchange method described herein can the various visual and quantification techniques of direct correlation.
By guaranteeing to improve interactivity to the tracking of vasculature part with to the operation fast that is visualized as of tracking vasculature part.Owing to only follow the trail of the finite part of vascular tree, it is faster that ground is carried out in feasible tracking.Can be at an arbitrary position with limited mutual execution based on the visual of path or tree with measure.By visual view feedback is provided, shown direct response to user action, this makes described instrument be more readily understood and grasp.The behavior of instrument may have more foresight to the user because the length of the vasculature part of following the trail of may with the location independent of mouse.This has provided when blood vessel moves at mouse and can predict and intelligible result.Technology described herein can for example be applied to the application at multi-modal 3D vessel analysis.These instruments provide senior observation to blood vessel, have cut apart, study and quantize.
Term 3D vascular quantizes to relate to the collection to using, and it uses different acquisition methods, is target with different vascular structures, but may be similar to the demand of expection measurement.The anatomical structure example is main artery, arteria carotis, coronary artery, peripheral leg artery and coronary artery.MR and CT and rotational x-ray are the examples that is used to gather the view data of described anatomical structure.The example of vascular research will be to seek blood vessel to broaden or blocked part.More specifically example will be to search pulmonary embolism in pulmonary artery.Vascular applications can be used the senior visualization technique as research purpose, such as curved surface MPR, makes visualization of blood vessels.These application can also be in view data some positions measure local vascular parameter such as blood vessel area and radius, thereby the degree of quantification of stenosis or aneurysmal size.Can pass through to carry out these measurements on the cross section of interested blood vessel, described cross section can automatically generate, or generates by the navigation along vessel centerline.Other application may concentrate on the surgery planning of getting involved such as the vascular of stent planning.
Fig. 1 shows the block scheme of the system that is used to carry out vessels analysis.This system can be made up by some chunks.Each chunk can partially or completely be realized with software.This system also comprises (not shown): processor, memory module, display, such as user's load module, the communication port of keyboard and/or mouse.Processor can be arranged to be used to carry out the instruction as the part of vessels analysis instrument.Memory module can comprise RAM, ROM and hard disk, such as the removable media of CD and DVD.Memory module can be used for the instruction of storage computation machine and/or be used to store medical image.Load module can be used to make the user can point out a position on tubular structure.In addition, load module can be used to make the user can control described system.Communication port can be used in another computer system and communicates by letter, for example server.Communication port can be arranged to be used to be connected to network, such as LAN (Local Area Network), wide area network and/or the Internet.For example, can visit other computer systems with retrieve image data, and transmit the vessels analysis report via described network.
Can provide display module 1 to show the 3-D view of representing tubular structure at least.This 3-D view obtains and can be stored in the memory module via communication port.Can use different types of visualization technique to show described image.For example, can show that direct volume presents, maximum intensity projection or slice view.Can in proper order or adopt simultaneously multiple presenting.Described system can allow image scaled, move and/or rotation so that obtain the suitable view of image.Described image can be represented tubular structure, and for example described image comprises the volume element (voxel) of indication at the gray-scale value at the concrete coordinate position place of volume.Can and discern blood vessel structure by their gray-scale value in some cases by the concrete shape of blood vessel structure.For example, 3-D view can also be expressed as surface model with vascular structure.
Can provide indicating module 2 so that the user can point out a position on blood vessel in tubular construction, with the position that obtains to point out.This indicating module can by make the user can the position in shown image on the location mouse pointer realize.Based on the position of mouse pointer, can calculate position on the blood vessel by mode known in this area.
Can provide identification module 3 to be positioned at the part of the tubular structure around the position discerned with identification, described part comprises branch arbitrarily, and extend from pointed position predetermined apart from size, thereby obtain the part discerned.For example, identification module 3 can be arranged to respond effectively the mouse moving event.
Display module 1 can be arranged to be used for being presented at the graphical annotations of shown 3-D view, the part of being discerned of its indication tubular structure.For example, illustrate vessel centerline with differentiable color.
Identification module 3 can comprise measurement module 4, is used to measure the preset distance size along tubular structure.For example, wave front (wavefront) algorithm is applied to following the trail of blood vessel, and each new ripple adds predetermined amount for the distance of measuring.Identification module 3 identifications reach the tubular structure of preset distance.A kind of alternatives is to measure Euclidean distance.Measurement module 4 can be arranged to be used to measure the length of the blood vessel of following the trail of from pointed position along described tubular structure.For example, the wave front algorithm is applied to following the trail of blood vessel, and each new ripple adds predetermined amount for the distance of measuring.Identification module 3 identifications reach the tubular structure of preset distance.Alternatively, measurement module 4 can be arranged to be used for determining from pointed position along the number of the branch that described tubular structure meets with.For example, according to the described blood vessel of pointed location recognition that arrives second branch.
Reorganization recombination module 5 can provide many curved surface reorganization (many curved surfaces MPR).The curvature of reorganization is corresponding to reaching from pointed position and away from the curvature of at least a portion of the blood vessel of the branch of tubular structure, and follows the curvature away from least a portion of two blood vessels of described branch.Can on display 1, show described reorganization.
Indicating module 2 can be arranged to be used to respond the mouse moving event.This is very easily, particularly when enough carrying out the blood vessel tracking apace.Can probe into blood vessel and need not clearly to click certain location.Except that of the response of described system, arrive from pointed position apart from the identification of the blood vessel structure of preset distance and the availability that note has also increased system to the mouse moving event.
3-D view can comprise volumetric image, and identification module 3 also comprises estimation module 6, is used for estimating local trace parameters in response to user's Location according to the volumetric image of pointed position.By directly estimating local trace parameters according to volumetric image (that is, the gray-scale value of voxel) in response to user's Location, it may be suitable for different vascular regions, and this allow to use described system and need not at first to select the blood vessel of concrete kind.
Estimation module 6 can be arranged to be used for the weighted mean of the key element of computation structure tensor, calculates described weighted mean around pointed position on the zone of intended size.This is effective, the universal mode of initialization tracing algorithm, and it is applicable to the vascular structure of numerous kinds.Can calculate this tensor fast.Be necessary to follow the trail of in advance vascular structure, this operation can be finished after the user has pointed out described position.
The system of Fig. 1 can be integrated into medical imaging workstation or be integrated in the control desk of medical imaging devices.It will be appreciated by those skilled in the art that medical imaging workstation and control desk comprise a processing unit and storer at least.It can also comprise many other unit, for example hard drive, RAM, buffer memory, display, data and instruction bus, display and keyboard.
Fig. 2 illustrates a kind of method of carrying out vessels analysis.This method comprises, in step 101, shows the 3-D view of representing tubular structure at least.In step 102, make the user can on blood vessel in tubular construction, point out a position, with the position that obtains to point out.In step 103, be identified in the part of pointed position tubular structure on every side.Described part can comprise branch arbitrarily and can extend predetermined distance scale from pointed position, with the part that obtains to be discerned.
For blood vessel is studied, the user can utilize the visualization technique that does not rely on path, geometry to make data visualization, and all many in this way plane reorganization of described technology (MPR), maximum intensity projection (MIP) or volume are drawn (VR).For for the visualization technique in path, need by angiocentric path.
A kind of algorithm that is used to follow the trail of blood vessel can be provided.Here, " tree " speech refers at least a portion of vascular structure, comprises branch and side shoot.Such algorithm can find the local tree topological structure that starts from the 3D optional position in the data (for example, user-selected position).In described tree topology structure, can encode to the relation between the position of coordinate, blood vessel and every branch of every blood vessel.Tree is extracted the restriction that may be subjected to maximum blood vessel trace lengths.This length can be along the vessel measurement of being followed the trail of but not along line measurement.
Mouse position on any visable representation of 3D medical image can be converted into { x, y, z} Data Position.For example, can use and corresponding { x, y, the z} position of visible blood vessel, mouse pointer below.Can use centering (centering) algorithm computation local vascular central point, and can use maximum length of vessel standard from this location tracking tree.Can use the tree followed the trail of to create based on tree or based on the visable representation in path and calculate measurement result.If the user is satisfied to the result, he can " stopping "/final definite result of tracking.After this, the user can determine to edit following the trail of the result.In case the user is satisfied, resulting image and corresponding measurement result can be sent to Reporting Tools.
In implementing example, " tree probe " can utilize on the function on the task panel by button activation.Activating " tree probe " afterwards, by rolling mouse above image automatically with the described tree of given track in range.Based on the local tree direction, aim at orthogonal view (orthoview) at the mouse position place.For example can be by using different color each branches of mark in image.Except that orthogonal view, can also create and show curved surface reorganization (in the position of vertical view).In addition, according to three vessel segments that depart from described branch, the view of many plane reorganization (MPR) can be aimed at immediate branch.
A kind of method can be provided, it utilizes the mutual fully and maximum repeatability in any visable representation of 3D medical image, and need not any click, follow the trail of the local tree topological structure from the mouse position of user's qualification and along given maximum length of vessel, wherein, the extracting position of each branch and blood vessel is known.Described tree topology structure can be used for alternatively display of visuallyization and represent that described visable representation presents etc. such as curved surface MPR, many curved surfaces MPR, the reorganization of stretching, orthogonal view, highlighted volume based on the tree topology structure.The vessel boundary detection method can be used for calculating quantization parameter along the local tree topological structure, and such as area, diameter, radius, and these parameters that are used to all can be obtained on the view are visual.Can propose to be used to quantize and at the automatic measuring position of other diseases based on the parameter of automatic derivation.Can only utilize a click just will send to Reporting Tools such as the information of image and measurement result.
The blood vessel tracing system can be carried out initialization by the preset value of choosing especially at concrete application.For example, preset value is specific for vascular regions (for example ' peripheral blood vessel ' or ' coronary vasodilator ') and/or image mode (for example CT).
Yet, also may add local parameter algorithm for estimating (preferably in real time), described algorithm carries out initialization to trace parameters (and randomly, upgrading described parameter during following the trail of).This allows to observe and need not to influence supposition to concrete blood vessel size.It allows tracer tools to follow the trail of the blood vessel of range of size broad.Similarly, the initialization of trace parameters allows to carry out on the MR image of intensity level the unknown of blood vessel before and observes (and may have a large amount of uneven pseudo-shadows).Provide the 3D source data volume of vascular structure, the user can select the visualization scheme that presents as MIP (maximum intensity projection), MPR (reorganization of many planes) or volume.For Location on blood vessel structure, the user can observe blood vessel by rolling mouse on image.Can take following steps blood vessel is carried out initialization: coordinate { x, y, the z} of pointed position in the retrieval 3D volume source in pointed position; Calculating local vascular orientation.Similarly, according to described local vascular orientation, can calculate the blood vessel profile.
The local vascular orientation can find by determining the picture structure orientation; The latter can directly calculate according to topography's GTG value by the utilization structure tensor.Can provide structure tensor by following formula
J = < g x g x > < g x g y > < g x g z > < g y g x > < g y g y > < g y g z > < g z g x > < g z g y > < g z g z >
Here, g iIt is image gradient along direction i.Parantheses＜〉be illustrated in the weighted mean on the zone of intended size.Can use Gauss (Gaussian) to blur and carry out the weighting operation.After the computation structure tensor, can computation of characteristic values { λ 0, λ 1,, λ 2And proper vector { v 0, v 1, v 2.Can be according to constraint condition λ 0≤ λ 1≤ λ 2Arrange eigenwert.Therefore, v 0Direction corresponding to the weighted product minimum of gradient.In tubular structure, this is corresponding to the local vascular direction.v 1And v 2Leap is perpendicular to the transversal plane of blood vessel.
Can determine the local vascular direction by the utilization structure tensor.Can determine estimation perpendicular to the intravascular space cross section of this vessel directions.For example, can determine the circular section, perhaps can carry out more detailed profile and detect.Can compare within the vascular cross-section and vascular cross-section outside image intensity.Such comparison can be used for initialization or upgrade the parameter of tracing algorithm.
Will be appreciated that the present invention also expands to computer program, specifically expand to the computer program on carrier or in the carrier, it is suitable for the present invention is tried out.Described program can be the form of source code, object code media supply and object code, such as the form of part compiling, or is suitable for any other forms of using in enforcement according to the present invention.It should also be appreciated that such program can have many different architecture designs.For example, the code of enforcement the method according to this invention or systemic-function can be subdivided into one or more subprocesss.Many different modes of distribution function are apparent to those skilled in the art between these subprocesss.Described subprocess can be stored together the program of self-contained to form (self-contained) with the form that can carry out document.The document carried out like this can comprise computer executable instructions, for example processor instruction and/or interpretive order (for example Java interpretive order).Alternatively, in the subprocess one, a plurality of or all be stored at least one external libraries document, and for example when operation, link statically or dynamically with master routine.Master routine comprises at least one at least once call in the subprocess.Similarly, subprocess can comprise the funcall to each other.An embodiment who relates to computer program, comprise with preceding method at least one treatment step in each corresponding computer executable instructions.These instructions can be subdivided into subprocess and/or be stored in can static state or one or more documents of dynamic link in.Another embodiment that relates to computer program comprise with aforementioned system and/or product at least one module in each corresponding computer executable instructions.These instructions can be subdivided into subprocess and/or be stored in can static state and/or one or more documents of dynamic link in.
The carrier of computer program can be any entity or the device of the program of can carrying.For example.Described carrier can comprise the storage medium such as ROM, for example CD ROM or semiconductor ROM; Or magnetic recording media, for example floppy disk or hard disk.In addition, described carrier can be the carrier that can transmit, and such as electricity or light signal, it can or pass through radio frequency or the transmission of other means via cable or optical cable.When described program is embedded in the sort signal, can be by cable or other devices or module appointment carrier.Alternatively, described carrier can be wherein to have embedded described procedure set to become circuit, and described integrated circuit is suitable for carrying out or is used to carry out relevant method.
Will be appreciated that the foregoing description is used for explanation but not is restriction the present invention, and those skilled in the art can design many alternatives under the situation that does not deviate from claim scope of the present invention.In the claims, place any reference marker between the parenthesis should not be interpreted as claim is construed as limiting." comprise " speech with and the use of combination do not get rid of element or step outside those elements described in the claim or the step.Indefinite article " one " or " one " do not get rid of the existence of a plurality of these class components.The present invention can implement by the hardware mode that comprises some discrete components, and the mode of the computing machine by suitable programming is implemented.Enumerated plurality of modules in the device claim, some in these modules can be by one in the hardware or identical content realization.Applied certain measures is not indicated the combination of these measures of use that can not be favourable in different mutually dependent claims.

Claims (12)

1. a system that is used to carry out vessels analysis comprises
Display module (1), it is used to show the 3-D view of representing tubular structure at least;
Indicating module (2), it is used to make the user can point out a position on blood vessel in tubular construction, to obtain pointed position; And
Identification module (3), it is used to discern the part that is positioned at pointed position described tubular structure on every side, and described part comprises branch arbitrarily, and extends predetermined distance scale from pointed position, with the part that obtains to be discerned.
2. system according to claim 1, described display module (1) is arranged to be used for the 3-D view display graphics note shown, and it indicates the part of being discerned of described tubular structure.
3. system according to claim 1, described identification module (3) comprises measurement module (4), it is used for measuring predetermined distance scale along described tubular structure.
4. system according to claim 3, described measurement module (4) is arranged to be used to measure the length of the blood vessel of following the trail of from pointed position along described tubular structure.
5. system according to claim 3, described measurement module (4) is arranged to determine from pointed position along the number of the branch that described tubular structure meets with.
6. system according to claim 1, also comprise recombination module (5), it is used to provide many curved surface reorganization, described many curved surfaces reorganization is followed the branch that reaches described tubular structure from pointed position and away from the curvature of at least a portion of the blood vessel of the branch of described tubular structure, and follows the curvature away from least a portion of two blood vessels of described branch.
7. system according to claim 1, described indicating module (2) is arranged to be used to respond the mouse moving event.
8. system according to claim 1, described 3-D view comprises volumetric image, described identification module (3) also comprises estimation module (6), and it is used for estimating local trace parameters in response to described user points out described position according to the described volumetric image of pointed position.
9. system according to claim 8, described estimation module (6) are arranged to be used for the weighted mean of the key element of computation structure tensor, calculate described weighted mean around pointed position above the zone of intended size.
10. a medical imaging workstation comprises system according to claim 1.
11. a method of carrying out vessels analysis comprises
The 3-D view of tubular structure is represented in demonstration (101) at least;
Make that (102) user can point out a position on blood vessel in tubular construction, to obtain pointed position; And
Identification (103) is positioned at the part of the described tubular structure around the pointed position, and described part comprises branch arbitrarily, and extends predetermined distance scale from pointed position, with the position that obtains to be discerned.
12. a computer program comprises being used to make processor system to carry out the instruction of method according to claim 11.
CN200980149367.7A 2008-12-10 2009-12-03 Vessel analysis Active CN102246206B (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP08171149.1 2008-12-10
EP08171149 2008-12-10
PCT/IB2009/055484 WO2010067276A1 (en) 2008-12-10 2009-12-03 Vessel analysis

Publications (2)

Publication Number Publication Date
CN102246206A true CN102246206A (en) 2011-11-16
CN102246206B CN102246206B (en) 2015-07-22

Family

ID=41718673

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200980149367.7A Active CN102246206B (en) 2008-12-10 2009-12-03 Vessel analysis

Country Status (7)

Country Link
US (1) US8611629B2 (en)
EP (1) EP2377095B1 (en)
JP (1) JP5715064B2 (en)
CN (1) CN102246206B (en)
BR (1) BRPI0917609A2 (en)
RU (1) RU2534948C2 (en)
WO (1) WO2010067276A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104185448A (en) * 2012-03-15 2014-12-03 富士胶片株式会社 Medical image display device, medical image display method and medical image display program
CN104427927A (en) * 2012-06-28 2015-03-18 皇家飞利浦有限公司 Fiber optic sensor guided navigation for vascular visualization and monitoring

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102010018460B4 (en) * 2010-04-27 2017-02-16 Siemens Healthcare Gmbh Method for determining at least one change of a tubular tissue structure of a living being, arithmetic unit and data carrier
US11412998B2 (en) 2011-02-10 2022-08-16 Karl Storz Imaging, Inc. Multi-source medical display
US10631712B2 (en) * 2011-02-10 2020-04-28 Karl Storz Imaging, Inc. Surgeon's aid for medical display
US10674968B2 (en) * 2011-02-10 2020-06-09 Karl Storz Imaging, Inc. Adjustable overlay patterns for medical display
EP2751726B1 (en) 2011-09-13 2019-10-09 Koninklijke Philips N.V. Vessel annotator
US8630764B2 (en) * 2011-12-02 2014-01-14 Continential Automotive Systems, Inc. System operation for telematics systems that use a battery
EP2863802B1 (en) * 2012-06-26 2020-11-04 Sync-RX, Ltd. Flow-related image processing in luminal organs
EP3503026A1 (en) * 2017-12-20 2019-06-26 Koninklijke Philips N.V. Device, system and method for interacting with vessel images
RU2684760C1 (en) * 2018-05-14 2019-04-12 Общество с ограниченной ответственностью "ЭНСИМ" Method and system for pre-operative modeling of medical procedure
US10813612B2 (en) 2019-01-25 2020-10-27 Cleerly, Inc. Systems and method of characterizing high risk plaques
KR102680085B1 (en) 2019-06-11 2024-06-28 서강대학교산학협력단 Vessel position tracking device and operating method thereof
US11501436B2 (en) 2020-01-07 2022-11-15 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US20220392065A1 (en) 2020-01-07 2022-12-08 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
US11969280B2 (en) 2020-01-07 2024-04-30 Cleerly, Inc. Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking
CN111681226B (en) * 2020-06-09 2024-07-12 上海联影医疗科技股份有限公司 Target tissue positioning method and device based on blood vessel identification
US20230289963A1 (en) 2022-03-10 2023-09-14 Cleerly, Inc. Systems, devices, and methods for non-invasive image-based plaque analysis and risk determination
CN115456950B (en) * 2022-08-15 2023-08-11 上海交通大学医学院附属第九人民医院 Automatic identification and positioning method, equipment and storage medium for puncture blood vessel

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010044576A1 (en) * 1994-10-27 2001-11-22 Vining David J. Method and system for producing interactive three-dimensional renderings of selected body organs having hollow lumens to enable simulated movement through the lumen
US20050110791A1 (en) * 2003-11-26 2005-05-26 Prabhu Krishnamoorthy Systems and methods for segmenting and displaying tubular vessels in volumetric imaging data
US20070031019A1 (en) * 2005-07-28 2007-02-08 David Lesage System and method for coronary artery segmentation of cardiac CT volumes
US20070160274A1 (en) * 2006-01-10 2007-07-12 Adi Mashiach System and method for segmenting structures in a series of images
CN101160602A (en) * 2005-04-12 2008-04-09 皇家飞利浦电子股份有限公司 A method, an apparatus and a computer program for segmenting an anatomic structure in a multi-dimensional dataset.
US20080095423A1 (en) * 2006-10-23 2008-04-24 Siemens Aktiengesellschaft Analysis method for data of an examination object, featuring correlation of a volume data record with a number of projection images

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04370882A (en) * 1991-06-20 1992-12-24 Hitachi Ltd Grid setting method
FR2708166A1 (en) * 1993-07-22 1995-01-27 Philips Laboratoire Electroniq A method of processing digitized images for the automatic detection of stenoses.
JP2776246B2 (en) * 1994-05-31 1998-07-16 日本電気株式会社 Mouse cursor tracking type enlarged display moving device
US6095976A (en) * 1997-06-19 2000-08-01 Medinol Ltd. Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen
JP2003511126A (en) * 1999-10-01 2003-03-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Analysis of object datasets in multidimensional space
AU2001229650A1 (en) * 2000-02-11 2001-08-20 The Government Of The United States Of America, As Represented By The Secretary Of The Department Of Health And Human Services Vessel delineation in magnetic resonance angiographic images
JP4109575B2 (en) * 2003-05-27 2008-07-02 アロカ株式会社 Ultrasonic image processing apparatus and three-dimensional data processing method
JP2008503303A (en) 2004-06-22 2008-02-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Method and apparatus for displaying a tracheobronchial tree
JP4503389B2 (en) * 2004-08-02 2010-07-14 株式会社日立メディコ Medical image display device
WO2006055031A2 (en) 2004-11-18 2006-05-26 Siemens Medical Solutions Usa, Inc. Method and system for local visualization for tubular structures
JP2006246941A (en) * 2005-03-08 2006-09-21 Toshiba Corp Image processing apparatus and vessel tracking method
US7379062B2 (en) * 2005-08-01 2008-05-27 Barco Nv Method for determining a path along a biological object with a lumen
US7711164B2 (en) * 2005-08-02 2010-05-04 Siemens Medical Solutiions Usa, Inc. System and method for automatic segmentation of vessels in breast MR sequences
WO2007129616A1 (en) * 2006-05-02 2007-11-15 National University Corporation Nagoya University Insertion assist system of endoscope and insertion assist method of endoscope
JP4783219B2 (en) * 2006-06-16 2011-09-28 株式会社トプコン Ophthalmic imaging equipment
US7860283B2 (en) 2006-10-25 2010-12-28 Rcadia Medical Imaging Ltd. Method and system for the presentation of blood vessel structures and identified pathologies
JP4545169B2 (en) * 2007-04-12 2010-09-15 富士フイルム株式会社 Image display method, apparatus and program

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010044576A1 (en) * 1994-10-27 2001-11-22 Vining David J. Method and system for producing interactive three-dimensional renderings of selected body organs having hollow lumens to enable simulated movement through the lumen
US20050110791A1 (en) * 2003-11-26 2005-05-26 Prabhu Krishnamoorthy Systems and methods for segmenting and displaying tubular vessels in volumetric imaging data
CN101160602A (en) * 2005-04-12 2008-04-09 皇家飞利浦电子股份有限公司 A method, an apparatus and a computer program for segmenting an anatomic structure in a multi-dimensional dataset.
US20070031019A1 (en) * 2005-07-28 2007-02-08 David Lesage System and method for coronary artery segmentation of cardiac CT volumes
US20070160274A1 (en) * 2006-01-10 2007-07-12 Adi Mashiach System and method for segmenting structures in a series of images
US20080095423A1 (en) * 2006-10-23 2008-04-24 Siemens Aktiengesellschaft Analysis method for data of an examination object, featuring correlation of a volume data record with a number of projection images

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104185448A (en) * 2012-03-15 2014-12-03 富士胶片株式会社 Medical image display device, medical image display method and medical image display program
US9536346B2 (en) 2012-03-15 2017-01-03 Fujifilm Corporation Medical image display apparatus, medical image display method, and medical image display program
CN104427927A (en) * 2012-06-28 2015-03-18 皇家飞利浦有限公司 Fiber optic sensor guided navigation for vascular visualization and monitoring

Also Published As

Publication number Publication date
CN102246206B (en) 2015-07-22
RU2534948C2 (en) 2014-12-10
JP5715064B2 (en) 2015-05-07
JP2012511380A (en) 2012-05-24
US8611629B2 (en) 2013-12-17
EP2377095B1 (en) 2016-05-25
EP2377095A1 (en) 2011-10-19
WO2010067276A1 (en) 2010-06-17
BRPI0917609A2 (en) 2019-10-15
US20110235891A1 (en) 2011-09-29
RU2011128347A (en) 2013-01-20

Similar Documents

Publication Publication Date Title
CN102246206B (en) Vessel analysis
CN101548296B (en) Automated hierarchical splitting of anatomical trees
CN102171724B (en) The selection of medical image sequences snapshot
Mairhofer et al. RooTrak: automated recovery of three-dimensional plant root architecture in soil from X-ray microcomputed tomography images using visual tracking
US7773791B2 (en) Analyzing lesions in a medical digital image
CN106157295B (en) For dividing the calculating and visual system and method for uncertainty in medical image
JP5802175B2 (en) Image processing device
US20140105474A1 (en) Apparatus and method for analyzing medical image
EP2846310A2 (en) Method and apparatus for registering medical images
CN103987337A (en) Distorsion fingerprinting for EM tracking compensation, detection and error correction
CN101953693A (en) The estimation based on model of complete or part pet decay pattern
CN105792738A (en) Local FFR estimation and visualization for improved functional stenosis analysis
CN103608843A (en) Image display apparatus
CN106062753A (en) Enabling a user to study image data
JP2017526440A (en) Medical imaging device
CN106073824A (en) For carrying out the method and system of the assessment automatically of modulus of shearing and viscosity according to shearing wave imaging
CN104463825B (en) Device and method for the test object in three-dimensional volumetric images
Vogel et al. Trajectory optimization for intra-operative nuclear tomographic imaging
CN103415873A (en) TNM classification using image overlays
RU2508056C2 (en) Method of composition and calculation of volume in system of ultrasound visualisation
CN102656586A (en) Associating acquired images with objects
JPH0889501A (en) Medical image diagnostic device
Gibbons et al. MRI-derived cardiac fat modelling for use in ultrasound tissue labelling and classification
US12070357B2 (en) System and method for automatic association and display of video loop subject matter for enhanced identification
Guo et al. 3D reconstruction of ultrasonic carotid artery based on fractal dimension and marching cubes

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant