Source code for bsrn.qc.k_index
"""
BSRN Level 5 checks - radiometric index tests (k-tests).
BSRN 5 级检查 - 辐射指数测试 (k-测试)。
"""
import numpy as np
import pandas as pd
[docs]
def kb_kt_test(ghi, bni, bni_extra, zenith):
"""
Check if beam transmittance ($k_b$) is less than clearness index ($k_t$) [1]_.
检查直射透射率 ($k_b$) 是否小于物理上限制的晴朗指数 ($k_t$)。
Parameters
----------
ghi : numeric or Series
Global horizontal irradiance ($G_h$). [W/m^2]
水平总辐照度 ($G_h$)。[瓦/平方米]
bni : numeric or Series
Beam normal irradiance ($B_n$). [W/m^2]
法向直接辐照度 ($B_n$)。[瓦/平方米]
bni_extra : numeric or Series
Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]
地外法向辐照度 ($E_{0n}$)。[瓦/平方米]
zenith : numeric or Series
Solar zenith angle ($Z$). [degrees]
太阳天顶角 ($Z$)。[度]
Returns
-------
flags : Series or ndarray
Boolean flags (True = Pass). [bool]
布尔标记（True = 通过）。[布尔值]
References
----------
.. [1] Forstinger, A., et al. (2021). Expert quality control of solar
radiation ground data sets. In SWC 2021: ISES Solar World Congress.
International Solar Energy Society.
"""
mu0 = np.cos(np.radians(zenith))
kb = bni / bni_extra
kt = ghi / (bni_extra * mu0)
# Domain: GHI > 50 and kb > 0 and kt > 0 / 适用范围: GHI > 50 且 kb > 0 且 kt > 0
in_domain = (ghi > 50) & (kb > 0) & (kt > 0)
condition_met = kb < kt
if hasattr(in_domain, 'iloc'):
return (~in_domain) | condition_met
else:
return (not in_domain) or condition_met
[docs]
def kb_limit_test(bni, bni_extra, elevation, ghi):
"""
Check if beam transmittance ($k_b$) stays within absolute physical limits based on elevation [1]_.
根据海拔检查直射透射率 ($k_b$) 是否在绝对物理限值内。
Parameters
----------
bni : numeric or Series
Beam normal irradiance ($B_n$). [W/m^2]
法向直接辐照度 ($B_n$)。[瓦/平方米]
bni_extra : numeric or Series
Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]
地外法向辐照度 ($E_{0n}$)。[瓦/平方米]
elevation : numeric
Site elevation. [m]
站点海拔。[米]
ghi : numeric or Series
Global horizontal irradiance ($G_h$). [W/m^2]
水平总辐照度 ($G_h$)。[瓦/平方米]
Returns
-------
flags : Series or ndarray
Boolean flags (True = Pass). [bool]
布尔标记（True = 通过）。[布尔值]
References
----------
.. [1] Forstinger, A., et al. (2021). Expert quality control of solar
radiation ground data sets. In SWC 2021: ISES Solar World Congress.
International Solar Energy Society.
"""
kb = bni / bni_extra
# Domain: GHI > 50 and kb > 0 / 适用范围: GHI > 50 且 kb > 0
in_domain = (ghi > 50) & (kb > 0)
condition_met = kb < (1100 + 0.03 * elevation) / bni_extra
if hasattr(in_domain, 'iloc'):
return (~in_domain) | condition_met
else:
return (not in_domain) or condition_met
[docs]
def kt_limit_test(ghi, bni_extra, zenith):
"""
Check if clearness index ($k_t$) is within physically possible limits [1]_.
检查晴朗指数 ($k_t$) 是否在物理可能范围内。
Parameters
----------
ghi : numeric or Series
Global horizontal irradiance ($G_h$). [W/m^2]
水平总辐照度 ($G_h$)。[瓦/平方米]
bni_extra : numeric or Series
Extraterrestrial beam normal irradiance ($E_{0n}$). [W/m^2]
地外法向辐照度 ($E_{0n}$)。[瓦/平方米]
zenith : numeric or Series
Solar zenith angle ($Z$). [degrees]
太阳天顶角 ($Z$)。[度]
Returns
-------
flags : Series or ndarray
Boolean flags (True = Pass). [bool]
布尔标记（True = 通过）。[布尔值]
References
----------
.. [1] Forstinger, A., et al. (2021). Expert quality control of solar
radiation ground data sets. In SWC 2021: ISES Solar World Congress.
International Solar Energy Society.
"""
mu0 = np.cos(np.radians(zenith))
kt = ghi / (bni_extra * mu0)
# Domain: GHI > 50 and kt > 0 / 适用范围: GHI > 50 且 kt > 0
in_domain = (ghi > 50) & (kt > 0)
condition_met = kt < 1.35
if hasattr(in_domain, 'iloc'):
return (~in_domain) | condition_met
else:
return (not in_domain) or condition_met
