* Even Rouault <[email protected]> [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
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
