[R-sig-Geo] sm.sphere

2008-08-18 Thread Galkowski, Jan
In the sm package, the Description for sm.sphere promises "The angle of view may be altered.", and the Value writeup says "... a list containing the value of the smoothing parameter and the rotation angles of the displayed plot." Yet, if I invoke this in combinations of TRUE, FALSE settings for "h

[R-sig-Geo] estimates for paths of travel given two time series, drawn from two corresponding spatial densities

2008-02-28 Thread Galkowski, Jan
This is less an R question and more a request to the community for literature references and the like. I'm interested in the following inferential question. What's the max likelihood estimate for a path of travel of a point across a plane having stable, smooth densities of two separate variables

Re: [R-sig-Geo] snippet of code for Voronoi edge presentation

2008-01-03 Thread Galkowski, Jan
Just to followup for others, don't need help any longer. I've mastered PBSmapping's PolySet, I believe, that that's solved a number of these problems, including the business about having to write out and back in. Same usage manipulates world maps to create appropriate duplicates. PolySet is docume

[R-sig-Geo] snippet of code for Voronoi edge presentation

2008-01-02 Thread Galkowski, Jan
I'm implementing STRIPACK, the Delaunay triangulator/Voronoi tessellator for spheres, working from FORTRAN 77 code available via the ACM Digital Library. I'm also interested in displaying my data, and I have "operationally toured" a bunch of what's available for R, including rgdal, sp, TeachingDem

Re: [R-sig-Geo] adding new maps to map, mapdata, mapproj; archive of R-sig-geo

2007-12-18 Thread Galkowski, Jan
From: Mulholland, Tom [mailto:[EMAIL PROTECTED] Sent: Monday, December 17, 2007 8:04 PM To: Galkowski, Jan Subject: RE: [R-sig-Geo] adding new maps to map, mapdata, mapproj;archive of R-sig-geo I have been playing with the iplots package and it uses the map class. So I wrote something to test the

[R-sig-Geo] adding new maps to map, mapdata, mapproj; archive of R-sig-geo

2007-12-17 Thread Galkowski, Jan
Hi, Sorry if this has been covered somewhere. Is there a writeup describing how new maps get added to map, mapdata, mapproj? What I'm trying to do is overlay polygons in latitude and longitude onto these maps, then project them. The associated things you'd like to do with such overlays apply,

Re: [R-sig-Geo] point.in.polygon() on massive datasets

2007-12-14 Thread Galkowski, Jan
Hi Markus, I have had success with a somewhat smaller point-in-polygon problem using a relational database (PostgreSQL) with suitable indexes to determine containment. The database serves as a large lookup table for this purpose, using RODBC. I use R to build a series of SQL INSERT statements to

[R-sig-Geo] voronoi/Delaunay/Dirichlet tessellation on sphere in R or S?

2007-11-27 Thread Galkowski, Jan
There's Renka's STRIPACK, and TRIPACK, respectively, ACM TOMS Algorithms 772 and 751, and there's the R package "deldir" which does the Delaunay for a plane, but does anyone have or know of the tessellation in R for a sphere? Also, is there a standard indexing scheme for Delaunay facets, and perha