[GRASS-dev] d.where in gis.m ?

[Moritz Lennert]

Michael Barton wrote:
> Moritz,
> 
> Display your vector, but select a raster layer in the gism layer tree. Then
> you'll get a much simpler output (because rasters only have e,n,cat,label)
> on a single line. Save the output window to a file and this might do what
> you want.


This is already a bit better, but besides the fact that it is a bit of a 
hassle to create a raster just to get nicer coordinate output, there is 
still quite a lot of clutter with the command lines:

r.what -f input=pop_poten east_north=151323.550767,160949.243195


151323.550767|160949.243195||36458.9120110746|


r.what -f input=pop_poten east_north=151323.550767,160949.243195


151323.550767|160949.243195||36458.9120110746|


I just think that it should'nt be too complicated to copy the 
information that already exists (as it is displayed on the bottom right 
of the screen) to the output window in a nice format. d.path does this, 
for example:

L: from  M: to R: quit
Node 392: 147162.820000 168433.076000
Node 6812: 152070.313000 168202.302000
Node 2062: 148133.692000 167472.842000
Node 6864: 150575.348000 169256.368000
Node 1785: 153084.305000 168910.834000
Node 7870: 148412.785000 165596.435000
Node 204: 149007.947000 168352.584000
Node 5011: 153025.424000 166352.285000
Node 4364: 147451.703000 168614.870000
Node 4773: 151505.647000 168907.431000

But d.path needs an X monitor...

I know the best answer would be to do it myself, but I am not a big 
expert in tcl/tk.

Moritz