2022-09-06 14:51:58,913 (natcap.invest.utils) utils.prepare_workspace(158) INFO Writing log messages to F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\InVEST-Urban-Cooling-log-2022-09-06--14_51_58.txt 2022-09-06 14:51:58,913 (natcap.invest.ui.model) model._logged_target(1644) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.urban_cooling_model 3.9.2: aoi_vector_path F:/PhD/UrbanCooling3.0/data/PRG_core/core_projected.shp avg_rel_humidity 68.17 biophysical_table_path F:/PhD/UrbanCooling3.0/data/biophysical_tables/biophysical_table_V_building_intensity_adjusted.csv cc_method factors cc_weight_albedo 0.2 cc_weight_eti 0.2 cc_weight_shade 0.6 do_energy_valuation False do_productivity_valuation True green_area_cooling_distance 450 lulc_raster_path F:/PhD/UrbanCooling3.0/data/LULC_improved/veg/LULC.tif n_workers -1 ref_eto_raster_path F:/PhD/UrbanCooling3.0/data/Evapotranspiraiton/PRG_summer_mean/smr_et_prg.tif results_suffix t_air_average_radius 500 t_ref 15.8 uhi_max 4.0 workspace_dir F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder 2022-09-06 14:51:58,913 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(305) INFO Starting Urban Cooling Model 2022-09-06 14:52:04,159 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 4.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:09,219 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 8.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:15,397 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 13.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:21,553 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 18.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:27,662 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 23.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:33,790 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 28.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:39,899 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 33.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:46,007 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 38.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:52,128 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 43.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:52:58,238 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 48.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:04,515 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 53.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:10,627 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 58.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:16,741 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 63.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:22,869 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 68.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:28,978 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 73.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:35,125 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 78.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:41,234 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 83.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:47,348 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 88.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:53,473 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 93.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:53:59,596 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 98.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:54:02,036 (pygeoprocessing.geoprocessing) geoprocessing.logger_callback(3320) INFO Warp 100.0% complete F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\lulc.tif 2022-09-06 14:54:02,218 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 2 aligned: lulc.tif 2022-09-06 14:54:06,554 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 2 of 2 aligned: ref_eto.tif 2022-09-06 14:54:06,554 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 2 rasters. 2022-09-06 14:54:06,610 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:54:06,612 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:54:06,612 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:54:07,781 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:54:07,782 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:54:08,240 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:54:08,241 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:54:08,242 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:54:09,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:54:09,381 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:54:09,734 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:54:09,735 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:54:09,735 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:54:10,875 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:54:10,875 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:54:11,263 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:54:11,264 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:54:11,264 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:54:12,416 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:54:12,417 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:54:12,857 (natcap.invest.urban_cooling_model) urban_cooling_model.convolve_2d_by_exponential(1399) INFO Starting a convolution over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\green_area.tif with a decay distance of 45 2022-09-06 14:54:12,932 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2678) INFO starting convolve 2022-09-06 14:54:12,935 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2696) DEBUG start fill work queue thread 2022-09-06 14:54:12,936 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2700) DEBUG fill work queue 2022-09-06 14:54:12,936 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2714) DEBUG start worker thread 2022-09-06 14:54:12,936 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2726) INFO 1344 sent to workers, wait for worker results 2022-09-06 14:54:17,942 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 15.0% complete on cc_park.tif 2022-09-06 14:54:22,943 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 29.7% complete on cc_park.tif 2022-09-06 14:54:27,955 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 44.3% complete on cc_park.tif 2022-09-06 14:54:32,968 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 58.9% complete on cc_park.tif 2022-09-06 14:54:37,984 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 72.4% complete on cc_park.tif 2022-09-06 14:54:43,028 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 86.1% complete on cc_park.tif 2022-09-06 14:54:47,462 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2705) DEBUG work queue full 2022-09-06 14:54:47,725 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2792) INFO convolution worker 100.0% complete on cc_park.tif 2022-09-06 14:54:48,971 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2800) INFO need to normalize result so nodata values are not included 2022-09-06 14:54:49,871 (pygeoprocessing.geoprocessing) geoprocessing.(2829) INFO convolution nodata normalizer approximately 0.3% complete on cc_park.tif 2022-09-06 14:54:51,827 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2839) INFO convolution nodata normalize 100.0% complete on cc_park.tif 2022-09-06 14:54:54,273 (natcap.invest.urban_cooling_model) urban_cooling_model.flat_disk_kernel(1220) INFO Creating a disk kernel of distance 45 at F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\area_kernel.tif 2022-09-06 14:54:54,337 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2678) INFO starting convolve 2022-09-06 14:54:54,339 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2696) DEBUG start fill work queue thread 2022-09-06 14:54:54,340 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2700) DEBUG fill work queue 2022-09-06 14:54:54,340 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2714) DEBUG start worker thread 2022-09-06 14:54:54,341 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2726) INFO 336 sent to workers, wait for worker results 2022-09-06 14:54:58,761 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2705) DEBUG work queue full 2022-09-06 14:54:58,902 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2792) INFO convolution worker 100.0% complete on green_area_sum.tif 2022-09-06 14:54:59,844 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2800) INFO need to normalize result so nodata values are not included 2022-09-06 14:55:00,443 (pygeoprocessing.geoprocessing) geoprocessing.(2829) INFO convolution nodata normalizer approximately 0.3% complete on green_area_sum.tif 2022-09-06 14:55:02,021 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2839) INFO convolution nodata normalize 100.0% complete on green_area_sum.tif 2022-09-06 14:55:03,594 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(458) INFO Calculating Cooling Coefficient from factors 2022-09-06 14:55:04,281 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:55:04,283 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:55:04,283 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:55:05,318 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:55:05,318 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:55:06,495 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:55:06,496 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:55:06,496 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:55:07,846 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:55:07,846 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:55:08,921 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:55:08,922 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:55:08,922 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:55:11,956 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:55:11,956 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:55:13,026 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:55:13,027 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:55:13,027 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:55:13,854 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:55:13,855 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:55:15,096 (natcap.invest.urban_cooling_model) urban_cooling_model.convolve_2d_by_exponential(1399) INFO Starting a convolution over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\T_air_nomix.tif with a decay distance of 50 2022-09-06 14:55:15,178 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2678) INFO starting convolve 2022-09-06 14:55:15,180 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2696) DEBUG start fill work queue thread 2022-09-06 14:55:15,181 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2700) DEBUG fill work queue 2022-09-06 14:55:15,181 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2714) DEBUG start worker thread 2022-09-06 14:55:15,182 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2726) INFO 1344 sent to workers, wait for worker results 2022-09-06 14:55:20,179 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 13.6% complete on T_air.tif 2022-09-06 14:55:25,207 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 27.3% complete on T_air.tif 2022-09-06 14:55:30,233 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 40.3% complete on T_air.tif 2022-09-06 14:55:35,239 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 53.2% complete on T_air.tif 2022-09-06 14:55:40,250 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 65.2% complete on T_air.tif 2022-09-06 14:55:45,287 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 78.1% complete on T_air.tif 2022-09-06 14:55:50,311 (pygeoprocessing.geoprocessing) geoprocessing.(2786) INFO convolution worker approximately 91.3% complete on T_air.tif 2022-09-06 14:55:53,230 (pygeoprocessing.geoprocessing) geoprocessing._fill_work_queue(2705) DEBUG work queue full 2022-09-06 14:55:53,533 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2792) INFO convolution worker 100.0% complete on T_air.tif 2022-09-06 14:55:54,768 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2800) INFO need to normalize result so nodata values are not included 2022-09-06 14:55:55,687 (pygeoprocessing.geoprocessing) geoprocessing.(2829) INFO convolution nodata normalizer approximately 0.3% complete on T_air.tif 2022-09-06 14:55:57,752 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2839) INFO convolution nodata normalize 100.0% complete on T_air.tif 2022-09-06 14:56:00,158 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2022-09-06 14:56:00,159 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:00,187 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1067) INFO Taking zonal statistics of F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\reprojected_aoi.shp over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\cc.tif 2022-09-06 14:56:02,103 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-09-06 14:56:02,103 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-09-06 14:56:02,123 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-09-06 14:56:02,123 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-09-06 14:56:02,124 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-09-06 14:56:02,131 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-09-06 14:56:02,132 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:02,132 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-09-06 14:56:02,133 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:02,133 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:02,146 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-09-06 14:56:02,147 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on reprojected_aoi.shp 2022-09-06 14:56:02,154 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:02,484 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:03,292 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-09-06 14:56:03,293 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.5834622e+06 1.0000000e+01 0.0000000e+00 6.4770605e+06 0.0000000e+00 -1.0000000e+01] for ('F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\\intermediate\\cc.tif', 1) 2022-09-06 14:56:03,293 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-09-06 14:56:03,293 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.shp 2022-09-06 14:56:03,330 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1067) INFO Taking zonal statistics of F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\reprojected_aoi.shp over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\T_air.tif 2022-09-06 14:56:07,096 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-09-06 14:56:07,097 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-09-06 14:56:07,102 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-09-06 14:56:07,103 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-09-06 14:56:07,103 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-09-06 14:56:07,108 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-09-06 14:56:07,109 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:07,110 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-09-06 14:56:07,110 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:07,110 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:07,129 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-09-06 14:56:07,129 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on reprojected_aoi.shp 2022-09-06 14:56:07,178 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:07,511 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:09,063 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-09-06 14:56:09,064 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.5834622e+06 1.0000000e+01 0.0000000e+00 6.4770605e+06 0.0000000e+00 -1.0000000e+01] for ('F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\\intermediate\\T_air.tif', 1) 2022-09-06 14:56:09,064 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-09-06 14:56:09,065 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.shp 2022-09-06 14:56:09,116 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(598) INFO Starting work productivity valuation 2022-09-06 14:56:09,118 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_wbgt(1180) INFO Calculating WBGT 2022-09-06 14:56:09,128 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:56:09,129 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:56:09,129 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:56:10,813 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:56:10,814 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:56:12,020 (natcap.invest.urban_cooling_model) urban_cooling_model.map_work_loss(1328) INFO Calculating work loss using thresholds: [31.5, 32.0, 32.5] 2022-09-06 14:56:12,032 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:56:12,034 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:56:12,034 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:56:12,823 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:56:12,823 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:56:12,924 (natcap.invest.urban_cooling_model) urban_cooling_model.map_work_loss(1328) INFO Calculating work loss using thresholds: [27.5, 29.5, 31.5] 2022-09-06 14:56:12,935 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2022-09-06 14:56:12,936 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 6148 2022-09-06 14:56:12,936 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2022-09-06 14:56:13,717 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2022-09-06 14:56:13,717 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2022-09-06 14:56:13,811 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1067) INFO Taking zonal statistics of F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\reprojected_aoi.shp over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\wbgt.tif 2022-09-06 14:56:15,800 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-09-06 14:56:15,801 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-09-06 14:56:15,811 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-09-06 14:56:15,812 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-09-06 14:56:15,812 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-09-06 14:56:15,818 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-09-06 14:56:15,819 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:15,820 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-09-06 14:56:15,820 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:15,820 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:15,833 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-09-06 14:56:15,834 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on reprojected_aoi.shp 2022-09-06 14:56:15,843 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:16,171 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:17,044 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-09-06 14:56:17,045 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.5834622e+06 1.0000000e+01 0.0000000e+00 6.4770605e+06 0.0000000e+00 -1.0000000e+01] for ('F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\\intermediate\\wbgt.tif', 1) 2022-09-06 14:56:17,045 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-09-06 14:56:17,045 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.shp 2022-09-06 14:56:17,102 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1067) INFO Taking zonal statistics of F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\reprojected_aoi.shp over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\light_work_loss_percent.tif 2022-09-06 14:56:17,602 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-09-06 14:56:17,602 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-09-06 14:56:17,611 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-09-06 14:56:17,612 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-09-06 14:56:17,613 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-09-06 14:56:17,618 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-09-06 14:56:17,619 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:17,619 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-09-06 14:56:17,619 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:17,620 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:17,632 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-09-06 14:56:17,633 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on reprojected_aoi.shp 2022-09-06 14:56:17,663 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:17,999 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:18,427 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-09-06 14:56:18,427 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.5834622e+06 1.0000000e+01 0.0000000e+00 6.4770605e+06 0.0000000e+00 -1.0000000e+01] for ('F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\\intermediate\\light_work_loss_percent.tif', 1) 2022-09-06 14:56:18,428 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-09-06 14:56:18,428 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.shp 2022-09-06 14:56:18,455 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1067) INFO Taking zonal statistics of F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\reprojected_aoi.shp over F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\intermediate\heavy_work_loss_percent.tif 2022-09-06 14:56:18,949 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(824) INFO 1 of 1 aligned: clipped_raster.tif 2022-09-06 14:56:18,950 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(828) INFO aligned all 1 rasters. 2022-09-06 14:56:18,958 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1275) INFO build a lookup of aggregate field value to FID 2022-09-06 14:56:18,958 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1282) INFO creating disjoint polygon set 2022-09-06 14:56:18,959 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2257) INFO build shapely polygon list 2022-09-06 14:56:18,964 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2280) INFO build shapely rtree index 2022-09-06 14:56:18,965 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2293) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:18,965 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2297) INFO build poly intersection lookup 2022-09-06 14:56:18,966 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2319) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:18,966 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2353) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2022-09-06 14:56:18,978 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1303) INFO processing 1 disjoint polygon sets 2022-09-06 14:56:18,979 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1337) INFO disjoint polygon set 1 of 1 100.0% processed on reprojected_aoi.shp 2022-09-06 14:56:18,988 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1345) INFO rasterizing disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:19,321 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1363) INFO summarizing rasterized disjoint polygon set 1 of 1 reprojected_aoi.shp 2022-09-06 14:56:19,748 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1407) DEBUG unset_fids: 0 of 1 2022-09-06 14:56:19,748 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1412) DEBUG gt [ 1.5834622e+06 1.0000000e+01 0.0000000e+00 6.4770605e+06 0.0000000e+00 -1.0000000e+01] for ('F:/PhD/UrbanCooling3.0/InVEST_runs/13_07_2018/hourly/veg/no_water/0/New folder\\intermediate\\heavy_work_loss_percent.tif', 1) 2022-09-06 14:56:19,748 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1491) DEBUG remaining unset_fids: 0 of 1 2022-09-06 14:56:19,748 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1499) INFO all done processing polygon sets for reprojected_aoi.shp 2022-09-06 14:56:19,777 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(779) INFO Calculate UHI summary results uhi_results.shp 2022-09-06 14:56:19,778 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(783) INFO Loading t_air_stats 2022-09-06 14:56:19,779 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(786) INFO Loading cc_stats 2022-09-06 14:56:19,779 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(792) INFO Loading wbgt_stats 2022-09-06 14:56:19,779 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(798) INFO Loading light_loss_stats 2022-09-06 14:56:19,779 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(805) INFO Loading heavy_loss_stats 2022-09-06 14:56:19,780 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(818) INFO Creating uhi_results.shp 2022-09-06 14:56:19,827 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(883) DEBUG Average light loss: 0.0 2022-09-06 14:56:19,827 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(891) DEBUG Average heavy loss: 0.0 2022-09-06 14:56:19,856 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(735) INFO Urban Cooling Model complete. 2022-09-06 14:56:19,858 (natcap.invest.utils) utils.prepare_workspace(166) INFO Elapsed time: 4m 20.94999999999999s 2022-09-06 14:56:19,858 (natcap.invest.utils) utils.prepare_workspace(169) INFO Execution finished