Re: [GRASS-user] LiDAR Data?

[Markus Neteler]

On Thu, Jun 10, 2010 at 12:14 AM, Rich Shepard  wrote:
>  I thought that I could use LiDAR data within grass,

You can :)

> but the 6.4 r.in.gdal
> man page does not include that format. Is this the wrong module for these
> data?

There is no format "Lidar" - in fact there are many formats related to
Lidar data:

- raw x,z,y,w points in ASCII tables (v.in.ascii)
- LAS format (LAS tools + v.in.ascii)
- other point based ASCII formats (v.in.ascii)
- points in 3D SHAPE format (vector, v.in.ogr)
- rasterized DSM in GeoTIFF or whatever (raster, r.in.gdal)
- ...

>  There are two topographic quadrangle files that cover the project drainage
> basin, and the data DVDs cost $200 each. I want to make high resolution DEMs
> of the project basin and try to determine flood elevations for different
> storm events from those. According to the data source, "All data are format
> specific to ESRI GIS format. Data must be viewed using specialty software
> capable of viewing .shp, geotif, and ESRI grid formats."

So GRASS can do that.

>  I suppose that r.in.gdal will import the .shp files.

No, it is v.in.ogr.

> I've not worked with
> geotif or ESRI grid formats so any experiences of you folks with these
> formats and LiDAR data will be much appreciated.

For GeoTIFF, use r.in.gdal.

To process Lidar point data, then use v.lidar.*.
The already mentioned wiki page may be helpful, too.

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Build error: 6.5 SVN code

[Markus Neteler]

On Thu, Jun 10, 2010 at 1:13 AM, Rich Shepard  wrote:
>  I've been building the latest SVN code for 6.4 and 6.5 here without any
> errors. Today I get this:
>
> Started compilation: Wed Jun  9 15:21:25 PDT 2010
> --
> Errors in:
> /usr4/grass65/gui/wxpython/scripts

...
>  Yet when I follow the directions I cannot run make from that directory:
>
> [rshep...@salmo /usr4/grass65]$ cd gui/wxpython/scripts/
> [rshep...@salmo /usr4/grass65/gui/wxpython/scripts]$ make
> make: *** No targets specified and no makefile found.  Stop.

Did you run "make" with -jX for parallelized compilation? I also get
an error when doing so while plain "make" works. There appears to
be a dependency problem in the Makefile.

please try
cd /usr4/grass65/gui/wxpython/
make

to see if that goes ok.

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR Data?

[Rich Shepard]

On Wed, 9 Jun 2010, Mark Seibel wrote:


This is a great resource.   http://grass.osgeo.org/wiki/LIDAR#Micro-tutorial


Mark,

  Thanks.

Rich
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR Data?

[Mark Seibel]

This is a great resource.   http://grass.osgeo.org/wiki/LIDAR#Micro-tutorial

Mark

On Wed, Jun 9, 2010 at 6:14 PM, Rich Shepard  wrote:
>  I thought that I could use LiDAR data within grass, but the 6.4 r.in.gdal
> man page does not include that format. Is this the wrong module for these
> data?
>
>  There are two topographic quadrangle files that cover the project drainage
> basin, and the data DVDs cost $200 each. I want to make high resolution DEMs
> of the project basin and try to determine flood elevations for different
> storm events from those. According to the data source, "All data are format
> specific to ESRI GIS format. Data must be viewed using specialty software
> capable of viewing .shp, geotif, and ESRI grid formats."
>
>  I suppose that r.in.gdal will import the .shp files. I've not worked with
> geotif or ESRI grid formats so any experiences of you folks with these
> formats and LiDAR data will be much appreciated.
>
> Thanks,
>
> Rich
> ___
> grass-user mailing list
> grass-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user
>
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Build error: 6.5 SVN code

[Nikos Alexandris]

Rich Shepard:
>I've been building the latest SVN code for 6.4 and 6.5 here without any
> errors. Today I get this:
> 
> Started compilation: Wed Jun  9 15:21:25 PDT 2010
> --
> Errors in:
> /usr4/grass65/gui/wxpython/scripts
> --
> In case of errors please change into the directory with error and run
> 'make'.
> If you get multiple errors, you need to deal with them in the order they
> appear in the error log. If you get an error building a library, you will
> also get errors from anything which uses the library.
> --
> Finished compilation: Wed Jun  9 16:06:33 PDT 2010
> make: *** [default] Error 1
> 
>Yet when I follow the directions I cannot run make from that directory:
> 
> [rshep...@salmo /usr4/grass65]$ cd gui/wxpython/scripts/
> [rshep...@salmo /usr4/grass65/gui/wxpython/scripts]$ make
> make: *** No targets specified and no makefile found.  Stop.
 
> Suggestions?

Not specifically to the "failing-to-compile" script. But if you don't need the 
GUI, I guess you can go ahead and compile from the root directory as always. 
GRASS should work and only "miss" these wxpython scripts.

In case you require the GUI... no idea.

Nikos
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR Data?

[Nikos Alexandris]

(This didn't hit the list before...)

On Thursday 10 of June 2010 00:14:29 Rich Shepard wrote:
>I thought that I could use LiDAR data within grass, but the 6.4
> r.in.gdal man page does not include that format. Is this the wrong module
> for these data?
> 
>There are two topographic quadrangle files that cover the project
> drainage basin, and the data DVDs cost $200 each. I want to make high
> resolution DEMs of the project basin and try to determine flood elevations
> for different storm events from those. According to the data source, "All
> data are format specific to ESRI GIS format. Data must be viewed using
> specialty software capable of viewing .shp, geotif, and ESRI grid
> formats."
 
> I suppose that r.in.gdal will import the .shp files. I've not worked
> with geotif or ESRI grid formats so any experiences of you folks with
> these formats and LiDAR data will be much appreciated.

Rich,

Geotiff is more than supported in the OSGeo world. And ESRI grid's I've 
imported (and exported after some trial and error) successfully sometimes, not 
worked frequently with such a format though.

See also [1][2][3]. But the grass-gis-lidar experts will chime in I bet...

My 2 drachmas [4] (kept some of the previous greek coins for any case... :-p), 
Nikos

---
[1] 
[2] 
[3] Relevant: 

(
[4] 
)
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR Data?

[Nikos Alexandris]

On Thursday 10 of June 2010 00:14:29 Rich Shepard wrote:
>I thought that I could use LiDAR data within grass, but the 6.4
> r.in.gdal man page does not include that format. Is this the wrong module
> for these data?
> 
>There are two topographic quadrangle files that cover the project
> drainage basin, and the data DVDs cost $200 each. I want to make high
> resolution DEMs of the project basin and try to determine flood elevations
> for different storm events from those. According to the data source, "All
> data are format specific to ESRI GIS format. Data must be viewed using
> specialty software capable of viewing .shp, geotif, and ESRI grid
> formats."
> 
>I suppose that r.in.gdal will import the .shp files. I've not worked
> with geotif or ESRI grid formats so any experiences of you folks with
> these formats and LiDAR data will be much appreciated.

I am sending you also off-line some work I've done and includes some notes on 
Geotiff's (quite useful for me at the time). Results not include (due to 
external forces :-p). Hope it is worth a little of bit of something to you as 
well...

Cheers, Nikos
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Build error: 6.5 SVN code

[Rich Shepard]

  I've been building the latest SVN code for 6.4 and 6.5 here without any
errors. Today I get this:

Started compilation: Wed Jun  9 15:21:25 PDT 2010
--
Errors in:
/usr4/grass65/gui/wxpython/scripts
--
In case of errors please change into the directory with error and run
'make'.
If you get multiple errors, you need to deal with them in the order they
appear in the error log. If you get an error building a library, you will
also get errors from anything which uses the library.
--
Finished compilation: Wed Jun  9 16:06:33 PDT 2010
make: *** [default] Error 1

  Yet when I follow the directions I cannot run make from that directory:

[rshep...@salmo /usr4/grass65]$ cd gui/wxpython/scripts/
[rshep...@salmo /usr4/grass65/gui/wxpython/scripts]$ make
make: *** No targets specified and no makefile found.  Stop.

Suggestions?

Rich
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] LiDAR Data?

[Rich Shepard]

On Thu, 10 Jun 2010, Nikos Alexandris wrote:


Geotiff is more than supported in the OSGeo world. And ESRI grid's I've
imported (and exported after some trial and error) successfully sometimes,
not worked frequently with such a format though.


Nikos,

  I've imported ESRI grid data as the 10m resolution DEMs for the state
here. But, I've not before worked with geotiff files. I'll read the
references you provided and take all the help I can get grom the GRASS-LiDAR
experts. By next Tuesday I need to have some high resolution
terrain/hydrologic models run to estimate flood elevations given channel
size (from the LiDAR DEM) and point rainfall data over a 2 week period.


My 2 drachmas [4] (kept some of the previous greek coins for any case...
:-p),


  Turn them into gold if the price drops. :-)

Many thanks,

Rich
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] LiDAR Data?

[Rich Shepard]

  I thought that I could use LiDAR data within grass, but the 6.4 r.in.gdal
man page does not include that format. Is this the wrong module for these data?

  There are two topographic quadrangle files that cover the project drainage
basin, and the data DVDs cost $200 each. I want to make high resolution DEMs
of the project basin and try to determine flood elevations for different
storm events from those. According to the data source, "All data are format
specific to ESRI GIS format. Data must be viewed using specialty software
capable of viewing .shp, geotif, and ESRI grid formats."

  I suppose that r.in.gdal will import the .shp files. I've not worked with
geotif or ESRI grid formats so any experiences of you folks with these
formats and LiDAR data will be much appreciated.

Thanks,

Rich
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Import GRIB files to GRASS

[Randy Bradley]

>> I'm planning on using some global meteo datasets that are
>> distributed in GRIB format. has anyone ever tried to parse
>> GRIB NCEP files to GRASS?

It is so interesting to see vastly different approaches to solve a
similar problem.  I'm downloading NOAA forecast data which is in GRIB
format.  I run the files through degrib and then parse the shape files with
a C program using the shapelib libraries and create a new shape file which I
import into GRASS with v.in.ogr.  I plot the values with d.vect.thematic.

The resulting maps can be viewed at:
http://www.ars.usda.gov/npa/marc/heatstress

That reminds me, I originally set up a location for each region with an
eye forward on multithreading the map generation phases.  I forgot the
reason in the middle of the project and created several regions in one
location.  I had to revert back to a location for each region when I
implemented multithreading but reduced the total map generation time from
several hours to 15 minutes.  :-)


Randy


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.neighbors erroneous posting

[Dave Roberts]

Friends,

Sorry about an erroneous posting with respect to running 
r.neighbors with large neighborhoods (51 in this case).  My problem was 
the current region, not r.neighbors which seesm to work fine with large 
windows.


Dave
--

David W. Roberts office 406-994-4548
Professor and Head  FAX 406-994-3190
Department of Ecology email drobe...@montana.edu
Montana State University
Bozeman, MT 59717-3460
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Using Tck/tk vector digitizer for training areas

[Micha Silver]

On 06/09/2010 07:48 PM, Jenny Turner wrote:

Greetings

I'm using WinGRASS (from OSGEO4W installer) to do some supervised 
classification. Since I'm using winGRASS i.class function is not 
available neither wxGUI's vector digitizer. So I'm using tcl/tk but 
I'm a bit confused on using it.
About the background images I used a discussion from last week 
regarding this but now.. How can I define a category? I mean, How can 
atribute a polygon to a certain category?
This indeed is not obvious. What I've done is change the category mode 
from "Next not used" to "Manual entry". This is set in the lower part of 
the v.digit window (in tcltk interface). Now each time you digitize a 
*centroid* for a training area, set its cat value to some standard value 
you choose. For example, you might decide:

train_areacat

forest 1
agriculture2
urban  9
water 10


Now for each training polygon which is a forest, set the cat value of 
its centroid to 1, agri areas will get value 2, and so on for all the 
centroids of all the training polygons. (Note that the cat values for 
the boundary lines are pretty much irrelevant in GRASS's vector model).


When you're done digitizing training areas, you might add another column 
or two to the attrib table for the training areas -

v.db.addcol train_areas col="train_label varchar(16), RGBCOLOR varchar(16)"
then give each centroid cat value a suitable label and an RGB color 
combo, etc.  Makes life easier down the line...




This mail was received via Mail-SeCure System.


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

This mail was received via Mail-SeCure System.





--
Micha Silver
Arava Development Co. +972-52-3665918
http://www.surfaces.co.il


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.neighbors

[Dave Roberts]

It appears that running r.neighbors with really large neighborhoods
results in mostly NULLs.  A neighborhood of 31 seemed to work, but 51 
returned only a tiny portion of the map.


Dave
--

David W. Roberts office 406-994-4548
Professor and Head  FAX 406-994-3190
Department of Ecology email drobe...@montana.edu
Montana State University
Bozeman, MT 59717-3460
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] r.neighbors with large neighborhoods

[Dave Roberts]

I'm working landform analysis using r.neighbors.  The help file says 
that "size" can go up to 25 (or +/1 12 pixels).  I need much bigger 
neighborhoods than that, so I just entered a large number and it seems 
to work (although it's slow).  Does anybody know what the true limits of 
the r.neighbors code really is?


Thanks, Dave
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Import GRIB files to GRASS

[António Rocha]

Hello Hamish.

Yest, I used a NOAA's tool to verify it and it's misplaced. I will take 
a look at those tickets that you mention.



About this: /There is a cell-registration bug in GDAL: (grid vs 
cell-center convention) cells are all offset to the east and south by 
half a cell. You can fix this by hand with the r.region 
<http://grass.osgeo.org/grass64/manuals/html64_user/r.region.html> module.

/How can I use r.region to fix this?


Hamish wrote:

António wrote:
  

I have one question regarding Hamish's link.
/There is a cell-registration bug in GDAL: (grid vs
cell-center convention) cells are all offset to the east and
south by half a cell. You can fix this by hand with the
r.region <http://grass.osgeo.org/grass64/manuals/html64_user/r.region.html>
module.

/Is this still happening, this error?/




well, the linked GDAL bug tickets are still open:
  https://trac.osgeo.org/gdal/ticket/2637
  https://trac.osgeo.org/gdal/ticket/2550

so I guess so. I'd encourage you to verify for yourself if your data is
offset or not, and post that verification (or not) to the appropriate
bug report with as much detail as you can; and if the patch attached
to the ticket fixes it for you (or not).


Hamish



  



__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5181 (20100608) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com




  




__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5185 (20100609) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Using Tck/tk vector digitizer for training areas

[Jenny Turner]

Greetings

I'm using WinGRASS (from OSGEO4W installer) to do some supervised
classification. Since I'm using winGRASS i.class function is not available
neither wxGUI's vector digitizer. So I'm using tcl/tk but I'm a bit confused
on using it.
About the background images I used a discussion from last week regarding
this but now.. How can I define a category? I mean, How can atribute a
polygon to a certain category?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Changing menudata.xml to include an direct processing

[Franz Schiller]

Greetings

I'm making some addaptions in menudata.xml in order to have a few functions
available more rapidly (for me). I wanted to have a menu-item (I know how to
add it) that runs directly without opening any panel. This would mean that
all parameterization would be in the xml. Is this possible?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Temporal Analysis

[Markus Neteler]

On Wed, Jun 9, 2010 at 6:30 PM, Nikos Dumakis  wrote:
> Now I'm a bit lost on this versions stuff :)
>>
>> > - Do you have that specific code available?
>>
>> See SVN :)
>>
> What is the name of the module? r.series?

Yes.

>> > About: "Find it in GRASS 6.5 and 7 (also 6.4.1 once 6.4.0 is out):
>> > http://grass.osgeo.org/grass70/manuals/html70_user/r.series.html ";
>> > - You are saying that r.series function at GRASS6.5 and 7 includes that
>> > possibility?
>>
>> Yes.
>>
>> > And can it be back-ported to 6.4?
>>
>> Yes, but only for 6.4.1.
>
> What is the 6.4.1 release? I only heard about 6.4.0RC6, GRASS6.5 and 7 never
> 6.4.1... Is it as stable as GRASS6.4.0RC6?

- 6.4.0RC6 is the 6th release candidate of 6.4.0.
- 6.4.1 will be the successor of 6.4.0 (likely there will be release
candidates, too).
- etc.

