windkit 2.0.0__tar.gz → 2.0.1__tar.gz

{windkit-2.0.0 → windkit-2.0.1}/CHANGELOG.md

@@ -2,14 +2,22 @@
 
 All major changes are listed here.
 
-## 2.0.0 (2025-07-21)
+## [2.0.1 (2026-01-29)](https://docs.wasp.dk/windkit/v2.0.1/)
 
 ### Highlights
 
-### Performance improvements
+#### Performance improvements
 * Significantly reduced import time by lazy loading some packages and data. Should be approx 1/3 of previous import time.
 
-### Sheer and veer functionality
+#### Documentation
+* Restructured user guide into focused, object-based documentation pages:
+  - Wind climates section with dedicated pages for TSWC, BWC, WWC, and GWC including I/O formats and usage examples
+  - Topography section with elevation, roughness, and landcover map documentation including format support tables
+  - Spatial section explaining point, stacked_point, raster, and cuboid structures with operation examples
+  - Wind farms section documenting wind turbine datasets and WTG objects
+  - New overview section with data structures, schemas, and tutorial data documentation
+
+#### Sheer and veer functionality
 
 * Added new function `wk.shear_extrapolate` to extrapolate wind speeds to different heights using a power law shear profile.
 * Added new function `wk.shear_exponent` to calculate a shear exponent from wind speeds
@@ -23,16 +31,24 @@ All major changes are listed here.
 
 * new function `create_geowc`, which creates a dummy geostrophic wind climate
 * functions `mean_wind_speed`, `mean_power_density`, `ws_freq_gt_mean`, `ws_cdf` and `mean_ws_moment` are now also working when a geostrophic wind climate is passed.
-* bug fix to interpolation, which would create invalid spatial structures when a dimension was in exclude dims but not present in the target structure (e.g. when doing interpolation from `stacked_point` to `point`).
 
 ### New features
 * Added new function `wk.extrapolate_to_height` to extrapolate wind climate objects to different heights. Currently, only works for time-series wind climates (TSWC).
-* Fixed bug where a point dataset with several horizontal locations would be converted to a single horizontal location when using the function `to_stacked_point`
 * `wk.roughness_map_to_file` now allows to write to a old fashioned `.map` file for the special case where these is a single polygon and the external roughness length is not specified (None). For this case we assume that the external roughness length is the same as the roughness length inside the polygon.
 * `wk.roughness_map_to_file` now allows users to write to a `.map` as long all displacement heights in the map are zero.
+* For `wv_count` objects (thus also geostrophic wind climates) if all of the bins have nan's, the result of function `mean_wind_speed`, `mean_power_density`, `ws_freq_gt_mean` and `mean_ws_moment` operating on such a function is also nan, whereas before it would be 0.
+* Added new function `wk.load_tutorial_data` to download and load tutorial data directly into memory as windkit-compatible objects (xarray.Dataset and geopandas.GeoDataFrame).
+* Changed `wk.get_tutorial_data` to return a `Path` to the extracted tutorial data folder, allowing users to load files manually using windkit readers or generic loaders. The function now includes detailed docstrings describing all available files for each dataset, including .nc, .omwc, .wtg, .map, and .csv.
+* Added `wk.read_wind_turbines` and `wk.wind_turbines_to_file` functions to read and write wind turbine datasets in various formats. Both only support CSV for now. The CSV should follow the windkit wind turbine format specification.
+
+
+### Bug fixes
+* Fixed pixel misalignment when downloading raster maps from Google Earth Engine. Downloads now use `crs_transform` instead of `scale` to preserve exact pixel alignment with source datasets (e.g., Globcover, SRTM, WorldCover). This fixes half-pixel offsets (~150m for Globcover) and ~12-15m offsets for elevation data.
 * Bugfix in `read_bwc`: it now opens `bwc` stored in netcdf correctly.
 but some of this functionality was only available in pywasp, so it has been moved there.
