windml.util.features.compute_highlevel_features(turbine)
month_power += np.array(month_power_m)
ramps_up += np.array(ramps_up_m)
ramps_down += np.array(ramps_down_m)
power_freq += np.array(power_freq_m)
month = [
'jan', 'feb', 'mar', 'apr', 'may', 'jun', 'jul', 'aug', 'sep', 'oct',
'nov', 'dec'
]
figure = plt.figure(figsize=(15, 10))
# plot 1
plot1 = plt.subplot(2, 2, 1)
show_coord_topo(windpark, "Park Topology", show=False)
# plot 2
plot2 = plt.subplot(2, 2, 2)
plt.title("Monthly Power")
plot2.bar(range(1, 13), month_power, color='#c4d8eb')
plot2.set_ylabel('Power', color='black')
plot2.set_xlabel('Month', color='black')
plot2.set_xticks([i + 0.5 for i in range(1, 13)])
plot2.set_xticklabels(month)
ramps_up = list(ramps_up)
ramps_down = list(ramps_down)
# plot 3
length = 2 * len(windml.util.features.rampheights) + 1
for turbine in windpark.turbines[1:]:
feat, month_power_m, ramps_up_m, ramps_down_m, power_freq_m =\
windml.util.features.compute_highlevel_features(turbine)
month_power+=np.array(month_power_m)
ramps_up+=np.array(ramps_up_m)
ramps_down+=np.array(ramps_down_m)
power_freq+=np.array(power_freq_m)
month = ['jan','feb','mar','apr','may','jun','jul','aug','sep','oct','nov','dec']
figure = plt.figure(figsize=(15, 10))
# plot 1
plot1 = plt.subplot(2, 2, 1)
show_coord_topo(windpark, "Park Topology", show = False)
# plot 2
plot2 = plt.subplot(2, 2, 2)
plt.title("Monthly Power")
plot2.bar(range(1,13), month_power, color='#c4d8eb')
plot2.set_ylabel('Power', color='black')
plot2.set_xlabel('Month', color='black')
plot2.set_xticks([i+0.5 for i in range(1,13)])
plot2.set_xticklabels(month)
ramps_up = list(ramps_up)
ramps_down = list(ramps_down)
# plot 3
length = 2*len(windml.util.features.rampheights)+1
"""
Location of AEMO Turbines
-------------------------------------------------------------------------
The black dots illustrate the locations of the turbines. The red dot
marks the target turbine.
"""
# Author: Nils A. Treiber <*****@*****.**>
# Author: Jendrik Poloczek <*****@*****.**>
# License: BSD 3 clause
from __future__ import print_function
from windml.datasets.aemo import AEMO
from windml.visualization.show_coord_topo import show_coord_topo
# see http://windfarmperformance.info/ for more information of the park
windpark = AEMO().get_windpark(AEMO.park_id['cathrock'], 3200)
print('Working on windpark around target turbine', str(windpark.get_target_idx()))
print('Plotting windpark ...')
title = 'Turbines of AEMO Data Set'
show_coord_topo(windpark, title)
for windmill in windpark.mills[1:]:
feat, month_power_m, ramps_up_m, ramps_down_m, power_freq_m = windml.util.features.compute_highlevel_features(windmill)
month_power+=np.array(month_power_m)
ramps_up+=np.array(ramps_up_m)
ramps_down+=np.array(ramps_down_m)
power_freq+=np.array(power_freq_m)
month = ['jan','feb','mar','apr','may','jun','jul','aug','sep','oct','nov','dec']
figure = plt.figure(figsize=(15, 10))
# plot 1
plot1 = plt.subplot(2, 2, 1)
plt.title("Park Topology")
show_coord_topo(windpark, show = False)
# plot 2
plot2 = plt.subplot(2, 2, 2)
plt.title("Monthly Power")
plot2.bar(range(1,13), month_power, color='#c4d8eb')
plot2.set_ylabel('Power', color='black')
plot2.set_xlabel('Month', color='black')
plot2.set_xticks([i+0.5 for i in range(1,13)])
plot2.set_xticklabels(month)
ramps_up = list(ramps_up)
ramps_down = list(ramps_down)
# plot 3
length = 2*len(windml.util.features.rampheights)+1
def test_coord_topo(self):
show_coord_topo(self.windpark, show=False)
""" 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)
def test_coord_topo(self):
show_coord_topo(self.windpark, show = False)
def test_coord_topo(self):
show_coord_topo(self.windpark, 'wind park topology', show = False)
