Search tips
Search criteria 


Logo of wtpaEurope PMCEurope PMC Funders GroupSubmit a Manuscript
Biomed Spectrosc Imaging. Author manuscript; available in PMC 2017 September 20.
Published in final edited form as:
PMCID: PMC5495158

Stroke Onset Time Determination Using MRI Relaxation Times without Non-Ischaemic Reference in A Rat Stroke Model



Objective timing of stroke in emergency departments is expected to improve patient stratification. Magnetic resonance imaging (MRI) relaxations times, T2 and T1ρ, in abnormal diffusion delineated ischaemic tissue were used as proxies of stroke time in a rat model.


Both ‘non-ischaemic reference’-dependent and -independent estimators were generated. Apparent diffusion coefficient (ADC), T2 and T1ρ, were sequentially quantified for up to 6 hours of stroke in rats (n = 8) at 4.7T. The ischaemic lesion was identified as a contiguous collection of voxels with low ADC. T2 and T1ρ in the ischaemic lesion and in the contralateral non-ischaemic brain tissue were determined. Differences in mean MRI relaxation times between ischaemic and non-ischaemic volumes were used to create reference-dependent estimator. For the reference-independent procedure, only the parameters associated with log-logistic fits to the T2 and T1ρ distributions within the ADC-delineated lesions were used for the onset time estimation.


The reference-independent estimators from T2 and T1ρ data provided stroke onset time with precisions of ±32 and ±27 minutes, respectively. The reference-dependent estimators yielded respective precisions of ±47 and ±54 minutes.


A ‘non-ischaemic anatomical reference’-independent estimator for stroke onset time from relaxometric MRI data is shown to yield greater timing precision than previously obtained through reference-dependent procedures.

Keywords: Ischaemic stroke, Magnetic Resonance Imaging, MRI relaxometry, stroke onset time, rat

1. Introduction

The current pharmacological treatments of an ischaemic stroke, including recombinant plasminogen activator, can only be administered within strict time limits [1]. This is due to two key reasons, as (i) irreversible damage to brain parenchyma advances with time [2] and (ii) the benefits of pharmacotherapy are outweighed by adverse effects [3]. Estimates of the number of neural cells and functionality lost per minute of ischaemia [4] stress the earliest possible administration of thrombolytic and/or anti-ischaemic therapy. Thus, stroke duration is seen as a crucial proxy of both brain state and a predictor of functional recovery, and consequently will play a pivotal role in the objective determination of onset time in the clinical settings.

A growing body of evidence indicates magnetic resonance imaging (MRI) protocols provide clinically acceptable timing of stroke onset [58]. The difference between MRI relaxation times T2 and T1ρ in ischaemic brain tissue, as defined by its low apparent diffusion coefficient (ADC), and contra-lateral non-ischaemic tissue has been observed to increase linearly with time, thus, providing the means for estimating the onset of ischaemic stroke [6, 9]. For example, an application of the T2 difference approach for stroke timing has yielded promising results in patients [5], and a similar approach using T2-FLAIR as a proxy of stroke duration [7] has been trialled in patients with variable success [10]. Of these, the latter apparently has an acceptable specificity for clinical exploitation, however, the sensitivity remains low [7] when compared with the use of quantitative T2 relaxation times [5].

The key difficulty in the application of MRI relaxation times in stroke timing clinically is that it is difficult to know the baseline value at each voxel of the ischaemic lesion before the onset of ischaemia, such that the deviation from baseline due to ischaemia cannot be calculated with certainty. For this reason, some authors have sought to select a set of non-ischaemic voxels as a ‘reference’ region, representative of the lesion prior to the onset of ischaemia [5]. However, this approach is confounded by the inherently varying MRI relaxation times across the brain [11], such that the ‘reference’ voxels may not faithfully represent those of the lesion prior to the onset of ischaemia. Given that the distribution of relaxation times in the human brain is broader than the change caused by ischaemia, the challenge of selecting an appropriate set of reference voxels, and the consequences of failing to do so, have presented a major barrier.

Motivated by these factors, we present here a mathematical procedure for ‘non-ischaemic reference’-independent stroke time estimation. Relaxometric MRI data from permanent focal ischaemia in rats was used to assess the capability of the mathematical model for timing of stroke onset.

