本文整理汇总了Python中windml.datasets.nrel.NREL类的典型用法代码示例。如果您正苦于以下问题:Python NREL类的具体用法?Python NREL怎么用?Python NREL使用的例子?那么恭喜您, 这里精选的类代码示例或许可以为您提供帮助。
在下文中一共展示了NREL类的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的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)
开发者ID:Bengt,项目名称:windml,代码行数:7,代码来源:preprocessing_test.py
示例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)
开发者ID:Bengt,项目名称:windml,代码行数:35,代码来源:preprocessing_test.py
示例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
开发者ID:DeeplearningMachineLearning,项目名称:windml,代码行数:28,代码来源:forecast_horizon.py
示例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)
开发者ID:tverrbjelke,项目名称:windml,代码行数:7,代码来源:windmills_in_radius.py
示例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)
开发者ID:Bengt,项目名称:windml,代码行数:7,代码来源:preprocessing_test.py
示例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
开发者ID:Bengt,项目名称:windml,代码行数:57,代码来源:compare_regressors_param.py
示例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)
开发者ID:Bengt,项目名称:windml,代码行数:13,代码来源:preprocessing_test.py
示例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)
开发者ID:Bengt,项目名称:windml,代码行数:18,代码来源:preprocessing_test.py
示例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)
开发者ID:Bengt,项目名称:windml,代码行数:24,代码来源:preprocessing_test.py
示例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()
开发者ID:tverrbjelke,项目名称:windml,代码行数:6,代码来源:mapping_test.py
示例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))
开发者ID:Bengt,项目名称:windml,代码行数:6,代码来源:preprocessing_test.py
示例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
开发者ID:cigroup-ol,项目名称:windml,代码行数:31,代码来源:sequence.py
示例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()
开发者ID:cigroup-ol,项目名称:windml,代码行数:20,代码来源:windspeed_histogram.py
示例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)
开发者ID:tverrbjelke,项目名称:windml,代码行数:19,代码来源:plot_timeseries.py
示例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]
开发者ID:Bengt,项目名称:windml,代码行数:31,代码来源:svr_response_curve.py
示例16: setUpClass
def setUpClass(cls):
ds = NREL()
cls.turbine = ds.get_turbine(NREL.park_id['tehachapi'], 2004)
cls.windpark = ds.get_windpark(NREL.park_id['tehachapi'], 3, 2004)
开发者ID:cigroup-ol,项目名称:windml,代码行数:4,代码来源:visualization_test.py
示例17: NREL
"""
Topography of a Wind Windpark Near Tehachapi
-------------------------------------------------------------------------
This example shows the topography of a wind park near Tehachapi. The red dots
illustrate the locations of wind mills.
"""
from windml.datasets.nrel import NREL
from windml.visualization.show_coord_topo import show_coord_topo
radius = 30
name = 'tehachapi'
windpark = NREL().get_windpark(NREL.park_id['tehachapi'], 30, 2004)
print "Working on windpark around target mill", str(windpark.get_target_idx())
print "Plotting windpark ..."
show_coord_topo(windpark)
开发者ID:tverrbjelke,项目名称:windml,代码行数:20,代码来源:show_coord_topo.py
示例18: NREL
"""
Topography of a Windpark
-------------------------------------------------------------------------
This example shows the topography of a wind park near Tehachapi. The black dots
illustrate the locations of turbines. The red dot is the target turbine.
"""
# Author: Nils A. Treiber <[email protected]>
# License: BSD 3 clause
from windml.datasets.nrel import NREL
from windml.visualization.show_coord_topo import show_coord_topo
radius = 30
name = 'tehachapi'
windpark = NREL().get_windpark(NREL.park_id['tehachapi'], 30, 2004)
print "Working on windpark around target turbine", str(windpark.get_target_idx())
print "Plotting windpark ..."
title = "Some Turbines of NREL Data Set"
show_coord_topo(windpark, title)
开发者ID:Bengt,项目名称:windml,代码行数:24,代码来源:show_coord_topo.py
示例19: test_get_windpark
def test_get_windpark(self):
ds = NREL()
windpark = ds.get_windpark(NREL.park_id['tehachapi'], 10, 2004, 2005)
assert(len(windpark.mills) == 66)
开发者ID:tverrbjelke,项目名称:windml,代码行数:4,代码来源:nrel_test.py
示例20: NREL
from windml.datasets.nrel import NREL
from windml.mapping.power_mapping import PowerMapping
from windml.preprocessing.preprocessing import destroy
from windml.preprocessing.preprocessing import interpolate
from windml.visualization.plot_timeseries import plot_timeseries
import matplotlib.pyplot as plt
import matplotlib.dates as md
from pylab import *
from numpy import array, zeros, float32, int32
# get windpark 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()
measurements = target.get_measurements()[300:1000]
damaged, indices = destroy(measurements, method="nmar", percentage=.80,\
min_length=10, max_length=100)
neighbors = windpark.get_turbines()[:-1]
nseries = [t.get_measurements()[300:1000] for t in neighbors]
tinterpolated = interpolate(damaged, method='mreg',\
timestep=600,\
neighbor_series = nseries,\
reg = 'linear_model')
d = array([m[0] for m in tinterpolated])
开发者ID:Bengt,项目名称:windml,代码行数:30,代码来源:mreg_lin_interpolation.py
注:本文中的windml.datasets.nrel.NREL类示例由纯净天空整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论