From a1c0ae36bb1a29c7d547b5c315360ae7c3741479 Mon Sep 17 00:00:00 2001 From: ftong Date: Tue, 9 Sep 2025 14:41:02 +0200 Subject: [PATCH] set a minimum number of computed grid values to trigger upscaling of grid image --- src/seismic_hazard_forecasting.py | 12 ++++++++---- 1 file changed, 8 insertions(+), 4 deletions(-) diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py index 0703b50..3cd03c6 100644 --- a/src/seismic_hazard_forecasting.py +++ b/src/seismic_hazard_forecasting.py @@ -519,10 +519,14 @@ verbose: {verbose}") dtype=np.float64) # this reduces values to 8 decimal places iml_grid_tmp = np.nan_to_num(iml_grid[j]) # change nans to zeroes - # upscale the grid - up_factor = 4 - iml_grid_hd = resize(iml_grid_tmp, (up_factor * len(iml_grid_tmp), up_factor * len(iml_grid_tmp)), - mode='reflect', anti_aliasing=False) + # upscale the grid if there are at least 10 grid values + if np.count_nonzero(iml_grid_tmp) >= 10: + up_factor = 4 + iml_grid_hd = resize(iml_grid_tmp, (up_factor * len(iml_grid_tmp), up_factor * len(iml_grid_tmp)), + mode='reflect', anti_aliasing=False) + else: + iml_grid_hd = iml_grid_tmp + iml_grid_hd[iml_grid_hd == 0.0] = np.nan # change zeroes back to nan # trim edges so the grid is not so blocky