APOG_public/SP2XR
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Remove warnings for log10(0) in histograms calcualtion

Browse Source
This commit is contained in:
bertoz_b
2025-11-17 16:44:14 +01:00
parent 5a52ca2fdd
commit fc9cf7c861
1 changed files with 13 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

20
src/sp2xr/calibration.py
View File

@@ -228,14 +228,20 @@ def calibrate_single_particle(
) )
ddf_cal["time_lag"] = ddf_cal["Incand Peak Time"] - ddf_cal["Scatter Peak Time"] ddf_cal["time_lag"] = ddf_cal["Incand Peak Time"] - ddf_cal["Scatter Peak Time"]
# Suppress divide by zero warning for log10(0) - expected for particles with no signal # Calculate ratio using map_partitions to properly suppress warnings in Dask workers
import warnings def calculate_ratio_with_suppression(df):
"""Calculate log10 ratio with warning suppression for zero values."""
with np.errstate(divide="ignore", invalid="ignore"):
df["ratio_inc_scatt"] = np.log10(df["Incand relPeak"]) / np.log10(
df["Scatter relPeak"]
)
return df
with warnings.catch_warnings(): meta_with_ratio = ddf_cal._meta.copy()
warnings.filterwarnings("ignore", message="divide by zero encountered in log10") meta_with_ratio["ratio_inc_scatt"] = pd.Series([], dtype="float64")
ddf_cal["ratio_inc_scatt"] = np.log10(ddf_cal["Incand relPeak"]) / np.log10( ddf_cal = ddf_cal.map_partitions(
ddf_cal["Scatter relPeak"] calculate_ratio_with_suppression, meta=meta_with_ratio
) )
ddf_cal = define_flags(ddf_cal, instr_config, run_config) ddf_cal = define_flags(ddf_cal, instr_config, run_config)