rasterix 0.1a1__tar.gz

rasterix-0.1a1/.pre-commit-config.yaml

@@ -0,0 +1,69 @@
+ci:
+  autoupdate_schedule: quarterly
+
+repos:
+  - repo: https://github.com/astral-sh/ruff-pre-commit
+    # Ruff version.
+    rev: "v0.11.4"
+    hooks:
+      - id: ruff
+        args: ["--fix", "--show-fixes"]
+      - id: ruff-format
+
+  - repo: https://github.com/pre-commit/mirrors-prettier
+    rev: "v4.0.0-alpha.8"
+    hooks:
+      - id: prettier
+
+  - repo: https://github.com/pre-commit/pre-commit-hooks
+    rev: v5.0.0
+    hooks:
+      - id: check-yaml
+      - id: trailing-whitespace
+      - id: end-of-file-fixer
+      - id: check-docstring-first
+
+  - repo: https://github.com/executablebooks/mdformat
+    rev: 0.7.22
+    hooks:
+      - id: mdformat
+        additional_dependencies:
+          - mdformat-black
+          - mdformat-myst
+
+  # - repo: https://github.com/kynan/nbstripout
+  #   rev: 0.8.1
+  #   hooks:
+  #     - id: nbstripout
+  #       args: [--extra-keys=metadata.kernelspec metadata.language_info.version]
+
+  - repo: https://github.com/codespell-project/codespell
+    rev: v2.4.1
+    hooks:
+      - id: codespell
+        additional_dependencies:
+          - tomli
+
+  - repo: https://github.com/abravalheri/validate-pyproject
+    rev: v0.24.1
+    hooks:
+      - id: validate-pyproject
+
+  - repo: https://github.com/rhysd/actionlint
+    rev: v1.7.7
+    hooks:
+      - id: actionlint
+        files: ".github/workflows/"
+        args:
+          [
+            "-ignore",
+            "SC1090",
+            "-ignore",
+            "SC2046",
+            "-ignore",
+            "SC2086",
+            "-ignore",
+            "SC2129",
+            "-ignore",
+            "SC2155",
+          ]

rasterix-0.1a1/LICENSE

