Frequency of Cardiovascular Involvement in Familial Amyloidotic Polyneuropathy in Brazilian Patients

Background:Familial amyloidotic polyneuropathy (FAP) is a rare disease diagnosed in Brazil and worldwide. The frequency of cardiovascular involvement in Brazilian FAP patients is unknown.Objective:Detect the frequency of cardiovascular involvement and correlate the cardiovascular findings with the modified polyneuropathy disability (PND) score.Methods:In a national reference center, 51 patients were evaluated with clinical examination, electrocardiography (ECG), echocardiography (ECHO), and 24-hour Holter. Patients were classified according to the modified PND score and divided into groups: PND 0, PND I, PND II, and PND > II (which included PND IIIa, IIIb, and IV). We chose the classification tree as the statistical method to analyze the association between findings in cardiac tests with the neurological classification (PND).Results:ECG abnormalities were present in almost 2/3 of the FAP patients, whereas ECHO abnormalities occurred in around 1/3 of them. All patients with abnormal ECHO also had abnormal ECG, but the opposite did not apply. The classification tree identified ECG and ECHO as relevant variables (p < 0.001 and p = 0.08, respectively). The probability of a patient to be allocated to the PND 0 group when having a normal ECG was over 80%. When both ECG and ECHO were abnormal, this probability was null.Conclusions:Brazilian patients with FAP have frequent ECG abnormalities. ECG is an appropriate test to discriminate asymptomatic carriers of the mutation from those who develop the disease, whereas ECHO contributes to this discrimination.

Geometric distortion correction in EPI at ultra high field

Magdeburg, Univ., Fak. für Naturwiss., Diss., 2012

3-T direct MR arthrography of the wrist: value of finger trap distraction to assess intrinsic ligament and triangular fibrocartilage complex tears.

PURPOSE: To determine the value of applying finger trap distraction during direct MR arthrography of the wrist to assess intrinsic ligament and triangular fibrocartilage complex (TFCC) tears. MATERIALS AND METHODS: Twenty consecutive patients were prospectively investigated by three-compartment wrist MR arthrography. Imaging was performed with 3-T scanners using a three-dimensional isotropic (0.4 mm) T1-weighted gradient-recalled echo sequence, with and without finger trap distraction (4 kg). In a blind and independent fashion, two musculoskeletal radiologists measured the width of the scapholunate (SL), lunotriquetral (LT) and ulna-TFC (UTFC) joint spaces. They evaluated the amount of contrast medium within these spaces using a four-point scale, and assessed SL, LT and TFCC tears, as well as the disruption of Gilula's carpal arcs. RESULTS: With finger trap distraction, both readers found a significant increase in width of the SL space (mean Δ = +0.1mm, p ≤ 0.040), and noticed more contrast medium therein (p ≤ 0.035). In contrast, the differences in width of the LT (mean Δ = +0.1 mm, p ≥ 0.057) and UTFC (mean Δ = 0mm, p ≥ 0.728) spaces, as well as the amount of contrast material within these spaces were not statistically significant (p = 0.607 and ≥ 0.157, respectively). Both readers detected more SL (Δ = +1, p = 0.157) and LT (Δ = +2, p = 0.223) tears, although statistical significance was not reached, and Gilula's carpal arcs were more frequently disrupted during finger trap distraction (Δ = +5, p = 0.025). CONCLUSION: The application of finger trap distraction during direct wrist MR arthrography may enhance both detection and characterisation of SL and LT ligament tears by widening the SL space and increasing the amount of contrast within the SL and LT joint spaces.

Zeichentheorie und BWL : Zur Anwendung semiotischer Erkentnisse im wirtschaftlichen Diskurs

