laplacelaplaceLaplacelaplaceLaplaceLaplace (Operator)

Name

laplacelaplaceLaplacelaplaceLaplaceLaplace — Calculate the Laplace operator by using finite differences.

Signature

Description

laplacelaplaceLaplacelaplaceLaplaceLaplace filters the input images ImageImageImageImageImageimage using a Laplace operator. Depending on the parameter FilterMaskFilterMaskFilterMaskFilterMaskFilterMaskfilterMask the following approximations of the Laplace operator are used:

'n_4'

                           1
                  1       -4      1
                           1

'n_8'

                  1       1       1
                  1       -8      1
                  1       1       1

'n_8_isotropic'

                  10      22     10
                  22     -128    22
                  10      22     10

For the three filter masks the following normalizations of the resulting gray values is applied, (i.e., by dividing the result by the given divisor): 'n_4' normalization by 1, 'n_8', normalization by 2 and for 'n_8_isotropic' normalization by 32.

For a Laplace operator with size 3x3, the corresponding filter is applied directly, while for larger filter sizes the input image is first smoothed using a Gaussian filter (see gauss_imagegauss_imageGaussImagegauss_imageGaussImageGaussImage) or a binomial filter (see binomial_filterbinomial_filterBinomialFilterbinomial_filterBinomialFilterBinomialFilter) of size MaskSizeMaskSizeMaskSizeMaskSizeMaskSizemaskSize-2. The Gaussian filter is selected for the above values of ResultTypeResultTypeResultTypeResultTypeResultTyperesultType. Here, MaskSizeMaskSizeMaskSizeMaskSizeMaskSizemaskSize = 5, 7, 9, 11, or 13 must be used. The binomial filter is selected by appending '_binomial'"_binomial""_binomial""_binomial""_binomial""_binomial" to the above values of ResultTypeResultTypeResultTypeResultTypeResultTyperesultType. Here, MaskSizeMaskSizeMaskSizeMaskSizeMaskSizemaskSize can be selected between 5 and 39. Furthermore, it is possible to select different amounts of smoothing for the column and row direction by passing two values in MaskSizeMaskSizeMaskSizeMaskSizeMaskSizemaskSize. Here, the first value of MaskSizeMaskSizeMaskSizeMaskSizeMaskSizemaskSize corresponds to the mask width (smoothing in the column direction), while the second value corresponds to the mask height (smoothing in the row direction) of the binomial filter. Therefore,

     laplace(O:R:'absolute',MaskSize,N:)

for MaskSize > 3 is equivalent to

     gauss_image(O:G:MaskSize-2:) >
     laplace(G:R:'absolute',3,N:).

and

     laplace(O:R:'absolute_binomial',MaskSize,N:)

is equivalent to

     binomial_filter(O:B:MaskSize-2,MaskSize-2:) >
     laplace(B:R:'absolute',3,N:).

laplacelaplaceLaplacelaplaceLaplaceLaplace either returns the absolute value of the Laplace filtered image (ResultTypeResultTypeResultTypeResultTypeResultTyperesultType = 'absolute'"absolute""absolute""absolute""absolute""absolute") in a byte or uint2 image or the signed result (ResultTypeResultTypeResultTypeResultTypeResultTyperesultType = 'signed'"signed""signed""signed""signed""signed" or 'signed_clipped'"signed_clipped""signed_clipped""signed_clipped""signed_clipped""signed_clipped"). Here, the output image type has the same number of bytes per pixel as the input image (i.e., int1 or int2) for 'signed_clipped'"signed_clipped""signed_clipped""signed_clipped""signed_clipped""signed_clipped", while the output image has the next larger number of pixels (i.e., int2 or int4) for 'signed'"signed""signed""signed""signed""signed".

Parallelization

• Multithreading type: reentrant (runs in parallel with non-exclusive operators).
• Multithreading scope: global (may be called from any thread).
• Automatically parallelized on tuple level.
• Automatically parallelized on channel level.
• Automatically parallelized on domain level.

Parameters

Example (C)

read_image(&Image,"mreut");
laplace(Image,&Laplace,"signed",3,"n_8_isotropic");
zero_crossing(Laplace,&ZeroCrossings);

Result

laplacelaplaceLaplacelaplaceLaplaceLaplace returns 2 (H_MSG_TRUE) if all parameters are correct. If the input is empty the behaviour can be set via set_system('no_object_result',<Result>)set_system("no_object_result",<Result>)SetSystem("no_object_result",<Result>)set_system("no_object_result",<Result>)SetSystem("no_object_result",<Result>)SetSystem("no_object_result",<Result>). If necessary, an exception is raised.

Possible Successors

zero_crossingzero_crossingZeroCrossingzero_crossingZeroCrossingZeroCrossing, dual_thresholddual_thresholdDualThresholddual_thresholdDualThresholdDualThreshold, thresholdthresholdThresholdthresholdThresholdThreshold

Alternatives

diff_of_gaussdiff_of_gaussDiffOfGaussdiff_of_gaussDiffOfGaussDiffOfGauss, laplace_of_gausslaplace_of_gaussLaplaceOfGausslaplace_of_gaussLaplaceOfGaussLaplaceOfGauss, derivate_gaussderivate_gaussDerivateGaussderivate_gaussDerivateGaussDerivateGauss

See also

highpass_imagehighpass_imageHighpassImagehighpass_imageHighpassImageHighpassImage, edges_imageedges_imageEdgesImageedges_imageEdgesImageEdgesImage

Module

Foundation