official-apps/SeismicHazardForecasting
11
0
Fork 2
Code Issues Pull Requests Activity

Update src/seismic_hazard_forecasting.py

Browse Source
This commit is contained in:
ftong
2026-06-23 10:56:54 +02:00
parent 3576b9eb96
commit 72dc427a0e
1 changed files with 7 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
src/seismic_hazard_forecasting.py
+7 -4
View File
@@ -501,8 +501,7 @@ verbose: {verbose}")
grid_y_max = int(ceil(y_max / grid_dim) * grid_dim) grid_y_max = int(ceil(y_max / grid_dim) * grid_dim)
grid_y_min = int(floor(y_min / grid_dim) * grid_dim) grid_y_min = int(floor(y_min / grid_dim) * grid_dim)
grid_lat_max, grid_lon_max = utm.to_latlon(grid_x_max, grid_y_max, utm_zone_number, utm_zone_letter)
grid_lat_min, grid_lon_min = utm.to_latlon(grid_x_min, grid_y_min, utm_zone_number, utm_zone_letter)
# rectangular grid # rectangular grid
nx = int((grid_x_max - grid_x_min) / grid_dim) + 1 nx = int((grid_x_max - grid_x_min) / grid_dim) + 1
@@ -517,6 +516,10 @@ verbose: {verbose}")
nx = ny nx = ny
grid_x_max = int(grid_x_min + (nx - 1) * grid_dim) grid_x_max = int(grid_x_min + (nx - 1) * grid_dim)
# update grid extent in lat/lon
grid_lat_max, grid_lon_max = utm.to_latlon(grid_x_max, grid_y_max, utm_zone_number, utm_zone_letter)
grid_lat_min, grid_lon_min = utm.to_latlon(grid_x_min, grid_y_min, utm_zone_number, utm_zone_letter)
# new x and y range # new x and y range
x_range = np.linspace(grid_x_min, grid_x_max, nx) x_range = np.linspace(grid_x_min, grid_x_max, nx)
y_range = np.linspace(grid_y_min, grid_y_max, ny) y_range = np.linspace(grid_y_min, grid_y_max, ny)
@@ -861,8 +864,8 @@ verbose: {verbose}")
vmax_hd = np.nanmax(iml_grid_hd) vmax_hd = np.nanmax(iml_grid_hd)
# generate image overlay # generate image overlay
north, south = lat.max(), lat.min() # Latitude range north, south = grid_lat_max, grid_lat_min # Latitude range
east, west = lon.max(), lon.min() # Longitude range east, west = grid_lon_max, grid_lon_min # Longitude range
bounds = [[south, west], [north, east]] bounds = [[south, west], [north, east]]
map_center = [np.mean([north, south]), np.mean([east, west])] map_center = [np.mean([north, south]), np.mean([east, west])]