-* For `wv_count` objects (thus also geostrophic wind climates) if all of the bins have nan's, the result of function `mean_wind_speed`, `mean_power_density`, `ws_freq_gt_mean` and `mean_ws_moment` operating on such a function is also nan, whereas before it would be 0.
+* Fixed bug where a point dataset with several horizontal locations would be converted to a single horizontal location when using the function `to_stacked_point`
+* bug fix to interpolation, which would create invalid spatial structures when a dimension was in exclude dims but not present in the target structure (e.g. when doing interpolation from `stacked_point` to `point`).
+
 
 ### Changes
 * Changed default value of `use_bounds` to `True` in `wk.spatial.BBox.reproject` method.

{windkit-2.0.0 → windkit-2.0.1}/PKG-INFO

@@ -1,50 +1,72 @@
 Metadata-Version: 2.4
 Name: windkit
-Version: 2.0.0
+Version: 2.0.1
 Summary: WindKit provides core functionalities for working with wind resource data
-Home-page: https://docs.wasp.dk/windkit/latest
 Author: DTU Wind Energy
-License: BSD-3-Clause
+License-Expression: BSD-3-Clause
 Classifier: Programming Language :: Python :: 3
 Classifier: Operating System :: OS Independent
 Requires-Python: >=3.11
 Description-Content-Type: text/markdown
 License-File: LICENSE
+Requires-Dist: geopandas>=1.0
+Requires-Dist: h5netcdf
 Requires-Dist: lxml>=5.0
 Requires-Dist: netcdf4>=1.6
 Requires-Dist: numpy>=1.26
 Requires-Dist: packaging>=23
 Requires-Dist: pandas>=2.1
+Requires-Dist: platformdirs>=4.3
+Requires-Dist: pydantic>=2.7
+Requires-Dist: pydantic-settings>=2.3
 Requires-Dist: pyproj>=3.6
+Requires-Dist: pyyaml
+Requires-Dist: rasterio>=1.3
+Requires-Dist: requests
+Requires-Dist: rioxarray>=0.17
 Requires-Dist: scipy>=1.12
-Requires-Dist: shapely>=2
+Requires-Dist: shapely>=2.0
+Requires-Dist: toml
 Requires-Dist: xarray>=2023.9
-Requires-Dist: rioxarray>=0.17
-Requires-Dist: rasterio>=1.3
-Requires-Dist: geopandas>=1.0
-Requires-Dist: h5netcdf>=1.3
-Requires-Dist: h5py>=3.9
-Requires-Dist: matplotlib>=3.8
-Requires-Dist: pydantic-settings>=2.3
-Requires-Dist: pydantic>=2.7
-Requires-Dist: platformdirs>=4.3.8
-Requires-Dist: plotly>=5
-Requires-Dist: kaleido>=0.2.1
-Requires-Dist: requests>=2.31
-Requires-Dist: toml>=0.10
 Provides-Extra: plots
-Requires-Dist: matplotlib>=3.6; extra == "plots"
-Requires-Dist: plotly>=5.12; extra == "plots"
-Requires-Dist: kaleido>=0.2.1; extra == "plots"
+Requires-Dist: matplotlib>=3.8; extra == "plots"
+Requires-Dist: plotly>=5.12.0; extra == "plots"
+Requires-Dist: kaleido; extra == "plots"
+Provides-Extra: performance
+Requires-Dist: numba>=0.57.0; extra == "performance"
 Provides-Extra: data-apis
-Requires-Dist: ee; extra == "data-apis"
-Requires-Dist: pystac_client; extra == "data-apis"
-Requires-Dist: planetary_computer; extra == "data-apis"
-Requires-Dist: fsspec; extra == "data-apis"
 Requires-Dist: adlfs; extra == "data-apis"
