2023-10-13 08:41:55,870 (natcap.invest.utils) utils.prepare_workspace(161) Level 100 Writing log messages to [/development/InVEST/Outputs/CoastalVulnerability/InVEST-natcap.invest.coastal_vulnerability-log-2023-10-13--08_41_55.txt] 2023-10-13 08:41:55,871 (__main__) cli.main(459) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.coastal_vulnerability 3.14.0: aoi_vector_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/aoi_grandbahama_utm.shp bathymetry_raster_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/bathymetry.tif dem_averaging_radius 900 dem_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/dem_srtm_grandbahama.tif geomorphology_fill_value 3 geomorphology_vector_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/geomorphology_grandbahama.shp habitat_table_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/GrandBahama_Habitats/Natural_Habitats.csv landmass_vector_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/landmass_polygon.shp max_fetch_distance 15000 model_resolution 1000 n_workers -1 population_radius 500 population_raster_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/population_grandbahama.tif results_suffix GB shelf_contour_vector_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/continental_shelf_polyline_global.shp slr_vector_path workspace_dir /development/InVEST/Outputs/CoastalVulnerability wwiii_vector_path /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/WaveWatchIII_global.shp 2023-10-13 08:41:56,160 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2037) INFO starting reprojection 2023-10-13 08:41:56,163 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2077) INFO reprojection 100.0% complete on clipped_projected_landmass_GB.gpkg 2023-10-13 08:41:56,169 (natcap.invest.coastal_vulnerability) coastal_vulnerability.warp_and_mask_bathymetry(2187) INFO Masking positive values from bathymetry 2023-10-13 08:41:56,189 (natcap.invest.coastal_vulnerability) coastal_vulnerability.clip_and_project_raster(3022) WARNING /development/InVEST/SampleData/FromAppDownload/CoastalVulnerability/bathymetry.tif is unprojected and will be warped to match the AOI and resampled to a pixel size of 1000.0 meters 2023-10-13 08:41:56,200 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-13 08:41:56,201 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-13 08:41:56,228 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-13 08:41:56,229 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-13 08:41:56,232 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(1100) INFO preparing landmass geometry and interpolating shore points 2023-10-13 08:41:56,272 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(1149) INFO indexing geometry of landmass and creating shore points 2023-10-13 08:41:57,264 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(1215) INFO Finished creating 620 shore points in AOI 2023-10-13 08:41:57,316 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(1249) INFO Building spatial index for Wave Watch III points 2023-10-13 08:42:00,647 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(1295) INFO Interpolating Wave Watch III data to shore points 2023-10-13 08:42:00,947 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(1352) INFO Finished interpolating Wave Watch III data to shore points 2023-10-13 08:42:00,949 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1500) INFO Calculating wind exposure 2023-10-13 08:42:02,601 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1562) INFO Casting rays and extracting bathymetry values 2023-10-13 08:42:07,608 (natcap.invest.coastal_vulnerability) coastal_vulnerability.logger_callback(3437) INFO Wind exposure 48.23% complete. 2023-10-13 08:42:12,126 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1703) INFO Finished calculating wind exposure 2023-10-13 08:42:12,137 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wave_exposure(1880) INFO Calculating wave exposure 2023-10-13 08:42:13,923 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wave_exposure(1987) INFO Finished calculating wave exposure 2023-10-13 08:42:13,925 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_surge_exposure(2006) INFO Calculating surge potential 2023-10-13 08:42:14,379 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_surge_exposure(2071) INFO Finished calculating surge potential 2023-10-13 08:42:14,381 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(2095) INFO Calculating relief exposure 2023-10-13 08:42:14,546 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(470) INFO starting stats_worker 2023-10-13 08:42:14,547 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(476) INFO started stats_worker 2023-10-13 08:42:14,563 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(555) INFO 100.0% complete 2023-10-13 08:42:14,563 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(558) INFO Waiting for raster stats worker result. 2023-10-13 08:42:14,613 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(2136) WARNING 1 shore points are missing values after aggregating /development/InVEST/Outputs/CoastalVulnerability/intermediate/relief/positive_dem_GB.tif. No valid pixels were found within the search radius (900.0 meters) around these points. 2023-10-13 08:42:14,613 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(2141) INFO Finished calculating relief exposure 2023-10-13 08:42:14,620 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2435) INFO Searching for coral within 2000.0 meters of shore points 2023-10-13 08:42:14,631 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Ring Self-intersection at or near point 126818.06030000125 3035614.6246999996 2023-10-13 08:42:14,634 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Ring Self-intersection at or near point 140468.06020000053 3030544.6245999997 2023-10-13 08:42:14,646 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Ring Self-intersection at or near point 114488.06190000083 2946034.6234000013 2023-10-13 08:42:14,704 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2540) INFO Finished searching for coral 2023-10-13 08:42:14,707 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2435) INFO Searching for coastal_forest within 1000.0 meters of shore points 2023-10-13 08:42:14,716 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Ring Self-intersection at or near point 134495.52229999978 2940707.8449000004 2023-10-13 08:42:14,768 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2540) INFO Finished searching for coastal_forest 2023-10-13 08:42:14,771 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_raster_habitat(2377) INFO Searching for seagrass within 500.0 meters of shore points 2023-10-13 08:42:14,788 (natcap.invest.coastal_vulnerability) coastal_vulnerability._aggregate_raster_values_in_radius(3211) WARNING search radius around point (101535.59551799361, 2958347.559996179) extends beyond bounds of /development/InVEST/Outputs/CoastalVulnerability/intermediate/habitats/clipped_projected_seagrass__GB.tif 2023-10-13 08:42:14,800 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_raster_habitat(2409) INFO Finished searching for seagrass 2023-10-13 08:42:14,802 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2435) INFO Searching for mangrove within 2000.0 meters of shore points 2023-10-13 08:42:14,837 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] Ring Self-intersection at or near point 151947.11389999912 2970620.4354999983 2023-10-13 08:42:15,179 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_vector_habitat(2540) INFO Finished searching for mangrove 2023-10-13 08:42:15,184 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_habitat_rank(2599) INFO coral_GB: found protecting 120 shore points 2023-10-13 08:42:15,185 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_habitat_rank(2599) INFO coastal_forest_GB: found protecting 110 shore points 2023-10-13 08:42:15,186 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_habitat_rank(2599) INFO seagrass_GB: found protecting 382 shore points 2023-10-13 08:42:15,186 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_habitat_rank(2599) INFO mangrove_GB: found protecting 610 shore points 2023-10-13 08:42:15,217 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2037) INFO starting reprojection 2023-10-13 08:42:15,246 (osgeo) utils._log_gdal_errors(97) WARNING [errno 1] A geometry of type MULTILINESTRING is inserted into layer geomorphology_grandbahama of geometry type LINESTRING, which is not normally allowed by the GeoPackage specification, but the driver will however do it. To create a conformant GeoPackage, if using ogr2ogr, the -nlt option can be used to override the layer geometry type. This warning will no longer be emitted for this combination of layer and feature geometry type. 2023-10-13 08:42:15,248 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(2077) INFO reprojection 100.0% complete on geomorphology_projected_GB.gpkg 2023-10-13 08:42:15,257 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2641) INFO Assigning geomorphology rank 2023-10-13 08:42:15,259 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2646) INFO Build spatial index of geomorphology segments 2023-10-13 08:42:15,293 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2687) INFO searching for geomorphology in point buffers 2023-10-13 08:42:19,068 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2716) WARNING ************************************************** 1 points inherited the geomorphology_fill_value (3) because no geomorphology segments were found within the search distance (500.00 meters) around those points. You may wish to add geomorphology segments near these points, or you may wish to view the geomorphology and landmass polygon inputs in a GIS and edit one or both layers to make them align more closely. This subset of points has been saved to /development/InVEST/Outputs/CoastalVulnerability/intermediate/geomorphology/shore_points_missing_geomorphology_GB.gpkg ************************************************** 2023-10-13 08:42:19,069 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2730) INFO Finished assigning geomorphology rank 2023-10-13 08:42:19,072 (natcap.invest.coastal_vulnerability) coastal_vulnerability.aggregate_population_density(2237) INFO Aggregating population data 2023-10-13 08:42:19,117 (natcap.invest.coastal_vulnerability) coastal_vulnerability.aggregate_population_density(2276) WARNING 434 shore points are missing values after aggregating /development/InVEST/Outputs/CoastalVulnerability/intermediate/population/clipped_projected_pop_GB.tif. No valid pixels were found within the search radius (500.0 meters) around these points. 2023-10-13 08:42:19,117 (natcap.invest.coastal_vulnerability) coastal_vulnerability.aggregate_population_density(2281) INFO Finished aggregating population data 2023-10-13 08:42:19,118 (natcap.invest.coastal_vulnerability) coastal_vulnerability.assemble_results_and_calculate_exposure(2769) INFO Assembling exposure variables 2023-10-13 08:42:19,249 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_final_risk(2883) INFO calculating final risk scores 2023-10-13 08:42:19,399 (natcap.invest.utils) utils.prepare_workspace(169) INFO Elapsed time: 23.53s 2023-10-13 08:42:19,399 (natcap.invest.utils) utils.prepare_workspace(172) INFO Execution finished