Re: [gdal-dev] How proceed with GDAL GeoTIFF files with multiply subdataset

2023-03-08 Thread Brad Hards
On Wednesday, 8 March 2023 11:54:39 PM AEDT Raivo Rebane wrote: > Hi > > I found from documentation following about subdataset: > > > Subdatasets > > Multi-page TIFF files are exposed as subdatasets. On opening, a > subdataset name

Re: [gdal-dev] How proceed with GDAL GeoTIFF files with multiply subdataset

2023-03-08 Thread Raivo Rebane
Hi now I am using following c++ commands to read pages from multisubdatasets Geotiff file :    dataset = static_cast(GDALOpen(TiffFullName.c_str(), GA_ReadOnly));    if (dataset == NULL) {        std::cout << "Failed to open subdataset nr = " << subsets << " Read command =  "

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Even Rouault
The file is indeed fine and was recognized in earlier GDAL versions (<= 3.4). Bug fixed per https://github.com/OSGeo/gdal/pull/7382 Le 08/03/2023 à 14:22, Micha Silver a écrit : Inline... On 08/03/2023 14:33, Even Rouault wrote: Ah, the issue here is that the GeoTIFF keys in the .tif part

Re: [gdal-dev] How proceed with GDAL GeoTIFF files with multiply subdataset

2023-03-08 Thread Rahkonen Jukka
Hi, Do you know that your source file is a multipage TIFF? Could you share a small test file? There is python test for reading a multipage TIFF at https://github.com/OSGeo/gdal/blob/master/autotest/gcore/tiff_read.py#L582, please have a look. -Jukka Rahkonen- Lähettäjä: gdal-dev Puolesta Rai

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Micha Silver
Inline... On 08/03/2023 14:33, Even Rouault wrote: Ah, the issue here is that the GeoTIFF keys in the .tif partially define the projection (or possibly the .tif.aux.xml contain this invalid WKT). But METHOD["unnamed"] means that

Re: [gdal-dev] question about python installation

2023-03-08 Thread Greg Troxel
Michael Sumner writes: > I've been trying with PYTHONPATH and sys.path, and I think I was doing it > wrong ... but the simplicity of the /usr prefix had not occurred to me, > I'll be doing that. > > If I come back around with further insights I'll share them, but this is > working for me now. >

[gdal-dev] How proceed with GDAL GeoTIFF files with multiply subdataset

2023-03-08 Thread Raivo Rebane
Hi I found from documentation following about subdataset: Subdatasets Multi-page TIFF files are exposed as subdatasets. On opening, a subdataset name is GTIFF_DIR:{index}:filename.tif, where {index} starts at 1. I am using foll

Re: [gdal-dev] Can Geotiff file contain multiply Ortofotos

2023-03-08 Thread Raivo Rebane
I found from documentation following about subdataset: Subdatasets Multi-page TIFF files are exposed as subdatasets. On opening, a subdataset name is GTIFF_DIR:{index}:filename.tif, where {index} starts at 1. I am using followin

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Even Rouault
Ah, the issue here is that the GeoTIFF keys in the .tif partially define the projection (or possibly the .tif.aux.xml contain this invalid WKT). But METHOD["unnamed"] means that the projection method is unknown, so it is logical that GDAL/PROJ can't compute geographic coordinates from projected

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Micha Silver
Hi Even: Here's the full gdalinfo, with the reported CRS: gdalinfo HiBar_21022022_A1.tif Driver: GTiff/GeoTIFF Files: HiBar_21022022_A1.tif   HiBar_21022022_A1.tif.aux.xml

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Javier Jimenez Shaw
gdalinfo is usually showing the coordinates of the corners and center in the CRS of the GeoTIFF, and also in WGS84 in degrees as conversion of the previous. Looks like your CRS does not support such conversion. If you paste the full output of gdalinfo there should be some clues about the CRS used.

Re: [gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Even Rouault
Micha, which CRS is reported and which PROJ version do you use ? It is presumably one projection method for which PROJ only implements the forward conversion, that is from geographic to projected coordinates, but not the reverse way Even Le 08/03/2023 à 09:37, Micha Silver a écrit : I

[gdal-dev] ERROR 1: No inverse operation

2023-03-08 Thread Micha Silver
I received some Geotiff files, and gdalinfo shows: Corner Coordinates: ERROR 1: No inverse operation Upper Left  (  205370.241,  418790.721)   ERROR 1: No inverse operation Lower Left  (  205370.241,  4156