appelmar / gdalcubes_cpp

Earth observation data cubes from GDAL image collections
MIT License
74 stars 7 forks source link

GDAL3 and PROJ6 issues #34

Closed appelmar closed 4 years ago

appelmar commented 4 years ago

Running gdalcubes with GDAL3 and proj6 currently shows some problems:

  1. The extent of images in collections has swapped x and y axes (due to axis order in EPSG:4326)
  2. Reading image collections with gdalwarp tends to be significantly slower (seems to be related to proj, maybe just some sqlite mutex waits, needs more testing).

This issue aims at collecting further related problems and developing ideas how to solve these.

appelmar commented 4 years ago

For now, the workaround for 1. is to set CPLSetConfigOption("OGR_CT_FORCE_TRADITIONAL_GIS_ORDER", "YES") (see c24da710107f2b580170252f23e9529c281c7327).

appelmar commented 4 years ago

For 2. a potential solution would be to add a wrapper for gdalwarp that calls GDALReprojectionTransform() only once for every unique source SRS of relevant images.

appelmar commented 4 years ago

Fixed in v0.3.0