|
| 1 | +--- |
| 2 | +jupyter: |
| 3 | + jupytext: |
| 4 | + formats: ipynb,md |
| 5 | + text_representation: |
| 6 | + extension: .md |
| 7 | + format_name: markdown |
| 8 | + format_version: '1.3' |
| 9 | + jupytext_version: 1.13.0 |
| 10 | + kernelspec: |
| 11 | + display_name: Python 3 (ipykernel) |
| 12 | + language: python |
| 13 | + name: python3 |
| 14 | +--- |
| 15 | + |
| 16 | +# Matplotlib interactive examples |
| 17 | + |
| 18 | +This is a small example of using Matplotlib for interactive data visualization. |
| 19 | + |
| 20 | +To get interactive visualization you must "run" the notebook cells (in order). You can hit the "play" triangle above, or the "play-all" double triangle to run all the cells, or select a cell and hit "Shift-Enter". Once active, you can zoom and pan in the provided data using the GUI toolbar at left of each figure. |
| 21 | + |
| 22 | +```python |
| 23 | +import numpy as np |
| 24 | +import matplotlib.pyplot as plt |
| 25 | + |
| 26 | +# this sets up the Matplotlib interactive windows: |
| 27 | +%matplotlib widget |
| 28 | + |
| 29 | +# this changes the default date converter for better interactive plotting of dates: |
| 30 | +plt.rcParams['date.converter'] = 'concise' |
| 31 | +``` |
| 32 | + |
| 33 | +## Time series example: temperature at Dinosaur, Co |
| 34 | + |
| 35 | +Data is from https://www.ncei.noaa.gov/pub/data/uscrn/products/hourly02/ and is hourly temperature and solar radiation at Dinosaur, Colorado. |
| 36 | + |
| 37 | +### Load the data |
| 38 | + |
| 39 | +```python |
| 40 | +# Note the use of datetimes in the file complicate loading a bit. |
| 41 | +# We recommend using pandas or xarray for more elegant solutions |
| 42 | +# to handling complex timeseries data. |
| 43 | +with open('data/small.txt', 'r') as f: |
| 44 | + data = np.genfromtxt(f, dtype='datetime64[s],f,f,f', |
| 45 | + names=['date', 'doy', 'temp', 'solar']) |
| 46 | +datetime = data['date'] |
| 47 | +dayofyear = data['doy'] |
| 48 | +temperature = data['temp'] |
| 49 | +solar = data['solar'] |
| 50 | + |
| 51 | +# make two-day smoothed versions: |
| 52 | +temp_low = np.convolve(temperature, np.ones(48)/48, mode='same') |
| 53 | +solar_low = np.convolve(solar, np.ones(48)/48, mode='same') |
| 54 | +``` |
| 55 | + |
| 56 | +### Plot timeseries |
| 57 | + |
| 58 | +Here we plot the time series, add labels, title, and a small annotation about where the data came from. Using the toolbar on the left you can zoom in on the data to see, for instance, the day/night warming and cooling. You can also pan along the time series using the icon with four arrows. Note that because we have said `sharex=True` when creating the axes, the time window remains connected in both time seroes. |
| 59 | + |
| 60 | +```python |
| 61 | +fig, (ax0, ax1) = plt.subplots(2, 1, sharex=True, constrained_layout=True) |
| 62 | + |
| 63 | +# temperature: |
| 64 | +ax0.plot(datetime, temperature, label='hourly') |
| 65 | +ax0.plot(datetime, temp_low, label='smoothed') |
| 66 | +ax0.legend(loc='upper right') |
| 67 | +ax0.set_ylabel('Temperature $[^oC]$') # note the use of TeX math formatting |
| 68 | + |
| 69 | +# solar-radiation: |
| 70 | +ax1.plot(datetime, solar, label='hourly') |
| 71 | +ax1.plot(datetime, solar_low, label='smoothed') |
| 72 | +ax1.legend(loc='upper right') |
| 73 | +ax1.set_ylabel('Solar radiation $[W\,m^{-2}]$') # note the use of TeX math formatting |
| 74 | + |
| 75 | +ax0.set_title('Observations: Dinosaur, Colorado', loc='left') |
| 76 | +ax0.text(0.03, 0.03, 'https://www.ncei.noaa.gov/pub/data/uscrn/products/hourly02/', |
| 77 | + fontsize='small', fontstyle='italic', transform=ax0.transAxes); |
| 78 | +``` |
| 79 | + |
| 80 | +### Plot property-versus-property |
| 81 | + |
| 82 | +Lets make a scatter plot of solar radiation versus temperature, colored by day of the year. Note that we do not share the axes this time because the dynamic range of the hourly and smoothed data is different. |
| 83 | + |
| 84 | +```python |
| 85 | +fig, [ax0, ax1] = plt.subplots(1, 2, constrained_layout=True, figsize=(8, 4)) |
| 86 | + |
| 87 | +ax0.set_facecolor('0.6') |
| 88 | +sc = ax0.scatter(solar, temperature, c=dayofyear, s=2, cmap='RdBu_r') |
| 89 | +ax0.set_xlabel('Solar Radiation $[W\,m^{-2}]$') |
| 90 | +ax0.set_ylabel('Temperature $[^oC]$') |
| 91 | +ax0.set_title('Hourly', loc='left') |
| 92 | + |
| 93 | + |
| 94 | +ax1.set_facecolor('0.6') |
| 95 | +sc = ax1.scatter(solar_low[::24], temp_low[::24], c=dayofyear[::24], s=4, cmap='RdBu_r') |
| 96 | +ax1.set_xlabel('Solar Radiation $[W\,m^{-2}]$') |
| 97 | +ax1.set_ylabel('Temperature $[^oC]$') |
| 98 | +ax1.set_title('Smoothed', loc='left') |
| 99 | + |
| 100 | +fig.colorbar(sc, ax=[ax0, ax1], shrink=0.6, label='day of year') |
| 101 | +``` |
| 102 | + |
| 103 | +## Image data: global land-surface temperature |
| 104 | + |
| 105 | +Matplotlib can also plot two-dimensional data as in this example of land surface temperature. Note that for geographic data a package like cartopy or metpy is recommended. |
| 106 | + |
| 107 | +### Load the data |
| 108 | + |
| 109 | +```python |
| 110 | +dat = np.genfromtxt('data/MOD_LSTD_E_2021-08-29_gs_720x360.CSV', delimiter=',') |
| 111 | +dat = np.where(dat<1000, dat, np.NaN) |
| 112 | +dat = dat[::-1, :] |
| 113 | +lon = np.arange(-180.0, 180.1, 0.5) |
| 114 | +lat = np.arange(-90.0, 90.1, 0.5) |
| 115 | +date = '2021-08-29 to 2021-09-05' |
| 116 | +source = 'https://neo.sci.gsfc.nasa.gov/view.php?datasetId=MOD_LSTD_E&date=2021-09-01' |
| 117 | +``` |
| 118 | + |
| 119 | +### Plot data |
| 120 | + |
| 121 | +```python |
| 122 | +fig, ax = plt.subplots(constrained_layout=True, figsize=(7, 4)) |
| 123 | +ax.set_facecolor('0.8') |
| 124 | +pc = ax.pcolormesh(lon, lat, dat, shading='auto', cmap='inferno') |
| 125 | +ax.set_aspect(1.3) |
| 126 | +ax.set_xlabel('Longitude $[^o E]$') |
| 127 | +ax.set_ylabel('Latitude $[^o N]$') |
| 128 | +fig.colorbar(pc, shrink=0.6, extend='both', label='land temperature $[^oC]$'); |
| 129 | +ax.set_title(f'source: {source}', loc='left', fontsize='x-small'); |
| 130 | +``` |
| 131 | + |
| 132 | +```python |
| 133 | + |
| 134 | +``` |
0 commit comments