|
15 | 15 |
|
16 | 16 | import numpy as np |
17 | 17 | import pytest |
| 18 | +import xarray as xr |
18 | 19 |
|
19 | 20 | from cmethods import adjust |
| 21 | +from cmethods.core import _add_cmethods_metadata |
20 | 22 | from cmethods.distribution import ( |
21 | 23 | detrended_quantile_mapping, |
22 | 24 | quantile_delta_mapping, |
@@ -130,3 +132,70 @@ def test_adjust_failing_no_group_for_distribution(datasets: dict) -> None: |
130 | 132 | n_quantiles=100, |
131 | 133 | group="time.month", |
132 | 134 | ) |
| 135 | + |
| 136 | + |
| 137 | +def test_add_cmethods_metadata_with_dataarray() -> None: |
| 138 | + """Test that _add_cmethods_metadata adds correct attributes to a DataArray""" |
| 139 | + data = xr.DataArray( |
| 140 | + np.array([1, 2, 3, 4, 5]), |
| 141 | + dims=["time"], |
| 142 | + coords={"time": np.arange(5)}, |
| 143 | + ) |
| 144 | + |
| 145 | + result = _add_cmethods_metadata( |
| 146 | + data, |
| 147 | + method="linear_scaling", |
| 148 | + kind="+", |
| 149 | + n_quantiles=100, |
| 150 | + group="time.month", |
| 151 | + ) |
| 152 | + |
| 153 | + assert "cmethods_version" in result.attrs |
| 154 | + assert "cmethods_method" in result.attrs |
| 155 | + assert "cmethods_timestamp" in result.attrs |
| 156 | + assert "cmethods_source" in result.attrs |
| 157 | + |
| 158 | + assert result.attrs["cmethods_method"] == "linear_scaling" |
| 159 | + assert result.attrs["cmethods_kind"] == "+" |
| 160 | + assert result.attrs["cmethods_n_quantiles"] == "100" |
| 161 | + assert result.attrs["cmethods_group"] == "time.month" |
| 162 | + assert result.attrs["cmethods_source"] == "https://github.com/btschwertfeger/python-cmethods" |
| 163 | + assert "UTC" in result.attrs["cmethods_timestamp"] |
| 164 | + |
| 165 | + |
| 166 | +def test_add_cmethods_metadata_with_dataset() -> None: |
| 167 | + """Test that _add_cmethods_metadata adds correct attributes to a Dataset""" |
| 168 | + data = xr.Dataset( |
| 169 | + { |
| 170 | + "temperature": xr.DataArray( |
| 171 | + np.array([1, 2, 3, 4, 5]), |
| 172 | + dims=["time"], |
| 173 | + coords={"time": np.arange(5)}, |
| 174 | + ), |
| 175 | + }, |
| 176 | + ) |
| 177 | + |
| 178 | + result = _add_cmethods_metadata(data, method="quantile_mapping") |
| 179 | + |
| 180 | + assert "cmethods_version" in result.attrs |
| 181 | + assert "cmethods_method" in result.attrs |
| 182 | + assert "cmethods_timestamp" in result.attrs |
| 183 | + assert "cmethods_source" in result.attrs |
| 184 | + assert result.attrs["cmethods_method"] == "quantile_mapping" |
| 185 | + |
| 186 | + |
| 187 | +def test_add_cmethods_metadata_optional_params() -> None: |
| 188 | + """Test that _add_cmethods_metadata handles optional parameters correctly""" |
| 189 | + data = xr.DataArray( |
| 190 | + np.array([1, 2, 3]), |
| 191 | + dims=["time"], |
| 192 | + coords={"time": np.arange(3)}, |
| 193 | + ) |
| 194 | + |
| 195 | + result = _add_cmethods_metadata(data, method="variance_scaling") |
| 196 | + |
| 197 | + assert "cmethods_method" in result.attrs |
| 198 | + assert result.attrs["cmethods_method"] == "variance_scaling" |
| 199 | + assert "cmethods_kind" not in result.attrs |
| 200 | + assert "cmethods_n_quantiles" not in result.attrs |
| 201 | + assert "cmethods_group" not in result.attrs |
0 commit comments