https://github.com/cran/fields
Tip revision: beb6c9118b4f355fc630943f00274a7df8714fe1 authored by Doug Nychka on 26 April 2007, 00:00:00 UTC
version 3.5
version 3.5
Tip revision: beb6c91
quilt.plot.R
"quilt.plot" <-
function(x,y,z,nrow=64, ncol=64, grid=NULL, add.legend=TRUE,add=FALSE,...){
x<- as.matrix(x)
if( ncol(x)==2){
z<- y}
if( ncol(x)==1){
x<- cbind( x,y)}
if( ncol(x)==3){
z<- x[,3]
x<- x[,1:2]
}
# at this point x should be a 2 column matrix of x-y locations
# z is a vector or one column matrix of the z values.
#discretize data
as.image( z, x=x, nrow=nrow, ncol=ncol, na.rm=TRUE)-> out.p
#plot it
if( add.legend){
image.plot( out.p,col=tim.colors(64),add=add,...)
}
else{
image(out.p, col=tim.colors(64),add=add,...)
}
}