-Requires-Dist: zarr; extra == "data-apis"
-Provides-Extra: performance
-Requires-Dist: numba; extra == "performance"
+Requires-Dist: earthengine-api; extra == "data-apis"
+Requires-Dist: fsspec; extra == "data-apis"
+Requires-Dist: planetary-computer; extra == "data-apis"
+Requires-Dist: pystac-client; extra == "data-apis"
+Requires-Dist: zarr>2.15; extra == "data-apis"
+Provides-Extra: dev
+Requires-Dist: dash; extra == "dev"
+Requires-Dist: ipython; extra == "dev"
+Requires-Dist: nbval; extra == "dev"
+Requires-Dist: pre-commit; extra == "dev"
+Requires-Dist: pytest; extra == "dev"
+Requires-Dist: pytest-cov; extra == "dev"
+Requires-Dist: pytest-xdist; extra == "dev"
+Requires-Dist: windIO; extra == "dev"
+Provides-Extra: doc
+Requires-Dist: docutils; extra == "doc"
+Requires-Dist: jupyter; extra == "doc"
+Requires-Dist: jupyter-dash; extra == "doc"
+Requires-Dist: jupyter-sphinx; extra == "doc"
+Requires-Dist: myst-parser; extra == "doc"
+Requires-Dist: nbsphinx; extra == "doc"
+Requires-Dist: nbsphinx-link; extra == "doc"
+Requires-Dist: pydata-sphinx-theme; extra == "doc"
+Requires-Dist: recommonmark; extra == "doc"
+Requires-Dist: sphinx; extra == "doc"
+Requires-Dist: sphinx-argparse; extra == "doc"
+Requires-Dist: sphinx-copybutton; extra == "doc"
+Requires-Dist: sphinx-design; extra == "doc"
+Requires-Dist: sphinx-gallery; extra == "doc"
+Requires-Dist: sphinxcontrib-bibtex; extra == "doc"
+Requires-Dist: sphinxext-opengraph; extra == "doc"
 Dynamic: license-file
 
 # WindKit

windkit-2.0.1/pyproject.toml