Semiotics is hardly known in German business literature and management practice, despite the fact that its methodological approaches to discourse analysis have already been broadly absorbed by management theory in France and the United States. The present contribution points out why the general echo of semiotic theory and its numerous applications to business administration has remained limited so far, especially when compared with its potential for describing and explaining management problems. It is important, then, to show what the object of semiotic research is, independent of any semiotic orientation or school. What tools have been developed so far, and which discourses in business administration were chosen to apply these tools? The problems limiting a broader use of semiotic instruments in business administration are explained in detail, and the research perspectives are illustrated. Die in der deutschsprachigen Betriebswirtschaftslehre (BWL) noch wenig bekannte Semiotik und die in ihrem Rahmen entwickelten methodologischen Ansätze zur Analyse von Diskursen haben insbesondere in Frankreich und in den USA schon eine weite Verbreitung gefunden. Das Echo, auf welches die Semiotik trifft, bleibt im Vergleich zum Beitrag, den sie hinsichtlich der Beschreibung und Erklärung betriebswirtschaftlicher Tatbestände leisten könnte, noch sehr beschränkt - und dies, obwohl ihre grundsätzliche wissenschaftliche Leistungsfähigkeit inzwischen an Hand zahlreicher Beispiele aus dem betriebswirtschaftlichen Bereich belegt werden konnte. Es ist deshalb wichtig und interessant zu zeigen, was der Forschungsgegenstand der Semiotik ist - und zwar unabhängig von den verschiedenen Schulen der Semiotik. Welche Instrumente hat sie entwickelt? Auf welche betriebswirtschaftlichen Diskurse sind diese Instrumente bereits angewandt worden? Zum Schluss werden die Probleme dargestellt, die einen breiteren Einsatz semiotischer Instrumente in der Betriebswirtschaftslehre behindern und dabei auch die verschiedenen Forschungsrichtungen erläutert.

Alveolar echinococcosis of the liver: Diffusion-weighted MRI findings and potential role in lesion characterisation

