[Insight-users] Fusing two volumes, how-to?

[Robert Manzke]

Hi Siavash, Christopher, all

I have progressed a bit and am able to add two identical images. However,
if I take two images which do not "occupy the same physical space", e.g.
due to different voxel spacing etc. it fails.

I think I need to
- calculate the maximum bounds (/dimensions) of the fused image containing
the two individual images
- resample both images to the combined one which is large enough to contain
both other images

Is there any advice how to go from here?

I have attached my code so far...

This task must have been solved many times since if you have registered two
images using any method, one may want to represent those two images in one
image (space)...?

Thanks much,
Bob


On Tue, Sep 17, 2013 at 11:05 PM, Siavash Khallaghi <siavashk at ece.ubc.ca>wrote:

>  Hi Bob,
>
> You can do this using itk filters in the following order:
>
> 1. Use an itk::HistogramMatchingImageFilter to do histogram equalization.
> 2. Use an itk::ResampleImageFilter to map one volume spatially to the
> other.
> 3. Use an itk::AddImageFilter to add the result of Step 2.
> 4. Use an itk::MultiplyByConstantImageFilter on the result of Step 3 to
> create the average (there may be a filter that combines Steps 3 and 4).
>
> Siavash
>
>
> On 17/09/2013 10:46 AM, Christopher Mullins wrote:
>
> Check out the TubeTK MergeAdjacentImages application [1].  Looks like this
> will do several of the tasks you're looking for, and it's built on top of
> ITK.
>
>  [1]
> https://github.com/TubeTK/TubeTK/tree/master/Applications/MergeAdjacentImages
>
>
> On Tue, Sep 17, 2013 at 1:41 PM, Robert Manzke <paradenker at gmail.com>wrote:
>
>> Hi,
>>
>>  I want to create one single volume from two .mhd files.
>> I also have a 4x4 transformation matrix which maps one volume into
>> another.
>>
>>  So basically, I need an interpolator and a merge function, i.e. average
>> volume 1 and 2, to create a single volume.
>> The single volume will have a larger size than the two individual ones,
>> if the transform is not identity.
>>
>>  It can happen that both images have different gray level histograms, so
>> histogram matching before fusion is also desirable.
>>
>>  Can anyone lead me into the right direction how to do this with itk?
>>  I would volunteer to post my final code for a small how-to, once it
>> works.
>>
>>  Thanks,
>>  Bob
>>
>> _____________________________________
>> Powered by www.kitware.com
>>
>> Visit other Kitware open-source projects at
>> http://www.kitware.com/opensource/opensource.html
>>
>> Kitware offers ITK Training Courses, for more information visit:
>> http://www.kitware.com/products/protraining.php
>>
>> Please keep messages on-topic and check the ITK FAQ at:
>> http://www.itk.org/Wiki/ITK_FAQ
>>
>> Follow this link to subscribe/unsubscribe:
>> http://www.itk.org/mailman/listinfo/insight-users
>>
>>
>
>
>  --
> Christopher Mullins
> R&D Engineer
> Kitware Inc.,
> 919.869.8871
>
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects athttp://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:http://www.kitware.com/products/protraining.php
>
> Please keep messages on-topic and check the ITK FAQ at:http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:http://www.itk.org/mailman/listinfo/insight-users
>
>
>
> --
> PhD Candidate
> Robotics and Control Laboratory, Room 3090
> Electrical and Computer Engineering, UBC
> 2332 Main Mall, Vancouver, BC, Canada, V6T 1Z4www.ece.ubc.ca/~siavashk/
> Tel: (604)822-9215
>
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.php
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.itk.org/pipermail/insight-users/attachments/20130924/ee7462b2/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MetaImageFusionTool.cpp
Type: text/x-c++src
Size: 2485 bytes
Desc: not available
URL: <http://www.itk.org/pipermail/insight-users/attachments/20130924/ee7462b2/attachment.cpp>