本文整理汇总了Python中windml.datasets.nrel.NREL类的典型用法代码示例。如果您正苦于以下问题:Python NREL类的具体用法?Python NREL怎么用?Python NREL使用的例子?那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了NREL类的15个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: test_nmar_destroyer
def test_nmar_destroyer(self):
turbine = NREL().get_turbine(NREL.park_id['tehachapi'], 2004)
timeseries = turbine.get_measurements()[:1000]
damaged, indices = NMARDestroyer().destroy(timeseries, percentage=.50,\
min_length=10, max_length=50)
misses = MissingDataFinder().find(damaged, 600)
assert(len(misses) > 0)
示例2: test_mreg_interpolation_multi
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' : 10, '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)
示例3: compute_mse
def compute_mse(regressor, horizon):
# get wind park and corresponding target.
windpark = NREL().get_windpark(NREL.park_id['tehachapi'], 3, 2004, 2005)
target = windpark.get_target()
# use power mapping for pattern-label mapping.
feature_window = 3
mapping = PowerMapping()
X = mapping.get_features_park(windpark, feature_window, horizon)
y = mapping.get_labels_turbine(target, feature_window, horizon)
# train roughly for the year 2004, test for 2005.
train_to = int(math.floor(len(X) * 0.5))
test_to = len(X)
train_step, test_step = 25, 25
X_train=X[:train_to:train_step]
y_train=y[:train_to:train_step]
X_test=X[train_to:test_to:test_step]
y_test=y[train_to:test_to:test_step]
if(regressor == 'svr'):
reg = SVR(kernel='rbf', epsilon=0.1, C = 100.0,\
gamma = 0.0001).fit(X_train,y_train)
mse = mean_squared_error(reg.predict(X_test),y_test)
elif(regressor == 'knn'):
reg = KNeighborsRegressor(10, 'uniform').fit(X_train,y_train)
mse = mean_squared_error(reg.predict(X_test),y_test)
return mse
示例4: amount_of_windmills
def amount_of_windmills(radius, park):
target = NREL.park_id[park]
ds = NREL()
windpark = ds.get_windpark(target, radius, 2004, 2005)
target = ds.get_windmill(target, 2004, 2005)
windmills = windpark.get_windmills()
return len(windmills)
示例5: test_marthres_destroyer
def test_marthres_destroyer(self):
turbine = NREL().get_turbine(NREL.park_id['tehachapi'], 2004)
timeseries = turbine.get_measurements()[:1000]
damaged, indices = MARThresDestroyer().destroy(timeseries, percentage=.50,\
lower_bound = 0, upper_bound = 20)
misses = MissingDataFinder().find(damaged, 600)
assert(len(misses) > 0)
示例6: compute_mse
def compute_mse(regressor, param):
# get wind park and corresponding target. forecast is for the target
# turbine
park_id = NREL.park_id['tehachapi']
windpark = NREL().get_windpark(park_id, 3, 2004)
target = windpark.get_target()
# use power mapping for pattern-label mapping. Feature window length
# is 3 time steps and time horizon (forecast) is 3 time steps.
feature_window = 6
horizon = 3
mapping = PowerMapping()
X = mapping.get_features_park(windpark, feature_window, horizon)
Y = mapping.get_labels_turbine(target, feature_window, horizon)
# train roughly for the year 2004.
train_to = int(math.floor(len(X) * 0.5))
# test roughly for the year 2005.
test_to = len(X)
# train and test only every fifth pattern, for performance.
train_step, test_step = 5, 5
if(regressor == 'rf'):
# random forest regressor
reg = RandomForestRegressor(n_estimators=param, criterion='mse')
reg = reg.fit(X[0:train_to:train_step], Y[0:train_to:train_step])
y_hat = reg.predict(X[train_to:test_to:test_step])
elif(regressor == 'knn'):
# TODO the regressor does not need to be newly trained in
# the case of KNN
reg = KNeighborsRegressor(param, 'uniform')
# fitting the pattern-label pairs
reg = reg.fit(X[0:train_to:train_step], Y[0:train_to:train_step])
y_hat = reg.predict(X[train_to:test_to:test_step])
else:
raise Exception("No regressor set.")
# naive is also known as persistence model.
naive_hat = zeros(len(y_hat), dtype = float32)
for i in range(0, len(y_hat)):
# naive label is the label as horizon time steps before.
# we have to consider to use only the fifth label here, too.
naive_hat[i] = Y[train_to + (i * test_step) - horizon]
# computing the mean squared errors of Linear and naive prediction.
mse_y_hat, mse_naive_hat = 0, 0
for i in range(0, len(y_hat)):
y = Y[train_to + (i * test_step)]
mse_y_hat += (y_hat[i] - y) ** 2
mse_naive_hat += (naive_hat[i] - y) ** 2
mse_y_hat /= float(len(y_hat))
mse_naive_hat /= float(len(y_hat))
return mse_y_hat, mse_naive_hat
示例7: test_backward_copy_interpolation
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)
示例8: test_mreg_interpolation
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' : 10, '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)
示例9: test_topological_interpolation
def test_topological_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 = NMARDestroyer().destroy(measurements, percentage=.80,\
min_length=10, max_length=100)
tloc = (target.longitude, target.latitude)
neighbors = windpark.get_turbines()[:-1]
nseries = [t.get_measurements()[300:500] for t in neighbors]
nlocs = [(t.longitude, t.latitude) for t in neighbors]
t_hat = TopologicInterpolation().interpolate(\
damaged, method="topologic",\
timestep=timestep, location=tloc,\
neighbor_series = nseries,\
neighbor_locations = nlocs)
misses = MissingDataFinder().find(t_hat, timestep)
assert(measurements.shape[0] == t_hat.shape[0])
assert(len(misses) < 1)
示例10: setUpClass
def setUpClass(cls):
ds = NREL()
cls.windmill = ds.get_windmill(NREL.park_id['tehachapi'], 2004, 2005)
cls.windpark = ds.get_windpark(NREL.park_id['tehachapi'], 3, 2004, 2005)
cls.pmapping = PowerMapping()
cls.pdmapping = PowerDiffMapping()
示例11: test_nrel_repair
def test_nrel_repair(self):
ds = NREL()
target = ds.get_turbine(NREL.park_id['tehachapi'], 2005)
measurements = target.get_measurements()[:43504]
measurements = NRELRepair().repair(measurements)
assert(NRELRepair().validate(measurements))
示例12: NREL
"""
# Author: Oliver Kramer <[email protected]>
# License: BSD 3 clause
from __future__ import print_function
import sklearn
import numpy as np
import pylab as plt
from sklearn import manifold, decomposition
from builtins import range
from windml.datasets.nrel import NREL
# load data and define parameters / training and test sequences
K = 30
ds = NREL()
windpark = ds.get_windpark(NREL.park_id['tehachapi'], 10, 2004)
X = np.array(windpark.get_powermatrix())
X_train = X[:2000]
X_test = X[2000:2000 + 200 * 4]
# computation of ISOMAP projection
print("computation of ISOMAP projection")
X_latent = manifold.Isomap(K, n_components=2).fit_transform(X_train)
# computation of sequence of closest embedded patterns
sequence = []
for x in X_test:
win = 0
示例13: NREL
"""
Histogram of Wind Speeds
-------------------------------------------------------------
Histograms of wind speeds of a turbine near Cheyenne in the year 2004.
"""
# Author: Jendrik Poloczek <[email protected]>
# License: BSD 3 clause
import matplotlib.pyplot as plt
from pylab import plt
from windml.datasets.nrel import NREL
ds = NREL()
turbine = ds.get_turbine(NREL.park_id['cheyenne'], 2004)
speeds = list(map(lambda x : x[2], turbine.measurements))
plt.hist(speeds, color="#c4d8eb", bins=10, normed = 1)
plt.show()
示例14: NREL
"""
Time-Series of Wind Speed and Power
--------------------------------------------------
This example plots a time-series of a single
wind mill in the wind park 'tehachapi'.
"""
from matplotlib import dates
import matplotlib.pylab as plt
import numpy as np
import datetime, time
from windml.datasets.nrel import NREL
from windml.visualization.plot_timeseries import plot_timeseries
ds = NREL()
mill = ds.get_windmill(NREL.park_id['tehachapi'], 2004)
plot_timeseries(mill)
示例15: NREL
import matplotlib.pylab as plt
import datetime, time
import numpy as np
from numpy import array, matrix
from sklearn.grid_search import GridSearchCV
from sklearn.cross_validation import KFold
from sklearn import __version__ as sklearn_version
from sklearn.svm import SVR
from sklearn.neighbors import KNeighborsRegressor
from windml.datasets.nrel import NREL
from windml.visualization.plot_response_curve import plot_response_curve
ds = NREL()
turbine = ds.get_turbine(NREL.park_id['palmsprings'], 2004, 2006)
timeseries = turbine.get_measurements()
max_speed = 40
skip = 1
# plot true values as blue points
speed = [m[2] for m in timeseries[::skip]]
score = [m[1] for m in timeseries[::skip]]
# Second Plot: KNN-Interpolation
# Built patterns und labels
X_train = speed[0:len(speed):1]
Y_train = score[0:len(score):1]