
<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hello,<div><br></div><div>I don't think there is going to be anything turnkey to do what you want.</div><div><br></div><div>I'd looking to writing a custom filter for this specialized task. You can either base it off the the HessianRecursiveGaussianImageFilter [1] , which utilizes filters per-component and adaptors.</div><div><br></div><div>Or do a discrite approach along the lines of something I have done here [2]. </div><div><br></div><div>Hope that helps,</div><div>Brad</div><div><br></div><div>[1] <a href="https://github.com/InsightSoftwareConsortium/ITK/blob/1ce0e5dd3b2d9c7a7a9353dfb8b728ada8d71320/Modules/Filtering/ImageFeature/include/itkHessianRecursiveGaussianImageFilter.hxx#L187-L368">https://github.com/InsightSoftwareConsortium/ITK/blob/1ce0e5dd3b2d9c7a7a9353dfb8b728ada8d71320/Modules/Filtering/ImageFeature/include/itkHessianRecursiveGaussianImageFilter.hxx#L187-L368</a></div><div>[2] <a href="https://github.com/blowekamp/itkLocalDiscreteHessian/blob/30aa775494c20eaf9e085781ae233e0511773758/include/itkHessianImageFilter.hxx#L132-L174">https://github.com/blowekamp/itkLocalDiscreteHessian/blob/30aa775494c20eaf9e085781ae233e0511773758/include/itkHessianImageFilter.hxx#L132-L174</a></div><div><br><div><div>On Mar 27, 2015, at 10:16 AM, Thomas Seidel <<a href="mailto:seidel@cvrti.utah.edu">seidel@cvrti.utah.edu</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; font-size: 14px; font-family: Calibri, sans-serif;"><div>Hi all,</div><div><br></div><div>I would like to calculate the hessian matrix and get an output of type SymmetricSecondRankTensor pixel from a gradient image.</div><div><br></div><div>The irk HessianRecursiveGaussianImageFilter uses as input an intensity image.</div><div>However, I would like to implement a method for vessel detection based on the following algorithm:</div><div><br></div><div>1. GradientImageFilter</div><div>2. GradientVectorFlowImageFilter</div><div>3. get the hessian Matrix from 2. (this is my problem)</div><div>4. Hessian3DToVesselnessImageFilter</div><div><br></div><div><br></div><div>How can I get from the vector field in 2. to the hessian matrix I need for 3.?</div><div><br></div><div>Thanks for your help,</div><div>Thomas</div></div>

_____________________________________<br>Powered by <a href="http://www.kitware.com">www.kitware.com</a><br><br>Visit other Kitware open-source projects at<br><a href="http://www.kitware.com/opensource/opensource.html">http://www.kitware.com/opensource/opensource.html</a><br><br>Kitware offers ITK Training Courses, for more information visit:<br>http://www.kitware.com/products/protraining.php<br><br>Please keep messages on-topic and check the ITK FAQ at:<br>http://www.itk.org/Wiki/ITK_FAQ<br><br>Follow this link to subscribe/unsubscribe:<br>http://public.kitware.com/mailman/listinfo/insight-users<br></blockquote></div><br></div></body></html>