Patent Grant
7560926
This invention relates generally to Magnetic resonance imaging (MRI) and more particularly to diffusion weighted MRI imaging (DWI). Still more particularly, the invention relates to a method for optimizing the b-values to be used in diffusion weighted magnetic resonance imaging.
As is known in the art, Diffusion weighted imaging (DWI) is a method used in magnetic resonance imaging (MRI) and offers a means to evaluate an area of anatomy in terms of the motion of water molecules. The level of water diffusivity can provide an indication of the structure of the tissue at the cellular level. This in turn can provide an indication of the status of the tissue. In terms of tumors within an organ such as, for example, the liver, the water diffusivity within the tumor will be less than that of the healthy tissue of the organ. This is because the cancerous cells are denser, with more impediments to water motion. In a tumor which has undergone treatment to kill off the cancerous cells, the necrotic tumor will show an increase in diffusivity compared to the viable tumor. This is because, in killing the cancerous cells, the membranes are broken down, allowing greater motion of water molecules. Thus, measurement of water diffusivity can serve as a surrogate marker to evaluate tumor treatment response.
In DWI, the MR pulse sequence is designed in such a way as to produce higher signal intensity in proportion to the water diffusivity. The sequence's b-value parameter, which is related to the MR magnetic field gradient amplitude, is a key factor in controlling the diffusion weighting. The higher the b-value, the stronger the diffusion weighting. Conversely, the signal intensity falls in an exponential manner as a function of b-value. Therefore, the higher the b-value, the lower the signal intensity, and the higher the noise. Thus, there is a tradeoff in choosing the b-value such that it gives a strong diffusion weighting, but yet still has a high signal to noise ratio (SNR).
For a certain set of tissue types, the water diffusivity difference between the tissues will be most visible using a certain b-value. While there are some common guidelines as to what b-values should be used for certain tissue types, the results will vary by patient and other scan parameters and conditions. In general, the optimum b-value is considered to be that at which the maximum contrast-to-noise ratio is achieved. Since the optimum b-value for a certain patient, anatomy and scan is not known, the level of diffusivity is usually reported in terms of the apparent diffusion coefficient (ADC), which describes the rate of signal decay as a function of b-value and can be calculated from two or more images with different b-values (considering the slope of the line described by ln(SI(b)/SI(bo)) where bo reflects no diffusion weighting); where SI is signal intensity and In is the natural logarithm. In this way it is not necessary that any single scan is using the optimum b-value. By performing the ADC calculation on a pixel by pixel basis the results are presented as a parametric image often referred to as an ADC map.
However, depending on the b-values used, and the uncertainty in signal intensity in each image at a particular b-value, the resulting ADC values will have some degree of uncertainty, which decreases their diagnostic value, particularly in terms of doing follow-up studies where a comparison is being made between ADC values pre and post treatment.
In DWI studies usually multiple slices are being acquired for each dataset. For many cases, it is also desirable that the complete dataset be captured during one breath-hold to avoid motion of the anatomy and the need for registration within a single dataset. For this reason, and because each scan takes a fair amount of time, current practice is to make the best judgment in terms of what b-values to use to obtain the best ADC maps, and keep the number of scans with different b-values to a minimum.
In accordance with the present invention, a method is provided for selecting a set of b-values for a diffusion weighted scan. The method includes: identifying a single reference slice within a volume of interest (VOI); progressively changing b-values over a plurality of scans of the reference slice within the VOI wherein each one of the scans has a different b-value to obtain a plurality of diffusion weighted images; calculating apparent diffusion coefficient (ADC) maps for combinations of b-values from the obtained plurality of diffusion weighted images; evaluating the calculated ADC maps; selecting from the evaluation optimal b-values; and using the selected b-values for subsequent scans.
In one embodiment, a method is provided for selecting a set of b-values for a diffusion weighted scan. The method includes: identifying a single reference slice having a position and orientation providing a representative view of the volume of interest (VOI); performing n scans, where n is an integer greater than 1, on the reference slice at progressively increasing b-values to provide a set of pre-scan reference images; interactively identifying a region of interest (ROI), or multiple ROI, within the reference slice and relevant regions outside of the ROIs for reference, with the ROIs identified in one image being propagated to the other images of different b-values; calculating ADC maps for combinations of b-values as such combinations are available from the set of pre-scan reference images, including calculating a subset of the maps based on only two b-values, and other maps being based on more than two b-values; evaluating each combination of b-values according to a plurality of criteria; and selecting the b-values based on the evaluation.
In one embodiment, the criterion includes one or more of: goodness of fit parameters in the ADC map calculations; error in predicting SI for each b-value image; SNR within the b-value image; SNR within the ROI; contrast between the ROI and surrounding regions; similarity measures in the ROI and surrounding regions.
In one embodiment, the evaluation includes a weighting scheme based on measures pertaining to complete ADC maps and b-value images: ADC goodness of fit, ADC error of prediction and b-value image SNR.
In one embodiment, the weighting scheme includes giving greater weight to pixels within the ROI and near to the ROI based on the level of SNR and contrast in the vicinity of the pixel in the b-value images.
The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
The single FIGURE is a flowchart of the method for selecting a set of b-values for a diffusion weighted scan according to the invention.
Like reference symbols in the various drawings indicate like elements.
Referring now to the single FIGURE, a flowchart is shown of a method for selecting a set of b-values for a diffusion weighted scan. The method includes: identifying a single reference slice within a volume of interest (VOI); progressively changing b-values over a plurality of scans of slices within the VOI wherein each one of the scans has a different b-value to obtain a plurality of diffusion weighted images; calculating ADC maps for combinations of b-values from the obtained plurality of diffusion weighted images; evaluating the calculated ADC maps; selecting from the evaluation an optimal set of b-values; and using the selected set of b-values for subsequent scans.
At the start of a DWI study, certain pre-scans are performed to provide input to the b-value optimization software application. After some user input in identifying the region(s) of interest, the application searches for the optimum b-values. After determining the optimum values as judged by certain evaluation criteria those values are automatically or manually set in preparation for the next scans in the protocol.
The method is described in greater detail as follows:
A key component of the method is the evaluation algorithm which scores the ADC maps. Within that algorithm, the weighting scheme helps to minimize the effects of the ROI segmentation sometimes being less than perfect.
Described below are various optimization algorithms:
Each combination of b-values will be evaluated according to multiple criteria. Some of these criteria will relate to the ADC calculation itself, such as the goodness of fit parameters and the error in predicting SI for each b-value. As the SNR and contrast in an ADC map are dependent on those characteristics in the different b-value images from which the ADC map is derived, the evaluation algorithm will incorporate a weighting scheme based on the measurements of SNR and contrast in the b-value images. Since ADC is not expected to change regionally from pixel to pixel, similarity measures for adjacent pixels could also be used as part of the evaluation.
For each combination of b values, the overall image evaluation score is:
Where weight_SNR and weight_CNR will be based on those measures in the b value images.
The overall image weighted score and the ROI weighted scores are combined into a global score for each combination of b-values, for example:
and the combination of b-values with the highest global weighted score is selected for subsequent scans.
A further adaptation of the algorithm in Example 1 is to set a range of acceptable global scores, and select the global score within that range that employs the smallest number of b-values to minimize subsequent scan time.
Alternatively, the number of b-values used for a particular ADC map could be added to the weighting scheme described above.
Alternatively linear optimization methods are employed to minimize the number of b-values used (to minimize scan time), while simultaneously maximizing one or more of the other computed metrics such as SNR.
Because the pre-scans are only being done on one slice, many different b-values can be explored in a relatively short time. Once the optimum b-values have been determined, they are then used for a complete multi-slice scan.
DWI studies will usually report quantitative results, such as the overall ADC value within an ROI, this being a key measure in evaluating the tissue and comparing its state between studies. The accuracy and reproducibility of these measures are of primary concern. This invention provides a quantitative approach for determining the b-values to be used to achieve the most accurate results.
A number of embodiments of the invention have been described. Nevertheless, it will be understood that various modifications may be made without departing from the spirit and scope of the invention. Accordingly, other embodiments are within the scope of the following claims.
This application claims priority from U.S. Provisional application No. 60/979,892 filed Oct. 15, 2007, the entire subject matter thereof being incorporated herein by reference.
| Number | Name | Date | Kind |
|---|---|---|---|
| 5786692 | Maier et al. | Jul 1998 | A |
| 6163152 | Bernstein et al. | Dec 2000 | A |
| 6670812 | Mock et al. | Dec 2003 | B1 |
| 7310548 | Van Den Brink | Dec 2007 | B2 |
| 7511494 | Wedeen | Mar 2009 | B2 |
| Number | Date | Country | |
|---|---|---|---|
| 20090096448 A1 | Apr 2009 | US |
| Number | Date | Country | |
|---|---|---|---|
| 60979892 | Oct 2007 | US |
