Question

This is a python matplotlib question. So it would be great if you could show me in python method.

Histogram of wind gust speeds As before the file akaroawindgusts.txt contains hourly maximum wind gusts speeds at the Akaroa Electronic weather station (EW) through February 2014. To plot the wind gusts you extracted column 3. If you extrad column 2 as well in a separate call to np. loadtxt) you'll have a second array giving the direction of each gust in degrees One might expect that the histogram of wind gusts might be different for different directions. To investigate that, you are to write a function plot wind_gust histogram(filename, sin angle, nax_angle) that reads a file like akaroawindgusts.txt and plots a histogram of all the measured wind gusts that lie in the range min angle to max ongle inclusive. The plot resulting from a call plot_wind_gust histogram("akaroawindgusts.txt 180, 270) should be

I have this loadtxt that asked to plot histogram of wind gusts(column 3) that lie in direction angle(column 2) from min angle to max angle inclusively. I don't know how to include min_angle and max_angle into my codes.

Answer 1

import numpy as np
import matplotlib.pyplot as plt
from math import sin,pi
import math

def plot_wind_gust_histogram(filename,min_angle,max_angle):
   data1=np.loadtxt(filename,delimiter=";",skiprows=1,usecols=2) #speed
   data2=np.loadtxt(filename,delimiter=";",skiprows=1,usecols=1) #direction   add this
  
   data1=data1[np.where(np.logical_and( data2[:]>=min_angle , data2[:]<=max_angle ))].astype(np.int)   #add this

   axes=plt.axes()
   axes.hist(data1, bins=range(0,141,10) , facecolor='blue')
   axes.grid(True)
   axes.set_title("Wind gusts from file akaroawindgusts.txt in direction 180 and 270")
   axes.set_xlabel("Wind gust speed (km/h)")
   axes.set_ylabel("Count")
   axes.set_xticks(np.arange(0,141,10))
   axes.set_yticks(np.arange(0,70,10))
   plt.show()

plot_wind_gust_histogram("akaroawindgusts.txt",180,270)

1.py import numpy as np import matplotlib.pyplot as plt from math import sin,pi import math 4 def plot_wind gust_histogram(filename,min_angle,max_angle): datal=np.loadtxt(filename, delimiter=";",skiprows=1, usecols=2) #speed data2=np.loadtxt(filename, delimiter=";", skiprows=1, usecols=1) #direction add this data1-datal[np.where(np.logical_and ( data2[:]>=min_angle , data2[:]<=max_angle ))].astype(np.int) #add this axes=plt.axes() axes.hist(datal, bins=range(0, 141,10), facecolor='blue'), axes.grid(True) axes.set_title("Wind gusts from file akaroawindgusts.txt in direction 180 and 270") axes.set_xlabel("Wind gust speed (km/h)") axes.set ylabel("Count")|| axes.set xticks(np.arange(0, 141,10)) axes.set yticks(np.arange(0,70,10)) plt.show() plot wind gust histogram( "akaroawindgusts.txt", 180,270)