Hi Stefan,<br><br>Thank you very much for your quick response. I set the default pixel value to 500 and those holes did indeed have a value of 500.<br><br>Is there a way to specify that data in certain slices should always remain in the same plane while data in other slices are allowed to deform in any direction? MovingImageDerivativeScales seems like it might help if there is a way to specify that parameter for specific parts of the data.<br>
<br>Alternatively, is there a way to specify fixed image landmark positions AND moving image landmark positions as inputs for registration? I was looking at the thin-plate splines transform which appears to only allow fixed landmarks as input. Ideally, I would want to continue using a BSpline transform but be able to provide these landmarks in order to help guide the registration.<br>
<br>Thank you,<br>Kristen<br><br><div class="gmail_quote">On Thu, Aug 11, 2011 at 6:00 AM,  <span dir="ltr"><<a href="mailto:elastix-request@bigr.nl" target="_blank">elastix-request@bigr.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


Send elastix mailing list submissions to<br>
        <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://lists.bigr.nl/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/mailman/listinfo/elastix</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:elastix-request@bigr.nl" target="_blank">elastix-request@bigr.nl</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:elastix-owner@bigr.nl" target="_blank">elastix-owner@bigr.nl</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of elastix digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Re: Elastix Ubuntu New User (Stefan Klein)<br>
   2. Re: Elastix Ubuntu New User (Robinson, Jamie)<br>
   3. 3D Bspline registration (Kristen Tgavalekos)<br>
   4. Re: 3D Bspline registration (Stefan Klein)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Wed, 10 Aug 2011 13:38:07 +0200<br>