PURPOSE: To report the diffusion-weighted MRI findings in alveolar echinococcosis (AE) of the liver and evaluate the potential role of apparent diffusion coefficients (ADCs) in the characterisation of lesions. MATERIALS AND METHODS: We retrospectively included 22 patients with 63 AE liver lesions (≥1cm), examined with 3-T liver MRI, including a free-breathing diffusion-weighted single-shot echo-planar imaging sequence (b-values=50, 300 and 600s/mm(2)). Two radiologists jointly assessed the following lesion features: size, location, presence of cystic and/or solid components (according to Kodama's classification system), relative contrast enhancement, and calcifications (on CT). The ADCtotal, ADCmin and ADCmax were measured in each lesion and the surrounding liver parenchyma. RESULTS: Three type 1, 19 type 2, 17 type 3, three type 4 and 21 type 5 lesions were identified. The mean (±SD) ADCtotal, ADCmin and ADCmax for all lesions were 1.73±0.50, 0.76±0.38 and 2.63±0.76×10(-3)mm(2)/s, respectively. The mean ADCtotal for type 1, type 2, type 3, type 4 and type 5 lesions were 1.97±1.01, 1.76±0.53, 1.73±0.41, 1.15±0.42 and 1.76±0.44×10(-3)mm(2)/s, respectively. No significant differences were found between the five lesion types, except for type 4 (p=0.0363). There was a significant correlation between the presence of a solid component and low ADCmin (r=0.39, p=0.0016), whereas an inverse correlation was found between the relative contrast enhancement and ADCtotal (r=-0.34, p=0.0072). CONCLUSION: The ADCs of AE lesions are relatively low compared to other cystic liver lesions, which may help in the differential diagnosis. Although ADCs are of little use to distinguish between the five lesion types, their low value reflects the underlying solid component.

Coronary magnetic resonance angiography for assessment of the stent lumen: a phantom study.

PURPOSE: To evaluate the feasibility of visualizing the stent lumen using coronary magnetic resonance angiography in vitro. MATERIAL AND METHODS: Nineteen different coronary stents were implanted in plastic tubes with an inner diameter of 3 mm. The tubes were positioned in a plastic container filled with gel and included in a closed flow circuit (constant flow 18 cm/sec). The magnetic resonance images were obtained with a dual inversion fast spin-echo sequence. For intraluminal stent imaging, subtraction images were calculated from scans with and without flow. Subsequently, intraluminal signal properties were objectively assessed and compared. RESULTS: As a function of the stent type, various degrees of in-stent signal attenuation were observed. Tantalum stents demonstrated minimal intraluminal signal attenuation. For nitinol stents, the stent lumen could be identified, but the intraluminal signal was markedly reduced. Steel stents resulted in the most pronounced intraluminal signal voids. CONCLUSIONS: With the present technique, radiofrequency penetration into the stents is strongly influenced by the stent material. Thesefindings may have important implicationsforfuture stent design and stent imaging strategies.

Determinants of renal tissue oxygenation as measured with BOLD-MRI in chronic kidney disease and hypertension in humans.

Experimentally renal tissue hypoxia appears to play an important role in the pathogenesis of chronic kidney disease (CKD) and arterial hypertension (AHT). In this study we measured renal tissue oxygenation and its determinants in humans using blood oxygenation level-dependent magnetic resonance imaging (BOLD-MRI) under standardized hydration conditions. Four coronal slices were selected, and a multi gradient echo sequence was used to acquire T2* weighted images. The mean cortical and medullary R2* values ( = 1/T2*) were calculated before and after administration of IV furosemide, a low R2* indicating a high tissue oxygenation. We studied 195 subjects (95 CKD, 58 treated AHT, and 42 healthy controls). Mean cortical R2 and medullary R2* were not significantly different between the groups at baseline. In stimulated conditions (furosemide injection), the decrease in R2* was significantly blunted in patients with CKD and AHT. In multivariate linear regression analyses, neither cortical nor medullary R2* were associated with eGFR or blood pressure, but cortical R2* correlated positively with male gender, blood glucose and uric acid levels. In conclusion, our data show that kidney oxygenation is tightly regulated in CKD and hypertensive patients at rest. However, the metabolic response to acute changes in sodium transport is altered in CKD and in AHT, despite preserved renal function in the latter group. This suggests the presence of early renal metabolic alterations in hypertension. The correlations between cortical R2* values, male gender, glycemia and uric acid levels suggest that these factors interfere with the regulation of renal tissue oxygenation.

Metallic renal artery MR imaging stent: artifact-free lumen visualization with projection and standard renal MR angiography.

A cardiac-triggered free-breathing three-dimensional (3D) balanced fast field-echo projection renal magnetic resonance (MR) angiographic sequence was investigated for in-stent lumen visualization of a dedicated metallic renal artery stent. Fourteen prototype stents were deployed in the renal arteries of six pigs (in two pigs, three stents were deployed). Projection renal MR angiography was compared with standard contrast material-enhanced 3D breath-hold MR angiography. Artifact-free in-stent lumen visualization was achieved with both projection MR angiography and contrast-enhanced MR angiography. These promising results warrant further studies for visualization of in-stent restenosis.

Diffusion-weighted MR imaging of the spine at 3-T: feasibility, optimization of b-value and utility to differentiate benign from pathological vertebral compression fractures

Purpose: To evaluate the feasibility, determine the optimal b-value, and assess the utility of 3-T diffusion-weighted MR imaging (DWI) of the spine in differentiating benign from pathologic vertebral compression fractures.Methods and Materials: Twenty patients with 38 vertebral compression fractures (24 benign, 14 pathologic) and 20 controls (total: 23 men, 17 women, mean age 56.2years) were included from December 2010 to May 2011 in this IRB-approved prospective study. MR imaging of the spine was performed on a 3-T unit with T1-w, fat-suppressed T2-w, gadolinium-enhanced fat-suppressed T1-w and zoomed-EPI (2D RF excitation pulse combined with reduced field-of-view single-shot echo-planar readout) diffusion-w (b-values: 0, 300, 500 and 700s/mm2) sequences. Two radiologists independently assessed zoomed-EPI image quality in random order using a 4-point scale: 1=excellent to 4=poor. They subsequently measured apparent diffusion coefficients (ADCs) in normal vertebral bodies and compression fractures, in consensus.Results: Lower b-values correlated with better image quality scores, with significant differences between b=300 (mean±SD=2.6±0.8), b=500 (3.0±0.7) and b=700 (3.6±0.6) (all p<0.001). Mean ADCs of normal vertebral bodies (n=162) were 0.23, 0.17 and 0.11×10-3mm2/s with b=300, 500 and 700s/mm2, respectively. In contrast, mean ADCs were 0.89, 0.70 and 0.59×10-3mm2/s for benign vertebral compression fractures and 0.79, 0.66 and 0.51×10-3mm2/s for pathologic fractures with b=300, 500 and 700s/mm2, respectively. No significant difference was found between ADCs of benign and pathologic fractures.Conclusion: 3-T DWI of the spine is feasible and lower b-values (300s/mm2) are recommended. However, our preliminary results show no advantage of DWI in differentiating benign from pathologic vertebral compression fractures.

A Modulated Closed Form solution for Quantitative Susceptibility Mapping - A thorough evaluation and comparison to iterative methods based on edge prior knowledge.

The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one over-determined multiple orientation method (COSMOS) and anewly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65mm isotropic) brain data acquired at 7T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corrMCF=0.95, r(2)MCF =0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The l-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures.

