|
| 1 | +import numpy as np |
| 2 | +import pandas as pd |
| 3 | +import geopandas as gpd |
| 4 | +import shapely.geometry as geom |
| 5 | +import pytest |
| 6 | + |
| 7 | +from cubedynamics import fire_time_hull as fth |
| 8 | +from cubedynamics.fire_time_hull import ( |
| 9 | + FireEventDaily, |
| 10 | + HullClimateSummary, |
| 11 | + TimeHull, |
| 12 | + plot_climate_filled_hull, |
| 13 | +) |
| 14 | + |
| 15 | + |
| 16 | +def _synthetic_hull(days: int = 3, verts_per_layer: int = 4) -> TimeHull: |
| 17 | + verts = [] |
| 18 | + t_days = [] |
| 19 | + for layer in range(days): |
| 20 | + z = 10.0 + layer # intentionally non-1-based to catch index/order bugs |
| 21 | + t_days.extend([z] * verts_per_layer) |
| 22 | + for j in range(verts_per_layer): |
| 23 | + angle = 2.0 * np.pi * j / verts_per_layer |
| 24 | + radius = 1.0 + 0.1 * layer |
| 25 | + verts.append([radius * np.cos(angle), radius * np.sin(angle), z]) |
| 26 | + verts = np.asarray(verts, dtype=float) |
| 27 | + |
| 28 | + tris = [] |
| 29 | + for i in range(days - 1): |
| 30 | + for j in range(verts_per_layer): |
| 31 | + jn = (j + 1) % verts_per_layer |
| 32 | + v1 = i * verts_per_layer + j |
| 33 | + v2 = i * verts_per_layer + jn |
| 34 | + v3 = (i + 1) * verts_per_layer + jn |
| 35 | + v4 = (i + 1) * verts_per_layer + j |
| 36 | + tris.append([v1, v2, v3]) |
| 37 | + tris.append([v1, v3, v4]) |
| 38 | + |
| 39 | + gdf = gpd.GeoDataFrame( |
| 40 | + {"id": [1], "date": [pd.Timestamp("2020-07-01")], "geometry": [geom.box(0, 0, 1, 1)]}, |
| 41 | + crs="EPSG:4326", |
| 42 | + ) |
| 43 | + event = FireEventDaily(1, gdf, pd.Timestamp("2020-07-01"), pd.Timestamp("2020-07-03"), 0.0, 0.0) |
| 44 | + return TimeHull( |
| 45 | + event=event, |
| 46 | + verts_km=verts, |
| 47 | + tris=np.asarray(tris, dtype=int), |
| 48 | + t_days_vert=np.asarray(t_days, dtype=float), |
| 49 | + t_norm_vert=np.linspace(0, 1, len(t_days)), |
| 50 | + metrics={"days": float(days), "scale_km": 1.0, "volume_km2_days": 1.0, "surface_km_day": 1.0}, |
| 51 | + ) |
| 52 | + |
| 53 | + |
| 54 | +@pytest.fixture |
| 55 | +def _plotly_stub(monkeypatch): |
| 56 | + class _Mesh3d: |
| 57 | + def __init__(self, **kwargs): |
| 58 | + self.__dict__.update(kwargs) |
| 59 | + |
| 60 | + class _Figure: |
| 61 | + def __init__(self, data): |
| 62 | + self.data = data |
| 63 | + self.layout = type("Layout", (), {})() |
| 64 | + self.layout.scene = type("Scene", (), {})() |
| 65 | + self.layout.scene.aspectratio = type("Aspect", (), {})() |
| 66 | + |
| 67 | + def update_layout(self, **kwargs): |
| 68 | + scene = kwargs.get("scene", {}) |
| 69 | + aspect = scene.get("aspectratio", {}) |
| 70 | + self.layout.scene.aspectratio.z = aspect.get("z") |
| 71 | + |
| 72 | + def write_image(self, *_args, **_kwargs): |
| 73 | + return None |
| 74 | + |
| 75 | + monkeypatch.setattr(fth.go, "Mesh3d", _Mesh3d, raising=False) |
| 76 | + monkeypatch.setattr(fth.go, "Figure", _Figure, raising=False) |
| 77 | + |
| 78 | + |
| 79 | +def test_plot_climate_filled_hull_attaches_scalars_by_layer_order(_plotly_stub): |
| 80 | + hull = _synthetic_hull(days=3, verts_per_layer=4) |
| 81 | + summary = HullClimateSummary( |
| 82 | + values_inside=np.array([1.0, 2.0, 3.0]), |
| 83 | + values_outside=np.array([0.0, 0.0]), |
| 84 | + per_day_mean=pd.Series([1.0, 5.0, 9.0], index=pd.date_range("2020-07-01", periods=3, freq="D")), |
| 85 | + ) |
| 86 | + |
| 87 | + fig = plot_climate_filled_hull(hull, summary, color_limits=None) |
| 88 | + got = np.asarray(fig.data[0].intensity, dtype=float) |
| 89 | + expected = np.repeat(np.array([1.0, 5.0, 9.0]), 4) |
| 90 | + np.testing.assert_allclose(got, expected) |
| 91 | + |
| 92 | + |
| 93 | +def test_plot_climate_filled_hull_debug_z_mode_and_z_exaggeration(_plotly_stub): |
| 94 | + hull = _synthetic_hull(days=3, verts_per_layer=4) |
| 95 | + summary = HullClimateSummary( |
| 96 | + values_inside=np.array([1.0]), |
| 97 | + values_outside=np.array([0.0]), |
| 98 | + per_day_mean=pd.Series([2.0, 2.1, 2.2], index=pd.date_range("2020-07-01", periods=3, freq="D")), |
| 99 | + ) |
| 100 | + |
| 101 | + fig = plot_climate_filled_hull( |
| 102 | + hull, |
| 103 | + summary, |
| 104 | + scalar_debug_mode="z", |
| 105 | + z_exaggeration=2.8, |
| 106 | + color_limits=None, |
| 107 | + ) |
| 108 | + got = np.asarray(fig.data[0].intensity, dtype=float) |
| 109 | + np.testing.assert_allclose(got, hull.verts_km[:, 2]) |
| 110 | + assert float(fig.layout.scene.aspectratio.z) == 2.8 |
0 commit comments