From: Stefan Klein <<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>><br>
Subject: Re: [elastix] Elastix Ubuntu New User<br>
To: "Robinson, Jamie" <<a href="mailto:j.robinson.10@aberdeen.ac.uk" target="_blank">j.robinson.10@aberdeen.ac.uk</a>><br>
Cc: "<a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a>" <<a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a>><br>
Message-ID: <<a href="mailto:4E426D9F.5000704@erasmusmc.nl" target="_blank">4E426D9F.5000704@erasmusmc.nl</a>><br>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed<br>
<br>
Hi Jamie,<br>
That probably explains the error. You can better split the image stacks,<br>
and register the anterior views separately and the posterior views<br>
separately.<br>
<br>
Right now, I don't think elastix reads the images correctly even, since<br>
you use a 2d parameter file, and give a multi-image stack as input<br>
(which could be considered as a 3d image probably). Not sure what<br>
happens then.<br>
<br>
Hope this solves the problem.<br>
kind regards,<br>
Stefan<br>
<br>
On 10/08/2011 10:20, Robinson, Jamie wrote:<br>
> Hi Stefan,<br>
><br>
> Both image are 256x1024 with voxel size 2.4x2.4x1 mm. Each file is an<br>
> image stack with 2 slices - an anterior and posterior view.<br>
><br>
> Thanks,<br>
> Jamie<br>
><br>
><br>
> ------------------------------------------------------------------------<br>
> *From:* <a href="mailto:elastix-bounces@bigr.nl" target="_blank">elastix-bounces@bigr.nl</a> [<a href="mailto:elastix-bounces@bigr.nl" target="_blank">elastix-bounces@bigr.nl</a>] on behalf of<br>
> Stefan Klein [<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>]<br>
> *Sent:* Tuesday, August 09, 2011 9:57 PM<br>
> *To:* <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
> *Subject:* Re: [elastix] Elastix Ubuntu New User<br>
><br>
> Hi Jamie,<br>
> Thanks for your interest in elastix. The warnings are not a problem.<br>
> They just inform you that some parameters were not specified and that<br>
> elastix takes a default value.<br>
> The error means that the fixed and moving image do not have a<br>
> sufficiently large initial overlap. Or maybe your moving image is much<br>
> smaller than your fixed image. Could you email us:<br>
> - The size of the images<br>
> - The voxel size<br>
> Kind regards,<br>
> Stefan<br>
><br>
><br>
> On 9-8-2011 19:55, Robinson, Jamie wrote:<br>
>> I am working on developing some software that compares two nuclear<br>
>> medicine whole body bone scans both intra and intersubject.<br>
>> Previously, I have used a registration package bUnwarpJ in ImageJ for<br>
>> this purpose. However, I was not satisfied with the results (it uses a<br>
>> simple sum of squared intensity differences as a metric) and would now<br>
>> like to try something a bit more sophisticated and stumbled across<br>
>> elastix.<br>
>><br>
>> I followed the user manual on how to install the package and<br>
>> downloaded the default parameter configuration files. However, upon<br>
>> running the command<br>
>><br>
>> elastix -f FS1.dcm -m FS2.dcm -p ParameterFiles/Parameters_Affine.txt<br>
>> -p ParameterFiles/Parameters_BSpline.txt -out .<br>
>><br>
>> I am met with a whole bunch of warnings before the software fails. I<br>
>> have attached the log file but have no idea why the software fails.<br>
>> Can anyone offer any help.<br>
>><br>
>> Thanks,<br>
>> Jamie<br>
>><br>
>><br>
>> _______________________________________________<br>
>> elastix mailing list<br>
>> <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
>> <a href="http://lists.bigr.nl/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/mailman/listinfo/elastix</a><br>
<br>
--<br>
Stefan Klein<br>
<a href="tel:%2B31%2010%207043049" value="+31107043049" target="_blank">+31 10 7043049</a><br>
<a href="http://www.bigr.nl/people/StefanKlein" target="_blank">http://www.bigr.nl/people/StefanKlein</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Wed, 10 Aug 2011 11:55:18 +0000<br>
From: "Robinson, Jamie" <<a href="mailto:j.robinson.10@aberdeen.ac.uk" target="_blank">j.robinson.10@aberdeen.ac.uk</a>><br>
Subject: Re: [elastix] Elastix Ubuntu New User<br>
To: Stefan Klein <<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>><br>
Cc: "<a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a>" <<a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a>><br>
Message-ID:<br>
        <<a href="mailto:AB45752CD13C724FB1AAA26FDE4F6E37019926EE@AMSPRD0104MB121.eurprd01.prod.exchangelabs.com" target="_blank">AB45752CD13C724FB1AAA26FDE4F6E37019926EE@AMSPRD0104MB121.eurprd01.prod.exchangelabs.com</a>><br>


<br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Hi Stefan,<br>
<br>
The first solution works perfectly. In fact I only need to register one of the views as both anterior and posterior views are acquired simultaneously. Therefore, once I find the deformation field that registers the anterior views I can quickly register the posterior views with the mirror transformation. I look forward to using Elastix to help me with my task.<br>



