<div dir="ltr">Thanks Matt, helped a lot.<div><br></div><div>I know there are somethings I could find inside the ITK code but I'd like to ask you a few questions:</div><div>Are the gaussians normalized during this calculus? </div><div>Can I change the maximum value? </div><div>The variances help me adjust the radius but can I actually limit it? </div><div>is the variance related to the actual radius in voxels? say: a variance of sigma would give me a 99.7% of the values under a radius of 3*sigma?</div><div>what if I use the variance sigma but I want to limit it to a 2*sigma raidus? </div><div>can I change the mean?</div><div><br></div><div>Thanks!</div><div><br></div><div>Leslie</div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-04-09 16:38 GMT+02:00 Matt McCormick <span dir="ltr"><<a href="mailto:matt.mccormick@kitware.com" target="_blank">matt.mccormick@kitware.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Thu, Apr 9, 2015 at 9:17 AM, Leslie Solorzano <<a href="mailto:lesolorzanov@gmail.com">lesolorzanov@gmail.com</a>> wrote:<br>
> Hello world.<br>
<br>
:-)<br>
<span class=""><br>
<br>
> I want to use a gaussianSpatialObject and a<br>
> SpatialObjectToImageStatisticsCalculator to filter an image and make a<br>
> gaussian convolution. Why not using a DiscreteGaussianImageFilter? easy, I<br>
> need a gaussian with different radii in each dimension and I need to be able<br>
> to set the size, the maximum value and the sigma, which only a<br>
> gaussianSpatialObject can offer me.<br>
<br>
</span><span class="">> is there anyway I can use DiscreteGaussianImageFilter with different radii?<br>
> can I set the maximum value?<br>
<br>
</span>Yes, DiscreteGaussianImageFilter will give much better performance.<br>
<br>
To set a different Gaussian width in each direction, something like:<br>
<br>
  typedef itk::DiscreteGaussianImageFilter< ImageType, ImageType><br>
GaussianFilterType;<br>
  typedef GaussianFilterType::ArrayType VariancesType;<br>
  VariancesType variances;<br>
  variances[0] = 1.1;<br>
  variances[1] = 0.5;<br>
  variances[2] = 7.7;<br>
  GaussianFilterType::Pointer gaussianFilter = GaussianFilterType::New();<br>
  gaussianFilter->SetVariance( variances );<br>
<br>
would set different variances in each direction for a 3D image.<br>
<br>
HTH,<br>
Matt<br>
</blockquote></div><br></div>