2020-04-16 13:51:06,322 utils.prepare_workspace(115) INFO Writing log messages to C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2\InVEST-Coastal-Vulnerability-log-2020-04-16--13_51_06.txt 2020-04-16 13:51:06,324 model._logged_target(1642) Level 100 Starting model with parameters: Arguments for InVEST natcap.invest.coastal_vulnerability 3.8.0: aoi_vector_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/Base_Data/MD_AOI.shp bathymetry_raster_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/TopoBathy_forCV_m.tif dem_averaging_radius 1000 dem_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/TopoBathy_forCV_m.tif geomorphology_fill_value 1 geomorphology_vector_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/MD_geomorph_Albers.shp habitat_table_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/NaturalHabitat/Natural_Habitats.csv landmass_vector_path C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/Shoreline_clip.shp max_fetch_distance 1000 model_resolution 1000 n_workers -1 population_radius population_raster_path results_suffix shelf_contour_vector_path C:/InVEST_3.8.0_x86/sample_data/CoastalVulnerability/continental_shelf_polyline_global.shp slr_field UNKNOWN slr_vector_path workspace_dir C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2 wwiii_vector_path C:/InVEST_3.8.0_x86/sample_data/CoastalVulnerability/WaveWatchIII_global.shp 2020-04-16 13:51:06,370 coastal_vulnerability.execute(285) INFO Validating arguments 2020-04-16 13:51:06,389 validation._wrapped_validate_func(839) DEBUG Using ARG_SPEC for validation 2020-04-16 13:51:06,389 validation._wrapped_validate_func(843) INFO Starting whole-model validation with ARGS_SPEC 2020-04-16 13:51:06,570 validation._wrapped_validate_func(887) DEBUG Validation warnings: [] 2020-04-16 13:51:06,908 utils._log_gdal_errors(72) WARNING [errno 1] A geometry of type MULTILINESTRING is inserted into layer tmp_clipped_landmass of geometry type MULTIPOLYGON, 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. 2020-04-16 13:51:06,910 utils._log_gdal_errors(72) WARNING [errno 1] A geometry of type LINESTRING is inserted into layer tmp_clipped_landmass of geometry type MULTIPOLYGON, 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. 2020-04-16 13:51:07,411 geoprocessing.reproject_vector(1640) INFO starting reprojection 2020-04-16 13:51:07,413 utils._log_gdal_errors(72) WARNING [errno 1] A geometry of type MULTILINESTRING is inserted into layer tmp_clipped_landmass of geometry type MULTIPOLYGON, 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. 2020-04-16 13:51:07,414 utils._log_gdal_errors(72) WARNING [errno 1] A geometry of type LINESTRING is inserted into layer tmp_clipped_landmass of geometry type MULTIPOLYGON, 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. 2020-04-16 13:51:07,418 geoprocessing.reproject_vector(1683) INFO reprojection 100.0% complete on clipped_projected_landmass.gpkg 2020-04-16 13:51:13,160 coastal_vulnerability.warp_and_mask_bathymetry(1749) INFO Masking positive values from bathymetry 2020-04-16 13:51:18,197 geoprocessing.logger_callback(3239) INFO Warp 82.0% complete C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2\intermediate\wind_wave\TopoBathy_forCV_m_clipped.tif 2020-04-16 13:51:19,331 geoprocessing.logger_callback(3239) INFO Warp 100.0% complete C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2\intermediate\wind_wave\TopoBathy_forCV_m_clipped.tif 2020-04-16 13:51:20,361 geoprocessing.warp_raster(1841) DEBUG transforming bounding box from [1599594.229508998, 1826497.4953873642, 1809204.229508998, 2039767.4953873642] 2020-04-16 13:51:20,362 geoprocessing.warp_raster(1846) DEBUG transforming bounding to [1599594.2295089976, 1826497.4953873628, 1809204.2295089988, 2039767.4953873653] 2020-04-16 13:51:25,461 geoprocessing.logger_callback(3239) INFO Warp 63.0% complete C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2\intermediate\wind_wave\clipped_projected_bathy.tif 2020-04-16 13:51:28,870 geoprocessing.logger_callback(3239) INFO Warp 100.0% complete C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD2\intermediate\wind_wave\clipped_projected_bathy.tif 2020-04-16 13:51:29,843 geoprocessing.raster_calculator(375) INFO starting stats_worker 2020-04-16 13:51:29,845 geoprocessing.raster_calculator(381) INFO started stats_worker 2020-04-16 13:51:32,784 geoprocessing.raster_calculator(454) INFO 100.0%% complete 2020-04-16 13:51:32,784 geoprocessing.raster_calculator(457) INFO signaling stats worker to terminate 2020-04-16 13:51:32,784 geoprocessing.raster_calculator(459) INFO Waiting for raster stats worker result. 2020-04-16 13:51:32,784 threading.run(870) DEBUG payload is None, terminating 2020-04-16 13:51:33,468 coastal_vulnerability.prepare_landmass_line_index(627) INFO preparing landmass geometry 2020-04-16 13:51:33,624 coastal_vulnerability.prepare_landmass_line_index(680) INFO indexing geometry of landmass 2020-04-16 13:51:33,638 coastal_vulnerability.interpolate_shore_points(711) INFO creating shore points along edge of landmass 2020-04-16 13:51:33,687 coastal_vulnerability.interpolate_shore_points(796) INFO Finished creating 0 shore points in AOI 2020-04-16 13:51:33,710 coastal_vulnerability.interpolate_wwiii_to_shore(826) INFO Building spatial index for Wave Watch III points 2020-04-16 13:51:42,107 coastal_vulnerability.interpolate_wwiii_to_shore(872) INFO Interpolating Wave Watch III data to shore points 2020-04-16 13:51:42,108 coastal_vulnerability.interpolate_wwiii_to_shore(929) INFO Finished interpolating Wave Watch III data to shore points 2020-04-16 13:51:42,295 coastal_vulnerability.calculate_wind_exposure(1075) INFO Calculating wind exposure 2020-04-16 13:51:42,547 coastal_vulnerability.calculate_wind_exposure(1137) INFO Casting rays and extracting bathymetry values 2020-04-16 13:51:42,558 coastal_vulnerability.calculate_wind_exposure(1278) INFO Finished calculating wind exposure 2020-04-16 13:51:42,573 coastal_vulnerability.calculate_wave_exposure(1454) INFO Calculating wave exposure 2020-04-16 13:51:42,663 coastal_vulnerability.calculate_wave_exposure(1561) INFO Finished calculating wave exposure 2020-04-16 13:51:42,735 coastal_vulnerability.calculate_surge_exposure(1580) INFO Calculating surge potential 2020-04-16 13:51:44,739 coastal_vulnerability.calculate_surge_exposure(1650) INFO Finished calculating surge potential 2020-04-16 13:51:50,291 coastal_vulnerability.calculate_relief_exposure(1673) INFO Calculating relief exposure 2020-04-16 13:51:50,313 geoprocessing.warp_raster(1841) DEBUG transforming bounding box from [-1000.0000000004324, -1000.0000000030727, 1009.9999999995676, 1009.9999999969273] 2020-04-16 13:51:50,314 geoprocessing.warp_raster(1846) DEBUG transforming bounding to [-1000.0000000010518, -1000.0000000073214, 1010.0000000000052, 1010.0000000014323] 2020-04-16 13:51:50,324 geoprocessing.raster_calculator(375) INFO starting stats_worker 2020-04-16 13:51:50,325 geoprocessing.raster_calculator(381) INFO started stats_worker 2020-04-16 13:51:50,327 geoprocessing.raster_calculator(454) INFO 100.0%% complete 2020-04-16 13:51:50,327 geoprocessing.raster_calculator(457) INFO signaling stats worker to terminate 2020-04-16 13:51:50,328 geoprocessing.raster_calculator(459) INFO Waiting for raster stats worker result. 2020-04-16 13:51:50,328 threading.run(870) DEBUG payload is None, terminating 2020-04-16 13:51:50,328 threading.run(870) WARNING No valid pixels were received, sending None. 2020-04-16 13:51:50,338 coastal_vulnerability.calculate_relief_exposure(1703) INFO Finished calculating relief exposure 2020-04-16 13:51:50,361 coastal_vulnerability.search_for_habitat(1897) INFO Searching for Coastal_Forest within 1000 meters of shore points 2020-04-16 13:51:50,770 coastal_vulnerability.search_for_habitat(1963) WARNING No valid features exist in C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/NaturalHabitat/MD_CoastalForests_Albers.shp 2020-04-16 13:51:50,776 coastal_vulnerability.search_for_habitat(2005) INFO Finished searching for Coastal_Forest in proximity to shore points 2020-04-16 13:51:50,900 coastal_vulnerability.search_for_habitat(1897) INFO Searching for marsh within 1000 meters of shore points 2020-04-16 13:51:55,509 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1730525.7159806297 1863951.0598547906 2020-04-16 13:51:55,938 coastal_vulnerability.logger_callback(2892) INFO transforming marsh 66.94% complete. 2020-04-16 13:51:58,360 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1779921.3189727543 1862108.2297313288 2020-04-16 13:51:58,784 coastal_vulnerability.search_for_habitat(1963) WARNING No valid features exist in C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/NaturalHabitat/MD_emergent_wetlands_Albers.shp 2020-04-16 13:51:58,789 coastal_vulnerability.search_for_habitat(2005) INFO Finished searching for marsh in proximity to shore points 2020-04-16 13:51:58,841 coastal_vulnerability.search_for_habitat(1897) INFO Searching for sav within 500 meters of shore points 2020-04-16 13:52:00,863 coastal_vulnerability.search_for_habitat(1963) WARNING No valid features exist in C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/NaturalHabitat/MD_SAV_Albers.shp 2020-04-16 13:52:00,868 coastal_vulnerability.search_for_habitat(2005) INFO Finished searching for sav in proximity to shore points 2020-04-16 13:52:00,893 coastal_vulnerability.search_for_habitat(1897) INFO Searching for oyster within 100 meters of shore points 2020-04-16 13:52:00,952 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1690428.6771922419 1936192.7033486385 2020-04-16 13:52:01,023 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1704271.9733599168 1909756.9800894614 2020-04-16 13:52:01,054 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1699157.3273012638 1844745.6416061055 2020-04-16 13:52:01,104 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1692744.2019274402 1845537.7745831264 2020-04-16 13:52:01,137 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1670432.4964605442 1970655.5038107559 2020-04-16 13:52:01,232 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1674141.0924226737 1848931.8912851764 2020-04-16 13:52:01,239 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1649680.2793778866 1857707.6938342899 2020-04-16 13:52:01,386 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1685587.0079652793 1960423.0540095763 2020-04-16 13:52:01,412 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1661412.6218901328 1953862.1305724094 2020-04-16 13:52:01,416 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1660404.3573847346 1938484.1320184851 2020-04-16 13:52:01,433 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1661367.3744638236 1954461.0140032424 2020-04-16 13:52:01,436 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1660236.4223919469 1953726.4893873446 2020-04-16 13:52:01,477 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1678863.220498272 1962035.6003939845 2020-04-16 13:52:01,486 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1676516.5845371392 1984387.3134887591 2020-04-16 13:52:01,489 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1669711.3252976693 1972801.63767666 2020-04-16 13:52:01,521 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1668671.2633056666 1959626.4745778916 2020-04-16 13:52:01,543 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1667531.4533009457 1927391.9105538493 2020-04-16 13:52:01,805 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1704410.064379493 1908860.0192193445 2020-04-16 13:52:01,859 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1690809.4328907479 1964024.0889319442 2020-04-16 13:52:01,894 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1668407.0121014337 1929359.9414511072 2020-04-16 13:52:02,116 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1699727.0061757127 1848313.4846059345 2020-04-16 13:52:02,126 utils._log_gdal_errors(72) WARNING [errno 1] Ring Self-intersection at or near point 1698314.4766255682 1849931.9444692507 2020-04-16 13:52:02,316 coastal_vulnerability.search_for_habitat(1963) WARNING No valid features exist in C:/Users/kjd15/Documents/GIS/USCAproject/CoastalVulnerability/MD/MD_Inputs/NaturalHabitat/MD_Oysters_Albers.shp 2020-04-16 13:52:02,323 coastal_vulnerability.search_for_habitat(2005) INFO Finished searching for oyster in proximity to shore points 2020-04-16 13:52:02,339 coastal_vulnerability.calculate_habitat_rank(2034) INFO Coastal_Forest: found protecting 0 shore points 2020-04-16 13:52:02,341 coastal_vulnerability.calculate_habitat_rank(2034) INFO marsh: found protecting 0 shore points 2020-04-16 13:52:02,343 coastal_vulnerability.calculate_habitat_rank(2034) INFO sav: found protecting 0 shore points 2020-04-16 13:52:02,344 coastal_vulnerability.calculate_habitat_rank(2034) INFO oyster: found protecting 0 shore points 2020-04-16 13:52:02,475 geoprocessing.reproject_vector(1640) INFO starting reprojection 2020-04-16 13:52:07,476 geoprocessing.(1646) INFO reprojection approximately 48.7% complete on geomorphology_projected.shp 2020-04-16 13:52:12,476 geoprocessing.(1646) INFO reprojection approximately 98.6% complete on geomorphology_projected.shp 2020-04-16 13:52:12,589 geoprocessing.reproject_vector(1683) INFO reprojection 100.0% complete on geomorphology_projected.shp 2020-04-16 13:52:12,761 coastal_vulnerability.calculate_geomorphology_exposure(2089) INFO Assigning geomorphology rank 2020-04-16 13:52:12,764 coastal_vulnerability.calculate_geomorphology_exposure(2094) INFO Build spatial index of geomorphology segments 2020-04-16 13:52:18,967 coastal_vulnerability.calculate_geomorphology_exposure(2134) INFO searching for geomorphology in point buffers 2020-04-16 13:52:18,978 coastal_vulnerability.calculate_geomorphology_exposure(2177) INFO Finished assigning geomorphology rank 2020-04-16 13:52:19,175 coastal_vulnerability.assemble_results_and_calculate_exposure(2216) INFO Assembling exposure variables 2020-04-16 13:52:19,264 warnings._showwarnmsg(110) WARNING C:\InVEST_3.8.0_x86\invest-3-x86\numpy\lib\nanfunctions.py:1368: RuntimeWarning: Mean of empty slice 2020-04-16 13:52:19,264 model._logged_target(1648) ERROR Exception while executing Traceback (most recent call last): File "site-packages\natcap\invest\ui\model.py", line 1645, in _logged_target File "site-packages\natcap\invest\coastal_vulnerability.py", line 602, in execute File "site-packages\natcap\invest\coastal_vulnerability.py", line 2267, in assemble_results_and_calculate_exposure File "site-packages\natcap\invest\coastal_vulnerability.py", line 2404, in _bin_values_to_percentiles File "<__array_function__ internals>", line 6, in searchsorted File "site-packages\numpy\core\fromnumeric.py", line 1341, in searchsorted File "site-packages\numpy\core\fromnumeric.py", line 61, in _wrapfunc ValueError: object of too small depth for desired array 2020-04-16 13:52:19,265 model._logged_target(1651) INFO Execution finished 2020-04-16 13:52:19,266 utils.prepare_workspace(121) INFO Elapsed time: 1m 12.939999999999998s