2023-07-13 12:19:00,157 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\InVEST-natcap.invest.urban_cooling_model-log-2023-07-13--12_19_00.txt] 2023-07-13 12:19:00,157 (__main__) cli.main(458) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.urban_cooling_model 3.13.0: aoi_vector_path C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Greater_Springfield.shp avg_rel_humidity biophysical_table_path C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Biophysical_Table_Springfield.csv building_vector_path cc_method factors cc_weight_albedo cc_weight_eti cc_weight_shade do_energy_valuation False do_productivity_valuation False energy_consumption_table_path green_area_cooling_distance 100 lulc_raster_path C:\Users\tomas\Documents\ArcGIS\Projects\Springfield_cooling\Clipped_LULC_P2R5.tif n_workers -1 ref_eto_raster_path C:\Users\tomas\Documents\ArcGIS\Projects\Springfield_cooling\ET_500m.tif results_suffix t_air_average_radius 2000 t_ref 0.59 uhi_max 2.09 workspace_dir C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results 2023-07-13 12:19:00,159 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(408) INFO Starting Urban Cooling Model 2023-07-13 12:19:00,408 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 2 aligned: lulc.tif 2023-07-13 12:19:00,610 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 2 of 2 aligned: ref_eto.tif 2023-07-13 12:19:00,611 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 2 rasters. 2023-07-13 12:19:00,641 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:00,656 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:00,692 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:00,692 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:00,722 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:00,743 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:00,746 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:00,746 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:00,775 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:00,819 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:00,824 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:00,824 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:00,849 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:00,919 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:00,923 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:00,923 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:00,932 (natcap.invest.urban_cooling_model) urban_cooling_model.convolve_2d_by_exponential(1494) INFO Starting a convolution over C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\green_area.tif with a decay distance of 2 2023-07-13 12:19:00,968 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2983) INFO starting convolve 2023-07-13 12:19:01,004 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3031) INFO 2 sent to workers, wait for worker results 2023-07-13 12:19:01,094 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3095) INFO convolution worker 100.0% complete on cc_park.tif 2023-07-13 12:19:01,105 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3103) INFO need to normalize result so nodata values are not included 2023-07-13 12:19:01,126 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3139) INFO convolution nodata normalize 100.0% complete on cc_park.tif 2023-07-13 12:19:01,168 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2983) INFO starting convolve 2023-07-13 12:19:01,302 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3031) INFO 2 sent to workers, wait for worker results 2023-07-13 12:19:01,523 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3095) INFO convolution worker 100.0% complete on green_area_sum.tif 2023-07-13 12:19:01,530 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3103) INFO need to normalize result so nodata values are not included 2023-07-13 12:19:01,541 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3139) INFO convolution nodata normalize 100.0% complete on green_area_sum.tif 2023-07-13 12:19:01,555 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(560) INFO Calculating Cooling Coefficient from factors 2023-07-13 12:19:01,589 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:01,600 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:01,607 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:01,607 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:01,640 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:01,641 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:01,645 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:01,645 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:01,685 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:01,688 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:01,697 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:01,697 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:01,738 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-07-13 12:19:01,749 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-07-13 12:19:01,753 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-07-13 12:19:01,754 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-07-13 12:19:01,766 (natcap.invest.urban_cooling_model) urban_cooling_model.convolve_2d_by_exponential(1494) INFO Starting a convolution over C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\T_air_nomix.tif with a decay distance of 43 2023-07-13 12:19:01,853 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(2983) INFO starting convolve 2023-07-13 12:19:01,872 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3031) INFO 8 sent to workers, wait for worker results 2023-07-13 12:19:04,175 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3095) INFO convolution worker 100.0% complete on T_air.tif 2023-07-13 12:19:04,189 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3103) INFO need to normalize result so nodata values are not included 2023-07-13 12:19:04,209 (pygeoprocessing.geoprocessing) geoprocessing.convolve_2d(3139) INFO convolution nodata normalize 100.0% complete on T_air.tif 2023-07-13 12:19:04,249 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1161) INFO Taking zonal statistics of C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\reprojected_aoi.shp over C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\cc.tif 2023-07-13 12:19:04,324 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-07-13 12:19:04,324 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-07-13 12:19:04,341 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-07-13 12:19:04,342 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-07-13 12:19:04,342 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-07-13 12:19:04,350 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-07-13 12:19:04,351 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,352 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-07-13 12:19:04,353 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,354 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,362 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 3 disjoint polygon sets 2023-07-13 12:19:04,362 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,362 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,368 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,373 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,374 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,376 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,377 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,378 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,380 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,382 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for reprojected_aoi.shp 2023-07-13 12:19:04,395 (natcap.invest.urban_cooling_model) urban_cooling_model.pickle_zonal_stats(1161) INFO Taking zonal statistics of C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\reprojected_aoi.shp over C:\Users\tomas\Documents\ArcGIS Pro 3.1\Springfield_Cooling\Results\intermediate\T_air.tif 2023-07-13 12:19:04,423 (osgeo) utils._log_gdal_errors(97) ERROR [errno 2] d:\bld\gdal-split_1654513788754\work\gcore\gdalrasterblock.cpp, 1125: cannot allocate 524288 bytes 2023-07-13 12:19:04,423 (osgeo) utils._log_gdal_errors(97) ERROR [errno 1] GetBlockRef failed at X block offset 0, Y block offset 0: d:\bld\gdal-split_1654513788754\work\gcore\gdalrasterblock.cpp, 1125: cannot allocate 524288 bytes 2023-07-13 12:19:04,426 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(996) INFO 1 of 1 aligned: clipped_raster.tif 2023-07-13 12:19:04,426 (pygeoprocessing.geoprocessing) geoprocessing.align_and_resize_raster_stack(1000) INFO aligned all 1 rasters. 2023-07-13 12:19:04,442 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1524) INFO build a lookup of aggregate field value to FID 2023-07-13 12:19:04,442 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1531) INFO creating disjoint polygon set 2023-07-13 12:19:04,443 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2554) INFO build shapely polygon list 2023-07-13 12:19:04,450 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2580) INFO build shapely rtree index 2023-07-13 12:19:04,450 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2593) INFO poly feature lookup 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,451 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2597) INFO build poly intersection lookup 2023-07-13 12:19:04,451 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2623) INFO poly intersection feature lookup 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,452 (pygeoprocessing.geoprocessing) geoprocessing.calculate_disjoint_polygon_set(2655) INFO maximal subset build 100.0% complete on reprojected_aoi.shp 2023-07-13 12:19:04,460 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1551) INFO processing 3 disjoint polygon sets 2023-07-13 12:19:04,461 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 1 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,461 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 1 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,466 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 1 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,471 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 2 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,471 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 2 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,474 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 2 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,475 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1581) INFO disjoint polygon set 3 of 3 100.0% processed on reprojected_aoi.shp 2023-07-13 12:19:04,475 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1589) INFO rasterizing disjoint polygon set 3 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,478 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1607) INFO summarizing rasterized disjoint polygon set 3 of 3 reprojected_aoi.shp 2023-07-13 12:19:04,479 (pygeoprocessing.geoprocessing) geoprocessing.zonal_statistics(1756) INFO all done processing polygon sets for reprojected_aoi.shp 2023-07-13 12:19:04,492 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(881) INFO Calculate UHI summary results for uhi_results.shp 2023-07-13 12:19:04,493 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(884) INFO Loading t_air_stats 2023-07-13 12:19:04,493 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(887) INFO Loading cc_stats 2023-07-13 12:19:04,495 (natcap.invest.urban_cooling_model) urban_cooling_model.calculate_uhi_result_vector(919) INFO Creating uhi_results.shp 2023-07-13 12:19:04,567 (natcap.invest.urban_cooling_model) urban_cooling_model.execute(837) INFO Urban Cooling Model complete. 2023-07-13 12:19:04,567 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 4.41s 2023-07-13 12:19:04,569 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished