[GRASS-dev] [bug #4898] (grass) g.proj -p segfaults because of missing datum name

[Moritz Lennert]

Request Tracker wrote:
> this bug's URL: http://intevation.de/rt/webrt?serial_num=4898
> -------------------------------------------------------------------------
> 
> Subject: g.proj -p segfaults because of missing datum name
> 
> Platform: GNU/Linux/x86
> grass obtained from: CVS
> grass binary for platform: Compiled from Sources
> GRASS Version: cvs_head_20060724
> 
> Trying to use the "Use georeferenced file" for new location definition in the start-up screen I get a segfault.
> Searching for the cause, I find that it's in the call to g.proj, as calling g.proj -p on the file gives the following:
> 
>> g.proj -p UrbAdm_Mu.shp
> Trying to open with GDAL...ERROR 4: `UrbAdm_Mu.shp' not recognised as a supported file format.
> 
> Trying to open with OGR......succeeded.
> WARNING: Datum 'unnamed' not recognised by GRASS and no parameters found.
>          Datum transformation will not be possible using this projection
>          information.
> Segmentation fault
> 
> The warning comes from GPJ_osr_to_grass() which is in lib/proj/convert.c. I don't know why it segfaults, though.
> 
> Do we absolutely need a datum name as the information contained in the .prj file seems to be sufficient to create a location ?
> 

I tried introducing the datum name into the prj file and now I get:

A datum name bel72 (Reseau_National_Belge_1972) was specified without 
transformation parameters.
Note that the GRASS default for bel72 is 
towgs84=-99.1,53.3,-112.5,0.419,-0.830,1.885,-1.00.
Segmentation fault

Moritz