Re: [R-sig-Geo] geoR questions

2009-06-12 Thread Rubén Roa Ureta
María Eugenia Utgés wrote: Original Message Subject: Presentation and a "few" Questions (geoR) From:María Eugenia Utgés Date:Wed, June 10, 2009 5:21 pm To: r-sig-geo@stat.math.ethz.ch -

Re: [R-sig-Geo] geoR questions

2009-06-12 Thread Paulo Justiniano Ribeiro Jr
María I will try to answer some of your questions. Ohters will depend on the intended resulkts and may not be suitable for discussions hers. Hello Everyone I am making my doctorate in biology at the University of Buenos Aires, Argentina. I am working on survival and dispersal of the south ame

[R-sig-Geo] geoR questions

2009-06-12 Thread María Eugenia Utgés
Original Message Subject: Presentation and a "few" Questions (geoR) From:María Eugenia Utgés Date:Wed, June 10, 2009 5:21 pm To: r-sig-geo@stat.math.ethz.ch --

Re: [R-sig-Geo] [SoC v.autokrige2.py] autoKrige() and projected data

2009-06-12 Thread Anne Ghisla
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Edzer Pebesma ha scritto: > If newdata is created automatically (and this makes sense), it should > recieve the projection information of data. hello all, @Paul: I added the projection of obj to new_data before return. Please test if it works as part

Re: [R-sig-Geo] convert point data to count data

2009-06-12 Thread Adrian.Baddeley
Sorry, that should have been df <- data.frame(x=as.vector(rasterx.im(Z)), y=as.vector(rastery.im(Z)), count=as.vector(as.matrix(Z))) Adrian ___ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https:/

[R-sig-Geo] convert point data to count data

2009-06-12 Thread Adrian.Baddeley
Karla Falk writes: > I have a data set of points that I would like to convert to grid cell counts. > Ultimately I would like to end up with a table that has an x-y coordinate for > each grid cell, > and the count of the # of points falling in each grid cell. Does anyone know > a way to do thi

[R-sig-Geo] [Report #3] Kriging with GRASS and R: v.autokrige port to wxPython and more

2009-06-12 Thread Anne Ghisla
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hello all, this weekly report will be very short, as yesterday I sent this email [0] and made little progress. I'm having some problems in using autoKrige() function work, AFAIK because of projection information handling [1]. best regards and happy

Re: [R-sig-Geo] [SoC v.autokrige2.py] autoKrige() and projected data

2009-06-12 Thread Edzer Pebesma
If newdata is created automatically (and this makes sense), it should recieve the projection information of data. On another issue: what happens if both are provided but have different projection attributes? And what happens if data has long/lat coordinates (is.projected(data) == FALSE)? In that

[R-sig-Geo] [SoC v.autokrige2.py] autoKrige() and projected data

2009-06-12 Thread Anne Ghisla
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hello Paul, Edzer, the integration of automap package in v.autokrige2 is quite completed [0], thanks Paul for this clean interface to gstat! Sadly, I'm getting an error when autoKrige() receives no new_data argument. Examples in documentation all invo

Re: [R-sig-Geo] [Report #2] Kriging with GRASS and R: v.autokrige port to wxPython and more

2009-06-12 Thread Anne Ghisla
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Thomas Adams ha scritto: > Anne, > > Will this be scriptable? I guess it's not a huge deal since it would be > using a script that included GRASS and R using automap… Sure it will! I'll have a look at how to add this functionality in the next weeks.

[R-sig-Geo] Vedr.: Re: plotting a polygon with holes

2009-06-12 Thread Bjarke Christensen
That works like a charm: > myDEM.pix <- as(myDEM, 'SpatialPixelsDataFrame') > idx <- overlay(myDEM.pix, islands) > plot(oceans, col="blue") > plot(islands, col="white, add=T) > image(myDEM.pix[!is.na(idx),], add=T) Thank you very much, once again. Bjarke Christensen.