Hamid Ghorbani writes:
> In spatstat and for calculating the maximum absolute difference
> between the empirical Ripley's K-function (\hat(K), isotropic edge
> correction) and the theoretical K-function for a simulated Poisson
> process in a fixed window (suppose 3D case and fix number of points
Dear researcher,
Doing CSR test (based on Ripley K-function) for an observed
homogeneous point process in a Window W with intensity \lambda(=n/a) ,
we accept the null hypothesis of CSR if sup|\sqrt(\hat{K(t)}/\pi
-t|)<=c for t<=t_0 (c is unknown) which yields
max(0,(t-c)^2<=\hat{K}(t)<=\pi(c+t)^
Dear all,
I have a problem converting a large binary sparse weight matrix (62025 *
62025) of class dgTMatrix to listw.
First I create a sparse weight matrix with the code:
###weightsparse###
Wmatrix<-as_dgRMatrix_listw(Reg255.listw)
Since I want to model spatial autocorrelation i
Sorry, it keeps going to my spam folder so I skipped the subject
Hi:
Could someone show me how to add state borders to a shapefile.
Please follow the bellow link to download it.Thanks
https://secure.filesanywhere.com/fs/v.aspx?v=897263875a6472a99baa
library(rgdal)
library(maptools)
dsn="C:/patht
On Tue, May 11, 2010 at 10:53 AM, Roger Bivand wrote:
> On Mon, 10 May 2010, Stephen G. Eick wrote:
>
>> Are there any easy function to parse WKT and convert it into a sp object?
>> E.g. to parse
>>
>
> Not yet, not in OGR in rgdal anyway. There is a draft pair of MULTIPOLYGON
> <-> SpatialPolygo
On Mon, 10 May 2010, Stephen G. Eick wrote:
Are there any easy function to parse WKT and convert it into a sp object?
E.g. to parse
Not yet, not in OGR in rgdal anyway. There is a draft pair of MULTIPOLYGON
<-> SpatialPolygons functions in rgeos on R-forge - maybe add a wish
there, or cons
Are there any easy function to parse WKT and convert it into a sp object?
E.g. to parse
"MULTILINESTRING((-88.06354 41.80337,-88.06415 41.80297,-88.06447
41.80288,-88.06566 41.80275,-88.06586 41.80277,-88.06596 41.80284,-88.06598
41.8029,-88.06592 41.80433,-88.06592 41.80449,-88.06598 41.8047
On Tue, 5 Jan 2010, laure velez wrote:
Dear list member,
Although the answer to this might be straightforward I'm struggling with this
issue for weeks now :
I have a shapefile (only one polygon) and a raster (.tif). I would like to know
the value of each pixel on the raster which overlap wi
Dear list member,
Although the answer to this might be straightforward I'm struggling with this
issue for weeks now :
I have a shapefile (only one polygon) and a raster (.tif). I would like to
know the value of each pixel on the raster which overlap with the polygon in
the shapefile.
I ha
Hi Wesley,
The proj.4 string definition is available at (bookmark it!):
http://trac.osgeo.org/proj/wiki/GenParms
and then a list of projections is available at:
http://www.remotesensing.org/geotiff/proj_list/
all the best from Split,
Tom Hengl
http://spatial-analyst.net/
Roger,
I got the package. This is great. Thank you very much, I will have a test
run very soon. Have a nice weekend!
X.W.
On Thu, Apr 9, 2009 at 2:42 PM, Roger Bivand wrote:
> On Thu, 9 Apr 2009, x wong wrote:
>
> I do use Windows, and thank you very much for the help!
>>
>> I know how to do t
On Thu, 9 Apr 2009, x wong wrote:
I do use Windows, and thank you very much for the help!
I know how to do this in GIS as Jonathan suggested; however, I haven't tried
it in R with the solution Rick recommended since i am still cleaning my
input data.
The problem of this task is very obvious: i
I do use Windows, and thank you very much for the help!
I know how to do this in GIS as Jonathan suggested; however, I haven't tried
it in R with the solution Rick recommended since i am still cleaning my
input data.
The problem of this task is very obvious: it will be super inefficient if
all th
On Wed, 8 Apr 2009, x wong wrote:
Thank you very much, that is very helpful.
x.w.
On Wed, Apr 8, 2009 at 12:47 PM, rick reeves wrote:
Assuming that this dataset fits within the confines of R, this could be
done with sp package: spdistsN1() function
See attached sample solution:
http://n
Thank you very much, that is very helpful.
x.w.
On Wed, Apr 8, 2009 at 12:47 PM, rick reeves wrote:
> Assuming that this dataset fits within the confines of R, this could be
> done with sp package: spdistsN1() function
> See attached sample solution:
>
>
> http://nceas.ucsb.edu/scicomp/GISSem
Assuming that this dataset fits within the confines of R, this could be
done with sp package: spdistsN1() function
See attached sample solution:
http://nceas.ucsb.edu/scicomp/GISSeminar/UseCases/AssignClosestPointsR/AssignClosestPointsR.html
Hope this helps! RR
Jonathan Greenberg wrote:
It se
It seems that a raster-based approach would make more sense, rather than
the hugely computationally inefficient approach you are suggesting --
how about using a least cost path or euclidean distance approach? Both
are available in many GIS packages (ArcMap, GRASS GIS, etc...)
--j
x wong wrot
dear all,
I have a raster GIS map converted to points. Now, I am trying to calculate
the nearest neighbor distance for every point to a cluster of points within
this map. Since the number of points is more than two millions, I am
wondering whether there is a computation efficient way to do this.
Markus Loecher <[EMAIL PROTECTED]> writes:
> I am looking for a function that
> indicates whether a line segment falls inside a rectangle (i.e. crosses
> the rectangle boundaries somewhere).
In the spatstat package: Suppose W is the rectangle (or other polygon),
stored as an object of class "owin
Markus Loecher <[EMAIL PROTECTED]> writes:
> I am trying to compute a spatial density estimator where the samples are
> not points but lines. The usual kernel density estimators will not work
> very well, because of their fixed xy orientation. I am imagining adding up
> many "ridge" like density c
Dear Roger,
Regarding the shapefile question I think it was a memory problem because I had
88000 points and wen I reduced the number of points to half it worked fine.
Regarding the second question: thanks for the help and your solution almost
solves my problems, but I forgot to tell you that I
Hello,
In our lab we have an actual mapping problem and maybe there is an R
solution to it.
Aim is to calculate a surface from calculated data of fine root
biomass. Therefore a number of single trees (approx. 3000) have been
sampled. With a function from the literature it is possible to
es
On Tue, 21 Aug 2007, Sam Field wrote:
> List,
>
> I am looking for ways of estimating spatial autoregression models that adjust
> for a known source of heteroskedaticity and the Waller and Gotway (2004) text
> outline how this can be done in the case of the SAR model. If I work at it, I
> think I
List,
I am looking for ways of estimating spatial autoregression models that adjust
for a known source of heteroskedaticity and the Waller and Gotway (2004) text
outline how this can be done in the case of the SAR model. If I work at it, I
think I can implement this myself in R, but I wanted t
Yong Li wrote:
> require(maptools)
>
> #read in shape file
> pts <- readShapePoints(system.file("shapes/baltim.shp",
> package="maptools")[1])
>
> #define the dimension of grids
> x.min <- summary(pts)$bbox[1]
> y.min <- summary(pts)$bbox[2]
> x.max <- summary(pts)$bbox[3]
> y.max <- summary(pts)
require(maptools)
#read in shape file
pts <- readShapePoints(system.file("shapes/baltim.shp",
package="maptools")[1])
#define the dimension of grids
x.min <- summary(pts)$bbox[1]
y.min <- summary(pts)$bbox[2]
x.max <- summary(pts)$bbox[3]
y.max <- summary(pts)$bbox[4]
cellsize <- 20
x.n <- ro
On Wed, 11 Apr 2007, pir2.jv wrote:
> Merci de me désabonner de cette liste
Please *do* use the unsubscribe instructions given in the link at
the foot of each message - scroll to the foot of the page:
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
It's not hard to do, and gives the maintainer
Merci de me désabonner de cette liste
Cordialement
Jacques Vernin
___
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
On Thu, 15 Mar 2007, Roger Bivand wrote:
> On Thu, 15 Mar 2007, Adrian Baddeley wrote:
>
> >
> > On Tue, 13 Mar 2007, Michael Friendly wrote:
> >
> > > I've read the documentation, but can't find an example of how
> > > to specify a circular window for the ppp object that would contain the
> >
On Thu, 15 Mar 2007, Adrian Baddeley wrote:
>
> On Tue, 13 Mar 2007, Michael Friendly wrote:
>
> > I've read the documentation, but can't find an example of how
> > to specify a circular window for the ppp object that would contain the
> > data and the window.
>
> To create a circular window i
On Tue, 13 Mar 2007, Michael Friendly wrote:
> I've read the documentation, but can't find an example of how
> to specify a circular window for the ppp object that would contain the
> data and the window.
To create a circular window in spatstat, you can use the
function disc(). See help(disc).
On Fri, 1 Dec 2006 [EMAIL PROTECTED] wrote:
> Hello R-sig-geo list,
>
>I've been trying to figure this out by myself but i am getting kind of
> lost so i was hoping if you could point me to the right direction.
>I have a shape file for map of berkeley--the one in
> california--formatted i
Hello R-sig-geo list,
I've been trying to figure this out by myself but i am getting kind of
lost so i was hoping if you could point me to the right direction.
I have a shape file for map of berkeley--the one in
california--formatted in long/lat coordinates. I also have set of
locations in l
CTED]
> Office -- 803-777-5729
> Fax -- 803-777-4972
> http://people.cas.sc.edu/hardistf/
>
>
> On 9/28/06, Frank Hardisty <[EMAIL PROTECTED]> wrote:
> >
> >
> > -- Forwarded message --
> > From: "pir2.jv" <[EMAIL PROTECTED]&g
EMAIL PROTECTED]
> Date: Thu, 28 Sep 2006 16:42:13 +0200
> Subject: Re: [R-sig-Geo] (no subject)
>
> Le 28 sept. 06 à 16:18, Roger Bivand a écrit :
>
> >>>
> >>> If you say ls(), do the names of any of the objects begin with plot?
> >>
> >>
Le 28 sept. 06 à 16:18, Roger Bivand a écrit :
>>>
>>> If you say ls(), do the names of any of the objects begin with plot?
>>
>>> ls(pat="plot")
>> [1] "BBplot" "Bplot" "draw.plotmath.cell"
>> [4] "plot" "plotG" "plotM"
>> [7] "plotx"
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> Le 28 sept. 06 à 15:06, Roger Bivand a écrit :
>
> >>> By the way, why the load()?
> >>
> >> The load is the system response when I click the R.data in the
> >> project directory
> >
> > I think that this is the problem, and that you have an alternative
Le 28 sept. 06 à 15:06, Roger Bivand a écrit :
>>> By the way, why the load()?
>>
>> The load is the system response when I click the R.data in the
>> project directory
>
> I think that this is the problem, and that you have an alternative
> local
> copy of plot that disturbs the method dispa
On Thu, 28 Sep 2006, Roger Bivand wrote:
> On Thu, 28 Sep 2006, pir2.jv wrote:
>
> >
> > Le 28 sept. 06 à 13:09, Roger Bivand a écrit :
> >
> > >
> > > OK, the file size at least looks correct, and the output of the
> > > summary
> > > looks good too. By the way, why the load()?
> >
> > The
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> Le 28 sept. 06 à 13:09, Roger Bivand a écrit :
>
> >
> > OK, the file size at least looks correct, and the output of the
> > summary
> > looks good too. By the way, why the load()?
>
> The load is the system response when I click the R.data in the
> p
Le 28 sept. 06 à 13:09, Roger Bivand a écrit :
>
> OK, the file size at least looks correct, and the output of the
> summary
> looks good too. By the way, why the load()?
The load is the system response when I click the R.data in the
project directory
>
> Could you now try plot(xx, axes=TRUE
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> Le 28 sept. 06 à 12:13, Roger Bivand a écrit :
>
> > (Reverting to list)
> >
> >
> > So possibly there is an endian issue, or a contaminated shapefile.
> >
> > In a fresh R session, please do:
> >
> > library(maptools)
> > sessionInfo()
> > fn <- system.fi
Le 28 sept. 06 à 12:13, Roger Bivand a écrit :
> (Reverting to list)
>
>
> So possibly there is an endian issue, or a contaminated shapefile.
>
> In a fresh R session, please do:
>
> library(maptools)
> sessionInfo()
> fn <- system.file("shapes/sids.shp", package="maptools")[1]
> file.info(fn)
>
(Reverting to list)
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> Le 28 sept. 06 à 11:03, Roger Bivand a écrit :
>
> > On Thu, 28 Sep 2006, pir2.jv wrote:
> >
> >>
> >> I have a misfit whan I want to see the gallery:
> >>
> >> demo(gallery)
> >> . all is OK at the beginning
> >
> > Nothing
On Thu, 28 Sep 2006, pir2.jv wrote:
>
> I have a misfit whan I want to see the gallery:
>
> demo(gallery)
> . all is OK at the beginning
Nothing is OK. You clearly have a malfunctioning system, at least one of
the packages is wrongly installed. You have not included version
informatio
I have a misfit whan I want to see the gallery:
demo(gallery)
. all is OK at the beginning
but:
> library(sp)
> library(maptools)
> nc2 <- readShapePoly(system.file("shapes/sids.shp",
package = "maptools")[1], proj4string = CRS("+proj=longlat
+datum=NAD27"))
> plot(nc2)
On Wed, 27 Sep 2006, Edzer J. Pebesma wrote:
> Jacques,
>
> the projection functions are now in package rgdal which is on CRAN so
> should have Mac binaries. Look for the function spTransform.
>
Sorry, no automatic OSX binaries, because of dependence on external
libraries PROJ.4 and GDAL. Mac
Jacques,
the projection functions are now in package rgdal which is on CRAN so
should have Mac binaries. Look for the function spTransform.
Also, an email I sent to you last night answering a question related to
sp bounced; please resend the question to this list.
--
Edzer
pir2.jv wrote:
> I t
I thinh spproj is not available for Mac Users ?
> rSpatial = "http://r-spatial.sourceforge.net/R";
> install.packages("spproj", repos=rSpatial)
Avis : unable to access index for repository http://r-
spatial.sourceforge.net/R/bin/macosx/powerpc/contrib/2.3
Avis dans download.packages(pkgs, dest
I have been asked whether this list is willing to accept postings with
information on vacant positions in spatial data analysis. Like the main r-help
list, I feel that short postings from list subscribers giving such information,
with links to further details on a website, are positive, in the s
Hi folks,
After using spdep to construct a neighbor list from a SpatialPolygons
object, I plotted the neighbor list and noted that I needed one more link
that wasn't caught by the poly2nb command. I found the edit.nb() function
and can successfully use the interactive graphical interface to add th
You'll notice that I wasn't volunteering.. :-)
I guess my main comment re massive datasets, etc. is that once I learned
about the power of generic programming (we've been adapting the Boost
Graph Library to work with GIS), I always feel a great sense of waste
when I see nice algorithms chained pit
> -Original Message-
> From: Edzer J. Pebesma [mailto:[EMAIL PROTECTED]
> Sent: Friday, April 07, 2006 5:30 AM
> To: Tim Keitt
> Cc: Andy Bunn; r-sig-geo@stat.math.ethz.ch; [EMAIL PROTECTED]
> Subject: Re: [R-sig-Geo] (no subject); Clipping/Subsetting Sp objects
On Fri, 7 Apr 2006, Edzer J. Pebesma wrote:
> Andy, Tim, the magic phrases (following your example) are:
>
> clip.sp = SpatialPolygons(list(Polygons(list(Polygon(clip)), ID="clip")))
> fullgrid(x) = FALSE
> x.clip = x[!is.na(overlay(x, clip.sp)),]
> image(x.clip,col="blue",add=T)
>
> the first c
Andy, Tim, the magic phrases (following your example) are:
clip.sp = SpatialPolygons(list(Polygons(list(Polygon(clip)), ID="clip")))
fullgrid(x) = FALSE
x.clip = x[!is.na(overlay(x, clip.sp)),]
image(x.clip,col="blue",add=T)
the first command creates a SpatialPolygons object from
the coordinates.
We need R wrappers for:
http://www.vividsolutions.com/jts/jtshome.htm
or
http://geos.refractions.net/
My preference would be the latter.
THK
On Thu, 2006-04-06 at 16:38 -0400, Andy Bunn wrote:
> List,
>
> I want to clip (subset) an sp object SpatialGridDataFrame using a polygon.
> For instan
List,
I want to clip (subset) an sp object SpatialGridDataFrame using a polygon.
For instance in this example I want to create a new object of class
SpatialGridDataFrame that is clipped to the area inside the polygon on the
map.
data(meuse.grid)
coordinates(meuse.grid) <- ~x+y
gridded(meuse
Supose you have a data-frame "foo" which is something like:
CoordX CoordY Var Real
1 1 10 1
1 2 11 1
2 1 91
2 2 81
1 1 11 2
1 2 13 2
1 3 81
then you can read it as a geodata with:
foo.gd <- as.geo
Thank you for your reply.
unfortunately I haven't understood how to put the data in
the matrix and the realisations in the vector.
Could you be so kind to give me a little example? (I've
done a lot of tries but I didn't succeed in doing it).
Thank you so much for your help,
Michela
_
59 matches
Mail list logo
