hurricanes_analysis.py

# names of hurricanes
names = ['Cuba I', 'San Felipe II Okeechobee', 'Bahamas', 'Cuba II', 'CubaBrownsville', 'Tampico', 'Labor Day', 'New England', 'Carol', 'Janet', 'Carla', 'Hattie', 'Beulah', 'Camille', 'Edith', 'Anita', 'David', 'Allen', 'Gilbert', 'Hugo', 'Andrew', 'Mitch', 'Isabel', 'Ivan', 'Emily', 'Katrina', 'Rita', 'Wilma', 'Dean', 'Felix', 'Matthew', 'Irma', 'Maria', 'Michael']

# months of hurricanes
months = ['October', 'September', 'September', 'November', 'August', 'September', 'September', 'September', 'September', 'September', 'September', 'October', 'September', 'August', 'September', 'September', 'August', 'August', 'September', 'September', 'August', 'October', 'September', 'September', 'July', 'August', 'September', 'October', 'August', 'September', 'October', 'September', 'September', 'October']

# years of hurricanes
years = [1924, 1928, 1932, 1932, 1933, 1933, 1935, 1938, 1953, 1955, 1961, 1961, 1967, 1969, 1971, 1977, 1979, 1980, 1988, 1989, 1992, 1998, 2003, 2004, 2005, 2005, 2005, 2005, 2007, 2007, 2016, 2017, 2017, 2018]

# maximum sustained winds (mph) of hurricanes
max_sustained_winds = [165, 160, 160, 175, 160, 160, 185, 160, 160, 175, 175, 160, 160, 175, 160, 175, 175, 190, 185, 160, 175, 180, 165, 165, 160, 175, 180, 185, 175, 175, 165, 180, 175, 160]

# areas affected by each hurricane
areas_affected = [['Central America', 'Mexico', 'Cuba', 'Florida', 'The Bahamas'], ['Lesser Antilles', 'The Bahamas', 'United States East Coast', 'Atlantic Canada'], ['The Bahamas', 'Northeastern United States'], ['Lesser Antilles', 'Jamaica', 'Cayman Islands', 'Cuba', 'The Bahamas', 'Bermuda'], ['The Bahamas', 'Cuba', 'Florida', 'Texas', 'Tamaulipas'], ['Jamaica', 'Yucatn Peninsula'], ['The Bahamas', 'Florida', 'Georgia', 'The Carolinas', 'Virginia'], ['Southeastern United States', 'Northeastern United States', 'Southwestern Quebec'], ['Bermuda', 'New England', 'Atlantic Canada'], ['Lesser Antilles', 'Central America'], ['Texas', 'Louisiana', 'Midwestern United States'], ['Central America'], ['The Caribbean', 'Mexico', 'Texas'], ['Cuba', 'United States Gulf Coast'], ['The Caribbean', 'Central America', 'Mexico', 'United States Gulf Coast'], ['Mexico'], ['The Caribbean', 'United States East coast'], ['The Caribbean', 'Yucatn Peninsula', 'Mexico', 'South Texas'], ['Jamaica', 'Venezuela', 'Central America', 'Hispaniola', 'Mexico'], ['The Caribbean', 'United States East Coast'], ['The Bahamas', 'Florida', 'United States Gulf Coast'], ['Central America', 'Yucatn Peninsula', 'South Florida'], ['Greater Antilles', 'Bahamas', 'Eastern United States', 'Ontario'], ['The Caribbean', 'Venezuela', 'United States Gulf Coast'], ['Windward Islands', 'Jamaica', 'Mexico', 'Texas'], ['Bahamas', 'United States Gulf Coast'], ['Cuba', 'United States Gulf Coast'], ['Greater Antilles', 'Central America', 'Florida'], ['The Caribbean', 'Central America'], ['Nicaragua', 'Honduras'], ['Antilles', 'Venezuela', 'Colombia', 'United States East Coast', 'Atlantic Canada'], ['Cape Verde', 'The Caribbean', 'British Virgin Islands', 'U.S. Virgin Islands', 'Cuba', 'Florida'], ['Lesser Antilles', 'Virgin Islands', 'Puerto Rico', 'Dominican Republic', 'Turks and Caicos Islands'], ['Central America', 'United States Gulf Coast (especially Florida Panhandle)']]

