What are you expecting to get - a point for every raster cell? If so,
use gridded(isrlc) <- FALSE to convert to a SPDF first - then the
write will work.
More complicated conversions are possible, such as with
as.SpatialPolygons.GridTopology.
Cheers, Mike.
On Tue, Mar 23, 2010 at 5:32 PM, Agustin Lobo <alobolis...@gmail.com> wrote:
> Hi!
>
> I read a grid file into R with
> isrlc <-
> readGDAL("/media/Transcend/MASTER_ICTA2007_2008/GEODATA_miniprojs/Kefah/ISR_cov/isr_cov.gri")
>
> While isrlc seems correct, I'm getting an empty shapefile when I attempt to
> write:
>
> writeOGR(isrlc,dsn="/media/Transcend/MASTER_ICTA2007_2008/GEODATA_miniprojs/Kefah",
> layer="ISR_cov2",driver="ESRI Shapefile")
>
> Data from
> http://www.diva-gis.org/data/cov/ISR_cov.zip
>
> Thanks
>
> Agus
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
