<font size="2"><font face="verdana,sans-serif">Maybe you should ask this question on insight-developers mailing list? Perhaps there is some difference.</font></font><div><font size="2"><font face="verdana,sans-serif"><br></font></font></div>
<div><font size="2"><font face="verdana,sans-serif">HTH<br></font></font><br><div class="gmail_quote">On Sun, Apr 3, 2011 at 06:36, Mauro Maiorca <span dir="ltr"><<a href="mailto:mauromaiorca@gmail.com">mauromaiorca@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">nobody? :-(<br>
<br>
what I've done so far is transforming the input image with variable<br>
voxel spacing into one with constant voxel spacing (by using<br>
interpolation), then applying anisotropic diffusion, finally<br>
interpolating back the image to the original (variable) spacing. Well,<br>
it works, and I can apply any filter without worrying about the<br>
variable voxel spacing.<br>
<br>
However, I'm afraid this (re)interpolation trick might affect the<br>
final quality of the image, furthermore I'm just interested to few<br>
anisotropic diffusion filters. Thus, I'm still looking forward to find<br>
a way to modify the anisotropic diffusion filters already implemented<br>
in ITK in order to tackle the problem without using interpolation. Itk<br>
source code is well documented (thanks for that), but it's a quite big<br>
package, and it would be great to have directions regarding where in<br>
the code I should get my hands dirty on and what's the cleanest way to<br>
do it. I reckon itkFiniteDifferenceImageFilter.txx is one of the first<br>
target files I should inspect, and I guess at some stage I should play<br>
with the function<br>
void FiniteDifferenceImageFilter<TI,TO>::initializeFunctionCoefficients()<br>
I would like to treat the "voxel spacing mask" as image, something like that:<br>
typedef itk::Image<double, Dimension> variableVoxelSpacingMaskType;<br>
variableVoxelSpaceMaskType m_variableVoxelSpaceMask;<br>
I guess I have also to ask the user to insert the voxel space "mask"<br>
with a sort of SetVoxelSpaceMask() somewhere.<br>
<br>
Unfortunately, as you can see I have only fragmentary ideas of<br>
how/where I should modify the code, so any input would be appreciated.<br>
Since I haven't received any answer so far, please let me know whether<br>
I didn't explain my problem clearly enough, or maybe whether you think<br>
it's too complicate to be easily solved within the ITK toolkit.<br>
<br>
Regards,<br>
<font color="#888888">Mauro<br>
</font><div><div></div><div class="h5"><br>
<br>
On Mon, Mar 28, 2011 at 6:57 PM, Mauro Maiorca <<a href="mailto:mauromaiorca@gmail.com">mauromaiorca@gmail.com</a>> wrote:<br>
> Hello Itk-users,<br>
><br>
> I'm successfully using anisotropic diffusion filters on n-dimension volumes.<br>
> Unfortunately, many of my volumes don't have uniform voxel spacing. In<br>
> fact, in the simplest 3D scenario (n=3, x1=x, x2=y, x3=z) I can<br>
> summarise the voxel spacing within my images with a function f:R^3->R,<br>
> where df/dx=0, df/dy=0 df/dz=k (k constant). I was wondering whether<br>
> there is a way to apply diffusion parameters independently for each<br>
> dimension, or according to some sort of mask.<br>
> If negative, can you please point me to the quickest/cleanest way to<br>
> implement it within ITK (i.e. which files I need to modify, etc)? I'm<br>
> more then happy to share the code (InsightJournal) if other people<br>
> might be interested on this.<br>
><br>
> thanks a lot,<br>
> Mauro<br>
><br>
_____________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
</div></div></blockquote></div><br></div>