[gdal-dev] GDAL COG - Problem with Overviews

[Rahkonen Jukka (MML)]

Hi,

Please use “reply to all” for including the list as a recipient.
More details would lead to less guessing. Please add gdalinfo about the source and result (georeferencing is not important), tell how did you get the image of the original (which bands were selected, was any color enhancement applied). QGIS and a look at Layer properties – Symbology is good for that. Test with gdal_translate before going to Python and send us the command if you face the same trouble. If you cannot solve the problem yourself send a link to test image. It does not need to be your real image, any technically similar image is good and usually even better if you cut it smaller. A 200x200 pixel sized piece of the original image without georeferencing is one alternative that probably would not reveal secrets.

-Jukka Rahkonen-

Lähettäjä: Nicole Kamp <niki.kamp at gmail.com>
Lähetetty: torstai 25. marraskuuta 2021 15.22
Vastaanottaja: Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi>
Aihe: Re: [gdal-dev] GDAL COG - Problem with Overviews

Hello!
Thanks. I tried it with 0,1,2 - it did not work.


On Thu, Nov 25, 2021 at 1:49 PM Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi<mailto:jukka.rahkonen at maanmittauslaitos.fi>> wrote:
Hi,

Maybe the band list be zero based? Try [0,1,2].

-Jukka Rahkonen-

Lähettäjä: gdal-dev <gdal-dev-bounces at lists.osgeo.org<mailto:gdal-dev-bounces at lists.osgeo.org>> Puolesta Nicole Kamp
Lähetetty: torstai 25. marraskuuta 2021 14.44
Vastaanottaja: gdal-dev at lists.osgeo.org<mailto:gdal-dev at lists.osgeo.org>
Aihe: [gdal-dev] GDAL COG - Problem with Overviews

Dear GDAL Developers,
I have an overview problem with my Cloud Optimized GeoTIFF.
I tried to translate a 16 bit, 4 bands ortho image to a COG 8bit, 3 bands image with JPEG compression.

That is the original image.<https://pasteboard.co/CLur6VdD0CwQ.jpg>
That is the result.<https://pasteboard.co/MJv3lOjlI8iM.jpg>

It works fine with LZW compression, but I need the JPEG compression. What am I doing wrong?

Thanks,
Niki


Here is my code.
        for raster in glob.glob(str(parameters[self.INPUT_Pfad])+'/*.tif'):
            fileInfo = QFileInfo(raster)
            baseName = fileInfo.baseName()
            file_name = str(parameters[self.INPUT_Pfad])+'/'+baseName+'.tif'
            feedback.pushInfo(file_name)
            output_file1 = str(parameters[self.OUTPUT_Pfad])+'/'+baseName+'.tif'
            output_file2 = str(parameters[self.OUTPUT_Pfad])+'/'+baseName+'_cog.tif'

            # Create COG
            x_size = 6250
            y_size = 5000
            num_bands = 4

            src_ds = gdal.Open(file_name)
            format = "GTiff"
            driver = gdal.GetDriverByName(format)

            # Translate
            translate_options = gdal.TranslateOptions(format = 'COG', outputType = gdal.GDT_Byte, bandList =                                                     [1,2,3],  scaleParams=[''],
                                          creationOptions = ['TILED=YES','COMPRESS=JPEG','BIGTIFF=IF_NEEDED',                                                      'BLOCKXSIZE=512', 'BLOCKYSIZE=512', 'RESAMPLING=NEAREST',                                                                'OVERVIEWS=IGNORE_EXISTING', 'INTERLEAVE=PIXEL']
                                          )
            gdal.Translate(output_file1, src_ds, options=translate_options)

            translate_options = None
            src_ds = None
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20211125/ca1f343f/attachment-0001.html>