CA2854844C - Correlated diffusion imaging system and method for identification of biological tissue of interest - Google Patents

Correlated diffusion imaging system and method for identification of biological tissue of interest Download PDF

Info

Publication number
CA2854844C
CA2854844C CA2854844A CA2854844A CA2854844C CA 2854844 C CA2854844 C CA 2854844C CA 2854844 A CA2854844 A CA 2854844A CA 2854844 A CA2854844 A CA 2854844A CA 2854844 C CA2854844 C CA 2854844C
Authority
CA
Canada
Prior art keywords
tissue
diffusion
interest
signal
imaging
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CA2854844A
Other languages
French (fr)
Other versions
CA2854844A1 (en
Inventor
Alexander Sheung Lai Wong
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CA2854844A priority Critical patent/CA2854844C/en
Publication of CA2854844A1 publication Critical patent/CA2854844A1/en
Application granted granted Critical
Publication of CA2854844C publication Critical patent/CA2854844C/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4375Detecting, measuring or recording for evaluating the reproductive systems for evaluating the male reproductive system
    • A61B5/4381Prostate evaluation or disorder diagnosis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • 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
    • G01R33/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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
    • G01R33/563Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341Diffusion imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475User input or interface means, e.g. keyboard, pointing device, joystick
    • A61B5/748Selection of a region of interest, e.g. using a graphics tablet
    • A61B5/7485Automatic selection of region of interest
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • 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
    • G01R33/5602Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Biophysics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Signal Processing (AREA)
  • Public Health (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Reproductive Health (AREA)
  • Gynecology & Obstetrics (AREA)
  • Vascular Medicine (AREA)

Abstract

There is disclosed a novel form of imaging referred to in this disclosure as "correlated diffusion imaging" or CDI in which the tissue being imaged is characterized by a joint correlation of diffusion signal attenuation across multiple gradient pulse strengths and timings. Advantageously, by taking into account signal attenuation at different water diffusion motion sensitivities, correlated diffusion imaging can provide significantly improved delineation between cancerous tissue and healthy tissue when compared to existing diffusion imaging modalities. In an embodiment, the method comprises performing quantitative evaluation using receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment to study correlated diffusion imaging for the task of identification of biological tissue of interest. In another embodiment, the method comprises comparing T2- weighted imaging results with that obtained using standard diffusion imaging (via the apparent diffusion coefficient (ADC)) and with that obtained using CDI for tissue characterization and analysis. In still another embodiment, the method comprises of a dual-stage signal mixing configuration of CDI that provides better visualization of anatomical information while preserving strong delineation between cancerous tissue and healthy tissue.

Description

CORRELATED DIFFUSION IMAGING SYSTEM AND METHOD
FOR IDENTIFICATION OF BIOLOGICAL TISSUE OF INTEREST
FIELD OF THE INVENTION
The present invention relates generally to the field of medical imaging, and more particularly to a system and method for identification of biological tissue of interest utilizing improved imaging technologies and image processing techniques.
BACKGROUND
Various imaging technologies may be used to create detailed images of internal organs and internal physical features. For example, ultrasound utilizes high frequency sound waves that are reflected by internal tissue at different intensities to produce images viewable in real-time on a computer display. Ultrasound is commonly used for safely imaging a developing foetus in a pregnant woman, for example. Thermography utilizes differences in regional temperatures to create an infrared digital image of tissue, and is sometimes used in screening breast tissue for identification of potentially cancerous breast tumours. Tomography utilizes x-rays to provide imaging of a single plane or slice of tissue. One type of tomography is a Computed Axial Tomography scan (more commonly known as a CAT scan) which utilizes helical tomography to produce a series of thin slice or section images of a body to generate a virtual three-dimensional image of the inside of a body. While effective for screening certain types of tissues inside the body, existing imaging technologies and image processing techniques are still limited when it comes to screening for certain types of tissues located deep within the body which are substantially similar to the surrounding tissue.
Detection and localization of prostate cancer is one such challenge. Prostate cancer is the most common form of cancer diagnosed in men, with roughly 241,740 new cases in 2012 in the United States alone. Furthermore, prostate cancer is the second leading cause of cancer death in males in the United States, with an estimated 28,170 deaths in 2012. Similar per capita rates of deaths from prostate cancer can be found in other jurisdictions. Given that the median patient survival time for metastatic prostate cancer ranges from about 12.2 to 21.7 months, early clinical , , diagnosis of prostate cancer is key to improving the treatment and longevity of patients affected by prostate cancer.
Presently, conventional clinical diagnosis of prostate cancer involves a prostate specific antigen (PSA) screening, where high PSA levels are considered indicative of possible signs of prostate cancer. However, PSA screening has resulted in significant over-diagnosis of men suspected of having prostate cancer but who do not actually require treatment. As a consequence, many men are over-treated with therapies that carry significant risks in themselves.
Furthermore, there is still no reliable, widely accepted method of diagnostic imaging for prostate cancer. Although transrectal ultrasound (TRUS) is used routinely as a guide for biopsy, it cannot be used to visualize cancer foci because many tumours in the prostate gland are isoechoic and cannot be easily differentiated from surrounding tissue, resulting in sensitivity and specificity in the range of only around 40-50%. Positron emission tomography (PET) have also been investigated as a potential imaging modality for prostate cancer detection, with a number of different tracers that have shown promise for identifying prostate cancer. However, the spatial resolution achieved using PET may not be adequate to properly localize and detect early stage prostate cancer. T2-weighted magnetic resonance imaging (MRI) has also been investigated for prostate cancer detection, but currently requires highly-qualified subspecialty radiologists to interpret the data due to its weak delineation between cancerous tissue and healthy tissue.
Furthermore, in the peripheral zone of the prostate gland, the low T2 signal intensity that is associated with prostate cancer may also be due to a number of noncancerous abnormal conditions such as inflammation and hemorrhaging.
If detected at an early stage, the prognosis for recovery from prostate cancer is excellent. Hence, early detection and the localization of prostate cancer is crucial for diagnosis, as well as treatment via targeted focal therapy. Therefore, what is needed are further improvements in identification of tissues of interest utilizing novel imaging technologies and image processing techniques.
2 SUMMARY
The present invention relates to a system and method for identification of biological tissue of interest utilizing improved imaging technologies and image processing techniques. In an aspect, a new form of imaging is introduced and referred to in this disclosure as "correlated diffusion imaging" or CDI. With CDI, the tissue being imaged is characterized by a joint correlation of diffusion signal attenuation across multiple gradient pulse strengths and timings.
Advantageously, by taking into account signal attenuation at different water diffusion motion sensitivities, correlated diffusion imaging can provide significantly improved delineation between cancerous tissue and healthy tissue when compared to existing diffusion imaging modalities.
In an embodiment, the method comprises performing quantitative evaluation using receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment to study correlated diffusion imaging for the task of identification of biological tissue of interest.
In another embodiment, the method comprises comparing T2-weighted imaging results with that obtained using standard diffusion imaging (via the apparent diffusion coefficient (ADC)) and with that obtained using CDI for tissue characterization and analysis.
Experimental results suggest that CDI can provide significantly improved delineation between healthy tissue and cancerous tissue based on signal attenuation of the tissue at different water diffusion motion sensitivities. This has significant potential as a diagnostic tool for cancer detection and localization, such as cancerous tissue in the prostate gland.
In still another embodiment, the method comprises of a dual-stage signal mixing configuration of CDI that provides better visualization of anatomical information while preserving strong delineation between cancerous tissue and healthy tissue.
In still another aspect, CDI may be implemented as a diagnostic medical device which may be utilized as an aid in the detection and localization of cancer.
In this respect, before explaining at least one embodiment of the invention in detail, it is to be understood that the invention is not limited in its application to the details of construction and to
3 the arrangements of the components set forth in the following description or the examples provided therein, or illustrated in the drawings. The invention is capable of other embodiments and of being practiced and carried out in various ways. Also, it is to be understood that the phraseology and terminology employed herein are for the purpose of description and should not be regarded as limiting.
BRIEF DESCRIPTION OF THE DRAWINGS
The present invention will be better understood and objects of the invention will become apparent when consideration is given to the following detailed description thereof. Such description makes reference to the annexed drawings wherein:
FIG. 1 shows an illustrative method of performing correlated diffusion imaging (CDI) in accordance with an embodiment;
FIG. 2 shows an illustrative echo-planar sequence with two gradient pulses in accordance with an embodiment;
FIG. 3 shows illustrative receiver operator characteristic (ROC) curves in accordance with an embodiment;
FIGS. 4(a) and 4(b) show illustrative results obtained from prior art digital imaging techniques, in comparison to FIG. 4(c) which shows an illustrative digital image of potentially cancerous tissue acquired through CDI in accordance with embodiment;
FIGS. 5(a) and 5(b) show another illustrative example of prior art digital imaging techniques, in comparison to FIG. 5(c) which shows an illustrative digital image acquired through CDI in accordance with an embodiment;
FIGS. 6(a) and 6(b) show another illustrative example of prior art digital imaging techniques, in comparison to FIGS. 6(c) which shows illustrative digital images in accordance with an embodiment. FIG. 6(d) shows an illustrative example of a digital image acquired through CDI
overlaid on a T2-weighted image, in which the potentially cancerous tissue is highlighted in a different color;
4 FIG. 7 is an illustrative example of a method of performing correlated diffusion imaging (CDI) in accordance with another embodiment.
FIGS. 8(a) to 8(c) show another illustrative examples of digital images acquired through the dual-stage signal mixing configuration of CDI in accordance with an embodiment in which potentially cancerous tissue is highlighted with arrows.
FIG. 9 shows an illustrative example of a generic computing device that may provide a suitable operating environment for one or more embodiments.
In the drawings, embodiments of the invention are illustrated by way of example. It is to be expressly understood that the description and drawings are only for the purpose of illustration and as an aid to understanding, and are not intended as a definition of the limits of the invention.
DETAILED DESCRIPTION
As noted above, the present invention relates to a system and method for identification of biological tissue of interest utilizing improved imaging technologies and image processing techniques, termed herein as correlated diffusion imaging.
Diffusion imaging is a promising imaging modality for diagnosing cancer has been developed in which pairs of opposing magnetic field gradient pulses are applied to obtain sensitivity to the Brownian motion of water molecules in tissues. The differences in diffusion characteristics between tissue types facilitate for tissue characterization. As such, given the presumed high cellular density of prostate cancer, the associated tissues should exhibit restricted diffusion characteristics (and as such should have lower apparent diffusion coefficient (ADC) values).
While diffusion imaging shows considerable promise, particularly when used in multi-parametric imaging scenarios, delineating between cancerous tissue and healthy tissue in the prostate gland remains a challenge, due partly to the necessity for fine-tuning the strength, duration, and timing of the applied diffusion gradient pulses. Hence, the characteristics between cancerous tissue and healthy tissue may appear to have substantial overlap depending on the way the gradient pulses are applied, thus making it difficult to detect and localize cancer. As such, an alternative form of magnetic resonance imaging that gets around this issue is highly desired.
5 As will now be described in detail, the presented system and method introduces a new form of diffusion magnetic resonance imaging referred to in the present disclosure as correlated diffusion imaging (CDI), which takes advantage of the joint correlation in signal attenuation across multiple gradient pulse strengths and timings to not only reduce the dependency on the way diffusion gradient pulses are applied, but also improve delineation between cancerous and healthy tissue. To the best of the inventor's knowledge, there are no previous imaging techniques that take this novel approach to cancer assessment. Illustrative embodiments of the system and method will now be described. First, the materials and methods underlying CDI
are described in the Method Section below.
Method In an illustrative example, the methodology behind correlated diffusion imaging (CDI) is summarized in FIG. 1. First, multiple signal acquisitions are conducted from a subject at different gradient pulse strengths and timings. Second, the acquired signals (and any interpolated and/or extrapolated signals) are then mixed together to produce a final signal that characterizes the tissue being imaged. A detailed description of the steps involved in the methodology is presented below.
Imaging Protocol In one embodiment of CDI as a specific example for illustrating and evaluating the effectiveness of CDI for prostate cancer diagnosis, twenty patient cases with known prostate cancer were used as part of the study. The patients ranged in age from 58-80 years, with a median age of 69 years.
Informed consent was obtained from all patients, and approval for the study was obtained from the ethics review board of Sunnybrook Health Sciences Centre, located in Toronto, Canada. All results were reviewed by an expert radiologist with 16 years of experience interpreting body MRI and 11 years of experience interpreting prostate MRI.
Examinations using CDI were performed using a Philips Achieva 3.0T machine at the Sunnybrook Health Sciences Centre. In this specific example embodiment, the axial echo-planar sequence was performed for CDI with the following imaging parameters: TR range from 3336-6178 ms with a median of 4890 ms, and TE ranged from 61-67 ms with a median of 61 ms. The
6 resolution of the signal acquisitions ranged from 1.36 x 1.36 mm2 to 1.67 x 1.67 mm2 with a median of 1.56 x 1.56 mm2. Slice thickness ranged from 3.0-4.0 mm with a median of 3.5 mm.
The display field of view (DFOV) ranged from 20 x 20 cm2 to 24 x 24 cm2 with a median of 24 x 24 cm2.The invention is capable of other embodiments and of being practiced with other MRI
machines.
For comparison purposes, apparent diffusion coefficient (ADC) maps were also obtained using the same axial echo-planar sequence with the same imaging parameters and E2 =
{0s/mm2, 100s/
s/mm2, 1000s/mm2}, as it is considered state-of-the-art for prostate cancer analysis in existing diffusion imaging. Finally, axial T2-weighted imaging acquisitions with the same slice locations as the CDI sequence were obtained as a baseline reference of comparison.
Examinations using T2-weighted imaging were performed using a Philips Achieva 3.0T machine with the following imaging parameters: TR range from 4688-7504 ms with a median of 6481 ms, and TE range from 110-120 ms with a median of 120 ms. Slice thickness ranged from 3.0-4.0 mm with a median of 3.5 mm. The display field of view (DFOV) ranged from 20 x 20 cm2 to 24 x 24 cm2 with a median of 24 x 24 cm2.
Signal Acquisition In an embodiment, as a first step of the CDI imaging process, axial single-shot echo-planar sequences with two gradient pulses of equal magnitude (one pulse in each side of the 180 pulse to dephase and rephase the spins, respectively), as shown in FIG. 2 are used to obtain multiple signal acquisitions using a set of different configurations of gradient pulse strengths and timings, which we will denote as E2 ={ q I i=1,...,/V}, where q, denotes the th sequence. It will be appreciated, however, that other sequences may be used in alternative embodiments, including for example multi-shot echo-planar sequences and interleaved echo-planar sequences.
Imperfect rephasing occurs due to motion of water molecules, leading to attenuation in the acquired signal and thus allowing for the study of water diffusion based on signal attenuation behavior. By varying the configuration of gradient pulse strengths and timings between signal acquisitions, each signal acquisition is sensitive to a different degree of Brownian motion of water molecules in tissues, thus providing unique information with respect to the water diffusion
7 characteristics of the tissue being imaged. The different configurations of gradient pulse strengths and timings can be defined by the following set of parameters:
qi = GE, Si, IN,), (1) where, for the ith sequence, G, denotes the gradient pulse strength, 5, denotes the gradient pulse duration, and A, denotes time between gradient pulses. By grouping the gradient terms, the configuration of gradient pulse strengths and timings used for a particular sequence q, can be simplified to = /72282 A
Y ki-xt (2) where y denotes the proton gyromagnetic ratio.
Signal Mixing As the second step of the CDI imaging process, the multiple signal acquisitions are mixed together to obtain the final signal that characterizes the tissue being imaged. Here, of interest is not the signal attenuation obtained using the individual configurations of gradient pulse strengths and timings, but rather in the local correlation of signal attenuation across the different configurations of gradient pulse strengths and timings within a local spatial sub-volume V to provide a better overall characterization of the water diffusion properties of the tissue being imaged. As such, it is desirable to mix all of the signal acquisitions together into a single quantitative signal characterizing the local signal attenuation correlation.
It will be appreciated that the implementation of the present system and method is not limited to mixing signal acquisitions together, as interpolated and/or extrapolated signals based on signal acquisitions may also be mixed with the signal acquisitions together in other embodiments.
To achieve this goal, the following signal mixing function C(x) is introduced for characterizing local signal attenuation correlation, which is parameterized by diffusion range defined by [qa, qp]
and is defined as
8 C17õ .9,j (A) Sõar" (Al f (S,õ(x). . .9õ(x), V (4) (14, (s) dS,10 (4), (3) where x denotes spatial location, S denotes the acquired signal, pa,...,pfl denotes scaling factors, f denotes a function that relates Sqa (x),..., &gig f_xl within the sub-volume around x (denoted by V(x)). In this embodiment, the function f is set as the conditional joint probability (e.g. f (Sqa (x),..., Sqg (x) I V (x))). It will be appreciated, however, that the implementation of the present system and method is not limited to this particular function f, and other functions may be used for f in other embodiments.
For the embodiment of CDI used in this study, [q q [3] was set at [0s/mm2, 2000s/mm2], and V
was defined as a 7 mm3 spatial sub-volume for assessment purposes as it was found to provide good tissue delineation.
Image Analysis and Interpretation In experimentation conducted by the inventor, the ADC maps and CDI images were reconstructed using the ProCanVAS (Prostate Cancer Visualization and Analysis System) platform developed at the University of Waterloo Vision and Image Processing research group, and were analyzed such that each modality was analyzed independent of other modalities. All visual assessments were made by an expert radiologist with 16 years of experience interpreting body MRI and 11 years of experience interpreting prostate MRI. The invention is capable of other embodiments and of being practiced with other image visualization platforms.
Statistical Analysis Two different analysis strategies were performed to quantify the potential of CDI as a tool for prostate cancer detection and localization. In the first analysis strategy, a receiver operating characteristic (ROC) curve analysis was performed using CDI to quantitatively assess prostate detection and localization. The ROC curves were estimated assuming bivariate normal data. For illustrative purposes, the ROC curves obtained from the pooled data of all patient cases was plotted. To provide a quantitative assessment of diagnostic accuracy, the area under the ROC
curve (Az) was obtained as a single metric of diagnostic accuracy. For comparison purposes,
9 ROC curve analysis was also performed using ADC map as the baseline reference method for assessing prostate cancer using diffusion imaging.
In the second analysis strategy, the inventor studied whether CDI would be a useful imaging modality for building computer-aided clinical decision support systems to assist in the prostate cancer detection and localization process. To quantify the usefulness of CDI
for the purpose of building such systems, leave-one-out cross-validation (LOOCV) trials were performed across all patient cases. For each trial, a two-class Maximum Likelihood (ML) classifier model is trained based on the CDI signal intensity statistics of the individual voxels within the prostate gland (one class characterizing cancerous tissue, with the other class characterizing healthy tissue) across the training patient cases. This learned two-class ML classifier model is then used to calculate sensitivity, specificity, and accuracy for the validation patient case. This process is repeated for a number of trials so that each patient case is used once as the validation patient case. The same was performed on ADC for comparative purposes.
Results To visualize the diagnostic performance for all patient cases, ROC curves for CDI and ADC map results from all patient cases are shown in FIG. 3. It can be observed that improved ROC
characteristics are exhibited by CDI when compared to an ADC map. Furthermore, the area under the ROC curve for CDI is higher with Az = 0.9789, compared to the ROC
curve for ADC
map with Az = 0.9183. The overall sensitivity, specificity, and accuracy results from the LOOCV trials are shown in Table I, below.
Table 1 Leave-one-out cross-validation (LOOCV) results Sensitivity Specificity Accuracy CL 4 0.8676 0.9 :4 A pc. 0.8236 0.7679 17691 It can be observed that the sensitivity, specificity, and accuracy are higher for CDI when compared to ADC, which indicates the potential usefulness of CDI as an imaging modality for building computer-aided clinical decision support systems to assist in the prostate cancer detection and localization process.

FIGS. 4(a)-4(c), 5(a)-5(c), and 6(a)-6(c) show illustrative examples of slices from T2-weighted imaging, ADC map, and CDI of three patient cases out of the twenty patient cases used in the ROC analysis, and a number of observations can be made. Note that example slices show cancerous regions within the prostate gland, not benign prostatic hyperplasia (BPH) nodules.
There is weak visual delineation between prostate cancer and healthy tissue in the prostate gland in the T2-weighted imaging, thus making it difficult even for highly-qualified subspecialty radiologists to interpret (particularly in FIGS. 4(a) and 4(b) and 5(a) and 5(b) where there is no decrease in signal in the cancerous region). The ADC map provides improved visual delineation compared to the T2-weighted imaging; however, it can be observed that there are some cases (e.g., FIGS. 6(a) and 6(b)) where the boundary delineation between tumor and healthy tissue is still difficult to assess. In contrast, the results obtained from using correlated diffusion imaging or CDI in accordance with the present invention is shown in FIG. 6(c). As can be seen, FIG. 6(c) provides a much clearer indication of the locations and boundaries of the prostate cancer compared to the ADC maps for all patient cases. Hence, these experimental results are very promising and motivating for the potential of CDI as a diagnostic tool for prostate cancer detection and localization.
Discussion In this disclosure, a new form of diffusion magnetic resonance imaging has been introduced called correlated diffusion imaging or CDI, which quantifies joint correlation in signal attenuation across multiple diffusion gradient pulse strengths and timings.
The experimental results in the study evaluating the performance of CDI show that CDI can be an effective tool for prostate cancer detection and localization.
However, it is also important to understand the merits of CDI in relation to practical aspects of clinical image acquisition, post-processing, and analysis. One of the attractive characteristics of CDI from a clinical image acquisition perspective is that the signal acquisition process of CDI
can be performed on existing clinical imaging systems without hardware modifications, although the signal mixing process of CDI can also be built in hardware for integration into clinical imaging systems. The signal mixing and reconstruction process of CDI can all be performed post-acquisition on a computer workstation using additional computer-aided clinical decision support software such as the ProCanVAS (Prostate Cancer Visualization and Analysis System) platform developed at the University of Waterloo Vision and Image Processing research group.
Once reconstructed, the CDI images can be viewed on any existing DICOM viewer software, making it easy to integrate into existing radiology workflows.
One possible explanation for CDI's potential to be a more effective tool for prostate cancer detection and localization when compared to standard clinical practice ADC
maps may be related to the highly restrictive water diffusion nature of prostate cancer.
While different gradient pulse strengths and timings may be more sensitive to different degrees of water diffusion motion, this highly restrictive diffusion nature results in signal attenuation that is similar or lower than healthy tissue at all degrees of sensitivity. Therefore, this results in consistently low signal attenuation of prostate cancer compared with healthy tissue irrespective of gradient pulse strengths and timings which, in combination with the possible higher water content of cancerous tissue compared to healthy tissue, may lead to improved cancer and healthy tissue delineation in CDI.
Dual-Stage Signal Mixing FIG. 7 is an illustrative example of a method of performing correlated diffusion imaging (CDI) utilizing a dual-stage signal mixing methodology in accordance with another embodiment. In this embodiment, a dual-stage signal mixing configuration involves an additional third step of the CDI imaging process. After a step of the CDI imaging process in which multiple signal acquisitions, and any interpolated and/or extrapolated signals based on signal acquisitions, are mixed together using the first signal mixing function C(.20, the output quantitative signal CI of this first signal mixing process is then mixed with a baseline signal Sb (obtained in one embodiment using baseline sequence qb such as setting qb = Os/mm2) via a second signal mixing process defined by the signal mixing function C2():
C 2(4 == = TbS b ) (4) exp (Tel (x.))f (Sb(2),Ci(g....) V (1)) dSb(Z)dCi(), f cap where x denotes spatial location, Sb denotes the acquired baseline signal, tb denotes the scale factor of Sb, id i denotes the scale factor of C1, f denotes a function that relates Sb (x), C1 (xj within the local sub-volume around x (denoted by V (1). In this embodiment, the function f is set as the conditional joint probability (e.g., f (S b (x), C (x) V (x))). It will be appreciated, however, that the implementation of the present system and method is not limited to this particular function f, and other functions may be used for in other embodiments.
While the dual-stage mixing example as shown in FIG. 7 and described above uses an MRI
signal as the baseline signal Sb, it will be appreciated that other types of image signals may be used, including for example computed tomography (CT) and ultrasound.
In addition, during the mixing stage, the signals that are mixed may be a combination of acquired signals, and interpolated or extrapolated signals based on acquired signals.
A key benefit of the dual-stage signal mixing configuration of CDI is that anatomical information is better visualized in the resulting output signal C2 while preserving strong delineation between cancerous tissue and healthy tissue. In comparison to the single-stage signal mixing configuration of CDI, the dual-stage signal mixing configuration of CDI
provides significantly better visualization of the anatomical information while providing similar delineation between cancerous tissue and healthy tissue.
Illustrative examples of digital images acquired through the dual-stage signal mixing configuration of CDI are shown in FIGS. 8(a) to 8(c), in which potentially cancerous tissue is again highlighted with arrows.
Now referring to FIG. 9, shown is a schematic block diagram of a generic computing device. A
suitably configured computer device, and associated communications networks, devices, software and firmware may provide a platform for enabling one or more embodiments as described above. By way of example, FIG.9 shows a generic computer device 1000 that may include a central processing unit ("CPU") 1002 connected to a storage unit 1004 and to a random access memory 1006. The CPU 1002 may process an operating system 1001, application program 1003, and data 1023. The operating system 1001, application program 1003, and data 1023 may be stored in storage unit 1004 and loaded into memory 1006, as may be required.
Computer device 1000 may further include a graphics processing unit (GPU) 1022 which is operatively connected to CPU 1002 and to memory 1006 to offload intensive image processing calculations from CPU 1002 and run these calculations in parallel with CPU
1002. An operator 10010 may interact with the computer device 1000 using a video display 1008 connected by a video interface 1005, and various input/output devices such as a keyboard 1010, pointer 1012, and storage 1014 connected by an I/O interface 1009. In known manner, the pointer 1012 may be configured to control movement of a cursor or pointer icon in the video display 1008, and to operate various graphical user interface (GUI) controls appearing in the video display 1008. The computer device 1000 may form part of a network via a network interface 1011, allowing the computer device 1000 to communicate with other suitably configured data processing systems or circuits. One or more different types of sensors 1030 connected via a sensor interface 1032 may be used to search for and sense input from various sources. The sensors 1030 may be built directly into the generic computer device 1000, or optionally configured as an attachment or accessory to the generic computer device 1000.
Thus, in an aspect, there is provided a computer-implemented method operable on a computer device for detection and localization of tissue of interest, comprising:
receiving multiple signal acquisitions of imaged tissue from a subject; mixing in a first stage the acquired signals to obtain a joint correlation of diffusion signal attenuation in dependence upon the acquired signals; and displaying on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
In an embodiment, the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
In another embodiment, the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.
In another embodiment, obtaining a joint correlation of diffusion imaging comprises performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.
In another embodiment, the method further comprises mixing in the first stage one or more interpolated or extrapolated signals based on the acquired signals.

In another embodiment, the method further comprises comparing T2-weighted imaging results with that obtained using standard diffusion imaging via an apparent diffusion coefficient ADC);
and utilizing the comparison to obtain a joint correlation of diffusion signal attenuation.
In another embodiment, the method further comprises mixing in a second stage the obtained joint correlation of diffusion signal attenuation with a baseline signal containing anatomical information; and displaying the tissue of interest on the digital display, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
In another embodiment, the method further comprises providing one or more user-operable controls to allow a user to interact with the computer device to modify one or more image parameters, whereby the delineation of the tissue of interest on the digital display against the background of anatomical information may be further improved.
In another aspect, there is provided a system for detection and localization of tissue of interest, the system adapted to: receive multiple signal acquisitions of imaged tissue from a subject; mix in a first stage the acquired signals to obtain a joint correlation of diffusion signal attenuation in dependence upon the acquired signals; and display on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
In an embodiment, the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
In another embodiment, the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.
In another embodiment, obtaining a joint correlation of diffusion imaging comprises performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.

In another embodiment, the system is further adapted to mix in the first stage one or more interpolated or extrapolated signals based on the acquired signals.
In another embodiment, the system is further adapted to: compare T2-weighted imaging results with that obtained using standard diffusion imaging via an apparent diffusion coefficient ADC);
and utilize the comparison to obtain a joint correlation of diffusion signal attenuation.
In another embodiment, the system is further adapted to: mix in a second stage the obtained joint correlation of diffusion signal attenuation with a baseline signal containing anatomical information; and display the tissue of interest on the digital display, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
In another embodiment, the system is further adapted to: provide one or more user-operable controls to allow a user to interact with the system to modify one or more image parameters, whereby the delineation of the tissue of interest on the digital display against the background of anatomical information may be further improved.
In another aspect, there is provided a computer program product comprising a computer-readable medium, the computer program product storing code executable on a computer device for detection and localization of tissue of interest, the computer program product comprising: code for receiving multiple signals acquisitions of imaged tissue from a subject;
code for mixing in a first stage the acquired signals to obtain a joint correlation of diffusion signal attenuation in dependence upon the acquired signals; and code for displaying on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
In an embodiment, the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
In another embodiment, the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.

In another embodiment, the computer program product further comprises code for obtaining a joint correlation of diffusion imaging by performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.
While illustrative embodiments have been described above by way of example, it will be appreciated that various changes and modifications may be made without departing from the scope of the invention, which is defined by the following claims.

Claims (20)

1. A computer-implemented method operable on a computer device for detection and localization of tissue of interest, comprising:
receiving multiple signal acquisitions of imaged tissue from a subject; mixing in a first stage the acquired signals;
comparing T2-weighted imaging results with that obtained using standard diffusion imaging via an apparent diffusion coefficient (ADC);
utilizing the comparison to obtain a joint correlation of diffusion signal attenuation; and displaying on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
2. The computer-implemented method of claim 1, wherein the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
3. The computer-implemented method of claim 1, wherein the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.
4. The computer-implemented method of claim 1, wherein obtaining a joint correlation of diffusion imaging comprises performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.
5. The computer-implemented method of claim 1, further comprising mixing in the first stage one or more interpolated or extrapolated signals based on the acquired signals.
6. The computer-implemented method of claim 1, further comprising:
mixing in a second stage the obtained joint correlation of diffusion signal attenuation with a baseline signal containing anatomical information; and displaying the tissue of interest on the digital display, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
7. The computer-implemented method of claim 6, further comprising:
providing one or more user-operable controls to allow a user to interact with the computer device to modify one or more image parameters, whereby the delineation of the tissue of interest on the digital display against the background of anatomical information is further improved.
8. A computer-implemented system for detection and localization of tissue of interest, the computer-implemented system having a processor and a memory adapted to;
receive multiple signal acquisitions of imaged tissue from a subject;
mix in a first stage the acquired signals;
compare T2-weighted imaging results with that obtained using standard diffusion imaging via an apparent diffusion coefficient (ADC);
utilize the comparison to obtain a joint correlation of diffusion signal attenuation; and display on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
9. The system of claim 8, wherein the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
10. The system of claim 8, wherein the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.
11. The system of claim 8, wherein obtaining a joint correlation of diffusion imaging comprises performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.
12. The system of claim 8, wherein the system is further adapted to mix in the first stage one or more interpolated or extrapolated signals based on the acquired signals.
13. The system of claim 8, wherein the system is further adapted to:
mix in a second stage the obtained joint correlation of diffusion signal attenuation with a baseline signal containing anatomical information; and display the tissue of interest on the digital display, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
14. The system of claim 13, wherein the system is further adapted to:
provide one or more user-operable controls to allow a user to interact with the system to modify one or more image parameters, whereby the delineation of the tissue of interest on the digital display against the background of anatomical information is further improved.
15. A computer program product comprising a non-transitory computer-readable medium, the computer program product storing code executable on a computer device for detection and localization of tissue of interest, the computer program product comprising;
code for receiving multiple signals acquisitions of imaged tissue from a subject;
code for mixing in a first stage the acquired signals; code for comparing T2-weighted imaging results with that obtained using standard diffusion imaging via an apparent diffusion coefficient (ADC);
code for utilizing the comparison to obtain a joint correlation of diffusion signal attenuation:
and code for displaying on a digital display a resulting delineation of the imaged tissue of interest and background tissue based on the joint correlation of diffusion signal attenuation.
16. The computer program product of claim 15, wherein the multiple signal acquisitions are obtained from one of a magnetic resonance imaging (MRI) device, a computed tomography device, and an ultrasound device.
17. The Computer program product of claim 15, wherein the multiple signal acquisitions measure signal attenuation of the imaged tissue at different water diffusion motion sensitivities across multiple gradient pulse strengths and timings.
18. The computer program product of claim 15, further comprising code for obtaining a joint correlation of diffusion imaging by performing quantitative evaluation of the multiple signals using one or more of receiver operating characteristic (ROC) curve analysis, tissue class separability analysis, and visual assessment.
19. A computer-implemented method operable on a computer device for detection and localization of tissue of interest, comprising:
receiving multiple signal acquisitions of imaged tissue from a subject; mixing in a first stage the acquired signals to obtain a processed diffusion signal;
mixing in a second stage the obtained processed diffusion signal with a baseline signal;
and displaying the tissue of interest on a digital display based on the mixed processed diffusion signal and baseline signal, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
20. A computer-implemented system for detection and localization of tissue of interest, the computer-implemented system having a processor and a memory adapted to:
receive multiple signal acquisitions of imaged tissue from a subject; mix in a first stage the acquired signals to obtain a processed diffusion signal;
mix in a second stage the obtained processed diffusion signal with a baseline signal; and display the tissue of interest on a digital display based on the mixed processed diffusion signal and baseline signal, whereby the tissue of interest is delineated against a background of anatomical information to provide a visualization of the localization of the tissue of interest relative to the anatomical features of the subject.
CA2854844A 2014-06-19 2014-06-19 Correlated diffusion imaging system and method for identification of biological tissue of interest Active CA2854844C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CA2854844A CA2854844C (en) 2014-06-19 2014-06-19 Correlated diffusion imaging system and method for identification of biological tissue of interest

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CA2854844A CA2854844C (en) 2014-06-19 2014-06-19 Correlated diffusion imaging system and method for identification of biological tissue of interest

Publications (2)

Publication Number Publication Date
CA2854844A1 CA2854844A1 (en) 2015-12-19
CA2854844C true CA2854844C (en) 2018-10-02

Family

ID=54851578

Family Applications (1)

Application Number Title Priority Date Filing Date
CA2854844A Active CA2854844C (en) 2014-06-19 2014-06-19 Correlated diffusion imaging system and method for identification of biological tissue of interest

Country Status (1)

Country Link
CA (1) CA2854844C (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102016202254B4 (en) * 2016-02-15 2017-11-30 Siemens Healthcare Gmbh Model-free determination of image areas with anomalous diffusion based on diffusion-weighted magnetic resonance image data
DE102016213062A1 (en) 2016-07-18 2018-01-18 Siemens Healthcare Gmbh Determination of diffusion-weighted image data with a magnetic resonance device

Also Published As

Publication number Publication date
CA2854844A1 (en) 2015-12-19

Similar Documents

Publication Publication Date Title
Khalvati et al. Automated prostate cancer detection via comprehensive multi-parametric magnetic resonance imaging texture feature models
Nakamura et al. The mean apparent diffusion coefficient value (ADCmean) on primary cervical cancer is a predictive marker for disease recurrence
Kozlowski et al. Combined prostate diffusion tensor imaging and dynamic contrast enhanced MRI at 3T—quantitative correlation with biopsy
Reig et al. Breast MRI for evaluation of response to neoadjuvant therapy
Meng et al. Comparison between population average and experimentally measured arterial input function in predicting biopsy results in prostate cancer
Shu et al. Prediction of efficacy of neoadjuvant chemoradiotherapy for rectal cancer: the value of texture analysis of magnetic resonance images
Fedorov et al. A comparison of two methods for estimating DCE-MRI parameters via individual and cohort based AIFs in prostate cancer: a step towards practical implementation
Fan et al. Diagnosis of suspicious breast lesions using an empirical mathematical model for dynamic contrast-enhanced MRI
EP3307173B1 (en) System for identifying cancerous tissue
Long et al. MRI-based traditional radiomics and computer-vision nomogram for predicting lymphovascular space invasion in endometrial carcinoma
Borren et al. Accurate prostate tumour detection with multiparametric magnetic resonance imaging: dependence on histological properties
Crombé et al. High‐grade soft‐tissue sarcomas: Can optimizing dynamic contrast‐enhanced MRI postprocessing improve prognostic radiomics models?
Thoeny Diffusion-weighted MRI in head and neck radiology: applications in oncology
Villers et al. MRI in addition to or as a substitute for prostate biopsy: the clinician's point of view
Westphalen et al. Prostate cancer: prediction of biochemical failure after external-beam radiation therapy—Kattan nomogram and endorectal MR imaging estimation of tumor volume
US9619882B2 (en) Correlated diffusion imaging system and method for identification of biological tissue of interest
Priola et al. Comparison of CT and chemical-shift MRI for differentiating thymoma from non-thymomatous conditions in myasthenia gravis: value of qualitative and quantitative assessment
Goto et al. Adding a model-free diffusion MRI marker to BI-RADS assessment improves specificity for diagnosing breast lesions
Wong et al. Correlated diffusion imaging
Lo et al. Effect of imaging parameter thresholds on MRI prediction of neoadjuvant chemotherapy response in breast cancer subtypes
Yuan et al. A five-colour colour-coded mapping method for DCE-MRI analysis of head and neck tumours
Lambregts et al. Optimal imaging staging of rectal cancer
Mani et al. Early prediction of the response of breast tumors to neoadjuvant chemotherapy using quantitative MRI and machine learning
Créhange et al. Tumor volume and metabolism of prostate cancer determined by proton magnetic resonance spectroscopic imaging at 3T without endorectal coil reveal potential clinical implications in the context of radiation oncology
Vos et al. Role of serial multiparametric magnetic resonance imaging in prostate cancer active surveillance

Legal Events

Date Code Title Description
EEER Examination request

Effective date: 20170901