2. Methods

2.1. Animal model

Male Wistar rats (N=8, approximately 300 g) underwent permanent middle cerebral artery occlusion (MCAO) to induce focal stroke [12]. Animals were euthanized by intravenous injection of saturated KCI at the end of the experiment. Animal procedures were conducted according to European Community Council Directives 86/609/EEC guidelines and approved by the Animal Care and Use Committee of the University of Eastern Finland.

2.2. Image acquisition

MRI data were acquired according to the methods described by Jokivarsi et al. [6]. The trace of the diffusion tensor (ADC= 1/3 ·Trace[D]), T2 and T1ρ MR images from a single 2 mm slice were obtained at 30-60 min intervals for 6 h after induction of focal brain ischaemia using a horizontal Varian UnityInova MRI system operating at 4.7T with a volume coil transmit/quadrature half-volume receive set-up (Rapid Biomedical GmbH, Rimpar, Germany). The ADC was determined from spin echo images with 4 bipolar gradients along each axis with 4 b-values ranging from 0 to 1370 s mm2 (field of view 2.56 x 2.56 cm2, 64 x 128 data points, echo time 55 ms). Fast spin-echo (FSE) readout data (64 x 128 pixels, echo spacing 10 ms, FOV 2.56 x 2.56 cm2) was used for T and T2 MRI. The on-resonance spin-lock T MRI was acquired with a continuous wave T approach. A contrast formation block of AHP-SL-AHP (AHP = adiabatic half passage, SL = spin-lock) segment was added before readout. Adiabatic spin-lock pulses ranging from 8 to 64 ms were used [SL amplitude (B1,SL) 0.4 G, time to repetition (TR) 2.5 s and time to echo (TE) 6 ms]. T2 images were collected using the above FSE readout in front of a preparation block consisting of adiabatic pulses (AHP-TE/4-AFP-TE/2-AFP-TE/4-reverse; five TEs of 10, 20, 40, 60 and 80 ms). The adiabatic refocussing pulses in a single slice mode was used to minimise contributions of non-spin echo coherences [13]. Total acquisition time for ADC, T2 and T1ρ was 24 min.

2.3. Image processing

All image post-processing was performed using in-house designed software written in the MATLAB programming environment (MathWorks, Natick, Massachusetts, USA).

T2 and T1ρ –weighted MR images were filtered with a squared Hamming window function in the phase and frequency encoded dimensions of k-space prior to computing parametric maps. T2 and T1ρ relaxation maps were computed from a voxel-wise monoexponential fit, where the data were first transformed to a logarithmic space then fitted via linear least squares. T2 and T1ρ weighted images were generated by summing all echo and spin-lock time images, respectively.

Lesion detection was performed on reciprocal ADC images where regions of hyperintensity are suggestive of an ischaemic assault. More formally, we have defined a lesion as the largest connected cluster of voxels whose values all exceed one median absolute deviation above the whole brain median value [14]. The output from the detection procedure was a binary mask, where ischaemic voxels were flagged with the value of one and all others set to zero.

2.4. Statistical models for estimating onset time

At present, the state-of-the-art estimators for time from stroke onset are based upon measuring ΔT2 or ΔT1ρ, i.e. the difference between the mean T2 or T1ρ within the lesion and some user-defined reference region, for example, the contra-laterally mirrored lesion [5, 8, 9]. The mean value is typically estimated via the average of the voxel values within a given region. These averages are then used to train a linear regression model that is the main tool through which estimation is performed.

One of the primary disadvantages associated with these estimators is the arbitrary choice of the reference region, which if not chosen appropriately, in terms of the anatomical and/or tissue type match with respect to the ischaemic lesion, can dramatically affect their reliability. Furthermore, their robustness is governed by the ability to accurately estimate the mean T2 or T1ρ values, which is challenging given that the number of ischaemic voxels can vary considerably between strokes and with time in stroked tissue. Also, voxel values can exhibit heavy-tailed empirical distributions for which the mean value is not well defined, making the average estimator potentially unreliable.