>> > And about "6.4.0 out": Isn't Out already? :D
>>
> Isn't it 6.4.0RC6 the GRASS6.4.0 ?

No, it is just a release candidate. We found bugs therein which have
been subsequently
fixed (available via SVN, will be in 6.4.0[final]). Still some blocker
bugs need to be
fixed before releasing 6.4.0.

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Temporal Analysis

[Nikos Dumakis]

Now I'm a bit lost on this versions stuff :)

>
> > - Do you have that specific code available?
>
> See SVN :)
>
> What is the name of the module? r.series?



> > About: "Find it in GRASS 6.5 and 7 (also 6.4.1 once 6.4.0 is out):
> > http://grass.osgeo.org/grass70/manuals/html70_user/r.series.html "
> > - You are saying that r.series function at GRASS6.5 and 7 includes that
> > possibility?
>
> Yes.
>
> > And can it be back-ported to 6.4?
>
> Yes, but only for 6.4.1.
>
What is the 6.4.1 release? I only heard about 6.4.0RC6, GRASS6.5 and 7 never
6.4.1... Is it as stable as GRASS6.4.0RC6?



>
> > And about "6.4.0 out": Isn't Out already? :D
>
> Isn't it 6.4.0RC6 the GRASS6.4.0 ?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Martin Landa]

[back to ML]

Hi,

2010/6/9 Rebecca Bennett :
> Ubuntu Hardy, I was entering the commands in the Output section of the GUI,
> but  have now entered it into the bash terminal and it worked. I did not
> realise that there were commands that were fussy about where they were
> input... thanks for your patience.

so it's a TCL/TK GUI bug. Please fill the ticket in the GRASS trac.

Also please consider usage of wxGUI [1] which is also available in
GRASS 6.4. More advanced wxGUI can be found in GRASS 6.5.

g.gui wxpython

Martin

[1] http://grass.osgeo.org/wiki/WxGUI

-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Martin Landa]

Hi,

2010/6/9 Rebecca Bennett :
> I have deleted the column I added and tried all command line permutations of
> v.db.addcol map=contour_points_c...@permanent layer=1 `col="rast_value
> DOUBLE PRECISION"`
> v.db.addcol map=contour_points_c...@permanent layer=1 `columns="rast_value
> DOUBLE PRECISION"`

ones more - the correct form is (without ``)

v.db.addcol map=contour_points_c...@permanent layer=1 col="rast_value
DOUBLE PRECISION"

Note that you can edit only vector maps from the current mapset.

Martin

-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Previewing column names in input ASCII file without messing g.parser up

[Luigi Ponti]

Hi there,

Working on windows 7, I was trying to run a GRASS bash shell script 
(written for the Cygwin version of GRASS) using the 
WinGRASS-6.4.SVN-r42329-1 installed from 
. After some tuning to adopt the 
direct rendering approach and to adjust to msys pathname conventions, 
all runs fine (winGRASS is in "Program Files" and no problems occur 
using a variety of grass and shell commands) except for one feature of 
the script that I am not able to reproduce in the native winGRASS version.


The shell script basically imports (v.in.ascii) an ASCII text file with 
columns including latitude, longitude, and some variables, and then maps 
a variable of choice. It calls a little Perl script that prints column 
order and names before g.parser is called. In the Cygwin version, this 
results in (1) the output of the Perl script to be printed to the 
GRASS/bash terminal, and (2) the script GUI to be launched as well. This 
way, the user can look at names/numbers of columns in the ASCII file, 
and input column numbers of choice in the GUI for v.in.ascii.


