Re: [R-sig-Geo] Availability of MacOS X binary for rgeos

2011-09-22 Thread Roger Bivand
On Fri, 23 Sep 2011, hi_ono2...@ybb.ne.jp wrote: Hello. Does anyone know this? As the package page on CRAN http://cran.r-project.org/web/packages/rgdal/index.html says: rgdal binaries for OSX are available from the CRAN (extras) repository which you access by setting: setRepositories(in

Re: [R-sig-Geo] Availability of MacOS X binary for rgeos

2011-09-22 Thread Ben Madin
If you install the geos available from http://www.kyngchaos.com/software/frameworks the rgeos from source code works OK as long as you direct the lib and include to the right places. You may want proj4 as well. Check the dependencies, I think you will need the UNIX Image IO frameworks as well

Re: [R-sig-Geo] Distance from coast

2011-09-22 Thread Michael Sumner
Here's a really simple example using the default low res longlat coastline in the maps package. It should work for any coastal dataset that can extract the raw coordinates from, though there's no attempts at efficiency. spDistsN1 is really fast and works well for this though, be careful to use long

[R-sig-Geo] Distance from coast

2011-09-22 Thread Sandeep Patil
Hello does anyone know of any R -script or function which will give the value of nearest distance from coast of a point in USA. I have a number of latitude and longitude values spread evenly across USA and i want to find the nearest distance from coast for each of these points. If not a function

[R-sig-Geo] Availability of MacOS X binary for rgeos

2011-09-22 Thread hi_ono2001
Hello. Does anyone know this? Thanks in advance. ___ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo

[R-sig-Geo] Spatial Error Model for panel data- large spatial weight matrix problem

2011-09-22 Thread Rafael Terra
Dear R users, I'm estimating  a spatial error model for panel data by maximum likelihood (library splm, function spml). I'm using R version 2.13.1 (2011-07-08) platform: x86_64-pc-mingw32/x64 (64-bit). The operational system is Windows 7 64-bit and my computer has 4Gb of RAM. My panel data has

Re: [R-sig-Geo] small data set for an example

2011-09-22 Thread juantomas . sayago
Sorry for cross posting You can check also the asdarbook website or the cran of different packages they have really nice data. Juan Sent on the Sprint® Now Network from my BlackBerry® -Original Message- From: Greg Snow Sender: r-sig-geo-bounces@r-project.orgDate: Thu, 22 Sep 2011 08:56:

Re: [R-sig-Geo] small data set for an example

2011-09-22 Thread Greg Snow
In the TeachingDemos package there are 3 datasets 'ccc', 'h2h', and 'towork' which are real, not overly large (but you can decide if they are small or not) which contain lat, long, and other variables. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@

Re: [R-sig-Geo] symbols in spplot

2011-09-22 Thread piero campa
Dear all, I would like to share a function that, as suggested by Roger, displays arrows as SpatialLinesDataFrame objects that can be plotted with the spplot command. Little reproducible example: > # Setup > data(meuse) > data(meuse.grid) > attach(meuse) > coordinates(meuse) <- ~x+y > coordinates(m

Re: [R-sig-Geo] about ST kriging

2011-09-22 Thread Pınar Aslantas Bostan
Edzer, I have 256 meteorological station and 39 year value. > dim(w) [1] 256 39 1 Before I tried only spatial kriging and I did not get memory.size error. Just a moment ago, I used another computer (64-bit Windows and 6 GB RAM), again I got the same error. -Original Message- From:

Re: [R-sig-Geo] [spatstat] problem with a psp Line Segment Pattern when plotting and pointsOnlines performing

2011-09-22 Thread Mathieu Rajerison
Ok, thanks a lot. I wonder how I could have digitized one segment of length 0 (???) Anyway, I thank you a lot for your answer. I think it will solve my problem for the moment. Best, Mathieu 2011/9/22 Michael Sumner > It seems it is because one of your line segments is of length zero, > usin

Re: [R-sig-Geo] [spatstat] problem with a psp Line Segment Pattern when plotting and pointsOnlines performing

2011-09-22 Thread Michael Sumner
It seems it is because one of your line segments is of length zero, using your attachment: load("track.RData") library(spatstat) which(!lengths.psp(track.psp) > 0) [1] 17 Dropping the offender: noZeroLengths <- track.psp[lengths.psp(track.psp) > 0] (pts <- pointsOnLines(noZeroLengths, eps=50))

Re: [R-sig-Geo] downscaling grid data; alternative for interp()?

2011-09-22 Thread Paul Hiemstra
Hi, A trick to get Nearest Neighbor (NN) is to use idw() from the gstat package with setting the nmax parameter to 1. This limits the amount of points used for the inverse distance weighted average to one, leading to NN. I'm not sure if this is much faster than interp, but you can take a look. r

