[GRASSLIST:5443] Still problems with vectors in grass6 (long msg)
Carlos Henrique Grohmann
guano at usp.br
Wed Jan 19 14:10:35 EST 2005
I'm still with problems with my vect maps, so now I'll post all msgs.
I have this USGS Open File Report (97-470D, 1999) with a shape file of the
political divisions of South America. I want to import it and make a vect map
of South America without divisions. So I've imported it like this:
>v.in.ogr dsn=/mnt/cdrw_dvdrom/SAM/VIEWS/GEO output=cont5 layer=POL6_G
min_area=0.0001 >type=boundary snap=-1 -o
>Over-riding projection check.
>Proceeding with import...
>Layer: POL6_G
>-----------------------------------------------------
>Building topology ...
>Registering lines:
>337 primitives registered
>Building areas:
>
>196 areas built
>177 isles built
>Attaching islands:
>Attaching centroids:
>Topology was built.
>Number of nodes : 319
>Number of primitives: 337
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 337
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 177
>Number of incorrect boundaries : 93
>Number of areas without centroid : 196
>-----------------------------------------------------
>
*** stuff ***
>
>Building topology ...
>Building areas:
>
>196 areas built
>182 isles built
>Attaching islands:
>Topology was built.
>Number of nodes : 319
>Number of primitives: 341
>Number of points : 0
>Number of lines : 4
>Number of boundaries: 337
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Layer: POL6_G
>-----------------------------------------------------
>-----------------------------------------------------
>Building topology ...
>Registering lines:
>337 primitives registered
>Building areas:
>
>196 areas built
>182 isles built
>Attaching islands:
>
And I got my map, with all countries in South America. Then I try to dissolve
all commom boundaries with:
>v.extract input=cont5 output=cont6 type=boundary layer=1 new=1 -d
>Building topology ...
>Registering lines:
>333 primitives registered
>Building areas:
>
>196 areas built
>182 isles built
>Attaching islands:
>Attaching centroids:
>Topology was built.
>Number of nodes : 319
>Number of primitives: 333
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 333
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Writing attributes ...
>Removing duplicate centroids ...Building topology ...
>Registering lines:
>
>333 primitives registered
>Building areas:
And I got the same map! No dissolved commom boundaries..
So I did:
>v.extract input=cont5 output=cont6 type=boundary layer=1 new=1 -d
>Building topology ...
>Registering lines:
>333 primitives registered
>Building areas:
>
>196 areas built
>182 isles built
>Attaching islands:
>Attaching centroids:
>Topology was built.
>Number of nodes : 319
>Number of primitives: 333
>Number of points : 0
>Number of lines : 0
>Number of boundaries: 333
>Number of centroids : 0
>Number of areas : 196
>Number of isles : 182
>Number of areas without centroid : 196
>
>Writing attributes ...
>Removing duplicate centroids ...Building topology ...
>Registering lines:
>
>333 primitives registered
>Building areas:
All the same, the boundaries are still there.
Once I get an error like this:
>bgerror failed to handle background error.
> Original error: too many nested calls to Tcl_EvalObj (infinite loop?)
> Error in bgerror: too many nested calls to Tcl_Eval (infinite loop?)
and then GRASS crashed.
Pleas help. I don't know what to do.
I'll leave the shp file in
http://www.igc.usp.br/pessoais/guano/sam.html
thanks
--
+----------------------------------------------------------+
Carlos Henrique Grohmann - Guano
Geologist M.Sc - Doctorate Student at IGc-USP - Brazil
Linux User #89721 - guano at usp dot br
+----------------------------------------------------------+
More information about the grass-user
mailing list