[R-sig-Geo] adding ids and df to hexagonal polys ?

2009-09-09 Thread Horacio Samaniego
Hi, I am trying to understand the different data types in the sp package. However, I do not seem to get the way that spdf are built. For example, let say that I can build a hexagonal grid. That work just fine using > regio <- readShapePoly("regions.shp") > hpts <- spsample(regio,type="hexagon

[R-sig-Geo] Classes for spatio-temporal data?

2009-09-09 Thread Blair Christian
Hi All, I'm starting to work with spatio temporal data, and wondered if there were any extensions to the sp package that include a temporal component. (for example, any classes that might contain 3 objects/data frames, one for the spatial information, one for the temporal information, and one for

[R-sig-Geo] Intersection of two coordinates (triangulation)

2009-09-09 Thread Breitbach, Nils
Hi Community, I am new to this forum and I am starting to analyse my telemetry data trying to be independent from commercial software. So far all I wanted to do was possible with R. Now I am looking for a function that is able to do triangulation. The structure of my data is as follows: - two

Re: [R-sig-Geo] get raster pixel x,y points from lat/lon?

2009-09-09 Thread milton ruser
Hi David, I think it is possivel, but if you will deal with large images, I hardly suggestion you handle this on GRASS, and if you need play with R, you use spgrass6() package to connect R on grass. good luck milton On Wed, Sep 9, 2009 at 1:19 PM, David Mireles wrote: > I'm new to rgdal and g

[R-sig-Geo] get raster pixel x,y points from lat/lon?

2009-09-09 Thread David Mireles
I'm new to rgdal and gdal-- I have a very large nitf image (too big to fit in memory) that I would like to extract subset areas of interest based on a center lat/long coordinate (and possibly an elevation). I would like to use the getRasterData function extract the AOI (which requires pixel x,y of

Re: [R-sig-Geo] decimate points from a Lidar file

2009-09-09 Thread Jonathan Greenberg
Michael: Typically Lidar files come "tiled", e.g. you get a set of point clouds pre-subsetted to, say, 1km x 1km window. These subsetted data are usually a LOT easier to write some batch scripts to deal with than trying to analyze the entire database at once. Have you seen if the source

Re: [R-sig-Geo] Errror with gmeta6()

2009-09-09 Thread Agustin Lobo
I've found the problem: I had both libgdal1-1.5.0-grass and libgdal1-1.6.0-grass If I keep libgdal1-1.6.0-grass only and reinstall rgdal and spgrass6, things go back to normality. The problem is that libgdal1-1.5.0-grass is automatically installed by the qgis grass plugin. So without libgdal1-1.

Re: [R-sig-Geo] Errror with gmeta6()

2009-09-09 Thread Roger Bivand
On Wed, 9 Sep 2009, Agustin Lobo wrote: Hi! I've started spgrass6() with no problems (I get: GRASS GIS interface loaded with GRASS version: 6.4.0RC4 and location: carlos but then get an error with gmeta6: gmeta6() g.region: error while loading shared libraries: libgrass_vect.so: cannot open

[R-sig-Geo] Errror with gmeta6()

2009-09-09 Thread Agustin Lobo
Hi! I've started spgrass6() with no problems (I get: GRASS GIS interface loaded with GRASS version: 6.4.0RC4 and location: carlos but then get an error with gmeta6: > gmeta6() g.region: error while loading shared libraries: libgrass_vect.so: cannot open shared object file: No such file or dire

[R-sig-Geo] introducing first order autocorrelation functions in GLMM to account for temporal autocorrelation

2009-09-09 Thread jofre carnicer cols
Dear All, I would like to ask which might be the best manner to introduce a first order autocorrelation function in a GLMM model in order to account for temporal autocorrelation effects in a temporal series dataset, using package Lme4 (function lmer, Laplace approach) or glmmPQL functions (penal

Re: [R-sig-Geo] Producing ENFA derived suitability maps with adehabitat

2009-09-09 Thread Tomislav Hengl
Dear Alexandre, I think that this should be doable. Once you run the ENFA to estimate the habitat model, copy the values of new rasters (future habitat) and replace the original values in the "" object. Then simply re-run the predict.enfa method e.g.: # prepare the dataset for ENFA: > beidata

[R-sig-Geo] creating a grid from a contour vector shape

2009-09-09 Thread Alexandre VILLERS
Thank you for your messages. I wasn't very clear in my post but both answers will be usefull... Best regards Alex -- Alexandre Villers PhD Student Team "Biodiversity" Centre d'Etudes Biologiques de Chizé-CNRS UPR1934 79360 Beauvoir sur Niort Phone +33 (0)5 49 09 96 13 Fax +33 (0)5 49 09 65 2