My two cents.
Some of the best interface features I've found in Isatis, one of the
best commercial geostatistical sw I've used, are:
- interaction with the variogram map to select direction and plot the
ESPs together differentiating the various selected dirs with colors.
It permits to reveal anis
> With your permission, I will try to remove your message from the list
> archives.
Absolutely, I agree with removing it.
I'm really sorry for this. I've received the direct link to the
download, trusted as an open access series from Wiley... I shouldn't
have trusted it!
_
gt;
>> > Gregoire Dubois (PhD)
>> >
>> > DG Joint Research Centre - European Commission > Institute for
>> > Environment and Sustainability > Global Environment Monitoring Unit >
>> > Monitoring Of Natural resources for DEvelopm
Maybe it could be interesting for someone...
http://www.ebooksz.cn/2009/02/geostatistics-for-environmental.html
Giovanni
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
I just share a link to a paper I've recently read, about the thread
object: http://dx.doi.org/10.1016/j.csda.2004.06.015
giovanni
2008/12/11 Kamran Safi :
> Hi Corrado,
>
> A useful reference is: Diniz-Filho J.A.F. & Bini L.M. (2005). Modelling
> geographical patterns in species richness using ei
A correction.
cp["Z"] should be get(cp)[z], to evaluate the string to a variable name...
2008/11/19 G. Allegri <[EMAIL PROTECTED]>:
> Alessandro,
> I suggest to read some startup tutorials before asking in ML.
> Anyway you can do this:
>
> for(i in 1:500){
>
Alessandro,
I suggest to read some startup tutorials before asking in ML.
Anyway you can do this:
for(i in 1:500){
cp <- paste(sep="","centre_plot",i)
coordinates(cp) =~X+Y
proj4string(cp) <- CRS("+init=epsg:26911")
writeOGR(cp["Z"], paste(sep="",cp,".shp"), cp,driver="ESRI Shapefile")
}
2008/
> The maptools package has a function for computing the centroids of
> polygons - install that and try help.search("centroid").
x <- read.shape(system.file("shapes/sids.shp", package="maptools")
mycentroid <- get.Pcent(x)
___
R-sig-Geo mailing list
R-s
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
I suggest you also the following pages:
http://cran.r-project.org/web/packages/spgrass6/index.html
http://grass.itc.it/statsgrass/index.php
http://grass.itc.it/statsgrass/grass_geostats.html
A nice plugin from Carson Farmer for QGIS:
http://www.ftools.ca/manageR.html
Giovanni
2008/10/5 E. Anthon
I've just received "Applied spatial data analysis with R" from Springer.
I want to thanks publicly all the authors for this jewel. It will be
on my job desk from now on!
Giovanni Allegri
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://sta
I faced the problem of collecting informations about ArcGIS Geospatial
extension while I was following the geostatics course at university. A
month looking for documentation about what was behind the scene, but
nothing... Just basics explanations about kriging. That's when I've
discovered gstat!
I
ª Pilar Tugores Ferrà
>
> Becaria FPI - PhD Student
>
> Acústica de Pesquerías
>
> I.E.O. - Centro Oceanográfico de Baleares
>
> Muelle de Poniente s/n
>
> 07015 Palma de Mallorca (España)
>
> Tel.: (34) 971 401561
>
>
> -----Mensaje original-
> De:
> You spoke of more sophisticated methods of automatically choosing between
> models, what kind of methods did you have in mind?
Here's a list I gathered some time ago. Many of them are just buzzwords form me!
‧ Adjusted R-squared (Wherry 1931)
‧ Bootstrap (Efron 1979)
‧ Cross-validation (Stone
> At this stage I do
> this by computing the sums of squares between the model and the sample
> variogram and choose the one with the smallest SS. This is a rather crude
> way of selecting between the models.
Thanks Paul for automap. I'm planning to try it in the next occasion.
What are the major
kpss.test() from tseries package could be useful.
http://cran.r-project.org/web/packages/tseries/index.html
2008/8/7 Alessandro <[EMAIL PROTECTED]>:
> Hi All,
>
>
>
> is there un code to do a stationary test in R? because I have a data-base
> and I wish to know there is or not a stationary in my d
An Open Source version of GapMinder (completely developed from scratch
with Flex + AS3) is almost complete and coming out.
We're working on it... and will let you know when it will be ready :)
Giovanni
2008/8/6 David S. Bieri <[EMAIL PROTECTED]>:
> A beta online version of this software can be fo
Maybe this could help?
http://sekhon.berkeley.edu/library/sp/html/overlay-methods.html
2008/7/23 Rick Reeves <[EMAIL PROTECTED]>:
> Hi:
>
> Check out the FindPolys method in the PBSMapping package. here is a code
> snippet from a larger example:
> (I will be extracting a simpler use case for distr
http://grass.itc.it/grass63/binary/mswindows/native/
2008/7/21 Dylan Beaudette <[EMAIL PROTECTED]>:
> On Monday 21 July 2008, milton ruser wrote:
>> Dear all,
>>
>> I have a very large BINARY image (60,000x50,000 pixels) and need to
>> generate the distance map from the 1 values.
>> When the image
Roger didn't say it but a good comprehensive book about spatial data
handling in R is coming out in the next 2-3 months:
"Applied Spatial Data Analysis with R"
http://www.springer.com/public+health/epidemiology/book/978-0-387-78170-9
I've already boked it...
giovanni
2008/6/26 Roger Bivand <[EM
Thanks Roger.
I'll try it
2008/5/9 Roger Bivand <[EMAIL PROTECTED]>:
> On Mon, 5 May 2008, Roger Bivand wrote:
>
>> On Mon, 5 May 2008, G. Allegri wrote:
>>
>>> Thank you Roger. I thought to buy it to realize a course on spatial data
>>> analy
Thank you Roger. I thought to buy it to realize a course on spatial data
analysis and to suggest it to the alumns. I think it could be a good trade
off between hard mathematical theory and an R manual...
I've tried to insert the token indicated in the Discount Info, but Springer
returns that it's e
A question to Dr. Roger Bivand:
I'm about to order you book "Applied Spatial Data Analysis with R", but I
would like to have a look to the table of contents. It's not available on
the Springer page, could it be shared on the ml?
Does the book cover some teoretical foundaments about the arguments tr
Thanks Edzer. It's a sparse grid, indeed. I created it as SpatialPixels.
I won't care... :-)
2008/4/7, Edzer Pebesma <[EMAIL PROTECTED]>:
> G. Allegri wrote:
> > Hi list,
> > what does this warning mean?
> >
> > "Warning in points2grid(point
Hi list,
what does this warning mean?
"Warning in points2grid(points,tolerance):
grid has empty columns/rows in dimension 1"
It gets out when doing an overlay between points and grid
[overlay(grid, pointsdataframe)], with sp package.
I receive a complete list of values for the points, uploaded f
I answer myswelf:
> lines(gridlines(x,easts=pretty(bbox(x)[1,],n=20),norths=pretty(bbox(x)[2,],n=20)))
to create 20 lines. I didn't know the pretty() function...
2008/4/6, G. Allegri <[EMAIL PROTECTED]>:
> Is it possible to manage the lines spacing returned from the
> g
Is it possible to manage the lines spacing returned from the
gridlines() inside sp package?
Giovanni
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Are you sure you have the data stored in AUX? ArcGIS Desktop uses an AUX
file for images that do not support storing coordinate system information
internally. The image header and AUX file tells ArcMap about the coordinate
system including units and projection.
It happens that AUX files are created
Intereseting work...
I haven't had the time yet to investigate on it but, as I'm not a C++
developer, I was thinking to couple CGAL and R through thei respective
Python bindings for alpha shapes processing. Obviously, it would be
better having a low-level binding... but there's no time to learn
eve
yet.
>
> Yes that is what CGAL stands for.
>
> Murray
>
>
>
>
>
>
> Hi Murray,
>
> Regarding the first question, you can read into R first and then output
> it again. You could also call ogr2ogr using a system call like:
>
> system("ogr2ogr norma
system command which would be a faster option I guess.
>
> Thanks for you input.
>
>
> Murray
>
>
>
>
> G. Allegri wrote:
> > Hi Murray.
> > About the first question. As far as I understand, you're working in R
> > on a *imported* PostGIS
Hi Murray.
About the first question. As far as I understand, you're working in R
on a *imported* PostGIS layer (I suppose you use rgdal to do it,
right?). Then you would need to export it to a shapefile (easily done
with writeOGR). Why are you considering the use of ogr2ogr? As you
probably know,
Hello everyone.
I'm writing a chapter of my thesis about open source software for
spatial analysis. I've looked for documented case histories and
scientific literature about the use of R in the fields of lanscape
analysis, soil science, environmental analysis, geology, etc.
Yet I haven't found anyt
Roger, I have to get experience yet... I thought that the right way to
access the datas in Spatial*DataFrame should be through the use of
slots.
Thanks for the hint. Now things get easier! :-)
2008/3/12, Roger Bivand <[EMAIL PROTECTED]>:
> On Wed, 12 Mar 2008, G. Allegri wrote:
>
>
I'm trying to change a column name of the data slot of a SpatialPixelsDataFrame.
I've added the column this way:
> slot(EceZKRIG0_3,"data")[,index]<-Sim_mean_approx[[2]]
The column added is V151 (as 151 is the actual index).
I want to rename it, and I've tried:
> names(slot(EceZKRIG0_3,"data")[i
spsample! I looked inside the sp manual a lot of times, without seeing it!
:-)
Thanks Paul
2008/3/6, Paul Hiemstra <[EMAIL PROTECTED]>:
>
> G. Allegri schreef:
>
> > Hello list,
> > I need to create a SpatialPixelDataFrame to do a gaussian simulation
> over an
> &
Hello list,
I need to create a SpatialPixelDataFrame to do a gaussian simulation over an
area delimited by a polygon. The object I want to produce is something
similar to the "meuse.grid" example.
I could do it in Grass, converting first the polygon to raster (with the
resolution of the Grass regio
lues to do back-t?
2008/2/26, G. Allegri <[EMAIL PROTECTED]>:
> I'm sorry for asking again, but I haven't found an exhaustive answer yet.
> Is there a package or a function in R to execute normal score
> transformation (and back?).
>
> I know GSLIB is an alternati
I'm sorry for asking again, but I haven't found an exhaustive answer yet.
Is there a package or a function in R to execute normal score
transformation (and back?).
I know GSLIB is an alternative, but I would prefare to make the whole
work inside R...
Thanks,
Giovanni
I think it's better to wrap R scripting and void passing unuseful objects
back and forth between Python and R. In your script it doesn't seem that you
need direct interaction between Python and R (at least in this part of your
script).
I usually adopt something like this:
r("""CairoPNG(file="/var/
I have a question about using GLS estimation within the Regression
Kriging framework.
In Rossiter and Hengl texts it is stated that it makes not so much
difference, in many practical situations, using OLS rather then GLS.
I'd like to test it in my work.
What's the more feasible way to adobt it in R
Dear Edzer,
I've "medidated" on the answer you gave to Jose. Two considerations have raise:
1 - when you say that the approach of GLM is a way to consider
spatial dependence. I'm not sure about this. GLM are a way to account
for link functions between the dependent variables and covariates (ex.
P
Thanks, a really interesting post.
May I ask a stupid question? Where are the following function documented?
getPolygonCoordsSlot(getPolygonsPolygonsSlot(getSpPpolygonsSlot(poly.shp
)[[1]])[[1]])
They are in the sp package but I can't find references but in some posts on
this ML.
Thanks,
Giovann
I've solved using Cairo library. For now it's enough, then I'll try to
solve with Xvfb solution...
Thanks to all!
Giovanni
2008/1/16, G. Allegri <[EMAIL PROTECTED]>:
> I've tried using bitmap(), but it returns this taceback:
>
> [Wed Jan 16 16:10:05 2
GDD... I had troubles compiling thi package, and I can't solve it. It
tries to create some objects it cannot, but headers and dependencies
semm allright.
I try again
2008/1/16, Don MacQueen <[EMAIL PROTECTED]>:
> Have you looked into the GDD package?
>
> -Don
>
> At
ee
> ?bitmap and ?Cairo.
>
> Hope this helps,
>
> Paul
>
> G. Allegri wrote:
> > Thanks Marcelo,
> > I've tried using the suexec module in apache2 (it permits to change
> > userid and groupid on the base of the scripts called), but from
> > documenta
anaging the list at
> [EMAIL PROTECTED]
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of R-sig-Geo digest..."
>
>
> Today's Topics:
>
>1. Re: Spatially Constrained Clustering (Elias T. Krainski)
>
egression line.
>
>
> see also:
> FITTING DISTRIBUTIONS WITH R (by Vito Ricci)
> http://cran.r-project.org/doc/contrib/Ricci-distributions-en.pdf
>
>
> Tom Hengl
> http://spatial-analyst.net
>
>
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:
Hi everyone.
I'm sorry for the question maybe OT.
I'm trying to use R and Python to run some scripts via web interface.
I've successfully setup mod_python for Apache and the rpy module.
R needs X11 to use png() and jpeg() devices, so I have installed Xvfb
(X virtual framebuffer). I works correctly:
I'm trying to realize e regression kriging with gstat package on my
soil samples data. The response variable (ECe measuere) and covariates
appear positvely skewed.
As Tomislav Hengl suggests in its "framework for RK" [1], a logistic
transformation is proposed as a generic way to reduce the skewenes
Well, Berry, maybe I understand now what you mean. I was looking the problem
from R side, but what you will is reading and writing inside Qgis, without
invoking R shell, is it right?
In this case there's lot of work to do! :-)
Giovanni
2007/12/4, G. Allegri <[EMAIL PROTECTED]>:
&g
ROTECTED]>:
>
> G. Allegri wrote:
> > http://cran.r-project.org/src/contrib/Descriptions/rgdal.html
> >
> > Isn't enough?
> >
>
> That's what I'd use to write GDAL/OGR data to read into QGis, but if
> you want tighter coupling between QGis
http://cran.r-project.org/src/contrib/Descriptions/rgdal.html
Isn't enough?
Giovanni
2007/12/4, Barry Rowlingson <[EMAIL PROTECTED]>:
>
> Agustin Lobo wrote:
> > Hi!
> >
> > Does anyone has experience at displaying spatial R objects
> > with QGis (perhaps using the python console + Rpy or RSpyth
53 matches
Mail list logo