In winGRASS, the execution of the Perl script that prints column 
names/numbers prevents the GUI of the main shell script from starting up 
when called in the Cmd> text box of the Python wxgui (started using the 
shortcut created by the installer), with the following error message in 
the command output tab:


   self.goutput.RunCmd(cmd, switchPage=False)
 File "C:\Program
   Files\GRASS-64-SVN\etc\wxpython\gui_modules\goutput.py",
   line 354, in RunCmd

   menuform.GUI().ParseCommand(cmdlist, parentframe=self)
 File "C:\Program
   Files\GRASS-64-SVN\etc\wxpython\gui_modules\menuform.py",
   line 1853, in ParseCommand

   handler)
 File "C:\OSGeo4W\apps\Python25\lib\xml\sax\__init__.py",
   line 49, in parseString

   parser.parse(inpsrc)
 File
   "C:\OSGeo4W\apps\Python25\lib\xml\sax\expatreader.py", line
   107, in parse

   xmlreader.IncrementalParser.parse(self, source)
 File "C:\OSGeo4W\apps\Python25\lib\xml\sax\xmlreader.py",
   line 123, in parse

   self.feed(buffer)
 File
   "C:\OSGeo4W\apps\Python25\lib\xml\sax\expatreader.py", line
   211, in feed

   self._err_handler.fatalError(exc)
 File "C:\OSGeo4W\apps\Python25\lib\xml\sax\handler.py",
   line 38, in fatalError

   raise exception
   xml.sax._exceptions
   .
   SAXParseException
   :
   :2:0: syntax error

However, when calling the shell script in winGRASS text mode (started 
using the shortcut created by the installer) the Cygwin behavior is 
reproduced, i.e. column numbers and names are printed to the DOS 
terminal by the Perl script, and the GUI of the main shell script is 
launched as well. Apparently, this is due to the wxpython not having a 
DOS prompt where to print. Also, it seems to be calling code from 
OSGeo4W although I have checked my %PATH% environment variable, and 
there does not seem to be any reference to the C:\OSGeo4W installation: 
typing


   path | FINDSTR OSGeo4W

at the DOS prompt returns nothing.

Does anybody see a way around the wxpython problem that would allow 
column number/names in the input ASCII file to be previewed without 
preventing the g.parser GUI to start up?


A side note: it seems like the /scripts/ folder is not accessible inside 
custom shell scripts so that I had to add


   @set PATH=%GISBASE%\scripts\;%PATH%

in the corresponding .bat file that lives inside GISBASE/bin/ and is 
used to call the custom script itself from wxpython gui. Not sure why 
this happens.


Finally, I wanted to thank authors of posts from which I have collected 
info that enabled me to port shell scripts to native winGRASS (e.g. 
Glynn Clements and Elena Mitasova).


Kind regards and sorry for a long post of somewhat mixed content,

Luigi




___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Rebecca Bennett]

Hi all,

I have deleted the column I added and tried all command line permutations of
v.db.addcol map=contour_points_c...@permanent layer=1 `col="rast_value DOUBLE 
PRECISION"`
v.db.addcol map=contour_points_c...@permanent layer=1 `columns="rast_value 
DOUBLE PRECISION"`
that have so far been given in this thread, again with quotes and without. 
Unfortunately it does not recognise this command, which I why I used the GUI

The GUI is tcltk on GRASS 6.4

Micha - I have checked the points are within the raster (they were created from 
a contour map of the raster). The region is consistent with both the raster and 
vector

Thanks!





From: Micha Silver 
To: Rebecca Bennett 
Cc: Martin Landa ; grass-user@lists.osgeo.org
Sent: Wed, 9 June, 2010 10:12:28
Subject: Re: [GRASS-user] Re: Problems with v.what.rast

 On 09/06/2010 11:40, Rebecca Bennett wrote: 
> 
>Hi Martin,
>
>>I tried the command that Micha suggested v.db.addcol
>map=contour_points_c...@permanent layer=1
>>columns="rast_value double precision" but this age the following error
>
>
The quotes " " around the "rast_value DOUBLE PRECISION" are necessary.
Just to check can we see the output of 
v.info -c Contour_points_cats

You might also check that the points are all within the raster region.
i.e.:
v.info -g Contour_points_cats
r.info -g Dartmoor1
g.region -p
 

v.db.addcol Contour_points_cats col=rast_value DOUBLE PRECISION
>>Sorry  is not a valid option
>>Sorry  is not a valid option
>
>>So through the GUI I added a column to the attribute table which gave
>the command line out put 
>>v.db.addcol map=contour_points_c...@permanent layer=1
>{columns=rast_value double precision}
>
>>cheers!
>
>
>
>

From: >Martin Landa 
>To: Rebecca Bennett
>
>Cc: Micha Silver
>; grass-user@lists.osgeo.org
>Sent: Wed, 9 June,
>2010 9:37:14
>Subject: Re:
>[GRASS-user] Re: Problems with v.what.rast
>
>>Hi,
>
>>2010/6/9 Rebecca Bennett :
>>> I added a column v.db.addcol map=contour_points_c...@permanent
>layer=1
>>> {columns=rast_value double precision} (as my raster is a float),
>then ran
>>> with error
>
>>how did you get this command?
>
>>It should be
>
>>v.db.addcol map=contour_points_c...@permanent layer=1
>>columns="rast_value double precision"
>
>>Martin
>
>>-- 
>>Martin Landa  * http://gama.fsv.cvut.cz/~landa
>
>
>
>>This mail was received via Mail-SeCure System.
>


-- 
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918



  ___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Micha Silver]

On 09/06/2010 11:40, Rebecca Bennett wrote:

  
  
  Hi Martin,
  
I tried the command that Micha suggested v.db.addcol
map=contour_points_c...@permanent layer=1
columns="rast_value double precision" but this age the following error
  
  
  

The quotes " " around the "rast_value DOUBLE PRECISION" are necessary.
Just to check can we see the output of 
v.info -c Contour_points_cats

You might also check that the points are all within the raster region.
i.e.:
v.info -g Contour_points_cats
r.info -g Dartmoor1
g.region -p
 

  
  v.db.addcol Contour_points_cats col=rast_value DOUBLE PRECISION
Sorry  is not a valid option
Sorry  is not a valid option
  
So through the GUI I added a column to the attribute table which gave
the command line out put 
v.db.addcol map=contour_points_c...@permanent layer=1
{columns=rast_value double precision}
  
cheers!
  
  
  
  From:
Martin Landa 
  To: Rebecca Bennett

  Cc: Micha Silver
; grass-user@lists.osgeo.org
  Sent: Wed, 9 June,
2010 9:37:14
  Subject: Re:
[GRASS-user] Re: Problems with v.what.rast
  
Hi,
  
2010/6/9 Rebecca Bennett :
> I added a column v.db.addcol map=contour_points_c...@permanent
layer=1
> {columns=rast_value double precision} (as my raster is a float),
then ran
> with error
  
how did you get this command?
  
It should be
  
v.db.addcol map=contour_points_c...@permanent layer=1
columns="rast_value double precision"
  
Martin
  
-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa
  
  
  
  
  
  
This mail was received via Mail-SeCure System.



-- 
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Martin Landa]

Hi,

2010/6/9 Rebecca Bennett :
> I tried the command that Micha suggested v.db.addcol
> map=contour_points_c...@permanent layer=1
> columns="rast_value double precision" but this age the following error
>
> v.db.addcol Contour_points_cats col=rast_value DOUBLE PRECISION
> Sorry  is not a valid option
> Sorry  is not a valid option

the quotes are required, so `col="rast_value DOUBLE PRECISION"`.

> So through the GUI I added a column to the attribute table which gave the
> command line out put
> v.db.addcol map=contour_points_c...@permanent layer=1 {columns=rast_value
> double precision}

Which GUI (TCL/TK or wxGUI), which version of GRASS?

Martin

-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Rebecca Bennett]

Hi Martin,

I tried the commandthat Micha suggested v.db.addcol 
map=contour_points_c...@permanent layer=1
columns="rast_value double precision" but this age the following error

v.db.addcol Contour_points_cats col=rast_value DOUBLE PRECISION
Sorry  is not a valid option
Sorry  is not a valid option

So through the GUI I added a column to the attribute table which gave the 
command line out put 
v.db.addcol map=contour_points_c...@permanent layer=1 {columns=rast_value 
double precision}

cheers!





From: Martin Landa 
To: Rebecca Bennett 
Cc: Micha Silver ; grass-user@lists.osgeo.org
Sent: Wed, 9 June, 2010 9:37:14
Subject: Re: [GRASS-user] Re: Problems with v.what.rast

Hi,

2010/6/9 Rebecca Bennett :
> I added a column v.db.addcol map=contour_points_c...@permanent layer=1
> {columns=rast_value double precision} (as my raster is a float), then ran
> with error

how did you get this command?

It should be

v.db.addcol map=contour_points_c...@permanent layer=1
columns="rast_value double precision"

Martin

-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa



  ___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Re: Problems with v.what.rast