Re: [R-sig-Geo] about ST kriging

2011-09-22 Thread Edzer Pebesma
Pinar, what are the dimensions (space and time) of your data, i.e. what does dim(w) return? Have you tried kriging with a (spatial and/or temporal) subset? On 09/22/2011 01:38 PM, Pınar Aslantas Bostan wrote: > Hi Edzer, > > Thank you very much for your mail. I found two records have identical

[R-sig-Geo] downscaling grid data; alternative for interp()?

2011-09-22 Thread Klaus Vormoor
Hi, I am struggling with this quite a while already, and I haven't found a good solution so far. I want to disaggregate precipitation data. For this I have a grid with daily precipitation in 1x1 km resolution, and additionally 24 grids with 10x10 km resolution giving percentage charge of hourly pr

Re: [R-sig-Geo] about ST kriging

2011-09-22 Thread Pınar Aslantas Bostan
Hi Edzer, Thank you very much for your mail. I found two records have identical location, but it is not my fault because I checked the original source data and it is erroneous. Now I have another problem; R gives another error message: Error: cannot allocate vector of size 95.4 Mb My computer

Re: [R-sig-Geo] [spatstat] problem with a psp Line Segment Pattern when plotting and pointsOnlines performing

2011-09-22 Thread Mathieu Rajerison
Thanks for the answer! I coerced my SpatialLinesDataFrame object to a psp object using a non-NA column. Plotting is OK but pointsOnLines still fails: > track.psp <- as(track[, 2], "psp") > plot(track.psp) > pts<-pointsOnLines(track.psp, eps=50) Erreur dans ppp(x = x, y = y, window = win) : all c

Re: [R-sig-Geo] Calculate the surface of polygons within gridcells with gIntersection

2011-09-22 Thread Eelke Folmer
Hi Holly, Could you provide some insight into your int object. Does it happen for all i? My first guess is that your matrix 'int' might not be what you expect; try: grid[int[,i],] and see whether it is a geometry. Eelke Op 8:59 PM, Whiteley,Holly schreef: Dear list, I was wondering if someon

Re: [R-sig-Geo] small data set for an example

2011-09-22 Thread Barry Rowlingson
On Thu, Sep 22, 2011 at 1:14 AM, Hodgess, Erin wrote: > Dear R Sig Geo People: > > Does anyone have a small "real" data set with the longitude and latitude, > along with a couple of data values per location that they would be willing to > share, please? > > > > I'm using simulated data but thoug

Re: [R-sig-Geo] [spatstat] problem with a psp Line Segment Pattern when plotting and pointsOnlines performing

2011-09-22 Thread Michael Sumner
The problem is that the first attribute on the track object, "ID" is NA - and when plot.psp is trying to draw a legend it has nothing to base the image "ribbon" on. This should avoid the error: plot(track.psp, ribbon = FALSE) Try subsetting the attributes in the coercion to psp, choose which col

Re: [R-sig-Geo] Calculate the surface of polygons within gridcells with gIntersection

2011-09-22 Thread Whiteley,Holly
Dear list, I was wondering if someone could help me with polygon-polygon spatial queries? Similar to the below request in a previous thread, I have a SpatialPolygonDataFrame for an environmental factor (seabed sediment), and a gridded SpatialPolygonDataFrame. I am attempting to calculate the

[R-sig-Geo] [spatstat] problem with a psp Line Segment Pattern when plotting and pointsOnlines performing

2011-09-22 Thread Mathieu Rajerison
Hi, I Have a problem with a psp object. I can't plot it, neither perform a pointsOnLines calculation on it. --- Here is what I obtain with my different operations: > track <- readOGR("voie.shp", "voie") OGR data source with driver: ESRI Shapefile Source: "voie.shp", layer: "voie" with 1 features

[R-sig-Geo] FW: GWAR - R

2011-09-22 Thread Inaki Sagarzazu
Hi, I am trying to use gwr in R but I'm having problems because some of my polygons have missing values (NA) for the variables of interest. Do you have any suggestions on how to solve this? thanks

Re: [R-sig-Geo] help to change charactwer name in a row of data.frame

2011-09-22 Thread Paul Hiemstra
On 09/21/2011 03:05 PM, Ndjido Ardo BAR wrote: > Dear Gianni, > > you can also convert your x-column to be a character and then assign any > value to the line you will, just as follows > > df <- data.frame(x=c("a","b","c"),y=c(1,2,3)) > > df$x<-as.character(df$x) > > df[which(df$x=="b"),'x']<-"b%"