Running recent version of 5.7. The problem seems to occur while the
vectors are being cleaned. I have attached the shape files. I did
read this vector in fine with an earlier version of 5.7.
SHPFILE="utm_grd"
DBFFILE="utm_grid_fin"
GRASS 5.7.cvs:~/grassalpac/setupfiles > v.in.ogr -o dsn="$SHPDIR"
output=shutin layer="$SHPFILE"
Over-riding projection check.
Proceeding with import...
WARNING: The vector 'shutin' already exists and will be overwritten.
Layer: utm_grd
-----------------------------------------------------
Building topology ...
1168 primitives registered
1168 areas built
1168 isles built
Topology was built.
Number of nodes : 1168
Number of primitives: 1168
Number of points : 0
Number of lines : 0
Number of boundaries: 1168
Number of centroids : 0
Number of areas : 1168
Number of isles : 1168
Number of areas without centroid : 1168
-----------------------------------------------------
WARNING: Cleaning polygons, result is not guaranteed!
Building topology ...
Topology was built.
Number of nodes : 1168
Number of primitives: 1168
Number of points : 0
Number of lines : 0
Number of boundaries: 1168
Number of centroids : 0
Number of areas : -
Number of isles : -
-----------------------------------------------------
Break polygons:
Registering points ... 4858
All points (vertices): 6026
Registered points (unique coordinates): 4858
Points marked for break: 1168
Breaks: 0
-----------------------------------------------------
Remove duplicates:
Duplicates: 0
-----------------------------------------------------
Break boundaries:
Intersections: 6958
-----------------------------------------------------
Remove duplicates:
Duplicates: 721
-----------------------------------------------------
Clean boundaries at nodes:
Modifications: 231ERROR: Attempt to read dead line 13127
GRASS 5.7.cvs:~/grassalpac/setupfiles >
(attachments)
myutm.tar.gz (113 KB)
It is a bug in v.in.ogr (Vect_clean_small_angles_at_nodes),
I don't have a chance to fix it until monday. For now, you can
use v.in.ogr with snap=0.1 (or other reasonable threshold).
Radim
On Wednesday 25 August 2004 22:03, Craig Aumann wrote:
Running recent version of 5.7. The problem seems to occur while the
vectors are being cleaned. I have attached the shape files. I did
read this vector in fine with an earlier version of 5.7.
SHPFILE="utm_grd"
DBFFILE="utm_grid_fin"
GRASS 5.7.cvs:~/grassalpac/setupfiles > v.in.ogr -o dsn="$SHPDIR"
output=shutin layer="$SHPFILE"
Over-riding projection check.
Proceeding with import...
WARNING: The vector 'shutin' already exists and will be overwritten.
Layer: utm_grd
-----------------------------------------------------
Building topology ...
1168 primitives registered
1168 areas built
1168 isles built
Topology was built.
Number of nodes : 1168
Number of primitives: 1168
Number of points : 0
Number of lines : 0
Number of boundaries: 1168
Number of centroids : 0
Number of areas : 1168
Number of isles : 1168
Number of areas without centroid : 1168
-----------------------------------------------------
WARNING: Cleaning polygons, result is not guaranteed!
Building topology ...
Topology was built.
Number of nodes : 1168
Number of primitives: 1168
Number of points : 0
Number of lines : 0
Number of boundaries: 1168
Number of centroids : 0
Number of areas : -
Number of isles : -
-----------------------------------------------------
Break polygons:
Registering points ... 4858
All points (vertices): 6026
Registered points (unique coordinates): 4858
Points marked for break: 1168
Breaks: 0
-----------------------------------------------------
Remove duplicates:
Duplicates: 0
-----------------------------------------------------
Break boundaries:
Intersections: 6958
-----------------------------------------------------
Remove duplicates:
Duplicates: 721
-----------------------------------------------------
Clean boundaries at nodes:
Modifications: 231ERROR: Attempt to read dead line 13127
GRASS 5.7.cvs:~/grassalpac/setupfiles >
No prob!
Thanks Radim!
On Thu, 2004-08-26 at 05:34, Radim Blazek wrote:
It is a bug in v.in.ogr (Vect_clean_small_angles_at_nodes),
I don't have a chance to fix it until monday. For now, you can
use v.in.ogr with snap=0.1 (or other reasonable threshold).
Radim
I have fixed few problems in lib and modules.
It runs better now but it is not perfect, one centroid
cannot be attached to an area because the area is too small.
Some check must be added, probably to Vect_get_point_in_poly_isl.
It is better to run it with snap threshold anyway, because otherwise
there are too many overlaps and holes.
Radim
On Thursday 26 August 2004 17:14, Craig Aumann wrote:
No prob!
Thanks Radim!
On Thu, 2004-08-26 at 05:34, Radim Blazek wrote:
> It is a bug in v.in.ogr (Vect_clean_small_angles_at_nodes),
> I don't have a chance to fix it until monday. For now, you can
> use v.in.ogr with snap=0.1 (or other reasonable threshold).
>
> Radim