DEM in R

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20090722/3ae292cf/attachment.pl>
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20090722/fb6ab52e/attachment.pl>
Dear Craig,

Plotting of DEMs and geomorphometric analysis is possible in R. For example, to derive slope/aspect,
you could use this script:

https://www.msu.edu/~ashton/classes/866/notes/lab12/slope.R
(https://stat.ethz.ch/pipermail/r-sig-geo/2009-July/006088.html)

Here is an example of a perspective plot:
http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=27

But you will soon find out that R can not really work well with large maps (and you do not want to
re-implement all algorithms for DEM analysis for scratch), so you should also consider using some
open source GIS (that can be used in combination with R):

"Analysis of DEMs in R+ILWIS/SAGA"
http://spatial-analyst.net/wiki/index.php?title=Analysis_of_DEMs_in_R%2BILWIS/SAGA 

"Geomorphometry in R+SAGA+ILWIS+GRASS"
http://geomorphometry.org/R.asp 

Welcome to R-sig-geo!

T. Hengl

PS: Your e-mail footer says that you are sending "confidential information", that we should "not
disseminate, copy or take any action in reliance on it" and that if this is true we should "notify
the Department of Agriculture and Food, WA immediately". Funny.
-----Original Message-----
From: r-sig-geo-bounces at stat.math.ethz.ch [mailto:r-sig-geo-bounces at stat.math.ethz.ch] On Behalf
Of milton ruser
Sent: Wednesday, July 22, 2009 6:57 AM
To: Scanlan, Craig
Cc: r-sig-geo at stat.math.ethz.ch
Subject: Re: [R-sig-Geo] DEM in R

Hi Craig,

I think you can use several ways of generate your DEM from xyz points.
One is you use interp function on akima package.

For derive slope and aspect, may be you can use R with GRASS (package
spgrass6) or R with SAGA (package RSAGA).

How about you give a look at:
http://spatial-analyst.net/wiki/index.php?title=Main_Page or
http://www.geomorphometry.org/

bests

milton
brazil=toronto

On Tue, Jul 21, 2009 at 8:27 PM, Scanlan, Craig <
craig.scanlan at agric.wa.gov.au> wrote:

Hi,

I need to create a matrix of elevation data to plot with persp(). I have
data in the form of vectors x,y and z. Also, I need to be able to derive
attributes like slope and aspect from the DEM. I have not been able to
find any examples of these problems. Are there any available? Any
suggestions?

Thanks

Craig.

This e-mail and files transmitted with it are privileged and confidential
information
intended for the use of the addressee. The confidentiality and/or privilege
in this e-mail is
not waived, lost or destroyed if it has been transmitted to you in error.
If you received this
e-mail in error you must
(a) not disseminate, copy or take any action in reliance on it;
(b) please notify the Department of Agriculture and Food, WA immediately by
return e-mail to the sender;
(c) please delete the original e-mail.

This email has been successfully scanned by
McAfee Anti-Virus software.
Department of Agriculture and Food WA

       [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

	[[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo