create polygon cut script

We want to cut

[D:\repositories\zonal_stats_toolkit\snappdata\public_lands_only_fixed.gpkg](file:///D:/repositories/zonal_stats_toolkit/snappdata/public_lands_only_fixed.gpkg)

against

[D:\Users\richp\Downloads\gz_2010_us_050_00_5m\gz_2010_us_050_00_5m\gz_2010_us_050_00_5m.shp](file:///D:/Users/richp/Downloads/gz_2010_us_050_00_5m/gz_2010_us_050_00_5m/gz_2010_us_050_00_5m.shp)