2024-12-08 16:59:46,248 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(536) INFO Waiting for raster stats worker result.
2024-12-08 16:59:47,099 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘FID_沿海’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,099 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘省’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,100 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘省代码’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,100 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘市’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,100 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘市代码’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,100 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘市类型’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,100 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘省类型’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,101 (osgeo) utils._log_gdal_errors(98) ERROR [errno 1] Failed to create field name ‘面积’: cannot convert to ISO-8859-1
2024-12-08 16:59:47,132 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] One or several characters couldn’t be converted correctly from UTF-8 to ISO-8859-1. This warning will not be emitted anymore.
2024-12-08 16:59:48,070 (osgeo) utils._log_gdal_errors(98) WARNING [errno 3] Cannot find tms_NZTM2000.json (GDAL_DATA is not defined)
2024-12-08 16:59:48,447 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type.
2024-12-08 16:59:50,073 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2264) INFO starting reprojection
2024-12-08 16:59:50,084 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer disjoint_vector of geometry type POLYGON, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type.
2024-12-08 16:59:50,148 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2304) INFO reprojection 100.0% complete on reprojected.gpkg
2024-12-08 16:59:50,188 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1798) INFO Clipping rasters to their intersection with the vector
@Cityhunter , welcome to the community. Please read the Posting Guidelines and add more context to your post, including an attachment with the entire logfile. Thank you,