[gdal-dev] The "Bertin 1953" projection

[Nikos Alexandris]

* Even Rouault <even.rouault at spatialys.com> [2018-11-14 12:13:50 +0100]:

>> So, `ogr2ogr` fails to pick up "Bertin 1953":
>> ```
>> ogr2ogr -t_srs +proj=bertin1953 TM_WORLD_BORDERS_SIMPL-0.3_Bertin1953.shp
>> TM_WORLD_BORDERS_SIMPL-0.3.shp ERROR 1: Failed to process SRS definition:
>> +proj=bertin1953
>> ```
>
>You need to add the +wktext parameter for use in GDAL:
>	  -t_srs "+proj=bertin1953 +wktext"

It works:

```
ogr2ogr -t_srs "+proj=bertin1953 +wktext" -f GPKG TM_WORLD_BORDERS_SIMPL-0.3_Bertin1953.gpkg TM_WORLD_BO
RDERS_SIMPL-0.3.shp
Warning 1: A geometry of type MULTIPOLYGON is inserted into layer TM_WORLD_BORDERS_SIMPL-0.3 of geometry t
ype POLYGON, which is not allowed. This warning will no longer be emitted for this combination of layer an
d feature geometry type.
```

and
```
ogrinfo TM_WORLD_BORDERS_SIMPL-0.3_Bertin1953.gpkg -al -so
INFO: Open of `TM_WORLD_BORDERS_SIMPL-0.3_Bertin1953.gpkg'
      using driver `GPKG' successful.

Layer name: TM_WORLD_BORDERS_SIMPL-0.3
Metadata:
  DBF_DATE_LAST_UPDATE=2008-07-30
Geometry: Polygon
Feature Count: 246
Extent: (-14385600.000000, -8557200.000000) - (16159500.000000, 8556070.000000)
Layer SRS WKT:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["unknown",
            SPHEROID["WGS84",6378137,298.257223563]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["custom_proj4"],
    EXTENSION["PROJ4","+proj=bertin1953 +wktext"]]
FID Column = fid
Geometry Column = geom
FIPS: String (2.0)
ISO2: String (2.0)
ISO3: String (3.0)
UN: Integer (0.0)
NAME: String (50.0)
AREA: Integer (0.0)
POP2005: Integer64 (0.0)
REGION: Integer (0.0)
SUBREGION: Integer (0.0)
LON: Real (0.0)
LAT: Real (0.0)
```

Merci Even!

Next would be to make this usable through GRASS GIS or QGIS, for
example.  I'll write in grass-user next for that matter.

Nikos