Skip to content

Instantly share code, notes, and snippets.

@simbamangu
Created January 17, 2012 07:22
Show Gist options
  • Save simbamangu/1625422 to your computer and use it in GitHub Desktop.
Save simbamangu/1625422 to your computer and use it in GitHub Desktop.
Test of raster formats
for (i in 1:10) {
kid.img <- system.time(dem <- raster("~/workspace/TEMP/kideporegion.img"))
kid.asc <- system.time(dem <- raster("~/workspace/TEMP/kideporegion.asc"))
kid.tiffdeflate <- system.time(dem <- raster("~/workspace/TEMP/kideporegion_deflate.tif"))
kid.tiffnone <- system.time(dem <- raster("~/workspace/TEMP/kideporegion_none.tif"))
kid.tif <- system.time(dem <- raster("~/workspace/TEMP/kideporegion.tif"))
kid.tifpack <- system.time(dem <- raster("~/workspace/TEMP/kideporegion_pack.tif"))
kid.times <- data.frame(as.list(kid.img))
kid.times[2,] <- as.list(kid.asc)
kid.times[3,] <- as.list(kid.tiffdeflate)
kid.times[4,] <- as.list(kid.tiffnone)
kid.times[5,] <- as.list(kid.tifpack)
kid.times$type <- c("IMG", "ASC", "tif-Deflate", "tif-none", "tif-pack")
kid.times$diff <- (kid.times$elapsed / min(kid.times$elapsed)) - 1
kid.times$run <- i
if (i == 1) { summ <- kid.times }
else summ <- rbind(summ, kid.times)
}
final <- data.frame(tapply(summ$elapsed, summ$type, mean))
final <- cbind(final, tapply(summ$elapsed, summ$type, sd))
names(final) <- c("mean", "sd")
final$cv <- final$sd / final$mean
final$diff <- (final$mean - min(final$mean))
final
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment