Nothing
CARsimu <- function(LEVEL=6, rho=0.2499, row2=0, col2=0, rc1=0, cr1=0, maindi=1, rajz=TRUE) {
#--------------------------------------------------------------
#
# TITLE: CARsimu()
# AUTHOR: FERKO CSILLAG AND SANDOR KABOS, MODIFIED: TARMO REMMEL
# DATE: 23 January 2020
# CALLS: NA
# CALLED BY: singlemap()
# NEEDS: NA
# NOTES: SIMULATES CONDITIONAL AUTOREGRESSIVE (CAR) IMAGES
# THAT CAN BE CUT LATER BY PROPORTION TO PRODUCE
# BINARY MAPS. THE RHO PARAMETERS MUST SUM TO <0.5
# THE RESULTS ARE STATIONARY AND ISOTROPIC.
#
# WAS PREVIOUSLY wujabki()
#--------------------------------------------------------------
# SAVE GRAPHIC PARAMETERS AND RESTATE THEM ON EXIT
opar <- par(no.readonly=TRUE)
on.exit(par(opar))
col1 <- rho
row1 <- rho
M <- 2^LEVEL
NN <- M^2
FNN <- NN/2
I1 <- array(rep(0, NN))
I1[1] <- 1
NR1 <- array(rep(0, NN))
NR1[2] <- 1
NR1[NN] <- 1
NC1 <- array(rep(0, NN))
NC1[M + 1] <- 1
NC1[NN - M + 1] <- 1
NR2 <- array(rep(0, NN))
NR2[3] <- 1
NR2[NN - 1] <- 1
NC2 <- array(rep(0, NN))
NC2[2 * M + 1] <- 1
NC2[NN - 2 * M + 1] <- 1
NRC1 <- array(rep(0, NN))
NRC1[M + 2] <- 1
NRC1[NN - M] <- 1
NCR1 <- array(rep(0, NN))
NCR1[M] <- 1
NCR1[NN - M + 2] <- 1
k <- maindi * I1 - row1 * NR1 - col1 * NC1 - row2 * NR2 - col2 * NC2 - rc1 * NRC1 - cr1 * NCR1
d2 <- fft(k, inverse = FALSE)
qse <- sqrt(Re(d2))
par(mfrow = c(1, 1))
nkar <- 0
while(nkar == 0) {
xe <- rnorm(NN)
x <- fft(xe)
for(i in 1:NN) {
x[i] <- x[i]/qse[i]
}
y <- (fft(x, inverse = TRUE))/M
yre <- Re(y)
A <- matrix(yre, byrow = TRUE, nrow = M, ncol = M)
nkar <- 3
}
A <- A/sum(A*A)
if(rajz) {
# INSERT POSSIBLE DRAWING FUNCTIONALITY HERE BUT BEST TO MOVE IT OUT OF THIS FUNCTION
}
return(A)
} # END FUNCTION: CARsimu
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.