From dc6a32b44c41fb8da884777bcb7a944ec06f655c Mon Sep 17 00:00:00 2001
From: amsnyder <asnyder@usgs.gov>
Date: Mon, 8 Apr 2024 12:52:30 -0500
Subject: [PATCH] fix crs var detection
---
catalog/CooperMcKenzie/collection.json | 2 +-
catalog/CooperUpperDeschutes/collection.json | 2 +-
.../collection.json | 26 +--
.../collection.json | 26 +--
.../collection.json | 26 +--
...McKenzie_create_collection_from_zarr.ipynb | 168 ++++++++++++++----
...eschutes_create_collection_from_zarr.ipynb | 164 +++++++++++++----
...mulation_create_collection_from_zarr.ipynb | 59 ++++--
...analysis_create_collection_from_zarr.ipynb | 55 ++++--
...mulation_create_collection_from_zarr.ipynb | 55 ++++--
10 files changed, 425 insertions(+), 158 deletions(-)
diff --git a/catalog/CooperMcKenzie/collection.json b/catalog/CooperMcKenzie/collection.json
index 576c3808..d93187d5 100644
--- a/catalog/CooperMcKenzie/collection.json
+++ b/catalog/CooperMcKenzie/collection.json
@@ -204,7 +204,7 @@
},
"crs": {
"dimensions": [],
- "type": "data",
+ "type": "auxiliary",
"description": null,
"unit": null
},
diff --git a/catalog/CooperUpperDeschutes/collection.json b/catalog/CooperUpperDeschutes/collection.json
index 41d38046..235464c5 100644
--- a/catalog/CooperUpperDeschutes/collection.json
+++ b/catalog/CooperUpperDeschutes/collection.json
@@ -304,7 +304,7 @@
},
"crs": {
"dimensions": [],
- "type": "data",
+ "type": "auxiliary",
"description": null,
"unit": null
},
diff --git a/catalog/alaska_et_2020_ccsm4_historical_simulation/collection.json b/catalog/alaska_et_2020_ccsm4_historical_simulation/collection.json
index b6d76ce2..3e7c0288 100644
--- a/catalog/alaska_et_2020_ccsm4_historical_simulation/collection.json
+++ b/catalog/alaska_et_2020_ccsm4_historical_simulation/collection.json
@@ -37,7 +37,7 @@
2610000.0
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
},
"y": {
"type": "spatial",
@@ -48,10 +48,20 @@
-193582.29059999995
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
}
},
"cube:variables": {
+ "et0": {
+ "dimensions": [
+ "time",
+ "y",
+ "x"
+ ],
+ "type": "data",
+ "description": "Reference Evapotranspiration",
+ "unit": "mm"
+ },
"lat": {
"dimensions": [
"y",
@@ -70,19 +80,9 @@
"description": "longitude",
"unit": "degrees_east"
},
- "et0": {
- "dimensions": [
- "time",
- "y",
- "x"
- ],
- "type": "data",
- "description": "Reference Evapotranspiration",
- "unit": "mm"
- },
"polar_stereographic": {
"dimensions": [],
- "type": "data",
+ "type": "auxiliary",
"description": "CRS definition",
"unit": null
}
diff --git a/catalog/alaska_et_2020_era-interim_reanalysis/collection.json b/catalog/alaska_et_2020_era-interim_reanalysis/collection.json
index 367f5578..c0ad8f97 100644
--- a/catalog/alaska_et_2020_era-interim_reanalysis/collection.json
+++ b/catalog/alaska_et_2020_era-interim_reanalysis/collection.json
@@ -37,7 +37,7 @@
2610000.0
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
},
"y": {
"type": "spatial",
@@ -48,10 +48,20 @@
-193582.29059999995
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
}
},
"cube:variables": {
+ "et0": {
+ "dimensions": [
+ "time",
+ "y",
+ "x"
+ ],
+ "type": "data",
+ "description": "Reference Evapotranspiration",
+ "unit": "mm"
+ },
"lat": {
"dimensions": [
"y",
@@ -70,19 +80,9 @@
"description": "longitude",
"unit": "degrees_east"
},
- "et0": {
- "dimensions": [
- "time",
- "y",
- "x"
- ],
- "type": "data",
- "description": "Reference Evapotranspiration",
- "unit": "mm"
- },
"polar_stereographic": {
"dimensions": [],
- "type": "data",
+ "type": "auxiliary",
"description": "CRS definition",
"unit": null
}
diff --git a/catalog/alaska_et_2020_gfdl_historical_simulation/collection.json b/catalog/alaska_et_2020_gfdl_historical_simulation/collection.json
index a1c33b75..c2da734b 100644
--- a/catalog/alaska_et_2020_gfdl_historical_simulation/collection.json
+++ b/catalog/alaska_et_2020_gfdl_historical_simulation/collection.json
@@ -37,7 +37,7 @@
2610000.0
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
},
"y": {
"type": "spatial",
@@ -48,10 +48,20 @@
-193582.29059999995
],
"step": 20000.0,
- "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"Conversion\",\"name\":\"PROJ-based coordinate operation\",\"method\":{\"name\":\"PROJ-based operation method: +proj=stere +lat_0=90 +lat_ts=64 +lon_0=-152 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs\"},\"parameters\":[]}"
+ "reference_system": "{\"$schema\":\"https://proj.org/schemas/v0.5/projjson.schema.json\",\"type\":\"ProjectedCRS\",\"name\":\"unnamed\",\"base_crs\":{\"name\":\"WGS 84\",\"datum\":{\"type\":\"GeodeticReferenceFrame\",\"name\":\"World Geodetic System 1984\",\"ellipsoid\":{\"name\":\"WGS 84\",\"semi_major_axis\":6378137,\"inverse_flattening\":298.257223563}},\"coordinate_system\":{\"subtype\":\"ellipsoidal\",\"axis\":[{\"name\":\"Geodetic latitude\",\"abbreviation\":\"Lat\",\"direction\":\"north\",\"unit\":\"degree\"},{\"name\":\"Geodetic longitude\",\"abbreviation\":\"Lon\",\"direction\":\"east\",\"unit\":\"degree\"}]},\"id\":{\"authority\":\"EPSG\",\"code\":4326}},\"conversion\":{\"name\":\"unnamed\",\"method\":{\"name\":\"Polar Stereographic (variant B)\",\"id\":{\"authority\":\"EPSG\",\"code\":9829}},\"parameters\":[{\"name\":\"Latitude of standard parallel\",\"value\":64,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8832}},{\"name\":\"Longitude of origin\",\"value\":-152,\"unit\":\"degree\",\"id\":{\"authority\":\"EPSG\",\"code\":8833}},{\"name\":\"False easting\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8806}},{\"name\":\"False northing\",\"value\":0,\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1},\"id\":{\"authority\":\"EPSG\",\"code\":8807}}]},\"coordinate_system\":{\"subtype\":\"Cartesian\",\"axis\":[{\"name\":\"Easting\",\"abbreviation\":\"E\",\"direction\":\"south\",\"meridian\":{\"longitude\":90},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}},{\"name\":\"Northing\",\"abbreviation\":\"N\",\"direction\":\"south\",\"meridian\":{\"longitude\":180},\"unit\":{\"type\":\"LinearUnit\",\"name\":\"Meter\",\"conversion_factor\":1}}]}}"
}
},
"cube:variables": {
+ "et0": {
+ "dimensions": [
+ "time",
+ "y",
+ "x"
+ ],
+ "type": "data",
+ "description": "Reference Evapotranspiration",
+ "unit": "mm"
+ },
"lat": {
"dimensions": [
"y",
@@ -70,19 +80,9 @@
"description": "longitude",
"unit": "degrees_east"
},
- "et0": {
- "dimensions": [
- "time",
- "y",
- "x"
- ],
- "type": "data",
- "description": "Reference Evapotranspiration",
- "unit": "mm"
- },
"polar_stereographic": {
"dimensions": [],
- "type": "data",
+ "type": "auxiliary",
"description": "CRS definition",
"unit": null
}
diff --git a/workflows/archive/CooperMcKenzie_create_collection_from_zarr.ipynb b/workflows/archive/CooperMcKenzie_create_collection_from_zarr.ipynb
index 3faef1b1..82e72e2d 100644
--- a/workflows/archive/CooperMcKenzie_create_collection_from_zarr.ipynb
+++ b/workflows/archive/CooperMcKenzie_create_collection_from_zarr.ipynb
@@ -19,7 +19,9 @@
"cell_type": "code",
"execution_count": null,
"id": "201e0945-de55-45ff-b095-c2af009a4e62",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"import pystac\n",
@@ -132,7 +134,9 @@
"cell_type": "code",
"execution_count": null,
"id": "708f2cf5-79ab-49af-8067-de31d0d13ee6",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# open and view zarr dataset\n",
@@ -154,7 +158,9 @@
"cell_type": "code",
"execution_count": null,
"id": "482d204d-b5b6-40e5-ac42-55b459be1097",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# description of STAC collection\n",
@@ -188,7 +194,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ab91268f-7200-4cb1-979a-c7d75531d2c0",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"dims_auto_extract = ['X', 'Y', 'T']\n",
@@ -206,14 +214,34 @@
"## Get crs info"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "7d3fab2f-61ab-4ada-a2a6-c07982795f1b",
+ "metadata": {
+ "tags": []
+ },
+ "outputs": [],
+ "source": [
+ "crs_var = 'crs'"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
"id": "b03d52f3-1367-4255-a561-52ee4fc9e92d",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "crs = pyproj.CRS.from_cf(ds.crs.attrs)"
+ "# use pyproj to automatically extract crs info\n",
+ "crs = pyproj.CRS.from_cf(ds[crs_var].attrs)\n",
+ "\n",
+ "# alternatively, create the appropriate cartopy projection\n",
+ "# crs = ccrs.LambertConformal(central_longitude=crs_info.longitude_of_central_meridian, \n",
+ "# central_latitude=crs_info.latitude_of_projection_origin,\n",
+ "# standard_parallels=crs_info.standard_parallel)"
]
},
{
@@ -228,17 +256,21 @@
"cell_type": "code",
"execution_count": null,
"id": "38492a33-861e-46a3-8f0f-ceb7b2b5f42a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "ds.crs"
+ "ds[crs_var]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "4ba255b8-b76b-4f8d-8907-cc5589221e66",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs.to_proj4()"
@@ -264,7 +296,9 @@
"cell_type": "code",
"execution_count": null,
"id": "d46805e0-8e94-4ebe-aa01-d9a2d7051459",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# pull out lat/lon bbox for data\n",
@@ -283,7 +317,9 @@
"cell_type": "code",
"execution_count": null,
"id": "c80c650f-6dc6-4e2d-869c-0e674f528520",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"XX, YY = np.meshgrid(ds[dim_names_dict['X']].data, ds[dim_names_dict['Y']].data)"
@@ -293,7 +329,9 @@
"cell_type": "code",
"execution_count": null,
"id": "97fb5d51-ee5c-46d9-b256-6a5a51a42495",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"transformer = Transformer.from_crs(crs, \"EPSG:4326\", always_xy=True)\n",
@@ -304,7 +342,9 @@
"cell_type": "code",
"execution_count": null,
"id": "4ae8c4e5-9069-4972-bcf6-0bb824ea83a3",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(f'lower left coordinates (WGS84): {min(lon)}, {min(lat)}')\n",
@@ -315,7 +355,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6a054422-6b84-44bd-b102-ac293745cda8",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# create a spatial extent object \n",
@@ -334,7 +376,9 @@
"cell_type": "code",
"execution_count": null,
"id": "41a84995-867c-4152-8c57-85e3758bbb77",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# pull out first and last timestamps\n",
@@ -355,7 +399,9 @@
"cell_type": "code",
"execution_count": null,
"id": "1b1e37c4-5348-46ad-abc9-e005b5d6c02b",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"collection_extent = pystac.Extent(spatial=spatial_extent, temporal=temporal_extent)"
@@ -375,7 +421,9 @@
"cell_type": "code",
"execution_count": null,
"id": "adf6c59d-58cd-48b1-a5fd-3bb205a3ef56",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# define folder location where your STAC catalog json file is\n",
@@ -388,7 +436,9 @@
"cell_type": "code",
"execution_count": null,
"id": "7e96811b-95ae-406a-9728-55fc429d4e1f",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"if catalog.get_child(collection_id):\n",
@@ -417,7 +467,9 @@
"cell_type": "code",
"execution_count": null,
"id": "094832af-d22b-4359-b0f6-cf687acce5cc",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"asset_id = \"zarr-s3-osn\"\n",
@@ -459,7 +511,9 @@
"cell_type": "code",
"execution_count": null,
"id": "fc00946d-2880-491d-9b3b-3aeeb4414d6c",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# instantiate extention on collection\n",
@@ -478,7 +532,9 @@
"cell_type": "code",
"execution_count": null,
"id": "120a4914-3302-44a5-a282-0308ac84f040",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# list out dataset dimensions\n",
@@ -501,7 +557,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ea452f62-5644-49b6-8a4e-7dc4f649fd1a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs"
@@ -511,7 +569,9 @@
"cell_type": "code",
"execution_count": null,
"id": "1b1d05ff-8e43-44a7-8343-178b112c4ad6",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # the datacube extension can accept reference_system information as a numerical EPSG code, \n",
@@ -546,7 +606,9 @@
"cell_type": "code",
"execution_count": null,
"id": "13967ca9-1920-40f8-81cf-639d63439d71",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"time_step = pd.Timedelta(stac_helpers.get_step(ds, dim_names_dict['T'], time_dim=True)).isoformat()\n",
@@ -557,7 +619,9 @@
"cell_type": "code",
"execution_count": null,
"id": "8640779d-5131-4973-aad9-455d4acffd7a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for time steps:\n",
@@ -570,7 +634,9 @@
"cell_type": "code",
"execution_count": null,
"id": "de0efee3-af6e-4451-87a1-264503e24948",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for time steps, cont:\n",
@@ -594,7 +660,9 @@
"cell_type": "code",
"execution_count": null,
"id": "73306fb1-9b51-42d0-86fa-8d5e72644ee1",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"#x_step = stac_helpers.get_step(ds, dim_names_dict['X'])\n",
@@ -608,7 +676,9 @@
"cell_type": "code",
"execution_count": null,
"id": "2afe1dd0-ceb0-4a25-856e-38aefdbecd0c",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for spatial steps:\n",
@@ -623,7 +693,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ba997858-c048-4047-b03b-fc011a0d93b1",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for spatial steps, cont:\n",
@@ -645,7 +717,9 @@
"cell_type": "code",
"execution_count": null,
"id": "c046328f-4953-4df6-ac0d-2bab43f86888",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"y_step = stac_helpers.get_step(ds, dim_names_dict['Y'])\n",
@@ -659,7 +733,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6bd290eb-fdb3-4e2e-b9e8-4ebcfa3f608b",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for spatial steps:\n",
@@ -674,7 +750,9 @@
"cell_type": "code",
"execution_count": null,
"id": "e27871aa-d33f-4af2-9f4e-96a2e54f03ca",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for spatial steps, cont:\n",
@@ -696,7 +774,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6cfd212f-9309-44a8-845a-b2b6f536a598",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# get x, y bounds for extent of those dimensions (required)\n",
@@ -720,7 +800,9 @@
"cell_type": "code",
"execution_count": null,
"id": "acd45d3c-7845-47e6-9b7d-e35627a7ca9a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(dims)"
@@ -730,7 +812,9 @@
"cell_type": "code",
"execution_count": null,
"id": "5a443497-67a9-4dce-a8e9-b08d31a88223",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# create a dictionary of datacube dimensions you would like to assign to this dataset\n",
@@ -755,7 +839,9 @@
"cell_type": "code",
"execution_count": null,
"id": "e9272931-fc0b-4f2a-9546-283033e9cde8",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# drop metpy_crs coordinate we have added\n",
@@ -778,7 +864,7 @@
"vars_dict={}\n",
"for v in vars:\n",
" unit = stac_helpers.get_unit(ds, v)\n",
- " var_type = stac_helpers.get_var_type(ds, v)\n",
+ " var_type = stac_helpers.get_var_type(ds, v, crs_var)\n",
" long_name = stac_helpers.get_long_name(ds, v)\n",
" vars_dict[v] = pystac.extensions.datacube.Variable({'dimensions':list(ds[v].dims), 'type': var_type, 'description': long_name, 'unit': unit})"
]
@@ -795,7 +881,9 @@
"cell_type": "code",
"execution_count": null,
"id": "10141fd4-91d6-491d-878b-02653720891d",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# add dimesions and variables to collection extension\n",
@@ -828,7 +916,9 @@
"cell_type": "code",
"execution_count": null,
"id": "4b75791b-6b2d-40be-b7c6-330a60888fb5",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"if catalog.get_child(collection_id):\n",
diff --git a/workflows/archive/CooperUpperDeschutes_create_collection_from_zarr.ipynb b/workflows/archive/CooperUpperDeschutes_create_collection_from_zarr.ipynb
index a06fe0c7..2683d81d 100644
--- a/workflows/archive/CooperUpperDeschutes_create_collection_from_zarr.ipynb
+++ b/workflows/archive/CooperUpperDeschutes_create_collection_from_zarr.ipynb
@@ -19,7 +19,9 @@
"cell_type": "code",
"execution_count": null,
"id": "201e0945-de55-45ff-b095-c2af009a4e62",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"import pystac\n",
@@ -132,7 +134,9 @@
"cell_type": "code",
"execution_count": null,
"id": "708f2cf5-79ab-49af-8067-de31d0d13ee6",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# open and view zarr dataset\n",
@@ -154,7 +158,9 @@
"cell_type": "code",
"execution_count": null,
"id": "482d204d-b5b6-40e5-ac42-55b459be1097",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# description of STAC collection\n",
@@ -188,7 +194,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ab91268f-7200-4cb1-979a-c7d75531d2c0",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"dims_auto_extract = ['X', 'Y', 'T']\n",
@@ -206,14 +214,34 @@
"## Get crs info"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "931dc7c6-be81-47f2-b8d9-98f0e4ddfd64",
+ "metadata": {
+ "tags": []
+ },
+ "outputs": [],
+ "source": [
+ "crs_var = 'crs'"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
"id": "b03d52f3-1367-4255-a561-52ee4fc9e92d",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "crs = pyproj.CRS.from_cf(ds.crs.attrs)"
+ "# use pyproj to automatically extract crs info\n",
+ "crs = pyproj.CRS.from_cf(ds[crs_var].attrs)\n",
+ "\n",
+ "# alternatively, create the appropriate cartopy projection\n",
+ "# crs = ccrs.LambertConformal(central_longitude=crs_info.longitude_of_central_meridian, \n",
+ "# central_latitude=crs_info.latitude_of_projection_origin,\n",
+ "# standard_parallels=crs_info.standard_parallel)"
]
},
{
@@ -228,17 +256,21 @@
"cell_type": "code",
"execution_count": null,
"id": "38492a33-861e-46a3-8f0f-ceb7b2b5f42a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "ds.crs"
+ "ds[crs_var]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "4ba255b8-b76b-4f8d-8907-cc5589221e66",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs.to_proj4()"
@@ -264,7 +296,9 @@
"cell_type": "code",
"execution_count": null,
"id": "d46805e0-8e94-4ebe-aa01-d9a2d7051459",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# pull out lat/lon bbox for data\n",
@@ -283,7 +317,9 @@
"cell_type": "code",
"execution_count": null,
"id": "c80c650f-6dc6-4e2d-869c-0e674f528520",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"XX, YY = np.meshgrid(ds[dim_names_dict['X']].data, ds[dim_names_dict['Y']].data)"
@@ -293,7 +329,9 @@
"cell_type": "code",
"execution_count": null,
"id": "97fb5d51-ee5c-46d9-b256-6a5a51a42495",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"transformer = Transformer.from_crs(crs, \"EPSG:4326\", always_xy=True)\n",
@@ -304,7 +342,9 @@
"cell_type": "code",
"execution_count": null,
"id": "4ae8c4e5-9069-4972-bcf6-0bb824ea83a3",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(f'lower left coordinates (WGS84): {min(lon)}, {min(lat)}')\n",
@@ -315,7 +355,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6a054422-6b84-44bd-b102-ac293745cda8",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# create a spatial extent object \n",
@@ -334,7 +376,9 @@
"cell_type": "code",
"execution_count": null,
"id": "41a84995-867c-4152-8c57-85e3758bbb77",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# pull out first and last timestamps\n",
@@ -355,7 +399,9 @@
"cell_type": "code",
"execution_count": null,
"id": "1b1e37c4-5348-46ad-abc9-e005b5d6c02b",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"collection_extent = pystac.Extent(spatial=spatial_extent, temporal=temporal_extent)"
@@ -375,7 +421,9 @@
"cell_type": "code",
"execution_count": null,
"id": "adf6c59d-58cd-48b1-a5fd-3bb205a3ef56",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# define folder location where your STAC catalog json file is\n",
@@ -388,7 +436,9 @@
"cell_type": "code",
"execution_count": null,
"id": "7e96811b-95ae-406a-9728-55fc429d4e1f",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"if catalog.get_child(collection_id):\n",
@@ -417,7 +467,9 @@
"cell_type": "code",
"execution_count": null,
"id": "094832af-d22b-4359-b0f6-cf687acce5cc",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"asset_id = \"zarr-s3-osn\"\n",
@@ -459,7 +511,9 @@
"cell_type": "code",
"execution_count": null,
"id": "fc00946d-2880-491d-9b3b-3aeeb4414d6c",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# instantiate extention on collection\n",
@@ -487,7 +541,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ea452f62-5644-49b6-8a4e-7dc4f649fd1a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs"
@@ -497,7 +553,9 @@
"cell_type": "code",
"execution_count": null,
"id": "1b1d05ff-8e43-44a7-8343-178b112c4ad6",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # the datacube extension can accept reference_system information as a numerical EPSG code, \n",
@@ -532,7 +590,9 @@
"cell_type": "code",
"execution_count": null,
"id": "13967ca9-1920-40f8-81cf-639d63439d71",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"time_step = pd.Timedelta(stac_helpers.get_step(ds, dim_names_dict['T'], time_dim=True)).isoformat()\n",
@@ -543,7 +603,9 @@
"cell_type": "code",
"execution_count": null,
"id": "8640779d-5131-4973-aad9-455d4acffd7a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for spatial steps:\n",
@@ -558,7 +620,9 @@
"cell_type": "code",
"execution_count": null,
"id": "de0efee3-af6e-4451-87a1-264503e24948",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for time steps, cont:\n",
@@ -580,7 +644,9 @@
"cell_type": "code",
"execution_count": null,
"id": "73306fb1-9b51-42d0-86fa-8d5e72644ee1",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"x_step = stac_helpers.get_step(ds, dim_names_dict['X'])\n",
@@ -594,7 +660,9 @@
"cell_type": "code",
"execution_count": null,
"id": "2afe1dd0-ceb0-4a25-856e-38aefdbecd0c",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for spatial steps:\n",
@@ -609,7 +677,9 @@
"cell_type": "code",
"execution_count": null,
"id": "ba997858-c048-4047-b03b-fc011a0d93b1",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for spatial steps, cont:\n",
@@ -631,7 +701,9 @@
"cell_type": "code",
"execution_count": null,
"id": "c046328f-4953-4df6-ac0d-2bab43f86888",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"y_step = stac_helpers.get_step(ds, dim_names_dict['Y'])\n",
@@ -645,7 +717,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6bd290eb-fdb3-4e2e-b9e8-4ebcfa3f608b",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # optional debugging for spatial steps:\n",
@@ -660,7 +734,9 @@
"cell_type": "code",
"execution_count": null,
"id": "e27871aa-d33f-4af2-9f4e-96a2e54f03ca",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# # debugging for spatial steps, cont:\n",
@@ -682,7 +758,9 @@
"cell_type": "code",
"execution_count": null,
"id": "8b635e48-8d4c-4e29-ba8a-d20ee1bd4492",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# get x, y bounds for extent of those dimensions (required)\n",
@@ -706,7 +784,9 @@
"cell_type": "code",
"execution_count": null,
"id": "acd45d3c-7845-47e6-9b7d-e35627a7ca9a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"dims = list(ds.dims)\n",
@@ -717,7 +797,9 @@
"cell_type": "code",
"execution_count": null,
"id": "5a443497-67a9-4dce-a8e9-b08d31a88223",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# create a dictionary of datacube dimensions you would like to assign to this dataset\n",
@@ -742,7 +824,9 @@
"cell_type": "code",
"execution_count": null,
"id": "e9272931-fc0b-4f2a-9546-283033e9cde8",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# drop metpy_crs coordinate we have added\n",
@@ -765,7 +849,7 @@
"vars_dict={}\n",
"for v in vars:\n",
" unit = stac_helpers.get_unit(ds, v)\n",
- " var_type = stac_helpers.get_var_type(ds, v)\n",
+ " var_type = stac_helpers.get_var_type(ds, v, crs_var)\n",
" long_name = stac_helpers.get_long_name(ds, v)\n",
" vars_dict[v] = pystac.extensions.datacube.Variable({'dimensions':list(ds[v].dims), 'type': var_type, 'description': long_name, 'unit': unit})"
]
@@ -782,7 +866,9 @@
"cell_type": "code",
"execution_count": null,
"id": "10141fd4-91d6-491d-878b-02653720891d",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# add dimesions and variables to collection extension\n",
@@ -815,7 +901,9 @@
"cell_type": "code",
"execution_count": null,
"id": "4b75791b-6b2d-40be-b7c6-330a60888fb5",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"if catalog.get_child(collection_id):\n",
diff --git a/workflows/archive/alaska_et_2020_ccsm4_historical_simulation_create_collection_from_zarr.ipynb b/workflows/archive/alaska_et_2020_ccsm4_historical_simulation_create_collection_from_zarr.ipynb
index 187746e2..423be8d4 100644
--- a/workflows/archive/alaska_et_2020_ccsm4_historical_simulation_create_collection_from_zarr.ipynb
+++ b/workflows/archive/alaska_et_2020_ccsm4_historical_simulation_create_collection_from_zarr.ipynb
@@ -35,12 +35,12 @@
"import pandas as pd\n",
"import json\n",
"import numpy as np\n",
- "import metpy\n",
+ "import pyproj\n",
+ "from pyproj import Transformer\n",
"import cartopy.crs as ccrs\n",
"import cfunits\n",
"import json\n",
"import sys\n",
- "from pyproj import Transformer\n",
"sys.path.insert(1, '..')\n",
"import stac_helpers"
]
@@ -214,6 +214,18 @@
"## Get crs info"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "f168c59f-9326-42a7-8db3-28dc4b50a416",
+ "metadata": {
+ "tags": []
+ },
+ "outputs": [],
+ "source": [
+ "crs_var = 'polar_stereographic'"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
@@ -223,8 +235,13 @@
},
"outputs": [],
"source": [
- "ds = ds.metpy.parse_cf()\n",
- "crs = ds[list(ds.keys())[0]].metpy.pyproj_proj"
+ "# use pyproj to automatically extract crs info\n",
+ "crs = pyproj.CRS.from_cf(ds[crs_var].attrs)\n",
+ "\n",
+ "# alternatively, create the appropriate cartopy projection\n",
+ "# crs = ccrs.LambertConformal(central_longitude=crs_info.longitude_of_central_meridian, \n",
+ "# central_latitude=crs_info.latitude_of_projection_origin,\n",
+ "# standard_parallels=crs_info.standard_parallel)"
]
},
{
@@ -239,17 +256,21 @@
"cell_type": "code",
"execution_count": null,
"id": "03b05ed7-dfd7-497e-abb2-0711691502e3",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "ds.polar_stereographic"
+ "ds[crs_var]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "700f39b4-bb45-40d2-ac0b-dfbf813007e3",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs.to_proj4()"
@@ -275,7 +296,9 @@
"cell_type": "code",
"execution_count": null,
"id": "aff4370e-205c-4a35-9c14-1d74be8d1f11",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# pull out lat/lon bbox for data\n",
@@ -294,7 +317,9 @@
"cell_type": "code",
"execution_count": null,
"id": "1a958e87-8dc0-4a4d-985f-302fa09e13fb",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"XX, YY = np.meshgrid(ds[dim_names_dict['X']].data, ds[dim_names_dict['Y']].data)"
@@ -304,10 +329,12 @@
"cell_type": "code",
"execution_count": null,
"id": "d14bd8b0-d11c-4927-a6b1-7a8fcb35f5b8",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "transformer = Transformer.from_crs(crs.crs, \"EPSG:4326\", always_xy=True)\n",
+ "transformer = Transformer.from_crs(crs, \"EPSG:4326\", always_xy=True)\n",
"lon, lat = transformer.transform(XX.ravel(), YY.ravel())"
]
},
@@ -315,7 +342,9 @@
"cell_type": "code",
"execution_count": null,
"id": "5bd4dc6e-45a0-4af3-a29a-336b4f2d2532",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(f'lower left coordinates (WGS84): {min(lon)}, {min(lat)}')\n",
@@ -326,7 +355,9 @@
"cell_type": "code",
"execution_count": null,
"id": "dd752c5a-778a-4821-9478-55905d7ec121",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"# create a spatial extent object \n",
@@ -721,7 +752,7 @@
"vars_dict={}\n",
"for v in vars:\n",
" unit = stac_helpers.get_unit(ds, v)\n",
- " var_type = stac_helpers.get_var_type(ds, v)\n",
+ " var_type = stac_helpers.get_var_type(ds, v, crs_var)\n",
" long_name = stac_helpers.get_long_name(ds, v)\n",
" vars_dict[v] = pystac.extensions.datacube.Variable({'dimensions':list(ds[v].dims), 'type': var_type, 'description': long_name, 'unit': unit})"
]
diff --git a/workflows/archive/alaska_et_2020_era-interim_reanalysis_create_collection_from_zarr.ipynb b/workflows/archive/alaska_et_2020_era-interim_reanalysis_create_collection_from_zarr.ipynb
index 69ffaefd..c1920c82 100644
--- a/workflows/archive/alaska_et_2020_era-interim_reanalysis_create_collection_from_zarr.ipynb
+++ b/workflows/archive/alaska_et_2020_era-interim_reanalysis_create_collection_from_zarr.ipynb
@@ -35,12 +35,12 @@
"import pandas as pd\n",
"import json\n",
"import numpy as np\n",
- "import metpy\n",
+ "import pyproj\n",
+ "from pyproj import Transformer\n",
"import cartopy.crs as ccrs\n",
"import cfunits\n",
"import json\n",
"import sys\n",
- "from pyproj import Transformer\n",
"sys.path.insert(1, '..')\n",
"import stac_helpers"
]
@@ -214,6 +214,18 @@
"## Get crs info"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "2b97dba9-9af3-4979-a8ac-c2782985346e",
+ "metadata": {
+ "tags": []
+ },
+ "outputs": [],
+ "source": [
+ "crs_var = 'polar_stereographic'"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
@@ -223,8 +235,13 @@
},
"outputs": [],
"source": [
- "ds = ds.metpy.parse_cf()\n",
- "crs = ds[list(ds.keys())[0]].metpy.pyproj_proj"
+ "# use pyproj to automatically extract crs info\n",
+ "crs = pyproj.CRS.from_cf(ds[crs_var].attrs)\n",
+ "\n",
+ "# alternatively, create the appropriate cartopy projection\n",
+ "# crs = ccrs.LambertConformal(central_longitude=crs_info.longitude_of_central_meridian, \n",
+ "# central_latitude=crs_info.latitude_of_projection_origin,\n",
+ "# standard_parallels=crs_info.standard_parallel)"
]
},
{
@@ -239,17 +256,21 @@
"cell_type": "code",
"execution_count": null,
"id": "bec46024-befc-4161-a91a-e4a38daddeff",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "ds.polar_stereographic"
+ "ds[crs_var]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "eaaac451-7262-4473-aebb-ed7ff59f4b37",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs.to_proj4()"
@@ -292,7 +313,9 @@
"cell_type": "code",
"execution_count": null,
"id": "018385b2-121c-4d9c-ac3e-405e25db7bef",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"XX, YY = np.meshgrid(ds[dim_names_dict['X']].data, ds[dim_names_dict['Y']].data)"
@@ -302,10 +325,12 @@
"cell_type": "code",
"execution_count": null,
"id": "9cfdfda6-2bb3-4647-b922-e5bdf09cb611",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "transformer = Transformer.from_crs(crs.crs, \"EPSG:4326\", always_xy=True)\n",
+ "transformer = Transformer.from_crs(crs, \"EPSG:4326\", always_xy=True)\n",
"lon, lat = transformer.transform(XX.ravel(), YY.ravel())"
]
},
@@ -313,7 +338,9 @@
"cell_type": "code",
"execution_count": null,
"id": "f3305ce4-0143-4b6a-8b14-c8b1817a0803",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(f'lower left coordinates (WGS84): {min(lon)}, {min(lat)}')\n",
@@ -324,7 +351,9 @@
"cell_type": "code",
"execution_count": null,
"id": "87de6aeb-f6cd-40ab-a747-de757c9d0af1",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"spatial_extent = pystac.SpatialExtent(bboxes=[[min(lon).item(), min(lat).item(), max(lon).item(), max(lat).item()]])"
@@ -719,7 +748,7 @@
"vars_dict={}\n",
"for v in vars:\n",
" unit = stac_helpers.get_unit(ds, v)\n",
- " var_type = stac_helpers.get_var_type(ds, v)\n",
+ " var_type = stac_helpers.get_var_type(ds, v, crs_var)\n",
" long_name = stac_helpers.get_long_name(ds, v)\n",
" vars_dict[v] = pystac.extensions.datacube.Variable({'dimensions':list(ds[v].dims), 'type': var_type, 'description': long_name, 'unit': unit})"
]
diff --git a/workflows/archive/alaska_et_2020_gfdl_historical_simulation_create_collection_from_zarr.ipynb b/workflows/archive/alaska_et_2020_gfdl_historical_simulation_create_collection_from_zarr.ipynb
index 7a5470f6..f4239647 100644
--- a/workflows/archive/alaska_et_2020_gfdl_historical_simulation_create_collection_from_zarr.ipynb
+++ b/workflows/archive/alaska_et_2020_gfdl_historical_simulation_create_collection_from_zarr.ipynb
@@ -35,12 +35,12 @@
"import pandas as pd\n",
"import json\n",
"import numpy as np\n",
- "import metpy\n",
+ "import pyproj\n",
+ "from pyproj import Transformer\n",
"import cartopy.crs as ccrs\n",
"import cfunits\n",
"import json\n",
"import sys\n",
- "from pyproj import Transformer\n",
"sys.path.insert(1, '..')\n",
"import stac_helpers"
]
@@ -214,6 +214,18 @@
"## Get crs info"
]
},
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "d58ae3a4-f3cc-433b-93ff-6dd6d3627ff3",
+ "metadata": {
+ "tags": []
+ },
+ "outputs": [],
+ "source": [
+ "crs_var = 'polar_stereographic'"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
@@ -223,8 +235,13 @@
},
"outputs": [],
"source": [
- "ds = ds.metpy.parse_cf()\n",
- "crs = ds[list(ds.keys())[0]].metpy.pyproj_proj"
+ "# use pyproj to automatically extract crs info\n",
+ "crs = pyproj.CRS.from_cf(ds[crs_var].attrs)\n",
+ "\n",
+ "# alternatively, create the appropriate cartopy projection\n",
+ "# crs = ccrs.LambertConformal(central_longitude=crs_info.longitude_of_central_meridian, \n",
+ "# central_latitude=crs_info.latitude_of_projection_origin,\n",
+ "# standard_parallels=crs_info.standard_parallel)"
]
},
{
@@ -239,17 +256,21 @@
"cell_type": "code",
"execution_count": null,
"id": "276e86a1-fa2b-4ebf-93a1-2f2984ade96a",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "ds.polar_stereographic"
+ "ds[crs_var]"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "4823953c-3a8a-42f0-b325-98da43386882",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"crs.to_proj4()"
@@ -292,7 +313,9 @@
"cell_type": "code",
"execution_count": null,
"id": "885f7514-f91d-4e8d-961d-8d710ee22eec",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"XX, YY = np.meshgrid(ds[dim_names_dict['X']].data, ds[dim_names_dict['Y']].data)"
@@ -302,10 +325,12 @@
"cell_type": "code",
"execution_count": null,
"id": "339ccd0d-a613-47c0-8fa0-83c608ab4b72",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
- "transformer = Transformer.from_crs(crs.crs, \"EPSG:4326\", always_xy=True)\n",
+ "transformer = Transformer.from_crs(crs, \"EPSG:4326\", always_xy=True)\n",
"lon, lat = transformer.transform(XX.ravel(), YY.ravel())"
]
},
@@ -313,7 +338,9 @@
"cell_type": "code",
"execution_count": null,
"id": "7d28e137-b744-407d-b347-33b2f1772803",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"print(f'lower left coordinates (WGS84): {min(lon)}, {min(lat)}')\n",
@@ -324,7 +351,9 @@
"cell_type": "code",
"execution_count": null,
"id": "6e3a9999-baa7-4f02-9ae7-89ebc64e7662",
- "metadata": {},
+ "metadata": {
+ "tags": []
+ },
"outputs": [],
"source": [
"spatial_extent = pystac.SpatialExtent(bboxes=[[min(lon).item(), min(lat).item(), max(lon).item(), max(lat).item()]])"
@@ -719,7 +748,7 @@
"vars_dict={}\n",
"for v in vars:\n",
" unit = stac_helpers.get_unit(ds, v)\n",
- " var_type = stac_helpers.get_var_type(ds, v)\n",
+ " var_type = stac_helpers.get_var_type(ds, v, crs_var)\n",
" long_name = stac_helpers.get_long_name(ds, v)\n",
" vars_dict[v] = pystac.extensions.datacube.Variable({'dimensions':list(ds[v].dims), 'type': var_type, 'description': long_name, 'unit': unit})"
]
--
GitLab