outLayer = (
"ogr:dbname=\'" + gpkgPath + "\' table=\"emodnet_bathymetry\" (geom)"
)
params = {
"INPUT": tempDir + "bathymetry_pts_osgb.shp",
"MASK": gpkgPath + "|layername=study_area_water",
"OUTPUT": outLayer
}
processing.run("gdal:clipvectorbypolygon", params, feedback=feedback)
Current workaround (using GeoPandas):
import geopandas as gpd
data = gpd.read_file(
tempDir + "bathymetry_pts_osgb.shp",
mask=gpd.read_file(gpkgPath, layer="study_area_water")
)
data.to_file(gpkgPath, layer="emodnet_bathymetry", driver="GPKG")