n
example with 'meuse'!
Very best,
David
--
David Holstius
PhD Student in Environmental Health Sciences
UC Berkeley School of Public Health
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Enviado desde mi HTC
- Mensaje original -
De: r-sig-geo-requ...@stat.math.ethz.ch
Enviado: Miércoles, 15 de Septiembre de 2010 05:30 a.m.
Para: r-sig-geo@stat.math.ethz.ch
Asunto: R-sig-Geo Digest, Vol 85, Issue 15
Send R-sig-Geo mailing list submissions to
r-sig-geo@stat.math.e
Enviado desde mi HTC
- Mensaje original -
De: r-sig-geo-requ...@stat.math.ethz.ch
Enviado: Miércoles, 18 de Agosto de 2010 05:30 a.m.
Para: r-sig-geo@stat.math.ethz.ch
Asunto: R-sig-Geo Digest, Vol 84, Issue 18
Send R-sig-Geo mailing list submissions to
r-sig-geo@stat.math.ethz.
ith presence/absence:
image(tppc01, col=c("white", "green"))
map('world', add=T)
points(coordinates(tr), t="l", col="blue")
Some ideas ?
Many thanks for your time.
David
--
***
David PINAUD
Ingéni
Or it is a problem of your GRIDS argument. How is it organized?
> >> > qp <- list.files(getwd(), pattern=".sgrd")
> >> > filelist <- list.files(getwd(), pattern=".sgrd")
> >> > filelist
> >> [1] "g1.sgrd" "g2.
tim.haer...@lwf.bayern.de
Date: Monday, July 19, 2010 11:16 am
Subject: RE: [R-sig-Geo] RSAGA - rsaga.geoprocessor crash
To: david.de...@queensu.ca, r-sig-geo@stat.math.ethz.ch
> Hi David,
>
> > I am using SAGA v2.0.2 and RSAGA interface with R 2.9.2.
>
> Maybe it is a problem w
experience with this process?
Thanks in advance,
David Depew
Postdoctoral Fellow
School of Environmental Studies
Queen's University
Kingston, Ontario
K7L 3N6
david.de...@queensu.ca
P: (613) 533-6000 x77831
[[alternative HTML version deleted]]
_
2.11.1
> get("SYS", envir = .GRASS_CACHE)
[1] "WinNat"
Am I missing something? Any idea?
Thank you very much for your help and patience.
David
--
***
David PINAUD
Ingénieur de Recherche "Analyses spatiales"
Hi,
What R package provides a way to determine if a longitude/latitude
coordinate is situated inside a closed, complex spherical polygon (i.e., a
point inside a polygon on Earth's surface)?
For example:
p <- rbind( c(-123.4447,48.5146), c(-123.3171,48.5128),
c(-123.2607,48.4108), c(-123.4018,48.
tes of the main polygon ("continental France")
pf <- list(x=cf[,1], y=cf[,2]) # list of coordinates, as dp() needs a
list and not a matrix or dataframe...
cf1 <- dp(pf, 0.1) # simplification, with a bandwith of 0.1 decimal degree
plot(fr) # to see the result
points(cf1, col=&qu
Hi Michael,
maybe you should try the function dp() in the package shapefiles that is
an implementation of the Douglas-Peucker polyLine simplification algorithm.
HTH
David
Michael Friendly a écrit :
The Guerry package contains two maps of france (gfrance, gfrance85)
which are quite detailed
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
usage. – Oikos 102: 367–377.
Unfortunately, I don't know a R implementation to do this directly. One
(quick and dirty) way could be to weight the output of the Kernel
estimation by a raster representing distance to coast (or ocean/continent)?
Hope it helps
David
Tim Clark a écrit :
Dear
of the polygon. Could you
suggestme any package that implements this function or that could help
me to compute 'l' and 'w'?
Thank you,
--
David March Morlà
Doctorando / PhD candidate
Fisheries-GIS specialist
Instituto Mediterraneo de Estudios Avanzados (UIB-CSIC)
C/Miquel Marqués,
Hi
I got the same error trying to do something similar. For me the
problem was with the polygon data. One of the shapes was corrupted. My
tip is to try and identify what polygon renders the error and confirm
that it is OK.
-David Wahlund
On Wed, Mar 18, 2009 at 16:01, Wesley Roberts wrote
As of the last time I compiled code like the following (17-June-2008),
it displayed an indicator map:
require(sp)
data(meuse)
# indicator: which observations have high Zn content?
meuse$ind <- (meuse$zinc > 300)
table(meuse$ind)
coordinates(meuse) <- ~x +y
class(meuse) # SpatialPointsDataFrame
sp
On Thu, Jan 15, 2009 at 12:07, Edzer Pebesma
wrote:
>
>
> David Wahlund wrote:
>>
>> Hi,
>> trying to extract by mask using sp package.
>> Was surprised that I couldn't find 'overlay(x =
>> "SpatialPixelsDataFrame", y = "SpatialPolyg
eem
to be getting the index values of the pixels or something. The raster
is float, as seen from the summary below. But the return from
overlay() is ordered integer.
Am I missing something?
Thank you
David Wahlund
> summary(asc)
Object of class SpatialPixelsDataFrame
Coordinates:
min max
f all the point pairs). Any suggestions would be appreciated.Dave Warner David WarnerResearch Fishery BiologistUSGS Great Lakes Science Center1451 Green RoadAnn Arbor MI 48105734.214.9392
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.
Michael
adding fill = TRUE (so that polygon is used) works for me.
David
David Maxwell
Cefas
Lowestoft Laboratory
Suffolk
NR33 OHT
UK
www.cefas.co.uk
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED
A beta online version of this software can be found here:
http://www.gapminder.org/world/#$majorMode=chart$
rgds,
David
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Agustin Lobo
Sent: Wednesday, August 06, 2008 4:55 PM
To: sig-geo; martí
-data:when_to_use
http://wiki.r-project.org/rwiki/doku.php?id=tips:stats-spatial:kriging
and see if you can replace the blank spaces with something useful
Cheers
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
Maybe this would be an appropriate time to mention
http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data:when_to_use
How would list members answer this question?
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
2008/7/3
The spRbind methods also complain about non-unique IDs. Is there any
other difference between them and rbind?
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
2008/7/1 Roger Bivand <[EMAIL PROTECTED]>:
> On Tue, 1 Jul 20
created.
Is there a simple way to give new IDs to these objects, and to know
how many I have to give?
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
___
R-sig-Geo mailing list
R-sig-Geo
Yes, and a list of links to external tools would be great also.
Work so far is at
http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data . It's
very incomplete so please feel free to improve!
David Hugh-Jones
PhD Candidate
Essex University Department of Government
e who know a lot can provide answers,
and if (like me) you know barely anything, then you'll know what
questions need answering.
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
2008/7/1 Roger Bivand <[EMAIL PROTECTED]>:
Roger
Thanks for that tip. As it happens, what I am doing is quick and dirty
enough that it doesn't matter, but it could have bitten me fairly
badly.
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
2008/7/1 Roger Bivand &l
r to get
at the actual raw data, and this is so complicated and confusing that
it would be easier just to work with the raw data.
OK, I'll stop venting. If there's anything I could do to improve this
situation, I would gladly try.
David Hugh-Jones
PhD Candidate
Essex University Departmen
Thanks David for his comment about dp.
Quick question: is there any reasonably comprehensible API
documentation for the "sp" package? I have just spent about an hour
trying to get a list of points from a SpatialPolygons object. I
eventually just printed everything out and found the da
maybe you can try the function dp() in the package "shapefiles", which
is an implementation of the Douglas-Peucker polyLine simplification
algorithm.
Hope it helps
David
David Hugh-Jones a écrit :
Hi all
I have a big dataset of points and am doing stuff on them that takes a
lot o
e
objects? (I assume that there's no way to convert _to_ Arcinfo, though
I know it's possible to read from it).
Cheers
David Hugh-Jones
PhD Candidate
Essex University Department of Government
http://davidhughjones.googlepages.com
___
R-
This might be of interest too:
GOOGLE CHARTS FROM R: MAPS
http://www.iq.harvard.edu/blog/sss/archives/2008/04/google_charts_f_1.shtml
rgds,
David
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Barry Rowlingson
Sent: Friday, June 13, 2008 11:08 AM
To
nt mapped outputs is limited I plan to calculate the width of the
confidence interval on the 0,1 scale and map this.
Thanks again, this list is an excellent catalyst for learning
David
[EMAIL PROTECTED]
-Original Message-
From: Edzer Pebesma [mailto:[EMAIL PROTECTED]
Sent: 09 April 200
prediction (s0) = glm prediction (s0) + kriged residual prediction (s0)
for location s0
Is assuming independence of the two components adequate?
var rk(s0) ~= var glm prediction (s0) + var kriged residual prediction (s0)
Thanks again,
David Maxwell & Vanessa Stelzenmüller
[EMAIL PROTECTED]
--
r didn't do correctly?
Is there another way to set this directory?
Thank you for your help^.
David
--
***
David PINAUD
Ingénieur de Recherche "Analyses spatiales"
Centre d'Etudes Biologiques de Chizé - CNRS UPR1934
79360 Vil
I am sure Edzer can answer this, but I post it to the list because it
may be of general interest.
If I compute an indicator:
> data(meuse)
> meuse$i.Pb.med <- meuse$lead < median(meuse$lead)
I can of course compute an indicator variogram and krige at new
locations:
> vi <- variogram(i.Pb.med ~1
n, but this
> is not possible for dim_station, as there are only 3 dimensions. Do I have
> to redefine 51 other dimensions? I tried to use list() to define the other
> dimensions, but that didn't work.
>
> Many thanks
> Thomas
>
>
>
>
>
> -Ursprünglich
are the pixels for my area
>
> Thanks for any help!
> Thomas
>
>
> __
>
> Thomas Szegvary
> Institute of Environmental Geosciences
> Department of Geosciences
> University of Basel
> Bernoullistrasse 30
> CH - 4056 Basel
>
> Tel. 41-61-267 04 82
> Fax. 41-61-267 04 79
> Email: [EMAIL PROTECTED]
> www.radon.unibas.ch
> www.unibas.ch/environment
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
---
David W. Pierce/ Climate Research Division
Scripps Inst. Oceanog. / (858) 534-8276 (voice)
[EMAIL PROTECTED] / (858) 534-8561 (fax)
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma wrote:
> David, you may also notice that the map produced does not have a cell
> centre or cell 'crossing' at (0,0). Of course you'd assume enough
> cleverness that this would be enforced automatically, but too bad. Try
>
> v <- variogram(Ni ~ 1,
I've found a quite puzzling behaviour in gstat's variogram method, with
the map=T argument. I believe the result should be a symmetric map
because all the semivariances are squared differences. This is indeed
the case if I try an example from the Meuse set:
> library(gstat)
> data(meuse)
> coordin
Imagine if everybody on R-help automatically sent "out of the office"
replies to the list. R as we know it would cease to exist!
> Today's Topics:
>
>1. Distmap: the way I could do it (Agustin Lobo)
>2. Crossing native ArcGis GRID with a XY coordinate table
> (Milton Cezar Ribeiro)
Dear Agus,
Maybe you can try this:
dist<-distmap(yourdata)
dist.grid<-as.SpatialGridDataFrame.im(dist)
Regards
--
David March Morlà
Doctorando / PhD student
Instituto Mediterraneo de Estudios Avanzados (UIB-CSIC)
C/Miquel Marqués, 21
07190 Esporles
Islas Baleares - Spain
Tel: +34 971 6
(http://www.faocopemed.org/vldocs/028/publi04.pdf). Have
you got any suggestion to perform these path distance maps with R?
Thanks,
--
David March Morlà
Doctorando / PhD student
Instituto Mediterraneo de Estudios Avanzados (UIB-CSIC)
C/Miquel Marqués, 21
07190 Esporles
Islas Baleares - Spain
e the
minimum", still work? I can't see why not but my intuition is not
clear.
Cheers
David
On 16/06/07, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Fri, 15 Jun 2007, David Hugh-Jones wrote:
>
> > Hi R-spatialists,
> >
> > I'm working with a big list of
extract lines from polygons?
Cheers
David
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
rend surfaces', with the
predict(loess(z ~x+y,...),expand.grid(...)...) functions.
Another thing you might watch out for is that atmospheric model data often
has staggered grids, with some features defined at the grid points, and
others at the element centers.
Dave
--
Dr. David Forr
.
Regards,
--
David
On 6/2/07, Dirk Eddelbuettel <[EMAIL PROTECTED]> wrote:
>
> Brian,
>
> On 2 June 2007 at 19:08, Prof Brian Ripley wrote:
> | There are quite a few 'db API' features in RODBC that are not in DBI.
>
> Random idea of the day: Would it be wortwhile
On Wed, 7 Feb 2007, Roger Bivand wrote:
> On Wed, 7 Feb 2007, David Forrest wrote:
>
>> Hi,
>>
>> I've got a mis-projected UTM zone 18 polygon shapefile, the coordinates of
>> which need a transformation of (+287.4395, -77.5542) added to them to be
>> in
e function that would add this offset to the component
polygon coordinates, bbox, and labpts data?
Dave
--
Dr. David Forrest
[EMAIL PROTECTED](804)684-7900w
[EMAIL PROTECTED] (804)642-0662h
On Mon, 5 Feb 2007, Roger Bivand wrote:
On Fri, 2 Feb 2007, Virgilio Gómez-Rubio wrote:
Hi David,
...
The reason is that Sourceforge split out projects into separate addresses
last year to aid administration. I have updated the information on the
website, and added notes where appropriate
s -z3 -d:pserver:[EMAIL PROTECTED]:/cvsroot/r-spatial co -P
sp
They seem inconsistent with the instructions at
http://sourceforge.net/cvs/?group_id=84357 which seem to work.
Dave
--
Dr. David Forrest
[EMAIL PROTECTED](804)684-7900w
[EMAIL PROT
On Thu, 18 Jan 2007, Edzer J Pebesma wrote:
> David Forrest wrote:
>
>> On Wed, 17 Jan 2007, Roger Bivand wrote:
...
>> Hmm. I was thinking that the spplot() was the best-practices method for
>> spatial data.
>>
> spplot lets you choose class intervals and colou
file. It seems to
be working pretty well, and I'll see if he thinks its ready for sharing.
Dave
--
Dr. David Forrest
[EMAIL PROTECTED](804)684-7900w
[EMAIL PROTECTED] (804)642-06
On Wed, 17 Jan 2007, David Forrest wrote:
> On Wed, 17 Jan 2007, Roger Bivand wrote:
...
> I feel really close with this, except for the key.
Answering my own question by looking through the source at
http://r-spatial.cvs.sourceforge.net/r-spatial/sp/R/spplot.R?view=markup I
see that
color argument to fig2dev, which
> needs investigation.
>
> Roger
>
>>
>> If anyone has lots of time to spare, a SpatialGridDataFrame device with
>> red, green, blue and alpha channels would be nice (or a GDAL MEM device,
>> or whatever).
>>
>> Roger
>
On Wed, 17 Jan 2007, Roger Bivand wrote:
> On Tue, 16 Jan 2007, David Forrest wrote:
>
>> On Thu, 11 Jan 2007, Roger Bivand wrote:
>>
>>> On Thu, 11 Jan 2007, David Forrest wrote:
>>>
>>>> Hi All,
>> ...
>>
>> I think I can make
On Thu, 11 Jan 2007, Roger Bivand wrote:
> On Thu, 11 Jan 2007, David Forrest wrote:
>
>> Hi All,
...
I think I can make marginless plots in base, lattice, and sp graphics now,
but I'm still missing something about the data ranges for spplot.
#
On Thu, 11 Jan 2007, Roger Bivand wrote:
> On Thu, 11 Jan 2007, David Forrest wrote:
>
>> Hi All,
>>
>> I'd like to make graphic files of plots with controlled pixel->coordinate
>> mapping. Ultimately, I'd like to get the graphs out into Google earth
margin/border around my graphic that is hard to
correlate to the bbox.
Any hints?
Dave
--
Dr. David Forrest
[EMAIL PROTECTED](804)684-7900w
[EMAIL PROTECTED] (804)642-0662h
On Wed, 10 Jan 2007, Stephane Dray wrote:
> Dear David,
> you can simply add the line :
>
> deb http://cran.R-project.org/bin/linux/debian stable/
>
> to your sources.list. You can then obtain newer R versions for stable Debian.
Perfect. That seems exactly like what I w
e is a production server, and I'd like to not mix the
Debian releases on it if I can avoid it. Its using some rusty old
scripts with library(maps) on it.
Would you advise a separate R 2.4.0 + install for using sp on these
machines?
Dave
--
Dr. David Forrest
[EMAIL PROTECTED]
Hello Roger,
Yes, that's very helpful.
Thanks very much,
Dave Martell
On Mon, 14 Aug 2006, Roger Bivand wrote:
> On Mon, 14 Aug 2006, David Martell wrote:
>
> >
> > I am using the lm.morantest procedure in spdep to test for spatial
> > autocorrelation in my l
Martell
_
David L. Martell, Professor
Faculty of Forestry, University of Toronto
33 Willcocks Street, Toronto, Ontario, Canada, M5S 3B3
E-mail: [EMAIL PROTECTED]
URL: http://www.firelab.utoronto.ca
Phone: (416) 978-6960
Fax
n the small example data set,
all that I really need is just the two numbers 4 and 11 (before I had
said 3, but I realized the third is unnecessary). From these I plan to
construct my classes as <4,>=4 and <11, >=11 when I create my
classification in MapServer where I am doing the plotting.
a<-c(11,1,3,3,4,4,4,4,4,4,5)
This gives me the vector (3,1,1,1,2,2,2,2,2,2,2)
ideally I want the output for this to be just (3,5,11)
Taking this with the min/max of my data I can create the classes
1-3,4-5,6-11
On 2/24/06, Roger Bivand <[EMAIL PROTECTED]> wrote:
> On Fri, 24 Feb 2006, D
after?
Any help is appreciated,
David
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
/Views/Spatial.html
or http://sal.uiuc.edu/csiss/Rgeo/ to demonstrate the utility of
R-sig-Geo?
Dave
--
Dr. David Forrest
[EMAIL PROTECTED](804)684-7900w
[EMAIL PROTECTED] (804)642-0662h
http
) # works
col.whitebg() returns a list of light-backgrounded trellis settings that
work together.
>
> Francisco
>
>
>
> >From: David Forrest <[EMAIL PROTECTED]>
> >Reply-To: David Forrest <[EMAIL PROTECTED]>
> >To: r-sig-geo@stat.math.ethz.ch
> >Subje
und?
Also, the later examples on http://r-spatial.sourceforge.net/gallery/ seem
to fail by not finding SpatialPolygons() Should those examples be
revised?
meuse.sr =
SpatialPolygons(list(Polygons(list(Polygon(meuse.riv)),"meuse.riv")))
Dave
--
Dr. David Forre
70 matches
Mail list logo
