Monica Pisica <pisicandru <at> hotmail.com> writes: > > > Hi, > > I have a matrix of data which i can vizualize as an image - for example. > I would like to save this image as a > geotiff file or at a tiff file with a world file which holds the projection > of my data (ultimately the data > represent a map of some sort). I know i can save the data as an ESRI > grid, but i am not interested in that. > > I wonder if anybody knows about any code which will help me do that.
See the writeGDAL() function in the rgdal package, and hints on how to turn your matrix into a SpatialGridDataFrame object in the sp package. For example: data(volcano) str(volcano) image(volcano) grd <- GridTopology(c(0.5, 0.5), c(1, 1), c(87, 61)) SGDF <- SpatialGridDataFrame(grd, data=data.frame(volcano=c(volcano))) image(SGDF, "volcano") See the proj4string= argument to insert the projection in valid PROJ.4 format. Then: writeGDAL(SGDF, "volcano.tif", drivername = "GTiff", ...) using the options= id needed, to pass through create options as on: http://www.gdal.org/frmt_gtiff.html For follow-ups, please consider the R-sig-geo list. Roger > > Thanks in advance, > > Monica > _________________________________________________________________ > > s. It's easy! > > aspx&mkt=en-us > [[alternative HTML version deleted]] > > ______________________________________________ > R-help <at> stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > > ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
