[GRASS-user] Problem with v.digit -- "region for current mapset is not set"

Mon Jan 14 16:58:14 PST 2013

On Mon, Jan 14, 2013 at 02:24:48PM -0700, we recorded a bogon-computron collision of the <russo at bogodyn.org> flavor, containing:
> I'm attempting to do something I haven't done since March of 2011 --- digitizing
> a trail off of a scanned/georeferenced paper map.  I'm using Grass 6.5 (the
> svn version of develbranch_6).  I'm having trouble with v.digit even though
> I'm doing *precisely* the same thing I have done over and over again for years,
> just not recently.  I believe that a change made to v.digit between the last 
> time I did this and today may be the cause, and I don't know what to do to get
> it working again.

Confirmed: I went ahead and undid the changes to v.digit between revisions 
49662 and 49668, and this problem went away (and of course, v.digit now once 
again does what that commit was meant to stop it doing, monkeying with the 
current region, but I'm OK with that).

At least I'm able to do what I needed to do tonight.  But I'm curious whether
anyone else has seen v.digit start to complain about region setting since
December 2011 in Grass 6.x SVN.

> The problem I'm having is that when I fire up 
> 'v.digit -n map=newmapname bgcmd="d.rast basemapname"' I *always* get an
> error message that "ERROR: region for current mapset is not set" followed by
> instruction to run g.region.  At this point, I'm stuck, because even though
> the tool bar and digitizing window show up, the base raster that I'm trying to
> digitize off of is not displayed.
> 
> Thing is, the region IS set --- right before running v.digit I have done
> a "g.region rast=basemapname".
> 
> Looking through SVN logs, it seems that a minor mod to v.digit was made back
> in December of 2011, some months after the last time I did this process (which
> was in March of 2011, right before a backcountry ski trip when I was 
> digitizing a published trail map so I could make a custom map for my GPS).  
> Since that change was precisely one that modified how v.digit handles its
> temporary region, I'm wondering if it's busted.
> 
> Here's what I find:
> 
> v.digit *has* created a tmp_vdigit_<pid> file in the "windows" directory
> of my database for the location/mapset.  The contents of this file precisely
> match the WIND file for the mapset.  Yet somehow, v.digit still thinks
> the region isn't set.  I am not sure if it's getting WIND_OVERRIDE set properly,
> which it is supposedly doing via a putenv call.  
> 
> I tried turning on debugging with g.gisenv, but I didn't see anything in
> the output that meant anything to me.  This is what I see:
> 
> D2/3: Variable = 0x805b560
> D2/3: cmd 0 : d.rast estincomp
> D2/3: bg_add(): cmd = d.rast estincomp
> D2/3: cmd 0 : d.rast estincomp
> D2/3: G__read_Cell_head
> D2/3: G__read_Cell_head_array
> D3/3: region item: proj:       1
> D3/3: region item: zone:       13
> D3/3: region item: north:      4374022
> D3/3: region item: south:      4366938
> D3/3: region item: east:       359460
> D3/3: region item: west:       354486
> D3/3: region item: cols:       2487
> D3/3: region item: rows:       3542
> D3/3: region item: e-w resol:  2
> D3/3: region item: n-s resol:  2
> D3/3: region item: top:        1.000000000000000
> D3/3: region item: bottom:     0.000000000000000
> D3/3: region item: cols3:      2487
> D3/3: region item: rows3:      3542
> D3/3: region item: depths:     1
> D3/3: region item: e-w resol3: 2
> D3/3: region item: n-s resol3: 2
> D3/3: region item: t-b resol:  1
> D1/3: Region: N = 4374022.000000 S = 4366938.000000 E = 359460.000000 W = 354486.000000
> D2/3: Vect_open_new(): name = tenth_mtn_trails4
> D1/3: Vect_set_thresh(): thresh = 0.000000
> D1/3: V1_open_new_nat(): name = tenth_mtn_trails4
> D1/3: Vect_get_thresh(): thresh = 0.000000
> D1/3: dig__write_head()
> D1/3: write coor size (0) to head
> D2/3: coor body offset 14
> D3/3: dig_init_plus()
> D1/3: dig_spidx_init()
> D3/3: dig_cidx_init()
> D3/3: Vect_build(): build = 4
> Building topology for vector map <tenth_mtn_trails4>...
> D2/3: dig_cidx_free()
> D3/3: dig_cidx_init()
> D3/3: Vect_build_nat() build = 4
> D1/3: Vect_Rewind(): name = tenth_mtn_trails4
> Registering primitives...
> D3/3: Vect_read_next_line()
> D3/3: V1_read_next_line_nat()
> D3/3: Vect__Read_line_nat: offset = 14
> 0 primitives registered
> 0 vertices registered
> Building areas...
> 0 areas built
> 0 isles built
> Attaching islands...
> Attaching centroids...
> D2/3: dig_cidx_sort()
> Number of nodes: 0
> Number of primitives: 0
> Number of points: 0
> Number of lines: 0
> Number of boundaries: 0
> Number of centroids: 0
> Number of areas: 0
> Number of isles: 0
> D1/3: Vect_close(): name = tenth_mtn_trails4, mapset = PERMANENT, format = 0, level = 2
> D1/3: get coor info: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/coor
> D1/3: Info->size = 14, Info->mtime = 1358198180
> D1/3: Vect_save_topo()
> D1/3: Open topo: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/topo
> D2/3: topo body offset 142
> D2/3: topo body offset 142
> D2/3: Vect_cidx_save()
> D2/3: Open cidx: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/cidx
> D3/3: dig_write_cidx_head()
> D3/3: cidx body offset 13
> D3/3: dig_write_cidx_head()
> D3/3: cidx body offset 13
> D1/3: close history file
> D1/3: V1_close_nat(): name = tenth_mtn_trails4 mapset= PERMANENT
> D1/3: get coor info: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/coor
> D1/3: Info->size = 14, Info->mtime = 1358198180
> D1/3: dig__write_head()
> D1/3: write coor size (14) to head
> D2/3: coor body offset 14
> D1/3: Vect_get_thresh(): thresh = 0.000000
> D1/3: Vect_write_dblinks(): map = tenth_mtn_trails4, mapset = PERMANENT
> D1/3: dbln file: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/dbln
> D1/3: Dblinks written
> New empty map created.
> D1/3: Vect_open_old(): name = tenth_mtn_trails4 mapset= PERMANENT update = 1
> D1/3: Vect_set_thresh(): thresh = 0.000000
> D3/3: dig_init_plus()
> D1/3: dig_spidx_init()
> D3/3: dig_cidx_init()
> D1/3: open format file: 'PERMANENT/vector/tenth_mtn_trails4/frmt'
> D1/3: Vector format: 0 (native)
> D1/3: Vect_set_thresh(): thresh = 0.000000
> D1/3: Vect__read_head(): vector = tenth_mtn_trails4 at PERMANENT
> D1/3: Vect_set_thresh(): thresh = 0.000000
> D1/3: Level request = 0
> D1/3: Vect_open_topo(): name = tenth_mtn_trails4 mapset= PERMANENT
> D1/3: get coor info: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/coor
> D1/3: Info->size = 14, Info->mtime = 1358198180
> D2/3: Topo header: file version 5.0 , supported from GRASS version 5.0
> D2/3:   byte order 0
> D2/3:   header size 142
> D2/3:   with_z 0
> D2/3:   coor size 14
> D1/3: Topo head: coor size = 14, coor mtime = 1358198180
> D1/3: dig_load_plus()
> D3/3: dig_init_plus()
> D1/3: dig_spidx_init()
> D3/3: dig_cidx_init()
> D2/3: Topo header: file version 5.0 , supported from GRASS version 5.0
> D2/3:   byte order 0
> D2/3:   header size 142
> D2/3:   with_z 0
> D2/3:   coor size 14
> D3/3: dig_alloc_isle():
> D2/3: Vect_cidx_open(): name = tenth_mtn_trails4 mapset= PERMANENT
> D3/3: dig_cidx_init()
> D3/3: dig_read_cidx()
> D3/3: dig_cidx_init()
> D3/3: Cidx header: file version 5.0 , supported from GRASS version 5.0
> D3/3:   byte order 0
> D3/3:   header size 9
> D1/3: V1_open_old_nat(): name = tenth_mtn_trails4 mapset = PERMANENT
> D2/3: Coor header: file version 5.1 , supported from GRASS version 5.1
> D2/3:   byte order 0
> D2/3:   header size 14
> D2/3:   with_z 0
> D2/3:   coor size 14
> D1/3: get coor info: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/coor
> D1/3: Info->size = 14, Info->mtime = 1358198180
> D1/3: coor size in head = 14, real coor file size= 14
> D1/3: Vect_open_old(): vector opened on level 2
> D1/3: Vect_read_dblinks(): map = tenth_mtn_trails4, mapset = PERMANENT
> D3/3: Searching for FID column in OGR DB
> D1/3: dbln file: /users/russo/GIS/Colorado/PERMANENT/vector/tenth_mtn_trails4/dbln
> D1/3: Dblinks read
> D1/3: Vect_Rewind(): name = tenth_mtn_trails4
> D3/3: Vect_build_sidx_from_topo()
> D1/3: dig_spidx_init()
> D3/3: Spatial index was built
> D3/3: Vect_hist_command()
> D3/3: G_recreate_command()
> D1/3: Map opened
> D2/3: cat_init()
> D2/3: symb_line_init()
> D2/3: symb_node_init()
> D3/3: Starting Tk_Main.
> D3/3: v.digit Tcl_AppInit (...)
> D3/3: Starting toolbox.tcl
> D3/3: c_tool_centre()
> D2/3: var_init
> D2/3: Variable = 0x805b560
> D2/3: cat_max_get() field = 1
> D2/3: display_bg()
> ERROR: region for current mapset is not set
>        run "g.region"
> D2/3: display_map()
> 
> 
> I am unable to proceed here.  What else should I be looking at?  Has anyone
> else been able to use v.digit in the 6.x development branch?
> 
> -- 
> Tom Russo    KM5VY   SAR502   DM64ux          http://www.swcp.com/~russo/
> Tijeras, NM  QRPL#1592 K2#398  SOC#236        http://kevan.org/brain.cgi?DDTNM
>  echo "prpv_a'rfg_cnf_har_cvcr" | sed -e 's/_/ /g' | tr [a-m][n-z] [n-z][a-m]
> 
>  
> 
> 
> _______________________________________________
> grass-user mailing list
> grass-user at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/grass-user

-- 
Tom Russo    KM5VY   SAR502   DM64ux          http://www.swcp.com/~russo/
Tijeras, NM  QRPL#1592 K2#398  SOC#236        http://kevan.org/brain.cgi?DDTNM
 echo "prpv_a'rfg_cnf_har_cvcr" | sed -e 's/_/ /g' | tr [a-m][n-z] [n-z][a-m]