Re: [R-sig-Geo] Plotting direction vectors from an aspect map

2011-09-16 Thread bart
For doing this within R also have a look at ggplot2, with a simple sin and cos one could use geom_segment http://had.co.nz/ggplot2/geom_segment.html Bart On 09/16/2011 04:22 AM, Dylan Beaudette wrote: This comes to mind: http://casoilresource.lawr.ucdavis.edu/drupal/node/989 Note shameless

Re: [R-sig-Geo] Plotting direction vectors from an aspect map

2011-09-16 Thread Carson Farmer
I've done something similar for my own work (based on polygon centroids, but the idea is the same). The following (messy) code should do what you want using base graphics. The inputs are the coordinates of your cell centers, the angles of your arrows, and then a few additional parameters which

[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

Re: [R-sig-Geo] Converting spatialPointsDataFrame into ppp

2011-09-16 Thread GodinA
Dear all, I am new to spatial and point process analysis (as well as R), however most of my Ph.D. will rely heavily on these tools *sighs :) I am starting a new chapter where I intend to model elasmobranch i.e., shark and skate catch in the Northwest Atlantic Canadian waters using point process

Re: [R-sig-Geo] Converting spatialPointsDataFrame into ppp

2011-09-16 Thread Roger Bivand
On Fri, 16 Sep 2011, GodinA wrote: Dear all, I am new to spatial and point process analysis (as well as R), however most of my Ph.D. will rely heavily on these tools *sighs :) I am starting a new chapter where I intend to model elasmobranch i.e., shark and skate catch in the Northwest Atlantic

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,

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 b.rowling...@lancaster.ac.uk 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

[R-sig-Geo] projectRaster, MAJORITY method, rather than ngb?

2011-09-16 Thread Pieter Beck
Dear all, I am using raster's projectRaster to project a raster of a categorical variable to a grid of coarser resolution (which I specify using another raster object as the 'to' parameter). Rather than assigning the new cells the value originally found closest to their centers (i.e. the 'ngb'

[R-sig-Geo] Jenks classification for raster representation

2011-09-16 Thread Arnaud Mosnier
Dear list, In order to allow others benefiting from my errors, see below a presentation of a problem and it's solution by Roger Bivand (Thanks !). ## I want to use natural breaks (jenks) method to find class intervals into raster's values in order to plot it

[R-sig-Geo] GIS Lecturer position at the University of Illinois at Urbana-Champaign

2011-09-16 Thread Jonathan Greenberg
[Colleagues: Please forward this to any interested parties. Thanks!] The Department of Geography in the School of Earth, Society and Environment at the University of Illinois at Urbana-Champaign is looking for a Lecturer to fill a key teaching position in our Geographic Information Science

[R-sig-Geo] Sending raster:::aggregate a function with multiple parameters?

2011-09-16 Thread Jonathan Greenberg
R-sig-geo'ers: I wanted to run aggregate on a *degree* based azimuth image, with a circ.mean resampling. Since you can't just average degrees and get the right answer (e.g. the arithmetic mean of 359 degrees and 1 degrees is 180, but the true mean direction is 0), I adapted the circ.mean