Description Usage Arguments Details Value See Also Examples
This function decomposes an image into frequency component and residue of twodimensional ensemble patch transform by sifting process.
1 2 3 4 
x, y 
locations of regular grid at which the values in image 
z 
matrix of an image observed at location 
type 
patch type of 
tau 
a size parameter for twodimensional ensemble patch transform: 
theta 
a degree of clockwise rotation of a patch. 
process 
specifies transform types for patch and ensemble processes:

pquantile 
quantiles for lower and upper envelopes of patch transform. When it is 
equantile 
quantiles for lower and upper envelopes of ensemble patch transform. 
gamma 
controls the amount of envelope magnitude. 
boundary 
specifies boundary condition from 
stoprule 
stopping rule of sifting. The 
tol 
tolerance for stopping rule of sifting. 
maxiter 
the maximum number of sifting. 
check 
specifies whether the sifting process is displayed. When 
This function decomposes an image into frequency component and residue of twodimensional ensemble patch transform by sifting process for a size parameter.
eptcomp 
list of ept (ensemble patch transform) component at each sifting step when 
FC 
frequency component of ensemble patch transform by sifting process. 
residue 
residue of ensemble patch transform by sifting process. 
parameters 
a list of input parameters of 
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24  #### example : composite of two components having different frequencies
nr < nc < 128; x < seq(0, 1, length=nr); y < seq(0, 1, length=nc)
coscomp1 < outer(cos(20 * pi * x), cos(20 * pi * y))
coscomp2 < outer(cos(5* pi * x), cos(5 * pi * y))
cosmeanf < coscomp1 + coscomp2
op < par(mfcol=c(3,1), mar=c(0,0.5,2,0.5))
image(cosmeanf, xlab="", ylab="", col=gray(0:100/100), axes=FALSE, main="a composite image")
image(coscomp1, xlab="", ylab="", col=gray(0:100/100), axes=FALSE, main="highfrequency component")
image(coscomp2, xlab="", ylab="", col=gray(0:100/100), axes=FALSE, main="lowfrequency component")
#### Decomposition by Ensemble Patch Transform
outcossift < eptdecomp2d(z=cosmeanf, tau=8)
#### Decomposition Result
par(mfrow=c(2,2), mar=c(2,2,2,1))
image(outcossift$FC, xlab="", ylab="", col=gray(0:100/100), axes=FALSE, main="Decomposed HF")
persp(outcossift$FC, theta = 30, phi = 45, col = "white", xlab="X", ylab="Y", main="Decomposed HF")
image(outcossift$residue, xlab="", ylab="", col=gray(0:100/100), axes=FALSE, main="Residue")
persp(outcossift$residue, theta = 30, phi = 45, col = "white", xlab="X", ylab="Y", main="Residue")
par(op)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.