[R-sig-Geo] Exporting spatstat 'pixel image' to Arcview 9.x as an e00 or ascii
r.turner at auckland.ac.nz
Tue Jul 12 05:30:11 CEST 2011
This is definitely a case of the blind leading the blind, but I thought
Let "Im" be your pixel image. Do:
xy <- as.matrix(as.data.frame(raster.xy(as.owin(Im))))
v <- lookup.im(Im,xy[,1],xy[,2])
A <- SpatialPixelsDataFrame(points=xy,data=data.frame(v))
B <- as(A,"SpatialGridDataFrame")
Then "B" is an object of class "SpatialGridDataFrame" which ***may*** be
you need. I can't really tell, since I don't understand
The foregoing is kludgey, almost surely sub-optimal, and quite possibly
But there's a chance it will be of some use to you.
On 12/07/11 10:54, Honey-Marie de la Giroday wrote:
> Dear All,
> I have a pixel image (im) I received as a .RData file. I need to
> export this pixel image from R in a file format that can be imported
> into Arcview 9.x (e.g., .e00 or ascii).
> Below is some of the code that I have been trying to use. I have
> reviewed the manuals for rgdal, maptools, sp, and spatstat packages
> but have been unable to find another approach to converting the pixel
> image to an Arcview raster dataset. Please note that I have also
> searched the R-sig-geo archives and was unable to find this covered
> "#load spatst package
> #load file with model data
> #Plot best model data for 2004
> #Examine attributes of pixel image
>  "v" "dim" "xrange" "yrange" "xstep" "ystep" "xcol" "yrow"
>  "type" "units"
>  "im"
> #Load sp package
> #Coerce bestmodel2004 (pixel image) to SpatialGridDataFrame
> #......how do I do this?......
> #Output to Acrview = write.asciigrid(x, fname, attr=1, na.value=-9999.....
> Thank you for your time, and please let me know if I need to provide
> further information.
> H.M.C de la Giroday
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
More information about the R-sig-Geo