[R-SIG-Mac] using RMySQL on OS X 10.6.6 R 2.12.2
Verbesselt, Jan
jan.verbesselt at wur.nl
Thu Mar 24 11:19:35 CET 2011
Hi,
I trying to export data from large images (satellite data) in to MySQL data via R. The code works fine but after some time it gives to following error:
"
*** caught bus error ***
address 0x100000012, cause 'non-existent physical address'
Traceback:
1: match.fun(FUN)
2: lapply(z, is.object)
3: unlist(lapply(z, is.object))
4: order(jseq)
5: `[<-.data.frame`(`*tmp*`, , value = list(gid = 78014, yr = 2009, dy = 49, ndvi = 6448, rel = 0))
6: `[<-`(`*tmp*`, , value = list(gid = 78014, yr = 2009, dy = 49, ndvi = 6448, rel = 0))
7: write.table(value[from:to, , drop = FALSE], file = conb, append = TRUE, quote = FALSE, sep = "\t", na = .MySQL.NA.string, row.names = FALSE, col.names = FALSE, eol = "\n", ...)
8: safe.write(value, file = fn)
9: mysqlWriteTable(conn, name, value, ...)
10: is(object, Cl)
11: is(object, Cl)
12: .valueClassTest(standardGeneric("dbWriteTable"), "logical", "dbWriteTable")
13: dbWriteTable(con, "datamod13q1", datamod13q1, append = T, row.names = 0)
aborting ...
Bus error
cube008:MODIS janv$
"
As an example I am adding the R script that I am using:
Which runs through each pixel and exports the data to add table:
# read all pixels within a certain boundary
gid <- 1
col <- 7:428
row <- 2557:2760
"
# Create first xy table from one image
for (i in row) {
for (j in col) {
rc <- c(i,j)-1 # line, sample !
# read data for specific row and column
ndvi <- getRasterTable(rndvi, offset=rc, region.dim=c(1, 1))
rel <- getRasterTable(rrel, offset=rc, region.dim=c(1, 1))
# prepare to insert into MySQL database
(datamod13q1 <- data.frame(gid=gid,yr=yr, dy=dy, ndvi=ndvi$band1/10000,rel=rel$band1))
dbWriteTable(con, "datamod13q1", datamod13q1, append = T, row.names = 0)
gid <- gid + 1
}
}
"
Thanks for your help
J
More information about the R-SIG-Mac
mailing list