@@ -0,0 +1,113 @@
+[build-system]
+requires = ["setuptools>=45", "setuptools_scm[toml]>=6.2"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "windkit"
+description = "WindKit provides core functionalities for working with wind resource data"
+readme = "README.md"
+license = "BSD-3-Clause"
+authors = [
+    {name = "DTU Wind Energy"}
+]
+classifiers = [
+    "Programming Language :: Python :: 3",
+    "Operating System :: OS Independent",
+]
+dynamic = ["version"]
+requires-python = ">=3.11"
+dependencies = [
+    "geopandas>=1.0",
+    "h5netcdf",
+    "lxml>=5.0",
+    "netcdf4>=1.6",
+    "numpy>=1.26",
+    "packaging>=23",
+    "pandas>=2.1",
+    "platformdirs>=4.3",
+    "pydantic>=2.7",
+    "pydantic-settings>=2.3",
+    "pyproj>=3.6",
+    "pyyaml",
+    "rasterio>=1.3",
+    "requests",
+    "rioxarray>=0.17",
+    "scipy>=1.12",
+    "shapely>=2.0",
+    "toml",
+    "xarray>=2023.9",
+]
+
+[project.optional-dependencies]
+plots = [
+    "matplotlib>=3.8",
+    "plotly>=5.12.0",
+    "kaleido",
+]
+performance = [
+    "numba>=0.57.0",
+]
+data_apis = [
+    "adlfs",
+    "earthengine-api",
+    "fsspec",
+    "planetary-computer",
+    "pystac-client",
+    "zarr>2.15",
+]
+dev = [
+    "dash",
+    "ipython",
+    "nbval",
+    "pre-commit",
+    "pytest",
+    "pytest-cov",
+    "pytest-xdist",
+    "windIO",
+]
+doc = [
+    "docutils",
+    "jupyter",
+    "jupyter-dash",
+    "jupyter-sphinx",
+    "myst-parser",
+    "nbsphinx",
+    "nbsphinx-link",
+    "pydata-sphinx-theme",
+    "recommonmark",
+    "sphinx",
+    "sphinx-argparse",
+    "sphinx-copybutton",
+    "sphinx-design",
+    "sphinx-gallery",
+    "sphinxcontrib-bibtex",
+    "sphinxext-opengraph",
+]
+
+[project.scripts]
+split_combomap = "windkit.topography.vector_map:_split_combo"
+
+[tool.setuptools_scm]
+write_to = "windkit/_version.py"
+
+[tool.setuptools]
+include-package-data = true
+
+[tool.setuptools.packages.find]
+namespaces = true
+
+[tool.setuptools.package-data]
+"windkit.data" = ["*.json"]
+"windkit.data.landcover_tables" = ["*.json"]
+"windkit.data.plot" = ["*.json"]
+"windkit.data.loss_tables" = ["*.csv"]
+"windkit.data.uncertainty_tables" = ["*.csv"]
+
+[tool.pytest.ini_options]
+minversion = "6.0"
+norecursedirs = "dist build .tox"
+testpaths = ["test",]
+markers = [
+    "slow: marks tests as slow (deselect with '-m \"not slow\"')",
+    "performance: marks tests as performance tests (run with '-m performance')",
+]

windkit-2.0.1/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

{windkit-2.0.0 → windkit-2.0.1}/windkit/__init__.py

@@ -105,6 +105,8 @@ __all__ = [
     "create_wind_turbines_from_dataframe",
     "create_wind_turbines_from_arrays",
     "wind_turbines_to_geodataframe",
+    "read_wind_turbines",
+    "wind_turbines_to_file",
     # workspace
     "Workspace",
     # weng_workspace
@@ -142,6 +144,7 @@ __all__ = [
     "total_uncertainty_factor",
     # tutorial_data
     "get_tutorial_data",
+    "load_tutorial_data",
 ]
 
 from windkit._version import version as __version__
@@ -238,6 +241,8 @@ from windkit.wind_farm.wind_turbines import (
     create_wind_turbines_from_dataframe,
     create_wind_turbines_from_arrays,
     wind_turbines_to_geodataframe,
+    read_wind_turbines,
+    wind_turbines_to_file,
 )
 from windkit.workspace import Workspace
 from windkit.weng_workspace import WengWorkspace
@@ -301,4 +306,4 @@ from windkit.xarray_structures.direction import (
 from windkit.xarray_structures.wsbin import (
     create_wsbin_coords,
 )
-from windkit.tutorial_data import get_tutorial_data
+from windkit.tutorial_data import get_tutorial_data, load_tutorial_data

{windkit-2.0.0 → windkit-2.0.1}/windkit/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '2.0.0'
-__version_tuple__ = version_tuple = (2, 0, 0)
+__version__ = version = '2.0.1'
+__version_tuple__ = version_tuple = (2, 0, 1)
 
 __commit_id__ = commit_id = None

windkit-2.0.1/windkit/integrations/_earth_engine.py

@@ -0,0 +1,337 @@
+# (c) 2022 DTU Wind Energy
+"""
+Module that downloads elevation and roughness maps
+
+Currently driven by the google earth engine, which is a cloud interface service to access
+landcover and elevation data, among other things.
+
+Installation & Setup
+--------------------
+
+The earth engine is an optional install, you can install with conda::
+
+    conda install earthengine-api google-cloud-sdk
+
+After installation you will have to do a one-time authentication
+step from the command line: ``earthengine authenticate``
+
+This will open a browser where you will have to allow google to use
+you google account to retrieve data from the google servers. If you are on a machine
+without the ability to use a browser (such as an HPC
+cluster), you will have to use ``earthengine authenticate --quiet``, which requires you to
+to manually copy the authentication code into the terminal.
+
+In addition, you will have to `sign up
+<https://signup.earthengine.google.com/#!/>`_ for the google earth engine and give
+a reason why you want to use the program. Please pay particular attention to
+their terms of service.
+
+Automated Datasets
+------------------
+
+Currently, the databases that have been added are the Copernicus Global Land Cover
+(CGLS-LC100), Copernicus CORINE land Cover (CORINE), MODIS Global Land Cover MCD12Q1
+(MODIS), Globcover and WorldCover
+(https://developers.google.com/earth-engine/datasets/catalog/ESA_WorldCover_v100) as landcover databases and
+NASA SRTM Digital Elevation 30m(SRTM), ALOS DSM: Global 30 (ALOS) and NASADEM
+(https://developers.google.com/earth-engine/datasets/catalog/NASA_NASADEM_HGT_001?hl=en) as elevation databases.
+
+The landcover databases have standard conversion tables that are included in
+``windkit/data/landcover_tables``.
+
+Google Earth Engine provides lists of
+`elevation
+<https://developers.google.com/earth-engine/datasets/tags/elevation>`_ and
+`land cover
+<https://developers.google.com/earth-engine/datasets/tags/landcover>`_
+data sources, which provide additional details about the various datasources used
+in this library.
+"""
+
+import logging
+import math
+import tempfile
+import urllib.request
+import warnings
+from pathlib import Path
+
+import rioxarray
+
+from ..spatial import BBox, clip
+
+CRS_GEOGRAPHIC = "EPSG:4326"
+
+# Dataset configuration: maps datasource name to (url, band, is_collection)
+_DATASET_CONFIG = {
+    "CGLS-LC100": (
+        "COPERNICUS/Landcover/100m/Proba-V/Global/2015",
+        "discrete_classification",
+        False,
+    ),
+    "CORINE": ("COPERNICUS/CORINE/V20/100m/2018", "landcover", False),
+    "MODIS": ("MODIS/006/MCD12Q1/2018_01_01", "LC_Type1", False),
+    "SRTM": ("USGS/SRTMGL1_003", "elevation", False),
+    "NASADEM": ("NASA/NASADEM_HGT/001", "elevation", False),
+    "ALOS": ("JAXA/ALOS/AW3D30/V2_2", "AVE_DSM", False),
+    "Globcover": ("ESA/GLOBCOVER_L4_200901_200912_V2_3", "landcover", False),
+    "WorldCover": ("ESA/WorldCover/v100", "Map", True),
+}
+
+ELEVATION_DATASETS = {"SRTM", "ALOS", "NASADEM"}
+LIST_DATA_SOURCES = list(_DATASET_CONFIG.keys())
+
+logger = logging.getLogger(__name__)
+
+# Initialize Google Earth Engine
+try:
+    import ee as earth_engine
+
+    try:
+        credentials_path = Path.home() / "earth_engine.json"
+        if credentials_path.exists():
+            service_account = "windkit@windkit.iam.gserviceaccount.com"
+            credentials = earth_engine.ServiceAccountCredentials(
+                service_account, str(credentials_path)
+            )
+            earth_engine.Initialize(credentials)
+        else:
+            earth_engine.Initialize()
+    except Exception:
+        warnings.warn(
+            "Could not initialize Google Earth Engine. Run 'earth_engine.Authenticate()' and try again."
+        )
+
+except ImportError:
+    earth_engine = None
+
+
+def _bbox_to_ee_geometry(bbox):
+    return earth_engine.Geometry(bbox.__geo_interface__, bbox.crs.to_string())
+
+
+def _get_image(datasource):
+    """Get image from Google Earth Engine for the specified datasource.
+
+    Parameters
+    ----------
+    datasource : str
+        One of: 'CGLS-LC100', 'CORINE', 'MODIS', 'Globcover', 'WorldCover',
+        'SRTM', 'ALOS', 'NASADEM'
+
+    Returns
+    -------
+    image : ee.Image
+        Google Earth Engine image for the specified band
+    """
+    if earth_engine is None:
+        raise ValueError(
+            "ee (earthengine-api) is required to get maps from Google Earth Engine"
+        )
+
+    if datasource not in _DATASET_CONFIG:
+        valid_sources = ", ".join(LIST_DATA_SOURCES)
+        raise ValueError(f"Please specify a valid data source from {valid_sources}")
+
+    url, band, is_collection = _DATASET_CONFIG[datasource]
+
+    if is_collection:
+        dataset = earth_engine.ImageCollection(url).first()
+    else:
+        dataset = earth_engine.Image(url)
+
+    return dataset.select(band)
+
+
+def _get_ee_map(lat, lon, buffer_dist=20000, source="SRTM", vector=False):
+    """Extract map from a given lat, lon
+
+    Extract the smallest square which fits a cirle with radius buffer_dist
+    around the coordinates lat,lon.
+
+    Parameters
+    ----------
+    lat : float
+        Center latitude from which we extract a map
+    lon : float
+        Center longitude from which we extract a map
+    buffer_dist : int, optional
+        Distance in meters from the given (lat,lon) where a map is extracted, by default 20000
+    source : str {"CGLS-LC100", "CORINE", "MODIS", "Globcover", "WorldCover", "SRTM", "ALOS", "NASADEM"}, optional
+        Landcover or elevation datasource, by default "SRTM"
+    vector:
+        If true, return the map in vector format else return a raster map
+    """
+    if vector:
+        raise NotImplementedError("This feature is not yet available.")
+
+    bbox = BBox.utm_bbox_from_geographic_coordinate(lon, lat, buffer_dist)
+    ras = _get_raster_map_from_earth_engine(bbox, dataset=source)
+    return ras
+
+
+def _compute_aligned_origin(bbox_min, bbox_max, scale, origin):
+    """Compute a new origin aligned to the native grid but near the bbox.
+
+    This prevents Earth Engine from computing pixels from a distant origin
+    (e.g., -180 degrees) to the bbox by finding the nearest grid-aligned
+    point just outside the bbox.
+
+    Parameters
+    ----------
+    bbox_min : float
+        Minimum bound of bbox in this dimension (west or south)
+    bbox_max : float
+        Maximum bound of bbox in this dimension (east or north)
+    scale : float
+        Pixel scale (positive or negative)
+    origin : float
+        Native grid origin
+
+    Returns
+    -------
+    float
+        New origin aligned to native grid, positioned just outside the bbox
+    """
+    # Choose reference bound based on scale direction
+    # Positive scale: origin at min side, negative scale: origin at max side
+    reference_bound = bbox_min if scale > 0 else bbox_max
+    n_pixels = math.floor((reference_bound - origin) / scale)
+    return origin + n_pixels * scale
+
+
+def _create_aligned_transform(native_transform, bbox):
+    """Create an adjusted transform aligned to the native grid near the bbox.
+
+    Parameters
+    ----------
+    native_transform : list
+        Earth Engine transform: [xScale, xShear, xOrigin, yShear, yScale, yOrigin]
+    bbox : BBox
+        Bounding box for the download
+
+    Returns
+    -------
+    list or None
+        Adjusted transform, or None if transform cannot be computed
+    """
+    x_scale, x_shear, x_origin, y_shear, y_scale, y_origin = native_transform
+
+    if x_scale == 0 or y_scale == 0:
+        return None
+
+    west, south, east, north = bbox.bounds()
+
+    new_x_origin = _compute_aligned_origin(west, east, x_scale, x_origin)
+    new_y_origin = _compute_aligned_origin(south, north, y_scale, y_origin)
+
+    return [x_scale, x_shear, new_x_origin, y_shear, y_scale, new_y_origin]
+
+
+def _get_download_url(ee_image, download_params, dataset):
+    """Get download URL, falling back to scale if crs_transform fails.
+
+    Some datasets (e.g., NASADEM) may not respect the adjusted transform origin
+    and fail with grid dimension errors.
+
+    Parameters
+    ----------
+    ee_image : ee.Image
+        Earth Engine image
+    download_params : dict
+        Download parameters (may be modified in place)
+    dataset : str
+        Dataset name for logging
+
+    Returns
+    -------
+    str
+        Download URL
+    """
+    try:
+        return ee_image.getDownloadURL(download_params)
+    except Exception as e:
+        if "crs_transform" in download_params and "grid dimensions" in str(e).lower():
+            logger.debug(
+                f"crs_transform failed for {dataset}, falling back to scale: {e}"
+            )
+            del download_params["crs_transform"]
+            download_params["scale"] = ee_image.projection().nominalScale().getInfo()
+            return ee_image.getDownloadURL(download_params)
+        raise
+
+
+def _get_raster_map_from_earth_engine(bbox, dataset="NASADEM", band=None):
+    """
+    Get map from Google Earth Engine.
+
+    Parameters
+    ----------
+    bbox : windkit.spatial.BBox
+        Bounding box of the map to download.
+
+    dataset : str, optional
+        Dataset to retrieve, by default "NASADEM"
+
+    band : str, optional
+        Band to retrieve, by default None
+
+    Returns
+    -------
+    da : xarray.DataArray
+        DataArray with the map
+
+    Notes
+    -----
+    When the requested CRS matches the dataset's native CRS, the download uses
+    `crs_transform` to preserve exact pixel alignment with the source data.
+    When the CRS differs, `scale` is used which may result in slight pixel
+    misalignment (typically sub-pixel shifts).
+    """
+    if earth_engine is None:
+        raise ValueError(
+            "ee (earthengine-api) is required to get maps from Google Earth Engine"
+        )
+
+    if not isinstance(bbox, BBox):
+        raise ValueError("bbox must be a BBox object or a windkit.spatial.BBox object.")
+
+    ee_image = _get_image(dataset)
+    projection_info = ee_image.projection().getInfo()
+    native_crs = projection_info.get("crs")
+    native_transform = projection_info.get("transform")
+    requested_crs = f"EPSG:{bbox.crs.to_epsg()}"
+
+    download_params = {
+        "region": _bbox_to_ee_geometry(bbox.reproject(CRS_GEOGRAPHIC)),
+        "format": "GEO_TIFF",
+        "crs": requested_crs,
+    }
+
+    # Use crs_transform for precise alignment when CRS matches
+    crs_matches = native_crs == requested_crs and native_transform is not None
+    adjusted_transform = (
+        _create_aligned_transform(native_transform, bbox) if crs_matches else None
+    )
+
+    if adjusted_transform is not None:
+        download_params["crs_transform"] = adjusted_transform
+    else:
+        download_params["scale"] = ee_image.projection().nominalScale().getInfo()
+
+    download_url = _get_download_url(ee_image, download_params, dataset)
+
+    tmpfile = tempfile.NamedTemporaryFile(delete=False)
+    try:
+        logger.debug(f"Downloading file {tmpfile.name}")
+        urllib.request.urlretrieve(download_url, tmpfile.name)
+
+        da = rioxarray.open_rasterio(tmpfile.name)
+        da = da.rename({"spatial_ref": "crs", "x": "west_east", "y": "south_north"})
+        da = da.drop_vars("band", errors="i")
+        da.name = "elevation" if dataset in ELEVATION_DATASETS else "landcover"
+        da = clip(da, bbox)
+        da = da.sortby(["south_north", "west_east"])
+        return da
+    finally:
+        tmpfile.close()
+        Path(tmpfile.name).unlink()

{windkit-2.0.0 → windkit-2.0.1}/windkit/spatial/__init__.py

@@ -97,6 +97,7 @@ from ._struct import (
     _is_vertical,
 )
 from .decorators import _stack_then_unstack
+from ._conversions import _to_spatial_struct
 from .spatial import (
     add_projected_wrf_coordinates,
     equal_spatial_shape,