{
"metadata": {
"name": "Amazon river discharge"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "heading",
"level": 1,
"metadata": {},
"source": [
"Amazon river discharge"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"[Nikolay Koldunov](http://koldunov.net)\n",
" "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"I want to show a small example of the work flow that is more or less typical during research process. Often you see some interesting feature in your data and want to investigate it in more detail. If you not lucky enough to work with the model data, this would require dealing with multiple data sources, and possibly multiple file formats. Having all data sets in one place in consistent format becomes very handy for this type of applications."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"I begin with the visualization of river discharge data, that I did for one of my projects. It is performed in [Tableau](http://www.tableausoftware.com/public), and use [Dai and Trenberth Global River Flow and Continental Discharge Dataset](http://www.cgd.ucar.edu/cas/catalog/surface/dai-runoff/index.html). This data set is also available on Marinexplore. "
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from IPython.display import HTML\n",
"HTML('')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
""
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 7,
"text": [
""
]
}
],
"prompt_number": 7
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"This visualization is interactive, so if you choose several rivers on the map it will show a bar chart with distribution of river discharges. "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"I also want to show some simple statistics about my data, so I create another visualization with only 20 biggest rivers:"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"HTML('')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
""
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 8,
"text": [
""
]
}
],
"prompt_number": 8
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Look at Amazon. I knew that it's a biggest river in the world, but I never though it's domination is so overwhelming. This huge river discharge should be visible somehow in salinity field of the ocean. In order to test this I decide to look at satellite salinity data."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"But first we would need to plot our river discharge data in pyhton. [Text file](http://www.cgd.ucar.edu/cas/catalog/surface/dai-runoff/coastal-stns-byVol-updated-oct2007.txt) with anual mean river discharge values from Dai and Trenberth have been converted in .csv format with vertical bar as separator and is available [here](https://dl.dropboxusercontent.com/u/107639/notebooks/amazon/rivers2.csv). We are going to open it with [pandas](http://pandas.pydata.org/), and plot with basemap. "
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"%pylab inline"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"\n",
"Welcome to pylab, a matplotlib-based Python environment [backend: module://IPython.kernel.zmq.pylab.backend_inline].\n",
"For more information, type 'help(pylab)'.\n"
]
}
],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import pandas as pd"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 44
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"df = pd.read_csv('rivers2.csv', sep='|', na_values=[' Null', ' Null'])"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 3
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Here is our data in pandas DataFrame:"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"df.head()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" No | \n",
" m2s_ra | \n",
" lonm | \n",
" latm | \n",
" area(km2) | \n",
" Vol | \n",
" nyr | \n",
" yrb | \n",
" yre | \n",
" elev(m) | \n",
" CT | \n",
" CN | \n",
" River_Name | \n",
" OCN | \n",
" Station_Name | \n",
"
\n",
" \n",
" \n",
" \n",
" | 0 | \n",
" 1 | \n",
" 12462 | \n",
" -51.75 | \n",
" -0.75 | \n",
" 4618750 | \n",
" 5389.537 | \n",
" 79.0 | \n",
" 1928 | \n",
" 2006 | \n",
" 37 | \n",
" BR | \n",
" SA | \n",
" Amazon | \n",
" ATL | \n",
" Obidos, Bra | \n",
"
\n",
" \n",
" | 1 | \n",
" 2 | \n",
" 10291 | \n",
" 12.75 | \n",
" -5.75 | \n",
" 3475000 | \n",
" 1270.203 | \n",
" 97.1 | \n",
" 1903 | \n",
" 2000 | \n",
" 270 | \n",
" CD | \n",
" AF | \n",
" Congo | \n",
" ATL | \n",
" Kinshasa, C | \n",
"
\n",
" \n",
" | 2 | \n",
" 3 | \n",
" 11474 | \n",
" -61.75 | \n",
" 9.25 | \n",
" 836000 | \n",
" 980.088 | \n",
" 75.8 | \n",
" 1923 | \n",
" 1999 | \n",
" -9999 | \n",
" VE | \n",
" SA | \n",
" Orinoco | \n",
" ATL | \n",
" Pte Angostu | \n",
"
\n",
" \n",
" | 3 | \n",
" 4 | \n",
" 10374 | \n",
" 120.75 | \n",
" 32.25 | \n",
" 1705383 | \n",
" 905.141 | \n",
" 99.6 | \n",
" 1900 | \n",
" 2000 | \n",
" 19 | \n",
" CN | \n",
" AS | \n",
" Changjiang | \n",
" PAC | \n",
" Datong, Chi | \n",
"
\n",
" \n",
" | 4 | \n",
" 5 | \n",
" 10245 | \n",
" 89.75 | \n",
" 24.25 | \n",
" 554542 | \n",
" 670.943 | \n",
" 44.3 | \n",
" 1956 | \n",
" 2000 | \n",
" 19 | \n",
" BD | \n",
" AS | \n",
" Brahmaputra | \n",
" IND | \n",
" Bahadurabad | \n",
"
\n",
" \n",
"
\n",
"