Re: Error on gdal raster info in VRT - PROJ: No GeoPixelScale tag in overviews

Michael Smith via gdal-dev <[email protected]>
Newsgroups gmane.comp.gis.gdal.devel
Message-ID <[email protected]> 
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&gt;> <http://www.spatialys.com> <http://www.spatialys.com&gt;> <http://www.spatialys.com&gt;> <http://www.spatialys.com&amp;gt;&gt;>
My software is free, but my time generally not.
























_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
lmpx.com only provides a reader for public news (NNTP) servers. It is not affiliated with the servers or forums shown here and is not responsible for the content of articles, which is written by their respective authors.