A method for estimation of a surface from a binary volume.
- This filter implements a surface-fitting method for estimation of a surface from a binary volume. This process can be used to reduce aliasing artifacts which result in visualization of binary partitioned surfaces.
- The binary volume (filter input) is used as a set of constraints in an iterative relaxation process of an estimated ND surface. The surface is described implicitly as the zero level set of a volume \( \phi \) and allowed to deform under curvature flow. A set of constraints is imposed on this movement as follows:
\[ u_{i,j,k}^{n+1} = \left\{ \begin{array}{ll} \mbox{max} (u_{i,j,k}^{n} + \Delta t H_{i,j,k}^{n}, 0) & \mbox{\f$B_{i,j,k} = 1\f$} \\ \mbox{min} (u_{i,j,k}^{n} + \Delta t H_{i,j,k}^{n}, 0) & \mbox{\f$B_{i,j,k} = -1\f$} \end{array}\right. \]
- where \( u_{i,j,k}^{n} \) is the value of \( \phi \) at discrete index \( (i,j,k) \) and iteration \( n \) , \( H \) is the gradient magnitude times mean curvature of \( \phi \) , and \( B \) is the binary input volume, with 1 denoting an inside pixel and -1 denoting an outside pixel.
- NOTES
- This implementation uses a sparse field level set solver instead of the narrow band implementation described in the reference below, which may introduce some differences in how fast and how accurately (in terms of RMS error) the solution converges.
- REFERENCES
- Whitaker, Ross. "Reducing Aliasing Artifacts In Iso-Surfaces of Binary
Volumes" IEEE Volume Visualization and Graphics Symposium, October 2000, pp.23-32.
- PARAMETERS
- The MaximumRMSChange parameter is used to determine when the solution has converged. A lower value will result in a tighter-fitting solution, but will require more computations. Too low a value could put the solver into an infinite loop. Values should always be less than 1.0. A value of 0.07 is a good starting estimate.
- The MaximumIterations parameter can be used to halt the solution after a specified number of iterations.
- INPUT
- The input is an N-dimensional image of any type. It is assumed to be a binary image. The filter will use an isosurface value that is halfway between the min and max values in the image. A signed data type is not necessary for the input.
- OUTPUT
- The filter will output a level set image of real, signed values. The zero crossings of this (N-dimensional) image represent the position of the isosurface value of interest. Values outside the zero level set are negative and values inside the zero level set are positive values.
- IMPORTANT!
- The output image type you use to instantiate this filter should be a real valued scalar type. In other words: doubles or floats.
- USING THIS FILTER
- The filter is relatively straightforward to use. Tests and examples exist to illustrate. The important thing is to understand the input and output types so you can properly interpret your results.
- In the common case, the only parameter that will need to be set is the MaximumRMSChange parameter, which determines when the solver halts.
- See also
- itk::simple::AntiAliasBinary for the procedural interface
-
itk::AntiAliasBinaryImageFilter for the Doxygen on the original ITK class.
Definition at line 91 of file sitkAntiAliasBinaryImageFilter.h.
|
| static bool | GetGlobalDefaultDebug () |
| |
| static void | GlobalDefaultDebugOff () |
| |
| static void | GlobalDefaultDebugOn () |
| |
| static void | SetGlobalDefaultDebug (bool debugFlag) |
| |
| static void | GlobalWarningDisplayOn () |
| |
| static void | GlobalWarningDisplayOff () |
| |
| static void | SetGlobalWarningDisplay (bool flag) |
| |
| static bool | GetGlobalWarningDisplay () |
| |
| static double | GetGlobalDefaultCoordinateTolerance () |
| | Access the global tolerance to determine congruent spaces. More...
|
| |
| static void | SetGlobalDefaultCoordinateTolerance (double) |
| | Access the global tolerance to determine congruent spaces. More...
|
| |
| static double | GetGlobalDefaultDirectionTolerance () |
| | Access the global tolerance to determine congruent spaces. More...
|
| |
| static void | SetGlobalDefaultDirectionTolerance (double) |
| | Access the global tolerance to determine congruent spaces. More...
|
| |
| static bool | SetGlobalDefaultThreader (const std::string &threader) |
| | Set/Get the default threader used for process objects. More...
|
| |
| static std::string | GetGlobalDefaultThreader () |
| | Set/Get the default threader used for process objects. More...
|
| |
| static void | SetGlobalDefaultNumberOfThreads (unsigned int n) |
| |
| static unsigned int | GetGlobalDefaultNumberOfThreads () |
| | Set/Get the default threader used for process objects. More...
|
| |
| void | CheckImageMatchingDimension (const Image &image1, const Image &image2, const std::string &image2Name) |
| |
| void | CheckImageMatchingPixelType (const Image &image1, const Image &image2, const std::string &image2Name) |
| |
| void | CheckImageMatchingSize (const Image &image1, const Image &image2, const std::string &image2Name) |
| |
| virtual unsigned long | AddITKObserver (const itk::EventObject &, itk::Command *) |
| |
| virtual itk::ProcessObject * | GetActiveProcess () |
| |
| virtual void | OnActiveProcessDelete () |
| |
| virtual void | onCommandDelete (const itk::simple::Command *cmd) noexcept |
| |
| virtual void | PreUpdate (itk::ProcessObject *p) |
| |
| virtual void | RemoveITKObserver (EventCommand &e) |
| |
| | NonCopyable ()=default |
| |
| | NonCopyable (const NonCopyable &)=delete |
| |
| NonCopyable & | operator= (const NonCopyable &)=delete |
| |
| template<class TImageType > |
| static void | FixNonZeroIndex (TImageType *img) |
| |
| template<class TImageType > |
| static TImageType::ConstPointer | CastImageToITK (const Image &img) |
| |
| template<class TPixelType , unsigned int VImageDimension, unsigned int VLength, template< typename, unsigned int > class TVector> |
| static Image | CastITKToImage (itk::Image< TVector< TPixelType, VLength >, VImageDimension > *img) |
| |
| template<unsigned int VImageDimension, unsigned int VLength, template< unsigned int > class TVector> |
| static Image | CastITKToImage (itk::Image< TVector< VLength >, VImageDimension > *img) |
| |
| template<class TImageType > |
| static Image | CastITKToImage (TImageType *img) |
| |
| static const itk::EventObject & | GetITKEventObject (EventEnum e) |
| |
| template<typename T > |
| static std::ostream & | ToStringHelper (std::ostream &os, const T &v) |
| |
| static std::ostream & | ToStringHelper (std::ostream &os, const char &v) |
| |
| static std::ostream & | ToStringHelper (std::ostream &os, const signed char &v) |
| |
| static std::ostream & | ToStringHelper (std::ostream &os, const unsigned char &v) |
| |
Inheritance diagram for itk::simple::AntiAliasBinaryImageFilter: