Xarray
\n", + "\n", + "\n", + "> *DS Python for GIS and Geoscience* \n", + "> *October, 2020*\n", + ">\n", + "> *© 2020, Joris Van den Bossche and Stijn Van Hoey. Licensed under [CC BY 4.0 Creative Commons](http://creativecommons.org/licenses/by/4.0/)*\n", + "\n", + "---" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt" + ] + }, + { + "cell_type": "markdown", + "metadata": { + "deletable": true, + "editable": true + }, + "source": [ + "## Introduction: meet `xarray`" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Let's start with reading the Sentinel RGB bands for Herstappe again:" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": { + "deletable": true, + "editable": true + }, + "outputs": [], + "source": [ + "import xarray as xr" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": { + "deletable": true, + "editable": true + }, + "outputs": [], + "source": [ + "data_file = \"./data/herstappe/raster/2020-09-17_Sentinel_2_L1C_True_color.tiff\"" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": { + "collapsed": false, + "deletable": true, + "editable": true, + "jupyter": { + "outputs_hidden": false + } + }, + "outputs": [ + { + "data": { + "text/html": [ + "<xarray.DataArray (band: 3, y: 227, x: 447)>\n", + "[304407 values with dtype=float32]\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[304407 values with dtype=float32]
array([1, 2, 3])
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (band: 3, y: 227, x: 447)>\n", + "[304407 values with dtype=float32]\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[304407 values with dtype=float32]
array([1, 2, 3])
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (band: 3, y: 227, x: 447)>\n", + "array([[[0.3375 , 0.2195 , ..., 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, ..., 0.21275, 0.27125],\n", + " ...,\n", + " [0.38625, 0.40475, ..., 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, ..., 0.388 , 0.38725]],\n", + "\n", + " [[0.311 , 0.2525 , ..., 0.2315 , 0.24775],\n", + " [0.3485 , 0.312 , ..., 0.25475, 0.29 ],\n", + " ...,\n", + " [0.34225, 0.3465 , ..., 0.338 , 0.31475],\n", + " [0.3405 , 0.35075, ..., 0.341 , 0.349 ]],\n", + "\n", + " [[0.3235 , 0.277 , ..., 0.266 , 0.27925],\n", + " [0.35775, 0.31775, ..., 0.28525, 0.31625],\n", + " ...,\n", + " [0.3485 , 0.346 , ..., 0.334 , 0.31725],\n", + " [0.338 , 0.34375, ..., 0.3375 , 0.33775]]], dtype=float32)\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[[0.3375 , 0.2195 , ..., 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, ..., 0.21275, 0.27125],\n", + " ...,\n", + " [0.38625, 0.40475, ..., 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, ..., 0.388 , 0.38725]],\n", + "\n", + " [[0.311 , 0.2525 , ..., 0.2315 , 0.24775],\n", + " [0.3485 , 0.312 , ..., 0.25475, 0.29 ],\n", + " ...,\n", + " [0.34225, 0.3465 , ..., 0.338 , 0.31475],\n", + " [0.3405 , 0.35075, ..., 0.341 , 0.349 ]],\n", + "\n", + " [[0.3235 , 0.277 , ..., 0.266 , 0.27925],\n", + " [0.35775, 0.31775, ..., 0.28525, 0.31625],\n", + " ...,\n", + " [0.3485 , 0.346 , ..., 0.334 , 0.31725],\n", + " [0.338 , 0.34375, ..., 0.3375 , 0.33775]]], dtype=float32)
array([1, 2, 3])
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (band: 3, y: 227, x: 447)>\n", + "array([[[0.3375 , 0.2195 , ..., 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, ..., 0.21275, 0.27125],\n", + " ...,\n", + " [0.38625, 0.40475, ..., 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, ..., 0.388 , 0.38725]],\n", + "\n", + " [[0.311 , 0.2525 , ..., 0.2315 , 0.24775],\n", + " [0.3485 , 0.312 , ..., 0.25475, 0.29 ],\n", + " ...,\n", + " [0.34225, 0.3465 , ..., 0.338 , 0.31475],\n", + " [0.3405 , 0.35075, ..., 0.341 , 0.349 ]],\n", + "\n", + " [[0.3235 , 0.277 , ..., 0.266 , 0.27925],\n", + " [0.35775, 0.31775, ..., 0.28525, 0.31625],\n", + " ...,\n", + " [0.3485 , 0.346 , ..., 0.334 , 0.31725],\n", + " [0.338 , 0.34375, ..., 0.3375 , 0.33775]]], dtype=float32)\n", + "Coordinates:\n", + " * band (band) <U1 'R' 'G' 'B'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[[0.3375 , 0.2195 , ..., 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, ..., 0.21275, 0.27125],\n", + " ...,\n", + " [0.38625, 0.40475, ..., 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, ..., 0.388 , 0.38725]],\n", + "\n", + " [[0.311 , 0.2525 , ..., 0.2315 , 0.24775],\n", + " [0.3485 , 0.312 , ..., 0.25475, 0.29 ],\n", + " ...,\n", + " [0.34225, 0.3465 , ..., 0.338 , 0.31475],\n", + " [0.3405 , 0.35075, ..., 0.341 , 0.349 ]],\n", + "\n", + " [[0.3235 , 0.277 , ..., 0.266 , 0.27925],\n", + " [0.35775, 0.31775, ..., 0.28525, 0.31625],\n", + " ...,\n", + " [0.3485 , 0.346 , ..., 0.334 , 0.31725],\n", + " [0.338 , 0.34375, ..., 0.3375 , 0.33775]]], dtype=float32)
array(['R', 'G', 'B'], dtype='<U1')
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "array([[0.3375 , 0.2195 , 0.17675, ..., 0.167 , 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, 0.222 , ..., 0.197 , 0.21275, 0.27125],\n", + " [0.30775, 0.26425, 0.25875, ..., 0.272 , 0.28975, 0.27425],\n", + " ...,\n", + " [0.39175, 0.38825, 0.4335 , ..., 0.3815 , 0.3135 , 0.2055 ],\n", + " [0.38625, 0.40475, 0.42125, ..., 0.39725, 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, 0.41775, ..., 0.386 , 0.388 , 0.38725]],\n", + " dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[0.3375 , 0.2195 , 0.17675, ..., 0.167 , 0.17625, 0.2485 ],\n", + " [0.362 , 0.29525, 0.222 , ..., 0.197 , 0.21275, 0.27125],\n", + " [0.30775, 0.26425, 0.25875, ..., 0.272 , 0.28975, 0.27425],\n", + " ...,\n", + " [0.39175, 0.38825, 0.4335 , ..., 0.3815 , 0.3135 , 0.2055 ],\n", + " [0.38625, 0.40475, 0.42125, ..., 0.39725, 0.3875 , 0.3235 ],\n", + " [0.3905 , 0.40375, 0.41775, ..., 0.386 , 0.388 , 0.38725]],\n", + " dtype=float32)
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "[101469 values with dtype=float32]\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[101469 values with dtype=float32]
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "[101469 values with dtype=float32]\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[101469 values with dtype=float32]
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 10, x: 10)>\n", + "array([[0.35 , 0.1795 , 0.18425, 0.28175, 0.168 , 0.30175, 0.40075,\n", + " 0.3285 , 0.24075, 0.32025],\n", + " [0.361 , 0.1615 , 0.1495 , 0.35125, 0.2855 , 0.182 , 0.16325,\n", + " 0.19125, 0.16225, 0.17125],\n", + " [0.4145 , 0.40825, 0.1885 , 0.1665 , 0.17575, 0.21025, 0.2005 ,\n", + " 0.37025, 0.16375, 0.41225],\n", + " [0.22225, 0.41 , 0.318 , 0.1625 , 0.336 , 0.33325, 0.1925 ,\n", + " 0.2195 , 0.17925, 0.24175],\n", + " [0.1955 , 0.263 , 0.216 , 0.3695 , 0.21975, 0.26475, 0.178 ,\n", + " 0.2425 , 0.27425, 0.37175],\n", + " [0.1805 , 0.40075, 0.17625, 0.16225, 0.414 , 0.1695 , 0.15875,\n", + " 0.19075, 0.20575, 0.24225],\n", + " [0.205 , 0.1875 , 0.1635 , 0.3485 , 0.21375, 0.23775, 0.32575,\n", + " 0.17825, 0.332 , 0.18325],\n", + " [0.21775, 0.16275, 0.177 , 0.16475, 0.44825, 0.44575, 0.221 ,\n", + " 0.1505 , 0.3005 , 0.18625],\n", + " [0.2065 , 0.22525, 0.15925, 0.16425, 0.16 , 0.31875, 0.30975,\n", + " 0.26425, 0.15025, 0.289 ],\n", + " [0.203 , 0.1985 , 0.44625, 0.50825, 0.36 , 0.3295 , 0.21725,\n", + " 0.156 , 0.344 , 0.18875]], dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.574e+06 6.573e+06 6.573e+06 ... 6.57e+06 6.57e+06\n", + " * x (x) float64 5.975e+05 5.979e+05 5.983e+05 ... 6.007e+05 6.011e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[0.35 , 0.1795 , 0.18425, 0.28175, 0.168 , 0.30175, 0.40075,\n", + " 0.3285 , 0.24075, 0.32025],\n", + " [0.361 , 0.1615 , 0.1495 , 0.35125, 0.2855 , 0.182 , 0.16325,\n", + " 0.19125, 0.16225, 0.17125],\n", + " [0.4145 , 0.40825, 0.1885 , 0.1665 , 0.17575, 0.21025, 0.2005 ,\n", + " 0.37025, 0.16375, 0.41225],\n", + " [0.22225, 0.41 , 0.318 , 0.1625 , 0.336 , 0.33325, 0.1925 ,\n", + " 0.2195 , 0.17925, 0.24175],\n", + " [0.1955 , 0.263 , 0.216 , 0.3695 , 0.21975, 0.26475, 0.178 ,\n", + " 0.2425 , 0.27425, 0.37175],\n", + " [0.1805 , 0.40075, 0.17625, 0.16225, 0.414 , 0.1695 , 0.15875,\n", + " 0.19075, 0.20575, 0.24225],\n", + " [0.205 , 0.1875 , 0.1635 , 0.3485 , 0.21375, 0.23775, 0.32575,\n", + " 0.17825, 0.332 , 0.18325],\n", + " [0.21775, 0.16275, 0.177 , 0.16475, 0.44825, 0.44575, 0.221 ,\n", + " 0.1505 , 0.3005 , 0.18625],\n", + " [0.2065 , 0.22525, 0.15925, 0.16425, 0.16 , 0.31875, 0.30975,\n", + " 0.26425, 0.15025, 0.289 ],\n", + " [0.203 , 0.1985 , 0.44625, 0.50825, 0.36 , 0.3295 , 0.21725,\n", + " 0.156 , 0.344 , 0.18875]], dtype=float32)
array('R', dtype='<U1')array([6573518.286716, 6573117.160681, 6572716.034646, 6572314.908611,\n", + " 6571913.782576, 6571512.656541, 6571111.530506, 6570710.404471,\n", + " 6570309.278436, 6569908.152401])
array([597504.049916, 597904.831314, 598305.612713, 598706.394111,\n", + " 599107.175509, 599507.956908, 599908.738306, 600309.519704,\n", + " 600710.301103, 601111.082501])
<xarray.DataArray (band: 3, y: 227, x: 447)>\n", + "array([[[ True, True, False, ..., False, False, True],\n", + " [ True, True, True, ..., False, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]],\n", + "\n", + " [[ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]],\n", + "\n", + " [[ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]]])\n", + "Coordinates:\n", + " * band (band) <U1 'R' 'G' 'B'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05
array([[[ True, True, False, ..., False, False, True],\n", + " [ True, True, True, ..., False, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]],\n", + "\n", + " [[ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]],\n", + "\n", + " [[ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " ...,\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True],\n", + " [ True, True, True, ..., True, True, True]]])
array(['R', 'G', 'B'], dtype='<U1')
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 317, x: 625)>\n", + "array([[65535, 65535, 65535, ..., 65535, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 65535, 14221, 13779],\n", + " [65535, 65535, 65535, ..., 14631, 65535, 65535],\n", + " ...,\n", + " [65535, 65535, 65535, ..., 65535, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 14516, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 13844, 65535, 65535]], dtype=uint16)\n", + "Coordinates:\n", + " band int64 2\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[65535, 65535, 65535, ..., 65535, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 65535, 14221, 13779],\n", + " [65535, 65535, 65535, ..., 14631, 65535, 65535],\n", + " ...,\n", + " [65535, 65535, 65535, ..., 65535, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 14516, 65535, 65535],\n", + " [65535, 65535, 65535, ..., 13844, 65535, 65535]], dtype=uint16)
array(2)
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
<xarray.DataArray ()>\n", + "array(30)
array(30)
<xarray.DataArray ()>\n", + "array(30)
array(30)
<xarray.DataArray ()>\n", + "array(0.27586532, dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'
array(0.27586532, dtype=float32)
array('R', dtype='<U1')<xarray.DataArray (y: 227)>\n", + "array([0.27257997, 0.2730302 , 0.27882323, 0.28289318, 0.28358892,\n", + " 0.2827528 , 0.28218287, 0.27802128, 0.27297595, 0.27357492,\n", + " 0.2750811 , 0.27421978, 0.27253968, 0.267533 , 0.26884004,\n", + " 0.26667392, 0.26950392, 0.27182943, 0.27148938, 0.27283612,\n", + " 0.2725464 , 0.27282214, 0.2732914 , 0.2809407 , 0.28457773,\n", + " 0.28962082, 0.29109505, 0.2883283 , 0.28721252, 0.2862478 ,\n", + " 0.28497875, 0.27839935, 0.27497035, 0.27145863, 0.27309564,\n", + " 0.2773652 , 0.2777612 , 0.27870917, 0.27950057, 0.28156993,\n", + " 0.28009897, 0.2784961 , 0.28004864, 0.28016555, 0.2799614 ,\n", + " 0.28029698, 0.27859282, 0.27868232, 0.27850223, 0.28422594,\n", + " 0.2747545 , 0.27626175, 0.2782226 , 0.27482104, 0.2757511 ,\n", + " 0.27877182, 0.27805147, 0.27441725, 0.27426678, 0.2760285 ,\n", + " 0.2784206 , 0.27596644, 0.27611804, 0.2747377 , 0.2743188 ,\n", + " 0.2724491 , 0.2672086 , 0.26089823, 0.2578473 , 0.2546292 ,\n", + " 0.2492981 , 0.2499944 , 0.24945639, 0.2477114 , 0.25369686,\n", + " 0.2589916 , 0.2613322 , 0.26193288, 0.2671068 , 0.26986575,\n", + " 0.2739463 , 0.27676624, 0.27886856, 0.27532104, 0.275915 ,\n", + " 0.2746516 , 0.27479023, 0.2731057 , 0.27007213, 0.27255148,\n", + " 0.27285233, 0.26989484, 0.27013814, 0.2779396 , 0.27601507,\n", + " 0.27753636, 0.28057772, 0.2804737 , 0.28125224, 0.28297094,\n", + "...\n", + " 0.28391612, 0.28284732, 0.28177014, 0.28556767, 0.28089765,\n", + " 0.27344462, 0.27452073, 0.27802572, 0.2784575 , 0.28234118,\n", + " 0.2862209 , 0.290736 , 0.28936243, 0.2865095 , 0.28146195,\n", + " 0.28202295, 0.27974552, 0.27544856, 0.26893288, 0.26566273,\n", + " 0.2634905 , 0.2642114 , 0.27001736, 0.27205536, 0.2740498 ,\n", + " 0.28036857, 0.28384957, 0.2881885 , 0.2849536 , 0.28080592,\n", + " 0.27967507, 0.27667618, 0.2811029 , 0.27860457, 0.27888033,\n", + " 0.28203857, 0.28413534, 0.28942174, 0.2876454 , 0.28944966,\n", + " 0.2917763 , 0.28778747, 0.2874709 , 0.28477237, 0.27712303,\n", + " 0.27817672, 0.27468625, 0.2745794 , 0.27197203, 0.2717567 ,\n", + " 0.2735403 , 0.276094 , 0.274481 , 0.27807832, 0.2762254 ,\n", + " 0.2777869 , 0.27896476, 0.28867394, 0.29002458, 0.29157272,\n", + " 0.29093847, 0.2968149 , 0.29821423, 0.29791832, 0.29627797,\n", + " 0.2927662 , 0.28936523, 0.28149047, 0.2772953 , 0.27569237,\n", + " 0.27091107, 0.2602349 , 0.26218903, 0.25784284, 0.25383276,\n", + " 0.24703188, 0.24943289, 0.2481499 , 0.24555984, 0.25293064,\n", + " 0.25616834, 0.25984564, 0.2595867 , 0.2658294 , 0.25919014,\n", + " 0.25744295, 0.25550896, 0.2572332 , 0.25469968, 0.26596195,\n", + " 0.2730632 , 0.2658529 , 0.2639066 , 0.26685235, 0.27544463,\n", + " 0.28415158, 0.28765494], dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06
array([0.27257997, 0.2730302 , 0.27882323, 0.28289318, 0.28358892,\n", + " 0.2827528 , 0.28218287, 0.27802128, 0.27297595, 0.27357492,\n", + " 0.2750811 , 0.27421978, 0.27253968, 0.267533 , 0.26884004,\n", + " 0.26667392, 0.26950392, 0.27182943, 0.27148938, 0.27283612,\n", + " 0.2725464 , 0.27282214, 0.2732914 , 0.2809407 , 0.28457773,\n", + " 0.28962082, 0.29109505, 0.2883283 , 0.28721252, 0.2862478 ,\n", + " 0.28497875, 0.27839935, 0.27497035, 0.27145863, 0.27309564,\n", + " 0.2773652 , 0.2777612 , 0.27870917, 0.27950057, 0.28156993,\n", + " 0.28009897, 0.2784961 , 0.28004864, 0.28016555, 0.2799614 ,\n", + " 0.28029698, 0.27859282, 0.27868232, 0.27850223, 0.28422594,\n", + " 0.2747545 , 0.27626175, 0.2782226 , 0.27482104, 0.2757511 ,\n", + " 0.27877182, 0.27805147, 0.27441725, 0.27426678, 0.2760285 ,\n", + " 0.2784206 , 0.27596644, 0.27611804, 0.2747377 , 0.2743188 ,\n", + " 0.2724491 , 0.2672086 , 0.26089823, 0.2578473 , 0.2546292 ,\n", + " 0.2492981 , 0.2499944 , 0.24945639, 0.2477114 , 0.25369686,\n", + " 0.2589916 , 0.2613322 , 0.26193288, 0.2671068 , 0.26986575,\n", + " 0.2739463 , 0.27676624, 0.27886856, 0.27532104, 0.275915 ,\n", + " 0.2746516 , 0.27479023, 0.2731057 , 0.27007213, 0.27255148,\n", + " 0.27285233, 0.26989484, 0.27013814, 0.2779396 , 0.27601507,\n", + " 0.27753636, 0.28057772, 0.2804737 , 0.28125224, 0.28297094,\n", + "...\n", + " 0.28391612, 0.28284732, 0.28177014, 0.28556767, 0.28089765,\n", + " 0.27344462, 0.27452073, 0.27802572, 0.2784575 , 0.28234118,\n", + " 0.2862209 , 0.290736 , 0.28936243, 0.2865095 , 0.28146195,\n", + " 0.28202295, 0.27974552, 0.27544856, 0.26893288, 0.26566273,\n", + " 0.2634905 , 0.2642114 , 0.27001736, 0.27205536, 0.2740498 ,\n", + " 0.28036857, 0.28384957, 0.2881885 , 0.2849536 , 0.28080592,\n", + " 0.27967507, 0.27667618, 0.2811029 , 0.27860457, 0.27888033,\n", + " 0.28203857, 0.28413534, 0.28942174, 0.2876454 , 0.28944966,\n", + " 0.2917763 , 0.28778747, 0.2874709 , 0.28477237, 0.27712303,\n", + " 0.27817672, 0.27468625, 0.2745794 , 0.27197203, 0.2717567 ,\n", + " 0.2735403 , 0.276094 , 0.274481 , 0.27807832, 0.2762254 ,\n", + " 0.2777869 , 0.27896476, 0.28867394, 0.29002458, 0.29157272,\n", + " 0.29093847, 0.2968149 , 0.29821423, 0.29791832, 0.29627797,\n", + " 0.2927662 , 0.28936523, 0.28149047, 0.2772953 , 0.27569237,\n", + " 0.27091107, 0.2602349 , 0.26218903, 0.25784284, 0.25383276,\n", + " 0.24703188, 0.24943289, 0.2481499 , 0.24555984, 0.25293064,\n", + " 0.25616834, 0.25984564, 0.2595867 , 0.2658294 , 0.25919014,\n", + " 0.25744295, 0.25550896, 0.2572332 , 0.25469968, 0.26596195,\n", + " 0.2730632 , 0.2658529 , 0.2639066 , 0.26685235, 0.27544463,\n", + " 0.28415158, 0.28765494], dtype=float32)
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
<xarray.DataArray (band: 3)>\n", + "array([0.27586532, 0.28978536, 0.29871756], dtype=float32)\n", + "Coordinates:\n", + " * band (band) <U1 'R' 'G' 'B'
array([0.27586532, 0.28978536, 0.29871756], dtype=float32)
array(['R', 'G', 'B'], dtype='<U1')
<xarray.DataArray (quantile: 3, band: 3)>\n", + "array([[0.17125 , 0.23675001, 0.25999999],\n", + " [0.2615 , 0.28600001, 0.29475001],\n", + " [0.39274999, 0.34474999, 0.34024999]])\n", + "Coordinates:\n", + " * band (band) <U1 'R' 'G' 'B'\n", + " * quantile (quantile) float64 0.1 0.5 0.9
array([[0.17125 , 0.23675001, 0.25999999],\n", + " [0.2615 , 0.28600001, 0.29475001],\n", + " [0.39274999, 0.34474999, 0.34024999]])
array(['R', 'G', 'B'], dtype='<U1')
array([0.1, 0.5, 0.9])
<xarray.DataArray (y: 227, x: 447)>\n", + "array([[3.375 , 2.1950002, 1.7675 , ..., 1.67 , 1.7624999,\n", + " 2.4850001],\n", + " [3.62 , 2.9524999, 2.22 , ..., 1.97 , 2.1275 ,\n", + " 2.7125 ],\n", + " [3.0774999, 2.6425002, 2.5874999, ..., 2.7200003, 2.8975 ,\n", + " 2.7425 ],\n", + " ...,\n", + " [3.9175 , 3.8825 , 4.335 , ..., 3.815 , 3.1349998,\n", + " 2.055 ],\n", + " [3.8625 , 4.0474997, 4.2124996, ..., 3.9724998, 3.875 ,\n", + " 3.2350001],\n", + " [3.9050002, 4.0375 , 4.1775 , ..., 3.8600001, 3.88 ,\n", + " 3.8725 ]], dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05
array([[3.375 , 2.1950002, 1.7675 , ..., 1.67 , 1.7624999,\n", + " 2.4850001],\n", + " [3.62 , 2.9524999, 2.22 , ..., 1.97 , 2.1275 ,\n", + " 2.7125 ],\n", + " [3.0774999, 2.6425002, 2.5874999, ..., 2.7200003, 2.8975 ,\n", + " 2.7425 ],\n", + " ...,\n", + " [3.9175 , 3.8825 , 4.335 , ..., 3.815 , 3.1349998,\n", + " 2.055 ],\n", + " [3.8625 , 4.0474997, 4.2124996, ..., 3.9724998, 3.875 ,\n", + " 3.2350001],\n", + " [3.9050002, 4.0375 , 4.1775 , ..., 3.8600001, 3.88 ,\n", + " 3.8725 ]], dtype=float32)
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "array([[-1.0861897 , -1.5164031 , -1.733019 , ..., -1.7897615 ,\n", + " -1.7358519 , -1.3923124 ],\n", + " [-1.0161111 , -1.2199328 , -1.5050778 , ..., -1.6245515 ,\n", + " -1.5476375 , -1.3047143 ],\n", + " [-1.1784675 , -1.3308597 , -1.351893 , ..., -1.3019532 ,\n", + " -1.2387367 , -1.2937152 ],\n", + " ...,\n", + " [-0.9371314 , -0.94610584, -0.8358635 , ..., -0.96364444,\n", + " -1.159956 , -1.5823092 ],\n", + " [-0.95127046, -0.9044857 , -0.86452883, ..., -0.9231895 ,\n", + " -0.9480395 , -1.1285561 ],\n", + " [-0.9403273 , -0.9069594 , -0.8728721 , ..., -0.9519179 ,\n", + " -0.9467499 , -0.9486848 ]], dtype=float32)\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05\n", + "Attributes:\n", + " transform: (40.07813983073409, 0.0, 593476.1968631798, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (40.07813983073409, 40.1126034998577)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[-1.0861897 , -1.5164031 , -1.733019 , ..., -1.7897615 ,\n", + " -1.7358519 , -1.3923124 ],\n", + " [-1.0161111 , -1.2199328 , -1.5050778 , ..., -1.6245515 ,\n", + " -1.5476375 , -1.3047143 ],\n", + " [-1.1784675 , -1.3308597 , -1.351893 , ..., -1.3019532 ,\n", + " -1.2387367 , -1.2937152 ],\n", + " ...,\n", + " [-0.9371314 , -0.94610584, -0.8358635 , ..., -0.96364444,\n", + " -1.159956 , -1.5823092 ],\n", + " [-0.95127046, -0.9044857 , -0.86452883, ..., -0.9231895 ,\n", + " -0.9480395 , -1.1285561 ],\n", + " [-0.9403273 , -0.9069594 , -0.8728721 , ..., -0.9519179 ,\n", + " -0.9467499 , -0.9486848 ]], dtype=float32)
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "array([[ 0.02650002, -0.03299999, -0.06224999, ..., -0.05825001,\n", + " -0.05525 , 0.00075001],\n", + " [ 0.01349998, -0.01675001, -0.03525001, ..., -0.03925 ,\n", + " -0.04200001, -0.01874998],\n", + " [-0.00600001, -0.01349998, -0.01899999, ..., -0.01525 ,\n", + " -0.01499999, 0.0045 ],\n", + " ...,\n", + " [ 0.04975 , 0.0465 , 0.06774998, ..., 0.04425001,\n", + " 0.00424999, -0.04525 ],\n", + " [ 0.044 , 0.05824998, 0.06449997, ..., 0.05574998,\n", + " 0.04949999, 0.00875002],\n", + " [ 0.05000001, 0.053 , 0.05849999, ..., 0.04975 ,\n", + " 0.04700002, 0.03825 ]], dtype=float32)\n", + "Coordinates:\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05
array([[ 0.02650002, -0.03299999, -0.06224999, ..., -0.05825001,\n", + " -0.05525 , 0.00075001],\n", + " [ 0.01349998, -0.01675001, -0.03525001, ..., -0.03925 ,\n", + " -0.04200001, -0.01874998],\n", + " [-0.00600001, -0.01349998, -0.01899999, ..., -0.01525 ,\n", + " -0.01499999, 0.0045 ],\n", + " ...,\n", + " [ 0.04975 , 0.0465 , 0.06774998, ..., 0.04425001,\n", + " 0.00424999, -0.04525 ],\n", + " [ 0.044 , 0.05824998, 0.06449997, ..., 0.05574998,\n", + " 0.04949999, 0.00875002],\n", + " [ 0.05000001, 0.053 , 0.05849999, ..., 0.04975 ,\n", + " 0.04700002, 0.03825 ]], dtype=float32)
array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (y: 227, x: 447)>\n", + "array([[593496.5734331 , 593536.53357293, 593576.56896276, ...,\n", + " 611291.09701794, 611331.18440777, 611371.33479761],\n", + " [593496.59793308, 593536.60932292, 593576.61421276, ...,\n", + " 611291.12701794, 611331.22090777, 611371.35754761],\n", + " [593496.54368308, 593536.57832294, 593576.65096275, ...,\n", + " 611291.20201796, 611331.29790778, 611371.3605476 ],\n", + " ...,\n", + " [593496.6276831 , 593536.70232292, 593576.82571275, ...,\n", + " 611291.31151795, 611331.32165776, 611371.29179761],\n", + " [593496.62218308, 593536.71882292, 593576.81346274, ...,\n", + " 611291.32726794, 611331.39565776, 611371.40979761],\n", + " [593496.6264331 , 593536.71782293, 593576.80996276, ...,\n", + " 611291.31601795, 611331.39615778, 611371.47354761]])\n", + "Coordinates:\n", + " band <U1 'R'\n", + " * y (y) float64 6.578e+06 6.577e+06 6.577e+06 ... 6.569e+06 6.568e+06\n", + " * x (x) float64 5.935e+05 5.935e+05 5.936e+05 ... 6.113e+05 6.114e+05
array([[593496.5734331 , 593536.53357293, 593576.56896276, ...,\n", + " 611291.09701794, 611331.18440777, 611371.33479761],\n", + " [593496.59793308, 593536.60932292, 593576.61421276, ...,\n", + " 611291.12701794, 611331.22090777, 611371.35754761],\n", + " [593496.54368308, 593536.57832294, 593576.65096275, ...,\n", + " 611291.20201796, 611331.29790778, 611371.3605476 ],\n", + " ...,\n", + " [593496.6276831 , 593536.70232292, 593576.82571275, ...,\n", + " 611291.31151795, 611331.32165776, 611371.29179761],\n", + " [593496.62218308, 593536.71882292, 593576.81346274, ...,\n", + " 611291.32726794, 611331.39565776, 611371.40979761],\n", + " [593496.6264331 , 593536.71782293, 593576.80996276, ...,\n", + " 611291.31601795, 611331.39615778, 611371.47354761]])
array('R', dtype='<U1')array([6577529.547066, 6577489.434462, 6577449.321859, ..., 6568544.323882,\n", + " 6568504.211278, 6568464.098675])
array([593496.235933, 593536.314073, 593576.392213, ..., 611290.930018,\n", + " 611331.008158, 611371.086298])
<xarray.DataArray (band: 3, y: 317, x: 625)>\n", + "[594375 values with dtype=uint16]\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[594375 values with dtype=uint16]
array([1, 2, 3])
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
Xarray Datasets
\n", + "\n", + "\n", + "> *DS Python for GIS and Geoscience* \n", + "> *October, 2020*\n", + ">\n", + "> *© 2020, Joris Van den Bossche and Stijn Van Hoey. Licensed under [CC BY 4.0 Creative Commons](http://creativecommons.org/licenses/by/4.0/)*\n", + "\n", + "---" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": { + "deletable": true, + "editable": true + }, + "outputs": [], + "source": [ + "import numpy as np\n", + "import matplotlib.pyplot as plt\n", + "\n", + "import xarray as xr" + ] + }, + { + "cell_type": "markdown", + "metadata": { + "deletable": true, + "editable": true + }, + "source": [ + "## Meet `xarray.Dataset`" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "We already know `xarray.DataArray`: a single array with dimensions (and each dimension can have a name and coordinate values):" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "<xarray.DataArray (band: 2, y: 317, x: 625)>\n", + "[396250 values with dtype=float64]\n", + "Coordinates:\n", + " * band (band) int64 1 2\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0, -57.4483...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan)\n", + " scales: (1.0, 1.0)\n", + " offsets: (0.0, 0.0)\n", + " AREA_OR_POINT: Area
[396250 values with dtype=float64]
array([1, 2])
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
<xarray.Dataset>\n", + "Dimensions: (x: 625, y: 317)\n", + "Coordinates:\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Data variables:\n", + " b4 (y, x) float64 ...\n", + " b8 (y, x) float64 ...\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0, -57.4483...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan)\n", + " scales: (1.0, 1.0)\n", + " offsets: (0.0, 0.0)\n", + " AREA_OR_POINT: Area
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
[198125 values with dtype=float64]
[198125 values with dtype=float64]
<xarray.DataArray (band: 2, y: 317, x: 625)>\n", + "[396250 values with dtype=float64]\n", + "Coordinates:\n", + " * band (band) int64 1 2\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0, -57.4483...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan)\n", + " scales: (1.0, 1.0)\n", + " offsets: (0.0, 0.0)\n", + " AREA_OR_POINT: Area
[396250 values with dtype=float64]
array([1, 2])
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
<xarray.Dataset>\n", + "Dimensions: (x: 625, y: 317)\n", + "Coordinates:\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Data variables:\n", + " b4 (y, x) float64 0.0782 0.0884 0.1751 ... 0.05769 0.0676 0.08299\n", + " b8 (y, x) float64 0.2736 0.257 0.2479 0.2257 ... 0.3324 0.3388 0.3722\n", + " ndvi (y, x) float64 0.5554 0.4881 0.1721 0.153 ... 0.7042 0.6673 0.6353\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0, -57.4483...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan)\n", + " scales: (1.0, 1.0)\n", + " offsets: (0.0, 0.0)\n", + " AREA_OR_POINT: Area
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
array([[0.078202, 0.088396, 0.175097, ..., 0.103395, 0.104005, 0.086397],\n", + " [0.102205, 0.075898, 0.085496, ..., 0.073198, 0.065705, 0.065507],\n", + " [0.101106, 0.080003, 0.076997, ..., 0.068803, 0.073304, 0.070497],\n", + " ...,\n", + " [0.072694, 0.074205, 0.068101, ..., 0.140902, 0.100801, 0.068299],\n", + " [0.070298, 0.069505, 0.067201, ..., 0.073396, 0.089006, 0.059907],\n", + " [0.070603, 0.0748 , 0.072404, ..., 0.057694, 0.067597, 0.082994]])
array([[0.273594, 0.256992, 0.247898, ..., 0.283696, 0.245197, 0.232105],\n", + " [0.255695, 0.231998, 0.259693, ..., 0.2551 , 0.223301, 0.231006],\n", + " [0.226596, 0.242405, 0.267903, ..., 0.271794, 0.281193, 0.309697],\n", + " ...,\n", + " [0.269596, 0.253407, 0.291997, ..., 0.245594, 0.291707, 0.308995],\n", + " [0.2822 , 0.2925 , 0.293599, ..., 0.213993, 0.262196, 0.379797],\n", + " [0.288502, 0.248096, 0.261906, ..., 0.332403, 0.338796, 0.372198]])
array([[0.55541097, 0.48813784, 0.17210779, ..., 0.46578366, 0.40432598,\n", + " 0.45748096],\n", + " [0.42886378, 0.50698781, 0.5046415 , ..., 0.55407855, 0.54530095,\n", + " 0.5581515 ],\n", + " [0.38293909, 0.50371527, 0.5535106 , ..., 0.59598584, 0.58643251,\n", + " 0.62915396],\n", + " ...,\n", + " [0.57524964, 0.54699581, 0.62176363, ..., 0.27087528, 0.48637406,\n", + " 0.63795195],\n", + " [0.60114281, 0.61600067, 0.62748996, ..., 0.48922162, 0.49313521,\n", + " 0.72751249],\n", + " [0.60678168, 0.53669486, 0.56684468, ..., 0.70420497, 0.66733 ,\n", + " 0.63534578]])
<xarray.Dataset>\n", + "Dimensions: (date: 75, level: 78)\n", + "Coordinates:\n", + " * level (level) int32 0 1 2 3 4 5 6 7 8 ... 69 70 71 72 73 74 75 76 77\n", + " * date (date) datetime64[ns] 2012-07-13T22:33:06.019200 ... 2014-07...\n", + " lon (date) float64 -39.13 -37.28 -36.9 ... -33.83 -34.11 -34.38\n", + " lat (date) float64 47.19 46.72 46.45 46.23 ... 42.6 42.46 42.38\n", + "Data variables:\n", + " salinity (level, date) float64 35.64 35.51 35.57 35.4 ... nan 34.94 nan\n", + " temperature (level, date) float64 18.97 18.44 19.1 19.79 ... nan 3.714 nan\n", + " pressure (level, date) float64 6.8 6.1 6.5 5.0 ... 2e+03 nan 2e+03 nan
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,\n", + " 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,\n", + " 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,\n", + " 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,\n", + " 72, 73, 74, 75, 76, 77], dtype=int32)
array(['2012-07-13T22:33:06.019200000', '2012-07-23T22:54:59.990400000',\n", + " '2012-08-02T22:55:52.003200000', '2012-08-12T23:08:59.971200000',\n", + " '2012-08-22T23:29:01.968000000', '2012-09-01T23:17:38.976000000',\n", + " '2012-09-12T02:59:18.960000000', '2012-09-21T23:18:37.036800000',\n", + " '2012-10-02T03:00:17.971200000', '2012-10-11T23:13:27.984000000',\n", + " '2012-10-22T02:50:32.006399999', '2012-10-31T23:36:39.974400000',\n", + " '2012-11-11T02:40:46.041600000', '2012-11-20T23:08:29.990400000',\n", + " '2012-12-01T02:47:51.993600000', '2012-12-10T23:23:16.972800000',\n", + " '2012-12-21T02:58:48.979200000', '2012-12-30T23:07:23.030400000',\n", + " '2013-01-10T02:56:43.008000000', '2013-01-19T23:24:26.956800000',\n", + " '2013-01-30T02:43:53.011199998', '2013-02-08T23:15:27.043200000',\n", + " '2013-02-19T01:12:50.976000000', '2013-02-28T23:07:13.008000000',\n", + " '2013-03-11T02:43:30.979200000', '2013-03-20T23:17:22.992000000',\n", + " '2013-03-31T01:50:38.025600000', '2013-04-09T23:19:07.968000000',\n", + " '2013-04-20T02:53:29.990400000', '2013-04-29T23:28:33.024000000',\n", + " '2013-05-10T02:50:18.009600000', '2013-05-19T23:21:05.990400000',\n", + " '2013-05-30T02:50:30.969600000', '2013-06-08T23:32:49.027200000',\n", + " '2013-06-19T03:42:51.004800000', '2013-06-28T23:32:16.022400000',\n", + " '2013-07-09T03:28:30.979200000', '2013-07-18T23:33:57.974400000',\n", + " '2013-07-29T03:15:42.019200000', '2013-08-07T23:25:02.035200004',\n", + " '2013-08-18T01:47:44.966400000', '2013-08-28T03:02:59.020799996',\n", + " '2013-09-07T03:03:51.984000000', '2013-09-16T23:32:07.987200000',\n", + " '2013-09-27T04:08:27.974400000', '2013-10-06T23:25:39.964800000',\n", + " '2013-10-17T02:55:50.995200000', '2013-10-27T03:45:47.001600000',\n", + " '2013-11-06T01:14:52.022400000', '2013-11-16T03:29:54.009600000',\n", + " '2013-11-26T03:03:56.995200000', '2013-12-05T23:33:59.011200000',\n", + " '2013-12-16T02:58:01.977600000', '2013-12-25T23:22:43.017600000',\n", + " '2014-01-05T02:52:06.009600000', '2014-01-14T23:41:18.009600000',\n", + " '2014-01-25T03:00:43.977600000', '2014-02-03T23:29:13.977600000',\n", + " '2014-02-14T02:50:11.961600000', '2014-02-23T23:03:21.974400000',\n", + " '2014-03-06T02:58:03.964800000', '2014-03-15T23:10:28.012800000',\n", + " '2014-03-26T02:51:22.032000000', '2014-04-04T23:25:58.972800000',\n", + " '2014-04-15T03:00:45.964800000', '2014-04-24T23:24:40.003200000',\n", + " '2014-05-05T02:56:22.012800000', '2014-05-15T00:10:06.009600000',\n", + " '2014-05-25T03:02:43.036800000', '2014-06-03T23:34:53.961600000',\n", + " '2014-06-14T03:01:23.980800000', '2014-06-23T23:24:31.968000000',\n", + " '2014-07-04T03:08:30.019200000', '2014-07-13T23:47:43.008000000',\n", + " '2014-07-24T03:02:33.014400000'], dtype='datetime64[ns]')
array([-39.13 , -37.282, -36.9 , -36.89 , -37.053, -36.658, -35.963,\n", + " -35.184, -34.462, -33.784, -32.972, -32.546, -32.428, -32.292,\n", + " -32.169, -31.998, -31.824, -31.624, -31.433, -31.312, -31.107,\n", + " -31.147, -31.044, -31.14 , -31.417, -31.882, -32.145, -32.487,\n", + " -32.537, -32.334, -32.042, -31.892, -31.861, -31.991, -31.883,\n", + " -31.89 , -31.941, -31.889, -31.724, -31.412, -31.786, -31.561,\n", + " -31.732, -31.553, -31.862, -32.389, -32.318, -32.19 , -32.224,\n", + " -32.368, -32.306, -32.305, -32.65 , -33.093, -33.263, -33.199,\n", + " -33.27 , -33.237, -33.221, -33.011, -32.844, -32.981, -32.784,\n", + " -32.607, -32.87 , -33.196, -33.524, -33.956, -33.944, -33.71 ,\n", + " -33.621, -33.552, -33.828, -34.11 , -34.38 ])
array([47.187, 46.716, 46.45 , 46.23 , 45.459, 44.833, 44.452, 44.839,\n", + " 44.956, 44.676, 44.13 , 43.644, 43.067, 42.662, 42.513, 42.454,\n", + " 42.396, 42.256, 42.089, 41.944, 41.712, 41.571, 41.596, 41.581,\n", + " 41.351, 41.032, 40.912, 40.792, 40.495, 40.383, 40.478, 40.672,\n", + " 41.032, 40.864, 40.651, 40.425, 40.228, 40.197, 40.483, 40.311,\n", + " 40.457, 40.463, 40.164, 40.047, 39.963, 40.122, 40.57 , 40.476,\n", + " 40.527, 40.589, 40.749, 40.993, 41.162, 41.237, 41.448, 41.65 ,\n", + " 42.053, 42.311, 42.096, 41.683, 41.661, 41.676, 42.018, 42.395,\n", + " 42.532, 42.558, 42.504, 42.63 , 42.934, 42.952, 42.777, 42.722,\n", + " 42.601, 42.457, 42.379])
array([[35.6389389 , 35.51495743, 35.57297134, ..., 35.82093811,\n", + " 35.77793884, 35.66891098],\n", + " [35.63393784, 35.5219574 , 35.57397079, ..., 35.81093216,\n", + " 35.58389664, 35.66791153],\n", + " [35.6819458 , 35.52595901, 35.57297134, ..., 35.79592896,\n", + " 35.66290665, 35.66591263],\n", + " ...,\n", + " [34.91585922, 34.92390442, 34.92390442, ..., 34.93481064,\n", + " 34.94081116, 34.94680786],\n", + " [34.91585922, 34.92390442, 34.92190552, ..., 34.93280792,\n", + " 34.93680954, 34.94380951],\n", + " [34.91785812, 34.92390442, 34.92390442, ..., nan,\n", + " 34.93680954, nan]])
array([[18.97400093, 18.43700027, 19.09900093, ..., 19.11300087,\n", + " 21.82299995, 20.13100052],\n", + " [18.74099922, 18.39999962, 19.08200073, ..., 18.47200012,\n", + " 19.45999908, 20.125 ],\n", + " [18.37000084, 18.37400055, 19.06500053, ..., 18.22999954,\n", + " 19.26199913, 20.07699966],\n", + " ...,\n", + " [ 3.79299998, 3.81399989, 3.80200005, ..., 3.80699992,\n", + " 3.81100011, 3.8599999 ],\n", + " [ 3.76399994, 3.77800012, 3.75699997, ..., 3.75399995,\n", + " 3.74600005, 3.80599999],\n", + " [ 3.74399996, 3.74600005, 3.7249999 , ..., nan,\n", + " 3.71399999, nan]])
array([[ 6.80000019, 6.0999999 , 6.5 , ..., 7.0999999 ,\n", + " 7.20000029, 6.5 ],\n", + " [ 10.69999981, 10.59999943, 10.39999962, ..., 10.79999924,\n", + " 11.09999943, 10.39999962],\n", + " [ 15.69999981, 14.09999943, 14.89999962, ..., 15.89999962,\n", + " 15.59999943, 15.89999962],\n", + " ...,\n", + " [1900.60009766, 1900. , 1900.19995117, ..., 1899.70007324,\n", + " 1900.40002441, 1899.80004883],\n", + " [1949.90002441, 1950. , 1949.89990234, ..., 1950.59997559,\n", + " 1950.20007324, 1949.70007324],\n", + " [1999.30004883, 1998. , 1998.5 , ..., nan,\n", + " 2000.40002441, nan]])
<xarray.Dataset>\n", + "Dimensions: (level: 78)\n", + "Coordinates:\n", + " * level (level) int32 0 1 2 3 4 5 6 7 8 ... 69 70 71 72 73 74 75 76 77\n", + "Data variables:\n", + " salinity (level) float64 35.91 35.9 35.9 35.9 ... 34.94 34.94 34.93\n", + " temperature (level) float64 17.6 17.57 17.51 17.42 ... 3.789 3.73 3.662\n", + " pressure (level) float64 6.435 10.57 15.54 ... 1.95e+03 1.999e+03
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,\n", + " 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,\n", + " 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,\n", + " 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,\n", + " 72, 73, 74, 75, 76, 77], dtype=int32)
array([35.9063218 , 35.90223138, 35.90313435, 35.90173139, 35.90544583,\n", + " 35.9100359 , 35.90946015, 35.91343146, 35.91967712, 35.92615988,\n", + " 35.93195456, 35.94055356, 35.94091596, 35.93905366, 35.93931069,\n", + " 35.93786745, 35.93525794, 35.93118039, 35.92534328, 35.91652257,\n", + " 35.90671895, 35.89617843, 35.88888019, 35.8789927 , 35.86946183,\n", + " 35.8598671 , 35.85061713, 35.84211978, 35.83150467, 35.81969395,\n", + " 35.80945061, 35.80092265, 35.79078674, 35.77886525, 35.76833627,\n", + " 35.75838795, 35.74923783, 35.73923559, 35.73000444, 35.71877237,\n", + " 35.68864513, 35.65607159, 35.62678265, 35.59231774, 35.56205662,\n", + " 35.45401408, 35.41392634, 35.3810557 , 35.34845245, 35.31531555,\n", + " 35.28392568, 35.26568334, 35.2389473 , 35.21583745, 35.19686081,\n", + " 35.18231257, 35.1648436 , 35.15073542, 35.12509338, 35.10155869,\n", + " 35.08199799, 35.06317012, 35.0490097 , 35.03678253, 35.02174266,\n", + " 35.01135579, 35.00212936, 34.99386297, 34.98810328, 34.98008094,\n", + " 34.97214884, 34.96517645, 34.95664983, 34.9507985 , 34.9465696 ,\n", + " 34.94198907, 34.93844852, 34.93290652])
array([17.60172602, 17.57223609, 17.5145833 , 17.42326395, 17.24943838,\n", + " 17.03730134, 16.76787661, 16.44609588, 16.17439195, 16.04501356,\n", + " 15.65827023, 15.4607296 , 15.26114862, 15.12489191, 14.99133783,\n", + " 14.90160808, 14.81990544, 14.74535139, 14.66822971, 14.585027 ,\n", + " 14.49732434, 14.41904053, 14.35412163, 14.27102702, 14.19081082,\n", + " 14.11487838, 14.04347293, 13.98067566, 13.90994595, 13.83274319,\n", + " 13.76139196, 13.69836479, 13.62335132, 13.54185131, 13.46647295,\n", + " 13.39395946, 13.32541891, 13.25205403, 13.18131082, 13.10233782,\n", + " 12.89268916, 12.67795943, 12.4649189 , 12.2178513 , 11.98270268,\n", + " 11.1281081 , 10.80430666, 10.49702667, 10.1749066 , 9.83453334,\n", + " 9.48625332, 9.19793334, 8.66010666, 8.12324001, 7.60221333,\n", + " 7.15289333, 6.74250667, 6.39543999, 6.04598667, 5.74538665,\n", + " 5.48913333, 5.26604001, 5.08768 , 4.93479998, 4.77769334,\n", + " 4.65368 , 4.54237334, 4.44274664, 4.35933333, 4.2653784 ,\n", + " 4.17290539, 4.08902703, 3.99864865, 3.92163514, 3.85617567,\n", + " 3.78916217, 3.72950001, 3.66207691])
array([ 6.43466671, 10.56891882, 15.54246568, 20.46301361,\n", + " 25.42567552, 30.44459441, 35.44324375, 40.4391894 ,\n", + " 45.40810832, 50.37837879, 60.47297323, 70.48378413,\n", + " 80.40270347, 90.48243311, 100.51216311, 110.46081151,\n", + " 120.52702795, 130.49459282, 140.51216064, 150.40540376,\n", + " 160.40810559, 170.36216035, 180.41080949, 190.4108097 ,\n", + " 200.39999761, 210.34729499, 220.32026858, 230.31351224,\n", + " 240.28918808, 250.41486297, 260.39999843, 270.36891752,\n", + " 280.42432136, 290.42297075, 300.4229691 , 310.46351087,\n", + " 320.50675346, 330.5297266 , 340.41891521, 350.49729383,\n", + " 375.41080867, 400.3797294 , 425.29864626, 450.38378205,\n", + " 475.30675403, 550.47703016, 575.68400146, 600.42400716,\n", + " 625.30800456, 650.34533773, 675.33333984, 700.37067546,\n", + " 750.42400716, 800.36666992, 850.38534017, 900.4613387 ,\n", + " 950.45067383, 1000.38534261, 1050.38534668, 1100.45734212,\n", + " 1150.45201335, 1200.40534505, 1250.25067383, 1300.49467773,\n", + " 1350.40268392, 1400.41734538, 1450.25734212, 1500.40267253,\n", + " 1550.46401367, 1600.44055011, 1650.34595056, 1700.42163251,\n", + " 1750.42299013, 1800.35677358, 1850.33379467, 1900.20676401,\n", + " 1950.20946771, 1999.28975423])
<xarray.DataArray 'month' (date: 75)>\n", + "array([ 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 10, 11, 11, 12, 12, 12,\n", + " 12, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6,\n", + " 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 11, 11, 11,\n", + " 12, 12, 12, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5,\n", + " 5, 6, 6, 6, 7, 7, 7])\n", + "Coordinates:\n", + " * date (date) datetime64[ns] 2012-07-13T22:33:06.019200 ... 2014-07-24T...\n", + " lon (date) float64 -39.13 -37.28 -36.9 -36.89 ... -33.83 -34.11 -34.38\n", + " lat (date) float64 47.19 46.72 46.45 46.23 ... 42.72 42.6 42.46 42.38
array([ 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 10, 11, 11, 12, 12, 12,\n", + " 12, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6,\n", + " 6, 6, 7, 7, 7, 8, 8, 8, 9, 9, 9, 10, 10, 10, 11, 11, 11,\n", + " 12, 12, 12, 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5,\n", + " 5, 6, 6, 6, 7, 7, 7])
array(['2012-07-13T22:33:06.019200000', '2012-07-23T22:54:59.990400000',\n", + " '2012-08-02T22:55:52.003200000', '2012-08-12T23:08:59.971200000',\n", + " '2012-08-22T23:29:01.968000000', '2012-09-01T23:17:38.976000000',\n", + " '2012-09-12T02:59:18.960000000', '2012-09-21T23:18:37.036800000',\n", + " '2012-10-02T03:00:17.971200000', '2012-10-11T23:13:27.984000000',\n", + " '2012-10-22T02:50:32.006399999', '2012-10-31T23:36:39.974400000',\n", + " '2012-11-11T02:40:46.041600000', '2012-11-20T23:08:29.990400000',\n", + " '2012-12-01T02:47:51.993600000', '2012-12-10T23:23:16.972800000',\n", + " '2012-12-21T02:58:48.979200000', '2012-12-30T23:07:23.030400000',\n", + " '2013-01-10T02:56:43.008000000', '2013-01-19T23:24:26.956800000',\n", + " '2013-01-30T02:43:53.011199998', '2013-02-08T23:15:27.043200000',\n", + " '2013-02-19T01:12:50.976000000', '2013-02-28T23:07:13.008000000',\n", + " '2013-03-11T02:43:30.979200000', '2013-03-20T23:17:22.992000000',\n", + " '2013-03-31T01:50:38.025600000', '2013-04-09T23:19:07.968000000',\n", + " '2013-04-20T02:53:29.990400000', '2013-04-29T23:28:33.024000000',\n", + " '2013-05-10T02:50:18.009600000', '2013-05-19T23:21:05.990400000',\n", + " '2013-05-30T02:50:30.969600000', '2013-06-08T23:32:49.027200000',\n", + " '2013-06-19T03:42:51.004800000', '2013-06-28T23:32:16.022400000',\n", + " '2013-07-09T03:28:30.979200000', '2013-07-18T23:33:57.974400000',\n", + " '2013-07-29T03:15:42.019200000', '2013-08-07T23:25:02.035200004',\n", + " '2013-08-18T01:47:44.966400000', '2013-08-28T03:02:59.020799996',\n", + " '2013-09-07T03:03:51.984000000', '2013-09-16T23:32:07.987200000',\n", + " '2013-09-27T04:08:27.974400000', '2013-10-06T23:25:39.964800000',\n", + " '2013-10-17T02:55:50.995200000', '2013-10-27T03:45:47.001600000',\n", + " '2013-11-06T01:14:52.022400000', '2013-11-16T03:29:54.009600000',\n", + " '2013-11-26T03:03:56.995200000', '2013-12-05T23:33:59.011200000',\n", + " '2013-12-16T02:58:01.977600000', '2013-12-25T23:22:43.017600000',\n", + " '2014-01-05T02:52:06.009600000', '2014-01-14T23:41:18.009600000',\n", + " '2014-01-25T03:00:43.977600000', '2014-02-03T23:29:13.977600000',\n", + " '2014-02-14T02:50:11.961600000', '2014-02-23T23:03:21.974400000',\n", + " '2014-03-06T02:58:03.964800000', '2014-03-15T23:10:28.012800000',\n", + " '2014-03-26T02:51:22.032000000', '2014-04-04T23:25:58.972800000',\n", + " '2014-04-15T03:00:45.964800000', '2014-04-24T23:24:40.003200000',\n", + " '2014-05-05T02:56:22.012800000', '2014-05-15T00:10:06.009600000',\n", + " '2014-05-25T03:02:43.036800000', '2014-06-03T23:34:53.961600000',\n", + " '2014-06-14T03:01:23.980800000', '2014-06-23T23:24:31.968000000',\n", + " '2014-07-04T03:08:30.019200000', '2014-07-13T23:47:43.008000000',\n", + " '2014-07-24T03:02:33.014400000'], dtype='datetime64[ns]')
array([-39.13 , -37.282, -36.9 , -36.89 , -37.053, -36.658, -35.963,\n", + " -35.184, -34.462, -33.784, -32.972, -32.546, -32.428, -32.292,\n", + " -32.169, -31.998, -31.824, -31.624, -31.433, -31.312, -31.107,\n", + " -31.147, -31.044, -31.14 , -31.417, -31.882, -32.145, -32.487,\n", + " -32.537, -32.334, -32.042, -31.892, -31.861, -31.991, -31.883,\n", + " -31.89 , -31.941, -31.889, -31.724, -31.412, -31.786, -31.561,\n", + " -31.732, -31.553, -31.862, -32.389, -32.318, -32.19 , -32.224,\n", + " -32.368, -32.306, -32.305, -32.65 , -33.093, -33.263, -33.199,\n", + " -33.27 , -33.237, -33.221, -33.011, -32.844, -32.981, -32.784,\n", + " -32.607, -32.87 , -33.196, -33.524, -33.956, -33.944, -33.71 ,\n", + " -33.621, -33.552, -33.828, -34.11 , -34.38 ])
array([47.187, 46.716, 46.45 , 46.23 , 45.459, 44.833, 44.452, 44.839,\n", + " 44.956, 44.676, 44.13 , 43.644, 43.067, 42.662, 42.513, 42.454,\n", + " 42.396, 42.256, 42.089, 41.944, 41.712, 41.571, 41.596, 41.581,\n", + " 41.351, 41.032, 40.912, 40.792, 40.495, 40.383, 40.478, 40.672,\n", + " 41.032, 40.864, 40.651, 40.425, 40.228, 40.197, 40.483, 40.311,\n", + " 40.457, 40.463, 40.164, 40.047, 39.963, 40.122, 40.57 , 40.476,\n", + " 40.527, 40.589, 40.749, 40.993, 41.162, 41.237, 41.448, 41.65 ,\n", + " 42.053, 42.311, 42.096, 41.683, 41.661, 41.676, 42.018, 42.395,\n", + " 42.532, 42.558, 42.504, 42.63 , 42.934, 42.952, 42.777, 42.722,\n", + " 42.601, 42.457, 42.379])
<xarray.Dataset>\n", + "Dimensions: (level: 78, month: 12)\n", + "Coordinates:\n", + " * level (level) int32 0 1 2 3 4 5 6 7 8 ... 69 70 71 72 73 74 75 76 77\n", + " * month (month) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", + "Data variables:\n", + " salinity (month, level) float64 35.95 35.95 35.95 ... 34.94 34.94 34.93\n", + " temperature (month, level) float64 15.06 15.06 15.06 ... 3.779 3.733 3.663\n", + " pressure (month, level) float64 6.6 10.42 15.32 ... 1.95e+03 1.998e+03
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,\n", + " 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,\n", + " 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,\n", + " 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,\n", + " 72, 73, 74, 75, 76, 77], dtype=int32)
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12])
array([[35.94653575, 35.94653511, 35.9463679 , 35.94636854, 35.944369 ,\n", + " 35.94387054, 35.94136937, 35.93703651, 35.93436941, 35.9337031 ,\n", + " 35.93253581, 35.92838605, 35.92570241, 35.92271932, 35.92321905,\n", + " 35.91988691, 35.90888468, 35.89755122, 35.89805158, 35.90155156,\n", + " 35.89838473, 35.89155134, 35.89071782, 35.88738505, 35.88840294,\n", + " 35.8868866 , 35.8815848 , 35.87410037, 35.86493174, 35.85942904,\n", + " 35.85442924, 35.84742546, 35.83709145, 35.82842255, 35.8182532 ,\n", + " 35.81408564, 35.80641619, 35.79924901, 35.79274686, 35.78724543,\n", + " 35.75673993, 35.71923192, 35.67922592, 35.63004875, 35.60821215,\n", + " 35.47102674, 35.42835172, 35.37917964, 35.34684308, 35.31467247,\n", + " 35.3026708 , 35.2681675 , 35.2131621 , 35.21466192, 35.22882652,\n", + " 35.21299235, 35.19999059, 35.19215457, 35.14898618, 35.12281736,\n", + " 35.09448306, 35.05331421, 35.02397982, 35.0021464 , 34.97881317,\n", + " 34.97431374, 34.96848106, 34.96398099, 34.96398163, 34.95931498,\n", + " 34.95614815, 34.95331637, 34.94731712, 34.94448344, 34.94148445,\n", + " 34.93915113, 34.93648402, 34.93235588],\n", + " [35.89972496, 35.90005811, 35.89989217, 35.90039253, 35.89972496,\n", + " 35.90005938, 35.90122541, 35.90122604, 35.90105947, 35.90022596,\n", + " 35.90022723, 35.89639155, 35.89455986, 35.89272499, 35.89322599,\n", + " 35.89139175, 35.88905907, 35.88672574, 35.88305791, 35.88005829,\n", + "...\n", + " 35.10666733, 35.07386475, 35.06346512, 35.05146484, 35.0374649 ,\n", + " 35.02506638, 35.01266632, 35.00106583, 34.98906479, 34.98726425,\n", + " 34.97906494, 34.97006531, 34.95726624, 34.95106506, 34.94786606,\n", + " 34.94466705, 34.94266663, 34.92985725],\n", + " [36.0026652 , 36.00316556, 36.00333277, 36.00283305, 36.01826804,\n", + " 36.01855414, 36.01869638, 36.0191258 , 36.01926858, 36.01898248,\n", + " 36.01998247, 36.0188397 , 36.01526751, 36.00183868, 36.00869533,\n", + " 36.01626587, 36.02469526, 36.01412255, 36.00483595, 35.9976932 ,\n", + " 35.98940768, 35.97297832, 35.96712058, 35.95684814, 35.93784768,\n", + " 35.92327554, 35.90417371, 35.89404188, 35.88547134, 35.87091119,\n", + " 35.86049543, 35.85292053, 35.84177617, 35.83077294, 35.81791251,\n", + " 35.80876814, 35.80062376, 35.79119274, 35.78061894, 35.77347401,\n", + " 35.74446869, 35.70917565, 35.66716712, 35.63859177, 35.61030143,\n", + " 35.51071439, 35.47527913, 35.4458471 , 35.41055679, 35.3728387 ,\n", + " 35.32597569, 35.30368696, 35.28382874, 35.24739511, 35.22581972,\n", + " 35.2198203 , 35.18624496, 35.20252991, 35.15709741, 35.11523819,\n", + " 35.09195164, 35.06366457, 35.03737804, 35.01394817, 34.99623544,\n", + " 34.98466219, 34.97823497, 34.97309276, 34.96837779, 34.96066448,\n", + " 34.95609283, 34.94966398, 34.94437844, 34.94252232, 34.94066674,\n", + " 34.93838228, 34.93666676, 34.93206024]])
array([[15.05599976, 15.05833324, 15.05766662, 15.06166681, 15.05516672,\n", + " 15.04283333, 15.03349972, 15.01066653, 14.99283346, 14.99166663,\n", + " 14.98616664, 14.96616666, 14.95250003, 14.93933344, 14.94133313,\n", + " 14.92733335, 14.86250003, 14.7711668 , 14.75533326, 14.7420001 ,\n", + " 14.70083332, 14.66266648, 14.64016644, 14.54766687, 14.49650002,\n", + " 14.4278334 , 14.36633317, 14.30000003, 14.22200012, 14.16349983,\n", + " 14.08966668, 14.02549998, 13.92916663, 13.85066668, 13.77833335,\n", + " 13.74049981, 13.66950003, 13.61616659, 13.55666685, 13.51116673,\n", + " 13.3191665 , 13.10066652, 12.84433301, 12.52416658, 12.3496666 ,\n", + " 11.38833348, 11.02249988, 10.60700003, 10.29350011, 9.93483321,\n", + " 9.60100015, 9.24199994, 8.55033318, 8.09533326, 7.72949998,\n", + " 7.28616667, 6.88616673, 6.58066662, 6.13900002, 5.8378334 ,\n", + " 5.54400007, 5.21483334, 4.96049992, 4.75650001, 4.53549997,\n", + " 4.44733326, 4.34666665, 4.27249996, 4.21516657, 4.14200004,\n", + " 4.07133325, 4.00899998, 3.94000002, 3.87983338, 3.82100002,\n", + " 3.76816662, 3.72016664, 3.64350009],\n", + " [14.43433348, 14.43516684, 14.43666681, 14.43650007, 14.43683322,\n", + " 14.4388334 , 14.43983332, 14.43833319, 14.43016656, 14.42850002,\n", + " 14.4236668 , 14.40166664, 14.39283339, 14.38700008, 14.38533322,\n", + " 14.37683344, 14.36866665, 14.35533333, 14.34200017, 14.32849995,\n", + "...\n", + " 5.66700001, 5.34379997, 5.17019997, 5.00860004, 4.86059999,\n", + " 4.7368001 , 4.61699991, 4.50999994, 4.38580008, 4.32800002,\n", + " 4.227 , 4.13199997, 4.01800003, 3.93739996, 3.88059998,\n", + " 3.81880002, 3.75939999, 3.64700007],\n", + " [16.3378334 , 16.33916648, 16.34100008, 16.33999983, 16.44442872,\n", + " 16.44557149, 16.44699996, 16.44842843, 16.449143 , 16.44514302,\n", + " 16.45071398, 16.44457122, 16.43214308, 16.29371452, 15.99871458,\n", + " 15.77285712, 15.58657128, 15.44085721, 15.34957136, 15.26457119,\n", + " 15.15414279, 15.03057126, 14.97799982, 14.85214288, 14.69985717,\n", + " 14.59171418, 14.46257142, 14.38071428, 14.32242843, 14.21785709,\n", + " 14.14185715, 14.08314283, 14.00599984, 13.9140001 , 13.82471439,\n", + " 13.74700001, 13.67157146, 13.59100001, 13.51314286, 13.45285715,\n", + " 13.25885705, 13.04057135, 12.76942853, 12.54028566, 12.32628577,\n", + " 11.52057144, 11.22199999, 10.95385715, 10.64757116, 10.29900006,\n", + " 9.86400005, 9.55828571, 9.06914289, 8.46785722, 7.93571431,\n", + " 7.50214291, 6.95328569, 6.74199997, 6.24899994, 5.85328565,\n", + " 5.56214292, 5.30028575, 5.0560001 , 4.83942849, 4.65957144,\n", + " 4.52271427, 4.42871428, 4.33642851, 4.26228571, 4.16785717,\n", + " 4.09142859, 4.01014277, 3.9365715 , 3.883 , 3.83199995,\n", + " 3.77928577, 3.73271431, 3.66299999]])
array([[ 6.59999998, 10.41666651, 15.3166666 , 20.29999955,\n", + " 25.29999987, 30.29999987, 35.30000051, 40.33333397,\n", + " 45.2833341 , 50.31666629, 60.20000013, 70.21666845,\n", + " 80.36666743, 90.63333384, 100.5 , 110.61666743,\n", + " 120.6166687 , 130.53333028, 140.48332723, 150.51666514,\n", + " 160.53333282, 170.43333181, 180.53333028, 190.63332876,\n", + " 200.43333181, 210.16666412, 220.31666565, 230.36666361,\n", + " 240.28333282, 250.04999797, 259.7666626 , 270.04999797,\n", + " 280.5333252 , 290.61666361, 300.6499939 , 310.66666158,\n", + " 320.58332316, 330.5333252 , 340.44999186, 350.41666158,\n", + " 375.31667074, 400.51666768, 425.41665649, 450.5 ,\n", + " 475.3833313 , 550.08333333, 575.68334961, 600.3666687 ,\n", + " 625.33334351, 650.41667684, 675.56667074, 700.45002238,\n", + " 750.58334351, 800.18332926, 850.65001424, 900.65000407,\n", + " 950.46668498, 1000.48332723, 1050.43332926, 1100.50002035,\n", + " 1150.5 , 1200.51668294, 1249.86667887, 1300.45001221,\n", + " 1350.51668294, 1400.38336182, 1450.45001221, 1500. ,\n", + " 1550.63334147, 1600.45003255, 1650.5 , 1700.55000814,\n", + " 1750.56669108, 1800.48335775, 1850.28334554, 1900.11665853,\n", + " 1950.05000814, 1999.95001221],\n", + "...\n", + " [ 6.42857143, 10.55714308, 15.33333349, 20.35000006,\n", + " 25.34285709, 30.35714258, 35.35714177, 40.35714231,\n", + " 45.35714286, 50.4857145 , 60.2714277 , 70.17142923,\n", + " 80.20000131, 90.7285723 , 100.81428637, 110.71428571,\n", + " 120.59999956, 130.62857274, 140.61428833, 150.51428877,\n", + " 160.51429095, 170.47143119, 180.51428659, 190.41429138,\n", + " 200.48571777, 210.62857274, 220.58571734, 230.32857186,\n", + " 239.81428746, 250.80000305, 260.55714199, 270.65714373,\n", + " 280.44285366, 290.2714277 , 300.59999302, 310.18571254,\n", + " 320.57142857, 330.68571254, 340.67142596, 350.5857108 ,\n", + " 375.30000087, 400.28571865, 425.05713763, 450.54285104,\n", + " 475.0857108 , 550.54284668, 575.49999128, 600.34284319,\n", + " 625.38570731, 650.41427176, 675.52856445, 700.49999128,\n", + " 750.60000174, 800.3999721 , 850.62855748, 900.54284668,\n", + " 950.64285714, 1000.29998779, 1050.52858189, 1100.44285366,\n", + " 1150.48573521, 1200.54288156, 1249.97141811, 1300.47143555,\n", + " 1350.50001744, 1400.15713937, 1450.29999651, 1500.5 ,\n", + " 1550.64285714, 1600.35714286, 1650.21428571, 1700.0857108 ,\n", + " 1750.51428223, 1800.40002441, 1850.51428223, 1900.41430664,\n", + " 1950.30001395, 1998.1499939 ]])
<xarray.Dataset>\n", + "Dimensions: (date: 25, level: 78)\n", + "Coordinates:\n", + " * date (date) datetime64[ns] 2012-07-31 2012-08-31 ... 2014-07-31\n", + " * level (level) int32 0 1 2 3 4 5 6 7 8 ... 69 70 71 72 73 74 75 76 77\n", + "Data variables:\n", + " salinity (date, level) float64 35.58 35.58 35.6 ... 34.94 34.94 34.94\n", + " temperature (date, level) float64 18.71 18.57 18.37 ... 3.811 3.754 3.714\n", + " pressure (date, level) float64 6.45 10.65 14.9 ... 1.95e+03 2e+03
array(['2012-07-31T00:00:00.000000000', '2012-08-31T00:00:00.000000000',\n", + " '2012-09-30T00:00:00.000000000', '2012-10-31T00:00:00.000000000',\n", + " '2012-11-30T00:00:00.000000000', '2012-12-31T00:00:00.000000000',\n", + " '2013-01-31T00:00:00.000000000', '2013-02-28T00:00:00.000000000',\n", + " '2013-03-31T00:00:00.000000000', '2013-04-30T00:00:00.000000000',\n", + " '2013-05-31T00:00:00.000000000', '2013-06-30T00:00:00.000000000',\n", + " '2013-07-31T00:00:00.000000000', '2013-08-31T00:00:00.000000000',\n", + " '2013-09-30T00:00:00.000000000', '2013-10-31T00:00:00.000000000',\n", + " '2013-11-30T00:00:00.000000000', '2013-12-31T00:00:00.000000000',\n", + " '2014-01-31T00:00:00.000000000', '2014-02-28T00:00:00.000000000',\n", + " '2014-03-31T00:00:00.000000000', '2014-04-30T00:00:00.000000000',\n", + " '2014-05-31T00:00:00.000000000', '2014-06-30T00:00:00.000000000',\n", + " '2014-07-31T00:00:00.000000000'], dtype='datetime64[ns]')
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,\n", + " 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,\n", + " 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,\n", + " 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,\n", + " 72, 73, 74, 75, 76, 77], dtype=int32)
array([[35.57694817, 35.57794762, 35.60395241, ..., 34.91988182,\n", + " 34.91988182, 34.92088127],\n", + " [35.45191574, 35.45191574, 35.45091629, ..., 34.92390442,\n", + " 34.92190552, 34.92290497],\n", + " [35.51696014, 35.51796341, 35.51696014, ..., 34.92086029,\n", + " 34.92086029, 34.91938019],\n", + " ...,\n", + " [35.90583801, 35.90784073, 35.90884018, ..., 34.95571518,\n", + " 34.94576263, 34.93976212],\n", + " [35.83697128, 35.83697128, 35.83496857, ..., 34.93985748,\n", + " 34.93885803, 34.93681145],\n", + " [35.77793884, 35.66791153, 35.66591263, ..., 34.94081116,\n", + " 34.93680954, 34.93680954]])
array([[18.7055006 , 18.57049942, 18.37200069, ..., 3.80349994,\n", + " 3.77100003, 3.745 ],\n", + " [19.14699936, 19.14599991, 19.14299965, ..., 3.75 ,\n", + " 3.71799994, 3.69199991],\n", + " [20.14500046, 20.14800072, 20.14500046, ..., 3.65700006,\n", + " 3.61199999, 3.56850004],\n", + " ...,\n", + " [15.48600006, 15.48799992, 15.48799992, ..., 3.91199994,\n", + " 3.81500006, 3.75399995],\n", + " [16.84399986, 16.84600067, 16.84600067, ..., 3.8440001 ,\n", + " 3.79399991, 3.74299991],\n", + " [20.13100052, 19.45999908, 19.26199913, ..., 3.81100011,\n", + " 3.75399995, 3.71399999]])
array([[ 6.45000005, 10.64999962, 14.89999962, ..., 1900.30004883,\n", + " 1949.95001221, 1998.65002441],\n", + " [ 5.9000001 , 10.39999962, 15.5 , ..., 1900.19995117,\n", + " 1949.89990234, 1997.30004883],\n", + " [ 6.29999971, 10.5 , 15.59999943, ..., 1899.89990234,\n", + " 1949.90002441, 1998.65002441],\n", + " ...,\n", + " [ 6.5 , 10.89999962, 15.80000019, ..., 1900.69995117,\n", + " 1950.09997559, 1998.5 ],\n", + " [ 6.80000019, 10.69999981, 15.90000057, ..., 1900.5 ,\n", + " 1950. , 2000.35003662],\n", + " [ 7.0999999 , 10.79999924, 15.89999962, ..., 1899.80004883,\n", + " 1950.20007324, 2000.40002441]])
<xarray.Dataset>\n", + "Dimensions: (date: 75, level: 78)\n", + "Coordinates:\n", + " * level (level) int32 0 1 2 3 4 5 6 7 8 ... 69 70 71 72 73 74 75 76 77\n", + " * date (date) datetime64[ns] 2012-07-13T22:33:06.019200 ... 2014-07...\n", + " lon (date) float64 -39.13 -37.28 -36.9 ... -33.83 -34.11 -34.38\n", + " lat (date) float64 47.19 46.72 46.45 46.23 ... 42.6 42.46 42.38\n", + "Data variables:\n", + " salinity (level, date) float64 nan nan nan nan nan ... nan nan nan nan\n", + " temperature (level, date) float64 nan nan nan nan nan ... nan nan nan nan\n", + " pressure (level, date) float64 nan nan nan nan nan ... nan nan nan nan
array([ 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17,\n", + " 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,\n", + " 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,\n", + " 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71,\n", + " 72, 73, 74, 75, 76, 77], dtype=int32)
array(['2012-07-13T22:33:06.019200000', '2012-07-23T22:54:59.990400000',\n", + " '2012-08-02T22:55:52.003200000', '2012-08-12T23:08:59.971200000',\n", + " '2012-08-22T23:29:01.968000000', '2012-09-01T23:17:38.976000000',\n", + " '2012-09-12T02:59:18.960000000', '2012-09-21T23:18:37.036800000',\n", + " '2012-10-02T03:00:17.971200000', '2012-10-11T23:13:27.984000000',\n", + " '2012-10-22T02:50:32.006399999', '2012-10-31T23:36:39.974400000',\n", + " '2012-11-11T02:40:46.041600000', '2012-11-20T23:08:29.990400000',\n", + " '2012-12-01T02:47:51.993600000', '2012-12-10T23:23:16.972800000',\n", + " '2012-12-21T02:58:48.979200000', '2012-12-30T23:07:23.030400000',\n", + " '2013-01-10T02:56:43.008000000', '2013-01-19T23:24:26.956800000',\n", + " '2013-01-30T02:43:53.011199998', '2013-02-08T23:15:27.043200000',\n", + " '2013-02-19T01:12:50.976000000', '2013-02-28T23:07:13.008000000',\n", + " '2013-03-11T02:43:30.979200000', '2013-03-20T23:17:22.992000000',\n", + " '2013-03-31T01:50:38.025600000', '2013-04-09T23:19:07.968000000',\n", + " '2013-04-20T02:53:29.990400000', '2013-04-29T23:28:33.024000000',\n", + " '2013-05-10T02:50:18.009600000', '2013-05-19T23:21:05.990400000',\n", + " '2013-05-30T02:50:30.969600000', '2013-06-08T23:32:49.027200000',\n", + " '2013-06-19T03:42:51.004800000', '2013-06-28T23:32:16.022400000',\n", + " '2013-07-09T03:28:30.979200000', '2013-07-18T23:33:57.974400000',\n", + " '2013-07-29T03:15:42.019200000', '2013-08-07T23:25:02.035200004',\n", + " '2013-08-18T01:47:44.966400000', '2013-08-28T03:02:59.020799996',\n", + " '2013-09-07T03:03:51.984000000', '2013-09-16T23:32:07.987200000',\n", + " '2013-09-27T04:08:27.974400000', '2013-10-06T23:25:39.964800000',\n", + " '2013-10-17T02:55:50.995200000', '2013-10-27T03:45:47.001600000',\n", + " '2013-11-06T01:14:52.022400000', '2013-11-16T03:29:54.009600000',\n", + " '2013-11-26T03:03:56.995200000', '2013-12-05T23:33:59.011200000',\n", + " '2013-12-16T02:58:01.977600000', '2013-12-25T23:22:43.017600000',\n", + " '2014-01-05T02:52:06.009600000', '2014-01-14T23:41:18.009600000',\n", + " '2014-01-25T03:00:43.977600000', '2014-02-03T23:29:13.977600000',\n", + " '2014-02-14T02:50:11.961600000', '2014-02-23T23:03:21.974400000',\n", + " '2014-03-06T02:58:03.964800000', '2014-03-15T23:10:28.012800000',\n", + " '2014-03-26T02:51:22.032000000', '2014-04-04T23:25:58.972800000',\n", + " '2014-04-15T03:00:45.964800000', '2014-04-24T23:24:40.003200000',\n", + " '2014-05-05T02:56:22.012800000', '2014-05-15T00:10:06.009600000',\n", + " '2014-05-25T03:02:43.036800000', '2014-06-03T23:34:53.961600000',\n", + " '2014-06-14T03:01:23.980800000', '2014-06-23T23:24:31.968000000',\n", + " '2014-07-04T03:08:30.019200000', '2014-07-13T23:47:43.008000000',\n", + " '2014-07-24T03:02:33.014400000'], dtype='datetime64[ns]')
array([-39.13 , -37.282, -36.9 , -36.89 , -37.053, -36.658, -35.963,\n", + " -35.184, -34.462, -33.784, -32.972, -32.546, -32.428, -32.292,\n", + " -32.169, -31.998, -31.824, -31.624, -31.433, -31.312, -31.107,\n", + " -31.147, -31.044, -31.14 , -31.417, -31.882, -32.145, -32.487,\n", + " -32.537, -32.334, -32.042, -31.892, -31.861, -31.991, -31.883,\n", + " -31.89 , -31.941, -31.889, -31.724, -31.412, -31.786, -31.561,\n", + " -31.732, -31.553, -31.862, -32.389, -32.318, -32.19 , -32.224,\n", + " -32.368, -32.306, -32.305, -32.65 , -33.093, -33.263, -33.199,\n", + " -33.27 , -33.237, -33.221, -33.011, -32.844, -32.981, -32.784,\n", + " -32.607, -32.87 , -33.196, -33.524, -33.956, -33.944, -33.71 ,\n", + " -33.621, -33.552, -33.828, -34.11 , -34.38 ])
array([47.187, 46.716, 46.45 , 46.23 , 45.459, 44.833, 44.452, 44.839,\n", + " 44.956, 44.676, 44.13 , 43.644, 43.067, 42.662, 42.513, 42.454,\n", + " 42.396, 42.256, 42.089, 41.944, 41.712, 41.571, 41.596, 41.581,\n", + " 41.351, 41.032, 40.912, 40.792, 40.495, 40.383, 40.478, 40.672,\n", + " 41.032, 40.864, 40.651, 40.425, 40.228, 40.197, 40.483, 40.311,\n", + " 40.457, 40.463, 40.164, 40.047, 39.963, 40.122, 40.57 , 40.476,\n", + " 40.527, 40.589, 40.749, 40.993, 41.162, 41.237, 41.448, 41.65 ,\n", + " 42.053, 42.311, 42.096, 41.683, 41.661, 41.676, 42.018, 42.395,\n", + " 42.532, 42.558, 42.504, 42.63 , 42.934, 42.952, 42.777, 42.722,\n", + " 42.601, 42.457, 42.379])
array([[nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " ...,\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan]])
array([[nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " ...,\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan]])
array([[nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " ...,\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan],\n", + " [nan, nan, nan, ..., nan, nan, nan]])
<xarray.DataArray (band: 3, y: 317, x: 625)>\n", + "[594375 values with dtype=uint16]\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
[594375 values with dtype=uint16]
array([1, 2, 3])
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])
<xarray.DataArray (band: 3, y: 317, x: 625)>\n", + "array([[[12812, 14483, 28688, ..., 16941, 17039, 14156],\n", + " [16744, 12435, 14008, ..., 11993, 10764, 10731],\n", + " [16564, 13107, 12615, ..., 11272, 12009, 11551],\n", + " ...,\n", + " [11911, 12157, 11157, ..., 23085, 16515, 11190],\n", + " [11518, 11387, 11010, ..., 12026, 14582, 9814],\n", + " [11567, 12255, 11862, ..., 9453, 11075, 13599]],\n", + "\n", + " [[15663, 16777, 27115, ..., 17826, 18186, 16482],\n", + " [16040, 15188, 16433, ..., 15663, 14221, 13779],\n", + " [16400, 15548, 15221, ..., 14631, 15778, 15974],\n", + " ...,\n", + " [15319, 15450, 15188, ..., 21610, 18841, 15810],\n", + " [15401, 15532, 15401, ..., 14516, 16859, 15859],\n", + " [15696, 15483, 15352, ..., 13844, 15810, 17334]],\n", + "\n", + " [[17776, 18645, 28475, ..., 19153, 19497, 18235],\n", + " [18563, 17842, 18596, ..., 17613, 17023, 17154],\n", + " [18628, 17563, 18088, ..., 17400, 18137, 17776],\n", + " ...,\n", + " [17645, 17744, 17432, ..., 22380, 20152, 17727],\n", + " [17580, 17383, 17252, ..., 17776, 18284, 17350],\n", + " [17891, 17678, 17662, ..., 16335, 17301, 18415]]], dtype=uint16)\n", + "Coordinates:\n", + " * band (band) int64 1 2 3\n", + " * y (y) float64 6.639e+06 6.639e+06 6.639e+06 ... 6.621e+06 6.621e+06\n", + " * x (x) float64 3.925e+05 3.926e+05 3.926e+05 ... 4.282e+05 4.283e+05\n", + "Attributes:\n", + " transform: (57.32777121388223, 0.0, 392465.9217302375, 0.0,...\n", + " crs: +init=epsg:3857\n", + " res: (57.32777121388223, 57.44833434997156)\n", + " is_tiled: 0\n", + " nodatavals: (nan, nan, nan)\n", + " scales: (1.0, 1.0, 1.0)\n", + " offsets: (0.0, 0.0, 0.0)\n", + " AREA_OR_POINT: Area\n", + " TIFFTAG_RESOLUTIONUNIT: 1 (unitless)\n", + " TIFFTAG_XRESOLUTION: 1\n", + " TIFFTAG_YRESOLUTION: 1
array([[[12812, 14483, 28688, ..., 16941, 17039, 14156],\n", + " [16744, 12435, 14008, ..., 11993, 10764, 10731],\n", + " [16564, 13107, 12615, ..., 11272, 12009, 11551],\n", + " ...,\n", + " [11911, 12157, 11157, ..., 23085, 16515, 11190],\n", + " [11518, 11387, 11010, ..., 12026, 14582, 9814],\n", + " [11567, 12255, 11862, ..., 9453, 11075, 13599]],\n", + "\n", + " [[15663, 16777, 27115, ..., 17826, 18186, 16482],\n", + " [16040, 15188, 16433, ..., 15663, 14221, 13779],\n", + " [16400, 15548, 15221, ..., 14631, 15778, 15974],\n", + " ...,\n", + " [15319, 15450, 15188, ..., 21610, 18841, 15810],\n", + " [15401, 15532, 15401, ..., 14516, 16859, 15859],\n", + " [15696, 15483, 15352, ..., 13844, 15810, 17334]],\n", + "\n", + " [[17776, 18645, 28475, ..., 19153, 19497, 18235],\n", + " [18563, 17842, 18596, ..., 17613, 17023, 17154],\n", + " [18628, 17563, 18088, ..., 17400, 18137, 17776],\n", + " ...,\n", + " [17645, 17744, 17432, ..., 22380, 20152, 17727],\n", + " [17580, 17383, 17252, ..., 17776, 18284, 17350],\n", + " [17891, 17678, 17662, ..., 16335, 17301, 18415]]], dtype=uint16)
array([1, 2, 3])
array([6639463.535997, 6639406.087663, 6639348.639328, ..., 6621424.759011,\n", + " 6621367.310677, 6621309.862342])
array([392494.585616, 392551.913387, 392609.241158, ..., 428152.459311,\n", + " 428209.787082, 428267.114853])