{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "b558712f-231e-47d3-a044-9064827132a3", "metadata": {}, "outputs": [], "source": [ "# command to auto-reload modules when they are edited (easier for testing and debugging)\n", "%load_ext autoreload\n", "%autoreload 2" ] }, { "cell_type": "code", "execution_count": 4, "id": "be010a0b", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from scipy.signal.windows import tukey\n", "import healpix_geo\n", "import sys\n", "sys.path.append(\"..\")\n", "from healpix_analyse.make_rectangle import make_healpix_rectangle_from_lonlat, intersect_with_nested_cell_ids\n", "from healpix_analyse.resample import resample_to_latlon_grid\n", "from healpix_analyse.ps import ps" ] }, { "cell_type": "markdown", "id": "de4d4713-b421-436e-9dad-6bcf35511bcf", "metadata": { "jp-MarkdownHeadingCollapsed": true }, "source": [ "## Data from tina notebook's S2" ] }, { "cell_type": "code", "execution_count": 10, "id": "738c1085-dcc1-4f61-a108-baf296dc6ba0", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset> Size: 168MB\n",
"Dimensions: (cells: 4194304)\n",
"Coordinates:\n",
" cell_ids (cells) uint64 34MB dask.array<chunksize=(1048576,), meta=np.ndarray>\n",
"Dimensions without coordinates: cells\n",
"Data variables:\n",
" b02 (cells) float64 34MB dask.array<chunksize=(1048576,), meta=np.ndarray>\n",
" b03 (cells) float64 34MB dask.array<chunksize=(1048576,), meta=np.ndarray>\n",
" b04 (cells) float64 34MB dask.array<chunksize=(1048576,), meta=np.ndarray>\n",
" b08 (cells) float64 34MB dask.array<chunksize=(1048576,), meta=np.ndarray><xarray.Dataset> Size: 73MB\n",
"Dimensions: (cells: 1835008)\n",
"Dimensions without coordinates: cells\n",
"Data variables:\n",
" b02 (cells) float64 15MB dask.array<chunksize=(262144,), meta=np.ndarray>\n",
" b03 (cells) float64 15MB dask.array<chunksize=(262144,), meta=np.ndarray>\n",
" b04 (cells) float64 15MB dask.array<chunksize=(262144,), meta=np.ndarray>\n",
" b08 (cells) float64 15MB dask.array<chunksize=(262144,), meta=np.ndarray>\n",
" cell_ids (cells) uint64 15MB dask.array<chunksize=(262144,), meta=np.ndarray>\n",
"Attributes:\n",
" zarr_conventions: [{'uuid': '7b255807-140c-42ca-97f6-7a1cfecdbc38', 'nam...\n",
" dggs: {'name': 'healpix', 'refinement_level': 20, 'indexing_...