CN101563707A - Method and system for fiber tracking - Google Patents
Method and system for fiber tracking Download PDFInfo
- Publication number
- CN101563707A CN101563707A CNA2006800450882A CN200680045088A CN101563707A CN 101563707 A CN101563707 A CN 101563707A CN A2006800450882 A CNA2006800450882 A CN A2006800450882A CN 200680045088 A CN200680045088 A CN 200680045088A CN 101563707 A CN101563707 A CN 101563707A
- Authority
- CN
- China
- Prior art keywords
- tensor
- value
- vector
- threshold
- equal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Abstract
One embodiment of the present invention includes a system comprising to determine a direction of tracking a fiber based on a vector corresponding to a largest value of a set of values for a tensor, if an anisotropy value of the tensor is greater than or equal to a first threshold; and to apply a weighted function to the vector of the tensor to select the direction of tracking the fiber, if the anisotropy value of the tensor is less than or equal to the first threshold and larger than or equal to a second threshold. In one embodiment, the system is to generate the weighted function comprises a linear interpolation to be performed on the vector. In one embodiment, a weight of the linear interpolation is in part dependent on the anisotropy value of the tensor, wherein the weight corresponds to the anisotropy value relative to the first threshold. One embodiment of the invention is also executable as a method.
Description
Related application
The application's requirement was submitted on November 30th, 2005, branch is equipped with the U.S. Provisional Application No. of the title of application number 60/741356 for " Method and System for Fiber Tracking ", and this provisional application all is included in this.
Background technology
It is just in the Research on development field that diffusion tensor imaging (DTI) develops.Scanner is collected data preferably all the time, and doctor and scientist find the new application of these data by persistence.The success of diffusion magnetic resonance imaging (MRI) comes from effective notion, this effective notion just molecule at random, during the displacement of diffusion-driven, molecular detection exceeds the institutional framework of the micro-scale of common image resolution ratio.Because diffusion is a three dimensional process, thus can be anisotropic in in-house molecular migration rate, as in the white matter of brain.
Provide the details of tissue microstructure, can extract, distinguish and obtain (exploit) diffusion anisotropic effect.It is the application of the fiber tracking in brain that such improvement is used, and this application can provide seeing clearly for the connectivity problem.Also use DTI that subtle anomalies in the multiple disease (comprising apoplexy, multiple sclerosis, dislexia and schizophrenia) is shown, and in the current part that becomes a lot of routine clinical schemes.
Yet, exist demand for more intuitive inputting interface so that allow the user specify interested fibrous bundle (tract), thereby make them become the part of surgery planning and subsequent navigation.
Summary of the invention
One embodiment of the present of invention comprise system, are used for: if the anisotropy value of tensor is more than or equal to first threshold, then based on the direction of determining to follow the tracks of fiber corresponding to the peaked vector in the class value of tensor; And, if the anisotropy value of tensor is less than or equal to first threshold and more than or equal to second threshold value, then the vector of tensor is used the direction that weighting function selects to follow the tracks of fiber.In one embodiment, this system is used to generate weighting function, and this weighting function comprises the linear interpolation that will carry out on vector.In one embodiment, the weight of linear interpolation partly depends on the anisotropy value of tensor, and wherein, this weight is corresponding to the anisotropy value relevant with second threshold value with first threshold.One embodiment of the present of invention can also be carried out as method.
Description of drawings
Patent or application documents comprise at least one color drawings.The disclosed copy with color drawings of this patent or patented claim will provide after the required expense by Patent Office's request of receiving and payment.
Fig. 1 a-b illustrates the inputting interface according to an embodiment.
Fig. 2 further illustrates an embodiment of compute tensor module.
Fig. 3 illustrates an embodiment of the visualization module of the development that is used for selecting the tensor that calculates.
Fig. 4 a-h illustrates the example via the reference volume of difference development according to an embodiment.
Fig. 5 illustrates according to an embodiment and identifies three dimensional area of interest.
Fig. 6 illustrates according to an embodiment and identifies a plurality of three dimensional area of interest.
Fig. 7 illustrates according to an embodiment and adds group of fibers to another group of fibers.
Fig. 8 illustrates according to an embodiment and deletes group of fibers.
Fig. 9 illustrates the color according to an embodiment rename group of fibers and/or change group of fibers.
Figure 10 illustrates according to an embodiment and identifies two-dimensional regions of interest.
Figure 11 presents the pretreated process flow diagram that be used to generate fibrous bundle of description according to an embodiment.
Figure 12 presents description according to the tracking of the embodiment process flow diagram with the process that generates fibrous bundle.
Figure 13 presents the process flow diagram of description according to the process of the generation tracking direction of an embodiment.
Embodiment
One embodiment of the present of invention comprise the DTI module interfaces of being made up of one or more submodules, " compute tensor ", " development " just of these submodules, " fiber tracking " and " fiber management ".In other embodiment, can use the difference collection of submodule.In one embodiment, can utilize the inputting interface of 3 spatial degrees of freedom that are provided for controlling input to operate the DTI module.The hand that inputting interface is provided in the reference volume effectively touches closely (handaccess).
The use compute tensor module is selected source DTI volume and is used for the parameter of compute tensor.Visualization module provides the set of visualizations of volume, and these set of visualizations can be assisted diagnosis.In other embodiment, can load source DTI volume, can come compute tensor automatically based on the preset parameter collection, and show via predetermined visualization.
Use fiber track module to follow the tracks of and visualize neuron fibers.In one embodiment, the area-of-interest in reference volume is designated as one or more three dimensional area of interest (for example cube).
Also provide fiber management module to be organized in the fiber that generates in the fiber track module in one embodiment.For example, fiber management module allows that the user adds, rename or deletion fiber.In one embodiment, also support paint tool to restain fiber and also visually distinguish multiple fiber.
Inputting interface
As mentioned above, in one embodiment, can utilize provides the inputting interface (being also referred to as " input interface " here) of at least 3 spatial degrees of freedom that are used for controlling input to operate the DTI module.Inputting interface provides three dimensions to handle alternately.Fig. 1 illustrates the inputting interface according to an embodiment.The hand that inputting interface is provided in the reference volume effectively touches closely.
As description, for (rendered) object or reference volume to 3 D rendering work, inputting interface comprises one or more hand-held apparatus 102a-b, such as contact pilotage.In one embodiment, this apparatus allows by the user freely with 3 or more freedom control instrument.The correspondence image of hand-held apparatus (being also referred to as " virtual tool " here) is shown with mutual with reference volume together with its motion.
In one embodiment, the manipulation of hand-held device produces the correspondence of shown reference volume and handles, and this reference volume can be shown as 3-D view.The result is to rotate and to move in different directions this reference volume from different perspectives.In one embodiment, the user activates this hand-held device (for example by the button that is pressed on the hand-held device 102b), so that the motion of hand-held device produces the motion of shown reference volume, comprises position and the orientation of handling this reference volume.
In one embodiment, utilize radio frequency (rf) tracker to follow the tracks of the motion of hand-held apparatus.Can select among the embodiment, inputting interface can comprise that being used to utilize one or more hand-held devices to provide imports the haptic device of controlling.
Can select to provide hand-held device among the embodiment more than two.For example, a plurality of users can be mutual with reference volume long-range or locally.In one embodiment, can use single hand-held device to carry out one group of above-mentioned behavior.
In one embodiment, mirror 103 is arranged between user and the computer screen 106.Reference volume and the virtual tool of mirror reflection as showing by computer screen.User's hand can the work space 108 behind mirror in motion and with undertaken alternately by the reference volume shown in the reflection.The result is that the user can utilize both hands can not blur reference volume to the reference volume effect.In addition, in one embodiment, inputting interface is provided with workstation1 12, and this workstation1 12 comprises that support 114 shelves for user's arm.
In one embodiment, can on screen, show reference volume three-dimensionally.In one embodiment, LCD (LCD) jalousie glass 120 is used to three-dimensional ground perception reference volume.The LCD jalousie glass allow by the notion of using the alternate frame ordering make light by with graphoscope on image synchronization.A plurality of observers can wear jalousie glass and side by side check and discuss reference volume.
Usually the button and the slide block that need frequent access control application and actuating operator scheme and instrument such as the complicated applications of neurosurgery planning.Thus, embodiment comprises that virtual tool panel provides the integrated of work space with at least 3 spatial degrees of freedom and application controls.
In one embodiment, shown in Fig. 1 b, the solid surface 116 below virtual tool panel and the work space 108 overlaps.In one embodiment, virtual tool panel can be in response to touching pedestal 116 and present with handheld tool 102a, and disappear after removing instrument away from pedestal.Can comprise other behaviors alternately between handheld tool and virtual tool panel, these behaviors comprise the one or more of pressing button, towing slide block, control curve etc.
In one embodiment, can utilize process or module to operate inputting interface, this process or module can generate multi-modality images and from the real-time volume drawing of cutting apart and the three dimensions surface rendering of one or more acquisitions of multi-modality images based on following technology one or more: computed tomography (CT), positron emission tomography (PET), single photon emission tomography (SPECT), magnetic resonance imaging (MRI), magnetic resonance angiography imaging (MRA), volume are ultrasonic.By using inputting interface together with this process, tridimensional virtual reality (VR) environment is provided in one embodiment, wherein, the user can " reach in the three-dimensional data " to come alternatively to work with three-dimensional data in real time by utilizing both hands.
In one embodiment, this process comprises one or more following features: have an X-rayed three-dimensional screen body and draw and surface rendering; Multi-modality images merges; Automaton registration and the checking that is registered object; Cut apart; The operation prospecting tools that is used for cutting, cutting, boring, recovery, clone, roaming, linearity and volumetric measurement; Default colour and/or transparent mapped with volume drawing; According to DICOM, and support a plurality of file layouts (for example TIFF); And utilize three-dimensional playback and video derivation performance to catch three dimensions and handle alternately.
Diffusion tensor imaging (DTI) module
As mentioned above, in one embodiment, the DTI module comprises one group of submodule.These submodules are compute tensor module, visualization module, fiber track module and fiber management module.But can use the choice set of module and not depart from the present invention.For example, in one embodiment, can load source DTI volume, and can come compute tensor automatically together with the predetermined visualization way of selecting in advance based on the preset parameter collection.In this embodiment, can not comprise that in the DTI module calculating and/or visualization module are for user interactions.
Compute tensor module
Fig. 2 further illustrates an embodiment of compute tensor module.Compute tensor module is used for calculating diffusion tensor according to the source DTI volume that comprises a group 6 or a plurality of diffusion weighted images (DWI).The panel 202 of compute tensor module allows the user to select the intensity threshold of source DTI volume, diffusion sensitizing parameter (being also referred to as the b value) and DWI.
Provide the Object Selection device 204 of panel 202 to select source DTI volume.In one embodiment, provide the quantity selector switch 206 in the panel to specify in the b value of determining during the DWI acquisition process.Use the b value parameter to come compute tensor.
When pressing compute tensor button 210, calculate the tensor of source DTI volume.In one embodiment, feature 212 letters that provide on panel 202 are shown and are used for the parameter of compute tensor.In one embodiment, feature 212 only occurs after compute tensor.Can select among the embodiment, the parts that provide on the panel 202 of compute tensor module can change and not depart from the present invention.
Visualization module
Fig. 3 illustrates an embodiment of visualization module, and this visualization module is used for selecting the development of institute's compute tensor.In one embodiment, one group of development is arranged on the panel 302 of visualization module.In one embodiment, develop comprise following one or more: FA (fractional anisotropy) volume 304, ADC (apparent diffusion coefficient) volume 306, FA color volume 308, SEC (shape coding colour) volume 310, LA (linear anisotropic) volume 312, PA (in-plane anisotropy) volume 314, SA (spherical anisotropy) volume 316 and direction (maximum dispersal direction) line 318.Can use selectable set of visualizations and not depart from the present invention.For example, visualization module embodiment can not comprise FA.
In one embodiment, by being chosen in the development on the panel and then pressing " computed volume " button 320 on panel, draw the DTI tensor of source DTI volume.
FA, ADC, LA, PA, SA volume are the grayscale volumes that the diffusion property of DTI tensor is shown.In the FA volume, has the diffusion anisotropy more of the voxel indication tensor in this voxel of higher-strength, as by shown in the example among Fig. 4 (a).The ADC volume is represented the average diffusion rate of brain.The average diffusion intensity of the DTI tensor of intensity indication in this voxel of the voxel of ADC volume is shown in example in passing through Fig. 4 (b).LA, PA and SA volume are represented linearity, plane, the spherical diffusion properties of DTI tensor respectively, shown in Fig. 4 (d-f).Higher linearity, plane and the spherical diffusion of the tensor of higher-strength indication in this voxel in the voxel in LA, PA and the SA volume.
In one embodiment, the FA color volume is the color-coded volume that the directional information of tensor is shown, shown in example in passing through Fig. 4 (f).In one embodiment, the SEC volume provides the color-coded volume coding of the shape information of DTI tensor, shown in example in passing through Fig. 4 (g).The shape of the diffusion tensor of color indication in this voxel of voxel.In the example of Fig. 4 (g), red indication diffusion tensor is an oblong, and yellow and white to indicate diffusion tensor respectively be oblate and spherical.Can select to use selectable color combination among the embodiment.
In one embodiment, the direction method generates and shows one group of straight line of indicating the maximum dispersal direction of all DTI tensors, shown in example in passing through Fig. 4 (h).
Fiber track module
In one embodiment, can the tensor that calculate be developed via fiber track module is fibrous bundle (being also referred to as fiber bundle here).An embodiment allows the user to identify two dimension and/or three dimensional area of interest on reference volume, and calculates the fiber by these area-of-interests.In one embodiment, as mentioned above, can utilize the inputting interface of at least 3 spatial degrees of freedom that are provided for controlling input, operate fiber track module.In one embodiment, as shown in Figure 5, use the virtual tool 502 of inputting interface in shown reference volume 506, to identify three dimensional area of interest 504.
In one embodiment, as shown in Fig. 5, in order to select three dimensional area of interest in reference volume, the three-dimensional cube 508 of user's manipulation at the place, tip of virtual tool 502 is identified at the area-of-interest 504 in the reference volume 506 as further.In one embodiment, after handling three-dimensional cube 508 above the zone of user in reference volume, the user can select the zone of (for example pressing the button at the sidepiece place of virtual tool) reference volume as three dimensional area of interest 504.In response, three-dimensional cube remains on the reference volume 506 of mark three dimensional area of interest 504, as shown in Figure 5.After this, the user can calculate/generate the fibrous bundle that passes through three dimensional area of interest 504 via the control panel request via computed push-buttom 510.
In one embodiment, calculating and demonstration are by the fibrous bundle of area-of-interest.Selectively, can generate at the fiber of whole volume but do not shown, and identify thus and show fibrous bundle by area-of-interest.Therefore, as described herein, in one embodiment, the calculating that relates to and/or generate fibrous bundle or fiber bundle can comprise: sign and display fibers bundle, calculate and the display fibers bundle, or calculate and display fibers bundle not.
In one embodiment, as shown in Figure 6, can on reference volume, select a plurality of three dimensional area of interest 604a-c.The user can select to generate the fiber (for example common factor pattern 622 on the control panel 620) that only passes through the three dimensional area of interest of all selections, or selects to find all fibres (for example integrated mode 624 on the panel 620) by at least one of three dimensional area of interest.
Further as shown in Figure 6, after the user specified area-of-interest and selects to calculate fiber button 610, newly-generated fibrous bundle was automatically added to effective group of fibers.In addition, the user can come the newly-generated fibrous bundle of deletion from current effective group of fibers by selecting recovery button.The user can be by repeatedly selecting recovery button, delete continuously generate at different time and be included in fibrous bundle in effective group of fibers.
In one embodiment, can regulate the size of three-dimensional cube in real time, thereby and change the size of the area-of-interest that will select.Can be via the input controlling features that is included on the control panel, or via the input controlling features on the hand-held device that is included in inputting interface, the size of regulating three-dimensional cube.
Can select among the embodiment, can use the object except three-dimensional cube and do not depart from the present invention.For example, can provide default three dimensional area of interest.Can define and arrange default three dimensional area of interest by probabilistic method.For example, can use the DTI atlas information of common registration to come automatic detection may comprise the zone of special fiber bundle.Can use inputting interface and the DTI module described to come modifier area then here from the size or the vpg connection of default three dimensional area of interest.
In one embodiment, as shown in Figure 7, add newly-generated group of fibers 732 to second group of fibers 730 via the interpolation button 742 on the control panel of fiber management module 700.In one embodiment, group of fibers is the container that is used for fibrous bundle.In one embodiment, at any time, have a group of fibers effective, this effective group of fibers is called as current group of fibers.As shown in Figure 8, can be from concentrated deletion first group of fibers of group of fibers.In one embodiment, it is empty up to the group of fibers change that the user can delete the fibrous bundle that generates at different time continuously.
As shown in Figure 9, can also in fiber management module, change into another Show Color to the group of fibers rename and with it.As default, come the fiber in the group of fibers is carried out coloud coding based on diffusion direction information.Also can carry out again painted to them to distinguish with other group of fibers.Example as shown in Figure 9 has two selections of fiber color, and one is monochromatic, and another is the direction color.
In one embodiment, the user can also delete the three dimensional area of interest of mark before of utilizing the three-dimensional cube sign that is arranged in the reference volume.In one embodiment, the user can select the delete button 626 on the control panel 620 of fiber track module, as shown in Figure 6.The user can handle virtual tool then near the three-dimensional cube that adds before.When virtual tool is near the three-dimensional cube that adds before predetermined, highlight is given prominence to this three-dimensional cube.Activate (for example pressing the button on virtual tool one side) in response to the user, with the outstanding three-dimensional cube of deletion highlight.
In one embodiment, can identify zone in the reference volume to be avoided and can not be the part that is used to generate the area-of-interest of fibrous bundle.In one embodiment, can in the three dimensional area of interest interface, (for example on the panel of fiber track module) provide additional buttons.Can arrange on reference volume that one or more cubes identify the area-of-interest that will avoid.In one embodiment, adopt different colors (or different shapes) corresponding to the cube of the area-of-interest that will avoid with respect to the cube that sign is used to generate the area-of-interest of fibrous bundle.When the user begins to calculate fiber, the fiber tracking algorithm will be considered the area-of-interest that these are avoided, and abandon those fibers of avoiding area-of-interest by these.
In one embodiment, as shown in figure 10, select two-dimensional regions of interest selector button 1060, will show viewer 1062 in response to the user.Section is handled at two dimension slicing and some interfaces of viewer demonstration reference volume.Utilize viewer, the user can select to check reference volume axially, radially and crown section, and amplify or dwindle two dimension slicing.The user can draw outline specifies the two-dimensional regions of interest in the section, or selectively specifies a plurality of two-dimensional regions of interest in the section.
In one embodiment, the reference volume of using with fiber track module can be the DTI volume, or CT, PET, SPECT, MRI, MRA, volume is ultrasonic or with other multi-modal volumes of the common registration of DTI volume.In addition, in one embodiment, can use from the split image of one or more acquisitions of the multi-modal volume of the common registration of DTI volume, as the reference volume.
In one embodiment, fiber tracking panel 620 comprises that also the interface regulates the stop condition that is used for fiber tracking.In one embodiment, stop condition comprises following one or more: FA threshold value, maximum length threshold, minimum length threshold and deviation angle threshold.By regulating above-mentioned threshold value, the user can obtain to have the fiber of difformity and length.Use step length slider to come the smoothness of controlling fiber.For example, if step-length is less, then fiber can be shown as more level and smooth and more accurate.Use these control, the user can find the parameter that obtains expectation fiber tracking result then.
Fiber tracking
Figure 11 and 12 illustrates description according to the generation of an embodiment and the process flow diagram of following the tracks of the process of fiber.Figure 11 illustrates and describes the pretreated process flow diagram that is used to generate fibrous bundle.Handling in 1102, load diffusion weighted images (DWI) (being also referred to as source DTI volume), and be the voxel compute tensor of DWI.Handling in 1104, using and handle to come level and smooth tensor by the noise that reduces data.In one embodiment, use gaussian kernel and come level and smooth tensor.Can select among the embodiment, can use other processing and do not depart from scope of the present invention.Handling in 1106, for tensor computation proper vector and eigenwert.
Figure 12 illustrates description according to the tracking of the embodiment process flow diagram with the process that generates fibrous bundle.Handling in 1202, as previously mentioned, the one or more area-of-interests of sign in reference volume.For example, can select a plurality of area-of-interests, and/or can select area-of-interest as three dimensional area of interest.In addition, in one embodiment, can select the area-of-interest collection, so that the fiber that the result obtains is by the area-of-interest (pattern of just occuring simultaneously) of all selections.
Handling in 1204, be used for the tensor of the voxel of area-of-interest corresponding to the proper vector of eigenvalue of maximum by sign.The tracking of respective fiber is carried out along the direction of this proper vector.
Handling in 1206, following the tracks of along current tracking direction and continue than short distance and reach new point.In one embodiment, to be called as step-length can be fixed value or fitness value in this short distance.If fixed step size then is fixed as step-length certain value.In one embodiment, the user can regulate this value by described step length slider before regulating.If step-length is adaptive, then distance changes at the anisotropy value according to the tensor on before during the fiber tracking process.
Handling in 1208, if new putting outside the border that is in respect to source DTI volume follows the tracks of stopping 1210.Especially, in one embodiment, reference volume has the bounding box of the size of indication source DTI volume.If new point is positioned at this bounding box outside, then follow the tracks of outside the border and follow the tracks of stopping.Otherwise, to handle in 1212, the use tri-linear interpolation process is calculated the new tensor at new some place, and is this new tensor calculated characteristics vector sum eigenwert.
Handling in 1214, if the fractional anisotropy (FA of this new tensor
New) be worth less than FA
1(predetermined threshold just) then followed the tracks of and stopped 1210.Handling in 1216, if FA
NewAt FA
1With FA
2Between (FA
2Be predetermined threshold and FA
2＞FA
1), then use independent processing to generate next alternative tracking direction.Otherwise, to handle in 1218, the tracking of each fiber is carried out along the direction corresponding to the proper vector of the eigenvalue of maximum of new tensor.Can select among the embodiment, can use the measurement except FA and do not depart from the present invention.
Handling in 1220, if the deviation angle between current tracking direction and next alternative tracking direction greater than predetermined threshold, is then followed the tracks of and stopped 1210.Otherwise, to handle in 1222, fiber tracking continues along next tracking direction, and handles at processing 1206 places and continue once more.In one embodiment, carry out processing 1202 to 1222 for one group of voxel of the region of interest domain identifier that utilizes sign.Can select among the embodiment, can remove some described processing, and can comprise additional processing, and not depart from scope of the present invention.
Figure 13 illustrates to describe and works as FA as reference process 1216
NewAt FA
1With FA
2Between the time generate the embodiment of process flow diagram of the processing of tracking direction.Handling in 1302, determine proper vector e1 corresponding to the eigenvalue of maximum of new tensor.Handling in 1304, use the proper vector (e1) and the current tracking direction (v that determine in 1302 handling
In) interpolation generate next tracking direction.In one embodiment, interpolation comprises following linear interpolation: vector (v)=(1-w) * v
In+ W*e
1, wherein, w is the weight of interpolation.Handling in 1306, the designated result vector is (v) as next tracking direction.Can select among the embodiment, can use the variation of interpolation to determine tracking direction and do not depart from scope of the present invention.
In one embodiment, the weight (w) that is used for linear interpolation is FA
NewFunction, FA wherein
New=f (w).In one embodiment, work as FA
NewHour, weight (w) is also less, and vice versa.
In addition, in one embodiment, regulate the variation of weight with respect to the variation of FA value.In one embodiment, in order to regulate the variation of weight, use second order parabola (f), this second order parabola (f) is defined as (f)=a*w*w+b*w+c.In one embodiment, coefficient a, b, c are based on that the predetermined correspondence of using FA value and weight determines.Can select among the embodiment, can use the variation of the weight that is identified for interpolation and do not depart from scope of the present invention.
Sum up
Very clear, can make the many modifications and variations of this embodiment by those skilled in the art and do not depart from the spirit of innovative techniques of the present disclosure.
Above-mentioned processing can be used as one group of pending instruction and is stored in the storer of computer system.In addition, the instruction of carrying out above-mentioned processing can also be stored on the other forms of machine readable medium, comprises disk and CD.For example, described processing can be stored on the machine readable medium, and such as disk or CD, these disks or CD can drive (or computer-readable medium driving) via dish and visit.In addition, can via data network instruction be downloaded in the calculation element with the form of compiling and linked version.
Selectively, the logic that is used for carrying out above-mentioned processing can realize in additional computer and/or machine readable medium, these machine readable mediums such as discrete hardware components such as large scale integrated circuit (LSI), specific use integrated circuit (ASIC), firmware such as electric Erasable Programmable Read Only Memory EPROM (EEPROM) and electric, optics, sound and other forms of transmitting signal (for example carrier wave, infrared signal, digital signal etc.) etc.
Very clearly can make the many modifications and variations of this embodiment by those skilled in the art and do not depart from the spirit of innovative techniques of the present disclosure.
Claims (20)
1. method comprises:
If the anisotropy value of first tensor is more than or equal to first threshold, then based on the direction of selecting to follow the tracks of fiber corresponding to peaked first vector of a class value of described first tensor; And
If the anisotropy value of described first tensor is less than or equal to described first threshold and more than or equal to second threshold value, then the vector of described first tensor is used the direction that weighting function selects to follow the tracks of fiber.
2. method according to claim 1, wherein, described vector is a proper vector, and the eigenwert that a described maximal value and a described class value are described first tensors.
3. method according to claim 1 wherein, is used described weighting function to described vector and is also comprised described vector application linear interpolation.
4. method according to claim 3, wherein, described linear interpolation is at least in part based on vector before.
5. method according to claim 4, wherein, the weight of described linear interpolation is at least in part based on the anisotropy value of described first tensor.
6. method according to claim 5 also comprises: if the anisotropy value of described first tensor is less than or equal to second threshold value, then stop the tracking of described fiber, described second threshold value is less than described first threshold.
7. method according to claim 5 also comprises: before vector that calculates described first tensor and value, described first tensor is used filtrator.
8. method according to claim 7, wherein, described filtrator comprises Gaussian function.
9. method according to claim 6 also comprises: after following the tracks of preset distance, carry out Tri linear interpolation and generate second tensor; And
If the anisotropy value of described second tensor more than or equal to first threshold, then based on corresponding to the peaked vector in the class value of described second tensor, determines to follow the tracks of the next direction of fiber; And
If the anisotropy value of described second tensor is less than or equal to described first threshold and more than or equal to described second threshold value, then the vector of described second tensor is used the direction that weighting function selects to follow the tracks of fiber.
10. method according to claim 5, wherein, the weight of described linear interpolation is partly calculated via the second order parabola function.
11. method according to claim 5, wherein, described anisotropy value is fractional anisotropy (FA) value.
12. method according to claim 5, wherein, weight is corresponding to the anisotropy value relevant with second threshold value with described first threshold.
13. a method comprises:
If the anisotropy value of tensor is more than or equal to first threshold, then based on corresponding to the peaked vector in the class value of the tensor of voxel, the direction of determining to follow the tracks of fiber; And
If the anisotropy value of described tensor less than described first threshold, then at least in part based on described vector and the combination of tracking direction before, is determined the direction of the described fiber of tracking.
14. store the machine readable medium of one group of instruction above one kind, these instructions are carried out when operation and are comprised following method:
If the anisotropy value of first tensor more than or equal to first threshold, then based on corresponding to the peaked vector in the class value of described first tensor, is determined the direction of tracking fiber; And
If the anisotropy value of described first tensor is less than or equal to described first threshold and more than or equal to second threshold value, then the vector of described first tensor is used the direction that weighting function selects to follow the tracks of described fiber.
15. machine readable medium according to claim 14 wherein, is used described weighting function to described vector and also comprised: the vector to described first tensor is used linear interpolation.
16. machine readable medium according to claim 14, wherein, described method also comprises:
After following the tracks of preset distance, carry out Tri linear interpolation and generate second tensor; And
If the anisotropy value of this second tensor more than or equal to described first threshold, then based on corresponding to the peaked vector in the class value of described second tensor, determines to follow the tracks of the next direction of described fiber; And
If the anisotropy value of described second tensor is less than or equal to described first threshold and more than or equal to described second threshold value, then the vector of described second tensor is used the direction that described weighting function selects to follow the tracks of described fiber.
17. machine readable medium according to claim 15, wherein, the weight of described linear interpolation partly depends on the anisotropy value of described first tensor.
18. a system comprises:
Device is if the anisotropy value that is used for tensor is more than or equal to first threshold, then based on the direction of determining to follow the tracks of fiber corresponding to the peaked vector in the class value of this tensor; And
Device, the anisotropy value of described tensor is less than or equal to described first threshold and more than or equal to second threshold value, then the vector of described tensor is used the direction that weighting function selects to follow the tracks of described fiber if be used for.
19. a system comprises:
The unit is if the anisotropy value that is used for tensor more than or equal to first threshold, then based on corresponding to the peaked vector in the class value of this tensor, is determined the direction of tracking fiber; And
The unit, the anisotropy value of described tensor is less than or equal to described first threshold and more than or equal to second threshold value, then the vector of described tensor is used the direction that weighting function selects to follow the tracks of described fiber if be used for.
20. system according to claim 19 wherein, is used for the unit of described vector application weighting function is used for the vector of described tensor is used linear interpolation.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US74135605P | 2005-11-30 | 2005-11-30 | |
| US60/741,356 | 2005-11-30 | ||
| US11/336,376 | 2006-01-20 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN101563707A true CN101563707A (en) | 2009-10-21 |
Family
ID=41221586
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CNA2006800450882A Pending CN101563707A (en) | 2005-11-30 | 2006-11-30 | Method and system for fiber tracking |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101563707A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110111429A (en) * | 2019-03-16 | 2019-08-09 | 哈尔滨理工大学 | A method of detection single pixel blood vessel |
| CN112489220A (en) * | 2020-10-23 | 2021-03-12 | 浙江工业大学 | Nerve fiber continuous tracking method based on flow field distribution |
-
2006
- 2006-11-30 CN CNA2006800450882A patent/CN101563707A/en active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110111429A (en) * | 2019-03-16 | 2019-08-09 | 哈尔滨理工大学 | A method of detection single pixel blood vessel |
| CN110111429B (en) * | 2019-03-16 | 2022-11-18 | 哈尔滨理工大学 | Method for detecting single-pixel blood vessel |
| CN112489220A (en) * | 2020-10-23 | 2021-03-12 | 浙江工业大学 | Nerve fiber continuous tracking method based on flow field distribution |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101535828A (en) | Method and system for diffusion tensor imaging | |
| US7411393B2 (en) | Method and system for fiber tracking | |
| US20070279435A1 (en) | Method and system for selective visualization and interaction with 3D image data | |
| CN1981717B (en) | Method and apparatus for integrating three-dimensional and two dimensional monitors with medical diagnostic imaging workstations | |
| US9014438B2 (en) | Method and apparatus featuring simple click style interactions according to a clinical task workflow | |
| US20070279436A1 (en) | Method and system for selective visualization and interaction with 3D image data, in a tunnel viewer | |
| Coffey et al. | Interactive slice WIM: Navigating and interrogating volume data sets using a multisurface, multitouch VR interface | |
| CN102246206B (en) | Vessel analysis | |
| CN101604458A (en) | The method that is used for the computer aided diagnosis results of display of pre-rendered | |
| US20090322748A1 (en) | Methods,systems, and computer program products for GPU-based point radiation for interactive volume sculpting and segmentation | |
| JP6434959B2 (en) | Enabling users to study image data | |
| Bornik et al. | A hybrid user interface for manipulation of volumetric medical data | |
| US7889894B2 (en) | Method of navigation in three-dimensional image data | |
| Bornik et al. | Computer-aided liver surgery planning: an augmented reality approach | |
| JP2009525058A (en) | Method and system for diffusion tensor imaging | |
| CN101118574A (en) | Systems and methods for rule-based volume rendition and navigation | |
| CN101563707A (en) | Method and system for fiber tracking | |
| CN103339660A (en) | A method, an apparatus and an arrangement for visualizing information | |
| JP2006516909A (en) | Medical image display method and apparatus | |
| LaViola et al. | Virtual reality-based interactive scientific visualization environments | |
| Beard et al. | Interacting with image hierarchies for fast and accurate object segmentation | |
| Abhari et al. | Use of a mixed-reality system to improve the planning of brain tumour resections: preliminary results | |
| Durisetti et al. | Evaluation of Tools Used for 3D Reconstruction of 2D Medical Images | |
| Camp et al. | A system for interactive volume analysis (SIVA) of 4-D biomedical images | |
| Fairfield et al. | Volume curtaining: a focus+ context effect for multimodal volume visualization |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20091021 |