@@ -0,0 +1,191 @@
+Apache License
+Version 2.0, January 2004
+http://www.apache.org/licenses/
+
+TERMS AND CONDITIONS FOR USE, REPRODUCTION, AND DISTRIBUTION
+
+1. Definitions.
+
+"License" shall mean the terms and conditions for use, reproduction, and
+distribution as defined by Sections 1 through 9 of this document.
+
+"Licensor" shall mean the copyright owner or entity authorized by the copyright
+owner that is granting the License.
+
+"Legal Entity" shall mean the union of the acting entity and all other entities
+that control, are controlled by, or are under common control with that entity.
+For the purposes of this definition, "control" means (i) the power, direct or
+indirect, to cause the direction or management of such entity, whether by
+contract or otherwise, or (ii) ownership of fifty percent (50%) or more of the
+outstanding shares, or (iii) beneficial ownership of such entity.
+
+"You" (or "Your") shall mean an individual or Legal Entity exercising
+permissions granted by this License.
+
+"Source" form shall mean the preferred form for making modifications, including
+but not limited to software source code, documentation source, and configuration
+files.
+
+"Object" form shall mean any form resulting from mechanical transformation or
+translation of a Source form, including but not limited to compiled object code,
+generated documentation, and conversions to other media types.
+
+"Work" shall mean the work of authorship, whether in Source or Object form, made
+available under the License, as indicated by a copyright notice that is included
+in or attached to the work (an example is provided in the Appendix below).
+
+"Derivative Works" shall mean any work, whether in Source or Object form, that
+is based on (or derived from) the Work and for which the editorial revisions,
+annotations, elaborations, or other modifications represent, as a whole, an
+original work of authorship. For the purposes of this License, Derivative Works
+shall not include works that remain separable from, or merely link (or bind by
+name) to the interfaces of, the Work and Derivative Works thereof.
+
+"Contribution" shall mean any work of authorship, including the original version
+of the Work and any modifications or additions to that Work or Derivative Works
+thereof, that is intentionally submitted to Licensor for inclusion in the Work
+by the copyright owner or by an individual or Legal Entity authorized to submit
+on behalf of the copyright owner. For the purposes of this definition,
+"submitted" means any form of electronic, verbal, or written communication sent
+to the Licensor or its representatives, including but not limited to
+communication on electronic mailing lists, source code control systems, and
+issue tracking systems that are managed by, or on behalf of, the Licensor for
+the purpose of discussing and improving the Work, but excluding communication
+that is conspicuously marked or otherwise designated in writing by the copyright
+owner as "Not a Contribution."
+
+"Contributor" shall mean Licensor and any individual or Legal Entity on behalf
+of whom a Contribution has been received by Licensor and subsequently
+incorporated within the Work.
+
+2. Grant of Copyright License.
+
+Subject to the terms and conditions of this License, each Contributor hereby
+grants to You a perpetual, worldwide, non-exclusive, no-charge, royalty-free,
+irrevocable copyright license to reproduce, prepare Derivative Works of,
+publicly display, publicly perform, sublicense, and distribute the Work and such
+Derivative Works in Source or Object form.
+
+3. Grant of Patent License.
+
+Subject to the terms and conditions of this License, each Contributor hereby
+grants to You a perpetual, worldwide, non-exclusive, no-charge, royalty-free,
+irrevocable (except as stated in this section) patent license to make, have
+made, use, offer to sell, sell, import, and otherwise transfer the Work, where
+such license applies only to those patent claims licensable by such Contributor
+that are necessarily infringed by their Contribution(s) alone or by combination
+of their Contribution(s) with the Work to which such Contribution(s) was
+submitted. If You institute patent litigation against any entity (including a
+cross-claim or counterclaim in a lawsuit) alleging that the Work or a
+Contribution incorporated within the Work constitutes direct or contributory
+patent infringement, then any patent licenses granted to You under this License
+for that Work shall terminate as of the date such litigation is filed.
+
+4. Redistribution.
+
+You may reproduce and distribute copies of the Work or Derivative Works thereof
+in any medium, with or without modifications, and in Source or Object form,
+provided that You meet the following conditions:
+
+You must give any other recipients of the Work or Derivative Works a copy of
+this License; and
+You must cause any modified files to carry prominent notices stating that You
+changed the files; and
+You must retain, in the Source form of any Derivative Works that You distribute,
+all copyright, patent, trademark, and attribution notices from the Source form
+of the Work, excluding those notices that do not pertain to any part of the
+Derivative Works; and
+If the Work includes a "NOTICE" text file as part of its distribution, then any
+Derivative Works that You distribute must include a readable copy of the
+attribution notices contained within such NOTICE file, excluding those notices
+that do not pertain to any part of the Derivative Works, in at least one of the
+following places: within a NOTICE text file distributed as part of the
+Derivative Works; within the Source form or documentation, if provided along
+with the Derivative Works; or, within a display generated by the Derivative
+Works, if and wherever such third-party notices normally appear. The contents of
+the NOTICE file are for informational purposes only and do not modify the
+License. You may add Your own attribution notices within Derivative Works that
+You distribute, alongside or as an addendum to the NOTICE text from the Work,
+provided that such additional attribution notices cannot be construed as
+modifying the License.
+You may add Your own copyright statement to Your modifications and may provide
+additional or different license terms and conditions for use, reproduction, or
+distribution of Your modifications, or for any such Derivative Works as a whole,
+provided Your use, reproduction, and distribution of the Work otherwise complies
+with the conditions stated in this License.
+
+5. Submission of Contributions.
+
+Unless You explicitly state otherwise, any Contribution intentionally submitted
+for inclusion in the Work by You to the Licensor shall be under the terms and
+conditions of this License, without any additional terms or conditions.
+Notwithstanding the above, nothing herein shall supersede or modify the terms of
+any separate license agreement you may have executed with Licensor regarding
+such Contributions.
+
+6. Trademarks.
+
+This License does not grant permission to use the trade names, trademarks,
+service marks, or product names of the Licensor, except as required for
+reasonable and customary use in describing the origin of the Work and
+reproducing the content of the NOTICE file.
+
+7. Disclaimer of Warranty.
+
+Unless required by applicable law or agreed to in writing, Licensor provides the
+Work (and each Contributor provides its Contributions) on an "AS IS" BASIS,
+WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied,
+including, without limitation, any warranties or conditions of TITLE,
+NON-INFRINGEMENT, MERCHANTABILITY, or FITNESS FOR A PARTICULAR PURPOSE. You are
+solely responsible for determining the appropriateness of using or
+redistributing the Work and assume any risks associated with Your exercise of
+permissions under this License.
+
+8. Limitation of Liability.
+
+In no event and under no legal theory, whether in tort (including negligence),
+contract, or otherwise, unless required by applicable law (such as deliberate
+and grossly negligent acts) or agreed to in writing, shall any Contributor be
+liable to You for damages, including any direct, indirect, special, incidental,
+or consequential damages of any character arising as a result of this License or
+out of the use or inability to use the Work (including but not limited to
+damages for loss of goodwill, work stoppage, computer failure or malfunction, or
+any and all other commercial damages or losses), even if such Contributor has
+been advised of the possibility of such damages.
+
+9. Accepting Warranty or Additional Liability.
+
+While redistributing the Work or Derivative Works thereof, You may choose to
+offer, and charge a fee for, acceptance of support, warranty, indemnity, or
+other liability obligations and/or rights consistent with this License. However,
+in accepting such obligations, You may act only on Your own behalf and on Your
+sole responsibility, not on behalf of any other Contributor, and only if You
+agree to indemnify, defend, and hold each Contributor harmless for any liability
+incurred by, or claims asserted against, such Contributor by reason of your
+accepting any such warranty or additional liability.
+
+END OF TERMS AND CONDITIONS
+
+APPENDIX: How to apply the Apache License to your work
+
+To apply the Apache License to your work, attach the following boilerplate
+notice, with the fields enclosed by brackets "[]" replaced with your own
+identifying information. (Don't include the brackets!) The text should be
+enclosed in the appropriate comment syntax for the file format. We also
+recommend that a file or class name and description of purpose be included on
+the same "printed page" as the copyright notice for easier identification within
+third-party archives.
+
+   Copyright 2014-2025 Deepak Cherian
+
+   Licensed under the Apache License, Version 2.0 (the "License");
+   you may not use this file except in compliance with the License.
+   You may obtain a copy of the License at
+
+     http://www.apache.org/licenses/LICENSE-2.0
+
+   Unless required by applicable law or agreed to in writing, software
+   distributed under the License is distributed on an "AS IS" BASIS,
+   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+   See the License for the specific language governing permissions and
+   limitations under the License.

rasterix-0.1a1/PKG-INFO

@@ -0,0 +1,26 @@
+Metadata-Version: 2.4
+Name: rasterix
+Version: 0.1a1
+Summary: Raster extensions for Xarray
+License-Expression: Apache-2.0
+Keywords: xarray
+Classifier: Development Status :: 4 - Beta
+Classifier: Natural Language :: English
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Programming Language :: Python :: 3.13
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: affine
+Requires-Dist: pandas>=2
+Requires-Dist: numpy>=2
+Requires-Dist: xarray>=2025
+Dynamic: license-file
+
+# rasterix: Raster tricks for Xarray
+
+<img src="rasterix.png" width="400">

rasterix-0.1a1/README.md

@@ -0,0 +1,3 @@
+# rasterix: Raster tricks for Xarray
+
+<img src="rasterix.png" width="400">

rasterix-0.1a1/design_notes/raster_index.md

@@ -0,0 +1,150 @@
+# Design for a RasterIndex
+
+## TL;DR
+
+1. We propose designing a RasterIndex that can handle many ways of expressing a raster → model space transformation.
+1. We propose that this RasterIndex _record_ the information for this transformation (e.g. `GeoTransform`) _internally_ and remove that information from the dataset when constructed.
+1. Since the information is recorded internally, a user must intentionally write the transformation back to the dataset, destroying the index while doing so, and then write to disk.
+
+## Goals:
+
+### UX goals:
+
+1. Make it easy to read a GeoTIFF with CRS and raster -> model space transformation information in to Xarray with appropriate indexes. There are at least two indexes: one associated with the CRS; and one with the transformation.
+1. The raster ↔ model transformation information can be ambiguous, so an explicit API should be provided.
+   1. [RPCs](http://geotiff.maptools.org/rpc_prop.html):
+      > The RPC model in a GeoTIFF file is supplementary to all other GeoTIFF tags and not directly related. That is, it is possible to have a conventional set of GeoTIFF tags (such as a tiepoint + pixel scale + projected coordinate system description) along with the RPCCoefficientTag. The RPCCoefficientTag is always describing a transformation to WGS84, regardless of what geographic coordinate system might be described in the coordinate system description tags of the GeoTIFF file. It is also possible to have only the RPCCoefficientTag tag and no other GeoTIFF tags.
+   1. [GeoTransform is not in the GeoTIFF standard](https://docs.ogc.org/is/19-008r4/19-008r4.html). Instead that uses ModelTiepointTag, ModelPixelScaleTag, ModelTransformationTag.
+   1. [GCPs are ambiguous](https://gdal.org/en/stable/user/raster_data_model.html#gcps):
+      > The GDAL data model does not imply a transformation mechanism that must be generated from the GCPs … this is left to the application. However 1st to 5th order polynomials are common.
+   1. [And for extra fun](https://gdal.org/en/stable/user/raster_data_model.html#gcps):
+      > Normally a dataset will contain either an affine geotransform, GCPs or neither. It is uncommon to have both, and it is undefined which is authoritative.
+
+### Index design goals:
+
+1. The CRS Index allows us to assert CRS compliance during alignment. This is provided by XProj.
+1. The transform index should allow us to:
+   1. Do accurate alignment in pixel space unaffected by floating-point inaccuracies in model-space;
+   1. All possible transforms have **offsets** which means they need to be kept up-to-date during slicing.
+   1. Allow extracting metadata necessary to accurately represent the information on disk.
+
+## Some complications
+
+### Handling projections
+
+1. There are at least 5 ways to record a raster ⇒ model space transformation.
+1. Information for these transforms may be stored in many places depending on the reading library:
+   1. `ds.spatial_ref.attrs`(rioxarray stores GeoTransform, gcps here)
+   1. `ds.band_data.attrs` (rioxarray stores TIEPOINTS here)
+   1. `ds.attrs` possibly since this would be the most direct way to map TIFF tags
+   1. It seems possible to store RPCs as either arrays or as attrs.
+
+We'd like a design that is extensible to handle all 5 (or more) cases; again suggesting an explicit lower-level API.
+
+### Composing with CRSIndex
+
+[One unanswered question so far is](https://github.com/benbovy/xproj/issues/22#issuecomment-2789459387)
+
+> I think the major decision is "who handles the spatial_ref / grid mapping variable". Should it be exclusively handled by xproj.CRSIndex? Or should it be bound to a geospatial index such as rasterix.RasterIndex, xvec.GeometryIndex, xoak.S2PointIndex, etc.?
+
+i.e. does the Index wrap CRSIndex too, or compose with it.
+Some points:
+
+1. reprojection requires handling both the transform and the CRS.
+1. the EDR-like selection API is similar.
+
+Importantly, GDAL chooses to write GeoTransform to the `grid_mapping` variable in netCDF which _also_ records CRS information.
+
+This gives us two options:
+
+1. RasterIndex wraps CRSIndex too and handles everything.
+1. RasterIndex extracts projection information, however it is stored (e.g. GeoTransform), and tracks it internally. Any functionality that relies on both the transformation and CRS will need to be built as an accessor layer.
+
+Below is a proposal for (2).
+
+## Proposal for transform index
+
+We design RasterIndex as a wrapper around **one** of many transform based indexes:
+
+1. AffineTransformIndex ↔ GeoTransform
+1. ModelTransformationIndex ↔ ModelTransformationTag
+1. ModelTiepointScaleIndex ↔ ModelTiepointTag + ModelPixelScaleTag
+1. GCPIndex ↔ Ground Control Points
+1. RPCIndex ↔ Rational Polynomial Coefficients
+1. Subsampled auxiliary coordinates, detailed in [CF section 8.3](https://cfconventions.org/Data/cf-conventions/cf-conventions-1.12/cf-conventions.html#compression-by-coordinate-subsampling) and equivalent to GDAL's [geolocation arrays](https://gdal.org/en/stable/development/rfc/rfc4_geolocate.html) with `PIXEL_STEP` and/or LINE_STEP\` > 1.
+
+Each of the wrapped index has an associated transform:
+
+```python
+@dataclass
+class RasterTransform:
+    rpc: RPC | None # rpcs
+    tiepoint_scale : ModelTiepointAndScale | None  # ModelTiepointTag, ModelPixelScaleTag
+    gcps: GroundControlPoints | None
+    transformation : ModelTransformation | None
+    geotransform : Affine | None # GeoTransform
+
+    def from_geotransform(attrs: dict) -> Self:
+        ...
+
+    def from_tiepoints(attrs: dict) -> Self:
+        ...
+
+    def from_gcps(gcps: ?) -> Self:
+        ...
+
+    def from_rpcs(?) -> Self:
+        ...
+
+    def from_geolocation_arrays(?) -> Self:
+        ...
+```
+
+### Read-time
+
+These transforms are constructed by **popping** the relevant information from a user-provided source.
+This is analogous to an "encode/decode" workflow we currently have in Xarray.
+
+```python
+transform = rasterix.RasterTransform.from_geotransform(ds.spatial_ref.attrs)
+# transform = rasterix.RasterTransform.from_tiepoints(ds.band_data.attrs)
+```
+
+By **popping** the information, the transformation is stored _only_ on the index, and must be rewritten back to the dataset by the user when writing to the disk. This is the RioXarray pattern.
+
+Once a transform is constructed, we could do
+
+```python
+    index = RasterIndex.from_transform(transform, dims_and_sizes=...)
+    ds = ds.assign_coords(xr.Coordinates.from_xindex(index))
+    ds = ds.set_xindex("spatial_ref", xproj.CRSIndex)
+```
+
+### Write-time
+
+Before write, we must write the transform (similar to rioxarray) and _destroy_ the RasterIndex instance.
+This seems almost required since there are many ways of recording the transformation information on the dataset; and many of these approaches might be used in the same dataset.
+
+```python
+    # encodes the internal RasterTransform in attributes of the `to` variable
+    # destroys the RasterIndex
+    ds = ds.rasterix.write_transform(to: Hashable | rasterix.SELF, formats=["geotransform", "tiepoint"])
+    ds.rio.to_raster()
+```
+
+Here:
+
+1. `SELF` could mean write to the object (Dataset | DataArray) attrs
+1. `formats` allow you to record the same information in multiple ways
+1. `.rio.write_transform` could just dispatch to this method.
+
+## Appendix
+
+### Encode/decode workflow for subsampled coordinates
+
+Taking the example 8.3 from [CF section 8.3](https://cfconventions.org/Data/cf-conventions/cf-conventions-1.12/cf-conventions.html#compression-by-coordinate-subsampling), the decode step may consist in:
+
+1. turn the tie point coordinate variables `lat(tp_yc, tp_xc)` and `lon(tp_yc, tp_xc)` into `lat(yc, xc)` and `lon(yc, xc)` in to Xarray coordinates associated with a custom transformation index that stores only the tie points. In other words, uncompress the dimensions of the `lat` & `lon` coordinates without uncompressing their data.
+1. also remove the tie point index variables and the interpolation variable, and track their data / metadata internally in the index
+
+The encode step would then consist in restoring the compressed tie point coordinate & index variables as well as the interpolation variable.

rasterix-0.1a1/notebooks/examples.ipynb

@@ -0,0 +1,191 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "id": "0",
+   "metadata": {},
+   "source": [
+    "# Raster dataset examples\n",
+    "\n",
+    "Written by ChatGPT.\n",
+    "\n",
+    "## GeoTiffs with GCPs and Tiepoints"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "1",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "import rasterio\n",
+    "from rasterio.control import GroundControlPoint\n",
+    "from rasterio.transform import Affine\n",
+    "\n",
+    "# Define file path for output GeoTIFF\n",
+    "output_tiff = \"output_with_gcps_tiepoints.tif\"\n",
+    "\n",
+    "# Create a dummy raster (100x100) with random values\n",
+    "width, height = 100, 100\n",
+    "data = np.random.randint(0, 255, (1, height, width), dtype=np.uint8)\n",
+    "\n",
+    "# Define Affine transformation (default identity)\n",
+    "transform = Affine.translation(0, 0) * Affine.scale(1, -1)\n",
+    "\n",
+    "# Define Ground Control Points (GCPs)\n",
+    "gcps = [\n",
+    "    GroundControlPoint(row=10, col=20, x=-122.123, y=37.456, z=50),\n",
+    "    GroundControlPoint(row=30, col=40, x=-122.456, y=37.789, z=60),\n",
+    "    GroundControlPoint(row=50, col=60, x=-122.789, y=38.123, z=55),\n",
+    "]\n",
+    "\n",
+    "# Define Tie Points (Model Tiepoints)\n",
+    "tie_points = [\n",
+    "    (0, 0, 0, -122.1, 37.4, 0),\n",
+    "    (50, 50, 0, -122.2, 37.5, 0),\n",
+    "]\n",
+    "\n",
+    "# Create a new GeoTIFF with GCPs\n",
+    "with rasterio.open(\n",
+    "    output_tiff,\n",
+    "    \"w\",\n",
+    "    driver=\"GTiff\",\n",
+    "    height=height,\n",
+    "    width=width,\n",
+    "    count=1,\n",
+    "    dtype=rasterio.uint8,\n",
+    "    crs=\"EPSG:4326\",  # Define coordinate reference system (WGS84)\n",
+    "    transform=transform,\n",
+    ") as dst:\n",
+    "    dst.write(data)\n",
+    "\n",
+    "    # Add Ground Control Points\n",
+    "    dst.gcps = (gcps, dst.crs)\n",
+    "\n",
+    "    # Add Tie Points as metadata\n",
+    "    dst.update_tags(TIEPOINTS=str(tie_points))\n",
+    "\n",
+    "print(f\"GeoTIFF with GCPs and Tie Points saved as {output_tiff}\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "2",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import xarray as xr\n",
+    "\n",
+    "xr.open_dataset(\"output_with_gcps_tiepoints.tif\", engine=\"rasterio\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "3",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import rasterio\n",
+    "\n",
+    "# Open the saved GeoTIFF\n",
+    "output_tiff = \"output_with_gcps_tiepoints.tif\"\n",
+    "\n",
+    "with rasterio.open(output_tiff) as dataset:\n",
+    "    print(\"GeoTransform (Affine Transformation):\")\n",
+    "    print(dataset.transform)\n",
+    "    print(dataset.gcps)\n",
+    "    print(dataset.tags())"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "4",
+   "metadata": {},
+   "source": [
+    "## GDAL netCDF with GeoTransform\n",
+    "\n",
+    "rasterio cannot write to netCDF"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "5",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import numpy as np\n",
+    "from osgeo import gdal, osr\n",
+    "\n",
+    "# Define the output netCDF file\n",
+    "output_file = \"output.nc\"\n",
+    "\n",
+    "# Define raster dimensions\n",
+    "cols = 10  # Number of columns\n",
+    "rows = 10  # Number of rows\n",
+    "bands = 1  # Number of bands\n",
+    "\n",
+    "# Define geotransform (Origin X, Pixel Width, Rotation, Origin Y, Rotation, Pixel Height)\n",
+    "geotransform = [100.0, 30.0, 0.0, 200.0, 0.0, -30.0]  # Example values\n",
+    "\n",
+    "# Define a spatial reference (WGS84 in this case)\n",
+    "srs = osr.SpatialReference()\n",
+    "srs.ImportFromEPSG(4326)  # WGS84\n",
+    "proj_wkt = srs.ExportToWkt()\n",
+    "\n",
+    "# Create a new netCDF file using GDAL\n",
+    "driver = gdal.GetDriverByName(\"netCDF\")\n",
+    "dataset = driver.Create(output_file, cols, rows, bands, gdal.GDT_Float32)\n",
+    "\n",
+    "# Set geotransform and projection\n",
+    "dataset.SetGeoTransform(geotransform)\n",
+    "dataset.SetProjection(proj_wkt)\n",
+    "\n",
+    "# Create some dummy data (a simple gradient)\n",
+    "data = np.arange(cols * rows, dtype=np.float32).reshape(rows, cols)\n",
+    "\n",
+    "# Write data to the first band\n",
+    "band = dataset.GetRasterBand(1)\n",
+    "band.WriteArray(data)\n",
+    "band.SetNoDataValue(-9999)  # Set a no-data value\n",
+    "\n",
+    "# Flush and close dataset\n",
+    "band.FlushCache()\n",
+    "dataset = None\n",
+    "\n",
+    "print(f\"NetCDF file '{output_file}' created with geotransform and projection.\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "6",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import xarray as xr\n",
+    "\n",
+    "xr.open_dataset(\"output.nc\")"
+   ]
+  }
+ ],
+ "metadata": {
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}