Sobel operator extended for 4D, legacy support for 3D, OK?

Algorithms
1

Please have a look at my pull request, which aims to improve our SobelOperator, and add ND support (including 4D):

I think it is ready, but there are a few things to consider:

  • The proposed new implementation of SobelOperator::GenerateCoefficients() is based on a function make_nd_sobel_kernels generated by AI, and presented by Hans (@hjmjohnson) at A ND Sobel implementation is needed. · Issue #5702 · InsightSoftwareConsortium/ITK · GitHub . I don’t fully understand this function, but it appears to produce the same numbers as the original implementation of itk::SobelOperator, for 2D. So I think it’s OK, but an extra review might be helpful.

  • The AI generated make_nd_sobel_kernels produces different coordinates for 3D than the original implementation of itk::SobelOperator (from ITK <= 5.4). The pull request proposes to add an option UseLegacyCoefficients, to get the coordinates of the original ITK implementation. I guess, legacy code should use the original coordinates for 3D, and new code should use the coordinates from the AI generated function, right?

  • The coordinates for 3D from the original implementation of itk::SobelOperator (from ITK <= 5.4) were based on “Irwin Sobel. An Isotropic 3x3x3 Volume Gradient Operator. Technical report, Hewlett-Packard Laboratories, April 1995.” I propose to remove the reference to this publication, because I cannot find it. But is this report still available somewhere?

1 Like
2

For what it’s worth, I just found a “Technical Reports” page of HP Labs from 1995, archived at

As well as at https://web.archive.org/web/20051219115308/http://hpl.hp.com/techreports/95/

Unfortunately, the technical report of Irwin Sobel, “An Isotropic 3x3x3 Volume Gradient Operator” (which is specifically about extending the operator to 3D) is not there!!!

This technical report was first mentioned in “itkSobelOperator.h” by Josh Cates (@joshcates) , June 24, 2003, commit FIX: Fixed 3D implementation. The Sobel operator now works correctly… · InsightSoftwareConsortium/ITK@b4cc496 · GitHub

Has anyone seen the technical report?

Update: I just read in CEASAR: A Smooth, Accurate and Robust Centerline Extraction Algorithm Ingmar Bitter et al., Stony Brook University:

I. Sobel. An Isotropic 3x3x3 Volume Gradient Operator. Hewlett-Packard’s Voxelator V-3 CD-ROM, Aug. 1996. Online Information at http://www.hp.com/info/voxelator (Paper is only on CD, not online)

Does anyone still have the CD?!?

3

I sent an email to Josh. Let’s see whether he responds, and whether he can find the report.

1 Like
4

Thanks @dzenanz

Looking at archive.org, it appears that HP’s Voxelator CD, containing the Sobel paper was distributed at the SIGGRAPH 1996 conference! Unfortunately I wasn’t there :smiley:

1 Like