Motivated by the disadvantages mentioned above, we have explored the idea of designing robust reference-independent estimators that exploit information from the full empirical distribution of the T2 or T1ρ data in ischaemic voxels only. For improved robustness with large as well as small regions, we underpin the linear regression estimators with a parametric statistical model of the empirical distribution that provides a low-dimensional structure with favourable regularity properties.

2.5. Distributional models for MR parameters

Appropriate parametric statistical models for T2 and T1ρ voxel values within ischaemic lesions were identified by using the MATLAB statistical distribution fitting toolbox. In particular, histograms of T2 and T1ρ within the lesion were computed, then area-normalised to represent probability densities. Individual (specified below) distributions were fitted to the data by maximum likelihood estimation and compared with their empirical distributions. The model that on average scored the best Akaike information criteria (AIC) value was then selected to design the regression method.

2.6. Training and testing of estimators

As with the state-of-the-art, our estimators were chosen to take the form of a multiple linear regression. Here, it is important to note that the assumption of linearity applies only to the coefficients of the estimators and does not necessarily imply that MR parameters vary linearly in time.

To mitigate the possibility of over-fitting we used a leave-one-out cross-validation approach to split the data used to train and to test the estimators. Precisely, our dataset of 8 rats (total of 45 sets of relaxometric data) was iteratively separated into a training set consisting of the scans of 7 rats, and a testing set consisting of the scans of the remaining rat. For each T2 scan in the training set, the parameters associated with the statistical model of ischaemic T2 were estimated by maximum likelihood and then used collectively to train a T2 polynomial regression estimator (this procedure was iterated over all rats, yielding a total of 8 T2 and 8 T1ρ reference-independent estimators). A similar approach was used for signal intensities in T2 and T weighted images. Finally, the performance of the regression estimates was assessed by measuring the root mean-squared-error (RMSE) with respect to the true onset times.

To serve as a comparison, we also computed linear regressions based on ΔT2 and ΔT1ρ predictors, where the reference region was defined as the contra-lateral mirrored lesion as described previously [9].

3. Results

3.1. Lesion detection

ADC images, including lesion mask overlays, from a typical rat brain acquired over multiple time-points are given in Figure 1A. As indicated in the images, the lesion masks form a single contiguous region and their volumes generally increase from time of stroke onset reaching 38±8% of the MRI slice area by 6 hours. ADC values of voxels in the lesion were 49.7±13.4% from those determined in the contralateral non-ischaemic brain parenchyma confirming presence of severe ischaemia [15].

Figure 1
(A) ADC images (green) with lesion masks (pink), (B) T2 and (C) T maps from a typical rat across multiple time points (separated by approximately 1-hour intervals). (D) Distribution of T2 within ADC lesion and (E) distribution of T ...

3.2. Statistical models for MR parameters

Histograms are given for T2 (Figure 1D) and T1ρ (Figure 1E) within the identified lesions for a typical rat (see Figure 1A for ADC, ADC,1B1B for T2 and and1C1C for T images). The histograms show that there was a right-shift in both T2 and T1ρ distributions within the lesion as time from stroke onset increased. It is precisely this information that is exploited in the current ΔT2 and ΔT1ρ estimators [6]. In addition to the value shift, the histograms also reveal that the general spread of the T2 and Tρ distributions increases with time from stroke onset. A similar investigation focusing now on the contra-lateral mirror reference was also performed and found that both T2 and T1ρ distributions remain approximately stationary across time-points (not shown).

As mentioned earlier, in order to identify a low-dimensional representation of the empirical distribution, several statistical distributions were fitted to the data. Namely, the normal, lognormal, log-logistic [16] and scaled Burr distributions [17]. Of these models the 2-parameter log-logistic distribution, i.e.


was found to be the most suitable for both T2 and T1ρ values within the lesion, as determined by obtaining, on average, the best AIC value (Figure 1D and E display the distribution fits in red). Here, the log-logistic parameters (μ, σ) respectively denote the log-median value and distributional shape for a given MR parameter within the lesion.

The log-logistic model was also found to be suitable for describing the distributions of T2 and T1ρ signal intensities within the lesion. However, very little difference was observed between time-points, i.e. little to no drift or spread (Supplementary Information, Figure S1).

3.3. Regression estimators

