[gdal-dev] Puzzled about gdal_translate return all white image...

2021-08-20 Thread Carl Godkin
Hi, I'm very puzzled about why gdal_translate (v. 3.3.1) returns an all-white image when I select a smaller area than with a working case with a larger area. These two command lines are exactly the same except that the failing case has a smaller area: gdal_translate " https://basemap.nationalmap.

Re: [gdal-dev] Custom vertical transformations

2021-08-20 Thread Even Rouault
Regarding your immediate issue, the workaround I see would be that you use PROJ to create a regular grid that contains the vertical correction, use gdalwarp to transform it to the exact extent and resolution of the result of your 2D-only gdalwarp, and finally use gdal_calc to sum the two.

Re: [gdal-dev] Custom vertical transformations

2021-08-20 Thread Even Rouault
Kristian, your question kind of intersects https://lists.osgeo.org/pipermail/gdal-dev/2021-August/054529.html Two part answer: A) I guess that there's an issue with axis order and units in your pipeline. In particular if your DHYMSEA_1km_6150_542.tif has a EPSG geographic CRS set, then the

[gdal-dev] Custom vertical transformations

2021-08-20 Thread Kristian Evers via gdal-dev
Hi list, I am trying to apply a custom vertical transformation to a grid using gdalwarp and am struggling a bit to get the results I want. It's entirely possible I am trying something not within the capabilities of gdalwarp but I am sure someone here can help me figure that out. I am devising

Re: [gdal-dev] Define and transform a GeoTiff with a compound and projected CRS

2021-08-20 Thread Even Rouault
Hannah, there are 2 issues. * First your WKT is indeed invalid. Instead of using a WKT, as there are EPSG codes for the horizontal and vertical part, you could just use gdal_translate -a_srs EPSG:6341+5703  (EPSG:6341 is "NAD83(2011) / UTM zone 12N", which is slightly more appropriate given th