weac.parser.snowscope_parser module

Parse SnowScope CSV profiles into a force–density signal for WEAC.

Density is derived from the SnowScope hardness/force signal F (kPa) with a semilog model D = a ln(F) + b (kg/m³):

  • French Alps SP2/SMP, Hagenmüller 2018 (HAGENMULLER2018, default): a = 71.4, b = 21.5 (published as ρ = 21.5 + 71.4 log(σ))

Pass semilog_slope / semilog_intercept on SnowScopeParser or SnowScopeParser.extract_profile() to use custom a and b instead of a preset (label CUSTOM on the returned profile).

class weac.parser.snowscope_parser.SnowScopeProfile(depth_mm, penetration_resistance_kPa, density_kg_m3, density_method)[source]

Bases: object

Aligned SnowScope penetration resistance and derived density.

Parameters:

  • depth_mm (ndarray)

  • penetration_resistance_kPa (ndarray)

  • density_kg_m3 (ndarray)

  • density_method (Literal['HAGENMULLER2018', 'CUSTOM'])

depth_mm: ndarray
penetration_resistance_kPa: ndarray
density_kg_m3: ndarray
density_method: Literal['HAGENMULLER2018', 'CUSTOM']
__init__(depth_mm, penetration_resistance_kPa, density_kg_m3, density_method)
Parameters:

  • depth_mm (ndarray)

  • penetration_resistance_kPa (ndarray)

  • density_kg_m3 (ndarray)

  • density_method (Literal['HAGENMULLER2018', 'CUSTOM'])

Return type:

None

class weac.parser.snowscope_parser.SnowScopeParser(file_path, *, density_method='HAGENMULLER2018', semilog_slope=None, semilog_intercept=None)[source]

Bases: object

Parser for SnowScope CSV files.

Parameters:

  • file_path (str)

  • density_method (DensityMethod)

  • semilog_slope (float | None)

  • semilog_intercept (float | None)

__init__(file_path, *, density_method='HAGENMULLER2018', semilog_slope=None, semilog_intercept=None)[source]
Parameters:

  • file_path (str)

  • density_method (Literal['HAGENMULLER2018'])

  • semilog_slope (float | None)

  • semilog_intercept (float | None)

extract_profile(*, density_method=None, semilog_slope=None, semilog_intercept=None)[source]

Return depth, penetration resistance, and density for one parameterization.

Converts penetration resistance (kPa) to density (kg/m³) via the semilog model D = a ln(F) + b, clipping sub-1 kg/m³ values.

Coefficients are chosen in order: call-level semilog_slope / semilog_intercept, then call-level density_method, then constructor custom semilog (if set), else the constructor density_method preset.

Parameters:

  • density_method (Literal['HAGENMULLER2018'] | None)

  • semilog_slope (float | None)

  • semilog_intercept (float | None)

Return type:

SnowScopeProfile

extract_layers(slope_angle_deg=0.0, *, method='bin', density_method=None, semilog_slope=None, semilog_intercept=None, layer_thickness_mm=None, gradient_threshold=12.0)[source]

Segment the SnowScope density profile into WEAC slab layers (top-down).

The SnowScope probe descends along the global vertical, so depth_mm is a plumb depth while WEAC Layer.h is slope-normal. As in the SnowPilot parser, thicknesses are converted plumb -> slope-normal by cos(phi) via plumb_to_slope_normal(). The SnowScope file records no slope, so slope_angle_deg defaults to 0 (no scaling); scale only when the angle is known (pass it explicitly).

Parameters:

  • slope_angle_deg (float) – Slope angle [deg from horizontal]. 0 leaves h as recorded; non-zero scales h by cos(phi).

  • method (Literal['bin', 'gradient']) – Layering mode. "bin" groups at a fixed thickness (or the native spacing); "gradient" cuts where the density gradient over a fixed 2.5 mm span exceeds gradient_threshold.

  • density_method (Literal['HAGENMULLER2018'] | None) – Optional preset override for this call only.

  • semilog_slope (float | None) – Optional semilog slope a in D = a ln(F) + b for this call; must be paired with semilog_intercept.

  • semilog_intercept (float | None) – Optional semilog intercept b; must be paired with semilog_slope. Takes precedence over density_method.

  • layer_thickness_mm (float | None) – Bin mode only; None groups at native sample spacing with a 1 mm floor, a value groups samples into ~that thickness. Ignored in gradient mode.

  • gradient_threshold (float) – Gradient mode only; cut threshold T in kg/m^3/mm over the 2.5 mm span. Default 12.

Returns:

(layers, density_methods) with layers ordered surface -> ground and density_methods = [parameterization] * len(layers).

Return type:

tuple[list[Layer], list[str]]