I am attempting to use the Bayesian kriging function in geoR by Ribeiro
and I am encountering a problem with the "image" command.
The program is not responsive to the "ylim" command, such that:
image(temp.bayes,loc=pred.grid,col=rainbow(25, start=.6,
end=.1),xlab="Longitude",ylab="Latitude",
ylim=c(-30,30),borders=NULL)
legend.krige(y.leg=c(-150,-140), x.leg=c(-150,150),
values=temp.bayes$predictive$mean.simulations,col=rainbow(25, start=.6,
end=.1))
map('worldHires', ylim=c(-30,30), add=TRUE,interior=FALSE)
title(expression(paste("a.)",Delta, "T", degree,C)))
fails to generate a figure with the y-axis from 30 degrees south to 30
degrees north. Interestingly, the image command is responsive to the
"xlim"
command
image(temp.bayes,loc=pred.grid,col=rainbow(25, start=.6,
end=.1),xlab="Longitude",ylab="Latitude",
xlim=c(-30,30),borders=NULL)
but this changes the scale of both the y and x axes. I have tried
exporting the gridded data to the base R package to use the standard
"image" command; however the data are in a very cumbersome format.
Please help if you have encountered such a problem.
Thank you,
Ashley