[Insight-developers] Vector field to Jacobian determinant filter

[Luis Ibanez]

Hi Torsten,


Thanks for contributing your code to ITK.

The new filter has been committed to the CVS repository.
The test was added under Testing/Code/BasicFilters

and the Example was added under

          Insight/Examples/Registration

since it is quite relevant for deformable registration.

We simplified the author credits by removing the addresss
and your email. In general it is better not to put your
email since it get harvested by spammers over time.

The class should appear tomorrow in the Doxygen documentation.

If you have a chance please take a look at it and let us
know if there is anything that needs to be changed.


    Thanks a lot.



       Luis





-------------------------
Torsten Rohlfing wrote:

> 
> Good morning Luis!
> 
> I apologize for the delay. Attached is my code in a tar archive. I was 
> guessing that it would make more sense to send to you individually 
> rather than to the mailing list.
> 
> The archive has the following files.
> 
> itkDeformationFieldJacobianDeterminantFilter.h
> itkDeformationFieldJacobianDeterminantFilter.txx
> -- The actual filter class files
> 
> itkDeformationFieldJacobianDeterminantFilterTest.cxx
> -- A test (I am guessing here, since I haven't written one before). This 
> doesn't really test correct function yet, only lack of crashes and 
> exceptions. I'll write something more sophisticated once I am more 
> comfortable with the whole test thing.
> 
> DeformationFieldJacobian.cxx
> CMakeLists.txt
> -- A command line test tool; reads a deformation field, e.g., from 
> DemonsRegistration, and outputs a floating point image with the local 
> Jacobian determinant values of the def. field.
> 
> CREDITS
> -- my name and affiliation. Feel free to shorten as you see fit.
> 
> That's that. Let me know if there's anything else that needs to be done 
> right now.
> 
> There's even better news by the way: I have pretty much decided to 
> abandon independent code development on new projects and go with ITK 
> instead. So you should be getting a whole lot more code from me in the 
> future. :)
> 
> Best,
>  Torsten
> 
> Luis Ibanez wrote:
> 
>>
>> Hi Torsten,
>>
>> This is great news !,
>>
>> If we hurry up, we may be able to include your new filter for
>> the upcoming release 1.10.  We are freezing the CVS repository
>> on Friday January 21st.
>>
>>
>> Could you please send us your source code along with a short
>> line for the credits (author, department, institution).
>>
>> It will be most helpful and appreciated if you include a test
>> for the filter. Having a test will expedite the inclusion of
>> your filter in the repository.
>>
>> We also need any TABS to be removed from the source code
>> and replaced with two spaces.
>>
>> Regarding the alternative option of using this filter
>> as a measure of gradient magnitude, we usually prefer
>> to create a new filter with very specific functionality.
>> This makes documentation and testing much simpler.
>>
>>
>>
>> Please let us know if you have any questions.
>>
>>
>> Thanks a lot for contributing your code to ITK.
>>
>>
>>
>>     Regards,
>>
>>
>>        Luis
>>
>>
>>
>>
>> ---------------------
>> Torsten Rohlfing wrote:
>>
>>>
>>> Hi --
>>>
>>> it seems to me that it would be nice to have a filter that takes a 
>>> vector field (i.e., an n-dimensional image of n-dimensional vectors) 
>>> and produces an image of the local Jacobian determinants (for 
>>> deformation morphometry etc). In fact, this was suggested a few weeks 
>>> back on the itk-users mailing list I believe. As proposed there, I 
>>> have taken the VectorGradientMagnitudeImageFilter and modified it 
>>> accordingly. I would be happy to contribute this code if there is any 
>>> interest.
>>>
>>> Beyond that, I was wondering: one may consider the Jacobian 
>>> determinant as an alternative measure of gradient magnitude, albeit 
>>> one that is only well-defined for cases where the dimensionality of 
>>> the image equals that of the vectors at each pixel. So maybe it would 
>>> make more sense to include this computation as an alternative 
>>> operating mode in the VectorGradientMagnitudeImageFilter itself 
>>> rather than create a new class? The actual code is actually fairly 
>>> small, thanks to vxl's determinant computation.
>>>
>>> Best,
>>>  Torsten
>>>
>>> _______________________________________________
>>> Insight-developers mailing list
>>> Insight-developers at itk.org
>>> http://www.itk.org/mailman/listinfo/insight-developers
>>>
>>>
>>
>>
>>
>