Choledocholithiasis: repetitive thick-slab single-shot projection magnetic resonance cholangiopancreaticography versus endoscopic ultrasonography.

This prospective study compares repetitive thick-slab single-shot projection magnetic resonance cholangiopancreatography (MRCP) with endoscopic ultrasonography (EUS) for the detection of choledocholithiasis. Fifty-seven consecutive patients (36 women, mean age 61) referred for suspected choledocholithiasis underwent MRCP, followed by EUS. Each procedure was performed by different operators blinded to the results of the other investigation. MR technique included a turbo spin-echo T2-weighted axial sequence with selective fat saturation (SPIR/TSE, TE=70 ms, TR=1,600 ms), followed by coronal dynamic MRCP. The same thick-slab slice was sequentially acquired 12 times as breath-hold single-shot projection imaging (SSh, TE=900 ms, TE=8,000 ms) centred on the common bile duct (CBD). Two experienced radiologists independently and blindly evaluated MR images for the detection of CBD stones. Their inter-observer agreement kappa was determined. Secondly, the two observers read MR images in consensus again. CBD stones were demonstrated in 18 out of 57 patients (31.6 %) and confirmed by endoscopic retrograde cholangiography (ERCP, n=17) or intraoperative cholangiography (n=1). Clinical follow-up served as the "gold standard" in patients with negative results without following invasive procedure (n=28). Sensitivity, specificity, accuracy, positive and negative predictive value for MRCP resulting from consensus reading were 94.9%, 94.4%, 94.7%, 97.4% and 89.5%, respectively. Corresponding values of EUS were 97.4%, 94.4%, 96.5%, 97.4% and 94.4%. Inter-observer agreement kappa was 0.81. Repetitive thick-slab single-shot projection MRCP is an accurate non-invasive imaging modality for suspected choledocholithiasis and should be increasingly used to select those patients who require a subsequent therapeutic procedure, namely ERCP.

Intravoxel incoherent motion perfusion imaging in acute stroke: initial clinical experience.

INTRODUCTION: Intravoxel incoherent motion (IVIM) imaging is an MRI perfusion technique that uses a diffusion-weighted sequence with multiple b values and a bi-compartmental signal model to measure the so-called pseudo-diffusion of blood caused by its passage through the microvascular network. The goal of the current study was to assess the feasibility of IVIM perfusion fraction imaging in patients with acute stroke. METHODS: Images were collected in 17 patients with acute stroke. Exclusion criteria were onset of symptoms to imaging >5 days, hemorrhagic transformation, infratentorial lesions, small lesions <0.5 cm in minimal diameter and hemodynamic instability. IVIM imaging was performed at 3 T, using a standard spin-echo Stejskal-Tanner pulsed gradients diffusion-weighted sequence, using 16 b values from 0 to 900 s/mm(2). Image quality was assessed by two radiologists, and quantitative analysis was performed in regions of interest placed in the stroke area, defined by thresholding the apparent diffusion coefficient maps, as well as in the contralateral region. RESULTS: IVIM perfusion fraction maps showed an area of decreased perfusion fraction f in the region of decreased apparent diffusion coefficient. Quantitative analysis showed a statistically significant decrease in both IVIM perfusion fraction f (0.026 ± 0.019 vs. 0.056 ± 0.025, p = 2.2 · 10(-6)) and diffusion coefficient D compared with the contralateral side (3.9 ± 0.79 · 10(-4) vs. 7.5 ± 0.86 · 10(-4) mm(2)/s, p = 1.3 · 10(-20)). CONCLUSION: IVIM perfusion fraction imaging is feasible in acute stroke. IVIM perfusion fraction is significantly reduced in the visible infarct. Further studies should evaluate the potential for IVIM to predict clinical outcome and treatment response.

Functional mapping of the human visual cortex with intravoxel incoherent motion MRI.

Functional imaging with intravoxel incoherent motion (IVIM) magnetic resonance imaging (MRI) is demonstrated. Images were acquired at 3 Tesla using a standard Stejskal-Tanner diffusion-weighted echo-planar imaging sequence with multiple b-values. Cerebro-spinal fluid signal, which is highly incoherent, was suppressed with an inversion recovery preparation pulse. IVIM microvascular perfusion parameters were calculated according to a two-compartment (vascular and non-vascular) diffusion model. The results obtained in 8 healthy human volunteers during visual stimulation are presented. The IVIM blood flow related parameter fD* increased 170% during stimulation in the visual cortex, and 70% in the underlying white matter.

Quantitative Measurement of Brain Perfusion with Intravoxel Incoherent Motion MR Imaging.

Purpose: To evaluate the sensitivity of the perfusion parameters derived from Intravoxel Incoherent Motion (IVIM) MR imaging to hypercapnia-induced vasodilatation and hyperoxygenation-induced vasoconstriction in the human brain. Materials and Methods: This study was approved by the local ethics committee and informed consent was obtained from all participants. Images were acquired with a standard pulsed-gradient spin-echo sequence (Stejskal-Tanner) in a clinical 3-T system by using 16 b values ranging from 0 to 900 sec/mm(2). Seven healthy volunteers were examined while they inhaled four different gas mixtures known to modify brain perfusion (pure oxygen, ambient air, 5% CO(2) in ambient air, and 8% CO(2) in ambient air). Diffusion coefficient (D), pseudodiffusion coefficient (D*), perfusion fraction (f), and blood flow-related parameter (fD*) maps were calculated on the basis of the IVIM biexponential model, and the parametric maps were compared among the four different gas mixtures. Paired, one-tailed Student t tests were performed to assess for statistically significant differences. Results: Signal decay curves were biexponential in the brain parenchyma of all volunteers. When compared with inhaled ambient air, the IVIM perfusion parameters D*, f, and fD* increased as the concentration of inhaled CO(2) was increased (for the entire brain, P = .01 for f, D*, and fD* for CO(2) 5%; P = .02 for f, and P = .01 for D* and fD* for CO(2) 8%), and a trend toward a reduction was observed when participants inhaled pure oxygen (although P > .05). D remained globally stable. Conclusion: The IVIM perfusion parameters were reactive to hyperoxygenation-induced vasoconstriction and hypercapnia-induced vasodilatation. Accordingly, IVIM imaging was found to be a valid and promising method to quantify brain perfusion in humans. © RSNA, 2012.

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.