WO2002064034A1 - Method and apparatus for the non-invasive imaging of anatomic tissue structures - Google Patents

Method and apparatus for the non-invasive imaging of anatomic tissue structures Download PDF

Info

Publication number
WO2002064034A1
WO2002064034A1 PCT/US2001/004479 US0104479W WO02064034A1 WO 2002064034 A1 WO2002064034 A1 WO 2002064034A1 US 0104479 W US0104479 W US 0104479W WO 02064034 A1 WO02064034 A1 WO 02064034A1
Authority
WO
WIPO (PCT)
Prior art keywords
imaging data
nerve
tract
data
program
Prior art date
Application number
PCT/US2001/004479
Other languages
French (fr)
Inventor
David Roberts
Original Assignee
David Roberts
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 David Roberts filed Critical David Roberts
Priority to US10/467,789 priority Critical patent/US7672708B2/en
Priority to CA002437822A priority patent/CA2437822A1/en
Priority to EP01910573A priority patent/EP1361820B1/en
Priority to PCT/US2001/004479 priority patent/WO2002064034A1/en
Priority to IL15734901A priority patent/IL157349A0/en
Priority to CNA018230741A priority patent/CN1505490A/en
Priority to JP2002563835A priority patent/JP2004528873A/en
Publication of WO2002064034A1 publication Critical patent/WO2002064034A1/en
Priority to IL157349A priority patent/IL157349A/en
Priority to AU2007203153A priority patent/AU2007203153B2/en
Priority to US12/716,069 priority patent/US8195270B2/en

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/149Segmentation; Edge detection involving deformable models, e.g. active contour models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20092Interactive image processing based on input by user
    • G06T2207/20096Interactive definition of curve of interest
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30101Blood vessel; Artery; Vein; Vascular

