Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-19 Thread Mathieu Rajerison
gt; >> 2. Set the line to NA or 0. >> >> 3. Use clump() with directions=4. >> Jacob >> >> >> >> From: Barry Rowlingson >> To: Mathieu Rajerison >> Cc: r-sig-...@stat.math.ethz.ch >> Sent: Friday, 16 Se

Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-17 Thread Robert J. Hijmans
acob > > > ____ > From: Barry Rowlingson > To: Mathieu Rajerison > Cc: r-sig-...@stat.math.ethz.ch > Sent: Friday, 16 September 2011, 17:51 > Subject: Re: [R-sig-Geo] [raster] a railroad, a raster with a different > value on each side of it >

Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-16 Thread Jacob van Etten
September 2011, 17:51 Subject: Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it If your railroad is just a single line feature running approx N-S then: Create a box polygon for your study area, make sure the railroad just crosses it at N and S edges Use

Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-16 Thread Barry Rowlingson
On Fri, Sep 16, 2011 at 4:51 PM, Barry Rowlingson wrote: > If your railroad is just a single line feature running approx N-S then: > >  Create a box polygon for your study area, make sure the railroad just > crosses it at N and S edges > >  Use rgeos functions overlaying the RR line with the box p

Re: [R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-16 Thread Barry Rowlingson
If your railroad is just a single line feature running approx N-S then: Create a box polygon for your study area, make sure the railroad just crosses it at N and S edges Use rgeos functions overlaying the RR line with the box polygon to create the E and W polygons Job done? On Fri, Sep 16, 2

[R-sig-Geo] [raster] a railroad, a raster with a different value on each side of it

2011-09-16 Thread Mathieu Rajerison
Hi List, I've got a railroad shapefile and a corine land cover raster. I'd like to perform a land cover comparative analysis between the west side of the railroad and the right side. To do that, I tried first to have a raster with values=1 when on the east and 0 when on the west. But my calcula