[R-sig-Geo] [Fwd: Re: rgdal: problem with grid data]

[Agustin Lobo]

-------- Original Message --------
Subject: Re: [R-sig-Geo] rgdal: problem with grid data
Date: Tue, 23 Mar 2010 07:53:38 +0100
From: Agustin Lobo <aloboaleu at gmail.com>
Reply-To: Agustin.Lobo at ija.csic.es
To: Michael Sumner <mdsumner at gmail.com>
CC: sig-geo <r-sig-geo at stat.math.ethz.ch>
References: <4BA86074.8010709 at gmail.com> 
<522664f81003222349g39cc0c4ya756547d80c1d2f6 at mail.gmail.com>

ok, I see, I have to convert to Spatial Polygons data frame first, right?
Ideally, I should be able of clumping those cells with the same value,
so this is a bit more complicated than I expected. But should be possible.
I'll try and let the list know.
Thanks

Agus

Michael Sumner wrote:
> 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 <alobolistas at 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 at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
> 

-- 
Dr. Agustin Lobo
Institut de Ciencies de la Terra "Jaume Almera" (CSIC)
LLuis Sole Sabaris s/n
08028 Barcelona
Spain
Tel. 34 934095410
Fax. 34 934110012
email: Agustin.Lobo at ija.csic.es
http://www.ija.csic.es/gt/obster