Duke-UNC Brain Imaging and Analysis Center
BIAC Forums | Profile | Register | Active Topics | Members | Search | FAQ
Username:
Password:
Save Password   Forgot your Password?
 All Forums
 Support Forums
 Analysis Software Support
 MatLab - Calculating a Difference Image
 New Topic  Reply to Topic
 Printer Friendly
Author Previous Topic Topic Next Topic  

jvbinzak
Starting Member

1 Posts

Posted - Jun 21 2011 :  5:48:57 PM  Show Profile  Reply with Quote
I'm new to data processing so for give me if this comes off poorly. In the instances where my subjects QA shows a substancial RF spike or head displacement, I have been instructed to "Calculate a Difference Image" using MATLAB. My understanding is that this will allow me to remove the spikes and movement from my data before the first level of analysis.

Does anyone have experience dealing with managing large RF spikes and moments of head movement in their data? If so, what all is possible fix the data, and how do I go about it?

Specifically, can someone outline a step by step procedure to calculate a difference image using MATLAB? If so please make it as overt/fool proof as possible.

John Binzak
Research Assistant
Center for Child and Family Policy

dvsmith
Advanced Member

USA
218 Posts

Posted - Jun 22 2011 :  10:39:53 AM  Show Profile  Visit dvsmith's Homepage  Reply with Quote
Calculating the difference image (timepoint-to-timepoint difference) would mainly just allow you to better visualize the spikes and artifacts. If you want to do that, it's pretty easy to do:

funcfile = readmr(funcfile);
funcfile.data = diff(funcfile.data, 1, 4);

If you want to quantify and model spikes, you should take a look at an FSL tool called fsl_motion_outliers. It's also very easy to use.
Usage: fsl_motion_outliers <unfiltered 4D image> <number of dummy scans> <output confound file>
Go to Top of Page
  Previous Topic Topic Next Topic  
 New Topic  Reply to Topic
 Printer Friendly
Jump To:
BIAC Forums © 2000-2010 Brain Imaging and Analysis Center Go To Top Of Page
This page was generated in 0.61 seconds. Snitz Forums 2000