I posted this originally on github, but have moved it here at Sean's request. Additional details of this question/issue are discussed in the `pyproj` repository

It may be that rasterio doesn't support multilayer raster datasets. I'll make a note to open an issue in the tracker about this.

Installed rasterio through conda: conda create -n gpkg-test -c conda-forge python=3.7 rasterio

Hi, Correct, those dataset creation options need to be passed into rasterio.open(). One way to do that is to add to your metadata dict below metadata["profile"] = "GeoTIFF" Any extra keyword

Thanks again Sean, and Evan too. It seems I just need to get the right combination of options. The options that you have shown Even ....how are they implemented? They don't seem to go in the Env like

Indeed. You can play with the GTIFF specific PROFILE creation option (see https://gdal.org/drivers/raster/gtiff.html#creation-issues) PROFILE=BASELINE: .aux.xml only created if presence of

If you are creating a GeoTIFF it is possible that no XML file will be written. As Even Rouault explains in https://lists.osgeo.org/pipermail/gdal-dev/2010-May/024522.html, the PAM file is only created

This didn't work unfortunately Sean. No errors .... but nothing else was wrote out as I was hoping. If you get a few minutes would you mind testing to ensure I'm not going crazy? Thanks,

Hi, Can you tell us how you installed rasterio? If you installed one of the binary wheels from PyPI (which I built), you may be out of luck. I haven't paid any attention to making sure that they have

Sorry, James, I messed up and used the GDAL vernacular. Replace YES with True and I think you'll be good to go.-- Sean Gillies

Thanks for the tip Sean, but that doesn't work for me. with rasterio.Env(TIFF_USE_OVR=True, GDAL_PAM_ENABLED=YES, ESRI_XML_PAM=YES): ----- name 'YES' is not defined If I change 'YES' to True then

Hi James, I have tried this myself, but I think you'll want to add both GDAL_PAM_ENABLED=YES and ESRI_XML_PAM=YES to the

Hello, I'm trying to use rasterio to manipulate a file and then write it out, along with accompanying ESRI-type files. So the output files will be: raster.tif, raster.tif.aux.xml raster.tif.ovr I

I'm trying to add a raster layer to an existing GPKG with rasterio: properties = { 'driver':'GPKG', 'width':array.shape[1], 'height':array.shape[0], 'count':1,

Dion, Jonas: There is a new report in the FIona tracker that may be related: https://github.com/Toblerity/Fiona/issues/761. There I advised turning on CPL_CURL_VERBOSE=YES to see more details about

Dear Dior, Thanks!, this was perfect. I did't get it the first time, but know it runs smoothly. Many thanks! Javier -- Javier Lopatin Institute of Geography and Geoecology Karlsruhe Institute of

As I mentioned, the canonical way to do this is to leverage functools.partial. For example: from functools import partial from concurrent.futures import ProcessPoolExecutor

Hi Dion, Thanks for the tip. I thought that that might be the problem. Any ideas of how to pass a function to executor.map that need several arguments? so far I'm only able to pass functions with a

Hey Javier, this is not really a rasterio problem, mostly a Python peculiarity. In general, you cannot pass functions defined in a local closure to a different process. Just

Hi Guys! I'm still working with the parallelization of raster processes with rasterio, and I came across a problem that you might know about. I've some raster stacks of time series data from which I