Hi Dazhi,
The way to do this is use a point-in-polygon test. The sp package (which
maptools has as a dependency) has the function point.in.polygon which
performs a point in polygon test.
Dan
On 10/22/2010 06:40 AM, Dazhi Jiao wrote:
Hi,
I am new to spatial data analysis so sorry to ask a d
Hi Barry,
The last I looked (which was fairly recently) the sld plug-in for QGIS
has a very limited feature set. If the goal is to create an sld to go
with a shapefile, I'd argue for using uDig (another open source desktop
GIS product), which has really good sld support:
http://udig.refractio
sDataFrame object. If the
projection
# of the sp object is geographic (i.e., lon/lat) the areas are in
decimal degrees^2.
# Requires the sp and rgdal packages
# Programmer: Dan Putler
# Created: 06Nov09
# Modified: 06Nov09
areaPolygons<- function(spPoly, proj4string = NULL) {
if(cla
Roger,
This won't really matter since they are very close, but given the data
source and a North America based company, my guess is the underlying
datum is NAD83.
Dan
On 08/24/2010 10:45 AM, Roger Bivand wrote:
On Tue, 24 Aug 2010, Alexandre Villers wrote:
Hey,
There is an ESRI code
data set up of to 4000 points.
>
> Can anyone suggest ways to increase my memory limit or to modify my code to
> make this possible? I don't know what else to try!! Ideally, I would like to
> be able to make nsim=200, but I am trying
work with it in subsequent scripts
> that operate on this kc.edges.ll SpatialLinesDataFrame object?
>
> R, QGis, or even GRASS tools would be much appreciated.
>
> Thanks
> -Don
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
to use it in this way (to save
> on your local computer and as input to further analyses).
>
> Hadley
>
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
ternative to clustering is some sort of density algorithm
> that allows for value - but I can't quite get my head around how this
> might work.
>
> Could someone point me in the right direction - what other keywords
> should I be looking out for? what R packages are worth a
what program I should use and some tips for it.
> Thanks
> Juan
>
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
y have!
>
> Best
> Nils
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
eciate any pointers, thank you.
>
> Mark Na
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Coulmont
>
> This type of work is most efficiently done in a standard GIS. I would suggest
> installing GRASS GIS or something like that. This is a standard vector
> operation, and should be very fast. For this GRASS approach, see the
> v.overlay command.
>
> Cheers,
&g
s -- once in a
> SpatialLineDataFrame -- if I need to work not on a national level but
> on a municipal level ?
>
> Is it possible to create a subset of a smaller zone ?
> Coordinates:
> min max
> r1 2.560776 6.403847
> r2 49.4
it is likely a couple of
months away).
The function follows:
# dropLongLinks takes an nb (neighbor list) object, the original
coordinates,
# and a maximum distance, then removes links from the nb object for
links that
# exceed the maximum distance.
# Programer: Dan Putler
# Created: 08
I'd go with the shapefiles, the E00 files are ArcInfo export format
files, and it is very likely you would convert to shapefile sets along
the way.
Dan
On Sun, 2009-11-15 at 23:27 -0500, tom sgouros wrote:
> tom sgouros wrote:
>
> >
> > Dan Putler wrote:
>
>
sible make all these steps in R, and generate great maps,
> > but the most simple way is above. Personally, I use R to make the final maps
> > only when I have to produce a long series of simple maps.
> >
> > I think the best way is the integration of GIS software with R, each one
1.6.0 as an advanced install option. Eventually the SpatialLite
read/write capability for GDAL/OGR will make it to OSGEO4W, so it is
something of a moot point now.
Dan
On Sat, 2009-10-24 at 23:52 +0200, Roger Bivand wrote:
> On Sat, 24 Oct 2009, Dan Putler wrote:
>
> > Hi Edzer,
ialLite for the file-based exchange -- as the way forward. Do
> "standard" gdal/ogr installations on windows (whatever that means) come
> with support for sqlite and postgis, so that ogr2ogr takes you from one
> to another?
> --
> Edzer
>
> Barry Rowlingson wrote:
&g
as browsing such a
> >> file will show), and best avoided in its entirety. Simply saving a few
> >> hundred point coordinates may run you up to hundreds of Kb.
> >>
> >> Roger
> >>
> >>
> >>
> >>> Many thanks, Mark
> >>>
>
te any help you might be able to provide with this.
>
> Many thanks, Mark
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@stat.math.ethz.ch
> https://stat.ethz.ch/mailman/li
e are
a number of, and often available in online in PDF (start with An
Introduction to R, http://cran.r-project.org/doc/manuals/R-intro.pdf).
Now the answer to your question:
BestResult_R <- newData[newData$result == max(newData$result), ]
Dan Putler
On Tue, 2009-10-13 at 16:59 -0700, gianni
Sorry, I sent my initial reply only to Jim, not the list.
Dan
Forwarded Message
From: Jim Burke
To: Dan Putler
Subject: Re: [R-sig-Geo] How to negate %in%
Date: Mon, 25 May 2009 14:01:26 -0500
Thanks Dan,
Both yours and Torleif's suggestions work w
> >
> > Many thanks,
> > Wesley
> >
> >
> >
> > Wesley Roberts MSc.
> > Researcher: Earth Observation (Ecosystems)
> > Natural Resources and the Environment
> > CSIR
> > Tel: +27 (21) 888-2490
> > Fax: +27 (21) 888-2693
> >
ixels with the same value (i.e.,
> burned in my simulations). Could anyone point me towards such an
> algorithm, or help me get started finding one?
>
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
ist=c("bx","by","exampDat","exampGG","exampNB"))
Dan
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Hi Jim,
What are the nature of the attributes of the items, are they simply
actually distances from one another, or is their some other measures of
similarity/dissimilarity (say differences in soil types)?
Dan Putler
On Tue, 2008-10-28 at 13:35 -0700, Jim Bouldin wrote:
> I'm wondering
gt; > driver) shows no problems either, I'm afraid.
> >
> > Roger
> >
>
> version number:11.0
> vendor string:The X.Org Foundation
> vendor release number:10400090
> X.Org version: 1.4.0.90
>
> Sigh. Debugging X is such a pain...
>
> THK
>
> >
> >
> > >
> > > Roger
> > >
> > >>
> > >> THK
> > >>
> > >>>>
> > >>>> THK
> > >>>>
> > >>>>
> > >>>>>
> > >>>>> Roger
> > >>>>>
> > >>>>>>
> > >>>>>> THK
> > >>>>>>
> > >>>>>>
> > >>>>>
> > >>>>> --
> > >>>>> Roger Bivand
> > >>>>> Economic Geography Section, Department of Economics, Norwegian
> > School of
> > >>>>> Economics and Business Administration, Helleveien 30, N-5045 Bergen,
> > >>>>> Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
> > >>>>> e-mail: [EMAIL PROTECTED]
> > >>>>>
> > >>>>>
> > >>>>
> > >>>>
> > >>>>
> > >>>>
> > >>>>
> > >>>> --
> > >>>> Timothy H. Keitt, Assistant Professor
> > >>>> University of Texas at Austin http://www.keittlab.org/
> > >>>> Contact info and schedule at http://www.keittlab.org/~tkeitt/
> > >>>> Please support open document formats http://opendocument.xml.org/
> > >>>>
> > >>>
> > >>>
> > >>>
> > >>> --
> > >>> Timothy H. Keitt, Assistant Professor
> > >>> University of Texas at Austin http://www.keittlab.org/
> > >>> Contact info and schedule at http://www.keittlab.org/~tkeitt/
> > >>> Please support open document formats http://opendocument.xml.org/
> > >>>
> > >>
> > >>
> > >>
> > >>
> > >
> > >
> >
> > --
> > Roger Bivand
> > Economic Geography Section, Department of Economics, Norwegian School of
> > Economics and Business Administration, Helleveien 30, N-5045 Bergen,
> > Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
> > e-mail: [EMAIL PROTECTED]
> >
> >
>
>
>
--
Dan Putler
Sauder School of Business
University of British Columbia
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
tically efficient when there are large
> numbers of points. Accordingly, if the
> number of points in X exceeds the threshold nlarge, then only
> the border correction will be com-
> puted. Setting nlarge=Inf will prevent this from happening.
>
ses)
pointcase <- ppp(cases[,1], cases[,2], window=W) #generate the ppp
object
kdensity<-density.ppp(pointcase, 3000)
plot(kdensity, col=myPal)
rm(list=c("guichi","W","cases","pointcase","kdensity","myPal"))
___
Dan
On 25-Ju
Hi All,
To add some detail to Roger's earlier, post density.ppp in the
spatstat seems to be a very good answer to the original post since it
is specifically designed to estimate a kernel density for a point
process pattern. This function use a bivariate Gaussian smoother that
lends itself
Hi Markus,
There is a command line program called PAGC that does postal address
geocoding, which is open source (written to C). A link to it is:
http://sourceforge.net/projects/pagc/
Dan
On 6-May-07, at 5:31 PM, Markus Loecher wrote:
> Hello,
> is there an R package which would let me geocode
Denis,
If you are trying to do this within the context of using PBSmapping,
you might want to check out the PBSmodelling package as well since it
includes facilities to write dialog boxes using simple text files.
Plus, it is a Geo related package so it makes it appropriate for this
list ;-
32 matches
Mail list logo