2023-09-24 21:13:36,103 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [D:\窟野河\INVEST\2020_result\InVEST-natcap.invest.habitat_quality-log-2023-09-24--21_13_36.txt] 2023-09-24 21:13:36,103 (__main__) cli.main(458) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.habitat_quality 3.13.0: access_vector_path half_saturation_constant 0.05 lulc_bas_path lulc_cur_path D:\窟野河\INVEST\2020\2020landuse2.tif lulc_fut_path n_workers -1 results_suffix sensitivity_table_path D:\窟野河\INVEST\2020\sensitivity3.csv threats_table_path D:\窟野河\INVEST\2020\threat4.csv workspace_dir D:\窟野河\INVEST\2020_result 2023-09-24 21:13:36,104 (natcap.invest.habitat_quality) habitat_quality.execute(362) INFO Starting execute of Habitat Quality model. 2023-09-24 21:13:36,104 (natcap.invest.habitat_quality) habitat_quality.execute(368) INFO Creating workspace 2023-09-24 21:13:36,116 (natcap.invest.habitat_quality) habitat_quality.execute(380) INFO Checking Threat and Sensitivity tables for compliance 2023-09-24 21:13:36,131 (natcap.invest.utils) utils.build_lookup_from_csv(666) WARNING Empty or NaN values were found in the table: D:\窟野河\INVEST\2020\sensitivity3.csv. 2023-09-24 21:13:36,132 (natcap.invest.habitat_quality) habitat_quality.execute(404) INFO Validate threat rasters and collect unique LULC codes 2023-09-24 21:13:36,746 (natcap.invest.habitat_quality) habitat_quality.execute(459) INFO Checking threat raster values are valid ( 0 <= x <= 1 ). 2023-09-24 21:13:36,746 (natcap.invest.habitat_quality) habitat_quality.execute(470) INFO Aligning, resizing, and reprojecting raster inputs to that of the current land cover. 2023-09-24 21:13:37,481 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 5 aligned: 2020landuse2_aligned.tif 2023-09-24 21:13:38,068 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 5 aligned: cultivatedland_c_aligned.tif 2023-09-24 21:13:38,638 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 5 aligned: village_c_aligned.tif 2023-09-24 21:13:39,195 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 4 of 5 aligned: urban_c_aligned.tif 2023-09-24 21:13:39,744 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 5 of 5 aligned: industry_c_aligned.tif 2023-09-24 21:13:39,744 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 5 rasters. 2023-09-24 21:13:39,750 (natcap.invest.habitat_quality) habitat_quality.execute(529) INFO Starting habitat_quality biophysical calculations 2023-09-24 21:13:40,030 (natcap.invest.habitat_quality) habitat_quality.execute(587) INFO Calculating habitat quality for landuse: D:\窟野河\INVEST\2020_result\intermediate\2020landuse2_aligned.tif 2023-09-24 21:13:40,045 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-24 21:13:40,046 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-24 21:13:40,049 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(571) ERROR exception encountered in raster_calculator Traceback (most recent call last): File "pygeoprocessing\geoprocessing.py", line 517, in raster_calculator File "pygeoprocessing\geoprocessing.py", line 2179, in _map_dataset_to_value_op pygeoprocessing.geoprocessing.ReclassificationMissingValuesError: The following 1 raster values [-128] from "D:\窟野河\INVEST\2020_result\intermediate\2020landuse2_aligned.tif" do not have corresponding entries in the value map: {11: 0.6, 12: 0.5, 21: 0.9, 22: 0.8, 23: 0.5, 24: 0.3, 31: 0.8, 32: 0.6, 33: 0.4, 41: 0.8, 42: 0.9, 43: 0.7, 46: 0.6, 51: 0.0, 52: 0.0, 53: 0.0, 61: 0.1, 63: 0.1, 64: 0.9, 65: 0.0, 0: 0.0}. 2023-09-24 21:13:40,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(586) INFO Waiting for raster stats worker result. 2023-09-24 21:13:40,052 (pygeoprocessing.geoprocessing_core) threading.run(953) WARNING No valid pixels were received, sending None. 2023-09-24 21:13:40,053 (taskgraph.Task) Task.add_task(706) ERROR Something went wrong when adding task habitat_raster_c (6), terminating taskgraph. Traceback (most recent call last): File "natcap\invest\utils.py", line 1047, in reclassify_raster File "pygeoprocessing\geoprocessing.py", line 2186, in reclassify_raster File "pygeoprocessing\geoprocessing.py", line 517, in raster_calculator File "pygeoprocessing\geoprocessing.py", line 2179, in _map_dataset_to_value_op pygeoprocessing.geoprocessing.ReclassificationMissingValuesError: The following 1 raster values [-128] from "D:\窟野河\INVEST\2020_result\intermediate\2020landuse2_aligned.tif" do not have corresponding entries in the value map: {11: 0.6, 12: 0.5, 21: 0.9, 22: 0.8, 23: 0.5, 24: 0.3, 31: 0.8, 32: 0.6, 33: 0.4, 41: 0.8, 42: 0.9, 43: 0.7, 46: 0.6, 51: 0.0, 52: 0.0, 53: 0.0, 61: 0.1, 63: 0.1, 64: 0.9, 65: 0.0, 0: 0.0}. During handling of the above exception, another exception occurred: Traceback (most recent call last): File "taskgraph\Task.py", line 674, in add_task File "taskgraph\Task.py", line 1093, in _call File "natcap\invest\utils.py", line 1057, in reclassify_raster ValueError: Values in the LULC_c raster were found that are not represented under the 'lucode' column of the Sensitivity table. The missing values found in the LULC_c raster but not the table are: [-128]. 2023-09-24 21:13:40,054 (natcap.invest.utils) utils.prepare_workspace(166) ERROR Exception while executing natcap.invest.habitat_quality Traceback (most recent call last): File "natcap\invest\utils.py", line 1047, in reclassify_raster File "pygeoprocessing\geoprocessing.py", line 2186, in reclassify_raster File "pygeoprocessing\geoprocessing.py", line 517, in raster_calculator File "pygeoprocessing\geoprocessing.py", line 2179, in _map_dataset_to_value_op pygeoprocessing.geoprocessing.ReclassificationMissingValuesError: The following 1 raster values [-128] from "D:\窟野河\INVEST\2020_result\intermediate\2020landuse2_aligned.tif" do not have corresponding entries in the value map: {11: 0.6, 12: 0.5, 21: 0.9, 22: 0.8, 23: 0.5, 24: 0.3, 31: 0.8, 32: 0.6, 33: 0.4, 41: 0.8, 42: 0.9, 43: 0.7, 46: 0.6, 51: 0.0, 52: 0.0, 53: 0.0, 61: 0.1, 63: 0.1, 64: 0.9, 65: 0.0, 0: 0.0}. During handling of the above exception, another exception occurred: Traceback (most recent call last): File "natcap\invest\utils.py", line 164, in prepare_workspace File "invest\cli.py", line 469, in main File "natcap\invest\habitat_quality.py", line 600, in execute File "taskgraph\Task.py", line 674, in add_task File "taskgraph\Task.py", line 1093, in _call File "natcap\invest\utils.py", line 1057, in reclassify_raster ValueError: Values in the LULC_c raster were found that are not represented under the 'lucode' column of the Sensitivity table. The missing values found in the LULC_c raster but not the table are: [-128]. 2023-09-24 21:13:40,055 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 3.95s 2023-09-24 21:13:40,055 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished