Re: [R-sig-Geo] opening a grd file (Surfer)

2008-10-09 Thread Alexander Brenning
Hi Anne, in case you are not mainly interested in importing the (potentially large) Surfer grid into R, but want to do raster calculations instead (e.g. using other SAGA modules), you might still be interested in the RSAGA way of converting it: You didn't specify the file name of the SAGA ou

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Roger Bivand
On Thu, 9 Oct 2008, Barry Rowlingson wrote: Hi All, I'm now very impressed with PostGIS and what it can do with spatial data. Given the easy setup with Ubuntu, the easy admin with pgAdmin3, and the great interface with Qgis, I'm won over totally. Of course I still need to read my spatial data

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Dylan Beaudette
On Thursday 09 October 2008, Edzer Pebesma wrote: > Barry Rowlingson wrote: > > 2008/10/9 Edzer Pebesma <[EMAIL PROTECTED]>: > >> Why is this better practice than having a single 100x61 table with > >> monthly rainfall in the usual columns and a geometry column? > > > > Because when I get another

Re: [R-sig-Geo] opening a grd file (Surfer)

2008-10-09 Thread Roger Bivand
On Thu, 9 Oct 2008, Anne Goarant wrote: Hello! Does anyone know how to open a .grd file (surfer type file) with R (or convert this kind of file in ascii)? I tried using the RSAGA library but failed. rsaga.geoprocessor("io_grid", 3, list( FILE = '10s70s0e150e.sgrd')) This are bathymetry data

Re: [R-sig-Geo] opening a grd file (Surfer)

2008-10-09 Thread Dylan Beaudette
On Thursday 09 October 2008, Anne Goarant wrote: > Hello! > > Does anyone know how to open a .grd file (surfer type file) with R (or > convert this kind of file in ascii)? I tried using the RSAGA library but > failed. > rsaga.geoprocessor("io_grid", 3, list( FILE = '10s70s0e150e.sgrd')) > This are

[R-sig-Geo] opening a grd file (Surfer)

