2023-09-23 20:32:58,112 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\HQ_major_final\InVEST-natcap.invest.habitat_quality-log-2023-09-23--20_32_58.txt] 2023-09-23 20:32:58,112 (__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.5 lulc_bas_path lulc_cur_path D:\postdoc research\TM\1 Eshetu\Abay_upper_watershed\HabitatQ\lulc_final.tif lulc_fut_path n_workers -1 results_suffix final sensitivity_table_path D:\postdoc research\TM\1 Eshetu\Abay_upper_watershed\HabitatQ\sensitivity_ubnb_final.csv threats_table_path D:\postdoc research\TM\1 Eshetu\Abay_upper_watershed\HabitatQ\threats_final.csv workspace_dir C:\HQ_major_final 2023-09-23 20:32:58,113 (natcap.invest.habitat_quality) habitat_quality.execute(362) INFO Starting execute of Habitat Quality model. 2023-09-23 20:32:58,113 (natcap.invest.habitat_quality) habitat_quality.execute(368) INFO Creating workspace 2023-09-23 20:32:58,150 (natcap.invest.habitat_quality) habitat_quality.execute(380) INFO Checking Threat and Sensitivity tables for compliance 2023-09-23 20:32:58,164 (natcap.invest.habitat_quality) habitat_quality.execute(404) INFO Validate threat rasters and collect unique LULC codes 2023-09-23 20:33:04,761 (natcap.invest.habitat_quality) habitat_quality.execute(459) INFO Checking threat raster values are valid ( 0 <= x <= 1 ). 2023-09-23 20:33:04,761 (natcap.invest.habitat_quality) habitat_quality.execute(470) INFO Aligning, resizing, and reprojecting raster inputs to that of the current land cover. 2023-09-23 20:33:09,977 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 50.0% complete C:\HQ_major_final\intermediate\lulc_final_aligned_final.tif 2023-09-23 20:33:14,542 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 100.0% complete C:\HQ_major_final\intermediate\lulc_final_aligned_final.tif 2023-09-23 20:33:15,855 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 4 aligned: lulc_final_aligned_final.tif 2023-09-23 20:33:20,893 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 80.0% complete C:\HQ_major_final\intermediate\bare_f_aligned_final.tif 2023-09-23 20:33:22,057 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 100.0% complete C:\HQ_major_final\intermediate\bare_f_aligned_final.tif 2023-09-23 20:33:22,694 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 4 aligned: bare_f_aligned_final.tif 2023-09-23 20:33:27,762 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 68.0% complete C:\HQ_major_final\intermediate\construction_f_aligned_final.tif 2023-09-23 20:33:30,312 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 100.0% complete C:\HQ_major_final\intermediate\construction_f_aligned_final.tif 2023-09-23 20:33:30,961 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 3 of 4 aligned: construction_f_aligned_final.tif 2023-09-23 20:33:36,074 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 75.0% complete C:\HQ_major_final\intermediate\crops_f_aligned_final.tif 2023-09-23 20:33:37,492 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3601) INFO Warp 100.0% complete C:\HQ_major_final\intermediate\crops_f_aligned_final.tif 2023-09-23 20:33:38,558 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 4 of 4 aligned: crops_f_aligned_final.tif 2023-09-23 20:33:38,559 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 4 rasters. 2023-09-23 20:33:38,601 (natcap.invest.habitat_quality) habitat_quality.execute(529) INFO Starting habitat_quality biophysical calculations 2023-09-23 20:33:42,550 (natcap.invest.habitat_quality) habitat_quality.execute(587) INFO Calculating habitat quality for landuse: C:\HQ_major_final\intermediate\lulc_final_aligned_final.tif 2023-09-23 20:33:42,572 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:33:42,574 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:33:47,562 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO habitat_c_final.tif 42.9% complete 2023-09-23 20:33:52,564 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO habitat_c_final.tif 67.8% complete 2023-09-23 20:33:56,575 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:33:56,575 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:34:01,800 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:34:17,911 (pygeoprocessing.geoprocessing_core) geoprocessing.distance_transform_edt(2738) INFO Distance Transform Phase 2 2023-09-23 20:35:27,889 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:35:27,890 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:35:32,875 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_bare_f_c_final.tif 32.5% complete 2023-09-23 20:35:37,878 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_bare_f_c_final.tif 59.5% complete 2023-09-23 20:35:42,879 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_bare_f_c_final.tif 89.4% complete 2023-09-23 20:35:44,553 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:35:44,553 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:35:45,115 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:35:45,116 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:35:50,106 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_bare_f_c_final.tif 42.1% complete 2023-09-23 20:35:55,107 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_bare_f_c_final.tif 65.5% complete 2023-09-23 20:36:00,035 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:36:00,035 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:36:06,007 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO region_mask_rasteraqbpufmy.tif 90.5% complete 2023-09-23 20:36:06,720 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:36:22,988 (pygeoprocessing.geoprocessing_core) geoprocessing.distance_transform_edt(2738) INFO Distance Transform Phase 2 2023-09-23 20:37:40,560 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:37:40,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:37:45,522 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_construction_f_c_final.tif 28.0% complete 2023-09-23 20:37:50,523 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_construction_f_c_final.tif 52.3% complete 2023-09-23 20:37:55,524 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_construction_f_c_final.tif 71.8% complete 2023-09-23 20:38:00,526 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_exponential_construction_f_c_final.tif 92.2% complete 2023-09-23 20:38:02,051 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:38:02,052 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:38:02,751 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:38:02,755 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:38:07,739 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_construction_f_c_final.tif 41.2% complete 2023-09-23 20:38:12,742 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_construction_f_c_final.tif 63.1% complete 2023-09-23 20:38:17,743 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_construction_f_c_final.tif 92.3% complete 2023-09-23 20:38:18,638 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:38:18,638 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:38:24,587 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO region_mask_raster_jrdk7qm.tif 83.3% complete 2023-09-23 20:38:26,031 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:38:41,172 (pygeoprocessing.geoprocessing_core) geoprocessing.distance_transform_edt(2738) INFO Distance Transform Phase 2 2023-09-23 20:39:57,801 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:39:57,803 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:40:02,786 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_linear_crops_f_c_final.tif 34.3% complete 2023-09-23 20:40:07,786 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_linear_crops_f_c_final.tif 58.4% complete 2023-09-23 20:40:12,789 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO filtered_linear_crops_f_c_final.tif 79.1% complete 2023-09-23 20:40:17,260 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:40:17,260 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:40:18,102 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:40:18,102 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:40:23,039 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_crops_f_c_final.tif 40.5% complete 2023-09-23 20:40:28,041 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_crops_f_c_final.tif 62.3% complete 2023-09-23 20:40:33,042 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO sens_crops_f_c_final.tif 97.4% complete 2023-09-23 20:40:33,387 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:40:33,388 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:40:34,300 (natcap.invest.habitat_quality) habitat_quality.execute(717) INFO Starting raster calculation on total degradation 2023-09-23 20:40:34,380 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:40:34,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:40:39,368 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 13.5% complete 2023-09-23 20:40:44,368 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 24.4% complete 2023-09-23 20:40:49,373 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 34.2% complete 2023-09-23 20:40:54,374 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 43.5% complete 2023-09-23 20:40:59,377 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 52.3% complete 2023-09-23 20:41:04,388 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 60.8% complete 2023-09-23 20:41:09,390 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 69.2% complete 2023-09-23 20:41:14,392 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 79.6% complete 2023-09-23 20:41:19,396 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO deg_sum_c_final.tif 93.6% complete 2023-09-23 20:41:21,731 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:41:21,731 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:41:21,934 (natcap.invest.habitat_quality) habitat_quality.execute(735) INFO Starting raster calculation on quality 2023-09-23 20:41:21,962 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-09-23 20:41:21,963 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-09-23 20:41:26,951 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO quality_c_final.tif 29.1% complete 2023-09-23 20:41:31,954 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO quality_c_final.tif 44.2% complete 2023-09-23 20:41:36,956 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO quality_c_final.tif 57.8% complete 2023-09-23 20:41:41,960 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO quality_c_final.tif 70.7% complete 2023-09-23 20:41:46,960 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(550) INFO quality_c_final.tif 91.8% complete 2023-09-23 20:41:48,338 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-09-23 20:41:48,338 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-09-23 20:41:48,872 (natcap.invest.habitat_quality) habitat_quality.execute(748) INFO Baseline not provided to compute Rarity 2023-09-23 20:41:48,872 (natcap.invest.habitat_quality) habitat_quality.execute(775) INFO Habitat Quality Model complete. 2023-09-23 20:41:48,893 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 8m 50.77999999999997s 2023-09-23 20:41:48,894 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished