more development notes, add xproj to tests

Author: scottyhq

README.md

@@ -41,3 +41,15 @@ pip install rasterix
    hatch env run --env test.py3.13 run-pytest  # Run the tests without coverage reports
    hatch env run --env test.py3.13 run-coverage-html   # Run the tests with an html coverage report
    ```
+
+#### Using hatch
+
+Start a shell
+```
+hatch env run --env test.py3.13 ipython
+```
+
+Expose kernel for notebooks in JupyterLab
+```
+hatch env run --env test.py3.13 "python -m ipykernel install --user --name=rasterix"
+```

pyproject.toml

@@ -31,6 +31,7 @@ dependencies = [
     "pandas>=2",
     "numpy>=2",
     "xarray>=2025",
+    "xproj",
 ]
 dynamic=["version"]
 [project.optional-dependencies]
@@ -42,13 +43,16 @@ rasterize = [
 ]
 exactextract = ["exactextract"]
 test = [
-  "geodatasets",
   "dask-geopandas",
+  "geodatasets",
+  "exactextract",
+  "ipykernel",
+  "ipython",
+  "netCDF4",
   "odc-geo",
   "rasterio",
   "rioxarray",
-  "exactextract",
-  "netCDF4"
+  "ruff",
 ]
 
 [tool.hatch]
@@ -82,8 +86,9 @@ run-coverage = "pytest --cov-config=pyproject.toml --cov=pkg --cov-report xml --
 run-coverage-html = "pytest --cov-config=pyproject.toml --cov=pkg --cov-report html --cov=src"
 run-pytest = "run-coverage --no-cov"
 run-verbose = "run-coverage --verbose"
-run-mypy = "mypy src"
 list-env = "pip list"
+install-kernel = "python -m ipykernel install --user --name=rasterix"
+lint = "ruff check"
 
 [tool.ruff.lint]
 # E402: module level import not at top of file

tests/test_raster_index.py

@@ -2,17 +2,36 @@
 import rioxarray  # noqa
 import xarray as xr
 from affine import Affine
+import pytest
+
+from xarray.structure.merge import MergeError
 
 from rasterix import RasterIndex, assign_index
+# TODO: hook up xproj to remove need for import?
+import xproj
+
+def _open_test_raster():
+    source = "/vsicurl/https://noaadata.apps.nsidc.org/NOAA/G02135/south/daily/geotiff/2024/01_Jan/S_20240101_concentration_v3.0.tif"
+    return xr.open_dataarray(source, engine="rasterio")
 
 
 def test_rectilinear():
-    source = "/vsicurl/https://noaadata.apps.nsidc.org/NOAA/G02135/south/daily/geotiff/2024/01_Jan/S_20240101_concentration_v3.0.tif"
-    da_no_raster_index = xr.open_dataarray(source, engine="rasterio")
+    da_no_raster_index = _open_test_raster()
     da_raster_index = assign_index(da_no_raster_index)
     assert da_raster_index.equals(da_no_raster_index)
 
 
+def test_different_crs_same_geotransform():
+    da = _open_test_raster()
+    da1 = assign_index(da)
+    da2 = da1.copy()
+    da1 = da1.proj.assign_crs(spatial_ref='EPSG:32610')
+    da2 = da2.proj.assign_crs(spatial_ref='EPSG:32611')
+    with pytest.raises(
+        MergeError, match="conflicting values/indexes on objects to be combined for coordinate"
+    ):
+        da1 + da2
+
 # TODO: parameterize over
 # 1. y points up;
 # 2. y points down