Jason,
I agree completely that obtaining the GRASS GIS book
(http://www.grassbook.org/ ) is the way to go. I don't agree that it's
necessary to learn GIS concepts first using ArcGIS or QGIS. There are
many books on R, including "Applied Spatial Data Analysis with R".
Tom
On 9/21/10 12:06 P
Martin,
Depending on what you are doing, you may want to consider the use of
GRASS GIS (http://grass.osgeo.org/) in combination with R. This has been
a successful strategy for me in the past, not using Geotiffs, however.
The advantage is that the Geotiffs can be imported into GRASS GIS format
ylan
On Thu, Apr 22, 2010 at 8:32 AM, Thomas Adams wrote:
All:
I am using spplot to automatically generate a 21-panel eps graphic within a
script where the titles for the individual panes are too large and they get
cut-off at the beginning and end of the string. I have not been able to find
how
All:
I am using spplot to automatically generate a 21-panel eps graphic
within a script where the titles for the individual panes are too large
and they get cut-off at the beginning and end of the string. I have not
been able to find how to fix the font size so the titles will fit. I've
tried
Pebesma wrote:
Thomas, perhaps:
data <- readRAST6(map_list[[1]])
as map_list is a list, and readRAST6 seems to want a character vector,
which is in the first list element of map_list.
Thomas Adams wrote:
All:
With sp & spgrass6 I have no problem reading a list of maps with
readRAST
All:
With sp & spgrass6 I have no problem reading a list of maps with
readRAST6 when I do the following:
data<-readRAST6(c("map1","map2","map3",…))
however, I need to do this dynamically within a script, for example:
> inputfile<-Sys.getenv("R_INPUT_FILE")
> map_list<-scan(inputfile,list("")
little experiments…
I think the moral of the story is to restrict the vector map to the area
one is interested in and not to be sloppy. I'm still trying to grasp the
subtleties of v.generalize…
Thanks again,
Tom
Dylan Beaudette wrote:
On Tuesday 13 April 2010, Thomas Adams wrote:
All:
I'm using GRASS 6.4 with R 2.9.2 and sp/spgrass6. When plotting a vector
file, read in from GRASS, using spplot, the drawing takes many minutes
to complete. Apart from the fact that the vector file is probably more
detailed than I need, I know time is spent drawing off-screen. I tried
sp
Adam,
Do you mean PRISM maps such as these: http://www.prism.oregonstate.edu/
? If yes, then the question would be, do you want to duplicate the PRISM
methodology or effectively get similar results? One could use universal
kriging using gstat or methods proposed by
Tom
Adam A wrote:
Hello
Maybe something like this:
http://www.ars.usda.gov/sp2UserFiles/ad_hoc/1200SpatialWorkshop/01VinyardOverview.pdf
- Original Message -
From: milton ruser
Date: Tuesday, August 25, 2009 7:50 pm
Subject: Re: [R-sig-Geo] Kriging
> Hi Kabeli,
>
> I never saw Brian Vinyard slides 40! :
eeds to put some
effort to bridge this gap (check also with Olaf Conrad).
Here is more discussion about the same topic from few months ago:
https://stat.ethz.ch/pipermail/r-sig-geo/2009-April/005418.html
cheers,
Tom Hengl
-Original Message-
From: r-sig-geo-boun...@stat.math.eth
Josh,
I use R & GRASS together on both MacOS X and Linux without problems.
It's very direct to write a shell script that calls both R & GRASS
within a single script. Writing a R script that calls GRASS may be more
tedious to do what you want. Writing a shell script that calls both may
be the
Josh,
I can not answer your specific question WRT R, but before you embark on
this, you may consider R with GRASS GIS — they play together pretty
nicely and the slope and aspect capabilities are part of GRASS.
Regards,
Tom
Josh London wrote:
Hello
We are looking to mimic, in R, the slope (
Anne,
Will this be scriptable? I guess it's not a huge deal since it would be using a
script that included GRASS and R using automap…
Tom
Anne Ghisla wrote:
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Hello all,
there are some important news about the module!
First of all, the interface
Guy,
I don't know the answer to your question. However, depending on what you
are trying to do, might it be possible to use GRASS GIS in place of what
you are doing with IDL in ENVI? GRASS plays quite nicely with R. You
probably can read ENVI image data into R, but you may need to go through
Paul,
Not to mention the huge cost of ESRI products and the fact that they run only
on the MS-Windows platform. R also integrates well with GRASS GIS. Personally,
I can not see any compelling reasons to go the ESRI route vs R/GRASS…
I feel a great debt of gratitude to the selfless developers of
Paul,
This is outstanding; thank you very much! BTW, I noticed that when I run
your example that I get:
Warning: singular model in variogram fit
…actually, 2 such warnings with the OK example, and 1 with the UK
example. Are you seeing these too? Otherwise, the results look fine…
Thanks again
Dylan,
I think a solution using GRASS can be found on pages 110-111 of "Open
Source GIS: A GRASS GIS Approach", 3rd Ed. The same material is covered
in the 2nd Ed. as well, where you use r.mapcalc to combine two rasters
and judicious use of MASKs; a conditional statement in r.mapcalc is the key
Dylan,
I think a solution using GRASS can be found on pages 110-11 of "Open
Source GIS: A GRASS GIS Approach", 3rd Ed. The same material is covered
in the 2nd Ed. as well, where you use r.mapcalc to combine two rasters
and judicious use of MASKs; a conditional statement in r.mapcalc is the key.
Thomas,
Thanks to you and Roger for replies to my problem. It looks that you may
be right about the /Developer/SDKs/MacOSX10.5.sdk. I -thought- I had
updated the developer tools, but apparently not. I don't have the DVD
install disk with me, so this evening I'll make the install of the
correct
List,
I am attempting to build spgrass6 on MacOS X 10.5.1 (Leopard) on a 2.33. GHz
Intel Core Duo MacBook Pro. Using the Mac OS X Cocoa GUI, I am trying to
install spgrass6 from source. I have gdal installed fine; when I run
gdal-config at a term window prompt, I get:
Macintosh-2:bin teaiii$ g
Edzer,
Not to deflate Cesar's project, but from our past email exchanges you
are aware that I have been very pleased with the standalone version of
your gstat program and feel that it interfaces wonderfully with GRASS
6.x; I would highly recommend this approach to anyone. Also, for those
who m
Markus,
It's *possible* something like the Grid Analysis and Display System
(GrADS) (http://www.iges.org/grads/) could do the trick if you tweaked
the data right — but that may be too much of a chore. Look at the image
here http://www.iges.org/grads/gadoc/tutorial.html down the page a bit.
Reg
Gabor,
You *may* be able to accomplish this using R easily (I'm sure it's
doable), but I'm do not know how to do the conversion. I would suggest
as an alternative to do your work using GRASS GIS
(http://grass.itc.it/). R interfaces very well with GRASS and making
your LL to UTM conversion (and
Thierry
I have used used both GSTAT and GSLIB a lot and from a usability point
of view, I much prefer GSTAT because of it's ability to interface with
with GRASS GIS. This is particularly true of the standalone version of
GSTAT. I began using GSLIB, and while I continue to make use of it, I
fin
Matt,
Sorry, I do not have a very good solution — I encountered this exact
problem and with the frustrations I had at the time, I reverted back to
R 2.2.1. Life is good again… BTW, I'm using Redhat Enterprise 3 Linux.
Tom
Matt Moehr wrote:
> I have always used the R spatial packages on my hom
ents explicitly, to be sure that
> they are being understood correctly inside the functions. Depending on the
> class of the data= arugment, as far as I recall different versions of the
> gstat package treat the locations= argument differently.
>
> Roger
>
>
> On Fri, 12 May
attributes:
ohrfc.dem
Min. : 0.0
1st Qu.: 183.0
Median : 244.0
Mean : 302.2
3rd Qu.: 335.0
Max. :1707.1
> names(dem)
[1] "ohrfc.dem"
Regards,
Tom
Roger Bivand wrote:
> On Fri, 28 Apr 2006, Thomas Adams wrote:
>
>
>> Roger,
>>
>> Your su
mps file also has 'z' values of elevations. So, is this what you
already understood? Converting 'dem' to a SpatialPixelsDataFrame seemed
to only leave me with the grid locations and not the elevation values —
is this right.
Thanks again for your help!
Regards,
Tom
Roger Bivand
files have very different structures.
I'm not sure where to go to from here.
Regards,
Tom
Roger Bivand wrote:
> On Thu, 27 Apr 2006, Thomas Adams wrote:
>
>
>> List:
>>
>> I can not seem to work out the syntax for using R/gstat within a GRASS
>> 6.1 se
age -
From: Roger Bivand <[EMAIL PROTECTED]>
Date: Thursday, April 27, 2006 4:21 pm
Subject: Re: [R-sig-Geo] Example of universal kriging with R/gstat in GRASS
needed
> On Thu, 27 Apr 2006, Thomas Adams wrote:
>
> > List:
> >
> > I can not seem to work out the
List:
I can not seem to work out the syntax for using R/gstat within a GRASS
6.1 session to do universal kriging. I have a DEM (elevation data on a
grid) and point data for temperature; theoretically, the temperatures
should relate to elevation. So, I am trying to spatially interpolate the
tem
er generator
>data(temp_20050801_08): WARNING: Vector used instead of sites.
>GRASS site list temp_20050801_08: 0 cat, 3 dim, 1 str, 4 dbl.
>Warning: gstat/grass: specify z coordinate if you want to use it
>gstat/grass: 117 sites read successfully.
>
>which looks good (to me).
>
>
>
List:
I have point temperature data in GRASS6 (as vectors); I'm using GSTAT
2.4.5 which I have compiled with GRASS support (and confirmed with gstat
-v). I have:
# file: gstat.cmd.temps
#
# gstat command file
#
data(temp_20050101_08): 'temp_20050101_08', x=1, y=2, v=1;
variogram(temp_20050101_0
34 matches
Mail list logo