2023-09-15 09:33:59,950 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [/Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/Results/InVEST-natcap.invest.crop_production_percentile-log-2023-09-15--09_33_59.txt] 2023-09-15 09:33:59,950 (__main__) cli.main(459) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.crop_production_percentile 3.14.0: aggregate_polygon_path /Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/BaciaCavado/BaciaCavado_fixgeo.shp landcover_raster_path /Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/COS2018_crops_6set.tif landcover_to_crop_table_path /Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/luctocodegrapeecorn.csv model_data_path /Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/model_data n_workers -1 results_suffix CP2018 workspace_dir /Users/mac/Documents/Tese-MBEAG/InVEST_5/Crop_Production/Results 2023-09-15 09:33:59,983 (py.warnings) warnings._showwarnmsg(109) WARNING invest/cli.py:470: DeprecationWarning: `product` is deprecated as of NumPy 1.25.0, and will be removed in NumPy 2.0. Please use `prod` instead. 2023-09-15 09:33:59,999 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(528) INFO Processing crop maize 2023-09-15 09:33:59,999 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(533) INFO Clipping global climate bin raster to landcover bounding box. 2023-09-15 09:34:00,016 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_25th to climate bins. 2023-09-15 09:34:00,027 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,027 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,030 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,033 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate maize yield_25th yield raster to landcover resolution. 2023-09-15 09:34:00,151 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for maize at yield_25th 2023-09-15 09:34:00,165 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,165 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,185 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,186 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,200 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_50th to climate bins. 2023-09-15 09:34:00,210 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,210 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,211 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,211 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,214 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate maize yield_50th yield raster to landcover resolution. 2023-09-15 09:34:00,322 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for maize at yield_50th 2023-09-15 09:34:00,336 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,337 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,356 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,357 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,371 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_75th to climate bins. 2023-09-15 09:34:00,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,382 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,382 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,385 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate maize yield_75th yield raster to landcover resolution. 2023-09-15 09:34:00,495 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for maize at yield_75th 2023-09-15 09:34:00,508 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,509 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,529 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,529 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,544 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_95th to climate bins. 2023-09-15 09:34:00,553 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,554 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,554 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,554 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,557 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate maize yield_95th yield raster to landcover resolution. 2023-09-15 09:34:00,665 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for maize at yield_95th 2023-09-15 09:34:00,678 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,678 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,697 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,697 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,712 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(635) INFO Calculate observed yield for maize 2023-09-15 09:34:00,727 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,727 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,728 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,731 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(678) INFO Interpolating observed maize raster to landcover. 2023-09-15 09:34:00,862 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,862 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,878 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,879 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,895 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(528) INFO Processing crop grape 2023-09-15 09:34:00,895 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(533) INFO Clipping global climate bin raster to landcover bounding box. 2023-09-15 09:34:00,910 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_25th to climate bins. 2023-09-15 09:34:00,919 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:00,920 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:00,920 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:00,921 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:00,924 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate grape yield_25th yield raster to landcover resolution. 2023-09-15 09:34:01,033 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for grape at yield_25th 2023-09-15 09:34:01,047 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,047 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,066 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,067 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,079 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_50th to climate bins. 2023-09-15 09:34:01,089 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,089 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,090 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,091 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,094 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate grape yield_50th yield raster to landcover resolution. 2023-09-15 09:34:01,207 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for grape at yield_50th 2023-09-15 09:34:01,221 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,221 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,239 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,239 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,252 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_75th to climate bins. 2023-09-15 09:34:01,261 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,261 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,262 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,262 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,265 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate grape yield_75th yield raster to landcover resolution. 2023-09-15 09:34:01,367 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for grape at yield_75th 2023-09-15 09:34:01,380 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,380 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,399 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,399 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,412 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(566) INFO Map yield_95th to climate bins. 2023-09-15 09:34:01,422 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,423 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,424 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,424 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,428 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(595) INFO Interpolate grape yield_95th yield raster to landcover resolution. 2023-09-15 09:34:01,543 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(612) INFO Calculate yield for grape at yield_95th 2023-09-15 09:34:01,558 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,558 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,580 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,580 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,594 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(635) INFO Calculate observed yield for grape 2023-09-15 09:34:01,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,611 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,612 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,612 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,616 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(678) INFO Interpolating observed grape raster to landcover. 2023-09-15 09:34:01,740 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-15 09:34:01,740 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-15 09:34:01,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-15 09:34:01,756 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-15 09:34:01,872 (natcap.invest.crop_production_percentile) crop_production_percentile.tabulate_results(888) INFO Generating report table 2023-09-15 09:34:01,873 (natcap.invest.crop_production_percentile) crop_production_percentile.tabulate_results(916) INFO Calculating production area and summing observed yield. 2023-09-15 09:34:01,899 (natcap.invest.crop_production_percentile) crop_production_percentile.tabulate_results(916) INFO Calculating production area and summing observed yield. 2023-09-15 09:34:01,941 (natcap.invest.crop_production_percentile) crop_production_percentile.execute(730) INFO aggregating result over query polygon 2023-09-15 09:34:02,041 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2037) INFO starting reprojection 2023-09-15 09:34:02,042 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Value 1581427853.98000002 of field Shape_Area of feature 0 not successfully written. Possibly due to too larger number with respect to field width 2023-09-15 09:34:02,042 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2077) INFO reprojection 100.0% complete on aggregate_vector_CP2018.shp 2023-09-15 09:34:02,043 (natcap.invest.crop_production_percentile) crop_production_percentile.aggregate_to_polygons(1033) INFO Calculating zonal stats for maize yield_25th 2023-09-15 09:34:02,090 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-09-15 09:34:02,090 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-09-15 09:34:02,096 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-09-15 09:34:02,096 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-09-15 09:34:02,096 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list