Hi Alex, OK, sounds like you have done the sensible things. Could you put your images and weighting volumes somewhere for us to download (web page or ftp)? We'll take a look at them at this end and see if there is anything that can be done. All the best, Mark On Fri, 15 Oct 2004, Fornito, Alexander wrote: > Hi, > Thanks for your help - it's clarified a few things. However, I'm still > having a bit of trouble. Maybe you can tell me where I'm going wrong. > I've done what you suggested, trying to create a weight in the input > space. > I drew a mask around there areas of signal loss (orbital frontal and > inferior temporal areas) and saved this out as a binary image file (in > analyze format), with the mask valued at 0 and the rest (including > non-brain) as 1. I then included this image as the input weight using > the flirt GUI when registering the epi to the subjects T1, but it > appears to have made little difference to the actual registration. > I then tried again, this time drawing a mask around the ventricles and > weighting them as 2, the signal dropout as 0, and the rest as 1, but > again, there appears to be little difference in the registration > results. > Am I missing something? > Thanks again, > Alex > > -----Original Message----- > From: Mark Jenkinson [mailto:[log in to unmask]] > Sent: Friday, October 01, 2004 10:27 PM > To: [log in to unmask] > Subject: Re: [FSL] Registration woes > > Hi Alex, > > I agree with Tim. Especially in that it is very important to > distinguish > masking from weighting - they are quite different operations. > > I would say that there isn't much advantage in smoothing the cost > weights > though, so you could easily ignore that step. > > Also, I never recommend drawing weights based on a non-weighted > registration > (using mricro's linked view or anything else) as it is only likely to be > much good > if the registration is good, which is not going to be true when you > really need > the weights the most. Either draw it directly on the functions where > you can > see the signal loss (and use the -inweight) or draw a general one in the > appropriate > areas on a structural or standard-space image (and use the -refweight > option). > The latter saves a lot of time if you are doing this many times on > different > scans. > > It is also worth noting that the larger the ignored area in the weights > (voxels > with zero or near-zero values) then the less robust the initial search > phase of > the registration becomes. So for weights with relative large areas > being set > near zero it is best to use "no search" either via the GUI pull-down > options > or by using -nosearch in the command line. This will then require that > the > initial positions of the images are no more than say 20-30 degrees out > of > alignment, but as long as the slicing (sagittal, axial, oblique, etc) is > the same > or can be made the same (or nearly the same) using avwswapdim then > there should be no problems. > > Hope this is also useful. > > All the best, > Mark > > > > Tim Behrens wrote: > > >Hi Alex, > > > >I'll have a go at answering these - Mark can correct me when I get it > >wrong.. > > > > > > > >>1 - Would I have to mask just the initial highres, or both the initial > >>highres and task epi? If the latter, would I just mask the middle > volume, > >>as this is what is used in the registration? If so, my 4D file > contains > >>160 volumes. Would the middle correspond to the 79th, or the 80th? > >> > >> > > > >the cost function mask is a separate 3D image in reference (-refweight) > or > >input (-inweight) space ( you can actually apply both weightings if you > >like, but they have to be in the right spaces). Note that Cost function > >weighting is different from masking the original images, as masking > >introduces artificial edges which will drive the registration - usually > >not a good idea. > > > > > > > > > >>2 - The CBU page recommends co-registering the epi to the structural > and > >>using the "yoke" function in MRIcro to identify points where there is > >>signal dropout in the epi relative to the structural. These are the > >>portions that should be masked out. If I register the unmasked epi to > the > >>structural, wouldn't this yield a (relatively) poor correspondance > between > >>the images, so that using this as the basis for identifying signal > dropout > >>would not be the best way to go about it? Are there any alternative > ways > >>for idetifying which areas to mask out? > >> > >> > >> > > > >If you use -inweight, your cost function weight mask will be directly > in > >input space so you can just downweight the areas with susceptibility > >artefact, and upweight the boundaries you trust ( lateral walls of the > >ventricles etc.. ). If you wanted to you could apply a _generic_ MNI > space > >mask to the reference space (-refweight) which would do the same thing > - > >downweight occipital pole etc and upweight the ventricle walls. I guess > >this would be slightly less accurate than doing it individually for > each > >input volume, but I would think you could downweight pretty big areas > >without really detracting from the registration, so I would think this > >should be fine. > > > > > > > > > >>3 - Following from 2, would there be any problems in doing a straight > >>image subtraction to identify the areas to be masked? It seems to me > that > >>this would be plagued by the same problems as in 2, but that it would > be a > >>quicker way to go about it? > >> > >> > > > >I think this is answered by the previous point (??) > > > > > > > >>4 - The CBU page recommends smoothing the mask image (8mm in > accordance > >>with the SPM default). Would I also need to do this with FEAT? Would I > use > >>5mm, as this is the FEAT default (or alternatively, whatever level I > set > >>it to)? > >> > >> > >> > > > >I would think smoothing the mask image would be a good idea. However, > the > >smoothness here is different from the smoothness in Feat - here we are > >smoothing because of uncertainty in the location of the structures of > >interest, not due to GRF theory etc. In this case, I would thing 8mm > >smoothing of the mask image would be fine. > > > > > > > >>5 - As an aside, can someone tell me why the average epi template in > SPM > >>would be a poor choice as an intial highres? Does the fact that it's > not > >>the same brain as your subject introduce more problems than the > advantage > >>of having the added contrast? > >> > >> > >> > > > >Usually going to an initial highres, you would want to be able to > choose a > >low number of degrees of freedom to get a really robust registration. > If > >you use a reference brain which is a different global shape to the > input > >brain, you will need a high number of degrees of freedom even to get > close > >to a good registration, so I would think it would not be particularly > >advantageous to choose an averaged brain as an initial highres. > > > >Hope this is useful (and accurate!!) > > > > > >T > > > > > > > > > > > > > > > >>Sorry if some of these questions seem a bit basic - still trying to > get my > >>head around all of this! > >>Thanks heaps, I really appreciate your ongoing help! > >>Alex > >> > >> > >> >