Friday, December 16, 2011

Forward from http://www.yalescientific.org/2010/12/fourier-transform-natures-way-of-analyzing-data/

Fourier Transform: Nature’s Way of Analyzing Data

Described as “nature’s way of analyzing data” by Yale professor Ronald Coifman, the Fourier Transform is arguably the most powerful analytical tool in modern mathematics. Professor Peter Moore, a Yale structural biologist and professor of biophysics, agrees. “To form an image on your retina, the lens in your eye performs Fourier transformations on the light that enters it,” he explains. This tool is truly ubiquitous in nature, as our eyes and ears have subconsciously performed the Fourier transform to interpret sound and light waves for millions of years. Hence it was only a matter of time until the human intellect caught up to our internal processing systems and was able to functionally describe this process. After years of research, French Baron Jean-Baptiste-Joseph Fourier uncovered this powerful tool in the early 1800s, naming it the Fourier transform.
Fourier, a French military scientist, became interested in heat transfer in the late 1790s. In fact, many of his guests often complained that he kept his home uncomfortably warm. During Napoleon’s expansion campaigns, Fourier served on the Institute of Egypt’s scientific body in 1800, and after the French left Egypt, he set his efforts on repairing France from the devastation of the 1789 French Revolution. During this time, his obsession with heat transfer drove him to derive an equation describing the conduction of heat in solid bodies. Within seven years, he invented the Fourier transform to solve this equation.
The question itself was complicated; Fourier wanted to solve his equation to describe the flow of heat around an iron ring that attaches a ship’s anchor to its chain. He proposed that the irregular distribution of temperature could be described by the frequencies of many component sinusoidal waves around the ring. His premise was that the maximum temperature and position of the harmonics of these sinusoidal components could be derived via the Fourier transform of the originally irregular distribution of temperatures.
With today’s conceptions of mathematics and physics, these claims seem natural. Professor Coifman explains, “The time vibrations of any mechanical system is a combination of sines and cosines.” However, during the early 1800s, Fourier’s claims were radical. He proposed that discontinuous functions, such as temperature distributions, could be described by combining many continuous functions; for example, an infinite number of sinusoids could represent any function, including one with multiple jumps or discontinuities. Not surprisingly, these claims were met with heavy scrutiny. In fact, during one of Fourier’s research presentations, a contemporary French mathematician Joseph Louis Lagrange reportedly exclaimed that his ideas were “nothing short of impossible.”
Although Lagrange himself made a variety of mathematical contributions during the 1800s that greatly aided modern studies of astronomy and economics, the ubiquitous power of the Fourier transform in the modern mathematical world indicates that his doubts were misguided. Today, the Fourier transform can be applied in two different ways. First, it can be used to describe continuous functions – functions providing values for every real number. In these cases, the original function is deconstructed into component sinusoidal functions at every frequency, and these are combined by the Fourier integral operation. A second type of function that the Fourier transform can be applied is one consisting of numerous discrete values, a common form of data obtained from scientific experimentation. In these cases, a Fourier series is calculated by the sum of a series of sine or cosine functions at discrete frequencies.
Even though Fourier derived this method of analysis in the early 1800s, its general applicability to solving scientific problems in an efficient manner was greatly limited until the advent of modern electronic computation. Before the 1960s, Fourier transforms for newly discovered functions that were based on experimental data or not found in reference tables necessitated an intimidating and tedious amount of arithmetic calculations: for every n data points (usually well over 1,000 for many studies), approximately n additions and n2 multiplications were needed to perform the transform. As stated by Professor Moore, “the sheer computational difficulties” and effort involved in performing a Fourier transform on paper restricted its widespread usage in the sciences.
However, in 1965, Princeton mathematics professors James Cooley and John Tukey developed the so-called “Fast Fourier Transform” (FFT) algorithm at IBM’s Watson Research center. This computational method was an integral discovery that greatly expanded the potential use of Fourier methods to solve scientific problems, and was arguably based on the work of German mathematician Carl Freidrich Gauss in the early 1800s. The FFT exponentially reduces the number of multiplication steps needed to analyze curves. For example, if a curve consisted of 8 samples, then 82, or 64, multiplication steps would be needed for analysis via traditional Fourier methods. The FFT breaks this curve into four irreducible sets of 2 sample points each. Then, these are recombined into two four-point transforms, and finally into the 8-point transform of interest. Since each stage requires 8 multiplications, the total number of steps required is 16, just one-fourth of the original 64. Therefore, the advent of the Fast Fourier Transform, and its counterparts, such as the related Hartley transform, have allowed for a much more widespread application to many scientific fields dealing with fluctuations or wave-like phenomena. Professor Moore agrees, “One of the reasons that the Fourier transform has become so pervasive today is because the computation has become routine.”
Peter Moore has been using Fourier transform methods to solve biological problems since his days as a graduate student, explaining that Fourier transforms are nothing short of “pervasive in the physical sciences.” He primarily deals with its applications to crystallography and nuclear magnetic resonance (NMR) imaging, noting that Fourier transforms are essentially “built into the physics of many scientific phenomena.” Further applications lie in the realm of imaging and spectroscopy. Modern NMR methods collect data in the form of electrical signals as a function of time, but display them as a function of frequency. A Fourier transform is used to proceed from the time domain to the frequency domain. “I use them all the time,” Moore notes, “I can’t even remember a time when I wasn’t using them.” Professor Coifman, whose research deals with inventing new complex transformations, agrees with the ubiquity of this method, “Nowadays, there isn’t a single electronic instrument that doesn’t use a Fourier transform.” Indeed, the Fourier transform today is vital to audio and video compression; without it, MP3 players and digital video cameras would not exist!
“The image formation carried out by any focusing lens is most accurately described by the Fourier transform,” concludes Professor Moore. From this, the true ubiquity of the Fourier transform reveals itself; not only does it possess fundamental applications in modern electronics, biology, chemistry, and medicine, it is rooted in stimuli processing mechanisms that we relied upon for millions of years before the work of Jean-Baptiste-Joseph Fourier. Professor’s Coifman’s concluding words are therefore undeniable: the Fourier transform is “one of the most fundamental mathematical tools in today’s world.”
Further Reading:
Bracewell, Ronald. The Fourier Transform & Its Applications. New York: McGraw-Hill, 1999.
Bracewell, R. (1989) “The Fourier Transform,” Scentific American, June: 62-69.
Brigham, Oran. The Fast Fourier transform and its applications. New York: Prentice-Hall:
1988.
About the Author:
Rohit Thummalapalli is a sophomore in Ezra Stiles College. An aspiring Molecular, Cellular, and Developmental Biology and Applied Mathematics major, Rohit serves as Subscription Manager for the Yale Scientific Magazine. He also involved in MATHCOUNTS Outreach and works in Assistant Professor Elke Stein’s neurobiology research lab in the MCDB department.
Acknowledgements:
The author would like to thank professors Ronald Coifman and Peter Moore.




Tuesday, December 13, 2011

FSL atlas for DTI !

JHU one:




Tract one:

Tuesday, November 29, 2011

All my papers ... need more ... citation is not high ... sign ...

Y Wang, J Xiang, R Kotecha, J Vannest, Y Liu… - Brain topography, 2008 - Springer
... &) 4 J. Xiang 4 R. Kotecha 4 J. Vannest 4 Y. Liu 4 D. Rose 4 M. Schapiro 4 T. Degrauw MEG
Center, Division of Neurology, Cincinnati Children's Hospital Medical Center, MLC 2015, 3333
Burnet Avenue, Cincinnati, OH 45229-3039, USA e-mail: yingying.wang@cchmc.org ...
Cited by 7 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 7 versions - Import into EndNote

Y Wang, J Xiang, J Vannest, T Holroyd… - Clinical …, 2011 - Elsevier
... Yingying Wang a , b , c , Corresponding Author Contact Information , E-mail The Corresponding
Author , Jing Xiang b , Jennifer Vannest b , c , Tom Holroyd d , Daria Narmoneva a , Paul Horn ...
d MEG Core Facility, National Institute of Mental Health, Bethesda, MD, United States. ...
Cited by 1 - Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 4 versions - Import into EndNote

Y Wang, J Xiang, DF Rose, T Holroyd… - 17th International …, 2010 - Springer
... We adopted accumulated spectrograms because the MEG data recorded from the brain commonly
have very strong low-frequency ... Author: Yingying Wang Institute: Cincinnati Children's Hospital
Medical Center Street: 3333 Burnet Avenue, MLC 2015 City: Cincinnati, OH 45229 ...
Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 3 versions - Import into EndNote
 

Second author paper (2)


J Xiang, Y Wang, Y Chen, Y Liu… - Journal of …, 2010 - thejns.org
... Jing Xiang, MD, Ph.D. 1, 2 , Yingying Wang, M.Sc ... 1 , Nat Hemasilpin, MSEE 1 , Ki Lee, MD 2 ,
Francesco T. Mangano, DO 3 , Blaise Jones, MD 4 , and Ton deGrauw, MD, Ph.D. 2. 1 MEG Center,
2 Division of Neurology, 3 Division of Neurosurgery, and 4 Department of Radiology ...
Cited by 4 - Related articles - OhioLINK OLinks - Check Library Holdings - All 3 versions - Import into EndNote


X Huo, Y Wang, R Kotecha, EG Kirtman… - Brain topography, 2011 - Springer
Abstract Recent studies in adults have found consistent contralateral high gamma activities
in the sensorimotor cortex during unilateral finger movement. However, no study has
reported on this same phenomenon in children. We hypothesized that contralateral high ...
Cited by 1 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 3 versions - Import into EndNote
 

Third author paper (8)


J Xiang, Y Liu, Y Wang, EG Kirtman, R Kotecha… - Epileptic Disord, 2009 - jle.com
Summary: Purpose. Invasive intracranial recordings have suggested that high-frequency
oscillation is involved in epileptogenesis and is highly localized to epileptogenic zones. The
aim of the present study is to characterize the frequency and spatial patterns of high- ...
Cited by 9 - Related articles - All 4 versions - Import into EndNote


R Zhang, T Wu, Y Wang, H Liu, Y Zou, W Liu… - Seizure, 2011 - Elsevier
... Rui Zhang a , Ting Wu b , Yingying Wang d , e , Hongyi Liu a , Yuanjie Zou a , Wen Liu c , Jing
Xiang d , Chaoyong Xiao c , Lu Yang b and Zhen Fu f , Corresponding Author Contact ... b MEG
Center, Brain Hospital Affiliated of Nanjing Medical University, Nanjing, China. ...
Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 2 versions - Import into EndNote

 

Gamma oscillations in the primary motor cortex studied with MEG

X Huo, J Xiang, Y Wang, EG Kirtman… - Brain and …, 2010 - Elsevier
In recent years, there has been a growing interest on the role of gamma band (> 30Hz)
neural oscillations in motor control, although the function of this activity in motor control is
unknown clearly. With the goal of discussing the high frequency sources non-invasively ...
Cited by 4 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 6 versions - Import into EndNote 


X Wang, J Xiang, Y Wang, M Pardos… - … : The Journal of …, 2010 - Wiley Online Library
... of motor cortex in the pediatric migraine was altered, this study provides pilot data for further
investigation of the cerebral mechanisms of migraine with MEG and advanced ... (c) Analysis and
Interpretation of DataXiaoshan Wang, Jing Xiang, Xiaolin Huo, Yingying Wang, Andrew D ...
Cited by 1 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 3 versions - Import into EndNote


