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.
>
>