#21 Skips outer ring
In rare cases the outer ring is missing when reading
the VMap0 coverage BND. The following example reads a
small piece of sea polygons:
ogdi_info -no-dict -u
gltp:/vrf/c:/geodata/vmap0/sasaus/sasaus -f Area -l
OCEANSEA@BND(*) -r 15.6 14.8 95.0 97.0 -dl
The coordinates for the real second outer ring is
missing. The returned second outer ring is actually
the first hole in second polygon. This has the effect
that all holes are outside what is said to be the
outer ring.
In the enclosed image the blue polygon in the middle
should be a hole in the non existing outer ring. The
outlined polygons should also be holes.
Reading the same data with VPFView 2.1 shows it
correctly.
I am using:
3.1alpha3: Sun Mar 11 17:20:41 EST 2001
running at W2000
Logged In: YES
user_id=720294
It seems like removing a test in vrf_get_ring_coords in
driver/vrf/feature.c handles this problem. The following test
looks like an optimization and when removed the missing
polygons appear.
if (edge_rec.start_node == edge_rec.end_node)
done = TRUE;
Logged In: YES
user_id=720294
It seems like removing a test in vrf_get_ring_coords in
driver/vrf/feature.c handles this problem. The following test
looks like an optimization and when removed the missing
polygons appear.
if (edge_rec.start_node == edge_rec.end_node)
done = TRUE;