J Xiang, Z Xiao, Y Wang, Y Feng, H Qiao… - International Congress …, 2007 - Elsevier
The present study aimed to investigate whether magnetoencephalography (MEG)
information could result in the detection of subtle anatomical abnormalities at re-review of
conventional magnetic resonance imaging (MRI) by a new MEG guided post-image ...
Related articles - Find it with OLinks - Full-Text @ Pratt Library - Import into EndNote

 

Modeling the developmental patterns of auditory evoked magnetic fields in children

R Kotecha, M Pardos, Y Wang, T Wu, P Horn… - PloS one, 2009 - dx.plos.org
... Rupesh Kotecha 1 , Maria Pardos 1 , Yingying Wang 1 , Ting Wu 1 , Paul Horn 1 , 2 , David Brown
3 , 4 , Douglas Rose 1 , Ton deGrauw 1 , Jing Xiang 1 *. 1 MEG Center, Department of Neurology,
Cincinnati Children's Hospital Medical Center, Cincinnati, Ohio, United States of ...
Cited by 9 - Related articles - Cached - Find it with OLinks - Full-Text @ Pratt Library - All 9 versions - Import into EndNote 


Y Liu, J Xiang, Y Wang, JJ Vannest, AW Byars… - Brain topography, 2008 - Springer
... Yinhong Liu Æ Jing Xiang Æ Yingying Wang Æ Jennifer J. Vannest Æ Anna W. Byars Æ Douglas
F. Rose ... spatial and frequency differences between recognizing concrete and abstract words using
a 275 channel whole head magnetoencephalography (MEG) system. ...
Cited by 5 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - BL Direct - All 3 versions - Import into EndNote



R Kotecha, J Xiang, Y Wang, X Huo… - International Journal of …, 2009 - Elsevier
... a MEG Center, Department of Neurology, Cincinnati Children's Hospital Medical Center, 3333
Burnet Avenue, Cincinnati, OH, 45220, USA. Received 23 September 2008; accepted 31 October
2008. Available online 14 November 2008. Abstract. ... 3.1. Physical MEG waveforms. ...
Cited by 9 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 3 versions - Import into EndNote 


4th, 6th, 7th, or 10th author paper (3+4)


…, SK Holland, J Xiang, Y Wang - … , Proceedings of the …, 2010 - ieeexplore.ieee.org
Abstract This paper introduces a source localization technique that exploits the high
temporal resolution of neuronal magnetoencephalography (MEG) data to locate the
originating sources within the head. A traditional frequency beamforming algorithm was ...
Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 3 versions - Import into EndNote
 

Aberrant high-gamma oscillations in the somatosensory cortex of children with cerebral palsy: A meg study

…, M Bryce, S Huang, X Huo, Y Wang… - Brain and …, 2011 - Elsevier
Objective: Our study is to investigate somatosensory dysfunction in children with spastic
cerebral palsy (CP) using magnetoencephalography (MEG) and synthetic aperture
magnetometry (SAM). Methods: Six children with spastic CP and six age-and gender- ...
Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 3 versions - Import into EndNote


…, Y Chen, L Meng, X Wang, Y Wang - 17th International …, 2010 - Springer
Our previous studies have demonstrated that Morlet-wavelet transform with an extra large
sigma value could precisely determine the frequency signatures of neuromagnetic signals.
Unfortunately, the increase of frequency sensitivity is associated with a decrease of ...
Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 3 versions - Import into EndNote
 
…, D Rose, A Byars, D Brown, JH Seo, Y Wang… - Epilepsy research, 2010 - Elsevier
... Methods. In this magnetoencephalography (MEG) study, 10 patients and 10 age- and
gender-matched healthy controls were investigated with the multi-feature mismatch negativity
(MMN) paradigm. ... Antiepileptic drugs tapered more than 24 h before MEG study, Side, ...
Cited by 3 - Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 4 versions - Import into EndNote
 
