[GRASS-stats] Loading a point-vector table with 466 columns

[Markus Neteler]

On Mon, May 25, 2009 at 8:08 AM, Hamish <hamish_b at yahoo.com> wrote:
>
> Roger wrote:
>> >> The plugin also creates a fictitious third dimension in
>> >> (point at least) data that has created havoc, and has led
>> >> to readVECT6() getting a pointDropZ= argument - that's why it
>> >> says that wkbPoint is 3 with the plugin and (correctly) 2
>> >> otherwise.
> ...
>> OK, looks like calls to Vect_is_3d(poMap) missing in
>> *OGRGRASSLayer::GetFeatureGeometry about lines after 851 in
>> ogf_frmts/grass/ogrgrasslayer.cpp; the problem exists for
>> all geometry types, just emitting a z value even if the vect is
>> 2D. I can't see that Vect_read_line returns z as NULL if not 3D -
>> Vect_is_3d() is not used much in lib/vector/*.
>
> filed as gdalbug #3009:  https://trac.osgeo.org/gdal/ticket/3009

According to GDAL-trac now fixed:

On Mon, May 25, 2009 at 10:50 PM, GDAL <trac at osgeo.org> wrote:
> #3009: GRASS OGR driver: 2D maps read as 3D
> --------------------+-------------------------------------------------------
>  Reporter:  hamish  |        Owner:  rouault
>     Type:  defect  |       Status:  closed
>  Priority:  normal  |    Milestone:  1.7.0
> Component:  OGR_SF  |      Version:  unspecified
>  Severity:  normal  |   Resolution:  fixed
>  Keywords:  grass   |
> --------------------+-------------------------------------------------------
> Changes (by rouault):
>
>  * status:  assigned => closed
>  * resolution:  => fixed
>  * milestone:  => 1.7.0
>
> Comment:
>
>  For the 2D geometries reported as 3D, I've commited a fix in trunk in
>  r17122.
>
> --
> Ticket URL: <http://trac.osgeo.org/gdal/ticket/3009#comment:4>
> GDAL <http://trac.osgeo.org/gdal/>
> Geospatial Data Abstraction Library is a translator library for raster and vector geospatial data formats.


Kudos to Even Rouault (and the bug hunters here of course),
Markus