# damages (USD($)) of hurricanes
damages = ['Damages not recorded', '100M', 'Damages not recorded', '40M', '27.9M', '5M', 'Damages not recorded', '306M', '2M', '65.8M', '326M', '60.3M', '208M', '1.42B', '25.4M', 'Damages not recorded', '1.54B', '1.24B', '7.1B', '10B', '26.5B', '6.2B', '5.37B', '23.3B', '1.01B', '125B', '12B', '29.4B', '1.76B', '720M', '15.1B', '64.8B', '91.6B', '25.1B']

# deaths for each hurricane
deaths = [90,4000,16,3103,179,184,408,682,5,1023,43,319,688,259,37,11,2068,269,318,107,65,19325,51,124,17,1836,125,87,45,133,603,138,3057,74]

# 1
# Update Recorded Damages
conversion = {"M": 1000000,
              "B": 1000000000}
def updated_damages(damages_lst):
    updated_damage = []
    for damage in damages:
        if "M" in damage:
            damage = damage[:-1]
            damage = float(damage)
            new_damage = 1000000 * damage
            updated_damage.append(new_damage)
        elif "B" in damage:
            damage = damage[:-1]
            damage = float(damage)
            new_damage = 1000000000 * damage
            updated_damage.append(new_damage)
        else:
            new_damage = damage
            updated_damage.append(new_damage)
    return updated_damage


# test function by updating damages
updated_damages = list(updated_damages(damages))


# 2 
# Create a Table

# Create and view the hurricanes dictionary
def hurricane_dictionary(Name, Month, Year, Max_Sustained_Wind, Areas_Affected, Damage, Death):

    dictionary = {}
    for index in range(len(Name)):
        dictionary[Name[index]] = {
            "Name": Name[index],
            "Month": Month[index],
            "Year": Year[index],
            "Max Sustained Wind": Max_Sustained_Wind[index],
            "Areas Affected": Areas_Affected[index],
            "Damages": Damage[index],
            "Deaths": Death[index]
        }
    return dictionary

hurricanes = hurricane_dictionary(names, months, years, max_sustained_winds, areas_affected, updated_damages, deaths)
#print(hurricanes)

# 3
# Organizing by Year

# create a new dictionary of hurricanes with year and key
def hurricane_dictionary_by_year(hurricanes):
    hurricanes_by_year = {}
    for index in hurricanes:
        current_year = hurricanes[index]["Year"]
        current_hurricanes = hurricanes[index]
        if current_year not in hurricanes_by_year:
            hurricanes_by_year[current_year] = [current_hurricanes]
        else:
            hurricanes_by_year[current_year].append(current_hurricanes)
    return hurricanes_by_year

hurricane_dictionary_by_year(hurricanes)

# 4
# Counting Damaged Areas
# create dictionary of areas to store the number of hurricanes involved in
def hurricanes_by_affected_areas(hurricanes):
    hurricanes_by_affected_area = {}
    for index in hurricanes:
        current_area = hurricanes[index]["Areas Affected"]
        #print(current_area)
        current_hurricane = hurricanes[index]["Name"]
        for area in current_area:
            #print("AREA: ", area)
            if area not in hurricanes_by_affected_area:

                hurricanes_by_affected_area[area] = 1
            else:

                hurricanes_by_affected_area[area] +=1

    return hurricanes_by_affected_area

num_hurricanes_by_area = hurricanes_by_affected_areas(hurricanes)

# 5 
# Calculating Maximum Hurricane Count

# find most frequently affected area and the number of hurricanes involved in

new_value = max(num_hurricanes_by_area, key=num_hurricanes_by_area.get)
number = num_hurricanes_by_area[new_value] 

#print("Most_commonly affected area is", new_value, "with", number, "hurricanes")

# 6
# Calculating the Deadliest Hurricane

