threeStageParSel {DRIP}R Documentation

image denoising/deblurring, bandwidth selection, bootstrap


Select the bandwidth value for the image restoration method implemented in the function threeStage


threeStageParSel(image, bandwidth, edge1, edge2, nboot, blur=FALSE)



A square matrix object of size n by n, no missing value allowed.


Bandwidth values to be chosen from. Each of these values need to be an positive integer which specifies the number of pixels used in the local smoothing.


A matrix of 0 and 1 of the same size as image represents detected step edge pixels.


A matrix of 0 and 1 of the same size as image represents detected roof/valley edge pixels.


Required when blur is TRUE. Unused when blur is FALSE. An positive integer to specify the number of bootstraps to perform. See Qiu and Kang (2015) Statistica Sinica for details.


TRUE if the image contains blur, FALSE otherwise. If TRUE, the hybrid selection method proposed in Qiu and Kang (2015) Statistica Sinica is used. If FALSE, the leave-one-out cross validation is used.


Returns a list of the selected bandwdith, and a matrix of CV values with each entry corresponding to each choice of bandwdith.


Qiu, P., and Kang, Y. "Blind Image Deblurring Using Jump Regression Analysis," Statistica Sinica, 25, 2015, 879-899.


data(peppers) # Peppers image is bundled with the package and it is a
          # standard test image in image processing literature.
# Not run
#step.edges = stepEdgeLLK(peppers, 9, 17) # Step edge detection
#roof.edges = roofEdge(peppers, 6, 3000, edge1=step.edges) # Roof edge detection
#parSel = threeStageParSel(image=peppers, edge1=step.edges,
#edge2=roof.edges, bandwidth=4:5) # Time consuming

[Package DRIP version 1.7 Index]