[Martin Landa]

Hi,

2010/6/9 Rebecca Bennett :
> I added a column v.db.addcol map=contour_points_c...@permanent layer=1
> {columns=rast_value double precision} (as my raster is a float), then ran
> with error

how did you get this command?

It should be

v.db.addcol map=contour_points_c...@permanent layer=1
columns="rast_value double precision"

Martin

-- 
Martin Landa  * http://gama.fsv.cvut.cz/~landa
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Re: Problems with v.what.rast

[Rebecca Bennett]

Hi Micha 

Thanks for your help. I followed your instructions to add add separate cat 
values for each point, so the results of v.category Contour_points_cats 
opt=report now looks like this 
v.category Contour_points_cats opt=report
Layer: 2
type   countminmax
point 644858  1 644858
line   0  0  0
boundary   0  0  0
centroid   0  0  0
area   0  0  0
all   644858  1 644858
Layer/table: 1/Contour_points_cats
type   countminmax
point 644858  1 644858
line   0  0  0
boundary   0  0  0
centroid   0  0  0
area   0  0  0
all   644858  1 644858

I added a column v.db.addcol map=contour_points_c...@permanent layer=1 
{columns=rast_value double precision} (as my raster is a float), then ran with 
error

v.what.rast vector=contour_points_c...@permanent raster=dartmo...@permanent 
layer=1 column=rast_value

No record for category 553 in table 

No record for category 263 in table 
No record for category 668 in table 

essentially the "No Record" error was repeated lots of times and the tool 
stopped working.


I'm a bit stumped by this - any suggestions? (I have double checked my region 
to make sure it incorporates the whole area of interest)

cheers,
Rebecca




From: Micha Silver 
To: Rebecca Bennett 
Cc: grass-user@lists.osgeo.org
Sent: Tue, 8 June, 2010 17:57:57
Subject: Re: [GRASS-user] Problems with v.what.rast

On 06/08/2010 05:39 PM, Rebecca Bennett wrote: 
> 
>Dear all,
>
>>I'm having some trouble copying the values from a raster image
>(Dartmoor1) to a series of vector points (Contour_points - created from
>a contour vector map). When I try to use v.what.rast I get the
>following error:
>
>>v.what.rast vector=contour_poi...@permanent raster=dartmo...@permanent
>layer=1 column=level
>>More points (644858) of category 1, value set to 'NULL'
>>1 categories loaded from table
>>1 categories loaded from vector
>>0 categories from vector missing in table
>>1 duplicate categories in vector
>>1 records updated
>>0 update errors
>
>
Looks like all the points in Contour_points have the same cat value.
Check this with:
v.category Contour_points opt=report

You can add new separate cat values for each point like this:
v.category Contour_points opt=del out=Contour_points_tmp
v.category Contour_points_tmp opt=add out=Contour_points_cats 
(get rid of the tmp vector: g.remove vect=Contour_points_tmp)

Now check:
v.category Contour_points_cats opt=report

Next setup a column for the rast values:
v.db.addcol Contour_points_cats col="rast_value double" 
(Is your raster an integer or float raster??)
And now you should be able to get the raster values into the points
attrib table with:

v.what.rast Contour_points_cats rast=Dartmoor1 col=rast_value 



I think that this is an attribute table issue. I have
>followed the instructions in Example B 
>http://grass.itc.it/grass62/manuals/html62_user/v.what.rast.html >but  get the 
>following message 
>
>>v.to.db map=Contour_points option=cat
>>Reading features...
>>Updating database...
>>1 categories read from vector map (layer 1)
>>1 records selected from table (layer 1)
>>1 categories read from vector map exist in selection from table
>>0 records updated/inserted (layer 1)
>>Record (cat 1) already exists (not inserted)
>
>>I have also tried creating and attaching a new table but get the same
>result for v.what.vect. Any ideas welcome, thanks for reading!
>
>>Rebecca
>
>
>
>>This mail was received via Mail-SeCure System.
>
>
>___
>grass-user mailing list
>grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
>
>This mail was received via Mail-SeCure System.
>
>
>  


-- 
Micha Silver
Arava Development Co. +972-52-3665918
http://surfaces.co.il 



  ___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user