2008-10-09 Thread Anne Goarant
Hello! Does anyone know how to open a .grd file (surfer type file) with R (or convert this kind of file in ascii)? I tried using the RSAGA library but failed. rsaga.geoprocessor("io_grid", 3, list( FILE = '10s70s0e150e.sgrd')) This are bathymetry data from the GEBCO website. Thanks Anne

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Edzer Pebesma
Barry Rowlingson wrote: 2008/10/9 Edzer Pebesma <[EMAIL PROTECTED]>: Why is this better practice than having a single 100x61 table with monthly rainfall in the usual columns and a geometry column? Because when I get another year's data, I have to change the table definition. SQL does

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Barry Rowlingson
2008/10/9 Edzer Pebesma <[EMAIL PROTECTED]>: > Why is this better practice than having a single 100x61 table with monthly > rainfall in the usual columns and a geometry column? Because when I get another year's data, I have to change the table definition. SQL doesn't like that[1]. Keep your orig

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Edzer Pebesma
Barry Rowlingson wrote: Hi All, I'm now very impressed with PostGIS and what it can do with spatial data. Given the easy setup with Ubuntu, the easy admin with pgAdmin3, and the great interface with Qgis, I'm won over totally. Of course I still need to read my spatial data into R (until I d

Re: [R-sig-Geo] GLS spatial model-- range error

2008-10-09 Thread jofre carnicer cols
Sorry, I have just solved the problem. There was a typographic error in the code. Sorry for this mistake, Jofre --- El jue, 9/10/08, jofre carnicer cols <[EMAIL PROTECTED]> escribió: > De: jofre carnicer cols <[EMAIL PROTECTED]> > Asunto: GLS spatial model-- range error > Para: r-sig-geo@sta

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Barry Rowlingson
2008/10/9 G. Allegri <[EMAIL PROTECTED]>: > Hi Barry, > just a curiosity. Do you simply need a cloropleth with rainfall > amount/rate coloring, or you have to do more complex operations? > Just simple colouring, but I will have to make sure the scales are consistent across the 60 maps, so I will

Re: [R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread G. Allegri
Hi Barry, just a curiosity. Do you simply need a cloropleth with rainfall amount/rate coloring, or you have to do more complex operations? 2008/10/9 Barry Rowlingson <[EMAIL PROTECTED]>: > Hi All, > > I'm now very impressed with PostGIS and what it can do with spatial > data. Given the easy setup

[R-sig-Geo] PostGIS and sp integration

2008-10-09 Thread Barry Rowlingson
Hi All, I'm now very impressed with PostGIS and what it can do with spatial data. Given the easy setup with Ubuntu, the easy admin with pgAdmin3, and the great interface with Qgis, I'm won over totally. Of course I still need to read my spatial data into R (until I do all my analysis in Python.

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Henk Sierdsema
Hi Rubén, Thanks for defining 'small' better. For me a dataset up to a thousand points is a small dataset, but for most people this is quite a big data set... Henk Henk Sierdsema SOVON Vogelonderzoek Nederland / SOVON Dutch Centre for Field Ornithology Rijksstraatweg 178 6573 DG Beek-Ubbe

Re: [R-sig-Geo] GLS spatial model-- range error

2008-10-09 Thread Edzer Pebesma
I don't I think the authors of nlme are on this list, but I looked at the help of corExp, and that suggests that using the argument correlation = corExp(form = ~x+y) will be more useful that what you did. -- Edzer jofre carnicer cols wrote: Hello, I am computing a gls spatial model: require

Re: [R-sig-Geo] plotting Ripley fonctions

2008-10-09 Thread Alexandre VILLERS
g the aggregation is or the attraction. thanks ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo __ Information from ESET Mail Security, version of virus signature database 3506 (20081009) __ The message was chec

[R-sig-Geo] GLS spatial model-- range error

2008-10-09 Thread jofre carnicer cols
Hello, I am computing a gls spatial model: require(nlme) gls.exp<-gls(S1 ~ X2TPall + tmngem + tmaxjlm + arbdenkm + prat_herm + arbclakm + conreukm, data=snouter.df, correlation=corExp(form~X+Y), na.action=na.omit) and obtain the following error: Error en Initialize.corSpatial(X[[1L]], ..

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Alexander Brenning
Hi, landcover classification with four classes... though I haven't experimented with spatial logistic regression, I would expect that there is more information in the spectral satellite data (feature space) that you are probably using as predictors than in the adjacency/distance information.

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Henk Sierdsema
Hi Ivan, I have uploaded two papers from Tom Hengl et al on this subject for you. This should give you a good idea how to solve the problem. www.sovon.nl/temp/Hengl1.pdf www.sovon.nl/temp/Hengl2.pdf Henk Henk Sierdsema SOVON Vogelonderzoek Nederland / SOVON Dutch Centre for Field Ornitholog

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Rubén Roa-Ureta
Henk Sierdsema wrote: Hi Ivan, Can you tell me what the purpose is of your modelling? Is it simply producing spatial predictions based on a logistic model or do you want to incorporate spatial autocorrelation in the models? Given your last mail it seems you want to incorporate spatial autocor

[R-sig-Geo] plotting Ripley fonctions

2008-10-09 Thread Lassana Touré
Hi. I need some help for plotting L(r) and Lcross so that I can have my confidence interval(CI) parallel to X axis. The plot I have obtain so far the graphic is not horizontal, then it is difficult to interpret strong the aggregation is or the attraction. thanks ___

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread ivan valencia
Hi HENK My data comes from land use classification, a grid with 1km2 resolution. I have a bynary classification with 4 types of land use, for each land use type have to run a logistic regression with available covariables also at the same raster format. I want to consider spatial logistic regressi

Re: [R-sig-Geo] CODE for spatial logistic regression

2008-10-09 Thread Henk Sierdsema
Hi Ivan, Can you tell me what the purpose is of your modelling? Is it simply producing spatial predictions based on a logistic model or do you want to incorporate spatial autocorrelation in the models? Given your last mail it seems you want to incorporate spatial autocorrelation despite the fac