derivate_gaussderivate_gaussDerivateGaussDerivateGaussderivate_gauss (算子)

名称

derivate_gaussderivate_gaussDerivateGaussDerivateGaussderivate_gauss — 用高斯导数卷积图像。

签名

描述

derivate_gaussderivate_gaussDerivateGaussDerivateGaussDerivateGaussderivate_gauss convolves an image with the derivatives of a Gaussian and calculates various features derived therefrom. SigmaSigmaSigmaSigmasigmasigma is the parameter of the Gaussian (i.e., the amount of smoothing). If one value is passed in SigmaSigmaSigmaSigmasigmasigma the amount of smoothing in the column and row direction is identical. If two values are passed in SigmaSigmaSigmaSigmasigmasigma the first value specifies the amount of smoothing in the column direction, while the second value specifies the amount of smoothing in the row direction. The possible values for ComponentComponentComponentComponentcomponentcomponent are:

'none'"none""none""none""none""none"：

Smoothing only.

'x'"x""x""x""x""x"：

First derivative along x.

'y'"y""y""y""y""y"：

First derivative along y.

'gradient'"gradient""gradient""gradient""gradient""gradient"：

Absolute value of the gradient.

'gradient_dir'"gradient_dir""gradient_dir""gradient_dir""gradient_dir""gradient_dir"：

Gradient direction in radians.

'xx'"xx""xx""xx""xx""xx"：

Second derivative along x.

'yy'"yy""yy""yy""yy""yy"：

Second derivative along y.

'xy'"xy""xy""xy""xy""xy"：

Second derivative along x and y.

'xxx'"xxx""xxx""xxx""xxx""xxx"：

Third derivative along x.

'yyy'"yyy""yyy""yyy""yyy""yyy"：

Third derivative along y.

'xxy'"xxy""xxy""xxy""xxy""xxy"：

Third derivative along x, x and y.

'xyy'"xyy""xyy""xyy""xyy""xyy"：

Third derivative along x, y and y.

'det'"det""det""det""det""det"：

Determinant of the Hessian matrix:

'laplace'"laplace""laplace""laplace""laplace""laplace"：

Laplace operator (trace of the Hessian matrix):

'mean_curvature'"mean_curvature""mean_curvature""mean_curvature""mean_curvature""mean_curvature"：

Mean curvature H

'gauss_curvature'"gauss_curvature""gauss_curvature""gauss_curvature""gauss_curvature""gauss_curvature"：

Gaussian curvature K

'area'"area""area""area""area""area"：

Differential Area A

'eigenvalue1'"eigenvalue1""eigenvalue1""eigenvalue1""eigenvalue1""eigenvalue1"：

First eigenvalue

'eigenvalue2'"eigenvalue2""eigenvalue2""eigenvalue2""eigenvalue2""eigenvalue2"：

Second eigenvalue

'eigenvec_dir'"eigenvec_dir""eigenvec_dir""eigenvec_dir""eigenvec_dir""eigenvec_dir"：

Direction of the eigenvector corresponding to the first eigenvalue in radians

'kitchen_rosenfeld'"kitchen_rosenfeld""kitchen_rosenfeld""kitchen_rosenfeld""kitchen_rosenfeld""kitchen_rosenfeld"：

Second derivative perpendicular to the gradient

'zuniga_haralick'"zuniga_haralick""zuniga_haralick""zuniga_haralick""zuniga_haralick""zuniga_haralick"：

Normalized second derivative perpendicular to the gradient

'2nd_ddg'"2nd_ddg""2nd_ddg""2nd_ddg""2nd_ddg""2nd_ddg"：

Second derivative along the gradient

'de_saint_venant'"de_saint_venant""de_saint_venant""de_saint_venant""de_saint_venant""de_saint_venant"：

Second derivative along and perpendicular to the gradient

注意

Besides the pure C version there are specific implementations of derivate_gaussderivate_gaussDerivateGaussDerivateGaussDerivateGaussderivate_gauss for speed up. Such an optimization is applied in case it is supported by the system and the respective system parameter *_enable is set to 'true'"true""true""true""true""true", see set_systemset_systemSetSystemSetSystemSetSystemset_system。The following optimizations are supported (listed according to their priority):

• using AVX512f instructions ('avx512f_enable'"avx512f_enable""avx512f_enable""avx512f_enable""avx512f_enable""avx512f_enable")

• using AVX instructions ('avx_enable'"avx_enable""avx_enable""avx_enable""avx_enable""avx_enable")

• using SSE2 instructions ('sse2_enable'"sse2_enable""sse2_enable""sse2_enable""sse2_enable""sse2_enable")

These implementations are slightly inaccurate compared to the pure C version due to numerical issues. For example, using SSE2 instructions the inaccuracy is in order of magnitude of 1.0e-5 for 'byte' images and ComponentComponentComponentComponentcomponentcomponent set to 'none'"none""none""none""none""none", 'x'"x""x""x""x""x", or 'y'"y""y""y""y""y".

In case accuracy is preferred over performance, set all corresponding system parameter to 'false'"false""false""false""false""false" (using set_systemset_systemSetSystemSetSystemSetSystemset_system) before calling derivate_gaussderivate_gaussDerivateGaussDerivateGaussDerivateGaussderivate_gauss. This way derivate_gaussderivate_gaussDerivateGaussDerivateGaussDerivateGaussderivate_gauss does not use the accelerations. Do not forget to set the parameter back to 'true'"true""true""true""true""true" afterwards.

derivate_gaussderivate_gaussDerivateGaussDerivateGaussDerivateGaussderivate_gauss is only executed on an OpenCL device if SigmaSigmaSigmaSigmasigmasigma induces a filter width respectively height of up to 129 pixels. This corresponds to a SigmaSigmaSigmaSigmasigmasigma of less than 20.7 for ComponentComponentComponentComponentcomponentcomponent = 'none'"none""none""none""none""none". The OpenCL implementation is slightly inaccurate compared to the pure C version due to numerical issues.

请注意，若使用域缩减后的图像作为输入，滤波器算子可能会返回意外结果。请参阅 滤波器 一章

执行信息

• 支持 OpenCL 计算设备。
• 多线程类型：可重入（与非独占算子并行运行）。
• 多线程作用域：全局（可从任何线程调用）。
• 在元组级别上自动并行化。
• 在通道级别上自动并行化。
• 在域级别上自动并行化。

参数

示例（C）

read_image(&Image,"mreut");
derivate_gauss(Image,&Gauss,3.0,"x");
zero_crossing(Gauss,&ZeroCrossings);

可能的后继

zero_crossingzero_crossingZeroCrossingZeroCrossingZeroCrossingzero_crossing, dual_thresholddual_thresholdDualThresholdDualThresholdDualThresholddual_threshold

替代

laplacelaplaceLaplaceLaplaceLaplacelaplace, laplace_of_gausslaplace_of_gaussLaplaceOfGaussLaplaceOfGaussLaplaceOfGausslaplace_of_gauss, binomial_filterbinomial_filterBinomialFilterBinomialFilterBinomialFilterbinomial_filter, gauss_filtergauss_filterGaussFilterGaussFilterGaussFiltergauss_filter, smooth_imagesmooth_imageSmoothImageSmoothImageSmoothImagesmooth_image, isotropic_diffusionisotropic_diffusionIsotropicDiffusionIsotropicDiffusionIsotropicDiffusionisotropic_diffusion

另见

zero_crossingzero_crossingZeroCrossingZeroCrossingZeroCrossingzero_crossing, dual_thresholddual_thresholdDualThresholdDualThresholdDualThresholddual_threshold

模块

基础