def detect_clear_days(self, th=0.1):
if self.filled_data_matrix is None:
print('Generate a filled data matrix first.')
return
clear_days = find_clear_days(self.filled_data_matrix, th=th)
self.daily_flags.flag_clear_cloudy(clear_days)
return
def test_find_clear_days(self):
filepath = Path(__file__).parent.parent
data_file_path = (
filepath
/ "fixtures"
/ "clear_day_detection"
/ "one_year_power_signals_1.csv"
)
with open(data_file_path) as file:
data = np.loadtxt(file, delimiter=",")
expected_data_file_path = (
filepath / "fixtures" / "clear_day_detection" / "one_year_weights_1.csv"
)
with open(expected_data_file_path) as file:
expected_output = np.loadtxt(file, delimiter=",")
expected_output = expected_output >= 1e-3
# Underling solar-data-tools uses MOSEK solver and
# if it's not used, try with ECOS solver.
# However, fails with ECOS solver and raises cvx.SolverError.
try:
actual_output = find_clear_days(data)
except (cvx.SolverError, ValueError):
self.skipTest(
"This test uses MOSEK solver"
+ "because default ECOS solver fails with large data. "
+ "Unless MOSEK is installed, this test fails."
)
else:
np.testing.assert_array_equal(expected_output, actual_output)
def test_clear_day_weights(self):
data_file_path = os.path.abspath(
os.path.join(
os.path.dirname(__file__),
"../fixtures/clear_day_detection/one_year_power_signals_1.csv")
)
with open(data_file_path) as file:
data = np.loadtxt(file, delimiter=',')
expected_data_file_path = os.path.abspath(
os.path.join(
os.path.dirname(__file__),
"../fixtures/clear_day_detection/one_year_weights_1.csv"))
with open(expected_data_file_path) as file:
expected_output = np.loadtxt(file, delimiter=',')
# Underling solar-data-tools uses MOSEK solver and
# if it's not used, try with ECOS solver.
# However, fails with ECOS solver and raises cvx.SolverError.
try:
actual_output = find_clear_days(data, boolean_out=False)
except (cvx.SolverError, ValueError):
self.skipTest(
"This test uses MOSEK solver" +
"because default ECOS solver fails with large data. " +
"Unless MOSEK is installed, this test fails.")
else:
np.testing.assert_array_almost_equal(expected_output,
actual_output, 4)
def obtain_weights(self, power_signals_d):
try:
from solardatatools.clear_day_detection import find_clear_days
except ImportError:
print('Weights not set!')
print('Please make sure you have solar-data-tools installed')
weights = np.ones(power_signals_d.shape[1])
else:
weights = find_clear_days(power_signals_d, boolean_out=False)
return weights
def test_clear_day_weights(self):
data_file_path = os.path.abspath(
os.path.join(
os.path.dirname(__file__),
"../fixtures/clear_day_detection/one_year_power_signals_1.csv")
)
with open(data_file_path) as file:
data = np.loadtxt(file, delimiter=',')
expected_data_file_path = os.path.abspath(
os.path.join(
os.path.dirname(__file__),
"../fixtures/clear_day_detection/one_year_weights_1.csv"))
with open(expected_data_file_path) as file:
expected_output = np.loadtxt(file, delimiter=',')
actual_output = find_clear_days(data, boolean_out=False)
np.testing.assert_array_equal(expected_output, actual_output)