Re: [R-sig-Geo] Problem with %IN% ==FALSE when data frames are equal

2009-07-13 Thread Roger Bivand
On Mon, 13 Jul 2009, Jim Burke wrote: Hi everyone I have a issue with %IN% using FALSE deep within my R code processing. It errors if the comparison data frames are equal. DESIRED WORK AROUND: I would like examine my %IN% data frame fields before submitting then to the erroring block of code.

[R-sig-Geo] Problem with %IN% ==FALSE when data frames are equal

2009-07-13 Thread Jim Burke
Hi everyone I have a issue with %IN% using FALSE deep within my R code processing. It errors if the comparison data frames are equal. DESIRED WORK AROUND: I would like examine my %IN% data frame fields before submitting then to the erroring block of code. Then appropriately trap this error. Any

Re: [R-sig-Geo] kml generation question

2009-07-13 Thread Matt Oliver
It turns out to be a pretty significant difference.try image(lon,lat,mcsst) image(lonr,latr,mcsst, add = T) On Mon, Jul 13, 2009 at 12:04 PM, Barry Rowlingson < b.rowling...@lancaster.ac.uk> wrote: > On Mon, Jul 13, 2009 at 4:32 PM, Matt Oliver wrote: > > I tried to polygon.you

Re: [R-sig-Geo] kml generation question

2009-07-13 Thread Roger Bivand
On Mon, 13 Jul 2009, Barry Rowlingson wrote: On Mon, Jul 13, 2009 at 3:26 PM, Matt Oliver wrote: Dear r-sig-geo, I am trying to generate a kml for an image I am producing using image(). After reading about kml generation, I'm unsure if I sould be using kmlOverlay() or writeOGR, or some other f

Re: [R-sig-Geo] kml generation question

2009-07-13 Thread Barry Rowlingson
On Mon, Jul 13, 2009 at 4:32 PM, Matt Oliver wrote: > I tried to polygon.you are right...too  slow > > I'm not sure how to resample to get a regular grid...using interp or interpp > would not be wise because of the size of the grid... A quick bodge would be to ignore the slight change in lat

Re: [R-sig-Geo] Slope and Aspect calculations in R

2009-07-13 Thread Ashton Shortridge
Hello Josh, and list, Figuring out how to link R with GIS may be extremely useful to you, not just for calculating slope and aspect, but also for other analyses in the future. There are also reasons you might not want to; I am one of those people who likes to write his own code. Consequently I

Re: [R-sig-Geo] kml generation question

2009-07-13 Thread Matt Oliver
I tried to polygon.you are right...too slow I'm not sure how to resample to get a regular grid...using interp or interpp would not be wise because of the size of the grid... On Mon, Jul 13, 2009 at 11:09 AM, Barry Rowlingson < b.rowling...@lancaster.ac.uk> wrote: > On Mon, Jul 13, 2009 at

Re: [R-sig-Geo] kml generation question

2009-07-13 Thread Barry Rowlingson
On Mon, Jul 13, 2009 at 3:26 PM, Matt Oliver wrote: > Dear r-sig-geo, > > I am trying to generate a kml for an image I am producing using image(). > After reading about kml generation, I'm unsure if I sould be using > kmlOverlay() or writeOGR, or some other function. > > This is the data > > http:/

[R-sig-Geo] kml generation question

2009-07-13 Thread Matt Oliver
Dear r-sig-geo, I am trying to generate a kml for an image I am producing using image(). After reading about kml generation, I'm unsure if I sould be using kmlOverlay() or writeOGR, or some other function. This is the data http://www.ceoe.udel.edu/cms/moliver/20071003.276.0237.n17.nc require(nc

Re: [R-sig-Geo] spTransform: wrong results??

2009-07-13 Thread Agustin Lobo
Roger, First at all, I think that the problem is solved. Note that I was setting the CRS in QGIS, thus it is unequivocally specified as EPSG 23031, then read into R through readOGR() Anyway, the problem is that I was setting the towgs84 field as: proj4string(wpcast300SPDF) <- CRS(paste(" +proj=u