Murray Richardson utoronto.ca> writes:
>
> It would be great to have CGAL bindings in R.
>
Agreed. This is something I've been contemplating as well.
THK
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r
On Sat, Mar 29, 2008 at 4:36 PM, Barry Rowlingson
<[EMAIL PROTECTED]> wrote:
> Tim Keitt wrote:
>
> > version number:11.0
> > vendor string:The X.Org Foundation
> > vendor release number:10400090
> > X.Org version: 1.4.0.90
>
into QGIS causes a crash.
THK
>
> Dan
>
>
>
> On Sat, 2008-03-29 at 12:36 -0500, Tim Keitt wrote:
> > On Sat, Mar 29, 2008 at 12:21 PM, Roger Bivand <[EMAIL PROTECTED]> wrote:
> > >
> > > On Sat, 29 Mar 2008, Roger Bivand wrote:
> > >
&
On Sat, Mar 29, 2008 at 12:21 PM, Roger Bivand <[EMAIL PROTECTED]> wrote:
>
> On Sat, 29 Mar 2008, Roger Bivand wrote:
>
> > On Sat, 29 Mar 2008, Tim Keitt wrote:
> >
> >> On Sat, Mar 29, 2008 at 10:25 AM, Tim Keitt <[EMAIL PROTECTED]> wrote:
> >
On Sat, Mar 29, 2008 at 10:25 AM, Tim Keitt <[EMAIL PROTECTED]> wrote:
>
> On Sat, Mar 29, 2008 at 10:01 AM, Tim Keitt <[EMAIL PROTECTED]> wrote:
> >
> > On Sat, Mar 29, 2008 at 8:41 AM, Roger Bivand <[EMAIL PROTECTED]> wrote:
> > > On Fri, 28 Mar
On Sat, Mar 29, 2008 at 10:01 AM, Tim Keitt <[EMAIL PROTECTED]> wrote:
>
> On Sat, Mar 29, 2008 at 8:41 AM, Roger Bivand <[EMAIL PROTECTED]> wrote:
> > On Fri, 28 Mar 2008, Tim Keitt wrote:
> >
> > > I'm trying to see if I can track down a bug in X
On Sat, Mar 29, 2008 at 8:41 AM, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Fri, 28 Mar 2008, Tim Keitt wrote:
>
> > I'm trying to see if I can track down a bug in X windows that causes
> > it to crash when plotting (I believe) very small filled polygons (or
>
I'm trying to see if I can track down a bug in X windows that causes
it to crash when plotting (I believe) very small filled polygons (or
it may have to do with large numbers of polygons). I first noticed
this running QGIS. Loading a PostGIS layer with ~3K polygons will
cause X to abort. I can load
On 10/25/07, Agustin Lobo <[EMAIL PROTECTED]> wrote:
> Hi,
>
> Is image() able to display a color composite?
The displayDataset() function should use the colormap. You will need
to open with GDAL.open().
THK
>
> I've made an 8 bit color composite (with another software), saved it as
> geotif and
Hi Peter,
I've been doing this in postgis. Once you get used to it, its not too
difficult. Load data with shp2pgsql (or use rgdal + RODBC or similar
if you want to stay in R) and then create a new table to hold the
results. After that you can do "insert into restab select
buffer(the_geom, dist) fr
On 8/14/07, Paul Hiemstra <[EMAIL PROTECTED]> wrote:
> Mike,
>
> I tried, and failed, to directly load the HDF5 files into R using rgdal.
> The next thing I tried was to convert the .h5 file to a GTiff using
> FWTools. This seemed to work, but when I tried to load the file into R
> through rgdal I
>
> > Date: Wed, 18 Jul 2007 10:09:34 -0500
> > From: "Tim Keitt" <[EMAIL PROTECTED]>
> > Subject: [R-sig-Geo] Group-wise error variance and autocorrelation
> > To: "R geostat list"
> > Message-ID:
> > <[EMAIL PROTECTED]>
>
I'm trying to specify a model where the fixed effects are
dep ~ group + group:indep
and the error variance and spatial autocorrelation vary by group. I've
tried lme but AFAICT it only supports heteroscedasticity and
autocorrelation within groups. Any pointers?
THK
--
Timothy H. Keitt, Universi
Any algebraic operation supported by R can easily be applied to a GDAL
image using rgdal -- size is virtually unlimited if you pull in data
as raster lines or tiles. The approach is to create an output image
and then use getRasterData and putRasterData to read chunks of data,
operate on them, and t
On 7/16/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Sun, 15 Jul 2007, Tim Keitt wrote:
>
> > On 7/14/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> >> On Sat, 14 Jul 2007, Michal Gallay wrote:
> >>
> >> > Dear R Users,
> >>
On 7/14/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Sat, 14 Jul 2007, Michal Gallay wrote:
>
> > Dear R Users,
> >
> > could you please advise me on doing map algebra with spatial grids? It's
> > the first time I am using spatial objects in R. I have imported an ascii
> > grid file and wanted
This is perhaps because the 'reg.finalizer' function in R did not
originally have the 'on.exit' option and so dangling handles would be
left at the end of an R session. I just checked and the 'on.exit'
option is now available, so we need to enable it in the code. That
will force closing of all GDA
I rather like locfit.
THK
On 6/24/07, zhijie zhang <[EMAIL PROTECTED]> wrote:
> Dear Friends,
> Except kernel2d(splancs) function, are there any other functions on kernel
> density estimation in point pattern analysis? I use the kernel2d(splancs)
> function to anayze my dataset, and the result s
? Was that what you wanted?
>
> All best
> Miha
>
>
> --- Tim Keitt <[EMAIL PROTECTED]> wrote:
>
> > Anyone know of code that will do bilinear interpolation from irregular
> > samples to another set of irregular locations (including points
> > pos
> On Tue, 5 Jun 2007, Tim Keitt wrote:
>
> > Anyone know of code that will do bilinear interpolation from irregular
> > samples to another set of irregular locations (including points
> > possibly outside the convex hull of the original points)? Akima
> > interpolates
Anyone know of code that will do bilinear interpolation from irregular
samples to another set of irregular locations (including points
possibly outside the convex hull of the original points)? Akima
interpolates to a regular grid AFAICT. I did find one package on CRAN
that does this will splines, b
On 6/2/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> On Sat, 2 Jun 2007, Tim Keitt wrote:
>
> > On 6/2/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> >> On Sat, 2 Jun 2007, Tim Keitt wrote:
>
> [...]
>
> >> > Edzer's
On 6/2/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> On Sat, 2 Jun 2007, Tim Keitt wrote:
>
> > My Rdbi (+Rdbi.PgSQL) package was written to be much simpler and
> > easier to use than either DBI or RODBC. Unfortunately, I do not have
> > time to maintain it, so
My Rdbi (+Rdbi.PgSQL) package was written to be much simpler and
easier to use than either DBI or RODBC. Unfortunately, I do not have
time to maintain it, so I'm not sure what state it is in. Either DBI
or RODBC should work fine for pushing tables back and forth. RODBC is
probably the most up-to-da
The getRasterData command in rgdal will return a subset of the image
given an rgdal object handle. I believe you can now also use simple
array [] operations as well.
THK
On 6/1/07, Andrew Niccolai <[EMAIL PROTECTED]> wrote:
> Greetings Spatial R users/developers,
>
> I was curious as to whether o
I would also recommend you take a look at RODBC for general purpose
database access. My impression is that it has received a lot more
maintenance attention lately.
I can't comment on the windows gdal binaries as I'm not too familiar
with that platform.
THK
On 5/30/07, Mike Leahy <[EMAIL PROTECTE
; > > > > Roger
> > > > >
> > > > >
> > > > > On Thu, 5 Apr 2007 [EMAIL PROTECTED] wrote:
> > > > >
> > > > > > Yes, for many uses that is my choice also. For the conterminous
> > US
> > &g
On 4/9/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Mon, 9 Apr 2007 [EMAIL PROTECTED] wrote:
>
> > For preview graphics and for large areas such as continents, large
> > countries, hemispheres, or the whole earth, spherical projections are
> > often adequate. I can provide some of the ones I h
> > > For our hemispherical application, because we were gridding the
> data, we
> > > wanted parallels of latitude to be parallel in the projected
> coordinate
> > > space, which we wouldn't get with the Lambert azimuthal.
> > >
> > > (See attached
Thanks. My application is not that demanding. Really, I just want it
to look reasonable. My plan is to lay out the postings in the
projected coordinates and then back transform into geographic
coordinates for analysis. I tried lots of projections and found
Lamberts Azimuthal Equal Area to be quite
Should have thought of that -- updating now.
THK
On 4/5/07, Edzer J. Pebesma <[EMAIL PROTECTED]> wrote:
> Tim Keitt wrote:
>
> >>spplot(beta, "1")
> >>
> >>
> >Error in `[.data.frame`(as(obj, "data.frame"), zcol) :
> >
> spplot(beta, "1")
Error in `[.data.frame`(as(obj, "data.frame"), zcol) :
undefined columns selected
I can't seem to fix this. Version: 0.9-10
THK
--
Timothy H. Keitt, University of Texas at Austin
Contact info and schedule at http://www.keittlab.org/tkeitt/
Reprints at http://
Anyone know of a particularly good map projection for showing all of
North and South America without too much distortion?
THK
--
Timothy H. Keitt, University of Texas at Austin
Contact info and schedule at http://www.keittlab.org/tkeitt/
Reprints at http://www.keittlab.org/tkeitt/papers/
ODF att
>> But this is a simpler problem. Are there cases with interior holes?
> >>
> >> I have been meaning to write something like this for hexbin for a
> >>
> > while.
> >
> >> There
> >> are many cases where it would be nice to find approxim
On 3/16/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote:
>
>
> > My OSGeo summer of code idea is to link the BGL to postgis. I wasn't
> > originally thinking about topology per se, but it would achieve much the
> > same thing.
>
>
> I've thought a bit about this project, too, as it butts up against
t this is a simpler problem. Are there cases with interior holes?
> >
> > I have been meaning to write something like this for hexbin for a
> while.
> > There
> > are many cases where it would be nice to find approximations to the
> > density contours
> > and a q
to write something like this for hexbin for a while.
> There
> are many cases where it would be nice to find approximations to the
> density contours
> and a quick and dirty way is to threshold the hexagon counts, find the
> hull and
> smooth the perimeter.
>
> Nicholas
>
> On
ve been out
> of the
> GIS world to long. It seems to me that this would be something easy to
> solve,
> just tedious iteration of the polygon coordinates and some
> triangulation.
>
> Nicholas
>
> > Date: Thu, 15 Mar 2007 10:49:23 -0500
> > From: "Tim K
Thanks for the pointer to unionSpatialPolygons. I solved hull problem,
but in a one-off hackish way. I should have looked through maptools a
couple of hours ago...
THK
On 3/15/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Thu, 15 Mar 2007, Tim Keitt wrote:
>
> > Is there an
-- Forwarded message --
From: Frank Warmerdam <[EMAIL PROTECTED]>
Date: Mar 15, 2007 11:29 AM
Subject: [OSGeo-Announce] OSGeo Accepted for Google Summer of Code 2007
To: [EMAIL PROTECTED]
OSGeo is pleased to announce that Google has accepted OSGeo as mentoring
organization for the
Is there an 'sp' function that takes a polygon as its argument and
returns a set of line objects corresponding to the arcs in the
polygon?
Or better yet, a function that given a set of polygons, returns the
hull? (ie the set of singleton arcs after applying the polys to arcs
function)
THK
--
Ti
Sounds like a good wiki-page. Is there a wiki for r-sig-geo?
THK
On 3/6/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> Not sure if this will become an issue but in case I need
> it to justify my toolset does anyone have a list of
> recognizable organizations that use R for geographic
> proce
5)
plot(xy, asp = 1, axes = F, pch = 19, xlab = NA, ylab = NA, cex = 0.5)
lapply(genPolyList(xy), polygon, xpd = NA, ypd = NA)
Now to figure out the sp stuff.
Cheers,
Tim
On 3/3/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Fri, 2 Mar 2007, Tim Keitt wrote:
>
> > Hi All,
>
Hi All,
I'd like to make a hexagonal lattice (stored in a shapefile, ogr
coverage or postgis geometry). Is the general approach to make a
series of independent hexagons and place them in the plane such that
adjacent edges sit directly on top of each other? Or can one build a
mesh out of connected
In addition to Rogers advise, you might take a peek at the functions
available in postgis. Once you have polys in the postgis database, you
can do quite a lot of geometry ops by eg "select intersection(poly1,
poly2) from tabname" and so on. There's some grunt work involved. Its
certainly not pointy
Rick,
Looks like GDAL supports HDF4/5 format so you should be able to read data
using rgdal. I've not tried it myself. Let me know what you find.
Cheers,
Tim
On 1/10/07, Rick Reeves <[EMAIL PROTECTED]> wrote:
>
>
> Greetings:
>
> This week I am working with HDF format files containing ocean colo
Yes, that looks right. AFAIK it is the same on Debian.
THK
On 12/17/06, Roger Bivand <[EMAIL PROTECTED]> wrote:
>
> On Sun, 17 Dec 2006, Patrick Giraudoux wrote:
>
> > Tim Keitt a écrit :
> > > Try installing libgdal1-1.3.1-dev which contains the header files and
ery* good shape, up to date, and have been
> this way in the last couple of years (thanks to the DebianGIS project).
> Anyway, we don't want to start a flame war, don't we? I was just
> suggesting the easiest way to get recent packages through a single
> command.
> pc
>
You mean Ubuntu lags behind Debian unstable. Debian unstable is constantly
being updated. This is either a good thing or a bad thing depending on your
perspective. I ran Debian unstable for several years and believe me the
quality of Debian packaging varies a lot. Most Debian devs are quite good,
b
Try installing libgdal1-1.3.1-dev which contains the header files and should
have gdal-config.
You may run into problems with incompatible library versions. Ubuntu updates
every 6 months, but the R world seems never to stop rolling along.
I've found the Ubuntu "prevu" utility extremely helpful as
Hi Folks,
Just though I'd let people know that there is an experimental branch
("THK") of rgdal in CVS at the sourceforge site. It passes "R CMD
check" on my machine.
Major changes are:
1) "new" is no longer used to open datasets. "openDataset" (or
equivalently "GDAL.open") is now the only way t
My best guess is that you need to add your projection to one of the
files in /usr/share/gdal or in the projection definition files shipped
with rgdal. Please send it upstream if you do.
THK
On 9/11/06, Zev Ross <[EMAIL PROTECTED]> wrote:
> Hi All,
>
> I think this is probably a simple one. I'm re
Jonathan,
This is what rgdal was designed to do. Here's an example,
ds1 <- GDAL.open("input.tif")
driver <- new('GDALDriver', 'GTiff')
ds2 <- new('GDALTransientDataset', driver, nrow(ds1), ncol(ds1), type
= "Float64")
for (row in nrow(ds1)) {
x <- getRasterData(ds1, offset = c(row - 1, 0), reg
I found the type size defs in Rinternals.h. They are
R_XDR_INTEGER_SIZE and R_XDR_DOUBLE_SIZE. I've added a check in rgdal
to ensure that type sizes match up.
THK
--
Timothy H. Keitt
http://www.keittlab.org/
___
R-sig-Geo mailing list
R-sig-Geo@stat.m
On 9/8/06, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Fri, 8 Sep 2006, Tim Keitt wrote:
>
> > One possibility is that the size of the GDAL data type is not the same
> > as the size of the R data type. This will definitely cause problems as
> > allocation is done b
One possibility is that the size of the GDAL data type is not the same
as the size of the R data type. This will definitely cause problems as
allocation is done by R and the pointer is cast to void when passed to
GDAL. My original C++ code simply assumes that both are 32-bit (or at
least equivalent
Ah, that is very useful. I'd always assumed the hard part is computing the
area of the triangles on the surface of a sphere/elipse, but perhaps that is
not so hard after all. Anyone have a reference to some simple formulas? The
other question is how much error is introduced if one computes areas af
Do you mean the flat-plane area or the area projected onto the surface of
the ellipsoid? In either case, the computation is likely non-trivial. It
makes much more sense to project to UTM and do the calculation there.
THK
On 5/31/06, Patrick Giraudoux <[EMAIL PROTECTED]> wrote:
>
> Dear Listers,
>
Andy,
Maybe this is confusing because the data really ought to be in a dataframe
with factors coding for treatments?
Here's one solution:
foo[,'z',] <- sapply(foo[,'z',], function(x) max(c(x, 0)))
THK
On 5/26/06, Andy Bunn <[EMAIL PROTECTED]> wrote:
>
> > > This is driving me nuts. I have a 3-
You'll notice that I wasn't volunteering.. :-)
I guess my main comment re massive datasets, etc. is that once I learned
about the power of generic programming (we've been adapting the Boost
Graph Library to work with GIS), I always feel a great sense of waste
when I see nice algorithms chained pit
We need R wrappers for:
http://www.vividsolutions.com/jts/jtshome.htm
or
http://geos.refractions.net/
My preference would be the latter.
THK
On Thu, 2006-04-06 at 16:38 -0400, Andy Bunn wrote:
> List,
>
> I want to clip (subset) an sp object SpatialGridDataFrame using a polygon.
> For instan
I believe this can be done in PROJ by simply dumping the points as ASCII
and then calling 'geod' on the command line.
THK
On Thu, 2006-03-09 at 06:49 -0600, Sebastian Luque wrote:
> Hello,
>
> I have searched the archives and other sources, but cannot find an R
> function that would calculate th
On Tue, 2006-03-07 at 17:39 +, Barry Rowlingson wrote:
> Tim Keitt wrote:
>
> > That's a rather odd aspect of S methods
>
> one of many :)
>
> I think this stems from the S3 way that what other languages call
> "methods of objects", R has as
On Tue, 2006-03-07 at 16:24 +, Barry Rowlingson wrote:
> Wouter Buytaert wrote:
>
> > The resulting R function doesn't seem very complicated. I'm a newbie
> > though, so all recommendations welcome:
> >
> > topidx <- function(map) {
> >
> > if(!(class(map) == "SpatialGridDataFrame"))
> >
Barry,
I hope you will share your code! We've been using QGIS and I have in
mind just such an interface to R.
We have ENVI in the lab and it is very powerful for image processing. It
would take many many person years to reproduce all of the advanced
algorithms in that package. (It does a lot of s
I'm holding out for the Boost Geospatial Library... :-)
THK
On Wed, 2006-02-08 at 21:57 +0200, Ari Jolma wrote:
> Tim Keitt kirjoitti:
> > On Tue, 2006-02-07 at 22:25 +0100, Roger Bivand wrote:
> >
> >> No trouble - I think this hits me when I'm hurrying and
On Tue, 2006-02-07 at 22:25 +0100, Roger Bivand wrote:
> No trouble - I think this hits me when I'm hurrying and I re-use a
> not-closed handle, filling it with another handle. I can't replicate now.
Hmmm... I just realized that creating a raster band object and then
closing the associated datase
I seem to remember writing a wrapper for the gdal copy function,
although none of that has be extensively tested. If you copy an existing
dataset to a writable dataset, you could then retrieve the data, operate
on it and write back to the new dataset. It should (?) then have all the
metadata copied
Its better than no package at all, right? I think this is a metadata
issue. If packages were keyworded a la Gentoo, users could mask out
packages that don't pass R package checks, etc. (Maybe that's already in
place. I haven't tracked closely recent changes in R.)
THK
On Tue, 2005-08-23 at 08:32
69 matches
Mail list logo