[R-sig-Geo] DEM in R

[Tomislav Hengl]

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