{
"cells": [
{
"cell_type": "markdown",
"metadata": {
"editable": true,
"slideshow": {
"slide_type": ""
},
"tags": []
},
"source": [
"# tstoolbox: Statistics\n",
"I will use the water level record for USGS site '02240000 OCKLAWAHA RIVER NEAR CONNER, FL'. Use tsgettoolbox to download the data we will be using."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"from tsgettoolbox import tsgettoolbox\n",
"from tstoolbox import tstoolbox"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"ock_flow = tsgettoolbox.nwis(\n",
" sites=\"02240000\", startDT=\"2008-01-01\", endDT=\"2016-01-01\", parameterCd=\"00060\"\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"
\n",
" \n",
" \n",
" | \n",
" USGS-02240000-00060 | \n",
"
\n",
" \n",
" Datetime | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" 2008-01-01 | \n",
" 541 | \n",
"
\n",
" \n",
" 2008-01-02 | \n",
" 533 | \n",
"
\n",
" \n",
" 2008-01-03 | \n",
" 520 | \n",
"
\n",
" \n",
" 2008-01-04 | \n",
" 508 | \n",
"
\n",
" \n",
" 2008-01-05 | \n",
" 516 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" USGS-02240000-00060\n",
"Datetime \n",
"2008-01-01 541\n",
"2008-01-02 533\n",
"2008-01-03 520\n",
"2008-01-04 508\n",
"2008-01-05 516"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ock_flow.head()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [],
"source": [
"ock_stage = tsgettoolbox.nwis(\n",
" sites=\"02240000\", startDT=\"2008-01-01\", endDT=\"2016-01-01\", parameterCd=\"00065\"\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" USGS-02240000-00065 | \n",
"
\n",
" \n",
" Datetime | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" 2008-01-01 | \n",
" 3.67 | \n",
"
\n",
" \n",
" 2008-01-02 | \n",
" 3.64 | \n",
"
\n",
" \n",
" 2008-01-03 | \n",
" 3.60 | \n",
"
\n",
" \n",
" 2008-01-04 | \n",
" 3.58 | \n",
"
\n",
" \n",
" 2008-01-05 | \n",
" 3.57 | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" USGS-02240000-00065\n",
"Datetime \n",
"2008-01-01 3.67\n",
"2008-01-02 3.64\n",
"2008-01-03 3.60\n",
"2008-01-04 3.58\n",
"2008-01-05 3.57"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ock_stage.head()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The tstoolbox.rolling_window calculates the average stage over the spans listed in the 'span' argument."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"r_ock_stage = tstoolbox.rolling_window(\n",
" input_ts=ock_stage, span=\"1,14,30,60,90,120,180,274,365\", statistic=\"mean\"\n",
")\n",
"a_ock_stage = tstoolbox.aggregate(\n",
" input_ts=r_ock_stage, agg_interval=\"A\", statistic=\"max\"\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" USGS-02240000-00065_1_mean_max | \n",
" USGS-02240000-00065_14_mean_max | \n",
" USGS-02240000-00065_30_mean_max | \n",
" USGS-02240000-00065_60_mean_max | \n",
" USGS-02240000-00065_90_mean_max | \n",
" USGS-02240000-00065_120_mean_max | \n",
" USGS-02240000-00065_180_mean_max | \n",
" USGS-02240000-00065_274_mean_max | \n",
" USGS-02240000-00065_365_mean_max | \n",
"
\n",
" \n",
" Datetime | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" 2008-12-31 | \n",
" 7.73 | \n",
" 6.055000 | \n",
" 5.638333 | \n",
" 4.383333 | \n",
" 4.049444 | \n",
" 3.958000 | \n",
" 3.985111 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 2009-12-31 | \n",
" 5.29 | \n",
" 4.949286 | \n",
" 4.257667 | \n",
" 4.179167 | \n",
" 4.272111 | \n",
" 4.154667 | \n",
" 3.953278 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 2010-12-31 | \n",
" 7.49 | \n",
" 6.492857 | \n",
" 5.919667 | \n",
" 5.079167 | \n",
" 4.780000 | \n",
" 4.788167 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 2011-12-31 | \n",
" 6.22 | \n",
" 5.357143 | \n",
" 5.227000 | \n",
" 5.189333 | \n",
" 5.175444 | \n",
" 5.147917 | \n",
" 4.967722 | \n",
" 4.632701 | \n",
" 4.345342 | \n",
"
\n",
" \n",
" 2012-12-31 | \n",
" 6.12 | \n",
" 5.473571 | \n",
" 5.334667 | \n",
" 5.308000 | \n",
" 5.233444 | \n",
" 5.122083 | \n",
" 4.809944 | \n",
" 4.554051 | \n",
" 4.559014 | \n",
"
\n",
" \n",
" 2013-12-31 | \n",
" 5.39 | \n",
" 5.248571 | \n",
" 5.083000 | \n",
" 4.775333 | \n",
" 4.458444 | \n",
" 4.192667 | \n",
" 3.812778 | \n",
" 3.410584 | \n",
" NaN | \n",
"
\n",
" \n",
" 2014-12-31 | \n",
" 5.86 | \n",
" 5.762857 | \n",
" 5.704667 | \n",
" 5.614333 | \n",
" 5.510222 | \n",
" 5.437917 | \n",
" 5.213722 | \n",
" 3.709489 | \n",
" 3.379014 | \n",
"
\n",
" \n",
" 2015-12-31 | \n",
" 6.85 | \n",
" 5.256429 | \n",
" 5.001333 | \n",
" 4.884167 | \n",
" 4.792333 | \n",
" 4.882583 | \n",
" 5.126222 | \n",
" 5.060730 | \n",
" 4.767342 | \n",
"
\n",
" \n",
" 2016-12-31 | \n",
" 3.41 | \n",
" 3.485000 | \n",
" 3.549667 | \n",
" 3.667167 | \n",
" 3.866889 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" USGS-02240000-00065_1_mean_max USGS-02240000-00065_14_mean_max \\\n",
"Datetime \n",
"2008-12-31 7.73 6.055000 \n",
"2009-12-31 5.29 4.949286 \n",
"2010-12-31 7.49 6.492857 \n",
"2011-12-31 6.22 5.357143 \n",
"2012-12-31 6.12 5.473571 \n",
"2013-12-31 5.39 5.248571 \n",
"2014-12-31 5.86 5.762857 \n",
"2015-12-31 6.85 5.256429 \n",
"2016-12-31 3.41 3.485000 \n",
"\n",
" USGS-02240000-00065_30_mean_max USGS-02240000-00065_60_mean_max \\\n",
"Datetime \n",
"2008-12-31 5.638333 4.383333 \n",
"2009-12-31 4.257667 4.179167 \n",
"2010-12-31 5.919667 5.079167 \n",
"2011-12-31 5.227000 5.189333 \n",
"2012-12-31 5.334667 5.308000 \n",
"2013-12-31 5.083000 4.775333 \n",
"2014-12-31 5.704667 5.614333 \n",
"2015-12-31 5.001333 4.884167 \n",
"2016-12-31 3.549667 3.667167 \n",
"\n",
" USGS-02240000-00065_90_mean_max USGS-02240000-00065_120_mean_max \\\n",
"Datetime \n",
"2008-12-31 4.049444 3.958000 \n",
"2009-12-31 4.272111 4.154667 \n",
"2010-12-31 4.780000 4.788167 \n",
"2011-12-31 5.175444 5.147917 \n",
"2012-12-31 5.233444 5.122083 \n",
"2013-12-31 4.458444 4.192667 \n",
"2014-12-31 5.510222 5.437917 \n",
"2015-12-31 4.792333 4.882583 \n",
"2016-12-31 3.866889 NaN \n",
"\n",
" USGS-02240000-00065_180_mean_max \\\n",
"Datetime \n",
"2008-12-31 3.985111 \n",
"2009-12-31 3.953278 \n",
"2010-12-31 NaN \n",
"2011-12-31 4.967722 \n",
"2012-12-31 4.809944 \n",
"2013-12-31 3.812778 \n",
"2014-12-31 5.213722 \n",
"2015-12-31 5.126222 \n",
"2016-12-31 NaN \n",
"\n",
" USGS-02240000-00065_274_mean_max USGS-02240000-00065_365_mean_max \n",
"Datetime \n",
"2008-12-31 NaN NaN \n",
"2009-12-31 NaN NaN \n",
"2010-12-31 NaN NaN \n",
"2011-12-31 4.632701 4.345342 \n",
"2012-12-31 4.554051 4.559014 \n",
"2013-12-31 3.410584 NaN \n",
"2014-12-31 3.709489 3.379014 \n",
"2015-12-31 5.060730 4.767342 \n",
"2016-12-31 NaN NaN "
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"a_ock_stage"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"scrolled": true
},
"outputs": [],
"source": [
"q_max_avg_obs = tstoolbox.calculate_fdc(input_ts=a_ock_stage, sort_order=\"descending\")"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"
\n",
" \n",
" \n",
" | \n",
" USGS-02240000-00065_1_mean_max | \n",
" USGS-02240000-00065_14_mean_max | \n",
" USGS-02240000-00065_30_mean_max | \n",
" USGS-02240000-00065_60_mean_max | \n",
" USGS-02240000-00065_90_mean_max | \n",
" USGS-02240000-00065_120_mean_max | \n",
" USGS-02240000-00065_180_mean_max | \n",
" USGS-02240000-00065_274_mean_max | \n",
" USGS-02240000-00065_365_mean_max | \n",
"
\n",
" \n",
" Plotting_position | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
" | \n",
"
\n",
" \n",
" \n",
" \n",
" 10.000000 | \n",
" 3.41 | \n",
" 3.485000 | \n",
" 3.549667 | \n",
" 3.667167 | \n",
" 3.866889 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 11.111111 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.958000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 12.500000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.812778 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 16.666667 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.410584 | \n",
" NaN | \n",
"
\n",
" \n",
" 20.000000 | \n",
" 5.29 | \n",
" 4.949286 | \n",
" 4.257667 | \n",
" 4.179167 | \n",
" 4.049444 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.379014 | \n",
"
\n",
" \n",
" 22.222222 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.154667 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 25.000000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.953278 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 30.000000 | \n",
" 5.39 | \n",
" 5.248571 | \n",
" 5.001333 | \n",
" 4.383333 | \n",
" 4.272111 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 33.333333 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.192667 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 33.333333 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.709489 | \n",
" NaN | \n",
"
\n",
" \n",
" 37.500000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 3.985111 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 40.000000 | \n",
" 5.86 | \n",
" 5.256429 | \n",
" 5.083000 | \n",
" 4.775333 | \n",
" 4.458444 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.345342 | \n",
"
\n",
" \n",
" 44.444444 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.788167 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 50.000000 | \n",
" 6.12 | \n",
" 5.357143 | \n",
" 5.227000 | \n",
" 4.884167 | \n",
" 4.780000 | \n",
" NaN | \n",
" 4.809944 | \n",
" 4.554051 | \n",
" NaN | \n",
"
\n",
" \n",
" 55.555556 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.882583 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 60.000000 | \n",
" 6.22 | \n",
" 5.473571 | \n",
" 5.334667 | \n",
" 5.079167 | \n",
" 4.792333 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 60.000000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.559014 | \n",
"
\n",
" \n",
" 62.500000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.967722 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 66.666667 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.122083 | \n",
" NaN | \n",
" 4.632701 | \n",
" NaN | \n",
"
\n",
" \n",
" 70.000000 | \n",
" 6.85 | \n",
" 5.762857 | \n",
" 5.638333 | \n",
" 5.189333 | \n",
" 5.175444 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 75.000000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.126222 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 77.777778 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.147917 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 80.000000 | \n",
" 7.49 | \n",
" 6.055000 | \n",
" 5.704667 | \n",
" 5.308000 | \n",
" 5.233444 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 4.767342 | \n",
"
\n",
" \n",
" 83.333333 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.060730 | \n",
" NaN | \n",
"
\n",
" \n",
" 87.500000 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.213722 | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 88.888889 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" 5.437917 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
" 90.000000 | \n",
" 7.73 | \n",
" 6.492857 | \n",
" 5.919667 | \n",
" 5.614333 | \n",
" 5.510222 | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
" NaN | \n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" USGS-02240000-00065_1_mean_max \\\n",
"Plotting_position \n",
"10.000000 3.41 \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 5.29 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 5.39 \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 5.86 \n",
"44.444444 NaN \n",
"50.000000 6.12 \n",
"55.555556 NaN \n",
"60.000000 6.22 \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 6.85 \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 7.49 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 7.73 \n",
"\n",
" USGS-02240000-00065_14_mean_max \\\n",
"Plotting_position \n",
"10.000000 3.485000 \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 4.949286 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 5.248571 \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 5.256429 \n",
"44.444444 NaN \n",
"50.000000 5.357143 \n",
"55.555556 NaN \n",
"60.000000 5.473571 \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 5.762857 \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 6.055000 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 6.492857 \n",
"\n",
" USGS-02240000-00065_30_mean_max \\\n",
"Plotting_position \n",
"10.000000 3.549667 \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 4.257667 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 5.001333 \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 5.083000 \n",
"44.444444 NaN \n",
"50.000000 5.227000 \n",
"55.555556 NaN \n",
"60.000000 5.334667 \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 5.638333 \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 5.704667 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 5.919667 \n",
"\n",
" USGS-02240000-00065_60_mean_max \\\n",
"Plotting_position \n",
"10.000000 3.667167 \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 4.179167 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 4.383333 \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 4.775333 \n",
"44.444444 NaN \n",
"50.000000 4.884167 \n",
"55.555556 NaN \n",
"60.000000 5.079167 \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 5.189333 \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 5.308000 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 5.614333 \n",
"\n",
" USGS-02240000-00065_90_mean_max \\\n",
"Plotting_position \n",
"10.000000 3.866889 \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 4.049444 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 4.272111 \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 4.458444 \n",
"44.444444 NaN \n",
"50.000000 4.780000 \n",
"55.555556 NaN \n",
"60.000000 4.792333 \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 5.175444 \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 5.233444 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 5.510222 \n",
"\n",
" USGS-02240000-00065_120_mean_max \\\n",
"Plotting_position \n",
"10.000000 NaN \n",
"11.111111 3.958000 \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 NaN \n",
"22.222222 4.154667 \n",
"25.000000 NaN \n",
"30.000000 NaN \n",
"33.333333 4.192667 \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 NaN \n",
"44.444444 4.788167 \n",
"50.000000 NaN \n",
"55.555556 4.882583 \n",
"60.000000 NaN \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 5.122083 \n",
"70.000000 NaN \n",
"75.000000 NaN \n",
"77.777778 5.147917 \n",
"80.000000 NaN \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 5.437917 \n",
"90.000000 NaN \n",
"\n",
" USGS-02240000-00065_180_mean_max \\\n",
"Plotting_position \n",
"10.000000 NaN \n",
"11.111111 NaN \n",
"12.500000 3.812778 \n",
"16.666667 NaN \n",
"20.000000 NaN \n",
"22.222222 NaN \n",
"25.000000 3.953278 \n",
"30.000000 NaN \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 3.985111 \n",
"40.000000 NaN \n",
"44.444444 NaN \n",
"50.000000 4.809944 \n",
"55.555556 NaN \n",
"60.000000 NaN \n",
"60.000000 NaN \n",
"62.500000 4.967722 \n",
"66.666667 NaN \n",
"70.000000 NaN \n",
"75.000000 5.126222 \n",
"77.777778 NaN \n",
"80.000000 NaN \n",
"83.333333 NaN \n",
"87.500000 5.213722 \n",
"88.888889 NaN \n",
"90.000000 NaN \n",
"\n",
" USGS-02240000-00065_274_mean_max \\\n",
"Plotting_position \n",
"10.000000 NaN \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 3.410584 \n",
"20.000000 NaN \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 NaN \n",
"33.333333 NaN \n",
"33.333333 3.709489 \n",
"37.500000 NaN \n",
"40.000000 NaN \n",
"44.444444 NaN \n",
"50.000000 4.554051 \n",
"55.555556 NaN \n",
"60.000000 NaN \n",
"60.000000 NaN \n",
"62.500000 NaN \n",
"66.666667 4.632701 \n",
"70.000000 NaN \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 NaN \n",
"83.333333 5.060730 \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 NaN \n",
"\n",
" USGS-02240000-00065_365_mean_max \n",
"Plotting_position \n",
"10.000000 NaN \n",
"11.111111 NaN \n",
"12.500000 NaN \n",
"16.666667 NaN \n",
"20.000000 3.379014 \n",
"22.222222 NaN \n",
"25.000000 NaN \n",
"30.000000 NaN \n",
"33.333333 NaN \n",
"33.333333 NaN \n",
"37.500000 NaN \n",
"40.000000 4.345342 \n",
"44.444444 NaN \n",
"50.000000 NaN \n",
"55.555556 NaN \n",
"60.000000 NaN \n",
"60.000000 4.559014 \n",
"62.500000 NaN \n",
"66.666667 NaN \n",
"70.000000 NaN \n",
"75.000000 NaN \n",
"77.777778 NaN \n",
"80.000000 4.767342 \n",
"83.333333 NaN \n",
"87.500000 NaN \n",
"88.888889 NaN \n",
"90.000000 NaN "
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"q_max_avg_obs"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"tstoolbox.plot(\n",
" input_ts=q_max_avg_obs,\n",
" ofilename=\"q_max_avg.png\",\n",
" type=\"norm_xaxis\",\n",
" xtitle=\"Annual Exceedence Probability (%)\",\n",
" ytitle=\"Stage (feet)\",\n",
" title=\"N-day Maximum Rolling Window Average Stage for 0224000\",\n",
" legend_names=\"1,14,30,60,90,120,180,274,365\",\n",
")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"![](q_max_avg.png)"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [],
"source": [
"combined = ock_flow.join(ock_stage)\n",
"tstoolbox.plot(\n",
" input_ts=combined,\n",
" type=\"xy\",\n",
" ofilename=\"stage_flow.png\",\n",
" legend=False,\n",
" title=\"Plot of Daily Stage vs Flow for 02240000\",\n",
" xtitle=\"Flow (cfs)\",\n",
" ytitle=\"Stage (ft)\",\n",
")"
]
},
{
"cell_type": "markdown",
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
}
},
"source": [
"![](stage_flow.png)"
]
},
{
"cell_type": "markdown",
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
}
},
"source": [
"Note the signficant effect of tail-water elevation."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true,
"jupyter": {
"outputs_hidden": true
}
},
"outputs": [],
"source": []
}
],
"metadata": {
"jupytext": {
"formats": "ipynb,py:percent"
},
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.11.7"
}
},
"nbformat": 4,
"nbformat_minor": 4
}