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)