Hi, Mark, Great thanks! You are right!!! Last tonight, I figured it out when I was nearly giving up. I suddenly realized the reason of the problem may be FOV. So, I make the bcot1 and roi the same FOV and voxel size using flirt -usesqform before applying *****. Thanks again! Best, Zhijang Wang 2014-04-08 14:15 GMT+08:00 Mark Jenkinson <[log in to unmask]>: > Hi, > > If I assume that the names bcoT1 (in the initial flirt commands) is > equivalent to nbT1 (in the header info section) and that roi (flirt calls) > is equivalent to ROImask (headers), then I think I see what the problem is. > These images have completely different voxel sizes, and now I suspect that > they also have different FOV. In order to apply one transform to another > image, the second image has to be in the same space as the original input > image. So that means that the voxel size and FOV should be the same for > "bcoT1" and "roi". You should be able to overlay these two images within > FSLView. If not, then you cannot apply the transformation as the image > "roi" is not in the correct space to start with. > > Normally I would expect that "roi" is just a copy of the image "bcoT1" > but the intensities contain zeros everywhere except for within the ROI, > where they would contain the value 1. > > If you describe what "roi" is, and how it relates to "bcoT1" then we > should be able to help. > > All the best, > Mark > > > > On 6 Apr 2014, at 10:09, Zhijiang Wang <[log in to unmask]> wrote: > > Dear Mark, > > Great thanks to your reply! > Please let me re-tell you my question: > ####################################################################### > flirt -in bcoT1 -ref B0 -cost corration -o > t12b0 -omat t_mat > > # t12b0 is OK. > > flirt -in roi -ref B0 -o b0_roi -init t_mat > -applyxfm > # But, when applying the transformation matrix to the roi > mask, the resulted b0_roi.nii can not be alighned with b0/ > > Why? > > I am so confused/ > > ############################################################################ > > I just send my .nii file to you and fsl-email-list, but received "Fail in > sending this email because the size of files is greater than 200K". > > So, I have to attach some head information (*.mat, sform, qform). > > Could anybody tell me where the problem is? > > Great thanks! > > > > B0.head.mat > > -2.20312500000000 0 0 143.203125000000 > 0 2.20312500000000 0 -124.132568359375 > 0 0 2.20000004768372 -53.7819283008575 > 0 0 0 1 > fslorient -getqform B0 > -2.20312 0 -0 141 > 0 2.20312 -0 -121.929 > 0 0 2.2 -51.5819 > 0 0 0 1 > fslorient -getsform B0 > -2.20312 0 0 141 > 0 2.20312 0 -121.929 > 0 0 2.2 -51.5819 > 0 0 0 1 > > bT1.head.mat > 1.32840895652771 0.0126800928264856 0.0472275912761688 -104.481479177251 > -0.0154902637004852 0.999688148498535 -0.0220899917185307 -113.658335920423 > -0.0631656423211098 0.0215135142207146 0.998639881610870 -120.418306723237 > 0 0 0 1 > fslorient -getsform nbT1 > 1.32841 0.0126801 0.0472276 -103.093 > -0.0154903 0.999688 -0.02209 -112.696 > -0.0631656 0.0215135 0.99864 -119.461 > 0 0 0 1 > fslorient -getqform nbT1 > 1.32873 0 0.0436194 -99.8026 > 0 1 0 -80.1566 > -0.0580138 0 0.999048 -101.699 > 0 0 0 1 > > ROImaskfile.head.mat > -3.12486195564270 -0.0109078036621213 0.0366514474153519 101.246426661499 > -0.0109082190319896 3.12498092651367 1.35525271560688e-20 -122.445417311974 > 0.0272702574729919 9.51912588789128e-05 4.19984006881714 -71.5327780273146 > 0 0 0 1 > > fslorient -getsform ROImask > -3.12486 -0.0109078 0.0366514 98.1473 > -0.0109082 3.12498 1.35525e-20 -119.331 > 0.0272703 9.51913e-05 4.19984 -67.3056 > 0 0 0 1 > fslorient -getqform ROImask > -3.12486 -0.0109078 0.0366216 98.1473 > -0.0109082 3.12498 -5.22123e-08 -119.331 > 0.0272481 9.51524e-05 4.19984 -67.3056 > 0 0 0 1 > > > > > > > 2014-03-10 19:08 GMT+08:00 Mark Jenkinson <[log in to unmask]>: > >> Hi, >> >> It is difficult to know what is going wrong from your description here. >> The commands, subject to minor typos (e.g. corration), look correct. >> Could the problem be related to the output of the second command not >> using a floating-point datatype (you can fix this with "-datatype float" in >> the second flirt call) and hence it truncating a lot of the mask values to >> zero? >> >> If that is not the case then please give us a more informative >> description of the problem. >> >> All the best, >> Mark >> >> >> >> On 4 Mar 2014, at 13:15, Zhijiang Wang <[log in to unmask]> wrote: >> >> Dear fslers, >> >> I want to register a ROI mask into B0, but I am confused the followings: >> >> flirt -in bcoT1 -ref B0 -cost corration -o t12b0 -omat >> t_mat >> >> # t12b0 is OK. >> >> flirt -in roi -ref B0 -o b0_roi -init t_mat -applyxfm >> # But, when applying the transformation matrix to the roi mask, the >> resulted b0_roi.nii can not be alighned with b0/ >> >> Why? >> >> I am so confused/ >> >> who can share with your suggestions and comments. >> >> Best, >> >> >> > >