<div dir="ltr">Hi<div>This is my first time using the mailing list, I apologize if this question is out of place.</div><div>I am using python 3.6, with gdal python api version 3.6.0</div><div><br></div><div>I am working with pcraster framework, and need to store the resulting simulation files in a geopackage. I found the best way to do this is by converting the pcr object output to numpy array, and store the numpy array directly to gpkg file using:</div><div> </div><div>    creationOptions = ['TILE_FORMAT=PNG', 'APPEND_SUBDATASET=YES', 'PRECISION=0', str('RASTER_TABLE='+str(raster_table))]<br>    outputType= gdal.GDT_Float32    <br></div><div>    driver = gdal.GetDriverByName('GPKG')<br>    outDs = driver.Create(outpath, cols, rows, 1, outputType, options = creationOptions)<br></div><div>    outBand = outDs.GetRasterBand(1)<br>    outBand.SetNoDataValue(-999)<br>    outBand.WriteArray(numpy_array)  # this is where I store my np arrays<br>    outDs.SetGeoTransform(trans)<br>    ds = None<br>    outDs = None<br></div><div><br></div><div>the rows and cols are based on a file I'm using as reference, the output type</div><div>Issues I found:</div><div>1) I'm using the option "Tile format = PNG" because I found this option to be more effective to compress the final output files, I'm working with about 6000 outputs / Is that actual right, or there is another way to have tiff files compressed effectively with output data format as Float32?</div><div>2) As I said I am working with more than 6000 simulation outputs.. I found that my final gpkg created file has a limit of 1000 layers... Is that a default creation option?? How can I make sure my final gpkg file store all the layers?</div><div>3) What is the best way to access a certain via gdal python api the gpkg layers? Or is best to access the individual layer using sqlite3 lib? I use the following script to view list the raster tables, but I want to access a certain layer individually</div><div>ds = gdal.OpenEx('my.gpkg')<br>x = ds.GetMetadata_List('SUBDATASETS')<br></div><div><br></div><div><br clear="all"><div>Thank you!</div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div dir="ltr"></div><br></div></div></div></div>