Skip to content

AFNI/NIfTI Server

Sections
Personal tools
You are here: Home » AFNI » Community » Message Board

History of AFNI updates   Suggested reading for beginners   Add yourself to the AFNI map!

Re: AlphaSim and 3dmerge twice

|
June 20, 2012 04:17PM
Two smoothing operations will be roughly additive. This description taken from the wikipedia article on Gaussian blurring describes this:

Applying multiple, successive gaussian blurs to an image has the same effect as applying a single, larger gaussian blur, whose radius is the square root of the sum of the squares of the blur radii that were actually applied. For example, applying successive gaussian blurs with radii of 6 and 8 gives the same results as applying a single gaussian blur of radius 10, since \sqrt{6^2 + 8^2} = 10.

So two successive 8mm blurs would be the same as applying sqrt(64+64) = 11.3 mm gaussian blur.

3dFWHMx is a program that was intended to be used to estimate smoothness from the noise of a time series. 3dFWHMx does have a -2difMAD option that has been used in the past for estimating smoothness from structural images, but it's not clear this method is effective or useful. You can search our message board for previous posting regarding this option. I haven't had much luck with it though on MRI structural images. 3dLocalstat with "-stat fwhm" gives me something a little more useful and gives somewhat reasonable numbers. This is a short script to play around with (just an ad hoc attempt to get an estimate for smoothness). You may also want to look at 3dBlurtoFWHM.

#!/bin/tcsh
# est_smooth.csh
#
# estimate smoothness for a structural image (if such a thing is possible)
#
# usage:
# est_smooth.csh dset nbhd_size
#
# example:
#
# 3dmerge -1blur_fwhm 8.0 -prefix anat_ns_blur2 anat_ns_blur1+orig
# tcsh est_smooth.csh anat_ns_blur1+orig 5
#

set dset = $1
set nbhd_size = $2
rm -f temp_localblur+orig.*

# estimate smoothness measure locally
3dLocalstat -stat FWHM -mask $dset -prefix temp_localblur -nbhd "RECT($nbhd_size,$nbhd_size, $nbhd_size)" $dset

set medx = `3dBrickStat -median -positive temp_localblur+orig'[0]'`
set medy = `3dBrickStat -median -positive temp_localblur+orig'[1]'`
set medz = `3dBrickStat -median -positive temp_localblur+orig'[2]'`

echo "Median blur estimates in x,y,z are $medx[2], $medy[2], $medz[2]"
Subject Author Posted

AlphaSim and 3dmerge twice

David Zhu June 19, 2012 03:01PM

Re: AlphaSim and 3dmerge twice

Daniel Glen June 19, 2012 04:04PM

Re: AlphaSim and 3dmerge twice

David Zhu June 20, 2012 11:36AM

Re: AlphaSim and 3dmerge twice

Daniel Glen June 20, 2012 04:17PM



Author:

Your Email:


Subject:


Spam prevention:
Please, enter the code that you see below in the input field. This is for blocking bots that try to post this form automatically. If the code is hard to read, then just try to guess it right. If you enter the wrong code, a new image is created and you get another chance to enter it right.
CAPTCHA
Message:

Powered by Plone

phorum.org

This site conforms to the following standards: