Global & Local Temperature

This example uses historical data since 1880 on average global surface temperature changes from NASA's GISS Surface Temperature Analysis (GISTEMP) (original file) That site has lots of other data sets in csv format, too.

Click the "play" icon or press shift+enter to execute each cell.


In [ ]:
# import the software packages needed
import pandas as pd
import numpy as np
%matplotlib inline
import matplotlib as mpl
import matplotlib.pyplot as plt
inline_rc = dict(mpl.rcParams)
import datetime as dt
import matplotlib.dates as mdates

Importing a data set

This next cell may take a little while to run if it's grabbing a pretty big data set. The cell label to the left will look like "In [*]" while it's still thinking and "In [2]" when it's finished.


In [ ]:
# Combined land and ocean temperature averages (LOTI: Land Ocean Temperature Index)
data1 = pd.read_csv('http://github.com/adamlamee/CODINGinK12-data/raw/master/LOTI.csv', header=1).replace(to_replace="***", value=np.NaN)
data_LOTI = data1.apply(lambda x: pd.to_numeric(x, errors='ignore'))

# Only land temperature averages
data2 = pd.read_csv('http://github.com/adamlamee/CODINGinK12-data/raw/master/LAND.csv', header=1).replace(to_replace="***", value=np.NaN)
data_LAND = data2.apply(lambda x: pd.to_numeric(x, errors='ignore'))

data3 = pd.read_csv('./Orlando Temps.csv')
#data3 = pd.read_csv(https://gist.githubusercontent.com/GoodmanSciences/9953c23049f7814c0d26b5a49a5375f3/raw/72423fbdf5c54e14083d30bbbfccea79553f8025/Orlando%2520Temps.csv)

dates = ['1/1/1975','1/2/1975','01/3/1975']
x = [dt.datetime.strptime(d,'%m/%d/%Y').date() for d in dates]
y = "range(len(x))"

PART 1: Local Temperature Changes & Climate

We can view the first few rows of the file we just imported.


In [ ]:
# The .head(n) command displays the first n rows of the local temperature data file.
data3.head(4)

In [ ]:
data3.columns

In [ ]:
ax = data3.plot('date', 'avg', title="Average Temperature at Orlando International Airport (1975 - Present)", legend=False)
ax.set(xlabel="Date", ylabel="Average Temperature (Fahrenheit)")


# Get current size
fig_size = plt.rcParams["figure.figsize"]
 
# Prints: [8.0, 6.0]
print( "Current size:", fig_size)
 
# Set figure width to 20 and height to 6.5
fig_size[0] = 20
fig_size[1] = 6.5
plt.rcParams["figure.figsize"] = fig_size

plt.show()

PART 2: Global Temperature Changes & Climate


In [ ]:
# The .head(n) command displays the first n rows of the file.
data_LAND.head(5)

In [ ]:
x1 = data_LOTI.Year
y1 = data_LOTI.JanDec

# plt.plot() makes a line graph, by default
fig = plt.figure(figsize=(10, 5))
plt.plot(x1, y1)
plt.title('Average land an docean temperature readings')
plt.xlabel('Year')
plt.ylabel('Percent temp change')

In [ ]:
x2 = data_LAND.Year
y2 = data_LAND.JanDec

# plt.plot() makes a line graph, by default
fig = plt.figure(figsize=(10, 5))
plt.plot(x2, y2)
plt.title('Land temperature readings')
plt.xlabel('Year')
plt.ylabel('Percent temp change')

In [ ]:
# Wow, this needs a title and axis labels!
fig = plt.figure(figsize=(10, 5))
plt.plot(x1, y1, label="Land and Ocean")
plt.plot(x2, y2, label="Land only")


plt.legend()
plt.show()

Try editing some code and re-running the cell to see the effects. For more information on formatting the markdown text in a cell like this one, go to Help > Markdown > Basic Writing and Formatting Text.

To save your work: go to File > Save and Checkpoint. That only saves your edits as long as you're working inthis notebook. All is lost after 10 minutes of inactivity.

To really save your work: go to File > Download as > iPython notebook (or save as pdf if you just want to show someone a snapshot of what your code and output look like.

References