<html><body bgcolor="#FFFFFF"><div>Well, of course I meant rigid pre-registration ...<br><br>Am 17.02.2009 um 12:24 schrieb Andreas Schuh <<a href="mailto:andreas.schuh.84@googlemail.com">andreas.schuh.84@googlemail.com</a>>:<br><br></div><div></div><blockquote type="cite"><div><div>Hi!</div><div><br></div><div>Well, that's more or less the same issue, when I was thinking about this overlap thing. Well, I got a solver from a fellow that computed an overlap cube of two other cubes that has minimal extent/volume, where the direction of this cube is the one of the coordinate system (this solver is based on a simplex approach I think). I've implemented an ITK filter that uses this solver to get an ImageRegion and Direction matrix (orientation of the overlap cube), where I set the direction to the direction of the fixed image. These can then be passes to ResampleImageFilter to resample both images, fixed and moving, within the overlap. Unfortunately, my current implementation doesn't work. However, I didn't focus on this any more as for me, it was enough to just resample both images on a grid with the same direction and extent as the fixed image. This is because non-parametric registratration requires non-rigid pre-registration in most cases. So the misalignment of the image regions themselves isn't that much (and the borders seldomly contain something of interest). After pre-registration you could also just resample the images on a grid with the direction of the fixed image but the extent large enough to cover at least both fixed and moving images.</div><div><br></div><div>I think I just have an idea how you can get the mask you want.</div><div><br></div><div>Resample the moving image on a grid that has the same direction and size as the fixed image. So, use ResampleImageFilter and set the fixed image as reference image. Set the outside value appropriately to reflect the "masked out meaning".</div><div><span class="Apple-style-span" style="-webkit-composition-fill-color: rgba(175, 192, 227, 0.231373); -webkit-composition-frame-color: rgba(77, 128, 180, 0.231373); ">Depending on the mask format your metric assumes, you may just use the resulting image as mask or apply a threshold to get a binary mask image. </span><br></div><div><br></div><div>--</div><div>regards </div><div>Andreas</div><div> <br>Am 17.02.2009 um 09:43 schrieb Erik Türke <<a href="mailto:erik_tuerke@hotmail.com"><a href="mailto:erik_tuerke@hotmail.com">erik_tuerke@hotmail.com</a></a>>:<br><br></div><div></div><blockquote type="cite"><div>
Hi!<br><br>Thanks for the response.<br><br>I´ll try to explain yout the issue. <br>It is all about a regsitration problem.<br>I have two images, a fixed and a moving image. Lets consider the fixed image is of dimension [192 130 192] and the moving image is of dimension [192 7 192].<br>The metric will give you a sample exception error because too many samples are outside the moving image. <br>Now, what i want to do is to define a mask, which contains the overlap of the two images and pass it to the metric->SetFixedImageMask() method, which expects a ImageMaskSpatialObject.<br><br>So now my problem is how to compute an overlap(image, spatialobject) containing the fixed and the moving images, or the fixed and the moving spatial objects, respectively. <br>And i assume, i am capable of creating an ImageMaskSpatialObject with the help of the obtained overlap image or object and pass it to the metric.<br><br><br>thanks and best regards!<br><br><br>> From: <a href="mailto:andreas.schuh.84@googlemail.com"></a><a href="mailto:andreas.schuh.84@googlemail.com"><a href="mailto:andreas.schuh.84@googlemail.com">andreas.schuh.84@googlemail.com</a></a><br>> To: <a href="mailto:luis.ibanez@kitware.com"></a><a href="mailto:luis.ibanez@kitware.com"><a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a></a><br>> Subject: Re: [Insight-users] Getting overlap region<br>> Date: Mon, 16 Feb 2009 12:11:05 +0100<br>> CC: <a href="mailto:erik_tuerke@hotmail.com"><a href="mailto:erik_tuerke@hotmail.com">erik_tuerke@hotmail.com</a></a>; <a href="mailto:insight-users@itk.org"></a><a href="mailto:insight-users@itk.org"><a href="mailto:insight-users@itk.org">insight-users@itk.org</a></a><br>> <br>> Hi,<br>> <br>> the difference of determining the overlap region of two single <br>> ImageRegions and of getting the overlap of two oriented images or <br>> spatial objects is, that in the former case the regions are not <br>> oriented (considered having the same orientation) whereas in the <br>> latter cases you may have to be corncerned with the orientation of <br>> your regions, too. In this case ImageRegion::Crop won't help if you <br>> would like to have the oriented overlap region of minimal volume. <br>> That's a more difficult task, where you will have to solve some <br>> optimization problem probably by some kind of simplex algorithm.<br>> <br>> Otherwise, as you Erik are interested in spatial objects and I assume <br>> that the bounding boxes of these (SpatialObjects::GetBoundingBox) are <br>> all oriented along the axes of the coordinate system, you have a <br>> simpler task which can be solved by using ImageRegion::Crop.<br>> <br>> --<br>> regards<br>> Andreas<br>> <br>> Am 15.02.2009 um 22:52 schrieb Luis Ibanez <<a href="mailto:luis.ibanez@kitware.com"><a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a></a>>:<br>> <br>> ><br>> > Hi Erik,<br>> ><br>> > The easiest method is among ImageRegions.<br>> ><br>> > In that case, simply use:<br>> ><br>> > itk::ImageRegion::Crop()<br>> ><br>> > If you have two image regions<br>> ><br>> > itk::ImageRegion regionA;<br>> > itk::ImageRegion regionB;<br>> ><br>> > Then you can find the overlap by doing<br>> ><br>> > regionA.Crop( regionB )<br>> ><br>> > The common region will be stored in the<br>> > variable regionA.<br>> ><br>> ><br>> > Regards,<br>> ><br>> ><br>> > Luis<br>> ><br>> > ------------------<br>> > Erik Türke wrote:<br>> >> Hi!<br>> >> Thanks for your response.<br>> >> Actually it doesnt matter, if i am doing this with an image, an <br>> >> image region, or spatial object, since you can convert these types <br>> >> amoung themselves. I guess ;-)<br>> >> I think, i would prefere the easiest method.<br>> >> But if there is no difference between the effort of these methods, <br>> >> i would take the SpatialObject method.<br>> >> Thanks!<br>> >> best regards!<br>> >> > Date: Thu, 12 Feb 2009 11:16:21 -0500<br>> >> > From: <a href="mailto:luis.ibanez@kitware.com"></a><a href="mailto:luis.ibanez@kitware.com"><a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a></a><br>> >> > To: <a href="mailto:erik_tuerke@hotmail.com"></a><a href="mailto:erik_tuerke@hotmail.com"><a href="mailto:erik_tuerke@hotmail.com">erik_tuerke@hotmail.com</a></a><br>> >> > CC: <a href="mailto:insight-users@itk.org"></a><a href="mailto:insight-users@itk.org"><a href="mailto:insight-users@itk.org">insight-users@itk.org</a></a><br>> >> > Subject: Re: [Insight-users] Getting overlap region<br>> >> ><br>> >> ><br>> >> > Hi Erik,<br>> >> ><br>> >> > You question is too general.<br>> >> ><br>> >> > Could you elaborate on it please ?<br>> >> ><br>> >> > Do you :<br>> >> ><br>> >> > a) Have two images and want to create a third<br>> >> > image containing the overlap of the two images ?<br>> >> ><br>> >> > OR<br>> >> ><br>> >> > b) Have two itkImageRegions and want to compute<br>> >> > the itkImageRegion that represents the overlap<br>> >> > of the two ?<br>> >> ><br>> >> > OF<br>> >> ><br>> >> > c) Have two itkSpatialObjects and you are trying to<br>> >> > compute the itkSpatialObject representing their<br>> >> > overlap ?<br>> >> ><br>> >> ><br>> >> ><br>> >> > Thanks<br>> >> ><br>> >> ><br>> >> > Luis<br>> >> ><br>> >> ><br>> >> > ------------------<br>> >> > Erik Türke wrote:<br>> >> > > Hi!<br>> >> > ><br>> >> > > I am asking myself if there is any possibility to get the <br>> >> region of an<br>> >> > > overlap of 2 images or a new object containing the overlap of 2 <br>> >> other<br>> >> > > objects, respectively.<br>> >> > ><br>> >> > > Thanks fr your help!<br>> >> > ><br>> >> > > Best regards!<br>> >> > ><br>> >> > ><br>> >> > ><br>> >> > > <br>> >> --- <br>> >> ---------------------------------------------------------------------<br>> >> > > MSN Mobile: Die neuesten Infos aus der Welt der Stars und <br>> >> Sternchen zum<br>> >> > > Mitnehmen!<br>> >> > > <<a href="http://redirect.gimas.net/?cat=hmtl&n=M0809MSNMobile&d=http://info.mobile.de.msn.com/pc/default.aspx"><a href="http://redirect.gimas.net/?cat=hmtl&n=M0809MSNMobile&d=http://info.mobile.de.msn.com/pc/default.aspx">http://redirect.gimas.net/?cat=hmtl&n=M0809MSNMobile&d=http://info.mobile.de.msn.com/pc/default.aspx</a></a> <br>> >> ><br>> >> > ><br>> >> > ><br>> >> > > <br>> >> --- <br>> >> ---------------------------------------------------------------------<br>> >> > ><br>> >> > > _____________________________________<br>> >> > > Powered by <a href="http://www.kitware.com"></a><a href="http://www.kitware.com"><a href="http://www.kitware.com">www.kitware.com</a></a><br>> >> > ><br>> >> > > Visit other Kitware open-source projects at<br>> >> > > <a href="http://www.kitware.com/opensource/opensource.html"></a><a href="http://www.kitware.com/opensource/opensource.html"><a href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a></a><br>> >> > ><br>> >> > > Please keep messages on-topic and check the ITK FAQ at: <a href="http://www.itk.org/Wiki/ITK_FAQ"></a><a href="http://www.itk.org/Wiki/ITK_FAQ"><a href="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</a></a><br>> >> > ><br>> >> > > Follow this link to subscribe/unsubscribe:<br>> >> > > <a href="http://www.itk.org/mailman/listinfo/insight-users"></a><a href="http://www.itk.org/mailman/listinfo/insight-users"><a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a></a><br>> >> --- <br>> >> ---------------------------------------------------------------------<br>> >> Messenger Online Treff: Spontan chatten, über Hobbies reden, sofor <br>> >> t Spass haben? Jetzt klicken! <<a href="http://redirect.gimas.net/?cat=hmtl&n=M0809MOT&d=http://messenger.live.de/community/messenger-online-treff.aspx"><a href="http://redirect.gimas.net/?cat=hmtl&n=M0809MOT&d=http://messenger.live.de/community/messenger-online-treff.aspx">http://redirect.gimas.net/?cat=hmtl&n=M0809MOT&d=http://messenger.live.de/community/messenger-online-treff.aspx</a></a> <br>> >> ><br>> > _____________________________________<br>> > Powered by <a href="http://www.kitware.com"></a><a href="http://www.kitware.com"><a href="http://www.kitware.com">www.kitware.com</a></a><br>> ><br>> > Visit other Kitware open-source projects at<br>> > <a href="http://www.kitware.com/opensource/opensource.html"></a><a href="http://www.kitware.com/opensource/opensource.html"><a href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a></a><br>> ><br>> > Please keep messages on-topic and check the ITK FAQ at: <a href="http://www.itk.org/Wiki/ITK_FAQ"></a><a href="http://www.itk.org/Wiki/ITK_FAQ"><a href="http://www.itk.org/Wiki/ITK_FAQ">http://www.itk.org/Wiki/ITK_FAQ</a></a><br>> ><br>> > Follow this link to subscribe/unsubscribe:<br>> > <a href="http://www.itk.org/mailman/listinfo/insight-users"></a><a href="http://www.itk.org/mailman/listinfo/insight-users"><a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a></a><br><br><hr>Speichern Sie Ihre Daten von überall! <a href="http://redirect.gimas.net/?cat=hmtl&n=M0809Sky&d=http://skydrive.live.com/?mkt=de-de" target="_new">Kostenlos: Die Festplatte im Netz mit 5 GB Speicher!</a>
</div></blockquote></div></blockquote></body></html>