{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Kish Basin statistics" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import os\n", "\n", "import cartopy.crs as ccrs\n", "import contextily as cx\n", "import matplotlib.pyplot as plt\n", "\n", "import matplotlib.patches as mpatches\n", "import seaborn as sns\n", "from cartopy.mpl.ticker import LatitudeFormatter, LongitudeFormatter\n", "from matplotlib_scalebar.scalebar import ScaleBar\n", "\n", "from h2ss import data as rd\n", "from h2ss import functions as fns" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# basemap cache directory\n", "cx.set_cache_dir(os.path.join(\"data\", \"basemaps\"))" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "plt.rcParams[\"xtick.major.size\"] = 0\n", "plt.rcParams[\"ytick.major.size\"] = 0" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read data layers" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "ds, extent = rd.kish_basin_data_depth_adjusted(\n", " dat_path=os.path.join(\"data\", \"kish-basin\"),\n", " bathymetry_path=os.path.join(\"data\", \"bathymetry\"),\n", ")" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "ds = fns.net_to_gross(ds)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset> Size: 15MB\n", "Dimensions: (halite: 4, y: 237, x: 218)\n", "Coordinates:\n", " * y (y) float64 2kB 5.936e+06 5.936e+06 ... 5.889e+06 5.889e+06\n", " * x (x) float64 2kB 6.966e+05 6.968e+05 ... 7.398e+05 7.4e+05\n", " spatial_ref int64 8B 0\n", " * halite (halite) <U8 128B 'Fylde' 'Mythop' 'Preesall' 'Rossall'\n", " crs int64 8B 0\n", "Data variables:\n", " BaseDepth (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " Thickness (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " TopDepth (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " TopTWT (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " TopDepthSeabed (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " BaseDepthSeabed (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " Bathymetry (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " NetToGross (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan\n", " ThicknessNet (halite, y, x) float64 2MB nan nan nan nan ... nan nan nan
\n", " | BaseDepth | \n", "Thickness | \n", "TopDepth | \n", "TopTWT | \n", "TopDepthSeabed | \n", "BaseDepthSeabed | \n", "Bathymetry | \n", "NetToGross | \n", "ThicknessNet | \n", "
---|---|---|---|---|---|---|---|---|---|
count | \n", "63868.000000 | \n", "63456.000000 | \n", "63868.000000 | \n", "80386.000000 | \n", "63868.000000 | \n", "63868.000000 | \n", "63868.000000 | \n", "63456.000000 | \n", "63456.000000 | \n", "
mean | \n", "1011.301702 | \n", "112.835875 | \n", "899.392295 | \n", "653.242487 | \n", "852.308746 | \n", "964.218153 | \n", "-47.083550 | \n", "0.365550 | \n", "52.781874 | \n", "
std | \n", "536.478279 | \n", "113.380006 | \n", "507.401263 | \n", "259.464713 | \n", "514.491483 | \n", "543.935914 | \n", "24.235226 | \n", "0.102332 | \n", "68.099619 | \n", "
min | \n", "120.324500 | \n", "0.000000 | \n", "120.324500 | \n", "136.616000 | \n", "69.237007 | \n", "73.091809 | \n", "-125.472397 | \n", "0.261677 | \n", "0.000000 | \n", "
25% | \n", "574.370900 | \n", "4.344600 | \n", "476.546675 | \n", "468.859575 | \n", "423.502166 | \n", "517.393110 | \n", "-66.155502 | \n", "0.265696 | \n", "1.154345 | \n", "
50% | \n", "952.930250 | \n", "90.634850 | \n", "818.911100 | \n", "631.917300 | \n", "766.899198 | \n", "906.093512 | \n", "-41.403992 | \n", "0.345530 | \n", "31.317073 | \n", "
75% | \n", "1334.956375 | \n", "178.472675 | \n", "1207.349725 | \n", "776.623200 | \n", "1162.222693 | \n", "1298.014215 | \n", "-26.248158 | \n", "0.426796 | \n", "76.171349 | \n", "
max | \n", "3512.489700 | \n", "1111.067000 | \n", "3273.030800 | \n", "1911.048300 | \n", "3234.948776 | \n", "3475.165874 | \n", "-2.961994 | \n", "0.750000 | \n", "833.300250 | \n", "