PyQGIS – Copying a Geopackage to Another Folder

geopackagepyqgis

I try to copy a GPKG to another folder. The result is the correct named copied GPKG but there is only one layer in the copy. It’s the layer with the same name as the GPKG. The other features (1x polygon, 1x no geometry) from the original are not copied. The original GPKG is open with QGIS when I start the copy-process. Is there another simple way to copy the whole GPKG or do I have to create a new GPKG and copy the layers of the original GPKG to it?

from shutil import copy
inputgpkg = r'C:\Testfolder\original.gpkg'
outputgpkg = r'C:\Testfolder2\copy.gpkg'
copy(inputgpkg,outputgpkg)

Best Answer

As user30184 stated, one needs to copy the -wal and -shm file as well:

from shutil import copy
inputgpkg = r'C:\Testfolder\original.gpkg'
outputgpkg = r'C:\Testfolder2\copy.gpkg'
copy(inputgpkg,outputgpkg)
copy(inputgpkg+'-wal',outputgpkg+'-wal')
copy(inputgpkg+'-shm',outputgpkg+'-shm')

Related Solutions

[GIS] Using folder wildcards to copy files with Python

Try this:

import os, glob, shutil
root_dir = "C:\Users\xxxx\Desktop\Test\\"
country_dir = "Country_"
grid_path = "C:\Users\xxxx\Desktop\Test\Grid\Grid.shp"

# Get all files that constitute the Grid Shapefile
gridShpFiles = glob.glob(grid_path[:-3]+"*")

for path,dirname,files in os.walk(root_dir):
  if country_dir in path:
    for f in gridShpFiles:
      shutil.copy(f, path)

Take your time to set paths (lines 2, 3, and 4), I tested the script on GNU/Linux, but Windows paths are always trickier.

As you notice, you can avoid the wildcard; the expression if country_dir in path: does the work.

Finally, as @atrwork21 mentioned, you need to take all Shapefile files into account. That's what the second block does.

[GIS] Remove a raster from Geopackage

You can remove raster from GeoPackage with ogrinfo but it is a bit tricky because each raster table has also references in several metadata tables.

EDIT

GeoPackage has foreign key constraints which can prevent user from deleting some data. There are two alternatives to make deletes to success:

Delete stuff in right order, or
Disable foreign key check from SQLite database before deletes.

Here are ogrinfo commands ordered so that foreign keys don't make any trouble.

ogrinfo -sql "delete from gpkg_tile_matrix_set where table_name='raster2'" droptest.gpkg
ogrinfo -sql "delete from gpkg_tile_matrix where table_name='raster2'" droptest.gpkg
ogrinfo -sql "delete from gpkg_metadata_reference where table_name='raster2'" droptest.gpkg
ogrinfo -sql "delete from gpkg_contents where table_name='raster2'" droptest.gpkg
ogrinfo -sql "drop table raster2" droptest.gpkg

In my original answer I deleted things from the database in alphabetic order and I had to disable the foreign key check with a specific GDAL config option.

Drop the tiles

ogrinfo -sql "drop table raster2" droptest.gpkg

Remove line from gpkg_contents

ogrinfo -sql "delete from gpkg_contents where table_name='raster2'" droptest.gpkg

Remove references from gpkg_metadata_reference table. Ogrinfo does not accept the command with the same syntax as above because the SQLite database is now in invalid state because of orphaned foreign key constraints and you must use a special GDAL configuration option to override the validity check.

ogrinfo -sql "delete from gpkg_metadata_reference where table_name='raster2'" droptest.gpkg --config OGR_GPKG_FOREIGN_KEY_CHECK OFF

Remove references from gpkg_tile_matrix table.

ogrinfo -sql "delete from gpkg_tile_matrix where table_name='raster2'" droptest.gpkg --config  OGR_GPKG_FOREIGN_KEY_CHECK OFF

Remove references from gpkg_tile_matrix_set table

ogrinfo -sql "delete from gpkg_tile_matrix_set where table_name='raster2'" droptest.gpkg --config OGR_GPKG_FOREIGN_KEY_CHECK OFF

Now the raster table and all the references have been deleted. As a final step it might be good to vacuum the database.

ogrinfo -sql "VACUUM" droptest.gpkg

However, if the database is very big VACUUM will be slow. It could actually be faster to copy all the raster layers except the one that you want to drop into a new GeoPackage with gdal_translate. If you have also lots of vector layers in your GeoPackage that option may not be so fascinating. For simplifying the deletion you can save all the ogrinfo commands into a batch file and give the name of the raster file that you want to drop as a variable.

Best Answer

Related Solutions

[GIS] Using folder wildcards to copy files with Python

[GIS] Remove a raster from Geopackage

Related Question