def test_mreg_interpolation_multi(self):
park_id = NREL.park_id['tehachapi']
windpark = NREL().get_windpark(park_id, 3, 2004)
target = windpark.get_target()
timestep = 600
measurements = target.get_measurements()[300:350]
damaged, indices = MARDestroyer().destroy(measurements, percentage=.50)
before_misses = MissingDataFinder().find(damaged, timestep)
neighbors = windpark.get_turbines()[:-1]
count_neighbors = len(neighbors)
reg = 'knn' # KNeighborsRegressor(10, 'uniform')
regargs = {'n' : 8, 'variant' : 'uniform'}
processed = 0
missed = {k : count_neighbors for k in indices}
exclude = []
damaged_nseries = []
for neighbor in neighbors:
nseries = neighbor.get_measurements()[300:350]
damaged, indices = MARDestroyer().destroy(nseries, percentage=.50, exclude=exclude)
for index in indices:
if(index not in missed.keys()):
missed[index] = count_neighbors
missed[index] -= 1
if(missed[index] == 1):
exclude.append(index) # exclude in next iterations
damaged_nseries.append(damaged)
t_hat = MRegInterpolation().interpolate(damaged, timestep=timestep,\
neighbor_series=damaged_nseries, reg=reg, regargs=regargs)
after_misses = MissingDataFinder().find(t_hat, timestep)
assert(len(after_misses) < 1)
def destroy(timeseries, method, **args):
methods = {
'mar': MARDestroyer().destroy,
'nmar': NMARDestroyer().destroy,
'mar_with_threshold': MARThresDestroyer().destroy
}
return methods[method](timeseries, **args)
def test_backward_copy_interpolation(self):
park_id = NREL.park_id['tehachapi']
windpark = NREL().get_windpark(park_id, 10, 2004)
target = windpark.get_target()
timestep = 600
measurements = target.get_measurements()[300:500]
damaged, indices = MARDestroyer().destroy(measurements, percentage=.50)
before_misses = MissingDataFinder().find(damaged, timestep)
t_hat = BackwardCopy().interpolate(measurements, timestep=timestep)
after_misses = MissingDataFinder().find(t_hat, timestep)
assert(measurements.shape[0] == t_hat.shape[0])
assert(len(after_misses) < 1)
def test_mreg_interpolation(self):
park_id = NREL.park_id['tehachapi']
windpark = NREL().get_windpark(park_id, 3, 2004)
target = windpark.get_target()
timestep = 600
measurements = target.get_measurements()[300:500]
damaged, indices = MARDestroyer().destroy(measurements, percentage=.50)
before_misses = MissingDataFinder().find(damaged, timestep)
neighbors = windpark.get_turbines()[:-1]
reg = 'knn' # KNeighborsRegressor(10, 'uniform')
regargs = {'n' : 8, 'variant' : 'uniform'}
nseries = [t.get_measurements()[300:500] for t in neighbors]
t_hat = MRegInterpolation().interpolate(damaged, timestep=timestep,\
neighbor_series=nseries, reg=reg, regargs=regargs)
after_misses = MissingDataFinder().find(t_hat, timestep)
assert(len(after_misses) < 1)
def test_mar_destroyer(self):
turbine = NREL().get_turbine(NREL.park_id['tehachapi'], 2004)
timeseries = turbine.get_measurements()[:1000]
damaged, indices = MARDestroyer().destroy(timeseries, percentage=.50)
misses = MissingDataFinder().find(damaged, 600)
assert(len(misses) > 0)