pcc_cxk_stat.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "_\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170519110333 0 neu 1.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170519*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170519-S102835-E120109.018308.V05A.HDF5'])\n",
      "New enigma\n",
      "/automount/radar/scans/2017/2017-05/2017-05-19/n_ppi_010deg/*20170519110333*.h5\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-05/2017-05-19/n_ppi_010deg/n_ppi_010deg_12345_20170519110333_00.h5'])\n",
      "________Beam Blockage______\n",
      "67.609375\n",
      "________CLUTTER______\n",
      "67.609375\n",
      "________offset______\n",
      "66.609375\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-05/2017-05-19/raa01-rx_10000-1705191105-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "75\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n",
      "satlib.py:895: RuntimeWarning: Mean of empty slice\n",
      "  gr_ipoli_data[i] = np.nanmean(gr_data[rr < radius])\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/numpy/lib/nanfunctions.py:1434: RuntimeWarning: Degrees of freedom <= 0 for slice.\n",
      "  keepdims=keepdims)\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/numpy/lib/function_base.py:4033: RuntimeWarning: All-NaN slice encountered\n",
      "  r = func(a, **kwargs)\n",
      "satlib.py:898: RuntimeWarning: All-NaN slice encountered\n",
      "  gr_ipoli_max[i] = np.nanmax(gr_data[rr < radius])\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 39.120423)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 128.52925399999998)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20150330233003 1 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20150330*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20150330-S055056-E072328.006155.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20150330-S224903-E002135.006166.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2015/2015-03/2015-03-30/n_ppi_045deg/2015-03-30--23:30:03,00.mvol'])\n",
      "________Beam Blockage______\n",
      "43.789062\n",
      "________CLUTTER______\n",
      "42.22461\n",
      "________offset______\n",
      "40.22461\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2015/2015-03/2015-03-30/raa01-rx_10000-1503302330-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 57.51465300000001)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 121.00957)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20151015203657 1 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151015*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151015-S115821-E133055.009255.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151015-S194118-E211352.009260.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2015/2015-10/2015-10-15/n_ppi_045deg/2015-10-15--20:36:57,00.mvol'])\n",
      "________Beam Blockage______\n",
      "50.00586\n",
      "________CLUTTER______\n",
      "49.566406\n",
      "________offset______\n",
      "47.566406\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2015/2015-10/2015-10-15/raa01-rx_10000-1510152035-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 41.60774199999997)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 113.74102699999997)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170909021332 0 neu 1.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170909*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170909-S013353-E030627.020060.V05A.HDF5'])\n",
      "New enigma\n",
      "/automount/radar/scans/2017/2017-09/2017-09-09/n_ppi_010deg/*20170909021332*.h5\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-09/2017-09-09/n_ppi_010deg/n_ppi_010deg_12345_20170909021332_00.h5'])\n",
      "________Beam Blockage______\n",
      "60.13867\n",
      "________CLUTTER______\n",
      "44.20117\n",
      "________offset______\n",
      "43.20117\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-09/2017-09-09/raa01-rx_10000-1709090215-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "75\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 41.51082299999996)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 125.4860940000001)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20151216024501 0 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151216*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151216-S014742-E032015.010213.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20151216-S184555-E201828.010224.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2015/2015-12/2015-12-16/n_ppi_045deg/2015-12-16--02:45:01,00.mvol'])\n",
      "________Beam Blockage______\n",
      "48.804688\n",
      "________CLUTTER______\n",
      "48.804688\n",
      "________offset______\n",
      "46.804688\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2015/2015-12/2015-12-16/raa01-rx_10000-1512160245-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 41.05662300000006)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 107.883197)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170211153000 0 alt 3.5\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211-S144911-E162144.016802.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-02/2017-02-11/n_ppi_045deg/2017-02-11--15:30:00,00.mvol'])\n",
      "________Beam Blockage______\n",
      "37.583984\n",
      "________CLUTTER______\n",
      "37.583984\n",
      "________offset______\n",
      "34.083984\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-02/2017-02-11/raa01-rx_10000-1702111530-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 47.45134100000007)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 123.45165899999995)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20160601175950 1 alt 3.5\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20160601*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20160601-S002009-E015243.012826.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20160601-S171836-E185110.012837.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2016/2016-06/2016-06-01/n_ppi_045deg/2016-06-01--17:59:50,00.mvol'])\n",
      "________Beam Blockage______\n",
      "60.30664\n",
      "________CLUTTER______\n",
      "60.30664\n",
      "________offset______\n",
      "56.80664\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2016/2016-06/2016-06-01/raa01-rx_10000-1606011755-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 45.64480800000001)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 121.20911799999999)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170321043000 0 alt 3.5\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170321*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170321-S035103-E052338.017386.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-03/2017-03-21/n_ppi_045deg/2017-03-21--04:30:00,00.mvol'])\n",
      "________Beam Blockage______\n",
      "44.07617\n",
      "________CLUTTER______\n",
      "44.07617\n",
      "________offset______\n",
      "40.57617\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-03/2017-03-21/raa01-rx_10000-1703210430-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 47.48588700000005)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 122.54592200000002)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170211153000 0 alt 3.5\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211-S144911-E162144.016802.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-02/2017-02-11/n_ppi_045deg/2017-02-11--15:30:00,00.mvol'])\n",
      "________Beam Blockage______\n",
      "37.583984\n",
      "________CLUTTER______\n",
      "37.583984\n",
      "________offset______\n",
      "34.083984\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-02/2017-02-11/raa01-rx_10000-1702111530-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 47.1224400000001)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 122.75464699999998)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170810110333 0 neu 1.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170810*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170810-S102212-E115444.019599.V05A.HDF5'])\n",
      "New enigma\n",
      "/automount/radar/scans/2017/2017-08/2017-08-10/n_ppi_010deg/*20170810110333*.h5\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-08/2017-08-10/n_ppi_010deg/n_ppi_010deg_12345_20170810110333_00.h5'])\n",
      "________Beam Blockage______\n",
      "65.61719\n",
      "________CLUTTER______\n",
      "61.134766\n",
      "________offset______\n",
      "60.134766\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-08/2017-08-10/raa01-rx_10000-1708101105-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "75\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 36.29935699999987)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 117.40411700000004)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20171126032331 0 neu 1.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20171126*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20171126-S024618-E041852.021274.V05A.HDF5'])\n",
      "New enigma\n",
      "/automount/radar/scans/2017/2017-11/2017-11-26/n_ppi_010deg/*20171126032331*.h5\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-11/2017-11-26/n_ppi_010deg/n_ppi_010deg_12345_20171126032331_00.h5'])\n",
      "________Beam Blockage______\n",
      "49.679688\n",
      "________CLUTTER______\n",
      "48.185547\n",
      "________offset______\n",
      "47.185547\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-11/2017-11-26/raa01-rx_10000-1711260325-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "75\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 35.90111999999999)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 118.488067)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20140826221000 1 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140826*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140826-S133117-E150350.002799.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140826-S211407-E224640.002804.V05A.HDF5'])\n",
      "['/automount/radar-archiv/scans/2014/2014-08/2014-08-26/n_ppi_045deg/2014-08-26--22:10:00,00.mvol']\n",
      "________Beam Blockage______\n",
      "42.819336\n",
      "________CLUTTER______\n",
      "42.819336\n",
      "________offset______\n",
      "40.819336\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2014/2014-08/2014-08-26/raa01-rx_10000-1408262210-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 32.875676999999996)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 84.15136700000016)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20140806203538 1 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140806*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140806-S025632-E042905.002481.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20140806-S195445-E212718.002492.V05A.HDF5'])\n",
      "['/automount/radar-archiv/scans/2014/2014-08/2014-08-06/n_ppi_045deg/2014-08-06--20:35:38,00.mvol']\n",
      "________Beam Blockage______\n",
      "33.335938\n",
      "________CLUTTER______\n",
      "33.335938\n",
      "________offset______\n",
      "31.335938\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2014/2014-08/2014-08-06/raa01-rx_10000-1408062035-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 47.83331799999996)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 122.03982100000007)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20141008094000 1 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20141008*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20141008-S010511-E023741.003460.V05A.HDF5', '/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20141008-S084747-E102017.003465.V05A.HDF5'])\n",
      "['/automount/radar-archiv/scans/2014/2014-10/2014-10-08/n_ppi_045deg/2014-10-08--09:40:00,00.mvol']\n",
      "________Beam Blockage______\n",
      "40.648438\n",
      "________CLUTTER______\n",
      "40.648438\n",
      "________offset______\n",
      "38.648438\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2014/2014-10/2014-10-08/raa01-rx_10000-1410080940-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 33.73136799999975)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 86.63731300000018)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20141007023744 0 alt 2.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20141007*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20141007-S015721-E032951.003445.V05A.HDF5'])\n",
      "['/automount/radar-archiv/scans/2014/2014-10/2014-10-07/n_ppi_045deg/2014-10-07--02:37:44,00.mvol']\n",
      "________Beam Blockage______\n",
      "45.054688\n",
      "________CLUTTER______\n",
      "45.054688\n",
      "________offset______\n",
      "43.054688\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2014/2014-10/2014-10-07/raa01-rx_10000-1410070235-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 46.86607700000013)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 118.70936600000005)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170725152833 0 neu 1.0\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170725*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170725-S144757-E162031.019353.V05A.HDF5'])\n",
      "New enigma\n",
      "/automount/radar/scans/2017/2017-07/2017-07-25/n_ppi_010deg/*20170725152833*.h5\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-07/2017-07-25/n_ppi_010deg/n_ppi_010deg_12345_20170725152833_00.h5'])\n",
      "________Beam Blockage______\n",
      "58.64453\n",
      "________CLUTTER______\n",
      "58.64453\n",
      "________offset______\n",
      "57.64453\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-07/2017-07-25/raa01-rx_10000-1707251530-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "75\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 32.91543200000024)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 116.50259099999994)\n",
      "<type 'str'> <type 'int'> <type 'str'> <type 'float'>\n",
      "20170211152500 0 alt 3.5\n",
      "/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211*.HDF5\n",
      "('GPM:', ['/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.20170211-S144911-E162144.016802.V05A.HDF5'])\n",
      "[]\n",
      "Old enigma\n",
      "('Boxpol: ', ['/automount/radar/scans/2017/2017-02/2017-02-11/n_ppi_045deg/2017-02-11--15:25:00,00.mvol'])\n",
      "________Beam Blockage______\n",
      "36.432617\n",
      "________CLUTTER______\n",
      "36.432617\n",
      "________offset______\n",
      "32.932617\n",
      "________ATTCORR______\n",
      "________DPR Threshold______\n",
      "('RADOLANPFAD: ', '/automount/radar/dwd/rx/2017/2017-02/2017-02-11/raa01-rx_10000-1702111525-dwd---bin.gz')\n",
      "('Max --------------------------------------->', 150.0)\n",
      "('Min --------------------------------------->', 15.0)\n",
      "100\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/__main__.py:297: DeprecatedWarning: polar2lonlatalt_n is deprecated as of 0.11.3 and will be removed in 1.0.0. Use :func:`spherical_to_proj` instead.\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('Interpolationsdauer Boxpol auf DPR:', 46.4529950000001)\n",
      "('Interpolationsdauer RADOLAN auf DPR:', 119.74704800000018)\n",
      "('Gesamtlaufzeit: ', 44.957320349999996, 'min')\n"
     ]
    }
   ],
   "source": [
    "#!/usr/bin/env python\n",
    "\n",
    "import h5py\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "import glob\n",
    "import math\n",
    "import pandas as pd\n",
    "import wradlib\n",
    "from scipy import stats\n",
    "import matplotlib.cm as cm\n",
    "my_cmap = cm.get_cmap('jet',40)\n",
    "my_cmap.set_under('lightgrey')\n",
    "my_cmap.set_over('darkred')\n",
    "from pcc import get_miub_cmap as my_cmap\n",
    "from pcc import plot_radar\n",
    "from pcc import boxpol_pos\n",
    "from pcc import plot_borders\n",
    "import wradlib as wrl\n",
    "from osgeo import osr\n",
    "from satlib import ipoli_radi_stat\n",
    "from satlib import corcor\n",
    "Pos = boxpol_pos()\n",
    "blon0, blat0 = Pos['lon_ppi'], Pos['lat_ppi']\n",
    "bbx, bby = Pos['gkx_ppi'], Pos['gky_ppi']\n",
    "from time import *\n",
    "from pcc import cut_the_swath\n",
    "from satlib import good_overpasses_dpr_boxpol as overpasses_dpr_boxpol\n",
    "\n",
    "print (\"_\")\n",
    "\n",
    "\n",
    "\n",
    "tstart = clock()\n",
    "\n",
    "#print (overpasses_dpr_boxpol)\n",
    "\n",
    "for ii in overpasses_dpr_boxpol.keys():\n",
    "    ZP = overpasses_dpr_boxpol[ii][0],\n",
    "    pfadnr = overpasses_dpr_boxpol[ii][1]\n",
    "    enigma = overpasses_dpr_boxpol[ii][2]\n",
    "    offset = overpasses_dpr_boxpol[ii][3]\n",
    "\n",
    "    ZP = ZP[0]\n",
    "\n",
    "    print type(ZP), type(pfadnr), type(enigma), type(offset)\n",
    "    print ZP, pfadnr, enigma, offset\n",
    "    #ZP = '20180125170330'; pfadnr=0; enigma='neu'; offset=1\n",
    "\n",
    "\n",
    "    # Pfade zu den Dateien festlegen\n",
    "    #-----------------------------#\n",
    "\n",
    "    year = ZP[0:4]; ye = ZP[2:4]; m = ZP[4:6]; d = ZP[6:8]; ht = ZP[8:10]; mt = ZP[10:12]; st = ZP[12:14]\n",
    "\n",
    "    print('/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.' + year + m + d + '*.HDF5')\n",
    "    #pfad_radar = glob.glob('/automount/ags/velibor/gpmdata/dpr/2A.GPM.DPR.V6-20160118.' + year + m + d + '*.HDF5')\n",
    "    pfad_radar = glob.glob('/automount/ags/velibor/gpmdata/dprV7/2A.GPM.DPR.V7-20170308.' + year + m + d + '*.HDF5')\n",
    "\n",
    "    print ('GPM:',pfad_radar)\n",
    "    pfad_radar = pfad_radar[pfadnr]\n",
    "    #pfad_radar_Ku = pfad_radar[0]\n",
    "\n",
    "    deg_scan =  [\"/ppi_1p5deg/\",\"/ppi_2p4deg/\",\"/ppi_3p4deg/\",\n",
    "                 \"/n_ppi_010deg/\",\"/n_ppi_045deg/\",\n",
    "                 \"/n_ppi_082deg/\",\"/n_ppi_110deg/\",\"/n_ppi_140deg/\",\n",
    "                 \"/n_ppi_180deg/\",\"/n_ppi_280deg/\",\"/n_vertical_scan/\"][4]\n",
    "\n",
    "\n",
    "    if enigma=='neu':\n",
    "        print ('New enigma')\n",
    "        deg_scan =  [\"/n_ppi_010deg/\"][0]\n",
    "\n",
    "        boxpolpath = '/automount/radar/scans/' + year+ \"/\" +year +\"-\"+ m + \"/\" + year+ \"-\" + m +\"-\"+ d +\\\n",
    "                                   deg_scan+\"*\"+year+m+d+ht+mt+st+\"*.h5\"\n",
    "        print (boxpolpath)\n",
    "        ppi_datapath=glob.glob(boxpolpath)\n",
    "\n",
    "        print ('Boxpol: ',ppi_datapath)\n",
    "        ppi_datapath = ppi_datapath[0]\n",
    "\n",
    "\n",
    "    else:\n",
    "        try:\n",
    "            ppi_datapath=glob.glob('/automount/radar-archiv/scans/' + year+ \"/\" +\n",
    "                                   year +\"-\"+ m + \"/\" + year+ \"-\" + m +\"-\"+ d +\n",
    "                                   deg_scan+ year + \"-\" + m +\"-\"+ d + \"--\" +ht +\n",
    "                                   \":\"+mt+\":\"+st+\",*.mvol\")\n",
    "            print ppi_datapath\n",
    "            ppi_datapath = ppi_datapath[0]\n",
    "\n",
    "        except:\n",
    "            ppi_datapath=glob.glob('/automount/radar/scans/' + year+ \"/\" +\n",
    "                                   year +\"-\"+ m + \"/\" + year+ \"-\" + m +\"-\"+\n",
    "                                   d + deg_scan+ year + \"-\" + m +\"-\"+ d +\n",
    "                                   \"--\" +ht +\":\"+mt+\":\"+st+\",*.mvol\")\n",
    "            print ('Old enigma')\n",
    "            print ('Boxpol: ',ppi_datapath)\n",
    "            ppi_datapath = ppi_datapath[0]\n",
    "\n",
    "\n",
    "    # Wichtige Parameter festlegen\n",
    "    TH = 15 #Threshold um Nullen fuer Niederschlag raus zu filtern\n",
    "\n",
    "\n",
    "    #################################################### PPI BoxPol Daten einlesen\n",
    "    #------------------------------------------------------------------------------\n",
    "\n",
    "    ppi=h5py.File(ppi_datapath,'r')\n",
    "    data, attrs = wradlib.io.read_GAMIC_hdf5(ppi_datapath)\n",
    "\n",
    "    ZH0 = data['SCAN0']['ZH']['data']\n",
    "    #ZV = data['SCAN0']['ZV']['data']\n",
    "    #ZH0 = (ZH + ZV)*0.5\n",
    "    PHIDP = data['SCAN0']['PHIDP']['data']\n",
    "    r = attrs['SCAN0']['r']\n",
    "    az = attrs['SCAN0']['az']\n",
    "    elevation=attrs['SCAN0']['elevation']\n",
    "    lon_ppi = attrs['VOL']['Longitude']\n",
    "    lat_ppi = attrs['VOL']['Latitude']\n",
    "    alt_ppi = attrs['VOL']['Height']\n",
    "    rho = data['SCAN0']['RHOHV']['data']\n",
    "\n",
    "    R = ZH0\n",
    "\n",
    "    print (\"________Beam Blockage______\")\n",
    "    R[151:165]=np.nan\n",
    "    print (np.nanmax(R))\n",
    "\n",
    "    print (\"________CLUTTER______\")\n",
    "    rho_th  = 0.85\n",
    "    R[rho<= rho_th] = np.nan\n",
    "    print (np.nanmax(R))\n",
    "\n",
    "    print (\"________offset______\")\n",
    "    R = R - offset\n",
    "    print (np.nanmax(R))\n",
    "    #?\n",
    "    print (\"________ATTCORR______\")\n",
    "\n",
    "    \"\"\"pia = wrl.atten.correctAttenuationHB(\n",
    "        R,\n",
    "        coefficients = dict(a=4.57e-5, b=0.731, gate_length=1.0),\n",
    "        mode=\"warn\",\n",
    "        thrs=59.)\n",
    "    pia[pia > 4.8] = 4.8\n",
    "\n",
    "    print (\"________ATTCORR2______\")\n",
    "    R = R + pia\n",
    "    print (np.nanmax([R,pia]))\"\"\"\n",
    "\n",
    "    print (\"________DPR Threshold______\")\n",
    "    Z_boxpol = R\n",
    "\n",
    "    ### Threshold for DPR sensitivity\n",
    "    Z_boxpol[Z_boxpol<TH]=np.nan\n",
    "\n",
    "\n",
    "    ################################################################# DPR Einlesen\n",
    "    # -----------------------------------------------------------------------------\n",
    "    scan = 'MS'\n",
    "    gpmku = h5py.File(pfad_radar, 'r')\n",
    "    gpmku_HS = gpmku[scan]['SLV']\n",
    "    dpr_lat_1 = np.array(gpmku[scan]['Latitude'])\n",
    "    dpr_lon_1 = np.array(gpmku[scan]['Longitude'])\n",
    "    Z_dpr = np.array(gpmku_HS['zFactorCorrectedNearSurface'])\n",
    "    Z_dpr[Z_dpr < TH] = np.nan\n",
    "\n",
    "    ## Einlesen von Phase un Raintype\n",
    "    dpr_raintype = np.array(gpmku[scan]['CSF']['typePrecip'], dtype=float)\n",
    "    dpr_phase = np.array(gpmku_HS['phaseNearSurface'], dtype=float)\n",
    "\n",
    "\n",
    "    ############################################################## RADOLAN einlesen\n",
    "    # -----------------------------------------------------------------------------\n",
    "    mtt = mt\n",
    "    mtt = str(int(round(float(mtt)/5.0)*5.0))\n",
    "\n",
    "    if mtt == '0':\n",
    "        mtt = '00'\n",
    "    if mtt == '5':\n",
    "        mtt = '05'\n",
    "    if mtt == '60':\n",
    "        mtt = '55'\n",
    "\n",
    "    r_pro = 'rx'\n",
    "\n",
    "    pfad = ('/automount/radar/dwd/'+ r_pro +'/'+str(year)+'/'+str(year)+'-'+\n",
    "            str(m)+'/'+ str(year)+'-'+str(m)+'-'+str(d)+'/raa01-'+r_pro+'_10000-'+\n",
    "            str(ye)+str(m)+ str(d)+str(ht)+str(mtt)+'-dwd---bin.gz')\n",
    "\n",
    "    print ('RADOLANPFAD: ', pfad)\n",
    "\n",
    "    pfad_radolan = pfad[:-3]\n",
    "\n",
    "    try:\n",
    "        rw_filename = wradlib.util.get_wradlib_data_file(pfad)\n",
    "    except EnvironmentError:\n",
    "        rw_filename = wradlib.util.get_wradlib_data_file(pfad_radolan)\n",
    "\n",
    "    rwdata, rwattrs = wradlib.io.read_RADOLAN_composite(rw_filename)\n",
    "\n",
    "    radolan_zeit = rwattrs['datetime'].strftime(\"%Y.%m.%d -- %H:%M:%S\")\n",
    "\n",
    "    radolan_grid_xy = wradlib.georef.get_radolan_grid(900,900)\n",
    "    x = radolan_grid_xy[:,:,0]\n",
    "    y = radolan_grid_xy[:,:,1]\n",
    "    rwdata = np.ma.masked_equal(rwdata, -9999) / 2 - 32.5\n",
    "\n",
    "    ### Threshold for DPR sensitivity\n",
    "    rwdata[rwdata<TH]=np.nan\n",
    "    Z_radolan = rwdata\n",
    "\n",
    "\n",
    "    ######################################################## Cut the Swath for Bonn\n",
    "    # -----------------------------------------------------------------------------\n",
    "\n",
    "    dpr_lon, dpr_lat, Z_dpr = cut_the_swath(dpr_lon_1,dpr_lat_1,Z_dpr, eu=0)\n",
    "    dpr_lon, dpr_lat, dpr_raintype = cut_the_swath(dpr_lon_1,dpr_lat_1,dpr_raintype, eu=0)\n",
    "    dpr_lon, dpr_lat, dpr_phase = cut_the_swath(dpr_lon_1,dpr_lat_1,dpr_phase, eu=0)\n",
    "\n",
    "    ######################################################## Koordinaten Projektion\n",
    "    # -----------------------------------------------------------------------------\n",
    "\n",
    "    proj_stereo = wrl.georef.create_osr(\"dwd-radolan\")\n",
    "    proj_wgs = osr.SpatialReference()\n",
    "    proj_wgs.ImportFromEPSG(4326)\n",
    "\n",
    "    dpr_lon, dpr_lat = wradlib.georef.reproject(dpr_lon, dpr_lat,\n",
    "                                                projection_target=proj_stereo ,\n",
    "                                                projection_source=proj_wgs)\n",
    "\n",
    "    blon, blat = wradlib.georef.reproject(blon0, blat0,\n",
    "                                          projection_target=proj_stereo ,\n",
    "                                          projection_source=proj_wgs)\n",
    "\n",
    "\n",
    "    ############################################################### Dpr zuschneiden\n",
    "    #------------------------------------------------------------------------------\n",
    "    print ('Max --------------------------------------->', r[-1]/1000.)\n",
    "    inner_r = 15. # in km!!!!\n",
    "\n",
    "    print ('Min --------------------------------------->', inner_r)\n",
    "\n",
    "    lon0, lat0, radius = blon, blat, r[-1]/1000.\n",
    "    rr = np.sqrt((dpr_lat - lat0)**2 + (dpr_lon - lon0)**2)\n",
    "    position = rr < radius\n",
    "\n",
    "    # Maximum Range\n",
    "    Z_dpr[np.where(rr > radius)] = np.nan\n",
    "    dpr_raintype[np.where(rr > radius)] = np.nan\n",
    "    dpr_phase[np.where(rr > radius)] = np.nan\n",
    "\n",
    "    # Minimum Range\n",
    "    Z_dpr[np.where(rr < inner_r)] = np.nan\n",
    "    dpr_raintype[np.where(rr < inner_r)] = np.nan\n",
    "    dpr_phase[np.where(rr < inner_r)] = np.nan\n",
    "\n",
    "\n",
    "    ########################################################### RADOLAN zuschneiden\n",
    "    #------------------------------------------------------------------------------\n",
    "\n",
    "    rr2 = np.sqrt((y - lat0)**2 + (x - lon0)**2)\n",
    "    position2 = rr2 < radius\n",
    "\n",
    "    # Maximum Range\n",
    "    Z_radolan[np.where(rr2 > radius)] = np.nan\n",
    "\n",
    "    # Minimum Range\n",
    "    Z_radolan[np.where(rr2 < inner_r)] = np.nan\n",
    "\n",
    "\n",
    "    ########################################################### BoXPol zuschneiden\n",
    "    #------------------------------------------------------------------------------\n",
    "    # Minimum Range\n",
    "    a = (inner_r*1000.)/(r[1]-r[0])# inner_r in m!!!!\n",
    "    a = a.astype(int) # bin bis zu welchem radius alles NaN gesetzt wird\n",
    "    Z_boxpol[:,0:a]=np.nan\n",
    "    print (a)\n",
    "\n",
    "\n",
    "    ############################################################### Inverse Dezibel\n",
    "    #------------------------------------------------------------------------------\n",
    "    Z_boxpol = wradlib.trafo.idecibel(Z_boxpol)\n",
    "    Z_dpr = wradlib.trafo.idecibel(Z_dpr)\n",
    "    Z_radolan = wradlib.trafo.idecibel(Z_radolan)\n",
    "\n",
    "    ################################################## BoXPol interpolieren auf DPR\n",
    "    #------------------------------------------------------------------------------\n",
    "    t1 = clock()\n",
    "    radar_location = (lon_ppi, lat_ppi, alt_ppi)\n",
    "    #elevation = 1.5\n",
    "    azimuths = az\n",
    "    ranges = r\n",
    "    polargrid = np.meshgrid(ranges, azimuths)\n",
    "    lon, lat, alt = wradlib.georef.polar2lonlatalt_n(polargrid[0], polargrid[1],\n",
    "                                                     elevation, radar_location)\n",
    "    lon, lat = wradlib.georef.reproject(lon, lat, projection_target=proj_stereo ,\n",
    "                                        projection_source=proj_wgs)\n",
    "\n",
    "    grid_xy = np.vstack((dpr_lon.ravel(), dpr_lat.ravel())).transpose()\n",
    "\n",
    "    xy=np.concatenate([lon.ravel()[:,None],lat.ravel()[:,None]], axis=1)\n",
    "\n",
    "\n",
    "\n",
    "    Z_boxpol_ipoli, Z_boxpol_ipoli_std, Z_boxpol_ipoli_median, Z_boxpol_ipoli_max, Z_boxpol_ipoli_min = ipoli_radi_stat(xy,Z_boxpol.ravel(),grid_xy,2.5)\n",
    "    \n",
    "    Z_boxpol_ipoli = Z_boxpol_ipoli.reshape(dpr_lon.shape)\n",
    "    Z_boxpol_ipoli_std = Z_boxpol_ipoli_std.reshape(dpr_lon.shape)\n",
    "    Z_boxpol_ipoli_median = Z_boxpol_ipoli_median.reshape(dpr_lon.shape)\n",
    "    Z_boxpol_ipoli_max = Z_boxpol_ipoli_max.reshape(dpr_lon.shape)\n",
    "    Z_boxpol_ipoli_min = Z_boxpol_ipoli_min.reshape(dpr_lon.shape)\n",
    "\n",
    "    t2 = clock()\n",
    "    print ('Interpolationsdauer Boxpol auf DPR:', t2 - t1)\n",
    "\n",
    "\n",
    "    ################################################# RADOLAN interpolieren auf DPR\n",
    "    #------------------------------------------------------------------------------\n",
    "\n",
    "    xy_rado = np.vstack((x.ravel(), y.ravel())).transpose()\n",
    "\n",
    "    Z_radolan_ipoli, Z_radolan_ipoli_std, Z_radolan_ipoli_median, Z_radolan_ipoli_max, Z_radolan_ipoli_min = ipoli_radi_stat(xy_rado,Z_radolan.ravel().filled(np.nan),grid_xy,2.5)\n",
    "    \n",
    "    Z_radolan_ipoli = Z_radolan_ipoli.reshape(dpr_lon.shape)\n",
    "    Z_radolan_ipoli_std = Z_radolan_ipoli_std.reshape(dpr_lon.shape)\n",
    "    Z_radolan_ipoli_median = Z_radolan_ipoli_median.reshape(dpr_lon.shape)\n",
    "    Z_radolan_ipoli_max = Z_radolan_ipoli_max.reshape(dpr_lon.shape)\n",
    "    Z_radolan_ipoli_min = Z_radolan_ipoli_min.reshape(dpr_lon.shape)\n",
    "\n",
    "    t3 = clock()\n",
    "    print ('Interpolationsdauer RADOLAN auf DPR:', t3 - t2)\n",
    "\n",
    "    ################################################################### in  Dezibel\n",
    "    #------------------------------------------------------------------------------\n",
    "    Z_boxpol = wradlib.trafo.decibel(Z_boxpol)\n",
    "    Z_dpr = wradlib.trafo.decibel(Z_dpr)\n",
    "    Z_radolan = wradlib.trafo.decibel(Z_radolan)\n",
    "    Z_boxpol_ipoli = wradlib.trafo.decibel(Z_boxpol_ipoli)\n",
    "    Z_radolan_ipoli = wradlib.trafo.decibel(Z_radolan_ipoli)\n",
    "    \n",
    "    Z_boxpol_ipoli_std = wradlib.trafo.decibel(Z_boxpol_ipoli_std)\n",
    "    Z_radolan_ipoli_std = wradlib.trafo.decibel(Z_radolan_ipoli_std)\n",
    "    Z_boxpol_ipoli_median = wradlib.trafo.decibel(Z_boxpol_ipoli_median)\n",
    "    Z_radolan_ipoli_median = wradlib.trafo.decibel(Z_radolan_ipoli_median)\n",
    "    Z_boxpol_ipoli_max = wradlib.trafo.decibel(Z_boxpol_ipoli_max)\n",
    "    Z_radolan_ipoli_max = wradlib.trafo.decibel(Z_radolan_ipoli_max)\n",
    "    Z_boxpol_ipoli_min = wradlib.trafo.decibel(Z_boxpol_ipoli_min)\n",
    "    Z_radolan_ipoli_min = wradlib.trafo.decibel(Z_radolan_ipoli_min)\n",
    "\n",
    "\n",
    "    #zwischenspeichern\n",
    "    np.save('/automount/ags/velibor/data/dpr_boxpol_radolan/stat/minus_bias/'+ZP+'.npy',\n",
    "            [Z_dpr, Z_boxpol, Z_boxpol_ipoli, Z_radolan, Z_radolan_ipoli, dpr_phase, dpr_raintype, dpr_lon, dpr_lat, \n",
    "             Z_boxpol_ipoli_std, Z_boxpol_ipoli_median, Z_boxpol_ipoli_max, Z_boxpol_ipoli_min,\n",
    "             Z_radolan_ipoli_std, Z_radolan_ipoli_median, Z_radolan_ipoli_max, Z_radolan_ipoli_min])\n",
    "\n",
    "tend = clock()\n",
    "print ('Gesamtlaufzeit: ',(tend-tstart)/60., 'min')\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:66: DeprecationWarning: metadata {'config': True} was set from the constructor.  Metadata should be set using the .tag() method, e.g., Int().tag(key1='value1', key2='value2')\n",
      "  inline backend.\"\"\"\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:71: DeprecationWarning: metadata {'config': True} was set from the constructor.  Metadata should be set using the .tag() method, e.g., Int().tag(key1='value1', key2='value2')\n",
      "  'retina', 'jpeg', 'svg', 'pdf'.\"\"\")\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:85: DeprecationWarning: metadata {'config': True} was set from the constructor.  Metadata should be set using the .tag() method, e.g., Int().tag(key1='value1', key2='value2')\n",
      "  use `figure_formats` instead)\"\"\")\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:95: DeprecationWarning: metadata {'config': True} was set from the constructor.  Metadata should be set using the .tag() method, e.g., Int().tag(key1='value1', key2='value2')\n",
      "  \"\"\"\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:114: DeprecationWarning: metadata {'config': True} was set from the constructor.  Metadata should be set using the .tag() method, e.g., Int().tag(key1='value1', key2='value2')\n",
      "  \"\"\")\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/ipykernel/pylab/config.py:44: DeprecationWarning: InlineBackend._config_changed is deprecated: use @observe and @unobserve instead.\n",
      "  def _config_changed(self, name, old, new):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/traitlets/traitlets.py:770: DeprecationWarning: A parent of InlineBackend._config_changed has adopted the new @observe(change) API\n",
      "  clsname, change_or_name), DeprecationWarning)\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:98: DeprecationWarning: DisplayFormatter._formatters_default is deprecated: use @default decorator instead.\n",
      "  def _formatters_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:677: DeprecationWarning: PlainTextFormatter._deferred_printers_default is deprecated: use @default decorator instead.\n",
      "  def _deferred_printers_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:669: DeprecationWarning: PlainTextFormatter._singleton_printers_default is deprecated: use @default decorator instead.\n",
      "  def _singleton_printers_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:672: DeprecationWarning: PlainTextFormatter._type_printers_default is deprecated: use @default decorator instead.\n",
      "  def _type_printers_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:672: DeprecationWarning: PlainTextFormatter._type_printers_default is deprecated: use @default decorator instead.\n",
      "  def _type_printers_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:677: DeprecationWarning: PlainTextFormatter._deferred_printers_default is deprecated: use @default decorator instead.\n",
      "  def _deferred_printers_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:92: DeprecationWarning: DisplayFormatter._ipython_display_formatter_default is deprecated: use @default decorator instead.\n",
      "  def _ipython_display_formatter_default(self):\n",
      "/home/velibor/anaconda2/envs/wradlib/lib/python2.7/site-packages/IPython/core/formatters.py:669: DeprecationWarning: PlainTextFormatter._singleton_printers_default is deprecated: use @default decorator instead.\n",
      "  def _singleton_printers_default(self):\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEACAYAAABWLgY0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+UVGed5/H3t7rTzY+GhIEJ2E1CTGfWBLWHeM7mx2JC\ngYDs4ELG9WQz0rNmRHHOTlCio8E4vQ12nDUeJsHJrH9Eo8RN4o/V3STIMgIhdTqJSUxWkgYlatpL\n23QLmI1CuvlNP/vHvdVUN1XQVXXr163P65w+VN/6cZ+b2/n208/zfb6POecQEZFoipW6ASIiUjgK\n8iIiEaYgLyISYQryIiIRpiAvIhJhCvIiIhE25iBvZvVm9qKZ7TKz3WbWHhxvN7P9Zvaz4GtJ4Zor\nIiLZsGzy5M1sgnPuqJnVAM8BnwT+PfCWc+6+ArVRRERylNVwjXPuaPCwHqgFkr8hLMxGiYhIOLIK\n8mYWM7NdwAFgu3PupeCpO8zsFTP7hpldHHorRUQkJ1kN1wy/yWwy8L+B1cDvgTecc87M7gHe5pxb\nGW4zRUQkFzkFeQAzawMGU8fizWwWsNk515Lm9SqSIyKSA+dczkPi2WTXTEsOxZjZeGAR8JqZzUh5\n2QeBPedpaGS/2tvbS94GXZ+uTdcXva981Wbx2rcBD5tZDP+Xw/ecc//HzL5tZnOAIWAf8Im8WyUi\nIqEYc5B3zu0G3pPm+H8OtUUiIgXgeR5tbW309fXR1NRER0cHb3/720vdrILLpicv5xGPx0vdhIKK\n8vVF+dpA1wd+gF+0aBHd3d3Dx1544QW2b98e+UCf88Rr1icyc8U6l4hIqtbWVh599NFzjq9YsYJH\nHnmkBC0aOzPDFWPiVUSkUvX19aU93t/fX+SWFJ+CvIhEXlNTU9rjjY2NRW5J8Wm4RkQiL92YfHNz\nc0WMyec7XKMgLyJVIZld09/fT2NjY8Vk1yjIi4hEmCZeRUQkIwV5EZEIU5AXEYkwBXkRkQhTkBcR\niTAFeRGRCFOQFxGJMAV5EZEIU5AXEYkwBXkRkQhTkBcRiTAFeRGRCBtzkDezejN70cx2mdluM2sP\njk8xs21m9ksz+7GZXVy45oqISDayqkJpZhOcc0fNrAZ4Dvgk8B+B/+ec+4qZ3QVMcc6tTfNeVaEU\nEclSUatQOueOBg/r8TcBd8By4OHg+MPALbk2RkREwpVVkDezmJntAg4A251zLwHTnXMHAZxzB4BL\nw2+miIjkItue/JBz7lpgJnCdmb0Tvzc/4mVhNU5ERPJTm8ubnHNHzCwBLAEOmtl059xBM5sBHMr0\nvnXr1g0/jsfjxOPxXE4vIhJZiUSCRCIR2ueNeeLVzKYBp5xzh81sPPBj4MvAPOBN59y9mngVkXKU\n3N+1r6+PpqamitnfFYq4x6uZvRt/YjUWfH3POfclM/sT4PvAZUAPcKtz7o9p3q8gLyJF19nZydKl\nSxkYGBg+1tzczPbt2ysi0GsjbxGRDDzPo6WlZUSAT1qxYgWPPPJICVqVHW3kLSKSwZo1a9IGeID+\n/v4it6Y0FORFJJI8z2Pbtm0Zn29sbCxia0pHQV5EIqmtrY3jx4+nfa6hoYGOjo4it6g0FORFJJL6\n+vrSHq+pqWHLli0VMekaBgV5EYmkpqamtMc/8IEPcPPNNxe5NaWj7BoRiSTP81i0aBHd3d3Dxyop\ndTJJKZQiIhkkF0H19/fT2NhYUYugkhTkRUQiTHnyIiKSkYK8iEiEKciLSCR5nkdrayvz58+ntbUV\nz/NK3aSS0Ji8iEROpRclS6UxeRGRFJ2dnbzvfe87p2ZNd3c3bW1tJWpV6SjIi0jF8zyPW265halT\npzJv3jxOnz6d9nXVUpQsVU47Q4mIlIvOzk6WLFnCsWPHLvjaailKlkpj8iJSsZJDM5l67qkaGhro\n6urSmLyISCXwPI+lS5eOKcDX1tZWVVGyVAryIlKRzrchSKqGhgaeeuqpqipKlkpBXkQqTmdnJz/6\n0Y/O+5r6+nqWLVtGV1dX1QZ4yCLIm9lMM9tpZj83s91mtjo43m5m+83sZ8HXksI1V0SqXXIcfmho\nKONrFi5cyN69e3niiSeqcogm1ZgnXs1sBjDDOfeKmTUA/xdYDvwn4C3n3H0XeL8mXkUkJ8lqknv2\n7KGrq4tMsaSmpoadO3dGquee78TrmFMonXMHgAPB4wEz2wskq/Ln3AARkXQ8z+POO+/k6aef5siR\nI2N6T7VtCDIWOaVQmtkVQAJ4F/AZ4HbgMPAy8Bnn3OE071FPXkQuyPM8Vq1axY4dO7J6X6WmSF5I\n0XryKSdsAH4AfCro0X8N+KJzzpnZPcB9wMp07123bt3w43g8Tjwez6XNIhJR2SxsShWlFMlEIkEi\nkQjt87LqyZtZLfAjYKtz7qtpnp8FbHbOtaR5Tj15ERkhOdb++uuv43kehw4dyvozamtrI50iWeye\n/DeBX6QGeDObEYzXA3wQ2JNrY0SkOiSHZHbu3HneLJkLaWhoYMuWLZEN8GHIJrtmLtAJ7AZc8HU3\n8GFgDjAE7AM+4Zw7mOb96smLVKnkJOozzzzD4cOHOXPmTM6fVVNTQ0tLC7Nnz67IPVuzpT1eRaQs\npAbyt956CzNj/Pjx1NTU8Oabb+b9+VOmTOGmm25i48aNkQ/sqRTkq0iP57GprY2hvj5iTU3c3tHB\nrCr6YZfyk0uaY7YWLlzIgw8+WFWBPZWCfJXo8TweWLSI9d3dTAQGgfbmZlZv365AL0WR2lMfHBzk\noosuYnBwMOPCpHxNmDCBrVu3Vv14u4J8lVjf2srfP/ooE1OODQIbVqyg/ZFHStUsqQK55q1ny8y4\n5JJLqKur4/rrr6+6YZlMip4nL6Ux1Nc3IsADTASGqnCnGymsMCdJxyIWi7FgwYKqHpIpJAX5ChFr\namIQzunJx6pwpxspjGRw37p1KydPnizouerq6pg0aRJz585Vj73AFOQrQI/n8ceBAVaPG8cDx4+P\nHJPv6Ch186RCFbrHXldXx/jx46mvr+fMmTPU1tZqGKYENCZf5p7r7OShpUt5YGCAN4BvAL8ZN47G\nxYu5Y+NGTbpK1go5xl6taY6FpDH5COvxPO5dupTvDAwwEX+opgMYPH6cDZMmKcDLBRUyI2by5MlM\nmjSJyy67jObm5qpYmFSJFOTL2Ka2NlqCAJ9KE66SlFr7Zf/+/QwMDHD06FGccwwNDZ1TMuDEiRN5\nnU+TpJVHQb6MDfX1cRFowlWGFTvzpba2llgspknSCqYgX8ZiTU3cCrQD62F4wnV1QwPtmnCtCsUO\n6klaiBQdmngtY8lVriu7u/k+cAroamjgri1bmKv/+SKtmOmMSUprLE9a8Rpxw/Vq+vuJNTaqXk2E\nFaMOzGgaYy9/yq6JMBUkqx657oiULWXEVB/15MtUan68CpJFU9g9dzNj0qRJHD9+fDirRpOmlU/D\nNRHT43l0rFpFz44dPM65WTUqSFYZRtdWT/ezf/r06bzPo3H06NNwTYT0eB5fmjePQ729XA/Kj68w\nyZz1PXv2sHv37ry2tUtH6YySCwX5MtHjedxx883M2b+fS0H58WVudE/9zJkzoQf1JKUzSj40XFNi\nPZ7HxjvvpH/rVmpOnuTPgBhwK/AQafLju7o0Jl9ko1eVHj58mIGBgYKfV3VgBIo4XGNmM4FvA9Px\nN+3+unPun81sCvA9YBb+Rt63OucO59qgapKcXG0aGODfBMeGgNvwA/xKYAN+fvzLtbW0bdmiAF9E\nyUJeO3fuLFgvPR313CVMsSxeexr4tHPuncCNwN+Z2dXAWmCHc+4dwE7g8+E3M3qSxcceGBgghn8j\nPgYcBL6OH+CTC6BenTCBtqee0gKoIvA8j1tuuYWLL76YK6+8kh07dhQlwNfV1TF16lSWLVvGnj17\nFOAlNGPuyTvnDgAHgscDZrYXmAksB+YFL3sYSOAHfjmP1OJjMfzfoNOALwAbg38nAn+cOZN/6exU\nD74AkuPqzz///HCmy5tvvlmQc9XWjvxfTROoUiw5Tbya2RXAHOAFYLpz7iD4vwjM7NLQWhdhqcXH\nbge+BLThlxK+Pzh+9+WXc08ioQAfkmLWgampqaGlpYXZs2drwZGUVNZB3swagB8Anwp69KNnUzPO\nrq5bt274cTweJx6PZ3v6yBhdfOwLwD/h/1k0fepUGufO5dPaFCRro+unjxs3jvr6eo4ePVrQyVKl\nN0pYEokEiUQitM/LKrvGzGqBHwFbnXNfDY7tBeLOuYNmNgN42jl3TZr3KrsmRXLS9bMDAyo+FoJC\n7naUKrmq9OTJk0ycOFEBXQqu2Iuhvgn8IhngA0/ijzjcC3wEeCLXxlQLBfhwFLOglwp5SaUac0/e\nzOYCncBu/CEZB9wN/BQ/EeQyoAc/hfKPad6vnjx+Vs3qlpbhLf2SVLLgXBcqDRBGWYALUa66lFrR\nevLOueeAmgxPL8y1AdVGW/pllrroyPM8Dh06VNTzT5s2DecctbW1XH/99QrsEgkqa1Bk2tLPV8yy\nAJmouJdUAwX5Iqu2Lf1G56LX1tYWPNMl1eTJkxk3bhwnTpzgxIkTmiyVqqMgX0Q9nkfvwYPcE4vx\nD0NDwyULkpOuUUqXLFa2SzqaJBU5SwXKiiRZRriht5ePM7JkwdqtWyOVVVOsXY5SaehFokqbhlSA\n1DLCa4n2RiCe59HS0lKQ4ZjU0gBaeCTVQpuGlLnnOjv55yVLmHTsGDGiuxFIcux969atnDx5Mu/P\nMzMuueQS6urqlOkikgcF+ZAl68O/9swzDL71FuNOneL64LkhopNVMzo75tSpU3l9nsoCiBSGhmtC\nlOy1Tz52jEn4wdzhLy74KH4Rsgb8ImTJrJq7L7+cT5dxEbJ0C5KGhoZyTnecNm0ap06dUqaLyBhp\nTL4MpG6+fUNwbC3wFc6WEV4LvIFfRriH8i4jnBrYwyq9u3DhQmW7iORAY/IlNnrz7eQuLMk68bfi\nlw4u9zLChawDs3z5ch5//PFQP1NExkY9+Tx9dvlyxj35JA5/8+1kNZVkz/0B/F2eHsTfG5FYjKYF\nC/jsgw+WRYAvdD775ZdfTiKRUA9eJEcarimhHs/j7tmzuer48RG9dgdMwu+5vwF8A3i9vp6Z738/\nd5RRjXjP85g3bx69vb2hfq5y1kXCoyBfIqm57wAf5uzm2w8Ce4ETwIwpU5h5001lEdxHlxgYHBzk\nxIkToXy2qjWKFIaCfAkkx+EHenv5b5zNmimnlayFKACmxUgixacgXwLrW1s5/eijQPlkzaT20o8f\nP562/nouVAdGpLSUXVMCQ319xPBz35PVJJNZM6vHjy9KgC/0ptQafhGJBgX5LPV4Hnv27WM2MA1Y\nDWzAX806BFy8eHHoAT51M439+/dz+PDhgpXqnTlzJp2dnQrsIhGh4Zos9HgeDyxaxMrubu6nsKtX\ni7l/aZLSHUXKj8bki+izy5ez7sknmYg//r4RP/fd6utpDjE9slCpjaNpuzuR8le0IG9mDwEfAA46\n51qCY+34SSXJzTjvds79a4b3V3SQf66zk6/On8/302SotM+fz/qdO3P63NFDMQMDAxw5ciSUSdNU\nKgAmUpmKOfH6LfwFnN8edfw+59x9uTag3I2uSxNmFcnOzk6WLl1akPH1GTNmcMUVV9Dc3ExHR4cC\nukiVGnOQd849a2az0jyV82+YcpesKnns2DGuBz5Gmr1Zx4/PaW9Wz/NCDfBaZSoi6YSRXXOHmf01\n8DLwGefc4RA+s+R6PI97ly7lz48dG65LE2Y2zZo1a3IK8GbGpEmTOHnypEr1isgF5RvkvwZ80Tnn\nzOwe4D78lf1prVu3bvhxPB4nHo/nefrC2dTWRsvAADHOVpNM9uLbCXrxDQ2033//mD4v37z2uro6\nlixZooAuEnGJRIJEIhHa52WVXRMM12xOTryO9bng+YqaeL3rxhuZ8MILnGZkXZpk2YKXa2tpe+qp\nMZUtyDdbRqmNItUr34nX2IVfMvJ8pIzBm9mMlOc+COzJtSHloMfzWN/ayso5c/jZiy9yK3AQ+Doj\nA/yrEyZkFeBvvvnmMQf4yZMnU1dXR11dHVOnTmXZsmUK8CKSs2xSKB8D4sBU/NjXDswH5uAPT+8D\nPuGcO5jh/WXZk0/uyfrK00/ztiNHaAPuAu5lZFXJffgXednChRlrwafbKu/06dPnvC4Tba4hIqNp\nMVQektUkXW8vB4Hv4E+qngLuwV/wtAk/uMeAP9xwAxuff37EZ4S1MrWhoYGuri712EVkBBUoy9Ho\nevDT8dMih/AzaQaBWfh/rhB8v6G5ecTiJc/zOHToUJpPz86ECRPYsmWLAryIhK4qe/Kp9eD/LDgW\nA/4evyd/K/5QTWo+/N8AP8T/JRCGWCzGlClTlAIpIuel4ZocJPdlTZWaQTM6k+YnwNMhnl/ZMiIy\nVgryWUrdl/Wj+Ls6JfdkTe7s9AbwSvD6A8DrIZxXPXcRyYWCfJZSe/Gpuzq9BpwBavCDehiBXaUG\nRCRfCvJZ6PE81l59NV8+eXJ4X9bUevArgCfyPId2VBKRMCnIj0GP57GprY3nN2/m3x45knZf1j2c\nHaIZKzPjkksuoa6uTvXYRaQglEJ5AT2ex33xOH/729/yK0ZWkkzuy7qSCwf42lr/P5VqsotIJYls\nkO/xPP7x4x9n986dbHeODcCVpK8k+WKGz5g6dSpLlixRPXYRqViRC/Ke53F7ayuTfvIT3gZcwdlF\nTqm9+GQlyZX4JQtGU5qjiERBZIK853msWrWKHTt2sAC4NuW5QfzFTpl68RqKEZGoisTEa2dnJ0uC\nHZzmA/8OPxUS4KP4exYmFzmlrmL9azNue+wxbr3ttoK0S0QkX1U98ZosDrZ582aGhoa4Cj/AXwQk\naz8me++bgmOtwABwuqmJjsce471jKBcsIlKpKrYn73keixYtoru7e/jYe4GF+LVn7ufsStbUXPhW\n4K++8x313kWkIlRlnnxyI479QQXJpKuAJxlZB34vcIKgsuS0adz7wx+q9y4iFaPqhmtSx99TzQSu\nwd/FKVmDZiJwEjhz441889FHNZEqIlWnonrynufR0tLCwMDAiOOLgXr8TT/e4OxGH0PAwPLl3K/d\nlkSkQlVVT76tre2cAH8tcB3++PvE4Ks95fn2PHZrEhGpdBUV5FMnWQE+BLwDPwc+hj+xOjHl+UEg\n1thYrOaJiJSd2FhfaGYPmdlBM+tKOTbFzLaZ2S/N7MdmdnFhmukP1ex+9VWuAN6Nn0WzCX9CdQg/\noya5ipXg39UNDdze0VGoJomIlL0xB3ngW8D7Rx1bC+xwzr0D2Al8PqyGjfaRFSuYd+wYC/FLFVyP\n32u/HTiIP+G6En816z8At02YwMotW5ilyVYRqWJjDvLOuWeBP4w6vBx4OHj8MHBLSO3C8zxaW1u5\n4YYb+NOpU6l9/nneA8wAWhi52fYX8Mfk7wZera/nxLJl/MuePcxVqqSIVLl8x+Qvdc4dBHDOHTCz\nS0No0zkLnZIrWZO/kWKcHZ5Zjx/o78Efnmnv6lLvXUQkEPbE63lzJNetWzf8OB6PE4/H076ura1t\nxCTrDEaWKriNswueNuBvtv1ybS1tGp4RkQqXSCRIJBKhfV5WefJmNgvY7JxrCb7fC8SdcwfNbAbw\ntHPumgzvHXOe/Pz580dc5HzgvzOyVEFywdMp4NUJE1i7dauGZ0QkcvLNk89m4hXAgq+kJ/HnPgE+\nQv5bpAIwefLk4cdX4Qf4h4A78fdl/QXwX4CfT56s8XcRkfMYc0/ezB4D4sBU/ISWduBx4H8Cl+Fv\nl3qrc+6PGd4/pp6853nMmzeP3t5ewC869kzw4Zvw0yVjwL7Zs/nWz38+praLiFSqoq14dc59OMNT\nC3M9eTpr1qwZDvAwMosmuZJ1ENhw7bXnvllEREYoq9o1nucxe/Zsjh8/DoysKpm62YeyaESkWkSq\nds2aNWuGAzz4WTXXMHLLvhgw+V3vUoAXERmDsgnynuexbdu2EccOkGGoprm5uI0TEalQ2WbXFExb\nW9uIXjz4efErGVmPpr25WfVoRETGqGx68n19fSO+/xB+Ns0bwJeB34wbR+PixazeuFFDNSIiY1Q2\nQT41N34BfoBP1ofvAAaPH2fDpEkK8CIiWSiL4RrP89i1axfgB/gbGVkXnuD7of7+IrdMRKSylUWQ\nb2tro7e3l/n4AT6ZG59KG4CIiGSvLIJ8X1/fcKXJi9AGICIiYSmLMfmmpiZOczbAq8KkiEg4ymLF\nq+d5/OW73sV3jh4dDvDJCpMv19bS9tRTKkAmIlWp2FUoQ+d5Hn+7ahVvO3p0eAu/1BLCCvAiIrkr\n2XCN53msWrWKHTt2MAd4Fj8nfhN++QIDrlq0SAFeRCQPJQnyqeWEF+PXLk7mxLenvK79yJFSNE9E\nJDKKHuQ9z2PBggX09vZyFXBdcHyQkbnxSpkUEclfUcfkkxt079u3D/CrTMaAj3FuyuQnamqUMiki\nkqei9uTXrFkzYoPuA/jj79MYWU54CJgYjytlUkQkT0VNoRw3btw5lSb/AngHfn2a5KYgn2ts5HPP\nPqsgLyJVL98UyqIG+dHHPgR8EXgQ2AcQi9G0YAGfffBBBXgREcpkZygz2wccxh9pOeWcu+7874D5\nnK00eX9wbHBoiA3TpyvAi4iEJKwx+SEg7pz7w1henKxTo0qTIiKFFVZ2jWXzWVeiSpMiIsUQVpB3\nwHYze8nMPn6+Fy7Hz41XpUkRkcILa7hmrnPud2b2p/jBfq9z7tnRL5oCzAbeAtYFX8lKk10NDdyl\nSpMiUuUSiQSJRCK0zws9u8bM2oG3nHP3jTru/hL4X8H3PZytU7N7+nTuf/55BXgRkVFKXoXSzCaY\nWUPweCKwGNiT7rXv5OzwzCz84ZrPAX++cKECvIhIAYQxJj8deNbMdgEvAJudc9vSvTBd+YLV48dr\nHF5EpEDyHpN3znnAnLG8Nl35gosXL1YvXkSkQIq64nU58Chnyxe0Nzezevt2BXkRkQwqqqzBjTfc\nQMOBA8yZMYOJzc3c3tGhAC8ich4VFeSLdS4Rkagoi9o1Y7W+tZWhvj5iTU3qxYuIFEFRe/IDaDxe\nRCQbJc+Tz8bElH/Xd3ezqa2tmKcXEak6RQ3yqVRxUkSk8EoW5FVxUkSk8Ioa5FNXurYHKZQiIlI4\nRc2u2bBiBUP9/cQaG1mt7BoRkYJTnryISBmrqOwaEREpLgV5EZEIU5AXEYkwBXkRkQhTkBcRiTAF\neRGRCFOQFxGJsFCCvJktMbPXzOxXZnZXGJ8pIiL5y3sxlJnFgF8B7wP6gZeA25xzr416nRZDiYhk\nqRwWQ10H/No51+OcOwV8F1gewueKiEiewgjyTUBvyvf7g2MiIlJimngVEYmwMKpQ9gGXp3w/Mzh2\njnXr1g0/jsfjxOPxEE4vIhIdiUSCRCIR2ueFMfFaA/wSf+L1d8BPgb9yzu0d9TpNvIqIZCnfide8\ne/LOuTNmdgewDX/456HRAV5EREpD9eRFRMpYOaRQjln7/Pmsb22lx/OKeVoRkapV3J48Z/d3Xb19\nu7b/ExG5gIrqyQNMBNZ3d7Opra3YpxYRqTolyZOfCAz195fi1CIiVaUkQX4QiDU2luLUIiJVpehB\nPjkmf3tHR7FPLSJSdYo68fpf588n1tjI7R0dmnQVERmDfCdelScvIlLGKi67RkREikdBXkQkwhTk\nRUQiTEFeRCTCFORFRCJMQV5EJMIU5EVEIkxBXkQkwhTkRUQiTEFeRCTCFORFRCIsryBvZu1mtt/M\nfhZ8LQmrYSIikr8wevL3OefeE3z9awifV5ESiUSpm1BQUb6+KF8b6PqqXRhBPufqaFES9R+0KF9f\nlK8NdH3VLowgf4eZvWJm3zCzi0P4PBERCckFg7yZbTezrpSv3cG//wH4GnClc24OcAC4r9ANFhGR\nsQtt0xAzmwVsds61ZHheO4aIiOQgn01DavM5sZnNcM4dCL79ILAn02vzaaSIiOQmryAPfMXM5gBD\nwD7gE3m3SEREQlO0PV5FRKT4Cr7i1cyWmNlrZvYrM7ur0OcrBjPbZ2avmtkuM/tpcGyKmW0zs1+a\n2Y8rKdPIzB4ys4Nm1pVyLOP1mNnnzezXZrbXzBaXptVjl+H6Mi7kq6TrM7OZZrbTzH4eJEV8Mjge\nifuX5vpWB8ejcv/qzezFIJbsNrP24Hh49885V7Av/F8irwOzgIuAV4CrC3nOYnwBvwGmjDp2L/C5\n4PFdwJdL3c4srue9wByg60LXA8wGduEP9V0R3F8r9TXkcH3twKfTvPaaSro+YAYwJ3jcAPwSuDoq\n9+881xeJ+xe0eULwbw3wAnBdmPev0D3564BfO+d6nHOngO8Cywt8zmIwzv0raDnwcPD4YeCWorYo\nD865Z4E/jDqc6XqWAd91zp12zu0Dfo1/n8tWhuuD9Av5llNB1+ecO+CceyV4PADsBWYSkfuX4fqa\ngqcr/v4BOOeOBg/r8YO3I8T7V+gg3wT0pny/n7M3qJI5YLuZvWRmHwuOTXfOHQT/BxO4tGStC8el\nGa5n9D3to3LvabqFfBV7fWZ2Bf5fLC+Q+ecxCtf3YnAoEvfPzGJmtgt/rdF259xLhHj/VIUyN3Od\nc+8B/gL4OzO7CT/wp4rajHbUrmf0Qr5/KnF78mJmDcAPgE8FPd5I/Tymub7I3D/n3JBz7lr8v8Cu\nM7N3EuL9K3SQ7wMuT/l+ZnCsojnnfhf8+3vgcfw/lw6a2XTw1w8Ah0rXwlBkup4+4LKU11XkPXXO\n/d4Fg5zA1zn7J2/FXZ+Z1eIHwP/hnHsiOByZ+5fu+qJ0/5Kcc0eABLCEEO9foYP8S8BVZjbLzOqA\n24AnC3zOgjKzCUGvAjObCCwGduNf1+3Byz4CPJH2A8qXMXKMM9P1PAncZmZ1ZvZ24Crgp8VqZB5G\nXF/wP05S6kK+Sry+bwK/cM59NeVYlO7fOdcXlftnZtOSQ01mNh5YhD/vEN79K8LM8RL8GfFfA2tL\nPZMdwvW8HT9LaBd+cF8bHP8TYEdwrduAS0rd1iyu6TGgHzgB/Bb4G2BKpusBPo8/q78XWFzq9ud4\nfd8GuoJuF5PPAAAAY0lEQVR7+Tj+GGjFXR8wFziT8jP5s+D/uYw/jxG5vqjcv3cH1/RKcD1fCI6H\ndv+0GEpEJMI08SoiEmEK8iIiEaYgLyISYQryIiIRpiAvIhJhCvIiIhGmIC8iEmEK8iIiEfb/AdZG\nyxJy6H5TAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fbfe0f20a90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "% matplotlib inline\n",
    "plt.plot(np.unique(Z_radolan_ipoli_std), 'ok')\n",
    "plt.plot(np.unique(Z_boxpol_ipoli_std),'or')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG3pJREFUeJzt3W9wXGd1x/HvUQXBiRo3kEEZK7K1qFNK+NOYDukft2ZV\nHBpoJ2H6IiWVKGoA86KETIAOiYPGMqKdpk3DpCm8cMElwaKTQkucdErBGXWbtjOClDhNqBPakaVY\nyI0NTZtYhjJVffpirzarq6vdu6u9e+/e/X1mNNm9uqs9kezjR+c5z/OYuyMiIp2vJ+0ARESkNZTQ\nRURyQgldRCQnlNBFRHJCCV1EJCeU0EVEciJ2QjezHjN7zMweDJ5fYmZfM7Nvm9lXzWxrcmGKiEg9\njYzQbwaOVz2/FXjY3V8NzAC3tTIwERFpTKyEbmaXA28HPlN1+Trg3uDxvcA7WhuaiIg0Iu4I/ZPA\n7wDVy0r73f00gLs/C7yyxbGJiEgD6iZ0M/sV4LS7Pw5YjVu1h4CISIp6Y9yzC7jWzN4ObAF+1Mw+\nDzxrZv3uftrMLgPORL3YzJToRUSa4O61BtHr1B2hu/s+d9/u7q8C3gnMuPu7gIeA8eC2dwNHanyN\nzH/s378/9RgUp2JUnIpz9aMZm+lD/33gajP7NvCW4LmIiKQkTsmlwt3/Hvj74PFzwJ4kghIRkcZp\npWigWCymHUIsirN1OiFGUJyt1ilxNsOardXEfgMzT/o9RETyxszwVk+KiohIZ1BCFxHJCSV0EZGc\nUEIXEckJJXQRkZxQQhcRyQkldBGRnFBCFxHJCSV0EZGcUEIXEckJJXQRkZxQQhcRyQkldBGRnFBC\nFxHJCSV0EZGcUEIXEckJJXQRkZxQQhcRyQkldBGRnFBCFxHJCSV0EZGcUEIXEckJJXQRkZxQQhcR\nyQkldBGRnFBCFxHJCSV0EZGcUEIXEckJJXQRkZxQQhcRyQkldBGRnFBCFxHJCSV0EZGcUEIXEcmJ\nugndzC4ws6+b2TEze9LM9gfX95vZd8zsseDjmuTDFRGRjZi717/J7EJ3/76Z/QjwT8AHgbcBZ939\nrjqv9TjvIQIwPz/PxMQES0tLDAwMMDU1RaFQSDsskbYzM9zdGnlNb5yb3P37wcMLgtesZuiG3kyk\nlvn5ea6++mrm5uYq12ZnZzl69KiSukgMsWroZtZjZseAZ4Gj7v5o8KkPmNnjZvYZM9uaWJTSFSYm\nJtYkc4C5uTkmJiZSikiks8QdoZ8HdprZxcCXzewK4NPAx93dzewTwF3Ae6JePzk5WXlcLBYpFoub\nDFvyaGlpKfL6qVOn2hyJSPuVSiVKpdKmvkasGvqaF5hNAOeqa+dmtgN4yN3fEHG/augSy9jYGNPT\n0+uuj46Ocvjw4RQiEklPMzX0OF0ul66WU8xsC3A18LSZXVZ1268B32rkjUXCpqamGB4eXnNteHiY\nqamplCIS6Sx1R+hm9nrgXsrJvwe4391/18zuA64EzgMLwPvd/XTE6zVCl9hWu1xOnTrFtm3b1OUi\nXauZEXrDJZdGKaGLiDQukZKLiIh0BiV0EZGcUEIXEckJJXQRkZxQQhcRyQkldBGRnFBCl1TNz88z\nNjbGyMgIY2NjzM/Ppx2SSMdSH7qkJmp3xeHhYe2uKIL60KXDaHdFkdZSQpfUaHdFkdZSQpfUDAwM\nRF7ftm1bmyMRyQfV0CU1qqGLbEybc0nH0e6KItGU0EVEckJdLiIiXUwJXUQkJ5TQRURyQgldRCQn\nlNBFRHJCCV1EJCeU0EVEckIJXUQkJ5TQRURyQgld2koHWogkR0v/pW20GZdIfFr6L5mmAy1EkqWE\nLm2jAy1EkqWELm2jAy1EkqWELomqngRdXl5m+/btaz4/PDzM1NRUStGJ5IsmRSUxUZOgg4OD7Ny5\nk7Nnz+pAC5EampkU7U0qGJGoSdDFxUV2797NkSNHUopKJL9UcpHEaBJUpL00QpeWWj0jdGlpiYWF\nhch7NAkqkgzV0KVlomrmvb29rKysVJ5rIZFIPKqhS6qiauYrKysMDQ1RKBQ0CSqSsLoJ3cwuAB4B\nXhrc/yV3P2BmlwD3AzuABeB6d38+wVgl4zaqmRcKBWZmZtocjUj3qTsp6u4/BEbcfSdwJfA2M7sK\nuBV42N1fDcwAtyUaqWSeFg6JpCtWl4u7fz94eAHlUboD1wH3BtfvBd7R8uiko0xNTTE8PLzmmhYO\nibRPrElRM+sBvgkMA59y99vM7L/c/ZKqe55z95dHvFaTol1ktcvl1KlTqpmLbEIzk6INdbmY2cXA\nl4EPAv9QncDN7D/d/RURr1FCFxFpUOJdLu7+gpmVgGuA02bW7+6nzewy4MxGr5ucnKw8LhaLFIvF\nRt5WRCT3SqUSpVJpU1+j7gjdzC4F/tfdnzezLcBXgd8H3gw85+53mNlHgUvc/daI12uEnlPVi4gG\nBgZUXhFpoURKLmb2esqTnj3Bx/3u/rtm9nLgL4BB4BnKbYv/HfF6JfQc0ulDIslKvIbeDCX0fBob\nG2N6enrd9dHRUQ4fPpxCRCL5oiPopG208ZZI9iihS1O0iEgke1Rykaaohi6SLNXQpa20iEgkOUro\nIiI5oUlREZEupoQusc3PzzM2NsbIyAhjY2PMz8+nHZKIVFHJRWLRJKhIe6nkIomJOo1obm6OiYmJ\nlCISkTAldIlFC4lEsk8JXWLRQiKR7FMNXWJRDV2kvdSHLonSQiKR9lFCFxHJCXW5iIh0MSV02ZAW\nEol0FpVcJJImQUXSpZKLtIwWEol0HiV0iaSFRCKdRwldImkhkUjnUQ1dIqmGLpIu9aFLS2khkUh6\nlNBFRHJCXS4iIl1MCV0qtJBIpLOp5CKAJkFFskYlF2maFhKJdD4ldAG0kEgkD3rTDkDSs9qWuLS0\nxMLCQuQ9Wkgk0jlUQ+9SUTXz3t5eVlZWKs9VQxdJTzM1dI3Qu1RUzXxlZYWhoSEKhYIWEok04Zn5\neT43McH5pSV6BgYYn5piRxv/Dimhd6mNauaFQoGZmZk2RyPS+Z6Zn+eeq6/mwNwcFwHngP2zs9x0\n9GjbkromRbuUNt8Saa3PTUxUkjnARcCBuTk+18ZOMY3Qu0T1BOjAwAB79+5ldnZ2Xd/51NRUilGK\ndK7zS0uVZL7qIuB8GzvFlNC7QNQE6OzsLIcOHeLgwYPafEukBXoGBjgHa5L6OaCnjb/1qsulC4yN\njTE9Pb3u+ujoKIcPH04hIpH8iayhDw83XUNPpMvFzC4H7gP6gfPAQXe/x8z2A+8DzgS37nP3v20w\nZmkDLRoSSd6OQoGbjh7lzokJzp86Rc+2bdyUwS6XFeBD7v64mfUB3zSzo8Hn7nL3u5ILT5qlRUMi\n7bejUGB/ir/11k3o7v4s8GzweNnMngJWWyQa+nVA2iPuoiFNgIrkS0Nti2Y2BFwJfD249AEze9zM\nPmNmW1scmzSp1qKhkZERRkdHtQJUJIdid7kE5ZYvATcHI/VPAx93dzezTwB3Ae9JKE5pgBYNiXSn\nWAndzHopJ/PPu/sRAHf/btUtfwo8tNHrJycnK4+LxSLFYrGJUCUuLRoS6TylUolSqbSprxGrbdHM\n7gO+5+4fqrp2WVBfx8xuAd7k7r8R8Vq1LbZB9STo1q1bOXbsGCdPnqx8XhttiXSWRA6JNrNdwCPA\nk4AHH/uA36BcTz8PLADvd/fTEa9XQk9Y1CTo4OAgO3fu5OzZs1o0JNKBEknom6WEnjwtHBLJHx1B\n16W0cEhEQAk9FzQJKiKgkksuRNXQNQkq0tlUQ+9iq10u2jlRJB+U0EVE2ijJI+eU0LtI+MAKjchF\nWisqWQOVay9s3cr5Y8f4vZMnW7JdbpgSeo5p4ZBI+0TtbX7L4CBbzCoJfAK4lfUHWtw5OtqSHRcT\n2Q9d0hc16Rk2NzfHxMSE+s5FWiDqfND+xcU1CbwHUj9yLkwJvQNE7Z4YRX3nIs2rLrHMHz++LlmH\nE3gPpH7kXJgSegfYaOFQmPrORZoTLrFMsD5Znw9dGw/umwquVWroKZ4zoISeQeEJz4svvrjua3Rg\nhUjzwiWW97I+WZ8eHGRfVQ39UmB5cJDJnTu58OzZVI6cC9OkaMZstNGWma2ZBNXmWyKbEy6x3Hfm\nzNrPA7f09/P6K66gZ9u2tV0uwZmhrWxTDFOXSw5stNHWddddR19fnxYOibRAVIklyY6VZqjLJQc2\nqpe/8MILPPDAA22ORiSf4pRY0q6HN0MJPWO00ZZIMmp1sewAbgbeVVViSbse3gwl9AwILxravn37\nukVDmvAUaV6cLpZLgZ/asye1EksrqIaeMp02JJKM6hH5txYWuG9hoZLAnwHuJqLE0qJl+62gGnoH\nilo0tLi4yO7duzly5EhKUYl0nuoEHt5n5WOQyxJLmBJ6ynTakMjm1etaeQn5LLGEKaGnoLpmvrCw\nEHmPJkFFXhTe+XDP3r08fPBg5fny2bNrulbCy/THgf3AATq7i6UeJfQ2i6qZ9/b2srKyUnmuSVCR\nF4VH308Bd9x/P59aWakk570ve1nNfVZ2AO8BfnNoiNcVCrkpsYQpobdZVM18ZWWFoaEhCoWCJkFF\nQsI9438BlWQO5aT9qv/5n7r7rHx2eJi7MjTpmQQl9DbbqGZeKBSYmZlpczQi2VSrZ/w867etfS9w\n05Yt3PODH2R2n5V2UEJvMy0cEqmtXs941La1lwJb3/pW7uzrq+yzcnsXJPAw9aG3gU4bEonvwNgY\nH5me3rBn/Cngjt7eNTX0rPWQt4I258ogLRwSqa/RnQ8rXS5t2PUwLVpYlEFaOCRSW7PL8nft3t3e\nQDtAT9oB5J0WDonUttHOh+eC56sllXG18talEXrCNAkqUtv5paWuWJbfDkroCZuammJ2dnZN2UUL\nh0Re1DMw0BXL8ttBk6JtsNrlotOGpBuFl+2vOcotYiOtvHatNEpdLiKSKeEJz3PALYODbKk6bHn1\n2taqRUB57FpplBK6iLRdrY2zwvuQQzbP78witS1mRPVCooGBAZVYJLfqbZwV3occ1u+ESPD8vDq/\nNk0JvcWiFhLNzs5qJajkRq2TgMIbZ0XtQ34+4to5oEedX5umPvQWi1pINDc3x8TEREoRibTO6oj8\nI9PTHCiVeE2onBLeOGuc8j7k1T3lpwcH2bd9u/rME1B3hG5mlwP3Af2Uf15/6u5/bGaXAPdTbhtd\nAK539+cTjLUjaCGR5E2tEXl4BB5nH/Lbg8R958REZem++sxbI07JZQX4kLs/bmZ9wDfN7GvAbwEP\nu/sfmNlHgdsoz3V0NS0kkk4WNcH55RtvrNTIwzXxcdaeBHQ98NuhjbM22odcE6Ct13CXi5k9APxJ\n8PFmdz9tZpcBJXf/yYj7u6rLJaqGrt0UJStqdaRE9YPf0NfHny8vV5L4AeAjrE3qTwEfqxqBd8PG\nWe2QeNuimQ0BJeB1wKK7X1L1uefc/eURr+mqhA5aSCTZFNmRUjWajmon/BjwieqvAdxDxNmcXb4I\nKAmJti0G5ZYvATe7+7KZhbP0hll7cnKy8rhYLFIsFhuJseMUCgUO69dJyZh6R7lFtROGa+TdcjZn\nGkqlEqVSaVNfI9YI3cx6gb8GvuLudwfXngKKVSWXv3P310S8NvcjdPWdSyfYPzLCgaqEsVr7XrVR\nOeUP+/q4Jyi7aETePkmO0A8Bx1eTeeBBynMidwDvBrpic+9w8t67dy833nij+s4l88KbYIU7UsaJ\nPlj5PYcOcWdVTVwj8uyqO0I3s13AI8CTlMsqDuwDvkH5t7ZByqW16939vyNen5sRetSEZ19fH8vL\ny+vuHR0dVdlFUlc9CRqe9Iw6yk17qmSH9nJJQPWIfGFhgYWFhVivGxkZYWZmJtngRGrYaGOs6oSt\njpTs0l4uLRY1Io9LfeeStvAk6EXAJxcXuXP3bvZXHX+oo9zyQwk9pJkRebjsogMsJAvCJwGBNsHK\nOyX0Ks2MyIeHhzl06BAHDx5U37mkLrxMX5tgdRfV0KuMjY0xPT1d976hoSEKhYKSt6Su0UlPtRx2\nDk2KNqG6xHL8+HHOnDlT834t45e01NtnJWqlZ3hZviY9O4cmRWOoTuBbt27l2LFjnDx5suZrNCKX\ntEV1rNxw5MiafVaiVnq+BnhdocABdVx1hVwn9DiLgOrRiFzapZGj3C4C3lCVzGH9QiFQzbzruHui\nH+W3aL8TJ0748PDw6kIoB7yvr2/N840++vv7fWRkxEdHR/3EiROpxC/dZeHECf/w8LAvgzv4cfB3\n9/ZWnt8e/Lf6YxIqn3fwBfBbqq4tg394eNgX9Ge4IwW5s6F8m9sRetTJQVErOqPs2bNHqzwlcZs9\nyu164KaqfVYuBZYHB5msWjikZfrdJbcJfaOTg+pRD7kkpVZHSvjgiI2OcqvetjZqn5XblcC7Wq4S\nenhRUJTwIqDBwUF27tzJ2bNnNekpLdNoR0ozR7mtjr610lNW5SahRy0K6u3tZWVlpfJci4CkVRo+\n+adOR8o4zR/lJrIqN33oGy0KUsuhtForTv7RUW5ST9f1oYcXBUUpFAra9VA2rZEJzDgn/4wTvfd4\neASucoo0oqMSejOLgrTrocTRSAml3gRmVD+4OlKkHTomoTe7cZY6ViSs3oTlU8Ad99+/YQml3gTm\nOPFO/lFHirRax9TQ426c1d/fzxVXXKGauVTUahc8B9zQ17dmwjJc3w6fvfkMa0++18k/koRc19Dj\n9pVrUZDUSuDh0XbUEvp6JZSoFsL37d2r0bekrmMS+sDAQN17VGLpTo0k8DgTlnFLKJrAlKzJdEIP\nT4Ju3759zSSoFgV1p80k8DgTluEecE1gSqfIbEKPmgQdHBzk2muvVQLvIo2uuKyXwMeJN2GpEop0\nosxOim40CTo6OqoaeY41OoEZnrAMT2g+A9zN2gSuCUvpBB0/KRpnodApHXDb0Rrp944zgVlvBB5V\nLtFoW/IqtYTe7OETWijUuSKXzNfo925mxaUSuHSzVBJ6VH38yJEjdfcrVxdLumqNrnsGBhgPfjZx\nT92pt2S+2RWXSuDSrVJJ6I0cPqGFQu3TSDkkPLo+B9zyyCNsMdvwnkaXzI+jFZcijUgloTdy+IQW\nCsVTb/RcbzRdL2GHyyHh0fVFQP/iYs17Gu33rtUuqJ5vkfVSSegbLRIKHz7RrSWWcHKuV8qol4zj\njKbrJexwOSQ8uo5zzzi19/xW+URkc9qW0OstEurmwyfqterVK2XUS8ZxRtP1knF4NB1V3z5f5x4t\nmRdJVlsSeiOLhHbn8Ffpzbbq1Stl1EvGcUbT9RL2OGvLIVEn6pweHGRf1T88cU/dUflEpDXaktCj\nJkEXFxfZvXs3R44caUcIsW22Ft1oOSROq95mR89xRtPj1E7YUeWQqNE1wJ0TExuuuNSSeZHktCWh\nbzQJmoVFQrXKHc3Uohsth8Rp1atXyhindjKOM5qOm7DDyThqdL0/NImtEbhIe7QloS8sLEReT3qR\nUL3JxXrljmZq0Y2WQ+K06tUrZcRJxnFG03ETtohkU9sSem9vLysrK5VrSXSwNDq5WK/c0UwtutFy\nyDjxVjpC7VJGnGQcZzQtIp2rbV0uKysrDA0NUSgUmu5gafXkYqMTg0mUQxpp1VMpQ0RqqZvQzeyz\nwK8Cp939DcG1/cD7gDPBbfvc/W/rfa3Lf/ADftGdHuDU4iL3tbDXupnJxSRq0c1OJmqiUEQ2q+72\nuWb2C8AycF8ooZ9197vqvoGZDwE/BUxT+wzGdZOLVfeEE3a9cx/Dnyfia1Rvrfoo8CbWb61a+Ucl\nSL7h52v+IapxT6sSdqlUolgstuRrJakT4uyEGEFxtlqnxJnI9rnu/o9mtiPq/eK+yRibXzaexORi\n9ej5G088wS/t2tWyWnRS5ZBO+cPYCXF2QoygOFutU+JsxmZq6B8ws3cB/wx82N2f3+jGViwbT2py\ncTWBT05Osn9yMt7/uYhIBjWb0D8NfNzd3cw+AdxFeVV3pKxPLoqI5EGsI+iCkstDqzX0uJ8LPp/s\nGXciIjmV1BF0RlXN3Mwuc/dng6e/BnyrVQGJiEhz4rQtfgEoAq8ws5OUG0pGzOxKypWSBeD9CcYo\nIiIxxCq5iIhI9vUk9YXN7Boze9rM/s3MPprU+zTDzD5rZqfN7Imqa5eY2dfM7Ntm9lUz25pyjJeb\n2YyZ/auZPWlmH8xonBeY2dfN7FgQ5/4sxhnE1GNmj5nZg1mNEcDMFszsX4Lv6TeCa5mK1cy2mtkX\nzeyp4M/oz2Qwxp8IvoePBf993sw+mLU4g1hvMbNvmdkTZjZtZi9tJs5EErqZ9QB/Avwy8FrgBjP7\nySTeq0l/Rjm2arcCD7v7q4EZ4La2R7XWCvAhd38t8HPAbwffw0zF6e4/BEbcfSdwJfA2M7uKjMUZ\nuBk4XvU8izFCuZRZdPed7n5VcC1rsd4N/I27v4byusGnyViM7v5vwffwjcBPU26c+zIZi9PMtgE3\nAW8Mmkt6gRtoJk53b/kH8LPAV6qe3wp8NIn32kSMO4Anqp4/DfQHjy8Dnk47xlC8DwB7shwncCHl\ndQlvylqcwOXAUcrzQQ9m+WcOzAOvCF3LTKzAxcBcxPXMxBgR21uBf8hinMA2ygvXLwmS+YPN/l1P\nquQyACxWPf9OcC3LXunupwG83MHzypTjqTCzIcqj31nKP+BMxRmUMo4BzwJH3f1RshfnJ4HfAaon\njbIW4yoHjprZo2b23uBalmItAN8zsz8LyhkHzezCjMUY9uvAF4LHmYrT3U8BfwScBJaA5939YZqI\nM7Eaeg5kYrbYzPqALwE3u/sy6+NKPU53P+/lksvlwFVm9loyFKeZ/QrlzeUep/aWFal/LwO7vFwm\neDvlUtsvkqHvJ+VR5BuBTwVxnqP8W3iWYqwws5cA1wJfDC5lKk4z+zHgOspVg23ARWY2GhFX3TiT\nSuhLwPaq55cH17LstJn1Q7nPnhd3kkyNmfVSTuafd/fVs/oyF+cqd38BKAHXkK04dwHXmtkJ4M+B\nXzKzzwPPZijGCnf/j+C/36VcaruKbH0/vwMsuvs/B8//knKCz1KM1d4GfNPdvxc8z1qce4AT7v6c\nu/8f5Tr/z9NEnEkl9EeBHzezHWb2UuCdlOtCWbJmsRTl+MaDx+8GsnDY6SHguLvfXXUtU3Ga2aWr\ns+9mtgW4mvJmmZmJ0933uft2d38V5T+LM+7+LuAhMhLjKjO7MPitDDO7iHLt90my9f08DSya2U8E\nl94C/CsZijHkBsr/kK/KWpwngZ81s5eZmVH+fh6nmTgTLPRfA3wb+Hfg1jQnHSJi+wJwCvhh8M38\nLcoTEg8HMX8N+LGUY9wF/B/wOHAMeCz4nr48Y3G+PojtceAJ4PbgeqbirIr3zbw4KZq5GCnXp1d/\n5k+u/t3JWqyUO1seDWL9K2Br1mIM4rwQ+C7wo1XXshjnfsoDoSeAe4GXNBOnFhaJiOSEJkVFRHJC\nCV1EJCeU0EVEckIJXUQkJ5TQRURyQgldRCQnlNBFRHJCCV1EJCf+H4gJK8HdU3pXAAAAAElFTkSu\nQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fbfe111b750>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.unique(Z_radolan_ipoli_median), 'ok')\n",
    "plt.plot(np.unique(Z_boxpol_ipoli_median),'or')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHJdJREFUeJzt3X2QXOV15/HvGcuxhcaWZRMPHjGIdqcSgzHGTsXrjRI8\nk4hd4mzJm6SKDZ7GnoDBlbJlVi5SvKkzmmpXKi8sLlbx/kEM4U12Qbwhgio7Zih54qyr5BAjBYiA\n7A4tkKVF2EvWRuPEFUUnf/Tt1p07d/pt+nbfe/v3qZpi7p1uzWEER1fnOed5zN0REZHsGxl0ACIi\n0htK6CIiOaGELiKSE0roIiI5oYQuIpITSugiIjnRdkI3sxEze9LMHgmuN5nZY2b2vJl93cw2Jhem\niIi00skT+vXA4dD1TcDj7v4zwH7g5l4GJiIinWkroZvZucCHgS+Gbn8EuDf4/F7gP/c2NBER6US7\nT+ifB34HCI+Vjrn7CQB3fxl4e49jExGRDrRM6Gb2q8AJdz8EWJOXag8BEZEBWtfGa7YC283sw8B6\n4E1mdj/wspmNufsJMzsHeCXuzWamRC8i0gV3b/YQvULLJ3R3v8Xdz3P3dwK/Cex396uAR4GZ4GUf\nB/Y1+TVS/zE7OzvwGBSnYlScirP+0Y219KH/PnCZmT0P/HJwLSIiA9JOyaXB3f8K+Kvg81eBbUkE\nJSIindOkaGBycnLQIbRFcfZOFmIExdlrWYmzG9Ztrabtb2DmSX8PEZG8MTO814uiIiKSDUroIiI5\noYQuIpITSugiIjmhhC4ikhNK6CIiOaGELiKSE0roIiI5oYQuIpITSugiIjmhhC4ikhNK6CIiOaGE\nLiKSE0roGVWtVimVSkxNTVEqlahWq4MOSUQGTNvnZlC1WuWyyy5jcXGxca9YLDI/P0+hUBhgZCLS\nK9o+d0iUy+VlyRxgcXGRcrk8oIhEJA2U0DPo2LFjsfePHz/e50hEJE2U0DNo8+bNsffHx8f7HImI\npIlq6BmkGrpI/nVTQ1dCz6hqtUq5XOb48eOMj49TqVSUzEVyRAldRCQn1OUiIjLElNAzQoNEItKK\nSi4ZoEVQkeGjkktOaZBIRNqhhJ4BGiQSkXYooWeABolEpB2qoWeAaugiw0d96DmmQSKR4aKELiKS\nE+pyEREZYkroKaQhIhHphkouKaMFUBEBlVxyQUNEItKtlgndzN5gZt82s4Nm9rSZzQb3Z83su2b2\nZPBxefLh5p+GiESkW+tavcDdf2xmU+7+IzN7HfAtM/ta8OXb3f32ZEMcLhoiEpFutVVycfcfBZ++\ngdofAvWieEf1HWmtUqlQLBaX3SsWi1QqlQFFJCJZ0VZCN7MRMzsIvAzMu/sTwZc+bWaHzOyLZrYx\nsSiHSKFQYH5+nunpaaamppientaCqIi0paMuFzN7M/AwsAP4HvB9d3cz+xzwDne/JuY9Pjs727ie\nnJxkcnJyrXGLiOTKwsICCwsLjeu5ubnkJ0XNrAwshWvnZrYFeNTdL455vdoWRUQ6lEjbopmdXS+n\nmNl64DLgOTM7J/SyXwee6eQbi4hIb7XscgHeAdxrZiPU/gB40N2/amb3mdklwGngCPDJ5MIUEZFW\nNCkqIpJCmhQVERliSugiIjmhhC4ikhNK6CIiOaGELiKSE0roIiI5oYSeAjqhSER6QX3oA6YTikQk\njvrQM0gnFIlIryihD5hOKBKRXlFCHzCdUCQivaIa+oCphi4icbqpoSuhp0C1WqVcLnP8+HHGx8ep\nVCpK5iJDTgldRCQn1OUiIjLElNAHQINEIpIElVz6TIugItIOlVwyQINEIpIUJfQ+0yCRiCRFCb3P\nNEgkIklRDb3PVEMXkXaoDz0jNEgkIq0ooYuI5IS6XFJKfeci0g96Qk+YauYi0g09oaeQ+s5FpF/W\nDTqAvFPfucjweLFa5Z5ymdPHjjGyeTMzlQpb+vg3cSX0hKnvXGQ4vFitsueyy5hbXGQDsATMHjjA\njvn5viV1lVwSVqlUKBaLy+4Vi0UqlcqAIhKRJNxTLjeSOcAGYG5xkXv6WF7VE3rCCoUC8/Pz6jsX\nybnTx441knndBuB0H8urSuh9UCgUeOCBBwYdhoj0WLhm/syRIyzBsqS+BIz0sbyqhC4i0oVozfxZ\n4FPr1vGFU6fO1NCLRXb0sbyqPvQE1Ef7jx07xubNm1ViEcmhuVKJG/buXfZE/iyw6/zzuahQYGR8\nfE1dLt30oesJvcfiBokOHDigQSKRnImrmV8AXFQoMLd//yBCUpdLr2mQSGQ4jGzezFLkXr9r5lEt\nE7qZvcHMvm1mB83saTObDe5vMrPHzOx5M/u6mW1MPtz00yCRyHCYqVSYLRYbSb1eM58ZYEtyy5KL\nu//YzKbc/Udm9jrgW2b2NeA3gMfd/Q/N7EbgZuCmhONNPQ0SiQyHLYUCO+bnua1c5vTx44yMj7Oj\nz5OhUR0tiprZWcA3gd8G7gc+5O4nzOwcYMHd3xXznqFaFNVmXCLSC4nth25mI8B3gCLwBXe/2cz+\n0d03hV7zqru/Nea9Q5XQQQdYiMjaJX7AhZm9GXgY+Azw1+EEbmb/z93fFvOeoUvoIpJP/dx8K/G2\nRXf/oZktAJcDJ8xsLFRyeWW19+3evbvx+eTkJJOTk518WxGRgUt6862FhQUWFhbW9Gu0fEI3s7OB\nf3H3H5jZeuDrwO8DHwJedfc/CBZFN7n7ikXRYXhC1yCRSP7FDRItAbdNTzObwNYeST2hvwO4N6ij\njwAPuvtXzewA8JCZXQ28CFzRccQ5oEEikeGQhs23WmmnbfFp4P0x918FtiURVJY0GyTShlwi+VEf\nJBrk5lutaFJ0jTRIJDIc0jhIFKW9XNZIg0Qi+RTX0ZK2QaIo7ba4RhokEsmf2I6WYrGvx8kl3ofe\njbwndNAgkUje9LujJY62zx0QnUgkkn3hEkv18OHUd7TEUUIXkaEUTuA/3LiR0wcP8nsvvcQGoAyp\n72iJo5JLhzREJJJ90Rp5mdpWsfUE/iJwB1AJ7mWlhq4n9A5oiEgku6IHOt935EgjgY+w/Gl8C3A9\ncNXYGO+58MJUdrTEUULvgIaIRLIp+kS+i+UJfISVJZazgfdu29a3RdBeUELvgIaIRNIr2je+7brr\nePzOO2OfyF/P8gQ+Q63ssqLEkqKhoXYooXdAQ0Qi6RR9An8W+IMHH+QLp07FPpHPALPAXHD/bODk\nxAS73/c+znrttcyUWKK0KNoBDRGJpFO0b3wOuAFWvYZa0t91/vlcVCgwMj6e6N7m3dCiaMIKhQLz\n8/MaIhJJmehOiKdp/kS+BNxVLHJ7H7tW+kEJvUMaIhJJn+hOiNFFzi3ANcDHQk/kWSyptKKSi4hk\nXmwNfd26Rg19EH3ka6W9XBKgQSKRbGh0uQQ7ITa6XILrtNXIW1FC7zEtgooMRrMWxLjrrCXrdiih\n91ipVGLv3r0r7k9PT6uOLpKQVuWTPJRT2tFNQteJRU1okEik/+4plxvJHOAhaCTvuOsNwNziIveU\ny/0ONXWU0JvQIJFI/7VqQYxeQza2tu0HJfQmKpUKxWJx2b1isUglY+PAIllSb0FsXEPTa8jG1rb9\noBp6CzqNSCR5zfYmVw29g/cooYvIIMWd37lzYoKNoX1Voi2IWW9JbIcSuohkThrO70wjdbn0QLVa\npVQqMTU1RalUolqtDjokkVyLLoKCFjm7pb1cQnQikUh/RE8PyuL5nWmkkkuIBolEei9u6vPhq6/O\n1b4rSdD2uWukQSKR3opb8Lxy3z6+fPJk44n8AuDGU6dyvxNiPyihh2iQSGTtmh3GvAG4OJTM6y4A\nLioUmNu/v7/B5owSekilUuHAgQMrNuPSIJHI6pr1kEePfoOV53mCaua9ohp6hAaJRJprlsDLwE20\nPvrtj0ZH2RM8qatmHk996CKSqGhNPJrA68e8NV4P7GH50W+zxSK/dvfduR8MWistiopIoqI7IY6w\n/Om7k6Pftl56af8CHxJDn9B1IpFI+6JDQNEEPkPtqb1Cvg9jTquhTugaJBLpTPQw5hmWJ/CzgZMT\nE+wO7cOiFsT+GeoaugaJRJaLDgHNBB1eqy2Cxm2kpXp4byRSQzezc4H7gDFqe8vf6e57zGwWuBZ4\nJXjpLe7+lx3GPFAaJJJh16xjZQnY+c1vst5sRQLfvX17I4HfqgSeGu2UXE4Bn3X3Q2Y2CnzHzOaD\nr93u7rcnF16yNEgkw6xVx8oGYOzo0RX3Pn/0KLddeimz+/b1P2hpquVui+7+srsfCj4/Sa2NtJ4J\nO/rrQNroRCIZZq06Vla7p50Q06uj7XPN7HzgEuDbwa1Pm9khM/uimW3scWyJKxQKzM/PMz09zdTU\nFNPT01oQlaGxWsfKstfE3NNUZ3q1vSgalFsWgIq77zOznwS+7+5uZp8D3uHu18S8L7WLoiLDptk+\nKy8Cd7C85XDnxMSKGrqmOvsjscEiM1sHfAW43933Abj790Iv+RPg0dXev3v37sbnk5OTTE5OdhKj\niPRAtGb+LPCp0La1cS2Htwblx9vK5cZUp9oQk7GwsMDCwsKafo22ntDN7D5qT+OfDd07x91fDj7f\nCfycu3805r2pekLXIJEMq7ij3p4FdoWmONVymB5JtS1uBaaBp83sIODALcBHzewSamW2I8AnO464\nzzRIJMMmXGKpHj6sbWtzrmVCd/dvAa+L+VKmes4ByuXysmQOsLi4SLlc1iCR5EKrnRC1bW2+DdXo\nvwaJJG862cr2E6zcZ2W2WGSH2nRzY6gSugaJJEvizuJ8/M47Vz2bM5rAoz3kW4DrgavGxnjPhRdq\ngTOHhmovl7gaerFYVA1dUqHZ03bcQcpXjo4uO5szuhd53OESS8Bt09PMqsSYetoPvYX6IJFOJJK0\naTWG/xA0kjnEn83Zzla2KrHk21AldKgldS2AStq0GsM/TeuzOWfQVrbDbugSukhaNGspjD5tR68B\nrgB2hM7mXG0wSAl8eOS6hq4hIkmrViWW6Bh+XA1dZ3Pmmw6JDtECqKRNN/uohA+OaHS5KHkPBSX0\nEJ1GJGkSfSLfBXwu+hpgZ6ilUAl7uKnLJURDRJIm0UXP6IIm1Grg7922TS2F0rWO9kPPEg0RSZpE\n9x6fodY3Xt9rvF4Tn1FLoaxBbp/QK5UKBw4cWFFD12lEMggjmzcveyLfAlwDfCy006FaCmWtcltD\nhzNdLhoikqR1OqavgyKkFS2KivRJp2P6ajGUTimhi/RBqx5y7aEivdBNQs/Vomi1WqVUKjE1NUWp\nVKJarQ46JMmhbsb0NwCn1WElCcvNoqhOI5J+iXastDOmr4MkpB9y84Te7DQikV6qd6zUzXDmNCCo\n7bHyqXXr1JIofZebJ3QNEkm/zFQqzB440Ci7xG2Kde1113FbaAFULYnSD7lJ6BokkiRF2xJ/7e67\nlyXsuF0Nt1566YCilWGVmy4XbcYlvdSsLVE95NIPQ9+2qEEi6YVWbYmgNkRJ3tBvzqXTiKRbzba2\njbYlgtoQJZ1yldBF4nQ6lr8LmrYlgtoQJZ0yXXLRiUTSSrR8EjeWf+XoKF8OHbgcnfSMO3xCNXRJ\n2lCVXDRIJKtpVj55CBrJHGoJ+uJQMoczW9vOocOWJVsym9CbDRKpjj5cmnWkRMsncWP50cMm4ra2\n1WHLkgWZTegaJBJo3ZESTdZx9fArgB2jo+wJntSXgLuKRW5XSUUyJrMJXYNEwyG6oFkfn2+3I2WG\n5eWT+lh+uIZ+V7HINZFBIZVUJIsyuyiqQaJ8ajXQs3NigvVmy0oq4cOW47aufRbYFSqfNLpctC+5\npNjQDRZpkCj7miXwuIGeVnuPqyNF8mKoulxAg0RZ1EkCjxvoaVVSUUeKDLNMJ3TJllYLmNFkHbeA\neRp1pIisJlMJXYNE2dbqpJ9oAp+hlvTD5ZMTExPcEqqhqyNF5IzMJHQNEmVfq5N+ZliewOPKJ7cG\nXS63lcvqSBGJyMyiaKlUYu/evSvuT09Pq46eEXOlEjfs3dt0AXPnxAQbQwlcHSgyrBJZFDWzc4H7\ngDFqJcw/cff/bmabgAeplTGPAFe4+w86jrpNGiTKvnZO+lH9W6R77ZRcTgGfdfdDZjYKfMfMHgN+\nC3jc3f/QzG4Ebqa2xpUIDRJlQ6udDds56UdEutNxycXM/gL44+DjQ+5+wszOARbc/V0xr+9JyUWD\nROnXzs6G6gkXaU/ig0Vmdj6wAFwEHHX3TaGvverub415T88GizRIlC7Rp/GTr73G7kceWXXoB3TS\nj0i7Eh0sCsotXwGud/eTZhbN0qtm7d27dzc+n5ycZHJyspMYGzRIlB7Rp/El4Lo3vrHlzoY66Uck\n3sLCAgsLC2v6NdpK6Ga2jloyv9/d9wW3T5jZWKjk8spq7w8ndMmHaE/5BuCd//zPLXc21Ek/IvGi\nD7tzc3Md/xojbb7ubuCwu98RuvcItdZhgI8D+6JvWqtqtUqpVGJqaopSqUS1Wu31t5AOvFitMlcq\nMTs1xeL8/Iqn708AO9avZym4ru9sWL+u19DrOyaKSG+1rKGb2Vbgm8DT1MoqDtwC/A21A2AmqLUU\nX+Hu/z/m/V3V0LUIOnidbpy1BOz6yEd4y+hoo4tFOxuKdCdXuy1qkGiwWu27ol0NRZKVq90WNUjU\nf83O4ozuu7IFuB64amyM91x4oUbwRVIgtQldg0T9FX0ij57FGbfAeTbw3m3b1IIokhKpLbmohp68\nZk/kOjhCZLByVUMHDRL1UtxI/sNXX73siTx8lNuLwB7OHByhjbNE+it3CV16I24I6MrRUb4cnHIP\n7Z3FqQQu0j+5WhSV5lptghW+jpZTNgAXh5I5rDzKTQdHiGRPqp7QdSLR6pr1hEc3wYpeR8spoCdy\nkbTL9BO6TiRaXaue8Iegkbzjrl/Pyg6VK4Ado6PsCZ7U9UQukn2pSejlcnlZMgdYXFykXC7nfpCo\nVfnk5GuvNT2LM7oJVvR6hvhyyjWRvcnVRy6SbalJ6MM6SBS7h/iDDy7bQzy6i2G0J7zV9RbgGuBj\noXJKPXlvvfTSZP8FRaRvUpPQh3WQKLprYbRcEreL4QzLD1Oub4JVf1/0WuUUkeGQmkXRPA8SNSup\nVA8f5r5Xzuw8XC+NLHs/MLd+PXv+6Z9W7QmPboKlTbFEsi3zfeh5HCRqdSxbdIFztVN+orsYKkGL\n5FvmE3oezZVK3LB3b9sj9TqHU0QgY22Lee45D5dYqocPN+1Aidu18NrrrlP3iYh0bCAJPc8953E9\n462OZYvbtVDdJyLSqXaPoOupZj3nWRftWvkEZ5I66Fg2EUnOQJ7Q89xzfvrYMZVURGQgBpLQs9xz\nHm1BrD9Zh/cVV0lFRAZhIF0uWeo5b7YpVr0ffL3ZqhtlqUtFRLqRqbbFLPSct9oUi1XuaddCEVmr\nTLUtFgqF1G+6FV3gjG6Ktdq9C4CLCgXm9u9POkQRkYbU7OWSFs16yONaDk/H3FsCRjKwHiAi+dK3\nhJ6WQaJm+6pEa+TRHvIZlm+KtQScmJjgllANvVEzVxuiiPRZX2roL7zwQioWQTvdVyXupPu4g5Ih\n6HLRPisi0iOpXRSdnp5m7969K742PT3d1zp6q31VVtvpcGeoh1zJWkT6IbWLooMaJIqWV34UWuCE\nlfuqtDuWLyKSRn1J6EkNEnVSD18Crhwdbbqvygwra+Sqh4tIVmSqht5syKdVPZzgNX8UOhg5bggo\nrkauEouI9Ftqa+juvuZBolZDPu3UwwH+6wc/yKZiUSf7iEiqpbaGDt0NEoWfyJ85coT7jhxZdcin\nnXr4ErCpWFxRD9e+KiKSB6kdLIo+ke+iecJWPVxEhl1f9kMvlUpUq9WO3hMdu389Z/YUhzMJe7V9\nxs8GTk5MsHv7dmanprhtelobZIlIrvWlhg6dL4LOTk0xt7DQuH4R2EOtLr7aAqbq4SKSF6muoddP\nJGq3jj6yefOyEsoW4BrgY6FdDG+NSdiqh4vIsGr5hG5mdwH/CTjh7hcH92aBa4FXgpfd4u5/ucr7\n/aeA/wP8wtgYv3TBBSt6xuMOiojrI9e+4iIyLLp5Qsfdm34AvwBcAjwVujcLfLbVe4PX+mHwj4Kf\nBPfg+uPr1jWuT4JfOzHhnznvvBX3bti+3X93asp3T0/7kRde8KR84xvfSOzX7iXF2TtZiNFdcfZa\nVuKspefWOTb80XJR1N3/F/CPcX+AtPuHxkPAnZwpnzwEjWEegvtjR482nsbr9z5/9Cijb3oTc/v3\nM/vAA4k+mS+E6vVppjh7JwsxguLstazE2Y21dLl82swOmdkXzWxjsxdGe8Sj1/VAovc2AKdzcHC0\niEg/dJvQ/wfwTne/BHgZuL3VN1lqcg1nDooI00ERIiLta6tt0cy2AI96sCja7teCryfbFykiklOe\nUNuiEaqZm9k57v5ycPnrwDO9CkhERLrTMqGb2ZeASeBtZvYStQ6XKTO7hFql5AjwyQRjFBGRNiQ+\nKSoiIv2R2F4uZna5mT1nZv9gZjcm9X26YWZ3mdkJM3sqdG+TmT1mZs+b2ddbde70IcZzzWy/mf29\nmT1tZp9JaZxvMLNvm9nBIM7ZNMYZxDRiZk+a2SNpjRHAzI6Y2d8FP9O/Ce6lKlYz22hmf2Zmzwb/\njf67FMb408HP8Mngnz8ws8+kLc4g1p1m9oyZPWVme83sJ7qJM5GEbmYjwB8D/xF4N3Clmb0rie/V\npT+lFlvYTcDj7v4zwH7g5r5HtdwpasNb7wb+PfCp4GeYqjjd/cfAlLu/j9oA2q+Y2QdIWZyB64HD\noes0xgi1Uuaku7/P3T8Q3EtbrHcAX3X3C4D3As+Rshjd/R+Cn+H7gZ+l1jj3MCmL08zGgR3A+4Pm\nknXAlXQTZ6eTSO18AB8Evha6vgm4MYnvtYYYt7B8+vU5YCz4/BzguUHHGIn3L4BtaY4TOAv4W+Dn\n0hYncC4wT2096JE0/54DVeBtkXupiRV4M7AYcz81McbE9h+Av05jnMA4tf0HNwXJ/JFu/19PquSy\nGTgauv5ucC/N3u7uJwC81sHz9gHH02Bm51N7+j1A7Tc4VXEGpYyD1GYS5t39CdIX5+eB3wHCi0Zp\ni7HOgXkze8LMPhHcS1OsBeD7ZvanQTnjTjM7K2UxRv0X4EvB56mK092PA/8NeAk4BvzA3R+nizj7\nsh96RqVitdjMRoGvANe7+0lWxjXwON39tNdKLucCHzCzd5OiOM3sV6ltLneI5ltWDPxnGdjqtTLB\nh6mV2n6RFP08qT1Fvh/4QhDnErW/hacpxgYzez2wHfiz4Faq4jSztwAfoVY1GAc2mNl0TFwt40wq\noR8DzgtdnxvcS7MTZjYGtT57zuwkOTBmto5aMr/f3fcFt1MXZ527/xBYAC4nXXFuBbab2QvAl4Ff\nMrP7gZdTFGODu//f4J/fo1Zq+wDp+nl+Fzjq7n8bXP9Pagk+TTGG/QrwHXf/fnCdtji3AS+4+6vu\n/q/U6vw/TxdxJpXQnwB+ysy2mNlPAL9JrS6UJsuGpajFNxN8/nFgX/QNA3A3cNjd7wjdS1WcZnZ2\nffXdzNYDlwHPkqI43f0Wdz/P3d9J7b/F/e5+FfAoKYmxzszOCv5WhpltoFb7fZp0/TxPAEfN7KeD\nW78M/D0pijHiSmp/kNelLc6XgA+a2RvNzKj9PA/TTZwJFvovB54H/jdw0yAXHWJi+xJwHPhx8MP8\nLWoLEo8HMT8GvGXAMW4F/hU4BBwEngx+pm9NWZzvCWI7BDwF3BrcT1WcoXg/xJlF0dTFSK0+Xf89\nf7r+/07aYqXW2fJEEOufAxvTFmMQ51nA94A3he6lMc5Zag9CTwH3Ujt1s+M4NVgkIpITWhQVEckJ\nJXQRkZxQQhcRyQkldBGRnFBCFxHJCSV0EZGcUEIXEckJJXQRkZz4N9Wx+Io9YkQ9AAAAAElFTkSu\nQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fbfe111b6d0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.unique(Z_radolan_ipoli_max), 'ok')\n",
    "plt.plot(np.unique(Z_boxpol_ipoli_max),'or')\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGklJREFUeJzt3X+QVXd5x/H3s651gS2wNQ0UQmC9HWsS8gMdUztYvBRC\nY+uEpp2JjWxkJQb/qBij6YSm2ewSHGeiTJg06h84wYASG7UlJB1/hCS9s6UzaKKkRiHq4G6IEDaJ\nhp8JjsjTP+7dm7tnz9099/c9535eMwx7D/fc+9zN8vDN832+36+5OyIiEn9tjQ5ARESqQwldRCQh\nlNBFRBJCCV1EJCGU0EVEEkIJXUQkISZN6Gb2FjP7vpntM7Nnzaw/d73LzB4zs5+Z2ffMbEbtwxUR\nkWIsSh+6mU1199fM7E3A/wKfAP4B+LW7f87MbgO63H19bcMVEZFiIpVc3P213JdvAdoBB1YC23LX\ntwF/V/XoREQkskgJ3czazGwfcBTY7e5PAbPcfQTA3Y8C59cuTBERmUzUEfo5d18EXABcaWaXkB2l\nj3latYMTEZHo2kt5srufMLMMcDUwYmaz3H3EzGYDL4XdY2ZK9CIiZXB3K+X5UbpczhvtYDGzKcBV\nwAHgEaA397TVwK4Jgortr/7+/obH0Krxxzl2xd/4X3GPvxxRRuh/Amwzszay/wA85O7fNrO9wDfM\nbA3wPHBdWRGIiEhVTJrQ3f1Z4J0h138DLK9FUCIiUjqtFJ1EOp1udAgViXP8cY4dFH+jxT3+ckRa\nWFTRG5h5rd9DRCRpzAyv9qSoiIjEgxK6iEhCKKGLiCSEErqISEIooYuIJIQSuohIQiihi4gkhBK6\niEhCKKGLiCSEErqISEIooYuIJIQSuohIQiihi4gkhBK6iEhCKKGLiCSEErqISEIooYuIJIQSuohI\nQiihi4gkhBK6iEhCKKGLiCREe6MDEBGpt+eHhnigr49zhw/TNncuvRs3Mr+7u9FhVczcvbZvYOa1\nfg8RkYkUJvATM2Zwbt8+PnvoENOA00B/KsW63bubKqmbGe5uJd2jhC4iSfb80BD3XXUVGw4eZBrQ\nB6wHphU85zSwadUq+r/2tYbEGKachK4auogk2gN9fflkDtmkNy3wnGnAuSNH6htYDaiGLiKxFVYL\nB8Zce60gmUM2oZ9m/Ai9bc6c+gVeI5OWXMzsAmA7MAs4B2xx9/vMrB+4CXgp99Tb3f27Ifer5CIi\nFQsm7+Vr17JzzZr86Ps0cMu8eUwxG1Mfv76zk6+fOpVP4M8D9wIbofVq6GY2G5jt7s+YWSfwQ2Al\n8EHgpLvfM8n9SugikhdlVB28FjaRGUzUEF4fPwB8vrOT+3LPHU38MxYtYurJk7TNmdOUXS7lJPRJ\nSy7ufhQ4mvv6lJkdAOaOvmfJUYpIywpOUJ4GbhkcHDeqDl4LJuppwGWBZA7h9fGLgOkLF7IpleLc\nkSO0zZnDvzZhAq+GkmroZrYAuAL4PvBe4ONmdgPwNPBpdz9e7QBFJB6ijLxPnTw5ZoJyGjDrhRfG\nJevgtbBE/WbG18LPhVw7DXSlUk3VwVIrkRN6rtzyLeDm3Ej9S8Bd7u5m9hngHuDGGsUpIk0s6sh7\nbUdHpFF18FrYROZ1wLpAKWVk3jxuD7xnfyrFutw/LkkXKaGbWTvZZP5Vd98F4O4vFzzly8Cjxe4f\nGBjIf51Op0mn02WEKiLNKtgaWGzk/bYzZyKNqoPXesmWXQonMu9Ppbhx61Y2bdkyppQCsKmvL39t\nXUzKK5lMhkwmU9FrRFpYZGbbgVfc/VMF12bn6uuY2S3Au939QyH3alJUJOYmK6cM7d/P9pdeGnNP\nP7Ah+DrAhilTuO/11yfsTCl2rdknMqupJpOiZrYYWAU8a2b7AAduBz5kZleQ/cd0GPhYyRGLSNOJ\n1B4YMmkZZeR9HjBjxQo2dXZOOKoudi3JCbwatPRfRPLCauFR2gPDervDRtnN2O/drGoyQheR1hFW\nC4/SHjgfuBm4YdYsLr344tjXs+NKCV2khQXLK8Fl8hC9PfA84PLly8e1B7ZCu2CzUEIXaSGTbSN7\nfWen2gNjTDV0kYSabHIz6jL5/lSKa7du5fGC9sAxXS4F11ROqR7thy7SAqKsyAzrTAlOboa1FQJ8\n8j3voatgmbwSdWNoUlQk4aKuyLx+164xyTtscrPYNrKtskw+iZTQRZpIWJnk8S1bSt4LJawzJTi5\n2cv41Zeqg8ebErpIkwiOvg8Adz/0EF88e7bkvVDCOlOCk5vnAafmzWOgYPWl2grjTQldpEGCo/Hg\n6PsbkE/mUNpeKGGdKcX2PlECTw4ldJEGCKuFB0ff5xg/8v4osC6wF0pYC2FY8h4dfS9esqQeH1Ea\nQAldpMrK3Rc8OPoOm7SMuheKkndrUtuiSAXKPedybUcHO86cGftajN2J8ABwd3v7mBq69kJpHWpb\nFKmxSVdahrQLRt0XPGz0fdPataFlE5EwSugiEQXr3pWccxlWC+9Ppfjk5s3jErbKJhKVErpIRMGd\nCCs55zJsNK7Rt1RKCV0konOHD1f1nMuw0bhIJZTQRYoYt7Xs9Oktd86lxIu6XERChO6ZonMupY60\n26JIBFH7xAceeWTcxlV3rFzJzIK6txK41IraFkUmEXW3wrA9U6YB00+coP/hh+sdtkgkSuiSKJON\nvn8yPMz24eGy+sRPA21z5tTpk4iUTgldEiPK6PsOxrcaltInrq1lpZkpoUtihJ1YHxx9q09ckkwJ\nXWIryon1wdF3L28cvaY+cUkaJXSJhSibYIWdWB8cfc8HbgQ+vGABC7u71ScuiaK2RWl6YbXx4IHH\nEH5ifVjvuHYslDhQ26IkwmQn+RTbBOsiYPrChWwqOLFeo29pJUro0lSinOQD4ZObE51Yr1PspRUo\noUvDFY7Iw/rEw3rCwzbBUluhtLpJE7qZXQBsB2aRnWP6srv/m5l1AQ+RnWcaBq5z9+M1jFUSKDgi\nD+sTD+sJn+jMTJFWNemkqJnNBma7+zNm1gn8EFgJfAT4tbt/zsxuA7rcfX3I/ZoUlbyw+njhnikb\ngFsZX0rRHirSauqyOZeZPQx8Iffrfe4+kkv6GXd/R8jzldAFKF4fLzxb83ngPsb2iasrRVpRzbtc\nzGwBcAWwF5jl7iMA7n7UzM4v5bWk9YSt5AzWx8P6xFVKEYkmckLPlVu+Bdzs7qfMLDjsLjoMHxgY\nyH+dTqdJp9OlRSmJEDzxB4rXx+/RiFxaTCaTIZPJVPQakUouZtYO/BfwHXe/N3ftAJAuKLn8t7tf\nFHKvSi4CwIaeHm7dsUP1cZEIalZDN7PtwCvu/qmCa3cDv3H3uzUpKlGE1dBVHxcJV5OEbmaLgUHg\nWbJlFQduB34AfAOYR3Yu6zp3PxZyvxK65OW7XDQaF5mQjqCTphJ22ISSt0g0SujSNFReEamMEro0\nzGQLhiCb1DetWqV9VUQi0G6L0hBRN9SaBpw7cqT+AYq0iLZGByDxN9GCoUI6ZFmktjRCl7IUlliG\n9u/XIcsiTUAJXSY12fFvfeiQZZFmoElRGSPq2Z2Fx789D9wLbEQbaolUi7pcpCJRz+68A/hM8F7g\nllmzuPTii7VgSKQK1OUiFQmb3Aw7uzPs+LfzgMuXL1dLokgDqctF8sJ2QxxN3oVGj38bvT5aYunV\nhKdIQ2mELnltc+dGOrtTx7+JNCfV0CWv2HL9a7du5fGC5K36uEjtaVJUKqbdEEWagxK6lES7IYo0\nLyV0iUy7IYo0NyV0KUq7IYrEi/rQJa8wgZ+YMYNz+/bx2UOHtBuiSIIpoSdAlL1W1kPobojBEbp2\nQxSJLyX0mAtdrr9r15jl+m2g3RBFWoASesyE1cInW67fhnZDFGkFSugxEvVkoOBeK71kyy7B3RA/\nuXmzErhIgiihN7Eoo/GwWnhwuf55wKl58xhYtIipJ09qNC6SUGpbbBJR9iFf29HBjjNnxt4HbAip\nhWu5vki8qQ89pqLuQx7sViH33DtWrmRmQS1cyVsk/tSHHiOFI/KfDA+zfXh40n3Ii3WmqBYuIqCE\n3hDBEfkdjG8rLHaIhDpTRKQYJfQGCJ4MFJa8w/Yh12hcRCaihN4AwZOBeoF+YANvtBXqEAkRKdWk\nCd3M7gc+AIy4+2W5a/3ATcBLuafd7u7frVmUCRM8GWg+cCPw4QULWNjdPSZ5L16ypHGBikisRBmh\nfwW4D9geuH6Pu99T/ZCSJ6wlsX/v3jFdLfenUtyjrWtFpAKTJnR332Nm80P+qKR2mlYVuu/43r1c\nq3KKiFRZJTX0j5vZDcDTwKfd/XiVYoq1KKs7Nxw8yKYtW7TvuIhUVbkJ/UvAXe7uZvYZ4B6yZeBQ\nAwMD+a/T6TTpdLrMt21uUfda0b7jIhKUyWTIZDIVvUaklaK5ksujo5OiUf8s9+cts1J0Q08Pt+7Y\nEWl1p04GEpGJ1HKlqFFQMzez2e5+NPfw74GflPKmSREsr7xWUFoZpX3HRaReorQtPgikgbea2SGy\nLdNLzewK4BwwDHyshjE2pWL7r2h1p4g0ijbnKlNYeeUA8PmQ1Z3r1I4oIiXS5lw1VlhiGdq/f1x5\n5SJg+sKFbEqlNBoXkbpTQo8oWGLpY/z+K6eBrlRKk50i0hBtjQ4gLoIban2UN5I6vFFe6dVkp4g0\niEboEQU31JoP3AzcMGsWl158scorItJwSugRBTfUgmwHy+XLl6vEIiJNQSWXiHo3bqQ/lYptiWVo\naIienh6WLl1KT08PQ0NDjQ5JRKpMbYslyHe5xOzszqGhIa666ioOHjyYv5ZKpdi9ezfdMYhfpBXp\nkGjJGxoaoq+vj8OHDzM8PMzw8PC456xatYqvqVwk0pTUhy5A+Ig8zBFtECaSKEroRQT3aWnm8krh\naHzu3LmcPHly0mQOMGfOnDpEJyL1opIL4ScK7VyzZuyhFE2yhD+YvNeuXcuaNWvGJPCOjg7OnDkz\n4euohi7S3FRDL0OxTba+ntuPZVQzbHkbVkrp7Ozk1KlTke5fsGAB3d3dzJkzh40bNyqZizQx1dDL\nEFwBOg24LJDMR6/X+1CKKKWUYsl8ypQpvP766/nHGpGLJF/LJ/TgClCANxO+T0tbHWvOYaPxjo6O\nyPevWLGCzs5Ojhw5ohG5SIto+YQetgL0OmBd2Da4NVxEFGU0XqwuHiy7pFIpNm/erAQu0mJUQw+p\nofenUly7dSuPb9lSl0VExUbjYQk8rJSydetWtmzZotG4SIJoUrRMjVgBGmXhT5iVK1eqlCLSApTQ\nYyLqwh9NbIq0LnW5xERfX1+khT+a2BSRUiih10FwwjNKMtfEpoiUSgm9xootBgqjhT8iUomWS+i1\n3qMl6mKgsFZD1cdFpBItldBDWxT37q3aHi2lLAZauHAhqVRK9XERqZqWSuhhy/w3HDzIpr6+quzR\nEjbZWWwxUCqV0l7kIlJViU7owfLKawXJfFQle7REnewMaz/cGJOj60QkPhKb0IvtolitPVpKmexU\n+6GI1ENiFxZt6Onh1h07xiTvA8Dnw/ZoiVhDj7K6U5OdIlINWlhUIGwXxYuA6QsXsimVyi/zXxex\nyyXq6k5NdopIo0ya0M3sfuADwIi7X5a71gU8BMwHhoHr3P14DeMsWdguiqeBrlSqrAnQqKs7Ndkp\nIo0yacnFzN4LnAK2FyT0u4Ffu/vnzOw2oMvd1xe5v24ll8JJ0BMzZnBu3z4+e+hQxeWV0QnPvXv3\nTniPyisiUi0125zLzOYDjxYk9OeA97n7iJnNBjLu/o4i99YloYdNgt4ybx4zFi1i6smTJe2iWMpR\nb1rdKSK1UM8a+vnuPgLg7kfN7PwyX6dqwnrMN7/wApuWLKF/166SXiusvKLVnSLS7Ko1KTrhEHxg\nYCD/dTqdJp1OV+lt3xA2CVpuj/nhw4dDr2vCU0RqJZPJkMlkKnqNchP6iJnNKii5vDTRkwsTeq0U\nmwSN0mMerJdPnz499Hma8BSRWgkOdjds2FDya0StoS8gW0O/NPf4buA37n53s0yKFjtKbrJJ0LB6\n+bx58zAzDh06lL+m8oqI1FNNJkXN7EEgDbwVGAH6gYeBbwLzgOfJti0eK3J//btcSjhKrqenhx07\ndoy7rqPeRKSRdARdGZYuXRpat1q6dClPPvlk/QMSEUErRSOJWi+fU8b+LiIijdRSI3TVy0UkLlRy\nmYTq5SISFyq5TKJYf/mJEyd4+OGH6xyNiEh1xTahRzkbVPVyEWklsSy5ROk5V71cROKsZWroYYdX\nnAY2rVqV3xpX9XIRibOWqaFH2bdF9XIRaTVtjQ6gHKP7thQK7tsyd+7c0HtVLxeRpIplQu/duJH+\nVCqf1Edr6MvWrqWnp4elS5dy6tQpLrzwwjH3pVIpNm7cWPd4RUTqIZY1dBi/b8uytWvpXbNm3CTo\nokWLOHnypOrlIhIrLTMpGqbYJOiqVau05a2IxE45CT2WJZcwxSZBj5RxwIWISBwlJqFrElREWl1i\nSi5hC4m0aEhE4qqlaujBZf2j3St9fX1aNCQisdcyCV2jcRFJupaZFO3r6xuTzAEOHjxIX19fgyIS\nEWm8WCZ0dbSIiIwXy4SujhYRkfFUQxcRaUItMykKb3S5qKNFRJIosQl9z+Agd61eTcerr3Kmq4s7\nt23jvUuWVClCEZHmk8iEvmdwkHuXLeOBs2fzpxP1trdz8xNPKKmLSGIlMqGv6O5m5/DwuNOJrl2w\ngMeGhiqOT0SkGSWyD73j1VdDTyfqOHasEeGIiDStpk/oZ7q6Qk8nOjNzZiPCERFpWk2f0O/cto3e\n9vYxpxP1trdz57ZtjQxLRKTpVFRDN7Nh4DhwDvidu18Z8pzqdbkcO8aZmTPV5SIiiVf3SVEz+yXw\nLnd/dYLn1GX7XBGRJGnEpKhV4TVERKQKKk3GDuw2s6fM7KZqBLRncJAV3d1cM3MmK7q72TM4WI2X\nFRFJvPYK71/s7i+a2R+TTewH3H1PuS82uoho5+giouPH6V22DLSISERkUhUldHd/Mff7y2a2E7gS\nGJfQBwYG8l+n02nS6XTo6921enU+mUO23/yBs2e5dvVqLSISkUTLZDJkMpmKXqPsSVEzmwq0ufsp\nM5sGPAZscPfHAs+LPCl6zcyZPHL8ePj1V4vOu4qIJE45k6KVjNBnATvNzHOvsyOYzEt1pquL08eP\nj1vmr0VEIiKTa6q9XLQRl4hIViI259IiIhGRhCR0ERFJ6G6LIiISjRK6iEhCKKGLiCSEErqISEIo\noYuIJETDE7o24xIRqY6Gti1qIZGISLjY9aGv6O5m5/DwuKX+1y5YoM24RKSlxa4PvePVV8ckc8ju\nsNhx7FgjwhERibWGJvQzXV35w59HaTMuEZHyNDSh37ltG73t7fmkPlpDv3PbtkaGJSISSw3fy0Wb\ncYmIjBe7SVEREQkXu0lRERGpHiV0EZGEUEIXEUkIJXQRkYRQQhcRSQgldBGRhFBCFxFJCCV0EZGE\nUEIXEUkIJXQRkYRQQhcRSQgldBGRhFBCFxFJiIoSupldbWbPmdnPzey2agUlIiKlKzuhm1kb8AXg\nr4FLgOvN7B3VCqxZZDKZRodQkTjHH+fYQfE3WtzjL0clI/QrgV+4+/Pu/jvg34GVE92wZ3CQFd3d\nXDNzJiu6u9kzOFjB29dH3H8o4hx/nGMHxd9ocY+/HO0V3DsXeKHg8a/IJvlQewYHuXfZMnaePcs0\n4PTx4/QuWwZPPKETikREqqBuk6J3rV7NA7lkDjANeODsWe5avbpeIYiIJFrZR9CZ2XuAAXe/Ovd4\nPeDufnfgeTp/TkSkDHU7U9TM3gT8DFgGvAj8ALje3Q+U9YIiIlKRsmvo7v57M/s48BjZ0s39SuYi\nIo1T9ghdRESaS80mReO46MjM7jezETP7ccG1LjN7zMx+ZmbfM7MZjYyxGDO7wMyeNLOfmtmzZvaJ\n3PW4xP8WM/u+me3Lxd+fux6L+CG7NsPMfmRmj+QexyZ2ADMbNrP/y/03+EHuWiw+g5nNMLNvmtmB\n3N+BP49R7G/Pfc9/lPv9uJl9opz4a5LQY7zo6CtkYy60Hnjc3f8MeBL4l7pHFc1Z4FPufgnwF8A/\n5b7nsYjf3X8LLHX3RcAVwPvN7EpiEn/OzcD+gsdxih3gHJB290XuPtqCHJfPcC/wbXe/CLgceI6Y\nxO7uP899z98JvAs4DeyknPjdveq/gPcA3yl4vB64rRbvVYPY5wM/Lnj8HDAr9/Vs4LlGxxjxczwM\nLI9j/MBU4Gng3XGJH7gA2A2kgUfi+LMDDAFvDVxr+s8ATAcOhlxv+thDYl4B/E+58deq5BK26Ghu\njd6r1s539xEAdz8KnN/geCZlZgvIjnL3kv2BiEX8uZLFPuAosNvdnyI+8W8G/hkonJSKS+yjHNht\nZk+Z2Udz1+LwGbqBV8zsK7myxRYzm0o8Yg/6IPBg7uuS49dui6Vr6llkM+sEvgXc7O6nGB9v08bv\n7uc8W3K5ALjSzC4hBvGb2d8CI+7+DDBR33DTxR6w2LP/2/83ZEt2f0kMvv9ku/XeCXwxF/9pslWB\nOMSeZ2ZvBq4Bvpm7VHL8tUroh4ELCx5fkLsWRyNmNgvAzGYDLzU4nqLMrJ1sMv+qu+/KXY5N/KPc\n/QSQAa4mHvEvBq4xs18CXwf+ysy+ChyNQex57v5i7veXyZbsriQe3/9fAS+4+9O5x/9BNsHHIfZC\n7wd+6O6v5B6XHH+tEvpTwJ+a2Xwz+wPgH4FHavRe1WaMHWU9AvTmvl4N7Are0ES2Avvd/d6Ca7GI\n38zOG53FN7MpwFXAAWIQv7vf7u4XuvvbyP6sP+nuNwCP0uSxjzKzqbn/u8PMppGt5T5LPL7/I8AL\nZvb23KVlwE+JQewB15MdEIwqPf4aFvevJruS9BfA+kZPNkSM+UHgCPBb4BDwEaALeDz3WR4DZjY6\nziKxLwZ+DzwD7AN+lPtv8Ecxif/SXMzPAD8G/jV3PRbxF3yO9/HGpGhsYidbhx792Xl29O9sXD4D\n2c6Wp3Kf4T+BGXGJPRf/VOBl4A8LrpUcvxYWiYgkhCZFRUQSQgldRCQhlNBFRBJCCV1EJCGU0EVE\nEkIJXUQkIZTQRUQSQgldRCQh/h8IU3cYV1N3dQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fbfdf6fa710>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(np.unique(Z_radolan_ipoli_min), 'ok')\n",
    "plt.plot(np.unique(Z_boxpol_ipoli_min),'or')\n",
    "#plt.scatter(Z_radolan_ipoli_min,Z_boxpol_ipoli_min)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "ename": "TypeError",
     "evalue": "'numpy.ndarray' object is not callable",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mTypeError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-7-b9ee479a4eea>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mZ_radolan_ipoli_median\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\",\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"Z_radolan_ipoli_max,\"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m\"Z_radolan_ipoli_min\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[1;31mTypeError\u001b[0m: 'numpy.ndarray' object is not callable"
     ]
    }
   ],
   "source": [
    ", Z_radolan_ipoli_median, Z_radolan_ipoli_max, Z_radolan_ipoli_min"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "a = np.array([np.nan,np.nan,1,2,3,4.4,5.6,7.8])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "len(a\n",
    "   )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "np.count_nonzero(~np.isnan(a))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2.718281828459045"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.exp(1./1.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}