# R and the Geostatistical Software Library data format

[This article was first published on

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

**Bart Rogiers - Sreigor**, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

The *.gslib file format originates from the Geostatistical Software Library, but is also used in the follow-up of that software, i.e. the Stanford Geostatistical Modelling Software (SGeMS). Since not all geostatistical algorithms are implemented in R yet (especially the multiple-point ones), it is useful to have these functions for interacting with SGeMS. I also used it for 3D visualisation of geological datasets that I produced with R.

I extended some code that can be found here for writing *.gslib files from an R dataframe, with an added variable to enable a numeric no data flag. A reading function is also included to be able to reload the files back into R. Not so hard to program, but it might be useful to somebody..

write.gslib <- function(dat, file, nodatavalue=-99999) { cat('GSLIB file created in R\n', file=file) cat(length(names(dat)), file=file, append=TRUE) cat('\n', file=file, append=TRUE) write(cbind(names(dat)), file=file, append=TRUE) for(i in 1:ncol(dat)) dat[is.na(dat[,i]),i] <- nodatavalue for(i in 1:ncol(dat)) dat[dat[,i]=='NA',i] <- nodatavalue for(i in 1:ncol(dat)) dat[is.infinite(dat[,i]),i] <- nodatavalue write.table(dat,file=file, append=TRUE, sep='\t', col.names=FALSE, row.names=FALSE) } read.gslib <- function(filename) { char <- scan(filename, what=character()) nums <- which(as.numeric(char)!='NaN') nrcol <- as.numeric(char[nums[1]]) begin <- nums[2] eind <- nums[length(nums)] nrrows <- (eind-begin+1)/nrcol indices <- seq(1,nrrows*nrcol-1,nrcol) mat <- matrix(ncol=nrcol, nrow=nrrows) i <- 1 while(i <= nrcol) { mat[,i] <- as.numeric(char[(indices+i+begin-2)]) i <- i+1 } mat <- as.data.frame(mat) names(mat) <- c(char[(nums[1]+c(1:nrcol))]) return(mat) }

To

**leave a comment**for the author, please follow the link and comment on their blog:**Bart Rogiers - Sreigor**.R-bloggers.com offers

**daily e-mail updates**about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.

Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.