[GRASS-user] problem with v.proj building topology

Hi List,

I have a problem using v.proj. In in the destination location the vector map is incomplete, it seems that the topology is not completely built. In particular, there are missing polygons. The command line:

v.proj input=suolo_paes_ext location=UTM_32_ED50 mapset=mimmo output=suolo_paes_ext --o

gives the following WARNING:

Building topology for vector map <suolo_paes_ext>…
Registering primitives…
5305 primitives registered
26646 vertices registered
Building areas…
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored
100%
462 areas built
35 isles built
Attaching islands…
100%
Attaching centroids…
100%
Number of nodes: 4665
Number of primitives: 5305
Number of points: 0
Number of lines: 0
Number of boundaries: 4934
Number of centroids: 371
Number of areas: 462
Number of isles: 35
Number of incorrect boundaries: 4116
Number of centroids outside area: 224
Number of duplicate centroids: 12
Number of areas without centroid: 327

I’m using Debian, and tried both with grass6.4 debian package and grass6.5 svn

I also tried to use v.build after that, without results.

Any hints?

Thanks

-Margherita

Hi Margherita,

On Sat, Dec 4, 2010 at 8:18 PM, Margherita Di Leo <diregola@gmail.com> wrote:

Hi List,

I have a problem using v.proj. In in the destination location the vector map
is incomplete, it seems that the topology is not completely built. In
particular, there are missing polygons. The command line:

v.proj input=suolo_paes_ext location=UTM_32_ED50 mapset=mimmo
output=suolo_paes_ext --o

gives the following WARNING:

Building topology for vector map <suolo_paes_ext>...
Registering primitives...
5305 primitives registered
26646 vertices registered
Building areas...
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored

...

You are reprojecting from UTM_32_ED50 into ...?

I also tried to use v.build after that, without results.

In any case, you may need to clean up the topology *before* running v.proj.
Perhaps the original map isn't clean.

Markus

Hi Markus,

Thanks for answering. I solved using:

ogr2ogr -s_srs EPSG:23032 -t_srs EPSG:32632 repr_suolo_paes.shp suolo_paes_ext.shp

and imported the reprojected map with v.in.ogr. All polygons are reprojected now.
But I noticed that using v.proj the map was reprojected with an important error .I have another map of the same area, when I try to overlay them they are in two different places, very far apart. But when I overlay the map reprojected with ogr, everythings looks OK.
Couldn’t it be a bug of v.proj? Has anyone had a similar problem using v.proj?

Thank you

Margherita

2010/12/5 Markus Neteler <neteler@osgeo.org>

Hi Margherita,

On Sat, Dec 4, 2010 at 8:18 PM, Margherita Di Leo <diregola@gmail.com> wrote:

Hi List,

I have a problem using v.proj. In in the destination location the vector map
is incomplete, it seems that the topology is not completely built. In
particular, there are missing polygons. The command line:

v.proj input=suolo_paes_ext location=UTM_32_ED50 mapset=mimmo
output=suolo_paes_ext --o

gives the following WARNING:

Building topology for vector map <suolo_paes_ext>…
Registering primitives…
5305 primitives registered
26646 vertices registered
Building areas…
WARNING: Area of size = 0.0 ignored
WARNING: Area of size = 0.0 ignored

You are reprojecting from UTM_32_ED50 into …?

I also tried to use v.build after that, without results.

In any case, you may need to clean up the topology before running v.proj.
Perhaps the original map isn’t clean.

Markus