The anatomical reference-independent regression estimators were designed using (functions of) the parameters associated with the log-logistic fits, i.e. the μ and σ parameters. For the T1ρ-based estimator, a 3-parameter regression (including one constant term) was found to be most accurate:


(R2 = 0.9229, p [double less-than sign] 0.001), with an associated RMSE of ±27mins.

For the T2-based estimator, a 7-parameter estimator based on a third-order polynomial was chosen:


(R2 = 0.8918, p [double less-than sign] 0.001), with an associated RMSE of ±32mins.

Figures 2C and D show the fits for the T2 and T1ρ estimators, respectively, where 95% prediction intervals are given by green curves. These results were notably more accurate than those associated with the 2-parameter regression based on ΔT1ρ (Figure 2A, p < 0.001, pairwise Bartlett test) and ΔT2 (Figure 2B, p = 0.01, pairwise Bartlett test), which only achieved RMSEs of ±54mins and ±47mins respectively. It is worth noting that uncertainties here are larger than previously reported [9]. This is primarily due to different uncertainty quantification approaches. Rogers et al [9] used an analytical uncertainty propagation method that is liberal to approximate precision at a given time point, whereas here we assess precision empirically by cross-validation and report RMSE across all time points. An interval-wise breakdown of the results for all estimators is given in Table 1.

Figure 2
Estimated time vs. actual time of onset for several regression estimators. (A) Reference-dependent (B) Reference-dependent (C) Reference independent ΔT1ρ multiple linear regression, (D) Reference-independent ΔT2 multiple polynomial ...
Table 1
Interval-wise breakdown of estimated vs. actual onset times. Values are displayed as mean ± SD. 95% confidence intervals for estimated means were computed from mean ± t0.95,N [center dot]SEM, where SEM is standard error of mean, and t0.95, ...

In addition to measuring the precision of reference-dependent estimators, we investigated their robustness through consideration of randomly-shifted reference regions (Supplementary Information, Figure S2 and S3). The results showed that overall precision is sensitive to reference placement, with RMSEs often larger than those reported for the mirrored contra-lateral reference.

Estimators based on T2 and T1ρ signal intensities, averaged over their respective echo and spin-lock times, were also considered. However, they performed poorly because the statistical distributions lacked any significant time-dependence (Supplementary Information, Figure S1). A reference-independent estimator that combined both T2 and T1ρ log-logistic parameters was also constructed. However, our results indicated that these were highly correlated, and therefore no significant improvement in precision was observed.

4. Discussion

The current study has established a stroke onset timing procedure using quantitative relaxometric MRI data only from the ischaemic lesion delineated by low ADC. The use of low ADC for selection of tissue for relaxometric data assures that only ischaemic volumes [15] are used for timing. Low ADC eliminates non-ischaemic regions from timing analyses, because reperfusion in the early hours of ischaemia normalises ADC [18, 19]. It should be noted that quantitative ADC alone provides no time information [20], neither does the rate at which ADC lesion expands [14] in the early hours of ischaemia in permanent MCAO. This ‘anatomical reference-independent’ procedure performs better in timing than the one using MRI data from a non-ischaemic reference region. As the quantitative T2 (which is commonly available in standard clinical MRI systems) provides excellent precision for stroke time estimates in the preclinical model, the presented procedure is thought to bear significant potential for stroke patient stratification.

The current 1H MRI-based timing estimation procedures rely on the differences in average values of distributions of relaxation times between ischaemic tissue and non-ischaemic reference, i.e. ΔT2 and ΔT1ρ. A substantial body of evidence shows that the time dependence of both ΔT2 and ΔT1ρ in the rat brain can be adequately described by a linear function [6, 8, 9] and ΔT2 similarly in human brain [5]. The ΔT2–based estimators have been reported to classify stroke onset within 3 hours with a precision of 0.850 in rats [21] and 0.757 in humans [5] indicative of clinical utility. However, the precision of timing by these reference-dependent estimators strictly depends on the selection of an anatomically identical non-ischaemic reference region for the ischaemic lesion. This is a non-trivial task even for the simple rat brain (see Supplementary Information Figure S3). In humans, reference selection will become even more problematic due to inherent relaxation time differences in GM and WM [11], variation of T2 in WM due to microstructure [22], as well as the lack of ‘perfect anatomical symmetry’ of hemispheres to choose a suitable ‘mirror reference’ over the midline. Furthermore, stroke lesions in humans vary in location and size.

