So no error if its in the source projection: GTiff: ScanDirectories() GTiff: Opened 4320x4320 overview. GTiff: Opened 2160x2160 overview. GTiff: Opened 1080x1080 overview. GTiff: Opened 540x540 overview. GTiff: Opened 270x270 overview. PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded PROJ: pj_open_lib(us_nga_egm08_25.tif): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/us_nga_egm08_25.tif) - succeeded PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded PROJ: pj_open_lib(us_nga_egm08_25.tif): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/us_nga_egm08_25.tif) - succeeded Corner Coordinates: Upper Left ( 351390.250, 3424709.750) (100d33'20.58"W, 30d56'47.38"N)
If if you warp to 4326 in a vrt, then the vrt has errors: GDAL: GDALOpen(/vsis3/ele-gv-o1/16936682/GRID0000/Buckeye/received/pushed_2025_12_05/DMZ2/BuckEye/FinishedProducts/USA/Texas/DesertIce26_202511/14n3513420_20251108/14n3513420_20251108_dsm_1_050_egm08.cog.tif, this=0x55824e090a90) succeeds as GTiff. GDAL: Computing area of interest: -100.556, 30.9075, -100.51, 30.947 PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded PROJ: pj_open_lib(us_nga_egm08_25.tif): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/us_nga_egm08_25.tif) - succeeded PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded PROJ: pj_open_lib(us_nga_egm08_25.tif): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/us_nga_egm08_25.tif) - succeeded PROJ: Using coordinate operation Inverse of UTM zone 14N + Inverse of WGS 84 to EGM2008 height (2) PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded Warning 1: PROJ: No GeoPixelScale tag PROJ: Using coordinate operation Inverse of UTM zone 14N + Inverse of WGS 84 to EGM2008 height (2) PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded Warning 1: PROJ: No GeoPixelScale tag PROJ: Using coordinate operation Inverse of UTM zone 14N + Inverse of WGS 84 to EGM2008 height (2) PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded PROJ: Using coordinate operation Inverse of UTM zone 14N + Inverse of WGS 84 to EGM2008 height (2) PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded Creating output file that is 9362P x 8071L. GDAL: GDALDriver::Create(VRT,test1.vrt,9362,8071,1,Float32,0x55824e2d2400) -- Michael Smith RSGIS Center – ERDC CRREL NH US Army Corps On 2/19/26, 3:34 PM, "Michael Smith" <[email protected] <mailto:[email protected]>> wrote: GDAL: GDALOverviewDataset(/u02/efs/dged5_error/raster21058922_EPSG:32638+3855.vrt, this=0x55613209db80) creation. PROJ: Using coordinate operation Inverse of UTM zone 38N + Inverse of WGS 84 to EGM2008 height (2) PROJ: pj_open_lib(Und_min1x1_egm2008_isw=82_WGS84_TideFree): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/Und_min1x1_egm2008_isw=82_WGS84_TideFree) - succeeded ERROR 1: PROJ: No GeoPixelScale tag Overviews: 5130x4980, 2565x2490, 1283x1245, 641x623, 320x311 If I remove that file, it hits the standard egm 2008 grid GDAL: GDALOverviewDataset(/u02/efs/dged5_error/raster21058922_EPSG:32638+3855.vrt, this=0x5641ed2c68b0) creation. PROJ: Using coordinate operation Inverse of UTM zone 38N + Inverse of WGS 84 to EGM2008 height (1) PROJ: pj_open_lib(us_nga_egm08_25.tif): call fopen(/home/gridusr/gridpixi/.pixi/envs/default/share/proj/us_nga_egm08_25.tif) - succeeded Overviews: 5130x4980, 2565x2490, 1283x1245, 641x623, 320x311 On 2/19/26, 3:24 PM, "Michael Smith" <[email protected] <mailto:[email protected]> <mailto:[email protected] <mailto:[email protected]>>> wrote: Hmm, that helped a lot. Looks like its tied to our having the Und_min1x1_egm2008_isw=82_WGS84_TideFree file in the PROJ_DATA location. Removing that and the issue goes away. Mike -- Michael Smith RSGIS Center – ERDC CRREL NH US Army Corps On 2/19/26, 3:07 PM, "Even Rouault" <[email protected] <mailto:[email protected]> <mailto:[email protected] <mailto:[email protected]>> <mailto:[email protected] <mailto:[email protected]> <mailto:[email protected] <mailto:[email protected]>>>> wrote: Mike, this error comes from PROJ when trying to use a grid that would lack this GeoPixelScale GeoTIFF tag. Is your EGM2008 grid OK ? Even Le 19/02/2026 à 20:15, Michael Smith via gdal-dev a écrit : > GDAL 3.12.2 "Chicoutimi", released 2026/02/03 - from Conda-forge (both Linux > and OSX) > > I have a base file (COG with overviews) that when I wrap in a VRT to > reproject, I get an odd error but only on Linux, not on OSX (haven’t tested > on windows): > > Specifically, I see: > > Center ( 42.4860100, 15.2722355) ( 42d29' 9.64"E, 15d16'20.05"N) > Band 1 Block=512x128 Type=Float32, ColorInterp=Gray > NoData Value=-32767 > GDAL: GDALOverviewDataset(raster21058922_EPSG:32638+3855.vrt, > this=0x56307ab73b40) creation. > ERROR 1: PROJ: No GeoPixelScale tag > GDAL: GDALOverviewDataset(raster21058922_EPSG:32638+3855.vrt, > this=0x56307ac1acc0) creation. > ERROR 1: PROJ: No GeoPixelScale tag > GDAL: GDALOverviewDataset(raster21058922_EPSG:32638+3855.vrt, > this=0x56307aeaf410) creation. > ERROR 1: PROJ: No GeoPixelScale tag > GDAL: GDALOverviewDataset(raster21058922_EPSG:32638+3855.vrt, > this=0x56307af40750) creation. > ERROR 1: PROJ: No GeoPixelScale tag > GDAL: GDALOverviewDataset(raster21058922_EPSG:32638+3855.vrt, > this=0x56307afc7300) creation. > ERROR 1: PROJ: No GeoPixelScale tag > Overviews: 5130x4980, 2565x2490, 1283x1245, 641x623, 320x311 > GDAL: GDALClose(raster21058922_EPSG:32638+3855.vrt, this=0x56307a786b70) > GDAL: GDALClose(raster26-21058922_EPSG:32638+3855.vrt, this=0x56307a740e50) > > > The vrt that doesn't have reprojection doesn't get this. The base projection > here is UTM+EGM2008. > > I have not seen this error before. > -- http://www.spatialys.com <http://www.spatialys.com> <http://www.spatialys.com> <http://www.spatialys.com>> <http://www.spatialys.com> <http://www.spatialys.com>> <http://www.spatialys.com>> <http://www.spatialys.com&gt;>> My software is free, but my time generally not. _______________________________________________ gdal-dev mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/gdal-dev
