<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'>
Greetings, <BR>
&nbsp;<BR>
I had posted a question a while back about intensity-based registration for low-contrast images. Thank you Luis for providing some pointers. You had requested I post example images so that you and the rest of the community might suggest what metrics are best suited to try to do automatic registration. I was waiting until I was&nbsp;able to create my own images before I posted them. The following links show two images, one which is rotated 5 degrees with respect to the other (rotation about the image center). <BR>
&nbsp;<BR>

<META content="Microsoft SafeHTML" name=Generator>
<STYLE>
.ExternalClass .EC_hmmessage P
{padding:0px;}
.ExternalClass body.EC_hmmessage
{font-size:10pt;font-family:Tahoma;}
</STYLE>
<A href="http://www.hotlinkfiles.com/view/full/1600023_zwyje" target=_blank>http://www.hotlinkfiles.com/view/full/1600023_zwyje</A><BR>&nbsp;<BR><A href="http://www.hotlinkfiles.com/view/full/1600024_ujxkr" target=_blank>http://www.hotlinkfiles.com/view/full/1600024_ujxkr</A><BR>
&nbsp;<BR>
These images are part of 3D-tomographic datasets. I chose the two corresponding slices that showed the most contrast to post as an example. All other images have even less contrast than shown here. <BR>
&nbsp;<BR>
Bearing in mind that I am new to ITK, I was wondering whether someone could suggest a metric that would be good for automatically registering two 3D&nbsp;datasets of such images. I would be happy to post the full image-sets if that would be more helpful. I am trying to determine the shifts in the three axes as well as rotations about each axis. <BR>
&nbsp;<BR>
Any help that anybody can provide will be much appreciated. <BR>
&nbsp;<BR>
Vik<BR><BR>&gt; <BR>&gt; Hi Vikren,<BR>&gt; <BR>&gt; <BR>&gt; 1) Yes, it is possible with ITK to compute the rigid transform that<BR>&gt; will map the coordinate system of your fixed image into the<BR>&gt; coordinate system of the moving image.<BR>&gt; <BR>&gt; <BR>&gt; 2) Could you post screeshots of your low contract images in a public<BR>&gt; web site ?<BR>&gt; <BR>&gt; Being able to see the images will help us a lot for recommending<BR>&gt; you an image metric to use.<BR>&gt; <BR>&gt; <BR>&gt; In the absence of being able to see your image, here are two<BR>&gt; suggestions that come to mind:<BR>&gt; <BR>&gt; A) Normalized Correlation metric<BR>&gt; <BR>&gt; B) Passing your images through edge detection filters<BR>&gt; and then using the results as the actual images to<BR>&gt; register. (you will find several edge-detection<BR>&gt; filters in ITK, you probably want to start with<BR>&gt; a simple gradient magnitude filter).<BR>&gt; <BR>&gt; <BR>&gt; C) Mutual information is still an option that you<BR>&gt; could reconsider, and for which you may have to<BR>&gt; adjust the parameters of number of bins and<BR>&gt; number of samples.<BR>&gt; <BR>&gt; <BR>&gt; <BR>&gt; In general, you should be careful discarding a particular technique<BR>&gt; based on having tried it once, with only a reduced set of parameters.<BR>&gt; <BR>&gt; It is not uncommon that for a fine-tunning exercise to be required<BR>&gt; before you get acceptable results with a given technique.<BR>&gt; <BR>&gt; <BR>&gt; <BR>&gt; Regards,<BR>&gt; <BR>&gt; <BR>&gt; Luis<BR>&gt; <BR>&gt; <BR>&gt; ---------------------<BR>&gt; Vikren Sarkar wrote:<BR>&gt; &gt; Greetings,<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; I am new to ITK and am trying to use the registration algorithms to<BR>&gt; &gt; register two 3D x-ray image datasets. The problem is that all the<BR>&gt; &gt; images have very low contrast and the intensity-based algorithms<BR>&gt; &gt; fail when used with the mean square metric (I also tried mutual<BR>&gt; &gt; information with no success). I need to perform a rigid registration<BR>&gt; &gt; and am hoping someone with more experience could suggest a metric<BR>&gt; &gt; which would be appropriate for use with low contrast images.<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; Alternatively, I am open to the idea of using control points which<BR>&gt; &gt; the user would input. I am most interested in the actual values of<BR>&gt; &gt; translation, rotation and scaling required for the registration than<BR>&gt; &gt; in the final registered datasets themselves. The reason is that we<BR>&gt; &gt; are investigating the use of those images in patient positioning<BR>&gt; &gt; prior to radiation delivery. Can someone please point me in the<BR>&gt; &gt; right direction as to whether this is possible using ITK? It is<BR>&gt; &gt; probably obvious that I am very new at this and I would appreciate<BR>&gt; &gt; any help and/or pointers. Thank you.<BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; <BR>&gt; &gt; Vik<BR>&gt; &gt; <BR><BR><BR><br /><hr />Time for vacation? WIN what you need. <a href='http://www.gowindowslive.com/summergiveaway/?ocid=tag_jlyhm' target='_new'>Enter Now!</a></body>
</html>