[GRASSLIST:8104] Re: r.profile & subsequent chart plot

Ian MacMillan ian_macmillan at umail.ucsb.edu
Sun Aug 28 17:46:16 EDT 2005


hi all, sounds silly, but I just make a graph with excel, and then you 
can save that graph as a postscript or whatever other format is useful.

-Ian

On Aug 28, 2005, at 2:37 AM, Miha Staut wrote:

> Good plotting programmes are R or Gnuplot.
> Example with the use of R:
>
> GRASS 6.1.cvs (loc):~ > r.profile -i in=yourrast null="NA" > profile
>> R   # start R
>> profile <- read.table("profile", sep=" ", head=F)
>> str(pr)
> `data.frame':   439 obs. of  2 variables:
>  $ V1: num  0 100 200 300 400 500 600 700 800 900 ...
>  $ V2: int  NA NA NA NA NA NA NA NA NA NA ...
>> plot(pr$V1,pr$V2,ty="l")
>
> # If you wish to get height vs. length in equal scale use argument 
> asp=1, e.g.
>> plot(V1, V2, xlim = c(min(V1, na.rm=T), max(V1, na.rm=T)), ylim = 
>> c(min(V2,
> na.rm=T), max(V2, na.rm=T)), asp = 1, ty = "l")
>
> Miha Staut
>
>
> --- M S <mseibel at gmail.com> wrote:
>
>> i managed to get a nice output from grass's r.profile to get a cross 
>> section.
>>
>> anyone recommend a good program to then plot it with to creat the
>> cross section graphic?
>
>
>
> 	
> 	
> 		
> ___________________________________________________________
> Yahoo! Messenger - NEW crystal clear PC to PC calling worldwide with 
> voicemail http://uk.messenger.yahoo.com
>




More information about the grass-user mailing list