RecursiveExponentialFilter3d

Smoothes a three-dimensional image with a recursive algorithm implementing an exponential filter.

Access to parameter description

This algorithm performs a 3D smoothing filter whose impulse response is an exponential function. The smoothing is achieved by a separable exponential decreasing filter with a variable scale factor.

• Response to horizontal impulse: $ f(x)=b(\alpha|x|+1)e^{-\alpha|x|} $
• Response to vertical impulse: $ ~~~~f(y)=b(\alpha|y|+1)e^{-\alpha|y|} $
• Response to depth impulse: $ ~~~~~~~f(z)=b(\alpha|z|+1)e^{-\alpha|z|} $

Where:

• $x$, $y$ and $z$ represent the offsets from the pixel to process.
• $ \alpha $ is a sharpness factor. Low values produce high level of smoothing.
• $ b=\frac{\alpha}{4} $

The lower the $\alpha$ value, the smoother the result. Practically, this filter response is mostly visible in the [0.3, 5.3] range.
For convenience, the scale factor is controlled by a spread value in order to convert it between 0 et 100. A value of 0 applies a very limited smoothing, while 100 induces a strong smoothing: $$ SpreadValue=\frac{(5.3-\alpha)}{5}\times 100 $$ This filter is very close to a Gaussian filter. The main difference is the alpha coefficient, which has an inverse behavior compared to a Gaussian standard deviation.

This implementation is based on an Infinite Impulse Response (IIR) algorithm. It computes the sum of one causal and one anti-causal filter where previous results are used to compute the next result. Using this mode, the computation time is independent of the spread parameter.

See also

• GaussianFilter3d
• GradientOperator3d
• ExponentialFilter2d
• RecursiveExponentialFilter2d

// Function prototype.
std::shared_ptr< iolink::ImageView >
recursiveExponentialFilter3d( std::shared_ptr< iolink::ImageView > inputImage,
                              int32_t spreadValue,
                              std::shared_ptr< iolink::ImageView > outputImage = NULL );

// Function prototype.
recursive_exponential_filter_3d( input_image, spread_value = 60, output_image = None )

// Function prototype.
public static IOLink.ImageView
RecursiveExponentialFilter3d( IOLink.ImageView inputImage,
                              Int32 spreadValue = 60,
                              IOLink.ImageView outputImage = null );

// Command constructor.
RecursiveExponentialFilter3d();

// Get method of the inputImage parameter.
std::shared_ptr< iolink::ImageView > inputImage() const;
// Set method of the inputImage parameter.
void setInputImage( std::shared_ptr< iolink::ImageView > inputImage );

// Get method of the spreadValue parameter.
int32_t spreadValue() const;
// Set method of the spreadValue parameter.
void setSpreadValue( const int32_t& spreadValue );

// Get method of the outputImage parameter.
std::shared_ptr< iolink::ImageView > outputImage() const;
// Set method of the outputImage parameter.
void setOutputImage( std::shared_ptr< iolink::ImageView > outputImage );

// Method to launch the command.
void execute();

// Property of the inputImage parameter.
RecursiveExponentialFilter3d.input_image
// Property of the spreadValue parameter.
RecursiveExponentialFilter3d.spread_value
// Property of the outputImage parameter.
RecursiveExponentialFilter3d.output_image

// Method to launch the command.
execute()

// Command constructor.
RecursiveExponentialFilter3d()

// Property of the inputImage parameter.
RecursiveExponentialFilter3d.inputImage
// Property of the spreadValue parameter.
RecursiveExponentialFilter3d.spreadValue
// Property of the outputImage parameter.
RecursiveExponentialFilter3d.outputImage

// Method to launch the command.
Execute()

auto foam = readVipImage( std::string( IMAGEDEVDATA_IMAGES_FOLDER ) + "foam.vip" );

RecursiveExponentialFilter3d recursiveExponentialFilter3dAlgo;
recursiveExponentialFilter3dAlgo.setInputImage( foam );
recursiveExponentialFilter3dAlgo.setSpreadValue( 60 );
recursiveExponentialFilter3dAlgo.execute();

std::cout << "outputImage:" << recursiveExponentialFilter3dAlgo.outputImage()->toString();

foam = imagedev.read_vip_image(imagedev_data.get_image_path("foam.vip"))

recursive_exponential_filter_3d_algo = imagedev.RecursiveExponentialFilter3d()
recursive_exponential_filter_3d_algo.input_image = foam
recursive_exponential_filter_3d_algo.spread_value = 60
recursive_exponential_filter_3d_algo.execute()

print( "output_image:", str( recursive_exponential_filter_3d_algo.output_image ) );

ImageView foam = Data.ReadVipImage( @"Data/images/foam.vip" );

RecursiveExponentialFilter3d recursiveExponentialFilter3dAlgo = new RecursiveExponentialFilter3d
{
    inputImage = foam,
    spreadValue = 60
};
recursiveExponentialFilter3dAlgo.Execute();

Console.WriteLine( "outputImage:" + recursiveExponentialFilter3dAlgo.outputImage.ToString() );

auto foam = readVipImage( std::string( IMAGEDEVDATA_IMAGES_FOLDER ) + "foam.vip" );

auto result = recursiveExponentialFilter3d( foam, 60 );

std::cout << "outputImage:" << result->toString();

foam = imagedev.read_vip_image(imagedev_data.get_image_path("foam.vip"))

result = imagedev.recursive_exponential_filter_3d( foam, 60 )

print( "output_image:", str( result ) );

ImageView foam = Data.ReadVipImage( @"Data/images/foam.vip" );

IOLink.ImageView result = Processing.RecursiveExponentialFilter3d( foam, 60 );

Console.WriteLine( "outputImage:" + result.ToString() );