2023-02-25 10:23:38,517 (natcap.invest.utils) utils.prepare_workspace(162) Level 100 Writing log messages to C:/Users/User/Desktop/All documents and folders/SampleData_CVModel/Output_25Feb2023\InVEST-Coastal-Vulnerability-log-2023-02-25--10_23_38.txt 2023-02-25 10:23:38,519 (natcap.invest.ui.model) model._logged_target(1643) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.coastal_vulnerability 3.10.2: aoi_vector_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\AOI_R.S.shp bathymetry_raster_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Bathymetry_R.S_UTM.tif dem_averaging_radius 10 dem_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Clipped DEM.tif geomorphology_fill_value 5 geomorphology_vector_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Geomorphology_Projected_270123.shp habitat_table_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Habitats_RS\CaostalHabitats_NomsyO.csv landmass_vector_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\ne_10m_land.shp max_fetch_distance 5000 model_resolution 500 n_workers -1 population_radius 500 population_raster_path results_suffix tif shelf_contour_vector_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\continental_shelf_polyline_global.shp slr_field UNKNOWN slr_vector_path workspace_dir C:/Users/User/Desktop/All documents and folders/SampleData_CVModel/Output_25Feb2023 wwiii_vector_path C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\WaveWatchIII_global.shp 2023-02-25 10:23:38,524 (natcap.invest.coastal_vulnerability) coastal_vulnerability.execute(386) INFO Validating arguments 2023-02-25 10:23:38,524 (natcap.invest.validation) validation._wrapped_validate_func(1026) DEBUG Using ARG_SPEC for validation 2023-02-25 10:23:38,525 (natcap.invest.validation) validation._wrapped_validate_func(1030) INFO Starting whole-model validation with ARGS_SPEC 2023-02-25 10:23:38,525 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,526 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,527 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,532 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,533 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,575 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,576 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,604 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,605 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,615 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,616 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,636 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,636 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,666 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,666 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,676 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,676 (natcap.invest.validation) validation.timeout(716) DEBUG Starting file checking thread with timeout=5 2023-02-25 10:23:38,700 (natcap.invest.validation) validation.timeout(727) DEBUG File checking thread completed. 2023-02-25 10:23:38,785 (natcap.invest.validation) validation._wrapped_validate_func(1067) DEBUG Validation warnings: [] 2023-02-25 10:23:39,197 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2023-02-25 10:23:39,201 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on clipped_projected_landmass_tif.gpkg 2023-02-25 10:23:39,221 (natcap.invest.coastal_vulnerability) coastal_vulnerability.warp_and_mask_bathymetry(1806) INFO Masking positive values from bathymetry 2023-02-25 10:23:39,271 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1980) DEBUG transforming bounding box from [239202.86820000038, 476528.3000000008, 339233.12844779296, 540899.6248816821] 2023-02-25 10:23:39,273 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1985) DEBUG transforming bounding to [239202.86820000038, 476528.30000000075, 339233.12844779296, 540899.6248816821] 2023-02-25 10:23:39,292 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-02-25 10:23:39,294 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 20804 2023-02-25 10:23:39,294 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-02-25 10:23:39,305 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-02-25 10:23:39,306 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-02-25 10:23:39,315 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(722) INFO preparing landmass geometry and interpolating shore points 2023-02-25 10:23:39,678 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(771) INFO indexing geometry of landmass and creating shore points 2023-02-25 10:23:39,695 (py.warnings) warnings._showwarnmsg(109) WARNING natcap\invest\coastal_vulnerability.py:791: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry. 2023-02-25 10:23:40,459 (natcap.invest.coastal_vulnerability) coastal_vulnerability.prepare_landmass_line_index_and_interpolate_shore_points(833) INFO Finished creating 714 shore points in AOI 2023-02-25 10:23:40,499 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(867) INFO Building spatial index for Wave Watch III points 2023-02-25 10:23:47,490 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(913) INFO Interpolating Wave Watch III data to shore points 2023-02-25 10:23:48,482 (natcap.invest.coastal_vulnerability) coastal_vulnerability.interpolate_wwiii_to_shore(970) INFO Finished interpolating Wave Watch III data to shore points 2023-02-25 10:23:48,490 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1116) INFO Calculating wind exposure 2023-02-25 10:23:57,072 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1178) INFO Casting rays and extracting bathymetry values 2023-02-25 10:24:02,084 (natcap.invest.coastal_vulnerability) coastal_vulnerability.logger_callback(3059) INFO Wind exposure 47.06% complete. 2023-02-25 10:24:06,514 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wind_exposure(1320) INFO Finished calculating wind exposure 2023-02-25 10:24:06,534 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wave_exposure(1499) INFO Calculating wave exposure 2023-02-25 10:24:15,483 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_wave_exposure(1606) INFO Finished calculating wave exposure 2023-02-25 10:24:15,490 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_surge_exposure(1625) INFO Calculating surge potential 2023-02-25 10:24:16,046 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_surge_exposure(1690) INFO Finished calculating surge potential 2023-02-25 10:24:16,062 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(1714) INFO Calculating relief exposure 2023-02-25 10:24:16,374 (natcap.invest.coastal_vulnerability) coastal_vulnerability.clip_and_project_raster(2636) WARNING C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Clipped DEM.tif is unprojected and will be warped to match the AOI and resampled to a pixel size of 500 meters 2023-02-25 10:24:16,377 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1980) DEBUG transforming bounding box from [6.720414980900341, 4.3584021407486935, 7.497637203122563, 4.8150688074153605] 2023-02-25 10:24:16,378 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1985) DEBUG transforming bounding to [247001.48131756246, 481909.8424278702, 333391.74805300194, 532644.705654721] 2023-02-25 10:24:16,509 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(396) INFO starting stats_worker 2023-02-25 10:24:16,511 (pygeoprocessing.geoprocessing_core) threading.run(870) DEBUG stats worker PID: 20804 2023-02-25 10:24:16,511 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(402) INFO started stats_worker 2023-02-25 10:24:16,513 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(482) INFO 100.0% complete 2023-02-25 10:24:16,513 (pygeoprocessing.geoprocessing) geoprocessing.raster_calculator(485) INFO Waiting for raster stats worker result. 2023-02-25 10:24:16,581 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(1755) WARNING 16 shore points are missing values after aggregating C:/Users/User/Desktop/All documents and folders/SampleData_CVModel/Output_25Feb2023\intermediate\relief\positive_dem_tif.tif. No valid pixels were found within the search radius (10.0 meters) around these points. 2023-02-25 10:24:16,582 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_relief_exposure(1760) INFO Finished calculating relief exposure 2023-02-25 10:24:16,598 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_raster_habitat(1991) INFO Searching for E10N10 within 5000 meters of shore points 2023-02-25 10:24:16,921 (natcap.invest.coastal_vulnerability) coastal_vulnerability.clip_and_project_raster(2636) WARNING C:\Users\User\Desktop\All documents and folders\SampleData_CVModel\Habitats_RS\E10N10.tif is unprojected and will be warped to match the AOI and resampled to a pixel size of 500 meters 2023-02-25 10:24:16,924 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1980) DEBUG transforming bounding box from [6.6752962051433355, 4.313158268043368, 7.542568932329335, 4.860249177079568] 2023-02-25 10:24:16,925 (pygeoprocessing.geoprocessing) geoprocessing.warp_raster(1985) DEBUG transforming bounding to [241975.92722669555, 476897.3461333212, 338386.36616378464, 537659.8294261484] 2023-02-25 10:24:17,095 (natcap.invest.coastal_vulnerability) coastal_vulnerability.search_for_raster_habitat(2023) INFO Finished searching for E10N10 2023-02-25 10:24:17,106 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_habitat_rank(2213) INFO E10N10_tif: found protecting 714 shore points 2023-02-25 10:24:17,249 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1741) INFO starting reprojection 2023-02-25 10:24:17,252 (osgeo) utils._log_gdal_errors(98) WARNING [errno 1] A geometry of type MULTIPOLYGON is inserted into layer Geomorphology_Projected_270123 of geometry type POLYGON, 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-02-25 10:24:17,264 (pygeoprocessing.geoprocessing) geoprocessing.reproject_vector(1783) INFO reprojection 100.0% complete on geomorphology_projected_tif.gpkg 2023-02-25 10:24:17,284 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2256) INFO Assigning geomorphology rank 2023-02-25 10:24:17,289 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2261) INFO Build spatial index of geomorphology segments 2023-02-25 10:24:17,372 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2303) INFO searching for geomorphology in point buffers 2023-02-25 10:24:35,560 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_geomorphology_exposure(2347) INFO Finished assigning geomorphology rank 2023-02-25 10:24:35,567 (natcap.invest.coastal_vulnerability) coastal_vulnerability.assemble_results_and_calculate_exposure(2386) INFO Assembling exposure variables 2023-02-25 10:24:36,312 (natcap.invest.coastal_vulnerability) coastal_vulnerability.calculate_final_risk(2497) INFO calculating final risk scores 2023-02-25 10:24:36,987 (natcap.invest.utils) utils.prepare_workspace(170) INFO Elapsed time: 58.47s 2023-02-25 10:24:36,988 (natcap.invest.utils) utils.prepare_workspace(173) INFO Execution finished