Here we go, I'm assuming each ring is its own simple object, and a few
other things.
id = rep(seq(1,3),each=12)
az = rep(seq(0,359,30),3)
## assume you meant to use all ids here?
d = data.frame(id,az)
d$dist = round(rnorm(36,10,3),1)
xp = function(azimuth, r){cos(((90-azimuth)*pi)/180)*r}
yp = f
I'm trying to access GRASS functionality via spgrass6 and a plr script.
The following works (when executed directly from within R)...
library(spgrass6)
initGRASS(gisBase = "/usr/local/grass-6.4.3svn/", home = tempdir(), gisDbase =
"/usr/local/share/grassdata/", location = "aust_dem_1sec_wgs84",
Great!
thanks!
Erin
Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodge...@uhd.edu
-Original Message-
From: Michael Sumner [mailto:mdsum...@gmail.com]
Sent: Sun 7/15/2012 9:28 PM
To: Hodgess, Erin
Cc:
Just in terms of R, you could check out support for the OGR GML driver
in rgdal: http://www.gdal.org/ogr/drv_gml.html and the package sos4R
which has GML support (though may be too specialized). There are links
to the OGC documentation on the GDAL/OGR page.
Cheers, Mike.
On Mon, Jul 16, 2012 at 1
Dear Rusers:
I am using gstat package and met a problem about using krige. I searched
the mailing list, but still cannot solve it.
My sample dataset "testdata" could be downloaded from
http://www.fdsr.cn/en/xzzq.aspx?KID=4 . It includes a screenshot from
ArcGIS and the sample data. My codes are l
Hello again.
I would like to investigate the plotKML package and was wondering if there is a
tutorial or article or vignette available, please.
Thanks,
Erin
Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodge.
Hello everyone!
Does anyone use GML, please? If so, could you recommend a good tutorial or
intro book, please?
I've used R to produce KML, but have never used GML.
thanks,
Erin
Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - D
G'day list,
I'm calculating cross-variograms (variables below are sumso2="sulphur dioxide
emissions" and intensity="rain intensity" - I have 2500 sample points) using
GSTAT in R:
g <- gstat(NULL, id = "sumso2", form = sumso2 ~ 1, data=d1)
g <- gstat(g, id = "intensity", form = intensity ~ 1, da
Hello:
I have a .txt file i'm reading in R with 3670016 nrows made up of (256
lat, 512 lon and 28 forecast time) and ncol=31(no of days in the month). I
read the lat/lon data and got that my area of interest is at point row=93
to row 128 for lat and row 30 to 1 & 257 to 286 for lon. My issue is how
If you sent code it did not survive the journey. See the vignette in
sp for how to build polygons from scratch, section 7.1.
library(sp)
vignette("intro_sp")
On Mon, Jul 16, 2012 at 10:12 AM, wrote:
>
>
> r-sig-geo listers:
>
> Can anyone provide code or suggest commands for c
It seems my data-generating code was cut off. Here's another try:
id = rep(seq(1,3),each=12)           # Generate polygon
IDs'
az = rep(seq(0,359,30),3)              # Generate
azimuth sequence
c = data.frame(id
[1],az)Â Â Â Â Â Â Â Â Â Â Â Â Â Â Â
r-sig-geo listers:
Can anyone provide code or suggest commands for converting a data
frame with UTMS to a SpatialPolygons (sp package) object? Below, I
provide code to generate a sample dataframe, with three ordered
sequences of twelve coordinate pairs each.
id
Dear list members,
I am running R 2.15.1 on a MAC OSX 10.7.
I receive an error message regarding gpclib (see below):
> library(rgdal)
> library(maptools)
> library(rgeos)
Loading required package: stringr
rgeos: (SVN revision 330)
GEOS runtime version: 3.3.2-CAPI-1.7.2
Polygon checking: TRUE
13 matches
Mail list logo