<br>
Thanks,<br>
Jamie<br>
________________________________________<br>
From: Stefan Klein [<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>]<br>
Sent: Wednesday, August 10, 2011 12:38 PM<br>
To: Robinson, Jamie<br>
Cc: <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
Subject: Re: [elastix] Elastix Ubuntu New User<br>
<br>
Hi Jamie,<br>
That probably explains the error. You can better split the image stacks,<br>
and register the anterior views separately and the posterior views<br>
separately.<br>
<br>
Right now, I don't think elastix reads the images correctly even, since<br>
you use a 2d parameter file, and give a multi-image stack as input<br>
(which could be considered as a 3d image probably). Not sure what<br>
happens then.<br>
<br>
Hope this solves the problem.<br>
kind regards,<br>
Stefan<br>
<br>
On 10/08/2011 10:20, Robinson, Jamie wrote:<br>
> Hi Stefan,<br>
><br>
> Both image are 256x1024 with voxel size 2.4x2.4x1 mm. Each file is an<br>
> image stack with 2 slices - an anterior and posterior view.<br>
><br>
> Thanks,<br>
> Jamie<br>
><br>
><br>
> ------------------------------------------------------------------------<br>
> *From:* <a href="mailto:elastix-bounces@bigr.nl" target="_blank">elastix-bounces@bigr.nl</a> [<a href="mailto:elastix-bounces@bigr.nl" target="_blank">elastix-bounces@bigr.nl</a>] on behalf of<br>
> Stefan Klein [<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>]<br>
> *Sent:* Tuesday, August 09, 2011 9:57 PM<br>
> *To:* <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
> *Subject:* Re: [elastix] Elastix Ubuntu New User<br>
><br>
> Hi Jamie,<br>
> Thanks for your interest in elastix. The warnings are not a problem.<br>
> They just inform you that some parameters were not specified and that<br>
> elastix takes a default value.<br>
> The error means that the fixed and moving image do not have a<br>
> sufficiently large initial overlap. Or maybe your moving image is much<br>
> smaller than your fixed image. Could you email us:<br>
> - The size of the images<br>
> - The voxel size<br>
> Kind regards,<br>
> Stefan<br>
><br>
><br>
> On 9-8-2011 19:55, Robinson, Jamie wrote:<br>
>> I am working on developing some software that compares two nuclear<br>
>> medicine whole body bone scans both intra and intersubject.<br>
>> Previously, I have used a registration package bUnwarpJ in ImageJ for<br>
>> this purpose. However, I was not satisfied with the results (it uses a<br>
>> simple sum of squared intensity differences as a metric) and would now<br>
>> like to try something a bit more sophisticated and stumbled across<br>
>> elastix.<br>
>><br>
>> I followed the user manual on how to install the package and<br>
>> downloaded the default parameter configuration files. However, upon<br>
>> running the command<br>
>><br>
>> elastix -f FS1.dcm -m FS2.dcm -p ParameterFiles/Parameters_Affine.txt<br>
>> -p ParameterFiles/Parameters_BSpline.txt -out .<br>
>><br>
>> I am met with a whole bunch of warnings before the software fails. I<br>
>> have attached the log file but have no idea why the software fails.<br>
>> Can anyone offer any help.<br>
>><br>
>> Thanks,<br>
>> Jamie<br>
>><br>
>><br>
>> _______________________________________________<br>
>> elastix mailing list<br>
>> <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
>> <a href="http://lists.bigr.nl/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/mailman/listinfo/elastix</a><br>
<br>
--<br>
Stefan Klein<br>
<a href="tel:%2B31%2010%207043049" value="+31107043049" target="_blank">+31 10 7043049</a><br>
<a href="http://www.bigr.nl/people/StefanKlein" target="_blank">http://www.bigr.nl/people/StefanKlein</a><br>
<br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Wed, 10 Aug 2011 10:29:30 -0400<br>
From: Kristen Tgavalekos <<a href="mailto:kristen.tgav@gmail.com" target="_blank">kristen.tgav@gmail.com</a>><br>
Subject: [elastix] 3D Bspline registration<br>
To: <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
Message-ID:<br>
        <<a href="mailto:CAH1hsYuG%2BBhYE-T0zom8ZFbJVhb1wBmWm0AMpovMsctRKait-g@mail.gmail.com" target="_blank">CAH1hsYuG+BhYE-T0zom8ZFbJVhb1wBmWm0AMpovMsctRKait-g@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi Elastix team,<br>
