<div dir="ltr">Hi Dan,<div><br></div><div>To add to question 3:</div><div>For the lung data we always used a 4D mask, which is the concatenation of the 3D masks for the separate time points.</div><div><br></div><div>Regards,</div><div>Coert</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Feb 13, 2015 at 2:01 PM, Floris Berendsen <span dir="ltr"><<a href="mailto:floris@isi.uu.nl" target="_blank">floris@isi.uu.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Dan,<br><br>That is an interesting application. I was not involved in the works you mention, however I can give my view based on experiences with these methods and the elastix framework.<br><br>Answer 1): <br>Unfortunately, a set of 3d images must be provided as a 4D volume with its last dimension the image number or time index. In the current framework of elastix the groupwise registration is sort of hacked in, therefore it has these restrictions.<br><br>Asnwer 2): <br>I think my answer is threefold.<br>2a): I would assume that more (independent) image data can provide a better registration (estimate) in general. A higher SNR is a bonus as well.<br>2b): In the current groupwise registration framework the resulting transformation is defined to point from a floating 'spatial average' image domain to each time frame domain. Since you would like to model all deformations in the space of the TLC image you would need to chain the resulting transformations to each time frame with the inverse transformation of the TLC. This can be done with some scripting.<br>As a downside of including the TLC image in the time stack you  probably should not use any time domain smoothness for the registration, since that doesn't make sense for such a hybrid stack.<br>2c): Ultimately, you might want to adapt the groupwise framework such that you have a the TLC image as fixed image and the time stack as moving image. In that way you have the TLC image as a reference domain for modelling and the image similarity metric is still calculated in a groupwise fashion. You would need to write an adapted metric and possibly transformation for that, but the current framework allows this, I think.<br><br>Answer 3): <br>Not sure about the details in the paper. I would say that the lung mask needs to be 4D (and possibly different for each 3d sub volume), such that the sampler knows what the valid samples are.<br><br>Answer 4): <br>External tools/labour are needed to generate these label images. Delmon et al. use "Automated Segmentation of a Motion Mask to Preserve Sliding Motion in Deformable Registration of Thoracic CT" for that. I worked on an approach similar to that of Delmon et al. My souce code is not yet part of elastix, but the method is published so far is as:<br><br>Registration of organs with sliding interfaces and changing topologies<br><a href="http://spie.org/Publications/Proceedings/Paper/10.1117/12.2043447" target="_blank">http://spie.org/Publications/Proceedings/Paper/10.1117/12.2043447</a><br><br>Combining sliding organs with a groupwise approach has come to my mind as well. I think it can be done, but requires bigger adaptations of the code and might be limited to certain setups. Both our slidng organ methods require a label image that is defined in the fixed domain. In principle this conflicts with the floating spatial domain, but 2c) would solve this.<br><br>Best,<br>Floris<br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Thu, Feb 12, 2015 at 9:42 PM, Einstein, Daniel R <span dir="ltr"><<a href="mailto:Daniel.Einstein@pnnl.gov" target="_blank">Daniel.Einstein@pnnl.gov</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">



<div style="word-wrap:break-word">
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
Hello,</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
I am interested in applying a number of ElastiX 4D lung registration approaches to 11 dynamic rat lung images that are described in:</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<a href="http://www.ncbi.nlm.nih.gov/pubmed/22087338" target="_blank">http://www.ncbi.nlm.nih.gov/pubmed/22087338</a></div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
and</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<a href="http://www.ncbi.nlm.nih.gov/pubmed/23799057" target="_blank">http://www.ncbi.nlm.nih.gov/pubmed/23799057</a></div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
The parameter file (Par0012) available at <a href="http://elastix.bigr.nl/wiki/index.php/Par0012" target="_blank">http://elastix.bigr.nl/wiki/index.php/Par0012</a> seems to be what I am after. </div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
Question 1): do the images input to the elastic command line need to be truly 4D, or can they be multiple 3D images?</div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
<div><font face="Calibri,sans-serif">elastix -f <dynamic nD+t image> -m <dynamic nD+t image> -p <par filename> -out <output dir></font></div>
<div><font face="Calibri,sans-serif"><br>
</font></div>
<div><font face="Calibri,sans-serif">Question 2): We have a higher-resolution, higher SNR static TLC image from the same animal upon which we would like to base a model geometry. Assuming both sets of images (4D + TLC) are up/down sampled
 to be equal, do you think there is any benefit to including the TLC image in the groupwise registration?</font></div>
<div><font face="Calibri,sans-serif"><br>
</font></div>
<div><font face="Calibri,sans-serif">Question 3): The text of the Metz et al 2011 paper mentions that</font> "the proposed method, registration was always performed on the complete 4D image and using a lung mask". The command line above
 does not call for a mask. Were different settings used in the Metz et al paper and was a single mask used for all images or was a 4D mask image used?</div>
<div><br>
</div>
<div>Question 4): Later work with ElastiX described in Delmon et al 2013, describes the use of direction dependent B-splines for relative motion between the lung and the thorax (<a href="http://elastix.bigr.nl/wiki/index.php/Par0016" target="_blank">http://elastix.bigr.nl/wiki/index.php/Par0016</a>).
 This exercise requires a label image, that appears not to be readily available. Has any attempt been made to combine the methods in Par0012 and Par0016?</div>
<div><br>
</div>
<div>Thank you for your patience as I get my feet wet with your code.</div>
<div><br>
</div>
<div>Kind Regards,</div>
<div>Dan</div>
<div><br>
</div>
<div>
<div>
<div>___________________________________________</div>
<div style="font-family:Calibri,sans-serif"><font color="#000000"><font face="Calibri"><span style="font-size:14px">Daniel R Einstein, PhD</span></font></font></div>
<div style="font-family:Calibri,sans-serif">Computational Biology & Bioinformatics</div>
<div style="font-family:Calibri,sans-serif">Pacific Northwest National Laboratory</div>
<div style="font-family:Calibri,sans-serif"><a href="http://www.pnnl.gov" target="_blank">www.pnnl.gov</a></div>
<div style="font-family:Calibri,sans-serif"><br>
</div>
<div style="font-family:Calibri,sans-serif">and</div>
<div style="font-family:Calibri,sans-serif"><br>
</div>
<div style="font-family:Calibri,sans-serif">Department of Mechanical Engineering</div>
<div style="font-family:Calibri,sans-serif">University of Washington, Seattle</div>
<div style="font-family:Calibri,sans-serif"><a href="https://www.me.washington.edu" target="_blank">https://www.me.washington.edu</a></div>
<div style="font-family:Calibri,sans-serif"><br>
</div>
<div style="font-family:Calibri,sans-serif"><br>
</div>
</div>
</div>
<div><font face="Calibri,sans-serif"><br>
</font></div>
<div><font face="Calibri,sans-serif"><br>
</font></div>
<div style="color:rgb(0,0,0);font-family:Calibri,sans-serif;font-size:14px">
<br>
</div>
</div>

<br></div></div>_______________________________________________<br>
Elastix mailing list<br>
<a href="mailto:Elastix@bigr.nl" target="_blank">Elastix@bigr.nl</a><br>
<a href="http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix</a><br>
<br></blockquote></div><br></div>
<br>_______________________________________________<br>
Elastix mailing list<br>
<a href="mailto:Elastix@bigr.nl">Elastix@bigr.nl</a><br>
<a href="http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix" target="_blank">http://lists.bigr.nl/cgi-bin/mailman/listinfo/elastix</a><br>
<br></blockquote></div><br></div>