# find highest mortality hurricane and the number of deaths
def highest_mortality_by_hurricane(hurricanes):
  new_deaths_by_hurricane = 0
  for index in hurricanes:

    hurricane_name = hurricanes[index]["Name"]
    deaths_by_hurricane = hurricanes[index]["Deaths"]
    if deaths_by_hurricane > new_deaths_by_hurricane:
      new_deaths_by_hurricane = deaths_by_hurricane
      causative_hurricane = hurricane_name
  return causative_hurricane, new_deaths_by_hurricane

deadliest_hurricane = highest_mortality_by_hurricane(hurricanes)
print("The hurricane that caused most deaths is ", deadliest_hurricane[0], " having caused ",deadliest_hurricane[1], " deaths" )



# 7
# Rating Hurricanes by Mortality


# categorize hurricanes in new dictionary with mortality severity as key
def hurricane_mortality(hurricanes):
  hurricanes_by_mortality = {0:[],1:[],2:[],3:[],4:[],5:[]}
  for index in hurricanes:

    hurricane_name = hurricanes[index]["Name"]
    deaths_by_hurricane = hurricanes[index]["Deaths"]
    if deaths_by_hurricane > 10000:
      hurricanes_by_mortality[5].append(hurricanes[index]) 
    elif deaths_by_hurricane > 1000:
      hurricanes_by_mortality[4].append(hurricanes[index])
    elif deaths_by_hurricane > 500:
      hurricanes_by_mortality[3].append(hurricanes[index])
    elif deaths_by_hurricane > 100:
      hurricanes_by_mortality[2].append(hurricanes[index])
    elif deaths_by_hurricane > 0:
      hurricanes_by_mortality[1].append(hurricanes[index])
    else:
      hurricanes_by_mortality[0].append(hurricanes[index])
  return hurricanes_by_mortality

hurricane_mortality_risk = hurricane_mortality(hurricanes)
#print(hurricane_mortality_risk[5])
  
 

# 8 Calculating Hurricane Maximum Damage

# find highest damage inducing hurricane and its total cost
def max_damage (hurricanes):
  maximum_damage = 0
  for index in hurricanes:

    hurricane_name = hurricanes[index]["Name"]
    damages_by_hurricane = hurricanes[index]["Damages"]
    if damages_by_hurricane == "Damages not recorded":
      continue
    else:
      damages_by_hurricance = int(damages_by_hurricane)
      if damages_by_hurricane > maximum_damage:
        maximum_damage = damages_by_hurricane
        causative_hurricane = hurricane_name
  return maximum_damage, causative_hurricane

damaging_hurricane = max_damage(hurricanes)
maximum_damage = damaging_hurricane[0]
causative_hurricane = damaging_hurricane[1]
print("The hurricane that caused most damage is ", causative_hurricane, " having caused ", maximum_damage," dollars worth of damage" )




# 9
# Rating Hurricanes by Damage
damage_scale = {0: 0,
                1: 100000000,
                2: 1000000000,
                3: 10000000000,
                4: 50000000000}

 
# categorize hurricanes in new dictionary with damage severity as key
def hurricane_damage_rating(hurricanes):
  hurricanes_by_damages = {0:[],1:[],2:[],3:[],4:[],5:[]}
  for index in hurricanes:
    name_hurricane = hurricanes[index]["Name"]
    damages_by_hurricane = hurricanes[index]["Damages"]
    if damages_by_hurricane == "Damages not recorded":
      hurricanes_by_damages[0].append(hurricanes[index])
    else:
      damages_by_hurricane = int(damages_by_hurricane)
      if damages_by_hurricane > 50000000000:
        hurricanes_by_damages[5].append(hurricanes[index])
      elif damages_by_hurricane > 10000000000:
        hurricanes_by_damages[4].append(hurricanes[index])
      elif damages_by_hurricane > 1000000000:
        hurricanes_by_damages[3].append(hurricanes[index])
      elif damages_by_hurricane > 100000000:
        hurricanes_by_damages[2].append(hurricanes[index])
      else:
        hurricanes_by_damages[2].append(hurricanes[index])
  return hurricanes_by_damages
      
hurricane_ratings_by_damage = hurricane_damage_rating(hurricanes)
#print(hurricane_ratings_by_damage[0])