CROSS REFERENCE TO RELATED APPLICATIONS

[00011 This patent disclosure claims the benefit of U.S. Provisional Application No. 61/565,578 filed on December 1 , 2011 , which is hereby wholly incorporated by reference.

BACKGROUND

[0002] MRI is a diagnostic imaging technique in which image contrasts may be generated. Different image contrasts may be sensitive to different pathologies. Two image contrasts for diagnostic medicine are Ti contrast and T ₂ contrast. T ₂ contrast corresponds to the rate at which magnetization disappears perpendicular to the main magnetic field. T ₂ contrast is useful in clinical imaging. T ₂ contrast is often used to facilitate identification of tumors, tissue affected by stroke and other lesions, pathologies, and other differences between normal and diseased tissue.

[0003] For example, damaged myelin can be distinguished from undamaged myelin based on differences in the T ₂ decay. One type of damage to myelin, referred to as demyelination, impairs the conduction of signal for an affected neuron. In addition to causing degradation of a neuron's signal, demylenation results in bundle migration. Bundle migration refers to the distance between nerve fibers in a nerve bundle. Some brain diseases including multiple sclerosis (MS), leukoencephalopathy (LE), are related to loss of myelin water (MW). Myelin water fraction (MWF) can serve as a direct indicator of myelin integrity and provide quantitative measurements of myelin structure and change due to white matter diseases.

[0004] T ₂ contrast is typically generated using a spin-echo (SE) or turbo-spin- echo (TSE) sequence. A single spin-echo is formed through a combination of a 90 degree pulse and a 180 degree pulse. The 90 degree pulse generates signal perpendicular to the main magnetic field and is followed by the 180 degree pulse that refocuses transverse magnetization to form an echo. A purely T ₂ weighted spin echo signal may be acquired at a point in time after the 180 degree pulse. The time at which the echo is acquired, the echo time (TE), occurs at a time after the 180 degree pulse that is equal to the spacing between the 90 degree pulse and the 180 degree pulse.

[0005] Conventionally, a T ₂ decay curve may be measured using a single- slice multi-echo Carr-Purcell-Meiboom-Gill (CPMG) sequence with nonselective composite refocusing pulses and large z-axis gradient crushers. However, the single-slice coverage, long acquisition time (e.g., 25 minutes) due to nonselective refocusing pulses, and increasing signal-to-noise ratio (SNR) are obstacles to practical clinical applications. To increase volume coverage, a conventional multi-slice CPMG sequence with slice-selective refocusing pulse may be used. However, due to intrinsic sensitivity to refocusing imperfection, the measured curve deviates from the true T ₂ decay curve, which may be pronounced for the first several points.

[00061 In another conventional method, a T ₂ preparation spiral imaging sequence is used. However, this method yielded MWF maps with low resolution and poor image quality. Multi-slice techniques based on T ₂ * signal decay have also been proposed to achieve larger volume coverage and shorter acquisition time. However, high sensitivity to local field inhomogeneity is a limitation of this technique for practical applications since local field gradient along a slice direction can introduce non-exponential signal decay. Furthermore, a signal void region is observed in MWF maps even when high-order shim is performed to minimize field inhomogeneity before data acquisition and the correction of field inhomogeneity is performed in the post-processing.

BRIEF DESCRIPTION OF THE DRAWINGS

|0007] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate various example systems, methods, and other example embodiments of various aspects of the invention. It will be appreciated that the illustrated element boundaries (e.g., boxes, groups of boxes, or other shapes) in the figures represent one example of the boundaries. One of ordinary skill in the art will appreciate that in some examples one element may be designed as multiple elements or that multiple elements may be designed as one element. In some examples, an element shown as an internal component of another element may be implemented as an external component and vice versa. Furthermore, elements may not be drawn to scale.

[0008] Figure 1 illustrates a method associated with the spectral analysis of T ₂ spectral data.

[0009] Figure 2 illustrates an example of a normal nerve bundle and an abnormal nerve bundle.

[0010] Figure 3 illustrates an example data acquisition scheme with respect to the spectral analysis of T ₂ spectral data.

[0011] Figure 4A illustrates example refocusing slice profiles associated with T ₂ spectral data acquisition.

[0012] Figure 4B illustrates example T ₂ decay curves representative of acquired T ₂ spectral data.

[0013] Figure 5 illustrates graphs associated with a T ₂ decay curve and T ₂ spectral peaks. [0014] Figure 6A illustrates T ₂ decay curves of two adjacent example pixels.

[0015] Figure 6B illustrates the smoothing of T ₂ spectral data with varying values of the regularization constant, μ.

[0016] Figure 7 illustrates example intervals associated with spectral peaks.

[0017] Figure 8 illustrates example T ₂ sampling points associated with determining myelin water fraction.

[0018] Figure 9 illustrates an image associated with a conventional refocusing pulse and an image associated with a refocusing pulse that is three times larger.

[0019] Figure 10 illustrates images associated with non-weighted regularized T ₂ spectral data and images associated with weighted regularized T ₂ spectral data.

[0020] Figure 11 illustrates images associated with a healthy brain and images associated with a brain affected by lesions.

[0021] Figure 12 illustrates an apparatus configured to perform spectral analysis of T ₂ spectral data.

DEFINITIONS

[0022] The following includes definitions of selected terms employed herein. The definitions include various examples and/or forms of components that fall within the scope of a term and that may be used for implementation. The examples are not intended to be limiting. Both singular and plural forms of terms may be within the definitions.

[0023] References to "one embodiment", "an embodiment", "one example", "an example", and so on, indicate that the embodiment(s) or example(s) so described may include a particular feature, structure, characteristic, property, element, or limitation, but that not every embodiment or example necessarily includes that particular feature, structure, characteristic, property, element or limitation. Furthermore, repeated use of the phrase "in one embodiment" does not necessarily refer to the same embodiment, though it may.

[0024] "Computer component", as used herein, refers to a computer-related entity (e.g., hardware, firmware, software, software in execution, combinations thereof). Computer components may include, for example, a process running on a processor, a processor, an object, an executable, a thread of execution, a program, and a computer. A computer component(s) may reside within a process and/or thread. A computer component may be localized on one computer and/or may be distributed between multiple computers.

[0025] "Computer-readable medium", as used herein, refers to a non- transitory medium that stores signals, instructions and/or data. A computer- readable medium may take forms, including, but not limited to, non-volatile media, and volatile media. Non-volatile media may include, for example, optical disks, magnetic disks, and so on. Volatile media may include, for example, semiconductor memories, dynamic memory, and so on. Common forms of a computer-readable medium may include, but are not limited to, a floppy disk, a flexible disk, a hard disk, a magnetic tape, other magnetic medium, an ASIC, a CD, other optical medium, a RAM, a ROM, a memory chip or card, a memory stick, and other media from which a computer, a processor or other electronic device can read.

[0026] "Logic", as used herein, includes but is not limited to hardware, firmware, software and/or combinations thereof to perform a function(s) or an action(s), and/or to cause a function or action from another logic, method, and/or system. Logic may include a software controlled microprocessor, discrete logic (e.g., application specific integrated circuit (ASIC)), an analog circuit, a digital circuit, a programmed logic device, a memory device containing instructions, and so on. Logic may include a gate(s), a combinations of gates, other circuit components, and so on. In some examples, logic may be fully embodied as software. Where multiple logical logics are described, it may be possible in some examples to incorporate the multiple logical logics into one physical logic. Similarly, where a single logical logic is described, it may be possible in some examples to distribute that single logical logic between multiple physical logics.

[0027] Operable connection", or a connection by which entities are "operably connected", is one in which signals, physical communications, and/or logical communications may be sent and/or received. An operable connection may include a physical interface, an electrical interface, and/or a data interface. An operable connection may include differing combinations of interfaces and/or connections sufficient to allow operable control. For example, two entities can be operably connected to communicate signals to each other directly or through one or more intermediate entities (e.g., processor, operating system, logic, software). Logical and/or physical communication channels can be used to create an operable connection.

[0028] "Software", as used herein, includes but is not limited to, one or more computer instructions and/or processor instructions that can be read, interpreted, compiled, and/or executed by a computer and/or processor. Software causes a computer, processor, or other electronic device to perform functions, actions and/or behave in a desired manner. Software may be embodied in various forms including routines, algorithms, modules, methods, threads, and/or programs. In different examples software may be embodied in separate applications and/or code from dynamically linked libraries. In different examples, software may be implemented in executable and/or loadable forms including, but not limited to, a stand-alone program, an object, a function (local and/or remote), a servelet, an applet, instructions stored in a memory, part of an operating system, and so on. In different examples, computer-readable and/or executable instructions may be located in one logic and/or distributed between multiple communicating, cooperating, and/or parallel processing logics and thus may be loaded and/or executed in serial, parallel, massively parallel and other manners.

[0029] Suitable software for implementing various components of example systems and methods described herein may be developed using programming languages and tools (e.g., Java, C, C#, C++, SQL, APIs, SDKs, assembler). Software, whether an entire system or a component of a system, may be embodied as an article of manufacture and maintained or provided as part of a machine-readable medium.

[00301 "Signal", as used herein, includes but is not limited to, electrical signals, optical signals, analog signals, digital signals, data, computer instructions, processor instructions, messages, a bit, a bit stream, or other means that can be received, transmitted and/or detected.

[0031] "User", as used herein, includes but is not limited to, one or more persons, software, computers or other devices, or combinations of these.

[0032] Some portions of the detailed descriptions that follow are presented in terms of algorithms and symbolic representations of operations on data bits within a memory. These algorithmic descriptions and representations are used by those skilled in the art to convey the substance of their work to others. An algorithm, here and generally, is conceived to be a sequence of operations that produce a result. The operations may include physical manipulations of physical quantities. Usually, though not necessarily, the physical quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a logic, and so on. The physical manipulations create a concrete, tangible, useful, real-world result.

[0033] It has proven convenient at times, principally for reasons of common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, and so on. It should be borne in mind, however, that these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise, it is appreciated that throughout the description, terms including processing, computing, determining, and so on, refer to actions and processes of a computer system, logic, processor, or similar electronic device that manipulates and transforms data represented as physical (electronic) quantities.

[0034] Example methods may be better appreciated with reference to flow diagrams. While for purposes of simplicity of explanation, the illustrated methodologies are shown and described as a series of blocks, it is to be appreciated that the methodologies are not limited by the order of the blocks, as some blocks can occur in different orders and/or concurrently with other blocks from that shown and described. Moreover, less than all the illustrated blocks may be required to implement an example methodology. Blocks may be combined or separated into multiple components. Furthermore, additional and/or alternative methodologies can employ additional, not illustrated blocks.

DETAILED DESCRIPTION

[0035] Myelin is a major substance of the white matter in human brain. Myelin forms an electrically insulating layer to support and protect the axon of a neuron for the proper functioning of the nervous system. Myelin water is quantified with a multi-interval analysis of a T ₂ spectrum from a T ₂ decay curve. The intervals in human brain may include (1) myelin water, having a T ₂ component with a T ₂ range from 15 to 40 ms at 3.0 T; (2) tissue water, having a T ₂ component with a T ₂ range from 40 to 200 ms; (3) long T ₂ tissue water, having a T ₂ component with a T ₂ range from 200 to 800 ms, and (4) cerebrospinal fluid, having a T ₂ component with a T ₂ range from 800 to 2000 ms. The ratios of the components to the summation of the components are termed myelin water fraction ( WF), tissue water fraction (TWF), long T ₂ tissue water fraction (LWF), and cerebral spinal fluid fraction (CSF).

[0036] Example apparatus and methods are associated with the spectral analysis of magnetic resonance imaging (MRI) T ₂ spectral data associated with myelin. This method is called T ₂ SPectrum Analysis using a weighted Regularized Approach and non-loCal mean filter (T ₂ SPARC). In one embodiment, a modified multi-slice CPMG sequence is employed with a refocusing slice whose thickness is three times larger than an excitation slice. The size of the refocusing slice is configured with respect to the excitation slice to diminish deviations from the curve. Using a refocusing slice thickness that is three times larger than the excitation slices facilitates achieving a balance between acquisition time and quality of T ₂ decay curves. Additionally, since SNR affects the quality of MWF maps, in one embodiment, a weighted regularized approach with a large regularization coefficient is used in combination with a nonlocal mean (NLM) de-noising method to reduce the effect of noise and lessen requirement for the quality of T ₂ decay curve. While a refocusing slice three times larger than an excitation slice is described, one of ordinary skill in the art will appreciate that other ratios (e.g., 3.1 to 1 , 2.9 to 1 ) between the refocusing slice and the excitation slice may be employed.

[0037] Figure 1 illustrates a method 100 associated with T ₂ spectral analysis for myelin water imaging. Method 100 includes, at 1 10, accessing T ₂ spectral data from an object prepared by a pulse sequence and applied by an MRI apparatus. In one embodiment, the object may be a brain. Ten axial brain slices may be acquired over two scans. The slice position of the second scan may be shifted relative to the first scan so that the first scan and the second scan are interleaved. While ten slices and two acquisitions are described, one of ordinary skill in the art will appreciate that a greater and/or lesser number of slices and a great and/or lesser number of scans may be employed. [0038] In conventional T ₂ imaging, a first portion of the T ₂ spectral data is misrepresented due to a portion of the slice not being subject to an adequate refocusing pulse. Deviations in the first portion of the T ₂ decay curve have a profound effect on the MWF because MWF is derived from the T ₂ spectral data in the first interval from 15 to 40ms. Using a modified refocusing pulse with a slice thickness three times larger than the excitation pulse facilitates acquiring initial T ₂ spectral data with greater accuracy.

[0039] Method 100 also includes, at 120, controlling the MRI apparatus to perform a spectral analysis of the T ₂ spectral data. The spectral analysis includes uniformly weighting the T ₂ spectral data, using a pre-determined regularization constant, and using T ₂ spectral component thresholds to define the multiple intervals of myelin water imaging. In one example, the spectral analysis of the T ₂ spectral data includes, at 121 , decomposing the T ₂ spectral data using multi-exponential functions to produce T ₂ spectra. The multi-exponential functions may be employed using, for example, a non-negative least squares (NNLS) approach.

[0040] To improve the reliability of the result of the multi-exponential functions, the spectral analysis at 120, includes, at 122, regularizing the T ₂ spectral amplitude by applying a regularized constant μ. Conventionally, a small regularization constant may be used. However, a sufficiently large regularized constant, μ, can smoothen the T ₂ spectra so adjacent pixels of the T ₂ spectral data have similar characteristics.

[0041] If the T ₂ spectrum is uniformly sampled, the regularization is uniformly weighted for all spectral amplitudes. However, if T ₂ spectrum is logarithmically sampled, the regularization puts more weight on spectral amplitudes that correspond to larger log-spaced time intervals. The logarithmically sampled weighting disparity is further aggravated the larger the regularized constant, μ.

[0042] To ameliorate the logarithmically sampled weighting disparity, the example spectral analysis at 120 includes, at 123, weighting the multi- exponential functions. Weighting is performed so that in one example the regularization is uniformly weighted for the spectral amplitudes. The multi- exponential functions may be weighted with an element representative of inverse of log-spaced T ₂ time intervals. For example, at 123, an inverse of log-spaced T ₂ time interval matrix may be applied to the multi-exponential functions.

[0043] A plurality of T ₂ spectra result from applying weighted regularized multi-exponential functions to the T2 spectral data. The example spectral analysis 120 includes, at 124, partitioning the T ₂ spectra into intervals by fixed thresholds. In different examples, the intervals may correspond to MWF, TWF, LWF, and CSF. An individual fraction may correspond to an interval. In another embodiment, the intervals may be defined by determining thresholds.

[0044] The example spectral analysis 120 includes, at 125, generating at least one parametric map. For an individual pixel of a plurality of pixels, the ratios of at least one of MW, TW, LW, and CS components to the summation of all components are computed to determine the at least one corresponding to MWF, TWF, LWF, and CSF. The computed value of MWF, TWF, LWF, or CSF for the plurality of pixels may be used to generate the corresponding parametric map.

[0045] The example spectral analysis at 120 includes, at 126, denoising the at least one parametric map. In one example, the parametric map may be de- noised with a Non Local Mean (NLM) filter. The NLM filter is applied after the T ₂ spectral analysis (e.g., wrNNLS). Method 100 also includes, at 130, providing the at least one parametric map.

[0046] The method 100, as illustrated in Figure 1 , may be used to measure myelin water in white matter. Figure 2 illustrates an example of normal white matter 200 and abnormal white matter 250. The white matter 200 comprises a plurality of nerve fibers. A nerve fiber 210 has an axon 220. The axon 220 is an extended fiber of a nerve cell that carries a nerve signal. The axon 220 is surrounded by myelin 230. In normal white matter 200 the nerve fibers are grouped closely together. In abnormal white matter 250, the nerve fibers, such as nerve fiber 260, are not grouped together. Instead, the nerve fibers are subject to bundle migration, meaning the nerve fibers are moving away from one another.

[0047] Figure 3 illustrates an example data acquisition scheme 300 associated with the acquisition of T ₂ spectral data. In one example, a total of ten brain axial slices from two scans may be acquired using a modified pulse sequence on a 3T scanner using 12-channel head coil. A first scan 310 that includes five axial brain slices 311 , 312, 313, 314, and 315 is performed. Subsequently, a second scan 320 that also includes five axial brains slices 321 , 322, 323, 324, and 325 is performed. In the example, the individual brain axial slices of the ten brain axial slices have a slice thickness of 5mm. The slice position of the second scan 320 is shifted 7.5 mm relative to the first scan 310 to form two interleaved scans of ten consecutive slices with 2.5 mm gap. One of ordinary skill in the art will appreciate that greater and/or lesser number of scans may be employed and that different slice shifts and different slice thicknesses may be employed.

[0048] The imaging protocol for the two scans in the example may include: FOV = 220x220 mm; slice thickness = 5 mm; acquisition matrix = 256x256; 5 slices with 200% slice gap; TEs = 10 N, (N= 1 , 2, 3, 32); TR = 3000 ms; receiver bandwidth = 300 Hz/pixel; GRAPPA reduction factor = 2 and the reference line = 24. In the example, the acquisition time for the first scan 310 or the second scan 320 may be 5:29 minutes with a total time of 1 1 minutes.

[0049] Figure 4A illustrates a graph 400 associated with the thickness of refocusing slice profiles. Once the number of scans and slices is established based on the imaging protocol, a pulse sequence is determined. Imperfections in the conventional refocusing slice profile 420 contribute to flawed T ₂ spectral data. A modified pulse sequence with a modified refocusing profile 430 reduces signal deviations as compared to conventional refocusing slice profile 420.

[0050] An excitation slice thickness is demarcated by dashed lines 410 and 415. In one embodiment the excitation slice thickness is 5mm. A conventional refocusing RF pulse has a conventional refocusing slice profile 420. Specifically, there is a significant amount of excitation slice, demarcated by dashed lines 410 and 415, with a refocusing pulse less than 180°. This imperfection leads to large deviations of the first several points in T ₂ decay curve. The first several points are significant to MWF data since the T ₂ spectrum corresponding to myelin water is measured in the first interval from 15 - 40ms.

[0051] A modified refocusing profile 430 increases the slice thickness of the refocusing RF to be larger than the excitation slices demarcated by dashed lines 410 and 415. When the slice thickness of the refocusing RF is larger then the excitation slice, then it may be possible to generate more accurate T ₂ spectral data. The increased slice thickness of the modified refocusing profile 430 reduces signal deviations. The slice thickness of the modified refocusing slice profile 430 may be three times larger than an excitation slice thickness.

[0052] Figure 4B illustrates a conventional measured T ₂ decay curve 460 as compared to a modified T ₂ decay curve 470. Due to intrinsic sensitivity to refocusing imperfections, a conventional measured T ₂ decay curve 460 deviates from the modified T ₂ decay curve 470, as illustrated for a first portion 465 of the conventional measured T ₂ decay curve 460. Previous attempts to ameliorate the conventional refocusing profile 420, as shown in Figure 4A, have included the addition of constant slice crushers around each refocusing pulse. However, initial values of the first portion 465 of the conventional measured T ₂ decay curve 460 were still underestimated.

[0053] The modified refocusing slice profile 430 leads to a substantially ameliorated T ₂ decay curve 470 with a corrected first portion 475. Accordingly, the first portion 475 of the ameliorated T ₂ decay curve 470 yields improved T ₂ spectrum corresponding to myelin water measured in the first interval from 15 - 40ms. Thus, the modified pulse sequence yields more accurate T ₂ spectral amplitude(s) corresponding to myelin water.

[0054] Figure 5 illustrates example graphs associated with a T ₂ decay curve and T ₂ spectral peaks. The T ₂ spectral data was acquired from two scans of five axial brain slices acquired using a pulse sequence with a modified refocusing slice profile. The data yields a T ₂ decay curve 510. Spectral peaks 531 , 532, 533, and 534 of the spectral graph 530 are determined using spectral analysis 520.

[0055] T ₂ decay signal, y(t _n ), of the T ₂ decay curve 510, can be described using a general integral of exponential functions: y(t„) = dT _z , n = ±, 2, - , JV (1 )

[0057] where t _n is the measured time point equal to n-TE, and TE is time of echo; S(T ₂ ) is the amplitude of the spectral component at relaxation time T ₂ . The range of T ₂ is from zeros to maximum T ₂ value T _max (e.g., 2 seconds). This integral equation can be discretized as a summation of piecewise amplitude constants over M small ranges of T ₂ : r n = 1,2, ·· · , /V, (2)

[0058] where M is the T ₂ sampling points and

S _Tm = S _m AT _2nL (3)

[0059] where AT ₂ is a logarithmic T ₂ time interval. Equation (4) is in a general matrix form.

Y = E S _T , E _nm = e ¾m Vi

[0060] The T ₂ spectrum can be solved from this equation. The T ₂ spectrum was partitioned into four intervals: myelin water, T ₂ =15 - 40 ms; tissue water, 40 - 200 ms; tissue water with long T ₂ component, 200 - 800 ms; cerebral spinal fluid, 800 - 2000 ms.

[0061] The spectral amplitudes S _Tm are nonnegative. Thus, an approach, for example, Non-Negative Least Squares (NNLS) approach may be used to perform T ₂ spectrum analysis 520. Regularization is used to improve the reliability of the solution. Therefore, a regularized approach (e.g., regularized Non-Negative Least Squares (rNNLS)) is used to solve equation (4) by minimizing:

X ¹ = min{|| E S _T - K|| ₂ + ju||5 _T || ₂ } , S _r > 0 (5)

[0062] where ΙΙΊΙ ₂ represents the Euclidean norm, and μ is a regularization constant. The larger the μ value, the more the T ₂ spectrum is smoothened. Conventionally, a small value of μ is used in a typical MWF study to satisfy the constraint of ί.02χΙ _ιη < χ ^ζ < 1.02S* „ where χ^ _(η is the minimum misfit of a solution (e.g., an NNLS solution). A smaller μ value has been conventionally used to avoid smoothing of the spectral peaks derived from the measured T ₂ decay curve because smoothing was viewed as distorting the T ₂ spectral data. However, due to use of a smaller μ, the regularized approach has been sensitive to noise. Therefore, conventionally, four averages in acquisition are used to double the SNR. Having to acquire four averages prolongs the total acquisition time to over 25 minutes for a single slice. This may be an unacceptably long acquisition time.

100631 The regularized term S _T in equation (5) includes the spectral amplitude 5 and T ₂ time interval as shown in equation (3). When T ₂ spectral data is logarithmically sampled, regularization puts more weight on spectral amplitudes that correspond to larger log-spaced time intervals. Therefore, in addition to being regularized, the T ₂ spectrum amplitude is weighted to ameliorate the problems caused by using logarithmically spaced T ₂ time intervals. The larger the T ₂ time interval, the less weighting. With a larger μ value, weighting can be applied as follows:

X* = min{||£ 5 _r - y|| _¾ - ii.|| »' 5 _T j|. _B }, if^ -J- , S _T ≥Q (6)

[0064] where W is a weighting matrix and its element is equal to the inverse of log-spaced T ₂ time intervals. Using equation (6), the spectral amplitudes are equally weighted in regularization even with logarithmically spaced T ₂ time intervals. [0065] In one example, after the weighted regularized multi-exponential functions are applied, a filter may be applied. In one embodiment, the filter is a non-local mean (NLM) filter. The noise is reduced by averaging a large number of similar pixels regardless of whether they are spatially close or not. The de- noising performance of an NLM filter is designed for the presence of similar image structure in contrast to other classical neighborhood filters such as anisotropic filter and Gaussian filter. Therefore, a NLM filter is suitable to de- noise parametric maps of MWF and TWF since they have a large area with similar pixels. In one embodiment, a NLM filter with degree of filtering and radiuses of search and similarity window equal to 10, 5 and 2, respectively was used to improve image quality after all MWF, TWF, LWF and CSF parametric maps were scaled by multiplying 1000.

[0066] Figure 6A illustrates T ₂ decay curves of two adjacent pixels acquired using an example data acquisition scheme. A first T ₂ decay curve 610 corresponds to a first adjacent pixel. A second T ₂ decay curve 620 corresponds to a second adjacent pixel. The first and second pixels are located in white matter of the brain and have similar characteristics due to their proximity. Therefore, two adjacent pixels should have approximately the same T ₂ spectra despite any differences in the amount of noise to which the adjacent pixels are subjected.

[0067] Figure 6B illustrates the smoothing of T ₂ spectral data with varying values of the regularization constant, μ. A μ value is determined for the first T ₂ decay curve 610 (shown in Figure 6A) of the first adjacent pixel and the of a second T ₂ decay curve 620 (shown in Figure 6A) of the second adjacent pixel. The first T ₂ decay curve 610 and the second T ₂ decay curve 620 are converted to T ₂ spectra using the weighted regularized approach (e.g., wrNNLS) described with respect to Figure 5.

[0068] The μ=0 graph 630 illustrates the spectral amplitudes of two adjacent pixels. The first pixel spectral peak corresponds to a small peak component 631 and a large peak component 632. Likewise, the second pixel spectral peak corresponds to a small peak component 638 and a large peak component 639. The first spectral peak and the second spectral peak and their components are distinct. Therefore, despite the first pixel and the second pixel being adjacent, without a regularization constant (e.g., μ=0) the first pixel spectral peak and the second pixel spectral peak of the μ=0 graph 630 have different spectral characteristics. The small peak components lead to MWF equal to 0% in one pixel and 9.8% in the other pixel. Thus, despite the proximity of the first adjacent pixel and the second adjacent pixel, the un-regularized results show different characteristics.

[0069] A regularization constant, μ, is applied to smooth the T ₂ spectra so that adjacent pixels of the T ₂ spectral data have similar spectral characteristics. In one example, a regularization constant sufficient to achieve a desired uniformity between the T ₂ spectra (e.g.,T ₂ spectra with 1 % deviation) may be applied. Alternatively, a regularization constant sufficient to avoid a threshold difference (e.g., a 1% difference) in a spectral characteristic (e.g., spectral amplitude) between the adjacent pixels may be applied. The μ=0.1 graph 640 illustrates the spectral amplitudes of the two adjacent pixels with a μ=0.1 applied by an example weighted regularized approach (e.g., wrNNLS). The first pixel spectral peak corresponds to a small peak component 641 and a large peak component 642. Likewise, the second pixel spectral peak corresponds to a small peak component 648 and a large peak component 649.

[0070] The regularization constant, μ, being equal to 0.1 moves the first spectral peak large component 642 closer to the second spectral peak large component 649. The first pixel spectral peak and the second pixel spectral peak are also associated with small peak components. The two small peak components 641 and 648 are located near their respective large peak components. The peaks became broader with a small regularization coefficient μ= 0.1 , but the first pixel spectral peak small component 641 and the second pixel spectral peak small component 648 are distinct. 10071] A regularized constant is determined to be sufficiently large to smooth the T ₂ spectra so that adjacent pixels of the T ₂ spectral data have similar characteristics. Accordingly, when μ is increased to 0.3, as illustrated in the μ=0.3 graph 650, the first pixel spectral peak 652 and the second pixel spectral peak 659 coincide as broad peaks. Furthermore, the small peak components of the first pixel spectral peak and the second pixel spectral peak are subsumed by the broader first pixel spectral peak 652 and the second pixel spectral peak 659. This is indicative of the spectra near white matter and tissue component becoming one broad peak and having at least more similar spectral characteristics.

[0072] In μ=1.8 graph 660 the first pixel spectral peak 662 and the second pixel spectral peak 669 are approximately the same. Therefore, when μ is increased to 1.8, the MWFs of both pixels were computed equal to 12.7% and 2.1 %, respectively. Thus, as the μ value increases, the sensitivity to measurement noise in the weighted regularized approach decreases. Additionally, the T ₂ spectrum becomes increasingly smoother. The smoothness (fidelity) of spectrum and the sensitivity to noise are compromised to achieve a result indicative of the first adjacent pixel and the second adjacent pixel having more similar characteristics. Therefore, μ is defined to be sufficiently large to smooth the T ₂ spectra so that adjacent pixels of the T ₂ spectral data have characteristics that differ by less than a threshold amount. For example, μ may be defined such that spectral peaks overlap by at least a certain amount (e.g., 85%, 90%, 99%). In another example, μ may be defined to achieve a desired stability of myelin water fraction.

[0073] Figure 7 illustrates the spectral components derived from a T ₂ decay curve. The peaks correspond to components of the T ₂ decay curve. Specifically, the weighted regularized approach with a predetermined μ value, results in spectral components corresponding to myelin water component 710, tissue water component 720, long T ₂ tissue water component 730, and cerebrospinal fluid component 740. The ratios of the components to the summation of the components are termed myelin water fraction (MWF), tissue water fraction (TWF), long T ₂ tissue water fraction (LWF), and cerebral spinal fluid fraction (CSF) respectively.

[0074] Figure 7 illustrates the boundaries of the spectral components. The first boundary 715 separates the myelin water component 710 and the tissue water component 720. The second boundary 725 separates the tissue water component 720 and the long T ₂ tissue water component 730. The third boundary 735 separates the long T ₂ tissue water component 730 and the cerebrospinal fluid component 740. Defining the first boundary 715, the second boundary 725, and the third boundary 735 determines which T ₂ spectral points are used to determine the myelin water component 710, the tissue water component 720, the long tissue water component 730, and the cerebrospinal fluid component 740. Therefore, inaccuracies in the first boundary 715, the second boundary 725, and the third boundary 735 may vary fractions of the myelin water component 710, the tissue water component 720, the long T ₂ tissue water component 730, and/or the cerebrospinal fluid component 740.

[0075] Figure 8 illustrates the variation of average MWF 820 caused by using different numbers of T ₂ sampling points 8 0. Sampling points 810 corresponded to thresholds closest to the specified myelin water T ₂ threshold 40 ms. Therefore, the number of sampling points and a time interval can be used to determine an appropriate boundary between spectral components. When T ₂ is logarithmically sampled, the actual upper T ₂ boundary of myelin water varies with the different number of sampling points 810. In one example, the upper T ₂ boundary change was consistent with average MWF 820. Thus, the variation of average MWF 820 is caused by using the different number of sampling points 810.

[0076] Even with a specified myelin water T ₂ threshold (e.g. 40 ms), the actual threshold may vary with the number of sampling points 810 and become relevant to the robustness of the wrNNLS approach when T ₂ is logarithmically sampled. The largest change of MWF was 2%. A change of 2% in MWF is significant since the average MWF is 12%. [0077] Figure 9 illustrates an image associated with conventional refocusing pulses and an image associated with a refocusing pulse that is three times larger. Conventional myelin water fraction (MWF) image 900 was generated from a conventional refocusing profile. Modified MWF image 950 was generated from a modified refocusing profile. Conventional MWF image 900 and modified MWF image 950 are illustrated in the same gray scale.

[0078] Conventional MWF image 900 and modified MWF image 950 were generated using the same data processing techniques (e.g., wrNNLS approach with μ=1.8 and M = 96) and NLM filter. MWF in the conventional MWF image 900 is substantially lower than that in the modified MWF image 950. The average MWF value in a region of interest (ROI) 910 for the conventional MWF image 900 is 5.7%. The average MWF value in a corresponding ROI 960 of the modified MWF image 950 is 1 1.9%. Therefore, the modified MWF image 950 is more accurate than the conventional MWF image 900.

[0079] Figure 10 illustrates images associated with non-weighted regularized spectral data and images associated with weighted regularized spectral data. Non-weighted parametric maps 1000 correspond to non-weighted MWF 1010, non-weighted TWF 1020, non-weighted LWF 1030, and non-weighted CSF 1040. The non-weighted parametric maps underestimate the non-weighted MWF 1010 near the edge of brain. Larger contamination in the non-weighted LWF map 1030 from the non-weighted CSF map 1040 is also visible because spectral amplitudes corresponding to larger T ₂ value are overweighed in the regularization when using logarithmically spaced sampling.

[0080] The underestimation and contamination increased as the regularization constant increased. Accordingly, image quality of the non-weighted parametric maps 1000 was degraded and became unsuitable for assessing lesions in white matter. Weighted parametric maps 1050 correspond to weighted MWF 1060, weighted TWF 1070, weighted LWF 1080, and weighted CSF 1090. The weighted T ₂ SPARC method generated more uniform parametric maps than the non-weighted parametric maps 1000 as shown in Figure 10. [0081] Figure 11 illustrates images associated a healthy brain and images associated with a brain having lesions. The healthy parametric maps 1100 correspond to healthy map 1105, healthy MWF 1110, healthy TWF 1115, healthy LWF 1120, and healthy CSF 1125. For comparison, diseased parametric maps 1150 are illustrated. The diseased parametric maps 1150 correspond to diseased map 1155, diseased MWF 1160, diseased TWF 1165, diseased LWF 1170, and diseased CSF 1175.

[0082] The difference in contrast in the diseased parametric maps 1150 to the healthy parametric maps 1100 indicates deviations in the corresponding fraction. For example, in the region of interest indicated by arrow 1161 in diseased MWF 1160 has an MWF of 1.9%. In the region of interest indicated by arrow 1162 in diseased MWF 160 has an MWF of 4.5%. In the healthy MWF 1 10, the corresponding regions of interest have an MWF of 12.3%. The differences in contrast of the healthy parametric maps 1100 and the diseased parametric maps 1150 make the regions of interest indicated by the arrows visually identifiable.

[0083] Figure 12 illustrates an example MRI apparatus 1200 configured with a spectral analysis logic 1285. The apparatus 1200 includes a basic field magnet(s) 1210 and a basic field magnet supply 1220. Ideally, the basic field magnets 1210 would produce a uniform B ₀ field. However, in practice, the B ₀ field may not be uniform, and may vary over an object being imaged by the MRI apparatus 1200. MRI apparatus 1200 may include gradient coils 1230 configured to emit gradient magnetic fields like G _s , G _P and G _R . The gradient coils 1230 may be controlled, at least in part, by a gradient coils supply 1240.

[0084] MRI apparatus 1200 may also include an RF antenna 1250 that is configured to generate RF pulses and to receive resulting magnetic resonance signals from an object to which the RF pulses are directed. In some examples, how the pulses are generated and how the resulting MR signals are received may be controlled and thus may be selectively adapted during an MRI procedure. In one example, separate RF transmission and reception coils can be employed. The RF antenna 1250 may be controlled, at least in part, by an RF transmission- reception unit 1260. The gradient coils supply 1240 and the RF transmission- reception unit 1260 may be controlled, at least in part, by a control computer 1270.

[0085] The magnetic resonance signals received from the RF antenna 1250 can be employed to generate an image, and thus may be subject to a transformation process like a two dimensional FFT that generates pixilated image data. The transformation can be performed by an image computer 1280 or other similar processing device. The spectral analysis logic 1285 is configured to perform the methods described herein. The spectral analysis logic 1285 may include a filter logic to de-noise the pixilated image. The image data may then be shown on a display 1290. While an MR apparatus 1200 is illustrated, it is to be appreciated that in some examples the spectral analysis of spectral data may be employed with other imaging apparatus and/or methods (e.g., spectroscopy).

[0086] While Figure 1200 illustrates an example MRI apparatus 1200 that includes various components connected in various ways, it is to be appreciated that other MRI apparatus may include other components connected in other ways. The spectral analysis logic 1285 may be configured with elements of example apparatus described to perform example method described herein.

[0087] In different examples, spectral analysis logic 1285 may be permanently and/or removably attached to an MRI apparatus. While spectral analysis logic 1285 is illustrated as a single logic connected to the image computer 1280, it is to be appreciated that spectral analysis logic 1285 may be distributed between and/or operably connected to other elements of apparatus 1200. Spectral analysis logic 1285 may execute portions of the methods described herein.

[0088] To the extent that the term "or" is employed in the detailed description or claims (e.g., A or B) it is intended to mean "A or B or both". The term "and/or" is used in the same manner, meaning "A or B or both". When the applicants intend to indicate "only A or B but not both" then the term "only A or B but not both" will be employed. Thus, use of the term "or" herein is the inclusive, and not the exclusive use. See, Bryan A. Garner, A Dictionary of Modern Legal Usage 624 (2d. Ed. 1995).

[0089] To the extent that the phrase "one or more of, A, B, and C" is employed herein, (e.g., a data store configured to store one or more of, A, B, and C) it is intended to convey the set of possibilities A, B, C, AB, AC, BC, and/or ABC (e.g., the data store may store only A, only B, only C, A&B, A&C, B&C, and/or A&B&C). It is not intended to require one of A, one of B, and one of C. When the applicants intend to indicate "at least one of A, at least one of B, and at least one of C", then the phrasing "at least one of A, at least one of B, and at least one of C" will be employed.