update with pyogrio

requirements.txt

@@ -1 +1,2 @@
 geopandas >= 0.11.0
+pyogrio >= 0.5.1

shp2gpkg.py

@@ -1,48 +1,64 @@
-#!/usr/bin/env python3
+#!/usr/bin/env python
 """Grovel directories containing ESRI ShapeFile to create GeoPackage layers"""
 
 from os import walk
+import os
 import re
 import argparse
 import pandas as pd
-import geopandas as gp
 
+os.environ["USE_PYGEOS"] = "0"
+from pyogrio import read_dataframe, write_dataframe
 
-pd.set_option('display.max_columns', None)
+pd.set_option("display.max_columns", None)
 
-PARSER = argparse.ArgumentParser(description='Create ESRI shp data as layers in a GeoPackage file')
-PARSER.add_argument('--output', type=str, default='geo-model.gpkg', nargs='?',
-                    help='output filename')
-PARSER.add_argument('--search', type=str, default='.', nargs='?', help='search path')
-PARSER.add_argument('--crs', type=str, default='EPSG:32630',
-                    nargs='?', help='coordinate reference system')
+PARSER = argparse.ArgumentParser(
+    description="Create ESRI shp data as layers in a GeoPackage file"
+)
+PARSER.add_argument(
+    "--output", type=str, default="geo-model.gpkg", nargs="?", help="output filename"
+)
+PARSER.add_argument("--search", type=str, default=".", nargs="?", help="search path")
+PARSER.add_argument(
+    "--crs",
+    type=str,
+    default="EPSG:32630",
+    nargs="?",
+    help="coordinate reference system",
+)
 
 ARGS, REST = PARSER.parse_known_intermixed_args()
 REST = (REST + [None] * 3)[:3]
 OUTPATH, FILEPATH, CRS = [(i or j) for i, j in zip(REST, vars(ARGS).values())]
 
+
 def list_files(filepath, match):
-    """find all filenames containing match in directories under filepath """
+    """find all filenames containing match in directories under filepath"""
     files = ()
-    for (d, _, filenames) in walk(filepath):
-        if 'geopandas/datasets' in d:
+    for d, _, filenames in walk(filepath):
+        if "geopandas/datasets" in d:
             continue
-        files = files + tuple(f'{d}/{f}' for f in filenames if match in f)
+        files = files + tuple(f"{d}/{f}" for f in filenames if match in f)
     return files
 
-PATTERNS = [re.compile(i, re.IGNORECASE) for i in ['^shapefile', 'shape$', 'file$']]
+
+PATTERNS = [re.compile(i, re.IGNORECASE) for i in ["shapefile", "shape$", "file$"]]
+
 
 def get_layername(filepath):
     """return layer name from shape-filepath"""
-    r = filepath.split('/')[-1]
-    r = r.replace('.shp', '')
+    r = filepath.split("/")[-1]
+    r = r.replace(".shp", "")
     for p in PATTERNS:
-        r = [i for i in p.split(r) if i != ''][0]
+        r = p.split(r)[0]
     return r
 
-FILES = [f for f in list_files(FILEPATH, 'shp') if f[-4:] == '.shp']
+
+FILES = [f for f in list_files(FILEPATH, "shp") if f[-4:] == ".shp"]
+
 
 for f in FILES:
-    gf = gp.read_file(f, engine='pyogrio')
+    gf = read_dataframe(f)
     layername = get_layername(f)
-    gf.to_crs(CRS).to_file(OUTPATH, driver='GPKG', layer=layername, engine='pyogrio')
+    print(layername)
+    write_dataframe(gf.to_crs(CRS), OUTPATH, layer=layername)