<br>
I'm working on registering 3D images with a BSpline transform. They are<br>
limited angle tomography images with different depths (a 10 slice<br>
difference) and the resolution is 0.1 x 0.1 x 1 mm. My problem is that the<br>
final registered image has blank holes in it. Typically the top and bottom<br>
slices of the volume have the most holes. When I use the same parameters for<br>
2D registration with slices from these images, the registered image is well<br>
aligned with the static image and there are no holes. Does any one know what<br>
could be the reason for the errors I'm getting with 3D registration? As a<br>
side note, I've had a similiar issue with affine transforms but instead of<br>
holes, the slices in the final image look as if they were cropped.<br>
<br>
Here's some more info about my parameters if it helps:<br>
Optimizer: advanced stochastic gradient descent with automatic parameter<br>
estimation<br>
Metric: advanced mattes mutual information<br>
Image sampler: random coordinate<br>
Number of iterations: 1000<br>
Number of resolutions: 10<br>
<br>
Thank you,<br>
Kristen<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.bigr.nl/pipermail/elastix/attachments/20110810/14f21a66/attachment.html" target="_blank">http://lists.bigr.nl/pipermail/elastix/attachments/20110810/14f21a66/attachment.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 10 Aug 2011 22:36:08 +0200<br>
From: Stefan Klein <<a href="mailto:s.klein@erasmusmc.nl" target="_blank">s.klein@erasmusmc.nl</a>><br>
Subject: Re: [elastix] 3D Bspline registration<br>
To: <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
Message-ID: <<a href="mailto:4E42EBB8.7000106@erasmusmc.nl" target="_blank">4E42EBB8.7000106@erasmusmc.nl</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Hi Kristen,<br>
<br>
This is probably because of motion in the z-direction. At the boundary<br>
of the image, the transformation is likely to map outside the moving<br>
image. In the result image this looks like a hole. To verify that this<br>
is the case, you could set the (DefaultValue ... ) entry in your<br>
parameter file to some specific value, for example 1000, and check if<br>
you get these values in the holes. The default value is the value that<br>
elastix uses for pixels that map outside the moving image domain (when<br>
generating the deformed moving image).<br>
<br>
If you think that motion out of plane is not appropriate in your<br>
application (due to the thick slice thickness for example), you could<br>
add the following line to your parameter file:<br>
(MovingImageDerivativeScales 1 1 0)<br>
<br>
Then elastix will not make any displacement in the z-direction.<br>
<br>
Kind regards,<br>
Stefan<br>
<br>
On 10-8-2011 16:29, Kristen Tgavalekos wrote:<br>
> Hi Elastix team,<br>
><br>
> I'm working on registering 3D images with a BSpline transform. They<br>
> are limited angle tomography images with different depths (a 10 slice<br>
> difference) and the resolution is 0.1 x 0.1 x 1 mm. My problem is that<br>
> the final registered image has blank holes in it. Typically the top<br>
> and bottom slices of the volume have the most holes. When I use the<br>
> same parameters for 2D registration with slices from these images, the<br>
> registered image is well aligned with the static image and there are<br>
> no holes. Does any one know what could be the reason for the errors<br>
> I'm getting with 3D registration? As a side note, I've had a similiar<br>
> issue with affine transforms but instead of holes, the slices in the<br>
> final image look as if they were cropped.<br>
><br>
> Here's some more info about my parameters if it helps:<br>
> Optimizer: advanced stochastic gradient descent with automatic<br>
> parameter estimation<br>
> Metric: advanced mattes mutual information<br>
> Image sampler: random coordinate<br>
> Number of iterations: 1000<br>
> Number of resolutions: 10<br>
><br>
> Thank you,<br>
> Kristen<br>
><br>
><br>
> _______________________________________________<br>
> elastix mailing list<br>
> <a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
> <a href="http://lists.bigr.nl/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/mailman/listinfo/elastix</a><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://lists.bigr.nl/pipermail/elastix/attachments/20110810/d3a81333/attachment-0001.htm" target="_blank">http://lists.bigr.nl/pipermail/elastix/attachments/20110810/d3a81333/attachment-0001.htm</a><br>
<br>
------------------------------<br>
<br>
_______________________________________________<br>
elastix mailing list<br>
<a href="mailto:elastix@bigr.nl" target="_blank">elastix@bigr.nl</a><br>
<a href="http://lists.bigr.nl/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/mailman/listinfo/elastix</a><br>
<br>
<br>
End of elastix Digest, Vol 26, Issue 5<br>
**************************************<br>
</blockquote></div><br>