Definitions

  • the present invention relates to the non- invasive mag- netic imaging of anatomic tissue structures in isolation from surrounding tissues, and, in particular, to the non-invasive magnetic imaging of nerve tissues.
  • the present invention more specifically relates to the non-invasive quantitative measurement of damage to nerve tissue, such as the damage that occurs upon compression.
  • Peripheral nerves are compound structures. Each nerve trunk consists of a very large number of nerve "fibers” (neurites) that are protoplasmic extensions of neuronal cell bodies located in the dorsal root ganglia (sensory neurons) and ventral horn of the spinal cord (motor neurons) . All nerve fibers (except fine pain fibers) are surrounded by a fatty sheath of myelin which serves to enhance nerve impulse conduction (neural flow) . The nerve fibers are grouped into fasciculi, bundles of fibers bound together by a fine connective tissue membrane (epineureum) . The bundles of fasciculi comprising the nerve trunk are bound by a thicker membrane, the perineureum. The nerve trunk may be surrounded by other layers of protective connective tissue or fat.
  • nerve "fibers” nerves
  • All nerve fibers are surrounded by a fatty sheath of myelin which serves to enhance nerve impulse conduction (neural flow) .
  • the nerve fibers are
  • Neuropathy in peripheral nerves may occur through the agency of traumatic distortion, including compression, stretch and breakage, and/or associated ischemia.
  • the trauma may be acute, especially when the surrounding tissues are injured, or it may be chronic, and due to pathologic mechanical relations between the nerve and surrounding tissues : regional compression is considered to be the principal etiology in entrapment syndromes. There is rapid recovery from transient compression, but prolonged compression extends the recovery period because the nerve fibers are damaged. Remye- lination and re-establishment of conduction must take place. Chronic and repetitive compression prevents recovery and results in severe neurologic impairment .
  • Nerve trunks become thinned at the site of the compression, but may swell proximally. Characteristically conduction is slowed and conduction blocks occur as long segments of myelin are damaged.
  • the nerve trunk proximally will contain a full complement of yelinated fibers within its fasciculi; but at the level of compression and distally the number of myelinated fibers may be reduced significantly.
  • the relative density of the myelinated fibers between proximal and distal portions of the nerve is a measure of loss of neurophysiologic function in its distal distribution.
  • Classic sites for chronic nerve compression include the brachial plexus - especially where the cords of the plexus pass between the rib cage and clavicle, the elbow - where the ulna nerve winds around the olecranon process of the ulna bone, the carpal tunnel - where the median nerve lies between deep flexor tendons and the flexor retinaculum, the sacro-iliac region - where the lumbo-sacral trunk and sacral plexus lie on the bony surfaces of the pelvis, the gluteal region - where the sciatic nerve passes between the powerful gluteal musculature, the neck of the fibula - where the common pero- neal nerve winds about the external surface of the bone, and the tarsal canal - where the tibial nerve winds beneath the tarsal bones to enter the sole of the foot .
  • Some nerves are vulnerable to acute trauma as a result of their intimate relationship to bones subject to fracture. These include that portion of the brachial plexus lying beneath the clavicle, the radial nerve - which winds about the shaft of the humerus, and the common peroneal nerve, which winds about the neck of the fibula.
  • Carpal Tunnel Syndrome is defined as a pathologic condition in which the tunnel size is reduced, causing pressure on the median nerve with resultant pain and slight sensory impairment in the digits supplied by the nerve, and sometimes accompanied by slight wasting of the thenar muscles.
  • the carpal groove is a deep concavity on the palmar surface of the carpus formed by lateral and medial projections of the carpal bones.
  • the medial projection is formed by the pisiform bone ' and the hamulus - a hook-like projection of the hamate bone.
  • the lateral projection is formed by the tubercles of the scaphoid and trapezium bones .
  • the carpal groove is converted into an osteofibrous carpal tunnel by a strong fibrous retinaculum attached to the lateral and medial margins of the carpal groove.
  • the carpal tunnel transmits the tendons of the deep and superficial digital flexor muscles which in this region form a firm bundle .
  • the median nerve lies within the tunnel between the flexor retinaculum and the flexor tendon bundle. It emerges from behind the flexor digitorum superficialis distally in the forearm about 5 cm proximal to the flexor retinaculum of the carpus. It lies superficial to the deep digital flexor tendons and between the tendons of the superficial flexor muscle and the flexor carpi radialis. The nerve then passes deep to the flexor retinaculum of the wrist where normally it lies superficial to the stout bundle of flexor tendons.
  • Overt morphology associated with CTS includes carpal dislocations, carpal arthritis and tenosynovitis of the long flexor tendons. These conditions clearly restrict the space within the tunnel. There is often no apparent morphologic cause for the syndrome, but as complete division of the flexor retinaculum results in relief from symptoms in most cases it is assumed that mechanical compression that is not obvious from the gross morphology of the tunnel or its contents causes ischemia or directly compresses the median nerve.
  • Clinical symptoms includes weakening of the thenar muscles, parasthesias of the lateral palm, thumb and first two fingers. Other symptoms may include stiffness of the wrist and some pain.
  • Clinical examination consists of identifying those symptoms known to be associated with nerve damage.
  • the anatomic distribution of the branches of peripheral nerves is constant and well known.
  • the distribution of abnormal sensation or inability to utilize one or more muscles is a s ⁇ re indicator of which nerve ⁇ is involved.
  • Electromyography is used as a supplement to clinical examination.
  • the technique uses needle electrodes inserted into suspect muscle.
  • the intensity of the potentials measured is an indirect measure of the level of muscle activity.
  • Neural flow stimulation also involves the use of needle electrodes.
  • the technique compares input and output values of electric potentials . None of these diagnostic techniques definitively identifies the focus of the neuropathy. In most cases, if the examination and electromyographic results are typical of those associated with a particular syndrome, the focus of the neu- ropathy can be estimated with some certainty based on experience; but false positive diagnoses are relatively common and may lead to unnecessary surgery or incomplete relief of the symptoms.
  • electro yography is a relatively complex technique requiring some expertise to obtain a satis- factory diagnosis.
  • the use of needle electrodes is invasive, and may prove painful. This can lead to false readings or patient rejection of the study.
  • post-scan data processing may be employed to resolve the problem of obtaining physiological information with significant clinical diagnostic value from magnetic imaging data.
  • imaging data By first segmenting the imaging data to isolate anatomic tissue structures from surrounding tissues, high resolution three-dimensional images are obtained in which tissue microstructures are visible and provide clinically significant information relating to phys- iology.
  • tissue microstructures are visible and provide clinically significant information relating to phys- iology.
  • neurophysiological information with significant clinical diagnostic value relating to the efficacy of nerve flow.
  • a method for the non-invasive imaging 5 of an anatomic tissue structure in isolation from surrounding tissues including the steps of:
  • the imaging data my be acquired by either magnetic resonance or magnetic transfer imaging.
  • Applications for the method of the present invention are not limited to the imag-
  • the method of the present invention is particularly well suited for the imaging of vascular tissue and tumor tissue.
  • the present invention makes possible high resolution images of tumor tissue boundaries that assist the surgeon in
  • an apparatus for preforming the method of the invention.
  • An apparatus in accordance with the present invention includes:
  • an input device for receiving magnetic imaging data from a patient
  • a processing unit having:
  • a processor for sequentially executing the first, second and third programs
  • an output device for displaying the image or a storage device for storing the image an output device for displaying the image or a storage device for storing the image.
  • the present invention includes apparatuses for calculating the presence or loss of myelin in a nerve trunk or tract .
  • Such apparatuses include devices that measure the myelin without first constructing an image of the nerve trunk or tract.
  • An apparatus in accordance with this aspect of the present invention includes:
  • an input device for receiving imaging data from a patient of a nerve trunk or tract and the surrounding tissues
  • a processing unit having:
  • a processor for sequentially executing the first, second and third programs
  • processor executes :
  • the second program using the nerve trunk or tract imaging data isolated by execution of the first program
  • the third program using the fasciculi data population obtained by execution of the second program; and an output device for displaying the myelin calculation or a storage device for storing the myelin calculation.
  • the method and apparatus of the present invention are particularly useful for segmenting nerve tissue imaging data, enabling three-dimensional visualization of individual nerve fasciculi, and quantification of the percentage of fascicular tissue present in cross-sections of the nerve trunk.
  • the numerical values obtained are represented in bivariate plots of percentage and cross-sectional location and directly indicate to the physician those parts of the nerve trunk where myelin content is reduced and neural flow is impaired, and the extent of the degeneration.
  • a method for the imaging of a nerve trunk or tract suspected of demyelination including:
  • the present invention incorporates the discovery that the imaging data for the fasciculi can be used to calculate the presence or loss of myelin in a nerve trunk or tract.
  • preferred methods according to this aspect of the present ' invention further include the steps of calculating from the imaging data for the fasciculi, the presence or loss of myelin in the nerve trunk or tract suspected of demyelin- ation, and storing or displaying the myelin calculation.
  • the presence or loss of myelin in a nerve trunk or tract may be calculated from fasciculi imaging data independently of the construction of the nerve trunk image.
  • the present invention therefore also includes a method for the quantitative measurement of a nerve trunk or tract suspected of demyelina- tion, including the steps of:
  • a method for the imaging of a median nerve suspected of demyelination within one of a patient's carpal tunnels includ- ing the steps of: receiving from an input device, magnetic imaging data from a patient including at least a portion of the median nerve passing through a carpel tunnel ,-
  • This method can also be employed to calculate the presence or loss of myelin in the median nerve within a carpal tunnel from the fasciculi imaging data.
  • Methods in accord- ance with this aspect of the present invention further include the steps of calculating from the fasciculi imaging the presence or loss of median nerve myelin within the carpal tunnel, and storing or displaying the myelin calculation.
  • the presence or loss of myelin in the median nerve ith- in a carpal tunnel can also be calculated from the fasciculi imaging data independently of the construction of a median nerve image .
  • the present invention therefore also includes a "' method ' " for the quantitative"measurement of a' median nerve suspected of demyelination within one of a patient's carpal tunnels, including the steps of:
  • magnetic imaging data from a patient including at least a portion of the median nerve passing through a carpel tunnel ;
  • segmenting the imaging data to isolate the median nerve imaging data from any imaging data for surrounding tissues separating the median nerve imaging data into data populations corresponding to the fasciculi and other micro- structures of the median nerve;
  • FIG. 1 is a block diagram of the components of an exemplary apparatus in accordance with the present invention ;
  • FIG. 2 is a flow chart illustrating an exemplary set of steps in accordance with the method of the present invention
  • FIG. 3 is a prior art axial cross-section of a wrist, including the median nerve, obtained by magnetic resonance imaging;
  • FIG. 4 is a magnetic resonance image of a healthy median nerve within a carpal tunnel of a patient, constructed by an apparatus of the present invention employing the method of the present invention.
  • FIG. 5 is a magnetic resonance image of a median nerve within a carpal tunnel of a patient suffering from CTS, constructed by the same apparatus and method employed to construct the image of FIG. 4.
  • an imaging device 10 which includes a processing unit 12 having processor 14, memory 16 and stored programs 18, including post-scan data processing programs 20.
  • the processing device also includes an input device 26 and an output device 28.
  • the present invention is based on the discovery that by first segmenting imaging data to isolate anatomic tissue structures from surronding tissues before histographically resegmenting the imaging data, it is possible to construct images of the tissue microstructures having exceptional clarity.
  • the new methods and apparatus disclosed herein are particularly adapted for use in the diagnosis of neuropathy, especially the neuropathy resulting from compression and the associated demyelination, such as occurs with CTS .
  • the processing unit 12 is an R1000 workstation from SGI having about 128 MB associated RAM memory and a 4 GB hard or fixed drive memory 16.
  • the processor operates using a UNIX operating system to run application software 18 implementing the disclosed imaging device 10 and method.
  • processor is merely exemplary. Any computer or processing system having sufficient processing ability and memory may be used.
  • Imaging data acquired by magnetic resonance or magnetic transfer imaging is provided by imaging source 24 through an input device 26 to the processor 14.
  • the data source may be any source for providing magnetic resonance or magnetic transfer imaging data to a computer, such as a floppy disk containing an imaging data file, or the processor of a magnetic resonance or magnetic transfer imaging device (not shown), directly linked to a input port of processor 14, or by means of a modem interface.
  • the input device 26 may be a floppy disk drive, modem, or other data reading or data input device, and also includes an interactive, semi-automatic or automatic contouring device for segmenting of anatomic structures.
  • anatomic structures are defined as the functional substructures within a gross organ, such as the functional microstructures within a nerve trunk, i.e., fasciculi, and the like.
  • Input device 26 may also include a keyboard or similar input device for supplying commands to processor 14.
  • the command-supplying input device may also be provided as a second input device (not shown) .
  • the processor 14 receives commands and imaging data from the imaging data source through the input device (s) .
  • the imaging data my be stored in memory 16 for further processing or used by the processor 14 to construct an image of a tissue target.
  • the constructed image is sent to an output device 28 such as a visual image display device, for example, a video display terminal.
  • a video terminal for displaying the original magnetic resonance or magnetic transfer image may also be included.
  • the constructed image may also be displayed numerically or graphically through an appropriate output device.
  • the imaging device 10 performs the applications programs described herein in conjunction with Figure 2, which axe implemented from compiled source code in the "C" programming language.
  • the imaging device 10 includes conventional imaging software, such as software for the display of original magnetic images.
  • the imaging device also includes post-scan data processing programs, such as
  • clipping software for discarding irrelevant data
  • software for operating interactive, semi-automatic or automatic contouring devices for segmentation of anatomic structures
  • the imaging device also includes post-scan data process- ing software for the creation of binary images of segmented structures and binary slices and software for the creation of three-dimensional binary masks and software for the segmentation of gray-level data representing the segmented structure via the binary mask.
  • Post-scan data processing software for the histographic analysis of the gray- level data is also included, which performs single or iterative segmentation of the gray-level anatomic tissue structure into the intrinsic tissue microstructures, as well as software for the definition of material opacity, emission and surface reflection.
  • the imaging device post-scan data processing software also includes software for the creation of three-dimensional images of individual intrinsic tissue microstructures, or compound images of complete tissue structures, and software for image smoothing (aliasing) .
  • Post-scan data processing software is also included for the calculation of the central trunk (tract) axis (axes) , cross-sections and tissue percentages and the graphical display of values.
  • Figure 2 depicts an applications program for the imaging of a nerve trunk.
  • magnetic imaging data is acqu- ired in 1 mm slices along the nerve trunk axial plane.
  • the data is ' transferred to a processor for post scan processing.
  • gray-level images of the magnetic imaging slices are constructed by the image display software.
  • step 34 the interactive clipping software identifies the nerve trunk in the gray-level slices and in step 36 clip . the first nerve trunk slice to eliminate the surrounding tissue imaging data. The rest of the slice stack is then automatically clipped, with modification of the clipping boundary as required.
  • step 38 the semi-interactive segmentation ⁇ software identifies the nerve trunk boundary to the computer in the initial one or two slices, after which the software automatically segments the remaining slices and then checks and modifies the identified boundaries.
  • step 40 the binary masking software creates two- dimensional binary masks for the segmented nerve boundaries.
  • step 42 the two-dimensional binary masks are used by the gray-level data segmentation software to segment the gray- level data representing the nerve trunk from the original magnetic resonance or magnetic transfer image slices.
  • step 44 the histographic analysis software resegments the nerve trunk gray-level data for the segmented nerve trunk into nerve trunk microstructure data populations, including the fascicular bundles, the epineurium and inter-fascicular fascia. Interpolation of the resegmented data populations is performed in step 46. Software for the definition of material opacity, emission and surface reflection adjusts these features for individual microstructures to enhance their visualization within compound tissue images .
  • Still and movie images are then constructed by the three-dimensional imaging software in step 48. Both individual nerve tissue microstructure images and compound images of the nerve trunk composed of the individual nerve tissue microstructures are constructed.
  • step 50 images are displayed, either as printed i a- ges or by viewing software as video displays.
  • the tissue calculation software determines in step 52, from the fasciculi imaging data, the amount of fascicular material in selected cross-sections, the values for which are presented graphically in step 54 by the display software.
  • Nerve trunks are readily visualized in high resolution TI or T2 weighted magnetic resonance (MR) images.
  • MR magnetic resonance
  • the fascicles appear as bright "spots" interspersed in a darker gray matrix which represent interstitial tissues, surrounded by a dark border of variable thickness that represents the perineureum and several layers of thin supporting fascia.
  • the relatively high amplitude signal from the fasciculi relates to the presence of myelin.
  • Myelin is composed of lipids and protein, as in the plasma membranes of other cells, but the percentage of lipid is higher; and in addition to the usual suite of lipids found in plasma membranes, myelin contains a characteristic glycolipid. Myelin also contains up to 30% bound water.
  • fasciculi can be partially detected by direct visual inspection of cross-sectional MR images. But determination of the number or percentage of fasciculi present in a cross-section, comparisons between different sections is not satisfactory via visual inspection alone. It is often not possible to resolve darker fasciculi and small changes in signal amplitude may go unobserved.
  • the method of the present invention is thus used to generate three-dimensional images of the fasciculi, which can be visualized in a static format, or which may be rotated about the (approximate) long axis of the trunk.
  • the latter enhances the visualization process. Because the three-dimensional image can be limited to the fasciculi, the local density of these structures within the trunk, and their brightness (signal amplitude) are now relatively easy to visualize directly.
  • Nerve tissue demyelination is determined by first calculating the volume of fascicular tissue as a percentage of the volume of the gross trunk (tract) structure. Next, the area of fascicular tissue as a percentage of the area of the gross trunk (tract) structure as seen in selected cross- sections of the trunk (tract) is calculated. The contiguous cross sectional areas of the trunk (tract) of specified thickness (i.e., 1 mm) perpendicular to the long axis is then calculated.
  • the percentage of volume of fascicular tissue as a percentage of total volume of the trunk (tract) in contiguous cross-sections is calculated next.
  • a graphical display of the percentage volumes obtained vs. location of the sections on the central axis of the trunk is then constructed.
  • the variation of myelin content as indicated by these percentage values is a measure of the neurophysiologic (flow efficiency) status of the trunk (nerve) and thus a diagnostically significant index.
  • This technique enables a comparison of the percentage of fascicular tissue present, for example, in the forearm (where entrapment and compression is unlikely) , with the value for fascicular tissue within the carpal tunnel, and in regions of the trunk distal to the carpal tunnel . It is also possible to compare fascicular percentages in homologous regions from different individuals (compare between patients) .
  • the method of the present invention thus allows for the fascicular percentages to be calculated for cross-sectional areas of a nerve trunk.
  • fascicular percentages for 1 mm thick contiguous sections normal to the long axis of the trunk may be calculated, and then represented as a bivariate plot of percentage of fascicular material present vs. location within the trunk.
  • This graphical representation provides the CTS diagnostician, with a simple and direct means of determining in which precise regions of the median nerve trunk a loss of fascicular material has occurred, and the amount of fascicular tissue lost .
  • Magnetic Transfer Imaging is a new technique that has two advantages over conventional MR imaging when applied to analysis of nerve trunks. First, the in-plane resolution is improved. Second, the data values obtained with this technique are not scanner dependant and consequently direct inter-patient comparisons are possible.
  • MTI is particularly well-suited for use with the methods and apparatus of the present invention.
  • the method and apparatus of the present invention are also well-suited for the diagnosis of nerve pathology.
  • a clinical examination is typically first performed to determine the likely location of nerve pathology.
  • a high-resolu- tion magnetic image of the region is then obtained.
  • the MR data comprises thin (1 mm) contiguous slices in a plane which is more or less normal to the axis of the nerve.
  • High resolution images can be obtained by using scan times up to ten minutes.
  • Post-scan data processing is then applied to the MR data.
  • MR hard-copy (acetate) is not required.
  • the three- dimensional MR data is loaded into the zhree-dimensional imaging software. Steps in the subsequent data processing are as follows:
  • the original MR slices are displayed.
  • the nerve is identified and image contrast adjusted eo maximize its visualization by the operator. This cor.trast adjustment does not affect the signal differential marking the boundary of the nerve, and is not a factor in subsequent segmentation.
  • the size of the data set is reduced by data clipping to discard all data except that in the immediate vicinity of the nerve. Once a region of interest has been defined, clipping is automatic; thus, the size of the region of interest relates to the position of the nerve on successive images.
  • the original voxels are prismatic; i.e., in the Z direction (thickness of image-slice) they are usually 1 mm.
  • the dimen- sions of X and Y (in the plane of the slice) depend on the imaging protocol. Typically X and Y are 0.5-0.7 mm. Cubile voxels are required for 3 -dimensional rendering; thus, it is necessary to interpolate the data set in the Z direction in order to obtain slices corresponding in thickness to the X,Y dimensions.
  • the interpolation process is not simple data averaging; it utilizes a sophisticated algorithm which identifies the "shape" of regions of consistent signal in contiguous slices and calculates intermediate "shapes" (assuming linearity in shape change between slices) .
  • Shape based in- terpolation has been shown to produce highly accurate information in the interpolated slices. Interpolation is usually performed after clipping and segmentation on the binary image produced. This reduces the time required to perform interactive contouring during the process of segmentation.
  • the operator With interactive segmentation the operator traces the boundary of the nerve trunk section, recording coordinates on the boundary in so doing.
  • the boundary of the nerve as seen in cross-section is identified by the operator and indicated to the computer.
  • the computer is thus able to identify subsequently contiguous voxels with the same boundary signal characteristics. Once the boundary has been identified, the operator can indicate a starting point on it, then move the cursor to a next point some distance away from the first.
  • the computer is able to track the boundary between the two points using the boundary information provided, and draws a colored line to mark it. Because the boundary signal varies, and the presence of adjacent non-boundary pixels with similar signal characteristics, the boundary selected automatically may vary from that which the operator judges to be correct. In this case the boundary is redefined at the point of divergence (i.e., the boundary-contour is corrected) and the technique of indicating points at intervals along the boundary is continued until the boundary is closed. Boundary identification during contouring is used to refine the suite of signal characteristics used by the computer to determine the boundary automatically. After one or two boundaries have been contoured it is usually only necessary to indicate two or three points to obtain a complete, accurate boundary.
  • This technique can be extended to automatic identifica- tion of the boundary in subsequent slices, which speeds up the procedure.
  • the operator examines the automatically generated boundary contours and can correct any errors.
  • the continually refined boundary characteristics are stored by the computer and can be utilized to segment other data sets obtained via similar imaging protocols.
  • Live-path contouring is similar to live-wire contouring, but the computer only seeks boundary characteristics coinciding with those defined for it in a "path" of limited width on either side of the cursor. It is possible to run the cursor rapidly around the region of interest, without following the boundary exactly, to describe a contour within the limits of the defined path, for a very rapid method of contouring.
  • Live-wire and live-path contouring techniques are normally carried out on all the original MR slices. However, these techniques may be developed so that contours are only defined interactively on the first two or three slices, after which the computer will search for the corresponding boundary in the remaining slices and contour them automatically.
  • the segmentation techniques output the surface boundary of the region of interest by a binary separation of the three dimensional data based on those voxels within the boundary and those without it. Three-dimensional imaging of the bin- ary data therefore yields a "solid" object.
  • the boundary surface can be used as a "mask” to extract MR data within the region of interest from the original three-dimensional MR data matrix.
  • the extracted MR data set is referred to as a gray-level image.
  • the gray-level image thus contains information about the tissue microstructures contained within the nerve trunk.
  • the gray-level values for nerve fasciculi, matrix, and perineur- eum comprise distinct and overlapping populations. In the region of signal overlap the percentage of signal belonging to each material can be ranges of signal values assigned accordingly. This allows the different tissue types to be resegmented automatically based on this histographic data.
  • the histographic resegmentation may be performed iteratively in order to resolve any signal overlap to obtain the best visualization of the nerve tissue microstructures (i.e., visualization of the very distinctive fasciculi) .
  • Opacity, emission and surface strength parameters optimal for visualization of the nerve microstructures are then determined by creation of three-dimensional images and examination of the anatomy. Once values have been determined for these parameters, they may be used for more automated imaging in subsequent data sets. Improvement of visualization by adjusting these parameters is a distinct process from segmentation and does not select the tissue concerned from its surrounding matrix.
  • the output is visual and numerical.
  • Visual results con- sist of three-dimensional renderings of the trunk, including a compound image of the gross structure, which includes the perineurium, fascicles and interfascicular tissue.
  • the opacity of the perineurium is set so that this tissue is translucent . This enables the content of the nerve trunk to be viewed in the context of the perineurium.
  • the opacity of the interfascicular tissue is set to be greater than that of the perineurium, but less than that of the fascicles. The result is to display the microstructures of the nerve in their correct anatomic relationships.
  • the image is rotated about its long axis and recorded as multiple "frames" which are stored in a "cine” file. This is performed using conventional MR imaging software. Images are also created containing each of the material sets separately, and cine files created for each. These images allow direct examination of each microstructure separately.
  • a compound image consisting of the perineureum and fas- ciculae is also created and rendered in a cine file. This image allows the texture and thickness of the perineureum to be examined in the context of the fasciculi, and regions of membrane thickening or thinning compared with regions of MR signal fluctuation in the fasciculi.
  • the numerical output includes the cross-sectional area value for the nerve trunk normal to its long axis, calculated continuously throughout its length and displayed in a bivari- ate plot of cross-sectional area vs. section location. A similar bivariate plot is calculated to show the density of fascicular material within the cross-section vs. location.
  • the parameters measured by the method and apparatus of the present invention for nerve tissue include the general morphology of the trunk (tract) , the cross-sectional morphology, the size and distribution of the nerve fasciculi and the density of the fasciculi within cross-sections of the nerve trunk. These parameters are associated with physiologic changes accompanying nerve damage and dysfunction.
  • Gross morphology and cross-sectional morphology are phy- sical reflections of the pressure environment. It is well known that in regions of pressure a peripheral nerve reorganizes by becoming flattened in a plane normal to the compress- ive force. The amount of flattening almost certainly relates to altered function, but the sensitivity of the relationship may not be adequate for the purpose of early diagnosis.
  • cross-sectional fascicular density is a measure of MR signal from the "fascicular material" present in the section.
  • the MR signal may be expected to vary with the his- tochemistry of the nerve fibers (neurites) and their interstitial environment. This includes changes in the structure and chemistry of the myelin sheaths (composed of complex lipids and a high percentage of bound water) and in the ionic environment of the neurite and interstitial tissues.
  • variation in the MR signal is an indicator of functional changes in the neurites, including those due to trauma.
  • Signal variation can be detected in two ways. One way is by direct examination of the final image. Regions of nerve trauma are likely to be represented by localized disappearance of the fasciculi. The extent of fascicular disappearance is an indication of the severity of the nerve damage.
  • the MR signal variation can also be quantified by measuring the density of the fascicular material in cross-sections of the nerve trunk and plotting this data against the locations of the cross-sections.
  • Neural flow studies can be made on the median nerve in the region of the carpal tunnel because of the very superfic- ial placement of the trunk in this region.
  • neuro- pathies occur in other regions, both within and outside the central nervous system where neural flow studies are impractical because of the depth of the trunk.
  • the methods and apparatus of the present invention provide an innovative means of diagnosing and assessing neuropathy in such deep peripheral nerves, and in portions of the central nervous system.
  • the neuropathy need not be associated with compression or other trauma, but may also result from other causes, such as multiple sclerosis and related diseases.
  • multiple sclerosis and related diseases such as multiple sclerosis and related diseases.
  • practitioners were forced to rely upon the subjective assessment of the patient. This was particularly inaccurate in most circumstances, in which patients want to believe that a treatment is working and relate their well-being from this perspective.
  • the methods and apparatus of the present invention may also be used to assess other tissues that are inaccessible because of their depth within the body, including vascular tissues and tumor tissues.
  • the present invention makes it possible to visualize vascular blockage, as well as tumor tissue boundaries.
  • the ability to visualize tumor tissue boundaries makes the present invention a particularly useful tool for the cancer surgeon.
  • FIG. 3 A prior art axial cross-section of a wrist obtained by magnetic resonance imaging is shown in Fig. 3.
  • the cross- sectional image includes the median nerve and associated tendons and carpal bones . This Figure illustrates how such images provided little information of clinical diagnostic value relating to neurophysiology.
  • FIG. 4 A magnetic resonance image of a healthy median nerve trunk constructed by the method and apparatus of the present invention is shown in Fig. 4.
  • the white streaks are bundles of myelinated nerve fibers, i.e., fasciculi. This is in contrast to the magnetic resonance image of a median nerve trunk of a patient afflicted with CTS, shown in Fig. 5, and also constructed by the method and apparatus of the present invention.
  • the nerve fiber demyelination is dramatically evident.
  • MR signal amplitude was reduced where the nerve trunk passes beneath the flexor retinaculum. Reduction of signal amplitude is modest in subjects with no clinical symptoms of CTS, but becomes increasingly significant with CTS symptoms, and is pronounced in subjects diagnosed with severe CTS.
  • the four cadaver wrists were also subjected to MR scanning. Evidence of demyelination in these wrists was partially confirmed by extirpation of the median nerve from the carpal tunnel and histologic sectioning and staining.

Abstract

A method for the non-invasive imaging of an anatomic tissue structure in isolation from surrounding tissues, including: receiving from an input device magnetic imaging data from a patient of the anatomic tissue structure and surrounding tissues; segmenting the imaging data to isolate the anatomic tissue structure imaging data from the imaging data for the surrounding tissues; separating the anatomic tissue structure imaging data into data populations corresponding to tissue microstructures; constructing an image from the imaging data for at least one of the tissue microstructures; and storing or displaying the image. An apparatus embodying the disclosed method is also described, as well as a method for the quantitative measurement of a nerve tissue suspected of demyelination.

Description

METHOD AND APPARATUS FOR THE NON- INVASIVE IMAGING OF
ANATOMIC TISSUE STRUCTURES
Background Of The Invention
The present invention relates to the non- invasive mag- netic imaging of anatomic tissue structures in isolation from surrounding tissues, and, in particular, to the non-invasive magnetic imaging of nerve tissues. The present invention more specifically relates to the non-invasive quantitative measurement of damage to nerve tissue, such as the damage that occurs upon compression.
Clinically it is possible to estimate which neural pathway is involved in trauma by determining the end result loss of motor function and either abberant sensory function or complete loss of it. It is very difficult to determine where the lesion is and to what extent it has caused deterioration of the nerves involved.
Peripheral nerves are compound structures. Each nerve trunk consists of a very large number of nerve "fibers" (neurites) that are protoplasmic extensions of neuronal cell bodies located in the dorsal root ganglia (sensory neurons) and ventral horn of the spinal cord (motor neurons) . All nerve fibers (except fine pain fibers) are surrounded by a fatty sheath of myelin which serves to enhance nerve impulse conduction (neural flow) . The nerve fibers are grouped into fasciculi, bundles of fibers bound together by a fine connective tissue membrane (epineureum) . The bundles of fasciculi comprising the nerve trunk are bound by a thicker membrane, the perineureum. The nerve trunk may be surrounded by other layers of protective connective tissue or fat.
Neuropathy in peripheral nerves may occur through the agency of traumatic distortion, including compression, stretch and breakage, and/or associated ischemia. The trauma may be acute, especially when the surrounding tissues are injured, or it may be chronic, and due to pathologic mechanical relations between the nerve and surrounding tissues : regional compression is considered to be the principal etiology in entrapment syndromes. There is rapid recovery from transient compression, but prolonged compression extends the recovery period because the nerve fibers are damaged. Remye- lination and re-establishment of conduction must take place. Chronic and repetitive compression prevents recovery and results in severe neurologic impairment .
In chronic and repetitive compressions, which are more common than acute compressions, the damage is more extensive. Nerve trunks become thinned at the site of the compression, but may swell proximally. Characteristically conduction is slowed and conduction blocks occur as long segments of myelin are damaged. Depending on the severity of the chronic compression, the nerve trunk proximally will contain a full complement of yelinated fibers within its fasciculi; but at the level of compression and distally the number of myelinated fibers may be reduced significantly. The relative density of the myelinated fibers between proximal and distal portions of the nerve is a measure of loss of neurophysiologic function in its distal distribution.
Although nerve damage can occur anywhere within the body, some nerves are predisposed to chronic damage because of their structural relationships with other tissues. Classic sites for chronic nerve compression include the brachial plexus - especially where the cords of the plexus pass between the rib cage and clavicle, the elbow - where the ulna nerve winds around the olecranon process of the ulna bone, the carpal tunnel - where the median nerve lies between deep flexor tendons and the flexor retinaculum, the sacro-iliac region - where the lumbo-sacral trunk and sacral plexus lie on the bony surfaces of the pelvis, the gluteal region - where the sciatic nerve passes between the powerful gluteal musculature, the neck of the fibula - where the common pero- neal nerve winds about the external surface of the bone, and the tarsal canal - where the tibial nerve winds beneath the tarsal bones to enter the sole of the foot .
Some nerves are vulnerable to acute trauma as a result of their intimate relationship to bones subject to fracture. These include that portion of the brachial plexus lying beneath the clavicle, the radial nerve - which winds about the shaft of the humerus, and the common peroneal nerve, which winds about the neck of the fibula.
Carpal Tunnel Syndrome (CTS) is defined as a pathologic condition in which the tunnel size is reduced, causing pressure on the median nerve with resultant pain and slight sensory impairment in the digits supplied by the nerve, and sometimes accompanied by slight wasting of the thenar muscles. The carpal groove is a deep concavity on the palmar surface of the carpus formed by lateral and medial projections of the carpal bones. The medial projection is formed by the pisiform bone' and the hamulus - a hook-like projection of the hamate bone. The lateral projection is formed by the tubercles of the scaphoid and trapezium bones . The carpal groove is converted into an osteofibrous carpal tunnel by a strong fibrous retinaculum attached to the lateral and medial margins of the carpal groove. The carpal tunnel transmits the tendons of the deep and superficial digital flexor muscles which in this region form a firm bundle .
The median nerve lies within the tunnel between the flexor retinaculum and the flexor tendon bundle. It emerges from behind the flexor digitorum superficialis distally in the forearm about 5 cm proximal to the flexor retinaculum of the carpus. It lies superficial to the deep digital flexor tendons and between the tendons of the superficial flexor muscle and the flexor carpi radialis. The nerve then passes deep to the flexor retinaculum of the wrist where normally it lies superficial to the stout bundle of flexor tendons.
Overt morphology associated with CTS includes carpal dislocations, carpal arthritis and tenosynovitis of the long flexor tendons. These conditions clearly restrict the space within the tunnel. There is often no apparent morphologic cause for the syndrome, but as complete division of the flexor retinaculum results in relief from symptoms in most cases it is assumed that mechanical compression that is not obvious from the gross morphology of the tunnel or its contents causes ischemia or directly compresses the median nerve.
Clinical symptoms includes weakening of the thenar muscles, parasthesias of the lateral palm, thumb and first two fingers. Other symptoms may include stiffness of the wrist and some pain.
Clinical examination consists of identifying those symptoms known to be associated with nerve damage. The anatomic distribution of the branches of peripheral nerves is constant and well known. The distribution of abnormal sensation or inability to utilize one or more muscles is a sιre indicator of which nerve~ is involved.
Electromyography is used as a supplement to clinical examination. The technique uses needle electrodes inserted into suspect muscle. The intensity of the potentials measured is an indirect measure of the level of muscle activity. Neural flow stimulation also involves the use of needle electrodes. The technique compares input and output values of electric potentials . None of these diagnostic techniques definitively identifies the focus of the neuropathy. In most cases, if the examination and electromyographic results are typical of those associated with a particular syndrome, the focus of the neu- ropathy can be estimated with some certainty based on experience; but false positive diagnoses are relatively common and may lead to unnecessary surgery or incomplete relief of the symptoms. In addition, electro yography is a relatively complex technique requiring some expertise to obtain a satis- factory diagnosis. Furthermore, the use of needle electrodes is invasive, and may prove painful. This can lead to false readings or patient rejection of the study.
A further problem is experienced with clinical examinations. Numerous patients attempt to simulate symptomatology in order to qualify for insurance compensation. Identifying these patients using an essentially subjective examination is difficult and time consuming for the diagnostician.
Despite conflicting reports in the literature, there is no convincing evidence that nerve morphology viewed in unas- sisted two-dimensional magnetic images is a clinically practical means of obtaining information about histological degeneration and impaired neural flow. At pι~esent, CTS patients are not referred for magnetic imaging scans as a diagnostic procedure .
Summary Of The Invention
It has now been discovered that post-scan data processing may be employed to resolve the problem of obtaining physiological information with significant clinical diagnostic value from magnetic imaging data. By first segmenting the imaging data to isolate anatomic tissue structures from surrounding tissues, high resolution three-dimensional images are obtained in which tissue microstructures are visible and provide clinically significant information relating to phys- iology. When imaging nerve tissues, it is possible to obtain neurophysiological information with significant clinical diagnostic value relating to the efficacy of nerve flow.
Therefore, according to one aspect of the present invention, a method is provided for the non-invasive imaging 5 of an anatomic tissue structure in isolation from surrounding tissues, including the steps of:
I
receiving from an input device, magnetic imaging data from a patient of an anatomic tissue structure and the su rounding tissues;
10 segmenting the imaging data to isolate the anatomic tissue structure imaging data from the imaging data for the surrounding tissues ;
separating the anatomic tissue structure imaging data 15 into data populations corresponding to tissue microstructures;
constructing an image from imaging data for at least one of the tissue microstructures; and
storing or displaying the image.
20 The imaging data my be acquired by either magnetic resonance or magnetic transfer imaging. Applications for the method of the present invention are not limited to the imag-
""""" ϊng of nerve tissue"""but'"appTy~tO" 'e'ssential y an tissue" — structure for which it is difficult to obtain a high resolu-
25 tion image that is not obscured by surrounding tissues.
Thus, the method of the present invention is particularly well suited for the imaging of vascular tissue and tumor tissue. The present invention makes possible high resolution images of tumor tissue boundaries that assist the surgeon in
30 achieving the goal of complete tumor removal.
Therefore, according to another aspect of the present invention, an apparatus is provided for preforming the method of the invention. An apparatus in accordance with the present invention includes:
an input device for receiving magnetic imaging data from a patient;
a processing unit having:
a first program for segmenting the magnetic imaging data to isolate imaging data for an anatomic tissue structure from imaging data for surrounding tissues ;
a second program for separating the anatomic tissue structure imaging data into data populations corresponding to tissue microstructures;
a third program for constructing an image from imaging data for at least one of the microstructures; and
a processor for sequentially executing the first, second and third programs;
wherein the processor executes:
the first program using the imaging data from the input device,-
the second program using the anatomic tissue ' ™ " """structure' imagi'n-g"~d'ata* ~isrϊated~"by""exe~cαt:iOn"Ox the"" first program; and
the third program using the data populations obtained by execution of the second program; and
an output device for displaying the image or a storage device for storing the image.
The present invention includes apparatuses for calculating the presence or loss of myelin in a nerve trunk or tract . Such apparatuses include devices that measure the myelin without first constructing an image of the nerve trunk or tract. An apparatus in accordance with this aspect of the present invention includes:
an input device for receiving imaging data from a patient of a nerve trunk or tract and the surrounding tissues;
a processing unit having:
a first program for segmenting the magnetic imaging data to isolate imaging data for the nerve trunk or tract from the imaging data for the surrounding tissues;
a second program for separating the nerve trunk or tract imaging data into data populations corresponding to fasciculi and other nerve trunk or tract microstructures;
a third program for calculating the presence or loss of myelin in the nerve trunk or tract from the fasciculi imaging data; and
a processor for sequentially executing the first, second and third programs;
wherein the processor executes :
the first program using imaging"""data' "from"'Ωϊe~~ input device;
the second program using the nerve trunk or tract imaging data isolated by execution of the first program; and
the third program using the fasciculi data population obtained by execution of the second program; and an output device for displaying the myelin calculation or a storage device for storing the myelin calculation.
The method and apparatus of the present invention are particularly useful for segmenting nerve tissue imaging data, enabling three-dimensional visualization of individual nerve fasciculi, and quantification of the percentage of fascicular tissue present in cross-sections of the nerve trunk. The numerical values obtained are represented in bivariate plots of percentage and cross-sectional location and directly indicate to the physician those parts of the nerve trunk where myelin content is reduced and neural flow is impaired, and the extent of the degeneration.
Therefore, according to another aspect of the present invention, a method is provided for the imaging of a nerve trunk or tract suspected of demyelination, including:
receiving from an input device, magnetic imaging data from a patient of a nerve trunk or tract suspected of demyelination and the surronding tissues;
segmenting the imaging data to separate the imaging data for the nerve trunk or tract from the imaging data for the surrounding tissues;
separating the imaging data for the nerve trunk or tract into data populations corresponding to the fasciculi and other microstructures of the erve trunk or" tract;' '
constructing an image from imaging data for at least the fasciculi of the nerve trunk or tract; and
storing or displaying the image.
The present invention incorporates the discovery that the imaging data for the fasciculi can be used to calculate the presence or loss of myelin in a nerve trunk or tract.
Therefore, preferred methods according to this aspect of the present' invention further include the steps of calculating from the imaging data for the fasciculi, the presence or loss of myelin in the nerve trunk or tract suspected of demyelin- ation, and storing or displaying the myelin calculation.
According to another aspect of the present invention, the presence or loss of myelin in a nerve trunk or tract may be calculated from fasciculi imaging data independently of the construction of the nerve trunk image. The present invention therefore also includes a method for the quantitative measurement of a nerve trunk or tract suspected of demyelina- tion, including the steps of:
receiving from an input device, magnetic imaging data from a patient of a nerve trunk or tract suspected of demye- lination and the surrounding tissues;
segmenting the imaging data to separate the imaging data for the nerve trunk or tract from the imaging data for the surrounding tissues ;
separating the imaging data for the nerve trunk or tract into data populations corresponding to the fasciculi and other microstructures of the nerve trunk or tract;
calculating from the fasciculi imaging data, the presence or loss of myelin in the nerve trunk or tract; and
storing or' displaying the myelin calculation. ' '
The method of the present invention for quantitatively measuring nerve tissue demyelination is particularly useful for the non-invasive diagnosis of CTS . Therefore, according to yet another aspect of the present invention, a method is provided for the imaging of a median nerve suspected of demyelination within one of a patient's carpal tunnels, includ- ing the steps of: receiving from an input device, magnetic imaging data from a patient including at least a portion of the median nerve passing through a carpel tunnel ,-
segmenting the imaging data to isolate the median nerve imaging data from any imaging data for surrounding tissues;
separating the median nerve imaging data into data populations corresponding to the fasciculi and other microstructures of the median nerve;
constructing an image from imaging data for at least the fasciculi of the median nerve; and
storing or displaying the image.
This method can also be employed to calculate the presence or loss of myelin in the median nerve within a carpal tunnel from the fasciculi imaging data. Methods in accord- ance with this aspect of the present invention further include the steps of calculating from the fasciculi imaging the presence or loss of median nerve myelin within the carpal tunnel, and storing or displaying the myelin calculation.
The presence or loss of myelin in the median nerve ith- in a carpal tunnel can also be calculated from the fasciculi imaging data independently of the construction of a median nerve image . The present invention therefore also includes a "'method' "for the quantitative"measurement of a' median nerve suspected of demyelination within one of a patient's carpal tunnels, including the steps of:
receiving from an input device, magnetic imaging data from a patient including at least a portion of the median nerve passing through a carpel tunnel ;
segmenting the imaging data to isolate the median nerve imaging data from any imaging data for surrounding tissues; separating the median nerve imaging data into data populations corresponding to the fasciculi and other micro- structures of the median nerve;
calculating from the imaging data for the median nerve fasciculi, the presence or loss of myelin in the median nerve within the carpal tunnel ; and
storing or displaying the myelin calculation.
A more complete appreciation of the invention and many more other intended advantages can be readily obtained by reference to the following detailed description of the preferred embodiments and claims, taken in conjunction with the accompanying drawings, which disclose the principles of the invention, and the best modes which are presently contemplated for carrying them out .
Brief Description Of The Drawings
FIG. 1 is a block diagram of the components of an exemplary apparatus in accordance with the present invention ;
FIG. 2 is a flow chart illustrating an exemplary set of steps in accordance with the method of the present invention;
FIG. 3 is a prior art axial cross-section of a wrist, including the median nerve, obtained by magnetic resonance imaging;
FIG. 4 is a magnetic resonance image of a healthy median nerve within a carpal tunnel of a patient, constructed by an apparatus of the present invention employing the method of the present invention; and
FIG. 5 is a magnetic resonance image of a median nerve within a carpal tunnel of a patient suffering from CTS, constructed by the same apparatus and method employed to construct the image of FIG. 4. Detailed Description Of The Preferred Embodiments
Referring to the drawings, in which like reference numerals identify similar or identical elements, as shown in Figure 1, the present disclosure describes a method and app- aratus implementing an imaging device 10, which includes a processing unit 12 having processor 14, memory 16 and stored programs 18, including post-scan data processing programs 20. The processing device also includes an input device 26 and an output device 28.
The present invention is based on the discovery that by first segmenting imaging data to isolate anatomic tissue structures from surronding tissues before histographically resegmenting the imaging data, it is possible to construct images of the tissue microstructures having exceptional clarity. The new methods and apparatus disclosed herein are particularly adapted for use in the diagnosis of neuropathy, especially the neuropathy resulting from compression and the associated demyelination, such as occurs with CTS .
' In an exemplary embodiment, the processing unit 12 is an R1000 workstation from SGI having about 128 MB associated RAM memory and a 4 GB hard or fixed drive memory 16. The processor operates using a UNIX operating system to run application software 18 implementing the disclosed imaging device 10 and method.
It should be understood that the above described processor is merely exemplary. Any computer or processing system having sufficient processing ability and memory may be used.
Imaging data acquired by magnetic resonance or magnetic transfer imaging is provided by imaging source 24 through an input device 26 to the processor 14. The data source may be any source for providing magnetic resonance or magnetic transfer imaging data to a computer, such as a floppy disk containing an imaging data file, or the processor of a magnetic resonance or magnetic transfer imaging device (not shown), directly linked to a input port of processor 14, or by means of a modem interface. The input device 26 may be a floppy disk drive, modem, or other data reading or data input device, and also includes an interactive, semi-automatic or automatic contouring device for segmenting of anatomic structures. For purposes of the present invention "anatomic structures" are defined as the functional substructures within a gross organ, such as the functional microstructures within a nerve trunk, i.e., fasciculi, and the like.
Input device 26 may also include a keyboard or similar input device for supplying commands to processor 14. The command-supplying input device may also be provided as a second input device (not shown) .
The processor 14 receives commands and imaging data from the imaging data source through the input device (s) . The imaging data my be stored in memory 16 for further processing or used by the processor 14 to construct an image of a tissue target. The constructed image is sent to an output device 28 such as a visual image display device, for example, a video display terminal. A video terminal for displaying the original magnetic resonance or magnetic transfer image may also be included. The constructed image may also be displayed numerically or graphically through an appropriate output device.
The imaging device 10 performs the applications programs described herein in conjunction with Figure 2, which axe implemented from compiled source code in the "C" programming language. Referring to Figure 2, the imaging device 10 includes conventional imaging software, such as software for the display of original magnetic images. The imaging device also includes post-scan data processing programs, such as
"clipping" software for discarding irrelevant data, and software for operating interactive, semi-automatic or automatic contouring devices for segmentation of anatomic structures.
The imaging device also includes post-scan data process- ing software for the creation of binary images of segmented structures and binary slices and software for the creation of three-dimensional binary masks and software for the segmentation of gray-level data representing the segmented structure via the binary mask. Post-scan data processing software for the histographic analysis of the gray- level data is also included, which performs single or iterative segmentation of the gray-level anatomic tissue structure into the intrinsic tissue microstructures, as well as software for the definition of material opacity, emission and surface reflection.
The imaging device post-scan data processing software also includes software for the creation of three-dimensional images of individual intrinsic tissue microstructures, or compound images of complete tissue structures, and software for image smoothing (aliasing) . Post-scan data processing software is also included for the calculation of the central trunk (tract) axis (axes) , cross-sections and tissue percentages and the graphical display of values.
Figure 2 depicts an applications program for the imaging of a nerve trunk. In step 28 magnetic imaging data is acqu- ired in 1 mm slices along the nerve trunk axial plane. In step 30 the data is' transferred to a processor for post scan processing. In step 32 gray-level images of the magnetic imaging slices are constructed by the image display software.
In step 34 the interactive clipping software identifies the nerve trunk in the gray-level slices and in step 36 clip . the first nerve trunk slice to eliminate the surrounding tissue imaging data. The rest of the slice stack is then automatically clipped, with modification of the clipping boundary as required. In step 38 the semi-interactive segmentation software identifies the nerve trunk boundary to the computer in the initial one or two slices, after which the software automatically segments the remaining slices and then checks and modifies the identified boundaries.
In step 40 the binary masking software creates two- dimensional binary masks for the segmented nerve boundaries. In step 42 the two-dimensional binary masks are used by the gray-level data segmentation software to segment the gray- level data representing the nerve trunk from the original magnetic resonance or magnetic transfer image slices.
In step 44 the histographic analysis software resegments the nerve trunk gray-level data for the segmented nerve trunk into nerve trunk microstructure data populations, including the fascicular bundles, the epineurium and inter-fascicular fascia. Interpolation of the resegmented data populations is performed in step 46. Software for the definition of material opacity, emission and surface reflection adjusts these features for individual microstructures to enhance their visualization within compound tissue images .
Still and movie images are then constructed by the three-dimensional imaging software in step 48. Both individual nerve tissue microstructure images and compound images of the nerve trunk composed of the individual nerve tissue microstructures are constructed.
In step 50 images are displayed, either as printed i a- ges or by viewing software as video displays. Optionally, the tissue calculation software determines in step 52, from the fasciculi imaging data, the amount of fascicular material in selected cross-sections, the values for which are presented graphically in step 54 by the display software.
Nerve trunks are readily visualized in high resolution TI or T2 weighted magnetic resonance (MR) images. In cross- sections of a nerve trunk the fascicles appear as bright "spots" interspersed in a darker gray matrix which represent interstitial tissues, surrounded by a dark border of variable thickness that represents the perineureum and several layers of thin supporting fascia.
The relatively high amplitude signal from the fasciculi relates to the presence of myelin. Myelin is composed of lipids and protein, as in the plasma membranes of other cells, but the percentage of lipid is higher; and in addition to the usual suite of lipids found in plasma membranes, myelin contains a characteristic glycolipid. Myelin also contains up to 30% bound water.
When demyelination occurs, the loss of lipid molecules and the free interstitial water (which will affect ionic exchange rates between the free and bound water) reduces the amplitude of the MR signal within the bounding perineureum. Specifically, a percentage of the bright "spots" will become no longer visible. Pathologic tissues will appear darker in MR images than do normal tissues, and the presence of fasciculi will be less apparent.
The presence of fasciculi can be partially detected by direct visual inspection of cross-sectional MR images. But determination of the number or percentage of fasciculi present in a cross-section, comparisons between different sections is not satisfactory via visual inspection alone. It is often not possible to resolve darker fasciculi and small changes in signal amplitude may go unobserved.
The method of the present invention is thus used to generate three-dimensional images of the fasciculi, which can be visualized in a static format, or which may be rotated about the (approximate) long axis of the trunk. The latter enhances the visualization process. Because the three-dimensional image can be limited to the fasciculi, the local density of these structures within the trunk, and their brightness (signal amplitude) are now relatively easy to visualize directly.
With the present invention, it is now possible to quantify nerve trunk demyelination values from the generated ima- ges, which were heretofore measured either invasively or subjectively. Nerve tissue demyelination is determined by first calculating the volume of fascicular tissue as a percentage of the volume of the gross trunk (tract) structure. Next, the area of fascicular tissue as a percentage of the area of the gross trunk (tract) structure as seen in selected cross- sections of the trunk (tract) is calculated. The contiguous cross sectional areas of the trunk (tract) of specified thickness (i.e., 1 mm) perpendicular to the long axis is then calculated.
The percentage of volume of fascicular tissue as a percentage of total volume of the trunk (tract) in contiguous cross-sections is calculated next. A graphical display of the percentage volumes obtained vs. location of the sections on the central axis of the trunk is then constructed. The variation of myelin content as indicated by these percentage values is a measure of the neurophysiologic (flow efficiency) status of the trunk (nerve) and thus a diagnostically significant index.
This technique enables a comparison of the percentage of fascicular tissue present, for example, in the forearm (where entrapment and compression is unlikely) , with the value for fascicular tissue within the carpal tunnel, and in regions of the trunk distal to the carpal tunnel . It is also possible to compare fascicular percentages in homologous regions from different individuals (compare between patients) .
The method of the present invention thus allows for the fascicular percentages to be calculated for cross-sectional areas of a nerve trunk. For the median nerve, fascicular percentages for 1 mm thick contiguous sections normal to the long axis of the trunk may be calculated, and then represented as a bivariate plot of percentage of fascicular material present vs. location within the trunk. This graphical representation provides the CTS diagnostician, with a simple and direct means of determining in which precise regions of the median nerve trunk a loss of fascicular material has occurred, and the amount of fascicular tissue lost .
A drawback in using MR images for inter-patient comparisons is that the scale of gray-level values representing different types of tissue is scanner dependant and consequently varies between scanners. Magnetic Transfer Imaging (MTI) is a new technique that has two advantages over conventional MR imaging when applied to analysis of nerve trunks. First, the in-plane resolution is improved. Second, the data values obtained with this technique are not scanner dependant and consequently direct inter-patient comparisons are possible.
The possibility of using MTI for the analysis of nerve tissue has been discussed in the literature. The loss of bound protons in the degenerating myelin sheath has the effect, via magnetic transfer, of depressing the MR signal from free protons in the immediate vicinity of the fasciculi. The Magnetic Transfer function can be calculated, and the resultant values used to create MT images. Thus, MTI is particularly well-suited for use with the methods and apparatus of the present invention.
The method and apparatus of the present invention are also well-suited for the diagnosis of nerve pathology. In diagnostic methods according to the present invention, a clinical examination is typically first performed to determine the likely location of nerve pathology. A high-resolu- tion magnetic image of the region is then obtained.
The MR data comprises thin (1 mm) contiguous slices in a plane which is more or less normal to the axis of the nerve. High resolution images can be obtained by using scan times up to ten minutes.
Post-scan data processing is then applied to the MR data. MR hard-copy (acetate) is not required. The three- dimensional MR data is loaded into the zhree-dimensional imaging software. Steps in the subsequent data processing are as follows:
The original MR slices are displayed. The nerve is identified and image contrast adjusted eo maximize its visualization by the operator. This cor.trast adjustment does not affect the signal differential marking the boundary of the nerve, and is not a factor in subsequent segmentation. The size of the data set is reduced by data clipping to discard all data except that in the immediate vicinity of the nerve. Once a region of interest has been defined, clipping is automatic; thus, the size of the region of interest relates to the position of the nerve on successive images.
Interpolation of the imaging data is then performed. The original voxels are prismatic; i.e., in the Z direction (thickness of image-slice) they are usually 1 mm. The dimen- sions of X and Y (in the plane of the slice) depend on the imaging protocol. Typically X and Y are 0.5-0.7 mm. Cubile voxels are required for 3 -dimensional rendering; thus, it is necessary to interpolate the data set in the Z direction in order to obtain slices corresponding in thickness to the X,Y dimensions. The interpolation process is not simple data averaging; it utilizes a sophisticated algorithm which identifies the "shape" of regions of consistent signal in contiguous slices and calculates intermediate "shapes" (assuming linearity in shape change between slices) . Shape based in- terpolation has been shown to produce highly accurate information in the interpolated slices. Interpolation is usually performed after clipping and segmentation on the binary image produced. This reduces the time required to perform interactive contouring during the process of segmentation.
Automatic segmentation of MR data is difficult because of the non-linear signal. At present several different techniques are used to contour the boundary of the region of interest on the original MR images. They are interactive, live-wire and live-path contouring.
With interactive segmentation the operator traces the boundary of the nerve trunk section, recording coordinates on the boundary in so doing. In live-wire segmentation, the boundary of the nerve as seen in cross-section is identified by the operator and indicated to the computer. The computer is thus able to identify subsequently contiguous voxels with the same boundary signal characteristics. Once the boundary has been identified, the operator can indicate a starting point on it, then move the cursor to a next point some distance away from the first.
The computer is able to track the boundary between the two points using the boundary information provided, and draws a colored line to mark it. Because the boundary signal varies, and the presence of adjacent non-boundary pixels with similar signal characteristics, the boundary selected automatically may vary from that which the operator judges to be correct. In this case the boundary is redefined at the point of divergence (i.e., the boundary-contour is corrected) and the technique of indicating points at intervals along the boundary is continued until the boundary is closed. Boundary identification during contouring is used to refine the suite of signal characteristics used by the computer to determine the boundary automatically. After one or two boundaries have been contoured it is usually only necessary to indicate two or three points to obtain a complete, accurate boundary.
This technique can be extended to automatic identifica- tion of the boundary in subsequent slices, which speeds up the procedure. Usually the operator examines the automatically generated boundary contours and can correct any errors. The continually refined boundary characteristics are stored by the computer and can be utilized to segment other data sets obtained via similar imaging protocols.
Live-path contouring is similar to live-wire contouring, but the computer only seeks boundary characteristics coinciding with those defined for it in a "path" of limited width on either side of the cursor. It is possible to run the cursor rapidly around the region of interest, without following the boundary exactly, to describe a contour within the limits of the defined path, for a very rapid method of contouring.
Live-wire and live-path contouring techniques are normally carried out on all the original MR slices. However, these techniques may be developed so that contours are only defined interactively on the first two or three slices, after which the computer will search for the corresponding boundary in the remaining slices and contour them automatically.
The segmentation techniques output the surface boundary of the region of interest by a binary separation of the three dimensional data based on those voxels within the boundary and those without it. Three-dimensional imaging of the bin- ary data therefore yields a "solid" object. However, the boundary surface can be used as a "mask" to extract MR data within the region of interest from the original three-dimensional MR data matrix. The extracted MR data set is referred to as a gray-level image.
The gray-level image thus contains information about the tissue microstructures contained within the nerve trunk. The gray-level values for nerve fasciculi, matrix, and perineur- eum comprise distinct and overlapping populations. In the region of signal overlap the percentage of signal belonging to each material can be ranges of signal values assigned accordingly. This allows the different tissue types to be resegmented automatically based on this histographic data.
The histographic resegmentation may be performed iteratively in order to resolve any signal overlap to obtain the best visualization of the nerve tissue microstructures (i.e., visualization of the very distinctive fasciculi) .
Opacity, emission and surface strength parameters optimal for visualization of the nerve microstructures are then determined by creation of three-dimensional images and examination of the anatomy. Once values have been determined for these parameters, they may be used for more automated imaging in subsequent data sets. Improvement of visualization by adjusting these parameters is a distinct process from segmentation and does not select the tissue concerned from its surrounding matrix.
The output is visual and numerical. Visual results con- sist of three-dimensional renderings of the trunk, including a compound image of the gross structure, which includes the perineurium, fascicles and interfascicular tissue. The opacity of the perineurium is set so that this tissue is translucent . This enables the content of the nerve trunk to be viewed in the context of the perineurium. Similarly, the opacity of the interfascicular tissue is set to be greater than that of the perineurium, but less than that of the fascicles. The result is to display the microstructures of the nerve in their correct anatomic relationships.
To enhance the three-dimensionality of the image, and to enable the anatomic relationships to be examined from all angles, the image is rotated about its long axis and recorded as multiple "frames" which are stored in a "cine" file. This is performed using conventional MR imaging software. Images are also created containing each of the material sets separately, and cine files created for each. These images allow direct examination of each microstructure separately.
A compound image consisting of the perineureum and fas- ciculae is also created and rendered in a cine file. This image allows the texture and thickness of the perineureum to be examined in the context of the fasciculi, and regions of membrane thickening or thinning compared with regions of MR signal fluctuation in the fasciculi.
The numerical output includes the cross-sectional area value for the nerve trunk normal to its long axis, calculated continuously throughout its length and displayed in a bivari- ate plot of cross-sectional area vs. section location. A similar bivariate plot is calculated to show the density of fascicular material within the cross-section vs. location.
The parameters measured by the method and apparatus of the present invention for nerve tissue include the general morphology of the trunk (tract) , the cross-sectional morphology, the size and distribution of the nerve fasciculi and the density of the fasciculi within cross-sections of the nerve trunk. These parameters are associated with physiologic changes accompanying nerve damage and dysfunction.
Gross morphology and cross-sectional morphology are phy- sical reflections of the pressure environment. It is well known that in regions of pressure a peripheral nerve reorganizes by becoming flattened in a plane normal to the compress- ive force. The amount of flattening almost certainly relates to altered function, but the sensitivity of the relationship may not be adequate for the purpose of early diagnosis.
In regions of compression the size of the fascicles decreases, but they increase in number, forming smaller bundles. The cross-section of the neurites can be measured.
However, cross-sectional fascicular density is a measure of MR signal from the "fascicular material" present in the section. The MR signal may be expected to vary with the his- tochemistry of the nerve fibers (neurites) and their interstitial environment. This includes changes in the structure and chemistry of the myelin sheaths (composed of complex lipids and a high percentage of bound water) and in the ionic environment of the neurite and interstitial tissues. Thus, variation in the MR signal is an indicator of functional changes in the neurites, including those due to trauma.
Signal variation can be detected in two ways. One way is by direct examination of the final image. Regions of nerve trauma are likely to be represented by localized disappearance of the fasciculi. The extent of fascicular disappearance is an indication of the severity of the nerve damage. The MR signal variation can also be quantified by measuring the density of the fascicular material in cross-sections of the nerve trunk and plotting this data against the locations of the cross-sections.
Neural flow studies can be made on the median nerve in the region of the carpal tunnel because of the very superfic- ial placement of the trunk in this region. However, neuro- pathies occur in other regions, both within and outside the central nervous system where neural flow studies are impractical because of the depth of the trunk. The methods and apparatus of the present invention provide an innovative means of diagnosing and assessing neuropathy in such deep peripheral nerves, and in portions of the central nervous system.
The neuropathy need not be associated with compression or other trauma, but may also result from other causes, such as multiple sclerosis and related diseases. Until now it was not possible to assess quantitatively the progression of a disease like multiple sclerosis, or the effectiveness of treatment. Instead, practitioners were forced to rely upon the subjective assessment of the patient. This was particularly inaccurate in most circumstances, in which patients want to believe that a treatment is working and relate their well-being from this perspective.
The methods and apparatus of the present invention may also be used to assess other tissues that are inaccessible because of their depth within the body, including vascular tissues and tumor tissues. The present invention makes it possible to visualize vascular blockage, as well as tumor tissue boundaries. The ability to visualize tumor tissue boundaries makes the present invention a particularly useful tool for the cancer surgeon.
Proof of the concept of method of the invention for the analysis of neuropathy was obtained via a feasibility study of the median nerve in patients diagnosed with CTS and non- symptomatic subjects. The feasibility study is presented to illustrate certain aspects of the invention, and should not be interpreted as limiting the scope of the invention as it is defined by the claims.
Twenty wrists, including four cadaver wrists, were scanned and analyzed using the technique described above. Subjects scanned in vivo included clinically diagnosed patients with severe CTS, less severe CTS, and control (asymptomatic) subjects.
Results showed:
a. Clear feasibility of "dissecting" and quantifying nerve trunk structures is possible using this analytic system. A prior art axial cross-section of a wrist obtained by magnetic resonance imaging is shown in Fig. 3. The cross- sectional image includes the median nerve and associated tendons and carpal bones . This Figure illustrates how such images provided little information of clinical diagnostic value relating to neurophysiology.
A magnetic resonance image of a healthy median nerve trunk constructed by the method and apparatus of the present invention is shown in Fig. 4. The white streaks are bundles of myelinated nerve fibers, i.e., fasciculi. This is in contrast to the magnetic resonance image of a median nerve trunk of a patient afflicted with CTS, shown in Fig. 5, and also constructed by the method and apparatus of the present invention. The nerve fiber demyelination is dramatically evident.
b. MR signal amplitude was reduced where the nerve trunk passes beneath the flexor retinaculum. Reduction of signal amplitude is modest in subjects with no clinical symptoms of CTS, but becomes increasingly significant with CTS symptoms, and is pronounced in subjects diagnosed with severe CTS.
The four cadaver wrists were also subjected to MR scanning. Evidence of demyelination in these wrists was partially confirmed by extirpation of the median nerve from the carpal tunnel and histologic sectioning and staining.
Conclusions .-
The results strongly supported the hypothesis that nerve degeneration can be detected and localized in CTS. The ex- tent of degeneration determined by this procedure in CTS patients appears to be associated with the severity of symptoms.
Thus, using the methods and apparatus of the present invention, it is possible to determine if demyelination has occurred, the percentage of myelin loss in the affected region compared with other "functionally normal" regions, and to precisely localize the neuropathy. This provides a non- invasive, objective and quantified diagnostic procedure.
As will be readily appreciated, numerous variations and combinations of the features set forth above can be utilized without departing from the pre-sent invention as set forth in the claims. Such variations are not regarded as a departure from the spirit and scope of the invention, and all such variations are intended to be included within the scope of the following claims .

Claims

WHAT IS CLAIMED IS:
1. A method for the non-invasive imaging of an anatomic tissue structure in isolation from surrounding tissues, comprising:
receiving from an input device, magnetic imaging data from a patient of said anatomic tissue structure and said surrounding tissues ;
segmenting said imaging data to isolate said anatomic tissue structure imaging data from said imaging data for said surrounding tissues ,-
separating said anatomic tissue structure imaging data into data populations corresponding to tissue microstructures ;
constructing an image from imaging data for at least one of said tissue microstructures; and
storing or displaying said image .
2. The method of claim 1, wherein said magnetic imaging data is acquired by magnetic resonance imaging.
3. The method of claim 1, wherein said magnetic imaging data is acquired by magnetic transfer imaging.
4. The method of claim 1, wherein said anatomic tissue structure imaging data is separated into said microstructure data populations by histographic resegmenting.
5. The method of claim 4, wherein said histographic resegmenting is performed iteratively to separate any overlapping microstructure data populations.
6. The method of claim 1, wherein said image is displayed as a visual image.
7. The method of claim 1, wherein said anatomic tissue structure is segmented by interactive contouring.
8. The method of claim 1, wherein said anatomic tissue structure is segmented by live-wire contouring.
9. The method of claim 1, wherein said anatomic tissue structure is segmented by live-path contouring.
10. The method of claim 1, wherein said constructing step comprises constructing a compound image of said anatomic tissue structure from imaging data for a plurality of said tissue microstructures.
11. The method of claim 10, further comprising the step of adjusting at least one of the material opacity, emission and surface reflection for at least one of said tissue micro- structures to enhance the visualization of said compound image.
12. The method of claim 1, wherein said target tissue is a nerve trunk or tract .
13. The method of claim 12, wherein said microstructures comprise the perineureum, fasciculi and inter-fascicular tissue.
14. The method of claim 13, further including the step of calculating from fasciculi imaging data the presence or loss of myelin in said nerve trunk or tract.
15. The method of claim 14, further including the step of comparing said myelin calculation to a myelin calculation obtained for another portion of said nerve trunk or tract .
16. The method of claim 12, wherein said nerve trunk or tract is the median nerve.
17. The method of claim 1, wherein said anatomic tissue is tumor tissue.
18. The method of claim 1, wherein said anatomic tissue is vascular tissue.
19. An apparatus comprising:
an input device for receiving magnetic imaging data from a patient;
a processing unit comprising:
a first program for segmenting magnetic imaging data to isolate imaging data for a anatomic tissue structure from imaging data for surrounding tissues ;
a second program for separating said anatomic tissue structure imaging data into data populations corresponding to tissue microstructures ;
a third program for constructing an image from imaging data for at least one of said microstructures; and
a processor for sequentially executing said first, second and third programs;
wherein said processor executes:
said first program using said imaging data from said input device ,-
said second program using said anatomic tissue structure imaging data isolated by execution of said first program; and
said third program using at least one of said data populations obtained by execution of said second program; and an output device for displaying said image or a storage device for storing said image .
20. The apparatus of claim 19, wherein said second program separates said anatomic tissue structure imaging data into said microstructure data populations by histographic resegmenting.
21. The apparatus of claim 20, wherein said second program is executed iteratively by said processor to separate any overlapping microstructure data populations.
22. The apparatus of claim 19, wherein said first program segments said imaging data by interactive contouring.
23. The apparatus of claim 19, wherein said first program segments said imaging data by live-wire contouring.
24. The apparatus of claim 19, wherein said first program segments said imaging data by live-path contouring.
25. The apparatus of claim 19, wherein said output device is a visual display device.
26. The apparatus of claim 20, further comprising an output device for displaying the histograph obtained by said histographic resegmenting.
27. The apparatus of claim 19, wherein said input device comprises a device for receiving imaging data of a nerve trunk or tract and surronding tissues,- said second program separates said nerve trunk or tract imaging data into data populations corresponding to fasciculi and other nerve trunk or tract microstructures; said processing unit further comprises a fourth program for calculating from fasciculi imaging data the presence or loss of myelin in a nerve trunk or tract, wherein said processor executes said fourth program after the execution of said second program using imaging data for nerve trunk fasciculi obtained by execution of said second program; and said output device comprises a device for displaying said myelin calculation or a storage device for storing said myelin calculation.
28. The apparatus of claim 27, wherein said processing unit further comprises a fifth program for comparing myelin calculations from different nerve trunks or tracts or different portions of the same nerve trunk or tract, wherein said processor executes said fifth program after the execution of said fourth program using the myelin calculations obtained by execution of said fourth program.
29. The apparatus of claim 26, wherein said output device is a visual display device.
30. A method for imaging demyelination of the median nerve within one of a patient's carpal tunnels, comprising:
receiving from an input device, magnetic imaging data from a patient comprising at least a portion of the median nerve passing through a carpal tunnel,-
segmenting said imaging data to isolate the median nerve imaging data from imaging data for surronding tissues;
separating said median nerve imaging data into data populations corresponding to the fasciculi and other microstructures of said median nerve ;
constructing an image from at least said imaging data for said fasciculi of said median nerve; and
storing or displaying said image.
31. The method of claim 30, wherein said method further includes the step of calculating from said imaging data for said median nerve fasciculi the presence or loss of myelin in said portion of said median nerve passing through said carpal tunnel, and storing or displaying said myelin calculation.
32. The method of claim 30, wherein said median nerve imaging data is separated into microstructure data populations by histographic resegmenting.
33. The method of claim 32, wherein said histographic resegmenting is performed iteratively to separate any overlapping data populations.
34. The method of claim 30, wherein said median nerve imaging data is segmented by interactive contouring.
35. The method of claim 30, wherein said median nerve imaging data is segmented by live-wire contouring.
36. The method of claim 30, wherein said median nerve imaging data is segmented by live-path contouring.
37. The method of claim 30, wherein said median nerve microstructures also include the perineurium and inter- fascicular tissue.
38. The method of claim 30, wherein said imaging data is acquired by magnetic resonance imaging.
39. The method of claim 30, wherein said imaging data is acquired by magnetic transfer imaging.
40. The method of claim 30, wherein said constructing step comprises constructing a compound image of said median nerve comprising imaging data from a plurality of said nerve tissue microstructures.
41. The method of claim 40, wherein at least one of the material opacity, emission and surface reflection for at least one of said microstructures of said compound image are adjusted to enhance the visualization of said median nerve.
42. The method of claim 41, wherein said compound image also includes the perineurium and inter- fascicular tissue, and the opacity of said inter-fascicular tissue is set to be greater than the opacity of said perineurium, but less than that of said fasciculi.
43. The method of claim 42, wherein said opacity of said perineurium is set to appear translucent .
44. The method of claim 30, further comprising the step of comparing said myelin calculation to a myelin calculation obtained for another portion of said median nerve.
45. The method of claim 30, further comprising the step of comparing said myelin calculation to a myelin calculation obtained for the median nerve within said patient's other carpal tunnel .
46. A method for the imaging a nerve trunk or tract suspected of demyelination comprising:
receiving from an input device, magnetic imaging data from a patient of a nerve trunk or tract suspected of demyelination and the surrounding tissues;
segmenting said imaging data to separate said imaging data for said nerve trunk or tract from said imaging data for said surrounding tissues ,-
separating said imaging data for said nerve trunk or ""tract into data populations' corresponding to the fasciculi and other microstructures of said nerve trunk or tract;
constructing an image from at least said imaging data for said fasciculi of said nerve trunk or tract; and
storing or displaying said image.
47. The method of claim 46, further including the step of calculating from said fasciculi imaging data the presence or loss of myelin in said nerve trunk or tract suspected of demyelination, and storing or displaying said myelin calculation.
48. The method of claim 46, wherein said imaging data for said nerve trunk or tract is separated into microstruc- ture data populations by histographic resegmenting.
49. The method of claim 48, wherein said histographic resegmenting is performed iteratively to separate any overlapping microstructure data populations.
50. The method of claim 46, wherein said imaging data for said nerve trunk or tract is segmented by interactive contouring .
51. The method of claim 46, wherein said imaging data for said nerve trunk or tract is segmented by live-wire contouring .
52. The method of claim 46, wherein said imaging data for said nerve trunk or tract is segmented by live-path contouring .
53. The method of claim 46, wherein said microstructures of said nerve trunk or tract also include the perineureum and inter-fascicular tissue.
54. The method of claim 46, wherein said imaging data is "acquired by magnetic resonance imaging.
55. The method of claim 46, wherein said imaging data is acquired by magnetic transfer imaging.
56. The method of claim 46, 'wherein said constructing step comprises constructing a compound image of said nerve trunk or tract comprising imaging data from a plurality of nerve tissue microstructures.
57. The method of claim 56, wherein at least one of the material opacity, emission and surface reflection for at least one of said microstructures of said compound image is adjusted to enhance the visualization of said nerve trunk or tract.
58. The method of claim 57, wherein said compound image also includes the perineurium and inter-fascicular tissue, and the opacity of said inter-fascicular tissue is set to be greater than the opacity of said perineurium, but less than that of said fasciculi.
59. The method of claim 58, wherein said opacity of said perineurium is set to appear translucent .
60. The method of claim 47. , further comprising the step of comparing said myelin calculation to a myelin calculation obtained for another portion of said nerve trunk or tract.
61. The method of claim 46, wherein said nerve trunk or tract is the median nerve .
62. An apparatus comprising:
an input device for receiving magnetic imaging data from a patient for a nerve trunk or tract and the surrounding tissues ;
a processing unit., comprising:.., _..
a first program for segmenting said magnetic imaging data for said nerve trunk or tract from said magnetic imaging data for said surrounding tissues ;
a second program for separating said nerve trunk or nerve trunk imaging data into data populations corresponding to fasciculi and other nerve trunk or tract microstructures ,- a third program for calculating the presence or loss of myelin in said nerve trunk or tract from said fasciculi imaging data; and
a processor for sequentially executing said first, second and third programs,-
wherein said processor executes :
said first program using imaging data from said input device ;
said second program using said nerve trunk or track imaging data isolated by execution of said first program; and
said third program using said fasciculi data population obtained by execution of said second progra ; and
an output device for displaying said myelin calculation or a storage device for storing said myelin calculation.
63. The apparatus of claim 62, wherein said second program separates said nerve trunk or tract imaging data into said microstructure data populations by histographic resegmenting.
64. The apparatus of claim 63, wherein said second program is executed iteratively by said processor to separate any overlapping microstructure data populations.
65. The apparatus of claim 62, wherein said first program segments said imaging data by interactive contouring.
66 . The apparatus of claim 62, wherein said first program segments said data by live-wire contouring.
67. The apparatus of claim 62, wherein said first program segments said data by live-path contouring.
68. The apparatus of claim 62, wherein said output device is a visual display device.
69. The apparatus of claim 63, further comprising an output device for displaying the histograph obtained by said histographic resegmenting.
70. The apparatus of claim 62, wherein said processing unit comprises a fourth program for comparing myelin calcula- tions obtained from different nerve trunks or tracts or different portions of the same nerve trunk or tract, wherein said processor executes said fourth program after the execution of said third program using myelin calculations obtained by execution of said third program.
71. A method for the quantitative measurement of a nerve trunk or tract suspected of demyelination, comprising:
receiving from an input device, magnetic imaging data from a patient of a nerve trunk or tract suspected of demyelination and the surronding tissues;
segmenting said imaging data to separate said imaging data for said nerve trunk or tract from said imaging data for said surrounding tissues;
separating said imaging data for said nerve trunk or tract into data populations corresponding to the fasciculi and other microstructures of said nerve trunk or tract;
calculating from said fasciculi imaging data the presence or loss of myelin in said nerve trunk or tract; and
storing or displaying said myelin calculation.
72. The method of claim 71, wherein said imaging data for said nerve trunk or tract is separated into microstructure data populations by histographic resegmenting.
73. The method of claim 72, wherein said histographic resegmenting is performed iteratively to separate any overlapping microstructure data populations.
74. The method of claim 71, wherein said imaging data for said nerve trunk or tract is segmented by interactive contouring .
75. The method of claim 71, wherein said imaging data for said nerve trunk or tract is obtained by live-wire contouring .
76. The method of claim 71, wherein said imaging data for said nerve trunk or tract is obtained by live-path contouring.
77. The method of claim 71, wherein said imaging data is acquired by magnetic resonance imaging.
78. The method of claim 71, wherein said imaging data is acquired by magnetic transfer imaging.
79. The method of claim 71, further comprising the step of comparing said myelin calculation to a myelin calculation for another portion of said nerve trunk or tract.
80. The method of claim 71, whereir. said nerve trunk or tract is the median nerve.
81. A method for the quantitative r.easurement of median nerve demyelination within one of a patient's carpal tunnels, comprising : receiving from an input device, magnetic imaging data from a patient including at least a portion of the median nerve passing through a carpal tunnel .-
segmenting said imaging data to isolate said median nerve imaging data from any imaging data for surrounding tissues;
separating said median nerve imaging data into data populations corresponding to the fasciculi and other micro- structures of said median nerve ;
calculating from said imaging data for said median nerve fasciculi, the presence or loss of myelin in said median nerve within said carpal tunnel; and
storing or displaying said myelin calculation.
82. The method of claim 81, wherein said median nerve imaging data is separated into microstructure data populations by histographic resegmenting.
83. The method of claim 82, wherein said histographic resegmenting is performed iteratively to separate any overlapping data populations .
84. The method of claim 81, wherein said median nerve imaging data is segmented by interactive contouring.
85. The method of claim 81, wherein said median nerve imaging data is segmented by live-wire contouring.
86. The method of claim 81, wherein said median nerve imaging data is segmented by live-path contouring.
87. The method of claim 81, wherein said imaging data is acquired by magnetic resonance imaging.
88. The method of claim 81, wherein said imaging data is acquired by magnetic transfer imaging.
89. The method of claim 81, further comprising the step of comparing said myelin calculation to a myelin calculation obtained for another portion of said median nerve.
90. The method of claim 81, further comprising the step of comparing said myelin calculation to a myelin calculation obtained for the median nerve within said patient's other carpal tunnel .
PCT/US2001/004479 1998-07-23 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures WO2002064034A1 (en)

Priority Applications (10)

Application Number Priority Date Filing Date Title
US10/467,789 US7672708B2 (en) 1998-07-23 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures
CA002437822A CA2437822A1 (en) 2001-02-12 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures
EP01910573A EP1361820B1 (en) 2001-02-12 2001-02-12 Method for the non-invasive imaging of anatomic tissue structures
PCT/US2001/004479 WO2002064034A1 (en) 2001-02-12 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures
IL15734901A IL157349A0 (en) 2001-02-12 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures
CNA018230741A CN1505490A (en) 2001-02-12 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures
JP2002563835A JP2004528873A (en) 2001-02-12 2001-02-12 Method and apparatus for obtaining non-invasive images of anatomical tissue structures
IL157349A IL157349A (en) 2001-02-12 2003-08-11 Method and apparatus for the non-invasive imaging of anatomic tissue structures
AU2007203153A AU2007203153B2 (en) 2001-02-12 2007-07-06 Method and apparatus for the non-invasive imaging of anatomic tissue structures
US12/716,069 US8195270B2 (en) 1998-07-23 2010-03-02 Method and apparatus for the non-invasive imaging of anatomic tissue structures

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/US2001/004479 WO2002064034A1 (en) 2001-02-12 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US09/121,351 Continuation US6205235B1 (en) 1998-07-23 1998-07-23 Method and apparatus for the non-invasive imaging of anatomic tissue structures

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US10467789 A-371-Of-International 2001-02-12
US12/716,069 Division US8195270B2 (en) 1998-07-23 2010-03-02 Method and apparatus for the non-invasive imaging of anatomic tissue structures

Publications (1)

Publication Number Publication Date
WO2002064034A1 true WO2002064034A1 (en) 2002-08-22

Family

ID=21742334

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2001/004479 WO2002064034A1 (en) 1998-07-23 2001-02-12 Method and apparatus for the non-invasive imaging of anatomic tissue structures

Country Status (8)

Country Link
US (2) US7672708B2 (en)
EP (1) EP1361820B1 (en)
JP (1) JP2004528873A (en)
CN (1) CN1505490A (en)
AU (1) AU2007203153B2 (en)
CA (1) CA2437822A1 (en)
IL (2) IL157349A0 (en)
WO (1) WO2002064034A1 (en)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8221321B2 (en) 2002-06-07 2012-07-17 Verathon Inc. Systems and methods for quantification and classification of fluids in human cavities in ultrasound images
US8221322B2 (en) * 2002-06-07 2012-07-17 Verathon Inc. Systems and methods to improve clarity in ultrasound images
GB2391625A (en) 2002-08-09 2004-02-11 Diagnostic Ultrasound Europ B Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams
US7520857B2 (en) * 2002-06-07 2009-04-21 Verathon Inc. 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume
US20060025689A1 (en) * 2002-06-07 2006-02-02 Vikram Chalana System and method to measure cardiac ejection fraction
US20040127797A1 (en) * 2002-06-07 2004-07-01 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
US20090112089A1 (en) * 2007-10-27 2009-04-30 Bill Barnard System and method for measuring bladder wall thickness and presenting a bladder virtual image
US20080262356A1 (en) * 2002-06-07 2008-10-23 Vikram Chalana Systems and methods for ultrasound imaging using an inertial reference unit
US7819806B2 (en) * 2002-06-07 2010-10-26 Verathon Inc. System and method to identify and measure organ wall boundaries
US20090062644A1 (en) * 2002-06-07 2009-03-05 Mcmorrow Gerald System and method for ultrasound harmonic imaging
US20100036252A1 (en) * 2002-06-07 2010-02-11 Vikram Chalana Ultrasound system and method for measuring bladder wall thickness and mass
CN101313333B (en) * 2005-11-21 2013-04-24 皇家飞利浦电子股份有限公司 Method for creating a model of a structure
US8133181B2 (en) * 2007-05-16 2012-03-13 Verathon Inc. Device, system and method to measure abdominal aortic aneurysm diameter
US8167803B2 (en) * 2007-05-16 2012-05-01 Verathon Inc. System and method for bladder detection using harmonic imaging
US8160344B2 (en) * 2008-04-22 2012-04-17 Siemens Medical Solutions Usa, Inc. Iterative segmentation of images for computer-aided detection
US8225998B2 (en) * 2008-07-11 2012-07-24 Es&S Innovations Llc Secure ballot box
US8588486B2 (en) * 2009-06-18 2013-11-19 General Electric Company Apparatus and method for isolating a region in an image
JP5606832B2 (en) * 2010-03-05 2014-10-15 富士フイルム株式会社 Image diagnosis support apparatus, method, and program
US10492717B2 (en) * 2013-11-04 2019-12-03 The Research Foundation Of State University Of New York Methods, systems, and devices for determining and visually indicating demyelinated pathways
WO2017090044A1 (en) * 2015-11-26 2017-06-01 Tpm Medical Systems Ltd. Systems and methods for mri based imaging of peripheral nerves
CN108022291B (en) * 2017-12-01 2020-08-28 中山大学附属第一医院 Visual three-dimensional reconstruction method for human body peripheral nerve bundle type structure

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5560360A (en) 1992-03-09 1996-10-01 University Of Washington Image neurography and diffusion anisotropy imaging
US5706813A (en) * 1992-03-09 1998-01-13 University Of Washington Focal neurographic magnetic resonance imaging system
US6205235B1 (en) * 1998-07-23 2001-03-20 David Roberts Method and apparatus for the non-invasive imaging of anatomic tissue structures

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5185809A (en) * 1987-08-14 1993-02-09 The General Hospital Corporation Morphometric analysis of anatomical tomographic data
US4945478A (en) * 1987-11-06 1990-07-31 Center For Innovative Technology Noninvasive medical imaging system and method for the identification and 3-D display of atherosclerosis and the like
JP3457310B2 (en) * 1992-05-05 2003-10-14 ユニバーシティ オブ ワシントン Image nerve recording and diffusion anisotropic image processing
JPH07371A (en) * 1992-10-14 1995-01-06 Ge Yokogawa Medical Syst Ltd Neural path imaging method and mri system therefor
US5757953A (en) * 1996-02-29 1998-05-26 Eastman Kodak Company Automated method and system for region decomposition in digital radiographic images
US7048716B1 (en) * 1997-05-15 2006-05-23 Stanford University MR-compatible devices
AU4424499A (en) * 1998-06-08 1999-12-30 Brown University Research Foundation Method and apparatus for automatic shape characterization
US6526305B1 (en) * 1998-11-25 2003-02-25 The Johns Hopkins University Method of fiber reconstruction employing data acquired by magnetic resonance imaging
US7187810B2 (en) * 1999-12-15 2007-03-06 Medispectra, Inc. Methods and systems for correcting image misalignment
US6463315B1 (en) * 2000-01-26 2002-10-08 The Board Of Trustees Of The Leland Stanford Junior University Analysis of cerebral white matter for prognosis and diagnosis of neurological disorders
US7095890B2 (en) * 2002-02-01 2006-08-22 Siemens Corporate Research, Inc. Integration of visual information, anatomic constraints and prior shape knowledge for medical segmentations
JP4248822B2 (en) * 2002-08-29 2009-04-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Fiber rendering method and fiber rendering device
EP1644869A4 (en) * 2003-07-15 2011-01-26 Medical Metrx Solutions Inc Using statistical process control (spc) to demonstrate solution convergence in a technician guided segmentation system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5560360A (en) 1992-03-09 1996-10-01 University Of Washington Image neurography and diffusion anisotropy imaging
US5706813A (en) * 1992-03-09 1998-01-13 University Of Washington Focal neurographic magnetic resonance imaging system
US6205235B1 (en) * 1998-07-23 2001-03-20 David Roberts Method and apparatus for the non-invasive imaging of anatomic tissue structures

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
AMARTUR S C ET AL: "TISSUE SEGMENTATION FOR THREE-DIMENSIONAL DISPLAY OF HUMAN SPINES", MEDICAL PHYSICS, AMERICAN INSTITUTE OF PHYSICS. NEW YORK, US, vol. 18, no. 2, 1 March 1991 (1991-03-01), pages 305 - 308, XP000228818, ISSN: 0094-2405 *
HOWE F A: "MAGNETIC RESONANCE NEUROGRAPHY", MAGNETIC RESONANCE IN MEDICINE, ACADEMIC PRESS, DULUTH, MN, US, vol. 28, no. 2, 1 December 1992 (1992-12-01), pages 328 - 338, XP000330806, ISSN: 0740-3194 *
KOTTKE D ET AL: "Adaptive segmentation of coronary angiograms", PROCEEDINGS OF THE FOURTEENTH ANNUAL NORTHEAST BIOENGINEERING CONFERENCE (IEEE CAT. NO.88-CH2666-6), DURHAM, NH, USA, 10-11 MARCH 1988, 1988, New York, NY, USA, IEEE, USA, pages 287 - 290, XP000092904 *
ORANGE C M ET AL: "Magnetic contour tracing", PROCEEDING IEEE WORKSHOP ON VISUALIZATION AND MACHINE VISION (CAT. NO.94TH0636-1), PROCEEDINGS OF WORKSHOP ON VISUALIZATION AND MACHINE VISION, SEATTLE, WA, USA, 24 JUNE 1994, 1994, Los Alamitos, CA, USA, IEEE Comput. Soc. Press, USA, pages 33 - 44, XP002179287, ISBN: 0-8186-5875-4 *
YU XIAOHAN ET AL: "Interactive surface segmentation for medical images", ICSP '96. 1996 3RD INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING PROCEEDINGS (CAT. NO.96TH8116), PROCEEDINGS OF THIRD INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING (ICSP'96), BEIJING, CHINA, 14-18 OCT. 1996, 1996, New York, NY, USA, IEEE, USA, pages 1178 - 1181 vol.2, XP002179286, ISBN: 0-7803-2912-0 *

Also Published As

Publication number Publication date
US20040076317A1 (en) 2004-04-22
AU2007203153B2 (en) 2008-12-18
CA2437822A1 (en) 2002-08-22
US7672708B2 (en) 2010-03-02
US20100222665A1 (en) 2010-09-02
AU2007203153A1 (en) 2007-07-26
JP2004528873A (en) 2004-09-24
EP1361820B1 (en) 2012-12-26
US8195270B2 (en) 2012-06-05
EP1361820A1 (en) 2003-11-19
CN1505490A (en) 2004-06-16
IL157349A (en) 2006-10-05
IL157349A0 (en) 2004-02-19

Similar Documents

Publication Publication Date Title
US8195270B2 (en) Method and apparatus for the non-invasive imaging of anatomic tissue structures
US6205235B1 (en) Method and apparatus for the non-invasive imaging of anatomic tissue structures
EP1315125B1 (en) Image processing method and system for disease detection
AU2002251559B9 (en) Three-dimensional joint structure measuring method
US8139836B2 (en) Automatic segmentation of the heart and aorta in medical 3-D scans without contrast media injections
US6754376B1 (en) Method for automatic segmentation of medical images
JP2004032684A (en) Automated method and apparatus for detecting mass or substantial tissue deformation in medical image using computer
EP3195255B1 (en) Method for displaying easy-to-understand medical images
CN1502310A (en) Method and system for measuring disease relevant tissue changes
KR20080027861A (en) Image processing
Liao et al. Tracing the deformed midline on brain CT
Michopoulou Image analysis for the diagnosis of MR images of the lumbar spine
JP2016504070A (en) Quantifying tissue surface roughness based on image data and determining the presence of disease based on it
AU2001238167A1 (en) Methos and apparatus for the non-invasive imaging of anatomic tissue structures
Karwoski et al. Processing of CT images for analysis of diffuse lung disease in the lung tissue research consortium
Umakoshi et al. Quantitative follow-up assessment of patients with interstitial lung disease by 3D-curved high-resolution CT imaging parallel to the chest wall
EP3712844B1 (en) Method and system for fast assessment of brain change normality
JP2005046527A (en) Method for quantitative diagnosis of liver tissue using ultrasonic liver image
Koompairojn et al. Semi-automatic segmentation and volume determination of brain mass-like lesion
Kao et al. Fully automatic abdominal fat segmentation system from a low resolution CT image
Sato et al. Development of comparative reading system using 1D SOM for brain dock examinations

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CR CU CZ DE DK DM DZ EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT TZ UA UG US UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2437822

Country of ref document: CA

WWE Wipo information: entry into national phase

Ref document number: 157349

Country of ref document: IL

Ref document number: 2002563835

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2001238167

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 2001910573

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 018230741

Country of ref document: CN

WWP Wipo information: published in national office

Ref document number: 2001910573

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 10467789

Country of ref document: US

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642