Description Usage Arguments Details Value Author(s) References See Also Examples
Detection of edge points by the difference of two rotated and asymmetric Kernel- or M-Kernel-Estimators.
1 2 3 4 |
data |
numerical matrix representation of the (noisy) image. |
h1n, h2n |
positive numbers. Bandwidth for the kernels. |
asteps |
optional positive integer. Number of different angles used. |
estimator |
optional string. Estimator used within the windows. Possible values are:
|
kernel |
optional string. Kernel function for
|
score |
optional string. Score function for M-Kernel-Estimators
if
|
sigma |
optional positiv number. Parameter for the score function
|
kernelfunc |
optional function taking two numbers as arguments
and returning a positive number. Used as kernel function given
|
margin |
Optional value. Results near the margin are in
general not very reasonable. Setting |
edgepoints
implements several versions of the RDKE method,
introduced by Qiu in 1997.
The original method, which uses kernel estimates, is a generalized
version which uses M-Kernel-Estimators and two test procedures. The
test procedures are multiple tests for different angles for the
hypothesis of equal means (or medians) in both windows.
All methods apply rotating and scaling in the correct order (see
Garlipp, 2004).
A list of two numerical matrices. The first matrix contains the maximal jump height for every pixel if the chosen estimator is not a test procedure, and p-values otherwise. The second matrix contains the angle which leads to the maximal jump height or minimal p-value.
Tim Garlipp, TimGarlipp@gmx.de
Garlipp, T. (2004), On Robust Jump Detection in Regression Surface with Applications to Image Analysis, Carl-von-Ossietzky-Universitaet Oldenburg, Dissertation
Qiu, P. (1997), Nonparametric Estimation of Jump Surface, The Indian Journal of Statistics, 59A, No.2, 268-294.
1 2 3 4 5 6 7 8 9 10 | ## produce a matrix representation of a simple
## noisy image showing a black rectangle
y = matrix(rep(0, 60 * 60), nrow = 60)
y[21:40, 21:40] = 1
y = y + matrix(rnorm(60 * 60, 0, 0.2), nrow = 60)
image(y, col = gray(seq(0, 1, 1/255)))
## find the rectangle's edge points
ye = edgepoints(y, 0.05, 0.05, estimator = "M_median", kernel = "gauss")
image(ye[[1]] > 0.7, col = gray(c(1,0)))
|
row: 0/60
row: 1/60
row: 2/60
row: 3/60
row: 4/60
row: 5/60
row: 6/60
row: 7/60
row: 8/60
row: 9/60
row: 10/60
row: 11/60
row: 12/60
row: 13/60
row: 14/60
row: 15/60
row: 16/60
row: 17/60
row: 18/60
row: 19/60
row: 20/60
row: 21/60
row: 22/60
row: 23/60
row: 24/60
row: 25/60
row: 26/60
row: 27/60
row: 28/60
row: 29/60
row: 30/60
row: 31/60
row: 32/60
row: 33/60
row: 34/60
row: 35/60
row: 36/60
row: 37/60
row: 38/60
row: 39/60
row: 40/60
row: 41/60
row: 42/60
row: 43/60
row: 44/60
row: 45/60
row: 46/60
row: 47/60
row: 48/60
row: 49/60
row: 50/60
row: 51/60
row: 52/60
row: 53/60
row: 54/60
row: 55/60
row: 56/60
row: 57/60
row: 58/60
row: 59/60
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.