The stroke time estimators introduced in this paper use only the T2 or T1ρ data from low ADC-defined ischaemic brain tissue. Our data indicate that information embedded within the distributions of T2 and T1ρ can be exploited to build multiple linear regression estimators independent of MRI data from a non-ischaemic reference, provided relaxation data from multiple time points are available. In other words, the time dependency of the parameters in regression estimators needs to be determined similar to the ΔT2 and ΔT1ρ time dependencies used in the recent estimators [5, 6]. The pathophysiological basis underpinning the estimators is likely to be that the volume of brain tissue with abnormally elevated MRI relaxation times increases linearly with time [8]. This is due to the effects of cytotoxic oedema on both T [23] and T2 relaxation times with macromolecular breakdown [14] followed by vasogenic oedema on hydrodynamics within ischaemic tissue [24]. In addition to cytotoxic and vasogenic oedema, altered chemical exchange processes may influence T relaxation in ischaemia [25]. The spread of effects of ischaemia in the brain tissue in the early hours of the insult has been proposed to be the critical factor underpinning the time-dependent increase in T2 relaxation time due to cytotoxic oedema [14]. We speculate that this spread is captured by the log-logistic fits of the relaxation time distributions to provide accurate stroke time estimation.

Regarding the utility of the current time estimator in stroke imaging, several points should be considered. Firstly, while both T and T2 MRI are in use in general radiology, T is not part of brain MRI protocol in clinic due to specific absorption rate of energy concerns. Therefore, the estimator based on T2 data may bear clinical potentials. Secondly, it should be stressed that the time estimators are generated from the quantitative MRI data obtained during early hours of permanent ischaemia. Thus, the current estimators provide timing only from the tissue where residual perfusion is below the ischaemia threshold [26]. A previous study showed that reperfusion of ischaemic tissue, but not of that with compromised blood flow and normal ADC, leads to prolongation of tissue T2 [27]. Against this, reperfusion of previously ischaemic tissue during MRI data acquisition would invalidate the established ΔT2 vs stroke onset time relationships [6, 8] and potentially also the T2 relaxation time distribution characteristics used by the currently log-logistic estimator. In clinical settings reperfusion, either spontaneous or pharmacologically induced, is most likely to take place when the cloth is small and/or the occlusions involves small vessels; instead, reperfusion of proximal branches of cerebral arteries is uncommon [28]. Thirdly, in the current data set the number of voxels in ADC lesions varied from 180 to 1,330. A human stroke lesion of 5 cm3 in volume would have >1,600 voxels in a T2 maps acquired at routine spatial resolution of 3 mm3, i.e. even small strokes would have enough T2 data for generation of reference-independent estimators. Finally, a matter to be considered is that the ischaemic rat tissue is almost entirely gray matter and our time estimator apparently holds only for this tissue type.

5. Conclusions

The current study indicates that relaxometric MRI data may provide stroke time estimators without the need for non-ischaemic reference data. Critical issues in the clinical implementation of such protocols include both scan time and data post-processing. A recent study demonstrated a comprehensive stroke MRI protocol can be acquired in six minutes [29]. Investigational MRI scan procedures, such as magnetic resonance finger printing [30], have been shown to further shorten scan times, thus paving the way for expedited clinical brain examination. As far as the post-processing is concerned, the current standard computers are capable of analysing MRI data sets in seconds. Overall, quantitative MRI-based stroke timing taking minutes is likely to become reality for clinical exploitation in the near future.

Supplementary Material



The expert technical assistance by Ms Maarit Pulkkinen is greatly acknowledged. MP holds a Marie Curie Intra-European Fellowship for Career Development. BLM is a recipient of EPSRC PhD studentship. MJK is funded by Elizabeth Blackwell Institute Early Career Fellowship and by the Wellcome Trust international strategic support fund ISSF2 : 105612/Z/14/Z. KTJ and OHJG are funded by Academy of Finland, Finnish Technology Funding (TEKES) and Biocenter Finland. RAK and MP are funded by the Dunhill Medical Trust.


1. Hacke W, Donnan G, Fieschi C, Kaste M, von Kummer R, Broderick JP, Brott T, Frankel M, Grotta JC, Haley EC, Jr, Kwiatkowski T, et al. Association of outcome with early stroke treatment: pooled analysis of ATLANTIS, ECASS, and NINDS rt-PA stroke trials. Lancet. 2004;363:768–74. [PubMed]
2. Saver JL. Time is brain--quantified. Stroke. 2006;37:263–266. [PubMed]
3. Emberson J, Lees KR, Lyden P, Blackwell L, Albers G, Bluhmki E, Brott T, Cohen G, Davis S, Donnan G, Grotta J, et al. Effect of treatment delay, age, and stroke severity on the effects of intravenous thrombolysis with alteplase for acute ischaemic stroke: a meta-analysis of individual patient data from randomised trials. Lancet. 2014;384:1929–1935. [PMC free article] [PubMed]
4. Meretoja A, Keshtkaran M, Saver JL, Tatlisumak T, Parsons MW, Kaste M, Davis SM, Donnan GA, Churilov L. Stroke thrombolysis: save a minute, save a day. Stroke. 2014;45:1053–1058. [PubMed]
5. Siemonsen S, Mouridsen K, Holst B, Ries T, Finsterbusch J, Thomalla G, Ostergaard L, Fiehler J. Quantitative t2 values predict time from symptom onset in acute stroke patients. Stroke. 2009;40:1612–1626. [PubMed]
6. Jokivarsi KT, Hiltunen Y, Grohn H, Tuunanen PI, Grohn OHJ, Kauppinen RA. Estimation of the onset time of cerebral ischemia using T1r and T2 MRI in rats. Stroke. 2010;41:2335–2340. [PubMed]
7. Thomalla G, Cheng B, Ebinger M, Hao Q, Tourdias T, Wu O, Kim JS, Breuer L, Singer OC, Warach S, Christensen S, et al. DWI-FLAIR mismatch for the identification of patients with acute ischaemic stroke within 4.5 h of symptom onset (PRE-FLAIR): a multicentre observational study. Lancet Neurol. 2011;10:978–986. [PubMed]
8. McGarry BL, Rogers HJ, Knight MJ, Jokivarsi KT, Sierra A, Grohn OH, Kauppinen RA. Stroke onset time estimation from multispectral quantitative magnetic resonance imaging in a rat model of focal permanent cerebral ischemia. Int J Stroke. 2016:677–682. [PubMed]
9. Rogers HJ, McGarry BL, Knight MJ, Jokivarsi KT, Grohn OHJ, Kauppinen RA. Timing the ischaemic stroke by 1H MRI: Improved accuracy using absolute relaxation times over signal intensities. NeuroReport. 2014;25:1180–1185. [PubMed]
10. Cheng B, Brinkmann M, Forkert ND, Treszl A, Ebinger M, Kohrmann M, Wu O, Kang DW, Liebeskind DS, Tourdias T, Singer OC, et al. Quantitative measurements of relative fluid-attenuated inversion recovery (FLAIR) signal intensities in acute stroke for the prediction of time from symptom onset. J Cereb Blood Flow Metab. 2013;33:76–84. [PMC free article] [PubMed]
11. Wansapura JP, Holland SK, Dunn RS, Ball WS., Jr NMR relaxation times in the human brain at 3.0 tesla. J Magn Reson Imaging. 1999;9:531–538. [PubMed]
12. Longa EZ, Weinstein PR, Carlson S, Cummins R. Reversible middle cerebral artery occlusion without craniectomy in rats. Stroke. 1989;20:84–91. [PubMed]
13. Mitsumori F, Watanabe H, Takaya N, Garwood M, Auerbach EJ, Michaeli S, Mangia S. Toward understanding transverse relaxation in human brain through its field dependence. Magn Reson Med. 2012;68:947–953. [PMC free article] [PubMed]
14. Knight MJ, McGarry BL, Rogers HJ, Jokivarsi KT, Grohn OH, Kauppinen RA. A spatiotemporal theory for MRI T2 relaxation time and apparent diffusion coefficient in the brain during acute ischaemia: Application and validation in a rat acute stroke model. J Cereb Blood Flow Metab. 2016;36:1232–1243. [PMC free article] [PubMed]
15. Hoehn-Berlage M, Norris DG, Kohno K, Mies G, Leibfritz D, Hossmann KA. Evolution of regional changes in apparent diffusion coefficient during focal ischemia of rat brain: the relationship of quantitative diffusion NMR imaging to reduction in cerebral blood flow and metabolic disturbances. J Cereb Blood Flow Metab. 1995;15:1002–1011. [PubMed]
16. Forbes C, Evans M, Hastings N, Peacock B. Staistical Distributions. 4th ed. John Wiley & Sons; 2011.
17. Burr IW. Cumulative frequency functions. Ann Math Sci. 1942;13:215–232.
18. Li F, Han SS, Tatlisumak T, Liu KF, Garcia JH, Sotak CH, Fisher M. Reversal of acute apparent diffusion coefficient abnormalities and delayed neuronal death following transient focal cerebral ischemia in rats. Ann Neurol. 1999;46:333–342. [PubMed]
19. Kettunen MI, Gröhn OHJ, Lukkarinen JA, Vainio P, Silvennoinen MJ, Kauppinen RA. Interrelations of T1 and diffusion of water in acute cerebral ischemia of the rat. Magn Reson Med. 2000;44:833–839. [PubMed]
20. Knight RA, Dereski MO, Helpern JA, Ordidge RJ, Chopp M. Magnetic resonance imaging assessment of evolving focal cerebral ischemia. Comparison with histopathology in rats. Stroke. 1994;25:1252–1261. [PubMed]
21. McGarry BL, Rogers HJ, Knight MJ, Jokivarsi K, Grohn H, Kauppinen RA. Determining Stroke Onset Time Using Quantitative MRI: High Accuracy, Sensitivity and Specificity Obtained from Magnetic Resonance Relaxation Times. Cerebrovasc Dis Extra. 2016;6:60–65.
22. Knight MJ, Wood B, Coulthard E, Kauppinen RA. Anisotropy of spin echo T2 relaxation in the human brain in vivo. Biomed Spectr Imag. 2015;4:299–310.
23. Kettunen MI, Grohn OH, Silvennoinen MJ, Penttonen M, Kauppinen RA. Effects of intracellular pH, blood, and tissue oxygen tension on T1rho relaxation in rat brain. Magn Reson Med. 2002;48:470–477. [PubMed]
24. Siesjö BK. Mechanisms of ischemic brain damage. Crit Care Med. 1988;16:954–963. [PubMed]
25. Jin T, Autio J, Obata T, Kim SG. Spin-locking versus chemical exchange saturation transfer MRI for investigating chemical exchange process between water and labile metabolite protons. Magn Reson Med. 2011;65:1448–1460. [PMC free article] [PubMed]
26. Busza AL, Allen KL, King MD, van Bruggen N, Williams SR, Gadian DG. Diffusion-weighted imaging studies of cerebral ischemia in gerbils. Potential relevance to energy failure. Stroke. 1992;23:1602–1612. [PubMed]
27. Kettunen MI, Gröhn OHJ, Silvennoinen MJ, Penttonen M, Kauppinen RA. Quantitative assessment of the balance between oxygen delivery and consumption in the rat brain after transient ischemia with T2-BOLD magnetic resonance imaging. J Ceber Blood Flow Metab. 2002;22:262–270. [PubMed]
28. Miller JB, Merck LH, Wira CR, Meurer WJ, Schrock JW, Nomura JT, Siket MS, Madsen TE, Wright DW, Panagos PD, Lewandowski C. The Advanced Reperfusion Era: Implications for Emergency Systems of Ischemic Stroke Care. Ann Emerg Med. 2016 in press. [PubMed]
29. Nael K, Khan R, Choudhary G, Meshksar A, Villablanca P, Tay J, Drake K, Coull BM, Kidwell CS. Six-minute magnetic resonance imaging protocol for evaluation of acute ischemic stroke: pushing the boundaries. Stroke. 2014;45:1985–1991. [PubMed]
30. Ma D, Gulani V, Seiberlich N, Liu K, Sunshine JL, Duerk JL, Griswold MA. Magnetic resonance fingerprinting. Nature. 2013;495:187–192. [PMC free article] [PubMed]