Skip to content

Commit

Permalink
add more lines to roads if statement
Browse files Browse the repository at this point in the history
  • Loading branch information
Santonia27 committed Jan 10, 2025
1 parent 0253250 commit 523202b
Showing 1 changed file with 9 additions and 9 deletions.
18 changes: 9 additions & 9 deletions flood_adapt/database_builder/create_database.py
Original file line number Diff line number Diff line change
Expand Up @@ -1087,16 +1087,16 @@ def update_fiat_elevation(self):
roads = gpd.read_file(roads_path).to_crs(dem.spatial_ref.crs_wkt)
roads["geometry"] = roads.geometry.centroid # get centroids

x_points = xr.DataArray(roads["geometry"].x, dims="points")
y_points = xr.DataArray(roads["geometry"].y, dims="points")
roads["elev"] = (
dem.sel(x=x_points, y=y_points, band=1, method="nearest").to_numpy()
* conversion_factor
)
x_points = xr.DataArray(roads["geometry"].x, dims="points")
y_points = xr.DataArray(roads["geometry"].y, dims="points")
roads["elev"] = (
dem.sel(x=x_points, y=y_points, band=1, method="nearest").to_numpy()
* conversion_factor
)

exposure.loc[
exposure["Primary Object Type"] == "road", "Ground Floor Height"
] = 0
exposure.loc[
exposure["Primary Object Type"] == "road", "Ground Floor Height"
] = 0
exposure = exposure.merge(
roads[["Object ID", "elev"]], on="Object ID", how="left"
)
Expand Down

0 comments on commit 523202b

Please sign in to comment.