From 3ed9a613d8da9874f45d2a069686f1fa38c72ad9 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Mon, 5 Oct 2020 15:40:20 +0200 Subject: Database: use extended description for extent/area of use, as done by epsg.org WKT export --- src/4D_api.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src/4D_api.cpp') diff --git a/src/4D_api.cpp b/src/4D_api.cpp index 75fa6d04..d7fe13e1 100644 --- a/src/4D_api.cpp +++ b/src/4D_api.cpp @@ -1233,7 +1233,7 @@ std::vector pj_create_prepared_operations(PJ_CONTEXT *ctx, &west_lon, &south_lat, &east_lon, &north_lat, &areaName) ) { const bool isOffshore = - areaName && strstr(areaName, "offshore"); + areaName && strstr(areaName, "- offshore"); if( west_lon <= east_lon ) { op = add_coord_op_to_list(i, op, -- cgit v1.2.3 From 67af1d5741da08d9546b907e0da2c21c54c61b27 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Thu, 8 Oct 2020 15:49:30 +0200 Subject: Make proj_create_crs_to_crs() / cs2cs work with CRS with DatumEnsemble --- src/4D_api.cpp | 34 ++++++++++++---------------------- 1 file changed, 12 insertions(+), 22 deletions(-) (limited to 'src/4D_api.cpp') diff --git a/src/4D_api.cpp b/src/4D_api.cpp index d7fe13e1..3c6ab802 100644 --- a/src/4D_api.cpp +++ b/src/4D_api.cpp @@ -1069,28 +1069,18 @@ static PJ* create_operation_to_geog_crs(PJ_CONTEXT* ctx, const PJ* crs) { geodetic_crs_type == PJ_TYPE_GEOGRAPHIC_3D_CRS ) { auto datum = proj_crs_get_datum(ctx, geodetic_crs); - if( datum ) - { - auto cs = proj_create_ellipsoidal_2D_cs( - ctx, PJ_ELLPS2D_LONGITUDE_LATITUDE, nullptr, 0); - auto ellps = proj_get_ellipsoid(ctx, datum); - proj_destroy(datum); - double semi_major_metre = 0; - double inv_flattening = 0; - proj_ellipsoid_get_parameters(ctx, ellps, &semi_major_metre, - nullptr, nullptr, &inv_flattening); - auto temp = proj_create_geographic_crs( - ctx, "unnamed crs", "unnamed datum", - proj_get_name(ellps), - semi_major_metre, inv_flattening, - "Reference prime meridian", 0, nullptr, 0, - cs); - proj_destroy(ellps); - proj_destroy(cs); - proj_destroy(geodetic_crs); - geodetic_crs = temp; - geodetic_crs_type = proj_get_type(geodetic_crs); - } + auto datum_ensemble = proj_crs_get_datum_ensemble(ctx, geodetic_crs); + auto cs = proj_create_ellipsoidal_2D_cs( + ctx, PJ_ELLPS2D_LONGITUDE_LATITUDE, nullptr, 0); + auto temp = proj_create_geographic_crs_from_datum( + ctx, "unnamed crs", datum ? datum : datum_ensemble, + cs); + proj_destroy(datum); + proj_destroy(datum_ensemble); + proj_destroy(cs); + proj_destroy(geodetic_crs); + geodetic_crs = temp; + geodetic_crs_type = proj_get_type(geodetic_crs); } if( geodetic_crs_type != PJ_TYPE_GEOGRAPHIC_2D_CRS ) { -- cgit v1.2.3