…, T Kujala, J Xiang, J Vannest, Y Wang… - 17th International …, 2010 - Springer
... Ki Heyeong Lee1, Hisako Fujiwara1, Teija Kujala3, Jing Xiang1, Jennifer Vannest1,2, Yingying
Wang1, Nat ... In this MEG study we aimed to register two types of event-related ... which gives possibility
to simultaneously register acoustically and visually presented words (Wang et al. ...
Cited by 1 - Related articles - Find Full-Text at UC - OhioLINK OLinks - Check Library Holdings - All 2 versions - Import into EndNote
 
 
…, D Brown, P Horn, Y Wang… - International …, 2011 - informahealthcare.com
... 1,6 David Brown, 3 Paul Horn, 1,4 Yingying Wang, 1 Hisako Fujiwara, 1 Jing Xiang, 1 Marielle
A. Kabbouche, 1,6 Scott W. Powers, 5,6 and Andrew D. Hershey, 1,6 ... We would like to thank Nat
Hemasilphin and Elliah Kirtman for their technical assistance in MEG recordings. ...
Related articles - Find Full-Text at UC - Check Library Holdings - All 4 versions - Import into EndNote

Y Chen, J Xiang, EG Kirtman, Y Wang… - Clinical …, 2010 - Elsevier
... children. Methods. Sixty healthy children and 20 adults were studied with a whole-head
magnetoencephalography (MEG) system. The adults were included to find out when
the markers stabilize. Visual ... 2.3. MEG recordings. The MEG ...
Related articles - Find it with OLinks - Full-Text @ Pratt Library - All 4 versions - Import into EndNote

Monday, November 28, 2011

Good point from Dr. Nicols.

Dear Torsten,

Randomise doesn't give you the critical cluster size threshold by default, but it's easy to obtain.  If you add the -N option to randomise, for each corrp image you'll get a .txt file that gives the permutation distribution of the maximum statistic (whether that's max voxel T, max cluster size, max TFCE score, whatever).  If you find the 95%ile of that distribution, that's the 5% critical threshold based on permutation.  Loading that in Matlab and finding this percentile is easy:

MaxC=load('permdist.txt');
Nperm=length(MaxC);
sMaxC=sort(MaxC);
Level=0.05;
CritC=sMaxC(ceil(Nperm*Crit))

but you're right, it's something that I always report in papers.  I'll try to get that printed out with, say, the -v flag, in an future version.


My problem with 3dcluster, alphasim, and any Monte Carlo based cluster inference tool, is that you have to believe in the Gaussian autocorrelation *and* stationarity that the simulations are based on.  VBM data are widely acknowledged to exhibit nonstationary smoothness, but whenever I've looked at a FWHM image from FMRI data I see hints of structure there too.  Randomise or any permutation-based procedure will automatically account for any nonstationarity in the data, and is not vulnerable to errors in estimated FWHM smoothness (even if the data were stationary).

Permutation is "exact", in that it guaranteed to control false positive risk with very weak assumptions, but it's not perfect: Parametric models can provide better power *when* all the assumptions are satisfied [1].  But if lots and lots of people find better results with the Monte Carlo method than with permutation, it might be that the Monte Carlo method is inflating significances.  The traditional way of comparing methods, with Monte Carlo simulations of homogeneous smooth Gaussian noise, won't reveal this (as the parametric assumptions *define* the Monte Carlo method, and permutation can't out-perform that).  A large body of null data with real (i.e. messy) spatial structured noise would be needed to tested to see if there is a substantial statistical inefficiency in permutation cluster size inference.

Hope this helps!

-Tom

[1] However, in all the standard settings, e.g. t-tests, permutation tests have asymtotic relative efficiency of 1, i.e. will be as powerful as parametric tests when larger and larger sample sizes are considered.
[2] Random Field Theory makes these assumptions too, and additionally approximations in the P-value formulas, but these are just more reasons not to use RFT---thought RFT does at least have a way of handling nonstationarity.