windkit 0.6.3__tar.gz → 1.0.2__tar.gz

windkit-1.0.2/CHANGELOG.md

@@ -0,0 +1,700 @@
+# Changelog
+
+All major changes are listed here.
+
+## 1.0.2 (2025-07-21)
+
+WindKit 1.0.2 is a patch release to ensure consistent config and app data paths on windows. Both the config and app data paths are now set by default to ``C:\\Users\\<user>\\AppData\\Roaming\\DTU Wind Energy\\windkit``.
+
+## 1.0.1 (2025-07-15)
+
+WindKit 1.0.1 is a patch release that changes the config migration that shipped with 1.0.0. The changed ensures that the legacy `.ini` configuration files are migrated to the new `.toml` format without automatically removing the old file. This ensures that users are not constantly prompted to update their configuration file if they use different versions of WindKit in different environments.
+
+## 1.0.0 (2025-06-26)
+
+This is a major release that introduces a significantly cleaner and more intuitive API. We have removed many functions from the public namespace, renamed others for consistency, and introduced powerful new features. While this version includes many breaking changes, the result is a more streamlined and user-friendly library. We strongly recommend all users to carefully read these release notes and update their code.
+
+### Highlights
+
+#### Simplified and Consistent API
+The `windkit` namespace has been significantly cleaned up to make the API more predictable and easier to learn. Many internal or redundant functions have been removed from the public API, and others have been renamed for better clarity and consistency.
+*   Functions for reading and writing different objects in windkit are now grouped by type, such as `read_tswc`, `read_bwc`, `read_wwc`, and `read_gwc` for reading wind climate datasets, and `bwc_to_file`, `wwc_to_file`, and `gwc_to_file` for writing them.
+*   Functions for creating empty datasets are now prefixed with `create_*` (e.g., `create_bwc` instead of `empty_bwc`).
+*   Functions for creating wind climate objects from other data types now follow a `[climate]_from_[source]` pattern (e.g., `bwc_from_tswc` instead of `bwc_from_timeseries`).
+*   Validation functions are now consistently named `validate_*` (e.g., `validate_bwc` instead of `bwc_validate`).
+
+#### Revamped Map Reading and Writing
+The API for handling spatial maps has been completely overhauled. The generic `read_raster_map`, `read_vector_map`, `raster_map_to_file`, and `vector_map_to_file` functions have been replaced by specific, format-aware functions for each map type:
+*   `read_elevation_map` and `elevation_map_to_file`
+*   `read_roughness_map` and `roughness_map_to_file`
+*   `read_landcover_map` and `landcover_map_to_file`
+
+These new functions automatically infer the file format from the suffix (e.g., `.map`, `.tif`, `.gpkg`) and provide a more robust and intuitive interface for working with elevation, roughness, and landcover data.
+
+#### Polygon Representation for Landcover and Roughness Maps
+A major improvement is the new default support for polygon-based landcover and roughness maps. The functions `read_landcover_map` and `read_roughness_map` now return `geopandas.GeoDataFrame` objects with polygons when reading vector data, which is a more standard and powerful representation for GIS data. This enables more robust handling and validation of maps. For backward compatibility, you can still work with the legacy line-based format by setting `polygons=False`.
+
+#### Unified API for Wind Climate Statistics
+Calculating statistics from wind climate objects is now much simpler. A new set of generic functions can operate on any wind climate type (BWC, WWC, etc.):
+*   `wk.mean_wind_speed`
+*   `wk.mean_power_density`
+*   `wk.ws_cdf` (Cumulative Distribution Function)
+*   `wk.ws_freq_gt_mean`
+*   `wk.mean_ws_moment`
+
+This unified API replaces the climate-specific functions (e.g., `bwc_mean_windspeed`), simplifying your code and making it more versatile.
+
+#### New Loss and Uncertainty Modeling
+WindKit now includes tools to support loss and uncertainty analysis, designed to integrate seamlessly with PyWAsP. You can now create, validate, and summarize standardized loss and uncertainty tables using:
+*   `get_loss_table`, `total_loss`, `loss_summary`
+*   `get_uncertainty_table`, `total_uncertainty`, `uncertainty_summary`
+
+This provides a structured framework for AEP calculations.
+
+### Breaking Changes and Migration Guide
+Due to the extensive API changes, you will need to update your existing code. Here is a guide to the most common changes:
+
+#### General Renaming and Removals
+*   **Namespace Cleanup**: Many functions like `bwc_validate_wrapper`, `update_history`, and various `spatial` utilities have been removed from the public `windkit` namespace. 
+*   **Function Renaming**:
+  *   `windkit.read_timeseries_from_pandas` is now `windkit.tswc_from_dataframe`.
+  *   `windkit.bwc_from_timeseries` is now `windkit.bwc_from_tswc`.
+  *   `windkit.read_timeseries_from_csv` is now `windkit.read_tswc`.
+  *   All `windkit.empty_*` functions are now `windkit.create_*` (e.g., `empty_bwc` -> `create_bwc`).
+  *   `lines2poly` is now `lines_to_poly`, and `poly2lines` is now `poly_to_lines`.
+  *   `add_crs` is now `set_crs`.
+*   **Argument Renaming**: Several function arguments have been renamed for consistency. For example, `nsec` is now `n_sectors` and `nwsbins` is now `n_wsbins`. Check the function signatures if you encounter errors.
+*   **Removed Classes**: The `windkit.WindTurbines` class has been removed. Use the functions in the `windkit.wind_turbines` module to work with `xarray.Dataset` representations of turbines.
+
+
+#### Wind Climate I/O 
+Wind climate reading and writing has been standardized to one reader and one writer for each type:
+  *   Use `read_tswc` to read time series wind climate datasets.
+  *   Use `read_bwc` and `bwc_to_file` for Binned Wind Climate datasets.
+  *   Use `read_wwc` and `wwc_to_file` for Weibull Wind Climate datasets.
+  *   Use `read_gwc` and `gwc_to_file` for Generalized Wind Climate datasets.
+
+  #### Configuration File and Environment Variables
+  The configuration file format has been switched from INI to TOML. Existing `.ini` configuration files will be migrated automatically.
+
+  Additionally, user details can now be configured using case-insensitive environment variables, which is ideal for containerized or CI/CD environments:
+  *   `WINDKIT_NAME`
+  *   `WINDKIT_EMAIL`
+  *   `WINDKIT_INSTITUTION`
+
+  If no configuration is found (either via file or environment variables), you will be prompted for the information upon initialization.
+
+
+
+### Full list of changes
+
+#### Bug Fixes
+*   `clip_with_margin` now returns a `ValueError` when the bbox crosses the dateline.
+*   Fixed an issue where `read_raster_map` was checking for too high roughness values twice.
+*   Duplicate heights passed to `create_dataset` with `struct=stacked_point` are now correctly removed.
+*   `gdf_to_ds` now correctly handles unique points to ensure round-tripping with `ds_to_gdf` works as expected.
+*   Fixed a bug where `get_vector_map` for CORINE data could fail with a missing column error.
+*   `to_stacked_point` now correctly retains the input order of coordinates.
+*   `nearest_points` now raises an error on duplicate points to be compatible with newer xarray versions.
+*   Fixed a bug in `wk.spatial.add_crs` where custom projection information could be overwritten by a detected EPSG code.
+*   Fixed `wk.spatial.interp_unstructured_like` to correctly broadcast single-point sources.
+*   Fixed a bug in `nearest_points` when `keep_duplicates=False` and `return_distance=True`.
+
+
+#### API Breaks
+* `windkit.read_timeseries_from_pandas` renamed as `windkit.tswc_from_dataframe`
+* `windkit.bwc_from_timeseries` renamed as `windkit.bwc_from_tswc`
+* `windkit.read_timeseries_from_csv` renamed as `windkit.read_tswc`
+* all `windkit.empty_*` functions renamed to `windkit.create_*`
+* renamed `lines2poly` to `lines_to_poly` and `poly2lines` to `poly_to_lines`
+* removed `wk.lib_string` from public namespace
+* `windkit.spatial.interp_unstructured_like` and `windkit.spatial.interp_structured_like` will not interpolate over 1D variables by default
+ * Removed `read_raster_map`,`raster_map_to_file`, and `read_vector_map` and `vector_map_to_file` (made private)
+ * Added `read_elevation_map` and `elevation_map_to_file` (works for vector and raster)
+ * Added `read_roughness_map` and `roughness_map_to_file` (works for vector and raster)
+ * Added `read_landcover_map`, and `landcover_map_to_file` (works for vector and raster)
+* Old .gpkg files with change lines in landcover_map_to_file
+and roughness_map_to_file can be written by using polygons=False
+ * File formats for the readers and writers described above are now inferred from the filename's suffix, e.g., `.map` for map files.
+ * The functions `read_landcover_map` and `read_roughness_map` by default return polygons instead of changes lines. You can set the argument `polygons=False` to open the file as old fashioned change lines. You can try this if your file is not being opened correctly, but you are responsible for checking it's correctness, because there is less error checking for these types of files.
+ * Added `roughness_to_landcover` to convert roughness maps to landcover maps (vector or raster)
+ * Added `landcover_to_roughness` to convert landcover maps to roughness maps (vector or raster)
+ * Renamed `add_crs` to `set_crs` and and deprecated `add_crs` (which calls `set_crs`)
+* `get_vector_map` now returns landcover data as polygons by default.
+ * Removed `wk.bwc_ws_moment` from public namespace, use `wk.mean_ws_moment` instead
+ * Removed `wk.bwc_mean_windspeed` from public namespace, use `wk.mean_wind_speed` instead
+ * Removed `wk.bwc_mean_windspeed3` from public namespace, use `wk.mean_wind_speed3` instead
+ * Removed `wk.bwc_power_density` from public namespace, use `wk.mean_power_density` instead
+ * Removed `wk.bwc_ws_cdf` from public namespace, use `wk.ws_cdf` instead
+ * Removed `wk.bwc_ws_freq_gt_mean` from public namespace, use `wk.ws_freq_gt_mean` instead
+ * Removed `wk.bwc_to_tabfile` from public namespace, use `wk.bwc_to_file` instead
+ * Removed `wk.bwc_validate_wrapper` from public namespace
+ * Removed `wk.create_time_attributes` from public namespace
+ * Removed `wk.count_to_ws_freq_by_sector` from public namespace
+ * Removed `wk.read_rsffile`  from public namespace use `wk.read_wwc` instead
+ * Removed `wk.read_wrgfile` from public namespace, use `wk.read_wwc` instead
+ * Removed `wk.read_grdfile` from public namespace, use `wk.read_wwc` instead
+ * Removed `wk.read_pwcfile` from public namespace, use `wk.read_wwc` instead
+ * Removed `wk.to_rsffile` from public namespace, use `wk.wwc_to_file` instead
+ * Removed `wk.to_wrgfile` from public namespace, use `wk.wwc_to_file` instead
+ * Removed `wk.wwc_validate_wrapper` from public namespace
+ * Removed `wk.wwc_validate_structure`
+ * Removed `wk.wwc_validate_structure_wrapper`
+ * Removed `wk.wwc_mean_windspeed` from public namespace, use `wk.mean_wind_speed` instead
+ * Removed `wk.wwc_power_density` from public namespace, use `wk.mean_power_density` instead
+ * Removed `wk.to_libfile` from public namespace, use `wk.gwc_to_file` instead
+ * Removed `wk.gwc_validate_wrapper` from public namespace
+ * Removed `wk.wtg_validate_wrapper` from public namespace
+ * Removed `wk.windturbines_validate_wrapper` from public namespace
+ * Modified `wk.read_tswc` to be a generic tswc reader (not just for CSV)
+ * Renamed `wk.is_windturbines` to `wk.is_wind_turbines`
+ * Renamed `wk.windturbines_validate` to `wk.validate_windturbines`
+ * Renamed `wk.wtg_validate` to `wk.validate_wtg`
+ * Renamed `wk.gwc_validate` to `wk.validate_gwc`
+ * Renamed `wk.wwc_validate` to `validate_wwc`
+ * Renamed `wk.bwc_validate` to `wk.validate_bwc`
+ * Removed `wk.get_map`
+ * Removed `wk.spatial.interpolate_to_grid`
+ * Removed `wk.WindTurbines`
+ * Removed `wk.Workspace.get_vector_map`
+ * Added `wk.Workspace.get_elevation_map`
+ * Added `wk.Workspace.get_roughness_map`
+ * Removed `wk.bwc_from_counts` (from public namespace)
+ * Removed `wk.create_wv_count` (from public namespace)
+ * Removed `wk.update_history` (from public namespace)
+ * Removed `wk.update_coords_attrs` (from public namespace)
+ * Removed `wk.update_var_attrs` (from public namespace)
+ * Removed `wk.ALL_VARS_META` (from public namespace)
+ * Removed `wk.spatial.spatial_stack` (from public namespace)
+ * Removed `wk.spatial.spatial_unstack` (from public namespace)
+ * Removed `wk.spatial.stack_then_unstack` (from public namespace)
+ * Removed `wk.spatial.get_latitude` (from public namespace)
+ * Removed `wk.spatial.is_vertical` (from public namespace)
+ * Removed `wk.import_optional_dependency` (from public namespace)
+ * Removed `wk.reproject_gwc` (from public namespace)
+ * Removed `wk.get_ee_map` (from public namespace)
+ * Removed `wk.create_coords` 
+ * Renamed `wk.create_ws_bin_coords` to `wk.create_wsbin_coords` and changed the arguments to `bins`, `width` and `start`. `bins` is an integer number of bins or a list of bin edges. `width` is the width of the bins, and `start` is the starting point of the first bin. `width` and `start` are only used if `bins` is an integer.
+ * Removed `wk.create_ws_bin_coords_from_values` 
+ * Changed `wk.create_sector_coords` arguments to `bins` and `start`. `bins` is an integer number of bins or a list of bin edges. `start` is the starting point of the first bin. By default, `start` is -360/n_bins/2.
+ * Removed `wk.create_sector_coords_from_edges` 
+ * Removed `wk.create_direction_coords` (from public namespace)
+ * `wk.create_vector_map` now has bbox as a first positional argument and argument `dist` is removed
+ * Renamed `wk.get_land_cover_table` to `wk.get_landcover_table`
+
+#### New features
+* `wk.read_wtg` can read wind turbine generators files in netCDF format.
+* `wk.mean_ws_moment` to calculate the mean of a wind speed moment of any wind climate
+* `wk.ws_cdf` to calculate the cumulative distribution function of a wind speed of any wind climate
+* `wk.ws_freq_gt_mean` to calculate the fraction of probability mass above the mean of a wind speed of any wind climate
+* `wk.mean_wind_speed` to calculate the mean wind speed of any wind climate
+* `wk.mean_wind_speed3` to calculate the mean wind speed cubed of any wind climate
+* `wk.mean_power_density` to calculate the mean power density of any wind climate
+* All of the above calculates the the variable of all wind climates (not all functions are supports all climates yet).
+* Added `read_wwc` to read a weibull wind climate dataset from a file
+* Added `read_mfwwc` to read a multi-filey weibull wind climate dataset
+*  Added `wwc_to_file` to write a weibull wind climate dataset to a file
+* Added `read_gwc` to read a generalized wind climate dataset from a file
+* Added `gwc_to_file` to write a generalized wind climate dataset to a file
+* Added `bwc_to_file` to write a binned wind climate dataset to a file
+* New function `add_projected_wrf_coordinates` that is useful to process `WRF` output data to be used with `windkit`.
+* Existing function `windkit.read_wtg` can optionally read the wind turbine manufacturer name from a `.wtg` file if exists.
+* New function `windkit.metadata.update_coords_attrs` to update attributes of a dataset coordinates.
+* New function `create_baro` to create a baro file with dummy values
+* New function `create_stab` to create a stab file with dummy values
+* New function `windkit.snap_to_layer` that transform a vector map with polyons with gaps, open areas and inconsistent vertices into a valid vector map to run `pywasp`.
+* New example notebook tutorial for `windkit.get_era5` to get ERA5 data from Microsoft Planetary Computer
+* `windkit.raster_map_to_file` allows extra arguments that are forwarded to `rasterio.open`. Useful to create e.g. compressed geoTIFF files with the "compress" option.
+* New class `windkit.WengWorkspace` to read WAsP Engineering proyect files in XML format and converting its content to `xarray.Dataset` and `xarray.DataArray` instances such as extreme wind climate datasets, location datasets and raster map data arrays.
+* Enable roughness maps to be written to `.map` files
+* Added kw-argument `return_lctable` to functions reading landcover maps, set to false by default, which means only the landcover map is returned
+* Added kw-argument `convert_to_landcover` to functions reading roughness maps to optionally convert the map to landcover, returning both the map and the table (two outputs)
+* It is now allowed to use `wk.landcover_map_to_file` to write a landcover map to a `.map` file as long as the displacements are all zero.
+* `lines_to_polygons` is a more robust version of `lines_to_poly` and now tries to make the polygons valid if they are not. It is also much faster.
+* `polygons_to_lines` is a more robust version of `poly_to_lines` and is much faster.
+* Added function `windkit.plot.landcover_map`, which can plot a landcover map with some nicer color schemes.
+* Added function `windkit.plot.elevation_map`, which can plot an elevation map with some nicer color schemes.
+* Added function `windkit.add_landcover_table` to add roughness and other information from a wk.LandCoverTable to a geopandas dataframe with only an `id` column
+* in `get_raster_map`, an option to download from GWA4 landcover `GWA4_landcover` and `GWA4_elevation` has been added
+* Added function `windkit.get_tutorial_data` to download tutorial data from the windkit zenodo repository. The data is downloaded to a local cache directory, which is created if it does not exist.
+* Configuration can now be defined using case-insensitive environment variables `windkit_name`, `windkit_email` and `windkit_institution`. If these are not set, and the configuration file is not found, you will be prompted on initialization.
+
+#### Added functionality to work with losses
+- Defined the data structure of loss tables: a pandas DataFrame with required columns: "loss_name" and "loss_percentage". It can also have optional columns "loss_category", "loss_upper_bound", "loss_lower_bound", "loss_default", and "description"
+- Added function `get_loss_table` to obtain a loss tables (DataFrame) from known name. Currently, only "dtu_default" exists.
+- Added function `validate_loss_table` to validate a loss table is formatted correctly
+- Added function `total_loss` to calculate the total loss in percentage from a loss table
+- Added function `loss_summary` to print a formatted summary of the loss table
+- Added function `total_loss_factor` to calculate the loss factor (the factor that can be directly multiplied to subtract the loss from an AEP number)
+
+#### Added functionality to work with uncertainties
+- Defined the data structure of uncertainty tables: a pandas DataFrame with required columns: "uncertainty_kind", "uncertainty_name" and "uncertainty_percentage". It can also have optional columns "uncertainty_category", "uncertainty_upper_bound", "uncertainty_lower_bound", "uncertainty_default", and "description"
+- Added function `get_uncertainty_table` to obtain an uncertainty table (DataFrame) from known name. Currently, only "dtu_default" exists.
+- Added function `validate_uncertainty_table` to validate an uncertainty table is formatted correctly
+- Added function `total_uncertainty` to calculate the total uncertainty in percentage from an uncertainty table. Following the GUM guideline on adding uncertainties.
+- Added function `uncertainty_summary` to print a formatted summary of the uncertainty table
+- Added function `total_uncertainty_factor` to calculate the uncertainty factor (the factor that can be directly multiplied to add the uncertainty to an AEP estiamtion)
+
+#### Bug fixes
+* `clip_with_margin` now returns a `ValueError` when the bbox crosses the dateline
+* `read_raster_map` was checking for too high roughness twice, one is enough
+* If duplicate heights are passed to `create_dataset` with `struct=stacked_point` it will remove the duplicates
+as required in the `stacked_point` structure.
+* in `gdf_to_ds` we must create the target dataset using unique (`west_east`,`south_north`) points for `stacked_point` and unique `west_east` and `south_north` for cuboids to be able to round trip `ds_to_gdf` and `gdf_to_ds`.
+* get_vector_map(bb, "CORINE", source="dtu") could sometimes fail with "Exception: Missing required column: id", this has been fixed
+* read_vector_map also reads .gpkg files where a geometry column is missing on the landcovertable layer
+* to_stacked_point now retains input order of west_east/south_north coordinates
+* Nearest_points will now raise an error if there are duplicate points resulting from the nearest neighbour lookup if the data is in stacked_point or cuboid format. This is required because xarray >2024.02 does not allow duplicates in MultiIndexes anymore. A new keyword argument keep_duplicates controls the behaviour, the default being True, i.e. when possible behaviour is the same as before.
+* In the function `wk.spatial.add_crs` assigning a projection with more detailed info given in a proj4 string or custom CRS was overwritten when an EPSG was found for the projection. This is now fixed and the full custom CRS is added.
+* Fixed `wk.spatial.interp_unstructured_like` to broadcast single-point sources correctly to target.
+* Fix bug in `nearest_points` when `keep_duplicates=False` and `return_distance=True`: it threw an error that the arrays were not the same size. 
+
+#### Breaking changes
+* `empty_wasp_site_factors` renamed to `empty_wasp_site_effects`
+* in `empty_z0meso` argument `site_factors` was renamed renamed to `site_effects`
+* in `empty_pwc` argument `site_factors` was renamed renamed to `site_effects`
+- Argument `nwsbins`, `nws`, `n_bins` has been renamed to `n_wsbins` or `wsbins` (in case of array-like argument) in the following functions:
+    - `bwc_from_tswc`
+    - `create_ws_bin_coords`
+    - `create_ws_bin_coords_from_values`
+- Argument `nsec`, `nsecs`, `nwd`, `nbins` has been renamed to `n_sectors` in the following functions:
+    - `create_sector_coords`
+    - `create_direction_coords`
+    - `create_sector_coords_from_edges`
+    - `wd_to_sector`
+    - `bwc_from_tswc`
+- Argument `num_args` has been renamed to `n_args` in the following functions:
+    - `ltc_validate`
+- Argument `ws_bin_width`has been renamed to `wsbin_width` in the following functions:
+    - `create_ws_bin_coords`
+    - `bwc_from_tswc`
+* Metadata variable name `aep_deficit` has been renamed to `potential_aep_deficit`
+
+* Functions `read_raster_map` and `raster_map_to_file` were removed in favour of the new functions `read_elevation_raster`, `write_elevation_raster`, `read_landcover_raster`, and `write_landcover_raster`.
+* Functions `read_vector_map` and `vector_map_to_file` were removed in favour of the new functions `read_elevation_vector`, `write_elevation_vector`, `read_landcover_vector`, and `write_landcover_vector`.
+* Changed most instances where "header" is used to "description"
+* `get_vector_map` returns polygons instead of lines.
+* Switched configuration file from ini to toml, old file will be automatically migrated
+
+
+## 0.8.0 (2024-06-10)
+
+### New features
+
+#### Enhanced wind turbines class
+
+* New functions in `windkit.wind_turbines` to create point `xarray.Dataset` representation of WTGs from `pandas.DataFrame`, `geopandas.DataFrame` and raw point arrays. Additionally, the module adds a function `check_wtg_keys` to validate wtg keys and `wind_turbines_to_geodataframe` to convert the dataset back to a `geopandas.GeoDataFrame`.
+
+
+#### New Long Term Correction module
+
+The new module `ẁindkit.ltc` contains tools to perform long term correction of wind measurements. It is based on the windkit `xarray` time series wind climate dataset, therefore it
+is compatible with all the `windkit` toolset.
+* There are two methodologies available:
+  - class `LinRegMCP` implements a linear regression.
+  - class `VarRatMCP` implements a variance ratio linear regression.
+* There is also a  `ltc.scores` submodule to calculate scores such as `r2`, `rmse`, `bias` and `wasserstein distance`.
+
+
+#### Jupyer notebook examples
+
+The windkit package includes a folder with examples of use. currently there is one jupyter notebook example to use the long term correction module.
+
+
+#### New functions to handle wind data
+* New function `windkit.wd_to_sector` can return sectors of fixed width or sectors of equal probability (quantiles=True)
+* New function `resample_wind_and_direction` to resample a time series wind climate `xarray.Dataset` to a different sampling frequency.
+* New function `windkit.sector.create_sector_coords_from_edges` that returns a data array with sector coordinates for any array of sector edges, not necessarily equidistant.
+
+
+#### New function to implement cross prediction
+When implementing a cross prediction methodology, we need to create a combination of pair of spatial points with the different locations to predict.
+The new function `windkit.get_cross_predictions` allows that, from a `xarray.Dataset` with several points, return two aligned datasets with the points where we want
+to predict from and where we want to predict to.
+
+#### New functions on empty module
+* New function `windkit.empty.empty_raster_map` generates a `xarray.Dataset` raster map instance filled with either `NaN` or meaningful values.
+* New function  `windkit.empty.empty_wv_count` returns random meaningful values and it is more suitable for testing
+
+#### Improvements to get raster and vector source data
+* `get_vector_map` gets data in vector format from DTU servers (currently only support `CORINE`)
+* `get_raster_map`: added new source dtu, with access to `Viewfinder` and `ESA_CCI`, also allows non-geographic bounding boxes
+* Allow user to put `earth_engine.json` in home directory for authenticating to usage of google earth engine in `get_raster_map`.
+
+#### New functionality to handle conversion between geopandas dataframes and xarray datasets
+* New function `gdf_to_ds`, which converts a geopandas dataframe to a xr.Dataset in the best fitting windkit spatial structure.
+* New function `ds_to_gdf`, which converts a xr.Dataset in windkit spatial structure to a geopandas dataframe.
+
+
+### Changes
+
+* Added `flow_sep_height` argument in the metadata. This variable is used in pywasp to denote the height where flow separation occurs.
+* Added geostrophic wind climate variables `geo_turn` and `geo_wv_freq` to metadata
+* Moved `gwc_interpolate` to pywasp, because it relies on some code from pywasp and is only useful in a pywasp context
+* `get_roughness_map` and `get_raster_map` issues a warning if a landcover table is attempted to be read as roughness.
+* Removed constraint on first `gen_height` to be equal to 10. This is not a requirement for pywasp to work.
+* GML files with an empty description (the default) written by windkit can now be opened in the map editor
+* `bwc_to_tabfile` now creates the parent directory if does not exist, and formats the tab filename based on the dataset coordinates.
+* `windkit.Workspace.read_wwh` allows WAsP workspaces with displacement height maps (i.e. GML files, WAsP >12.7) to be imported.
+* New function `bwc_validate_structure` which checks only if the structure is OK for a `bwc` dataset, without checking for the actual values.
+* `windkit.bwc_to_tabfile` writes a default header "No Header" in the tab file if the dataset does not have the "wasp_header" attribute.
+* Added 'use_bounds' argument to `windkit.spatial.BBox.reproject` to allow the reprojection to respect the concave/convex shape of the bounding box due to reprojection.
+* The order of stacked point's is now preserved when round-trip via `windkit.spatial.spatial_stack` and `windkit.spatial.spatial_unstack`.
+* Added additional checks for consistency with PyWAsP use when reading in a roughness map with `windkit.read_raster_map`
+* Added method `windkit.WindTurbines.to_wtg_dict_and_turbines_ds` to convert the deprecated `WindTurbines` object to a dictionary and a `xarray.Dataset` of turbines.
+
+### Bug fixes
+
+* Fixed a bug in `windkit.read_bwc` , where it did not check correctly whether the requested crs matches the dataset crs.
+* `get_raster_map` no longer gives -9999 for areas where no data is available.
+* Fixed a bug in `windkit.bwc_to_tabfile`, where it allowed to write datasets with projected coordinates to a tab file, which is not allowed.
+* Fixed a bug in `windkit.bwc_from_timeseries` where it returned a point structure when the input was a stacked_point.
+* Fixed a bug in `LandCoverTable` that did not allow to instantiate an object when colors or labels are not `None`
+
+### Deprecations
+
+* `count_to_ws_freq_by_sector` is deprecated in favor of `bwc_from_counts`.
+* `binned_wind_climate.bwc_from_weibull` is removed in favor of `weibull_wind_climate.wwc_to_bwc` with simpler API passing in a `xr.Dataset` instead of weibull parameters one by one.
+* `WindTurbines` class is deprecated in favor of `wind_turbines.create_wind_turbines_from_*` functions to create `xr.Dataset` representation of WTGs.
+* `get_map` and `get_ee_map` are deprecated in favour of `get_raster_map`
+
+
+
+## 0.7.0 (2023-11-13)
+
+### New Features
+
+#### Additional interpolation functionality
+
+The interpolation functionality of windkit was expanded with the release, several low-level interpolation routines were added, allowing interpolation to and from all of windkit's spatial data structures.
+
+* Added `windkit.spatial.interp_structured_like` for interpolation from a "cuboid" to any spatial structure, using any of the interpolation methods available in `xarray.Dataset.interp`.
+* Added `windkit.spatial.interp_unstructured_like` and `windkit.spatial.interp_unstructured` and for interpolation between any spatial structures, using `nearest`, `linear`, or `cubic` methods. Additionally, the `natural` method can be used peform natural neighbor interpolation for 2D `point` structured target structures, using `metpy's implementation <https://unidata.github.io/MetPy/latest/examples/gridding/Natural_Neighbor_Verification.html>`.
+
+Additionally a few high level routines were added for interpolation specific fields.
+
+* Added `windkit.interpolate_gwc` to interpolate a generalized wind climate dataset.
+* Added `windkit.vinterp_wind_speed` for vertical interpolation of wind speeds
+* Added `windkit.vinterp_wind_direction` for vertical interpolation of wind directions
+
+#### Improvements to the map_conversion routines
+
+* Added `windkit.map_conversion.find_duplicate_lines` to identify duplicated lines in a GeoDataFrame, even if the line is reversed.
+* `windkit.poly_to_lines` now correctly sets the CRS to match the input data. Before it was `None`.
+* `windkit.poly_to_lines` and `windkit.lines2poly` now have an absolute tolerance argument, `atol`, which can be used to identify nearby points, which should be joined.
+
+#### Access to data from Microsoft's Planetary Computer platform
+
+* Elevation and landcover data obtained using `windkit.get_map` can now use either Microsoft's Planetary Computer platform or Google's Earth Engine.
+* Added `windkit.get_era5` to download ERA5 data from Microsoft's Planetary Computer platform.
+
+#### Calculate additional derived quantities from Weibull and Binned Wind climates
+
+* Added new functions `windkit.fit_weibull_wasp_m1_m3_fgtm`, `windkit.fit_weibull_wasp_m1_m3`, and `windkit.fit_weibull_k_sumlogm` for fitting weibull distributions to wind speed moments.
+* Added `windkit.weibull_cdf` and `windkit.weibull_pdf` calculate the cdf and pdf of a weibull distribution.
+* Added `windkit.weibull_freq_gt_mean` to calculate the fraction of probability mass above the mean of a weibull distribution.
+* Added `windkit.bwc_ws_cdf` and `windkit.bwc_freq_gt_mean` are similar functions for binned wind climates.
+
+#### Convenience functions to access standard land cover tables for common datasets
+* Added `windkit.LandCoverTable.get_table` and it's top-level alias `windkit.get_land_cover_table` to get a land cover table for a given land cover dataset. Landcover tables are currently provided for CGLS-LC100, CORINE, GlobCover, MODIS, WorldCover, and ESA_CCI.
+
+#### Even more spatial functionality
+
+* Added `windkit.spatial.clip_with_margin` to clip a raster-like dataset to the bounds of another dataset plus a margin. This is useful for clipping a dataset before doing other operations with it (e.g. interpolation) to reduce computation time.
+* Added `windkit.spatial.equal_spatial_shape`, to check if two xarray objects have the same struct and spatial shape. This differs from `windkit.spatial.are_spatially_equal`, which requires the values to also be the same.
+* Added `windkit.spatial.covers` to check if a cuboid dataset covers another spatial dataset (point, stacked_point,  or cuboid, or polygon).
+
+#### Further improvements to wrg and rsf file handling
+
+* `windkit.read_wrgfile`, `windkit.read_rsffile` have a new boolean argument `use_production`, where if it is True, the power colum in the file is interpreted as gross AEP, and when False, the power column in the file is interpreted as power density.
+* `windkit.read_wrgfile`, `windkit.read_rsffile` have a new boolean argument `use_production` where if it is True, the `gross_aep` data variable is written to the file, and if False, the `power_density` column is written to the file.
+
+
+### Bug fixes
+
+* `windkit.LandCoverTable.from_dict_ora`: Sets the correct values of `z0frac` and `dfrac` in the LandCoverTable description.
+* `windkit.BBox.buffer` now works with `shapely 2.0`.
+* Changed the default `nodata` value for `windkit.spatial.warp` from `np.nan` to `None`, which results in a value of `-9999` for integer fields, and `np.nan` for float fields. This fixes an issue with integer data, for which `np.nan` is undefined.
+* `windkit.read_rsffile` no longer returns a unneeded `point` coordinate, when `to_cuboid=False`.
+* `empty.empty_gwc` uses the same random number generator for all data variables, ensuring the same seed gives the same numbers.
+
+### Changes
+* `windkit.generalized_wind_climate.read_gwc` performs validations checks for all formats, not only `netCDF`
+* `windkit.generalized_wind_climate.to_libfile` detects if path is a file or a directory and creates subdirectories if they do not exist.
+* `windkit.generalized_wind_climate.to_libfile` formats `.lib` filename using the coordinates if no name is given, similar to what is done with multiple files.
+* function `gwc_validate` and decorator`gwc_validate_wrapper` perform extra validations on generalized wind climate datasets. The first generalized height must be 10.0, the first generalized roughness must be 0.0, and both must have at least 2 entries.
+* `windkit.read_cfdres`, now always returns a dataset with the coordinates `south_north` and `west_east` sorted from low to high.
+* `windkit.read_vector_map` will raise `RuntimeError` when trying to read a `.map` file with displacement height, which is not supported.
+* `wk.read_rsffile` now allows .rsf files with a single-point to be converted to cuboid.
+* Added `windkit.rotor_equivalent_wind_speed` function for calculating REWS
+* Added `windkit.read_wwc` that creates and validates a weibull wind climate `xarray.Dataset` from a file.
+* Added `datum_name` argument to `windkit.spatial.BBox.utm_bbox_from_geographic_coordinate` to allow the user to specify the datum name for the UTM projection.
+* `empty.empty_wasp_site_factors` and `empty.empty_met_fields` are filled with random numbers.
+* `empty.empty_pwc` can include any desired combination of met field and site factors. Defaults to `A_combined`, `k_combined`, `power_density` for met fields and `site_elev` for site factors.
+* `empty.empty_gwc` has a new boolean argument `wdfreq_constant`, which allows you to create a generalized wind climate object where `wdfreq` is constant along the `gen_height`. This is used to create datasets that can be written to the `.lib` file format.
+* `empty.empty_gwc` retains the shape of scalar structured datasets, i.e. a single point with no point dimension, only coordinate values. This structure is common when selecting a single location using `.sel` or `.isel`.
+* Added `windkit.empty.empty_tswc` that generates a time series wind climate `xarray.Dataset` with random but meaningful values.
+* Added `split_combomap` CLI utility, which reads a combo `.map` file and returns two `.gpkg` files, one for elevation and one for landcover.
+
+
+### Deprecations
+* `windkit.spatial.interpolate_to_grid` deprecated in favor of `windkit.spatial.interp_structured_like` and `windkit.spatial.interp_unstructured_like`
+
+### Breaking changes
+* `windkit.poly_to_lines` no longer has a `closed_map` argument. Instead, when `background_lc_id` is set, all empty space is filled with the provided id. This allows for maps that are no longer rectangular, e.g. due to reprojection, to be properly converted to line maps.
+* Removed accessors, which were previously registered for several wind climate data structures. You will need to use the functions instead.
+
+## 0.6.3 (2023-07-18)
+
+### New Features
+
+#### Enhanced bounding box object
+
+Several new methods were added to the BBox object to make it easier to work with.
+
+* Create BBox objects in more ways:
+    * `BBox.from_bounds` as alias to `BBox.from_cornerpts` and set default crs to `EPSG:4326`
+    * `BBox.from_point_and_buffer` to create a bbox from a x,y-point, a buffer size, and a crs
+    * `BBox.utm_bbox_from_geographic_coordinate` to create a bbox from a geographic coordinate (point) and a buffer in meters, projecting to the appropriate UTM coordinates.
+* Simplified reproject to geographic and UTM coordinates
+    * Added method `BBox.reproject_to_geographic()` to reproject a bbox to geographic coordinates
+    * Added method `BBox.reproject_to_utm()` to reproject a bbox to UTM coordinates
+* Return data in common GIS formats using:
+    * Added method `BBox.polygon` property to return a shapely polygon
+    * Added method `BBox.envelope()` to return an enveloped bbox in the current crs
+* Create a simple plot of your BBox using `BBox.plot()`
+
+
+#### Formatting improvements when writing .wrg/.rsf files
+
+* `windkit.to_rsffile` and `windkit.to_wrgfile` now attempts to match the WAsP formatting more closely and a new `formatters` argument allows the user to specify the formatting for each column. This allows for more control over the formatting of the output file. Additional formatting checks have also been added to ensure that the formatting is correct.
+
+#### Require regulation_type for wtg data structures
+
+* The WindTurbineGenerator now requires a regulation_type variable. When reading in a WTG using `windkit.read_wtg`, the default mode is pitch. You can set it to stall or have it autodetected using the new `windkit.wind_turbine.estimate_regulation_type` function.
+* Validation functions have been added to the WTG object type, so you can now ensure that your Dataset has the correct format before using the WTG.
+
+#### Allow use of Polygons for defining landcover data
+
+* `windkit.read_vector_map` Now can read polygon based landcover vector map files in `gpkg` format, which is commonly used in GIS tools. The function has a new parameter `background_lc_id` that is needed to setup the background for polygon based vector maps.
+
+### Changes
+* `windkit.read_vector_map`: For landcover data in `gpkg` format,`id` is now the default name used to identify landcover ID's in the landcover table. For old `gpkg` files where the name `index` is used, `windkit.read_vector_map` will print a `FutureWarning`.
+* Remove dependency on `cf-xarray`, as we have our own metadata vocabulary. This fixes an issue where Pooch was an optional dependency of `cf-xarray`, but required for the functionality used by WindKit.
+
+### BugFixes
+
+* Prevent exception when running in the python interpreter and calling `windkit.metadata.update_history`. Now a the generic message "Unknown python interpreter command." is added to the object history in this scenario. To get a correct history, it is recommended to use iPython for interactive sessions.
+* Fixed a bug in `BBox.buffer`, which would cause negative buffering to not preserve bbox shape
+
+## 0.6.2 (2023-06-20)
+
+### New Features
+
+#### New file reading functions
+
+* `windkit.read_timeseries_from_pandas` and `windkit.read_timeseries_from_csv` build a time-series wind climate `xarray.Dataset` from time-series data in tabular form, either a `pandas` Dataframe or a `csv` file. The functions attempt to auto-detect the wind speed and wind direction columns based on a standard list of names if they are not provided.
+* `windkit.read_cfdres` reads a WAsP CFD `.cfdres` file into an `xarray.Dataset` structure.
+
+#### Detection of CRS from WAsP .map files
+
+* `windkit.read_vector_map` and by extension `windkit.read_elevation_map` and `windkit.read_roughness_map` now try to detect CRS information from the header of WAsP `.map` files, which can be added using recent versions of the WAsP MapEditor. In case there is a problem (e.g. wrong format) it will raise an error.
+
+#### Updates to vertical profile plotting
+
+* `windkit.plot.vertical_profile` Function now works for data arrays without their attributes set. If the height coordinate does not have attributes, `Height [m]` will be used for the y-axis. If the data variable to plot does not have attributes, `DataArray.name` will be used for the x-axis.
+* `windkit.plot.vertical_profile` no longer allows plotting of objects with a `time` coordinate. This is due to performance limitations when plotting potentially thousands of points when a time-series object is passed.
+
+### Breaking changes
+
+* When creating a file using `windkit.to_rsffile` and `windkit.to_wrgfile` or reading a file using `wk.read_rsffile` and `wk.read_wrgfile`, the elevation variable has been renamed from `elevation` to `site_elev`. This now matches the variable name used in PyWAsP.
+* class `windkit.WindTurbines` does not allow wind turbine `xarray.Dataset` with more than one mode.
+
+### Changes
+
+* `windkit.spatial.nearest_points` will issue a warning instead of a `ValueError` when it is used on datasets with geographical coordinates.
+* `windkit.to_rsffile` and `windkit.to_wrgfile` calculate necessary parameters met_fields, "power_density", "combined_A", "combined_k", and "wspd", from a wwc if they aren't included.
+
+### Bug Fixes
+
+* Fixed a formatting bug in `wk.to_rsffile` and `wk.to_wrgfile` that resulted in incorrect column widths.
+
+## 0.6.1 (2023-03-16)
+
+### New Features
+* `windkit.wtg_cp` function returns the power_coefficient for given wind speeds from an wind turbine generator.
+* Added lookup table to convert WAsP Map Editor named projections to EPSG codes. This allows for projection to be autodetected from some GWC files.
+
+### Breaking changes
+
+* Removed `coords` from `windkit.read_gwc`. Instead you should use the separate `west_east`, `south_north`, and `height` arguments. Also removed `**kwargs`, as there was only a single option `header`, which is now its own argument.
+
+### Changes
+
+* Updated pyproj dependency to minimum 3.0
+* Added `gen_heights` and `gen_roughnesses` arguments to `windkit.empty.empty_gwc` to allow users to specify these coordinate values.
+* Removed the `mode` argument from `windkit.wtg_power`, `windkit.wtg_ct`, and `windkit.wtg_cp`, as they now work with any number of modes in the WTG. If you wish to use a single mode from a wind turbine object, you should select using the `isel` or `sel` method.
+
+### Bug Fixes
+
+* Allow projections that don't have an EPSG code to be added using `windkit.spatial.add_crs`.
+* Scalar spatial datasets would fail in `windkit.spatial.count_spatial_points`.
+* Ensured that `wk.Workspace.get_gwc` could parse GWC objects from WAsP 12.7 and 12.8 workspaces.
+* Fixed bug in `windkit.spatial.to_stacked_point`, where the order of the `west_east` and `south_north` coordinates would be flipped from expected. This caused an issue when converting to `point` and then back to `stacked_point`.
+* Prevent creation of NaN values for some curves when using `windkit.read_wtg` to read ``.wtg`` files with different modes and different wind speed tables.
+* Bug fix in `wind_climate._is_bwc` where an exception was not properly caught.
+
+## 0.6.0 (2023-01-29)
+
+Several new features in this release including support for Wind Turbines, the ability to convert roughness lines to roughness polygons, and the ability to read PWC files. Several important bug fixes also added.
+
+### New Features
+
+#### Support for wind turbines and AEP data
+
+* New metadata attributes for annual energy production and wind farm flow map
+* New class `windkit.WindTurbines` to hold collections of turbine locations and WTGs (power/ct curves and other details).
+* Added metadata for wind turbine generator `xarray.Dataset`.
+* New plotting function `windkit.plot.plot_wind_turbine_locations`.
+
+#### Ability to convert roughness lines to polygons and back
+
+* New module `windkit.map_conversion` That allows the user to both create polygons from
+  lines and convert polygons to lines, based on `geopandas.GeoDataFrame`. It contains
+  the following:
+  - `LineMap` and `PolyMap` classes: These classes enable conversion and nice plotting
+    of their types.
+  - `lines2poly` and `poly_to_lines` helper functions that perform the conversion from
+    lines to polygons and from polygons to line respectively, using the library
+    `geopandas.GeoDataFrame`.
+* New module `windkit.plot.landcover_map` to plot the `geopandas.GeoDataFrame` polygons
+  or lines as a map.
+* New module `windkit.plot.color` to provide color functions for maps with a focus on
+  roughness coloring.
+* Class`windkit.LandCoverTable` has a new method `add_colors_to_table` to store colors to be used in map plots.
+
+#### New spatial functionality
+
+* New function `windkit.spatial.interpolate_to_grid` to resample data to a target grid.
+* New function `windkit.spatial._get_latitude` that returns data-array of latitude values in wgs84.
+* New function `windkit.spatial.count_spatial_points` counts how many points are in a windkit dataset or dataarray.
+
+#### Enhancements to `windkit.get_map.get_ee_map`
+* Added NASADEM as elevation dataset
+* Added worldcover as landcover dataset
+* Added option to convert to vector in the cloud.
+
+#### Ability to read PWC files
+
+* `windkit.read_pwcfile` reads WAsP `.pwc` predicted wind climate file in XML format and returns a weibull wind climate `xarray.Dataset`
+
+### Breaking changes
+
+* Requires `geopandas>=0.11`.
+* Requires `rioxarray` for reading raster GIS files.
+
+### Deprecations
+
+### Changes
+
+* To support `pyproj` 3.4.0, changed from `pyproj.crs.CRS.is_exact_same` to `pyproj.crs.CRS.equals`, and dropped deprecated `skip_equivalent` kwargs.
+* Drop `south_north` and `west_east` coords before reassigning them as required in [future `xarray` version](https://github.com/pydata/xarray/issues/6505).
+* Added height dependent wind direction frequencies to GWC I/O routines
+  - For the .lib-file reader frequencies are duplicated across all gen_heights
+  - For the .lib-file writer a check is made to ensure all frequencies are these same across heights. If they are not, a gen_height is required to be set by the user.
+
+### Improvements
+ * Sped up creation of .lib strings for GWC objects. In testing 15 sites, went from 11s to <1.
+ * Retain meta data of `wsbin` coordinates now also present after combining `bwc`'s
+
+### Bug Fixes
+* `Windkit.wwc_mean_windspeed`  and `windkit.wwc_power_density` previously calculated emergent wind case using all-sector A&k values. Now it is correctly calculated using sectorwise, and an error is thrown if both `emergent` and `bysector` kwargs are `True`.
+* Fixed reading `tif` and `grd` files with negative `dx`.
+* Catch exception in `wind_climate._is_bwc`.
+* Prevent modification of input dataset when calling `weibull_wind_climate.weibull_combined`.
+* Return point or stacked-point dataset for when calling `windkit.spatial.create_dataset` with duplicated `west_east` and `south_north` values. Previously a cuboid was returned with duplicated points collapsed.
+* Use fixed column parsing in `read_rsf/read_wrg` rather than space splitting, when reading the number of sectors.
+* Ensure temporary file created when using `windkit/get_map.py` is always removed.
+* Create correct BWC when reading data with with irregular bins (i.e. as exported from windpro).
+
+## 0.5.1
+
+* Changed name of author's attributes from name/email to author/author_email
+* Fixed bug for 1d nearest point estimation. The KDtree throws an error if we do not expand the 2d point [x,y] to [[x],[y]]
+* Removed unneeded dependencies
+
+##  0.5.0
+
+### New Features
+* Added functions `read_wrgfile` and `to_wrgfile` to read and write .wrf resource grid files
+* Added functions for simple wind-related calculations:
+   - `windkit.wind_speed_and_direction` calculates wind speed and direction from horizontal wind vectors U and V
+   - `windkit.wind_vectors` calculates horizontal wind vectors U and V from wind speed and direction
+   - `windkit.wind_direction_difference` calculates the smallest (signed) circular difference between wind directions in degrees.
+   - `windkit.wind_speed` calculates horizontal wind speed from vectors U and V
+   - `windkit.wind_direction` calculates horizontal wind direction from vectors U and V
+   - `windkit.wd_to_sector` to calculate the bin-index for wind directions for a given number of bins/sectors
+* Support for new Xarray-based Wind Turbine Generator (WTG) object added
+*  New wind climate dataset validator function and decorator.
+   - The function receives an `xarray.Dataset` and returns `True` or `False`.
+   - The decorator wraps a function that receives a wind climate dataset and raises a `WindClimateValidationError` if the dataset is not valid.
+   - It checks dimensions, attributes and data variables.
+   - The functions available are:
+     - `bwc_validate`,`bwc_validate_wrapper` for binned wind climate dataset, imported from `windkit.binned_wind_climate`
+     - `wwc_validate`,`wwc_validate_wrapper` for weibull wind climate dataset, imported from `windkit.weibull_wind_climate`
+     - `gwc_validate`,`gwc_validate_wrapper` for generalized wind climate dataset, imported from `windkit.generalized_wind_climate`
+     - `ts_validate`,`ts_validate_wrapper` for time series wind climate dataset, imported from `windkit.time_series_wind_climate`
+
+* New `WindClimateValidationError`, raised by the newly introduced wind climate dataset validators.
+* New method `BBox.to_geoseries()` that converts Bounding box to `geopandas.Geoseries`
+* New `empty` module that generates `xarray.Dataset` instances filled with either `NaN` or meaningful values.
+   -  Suitable for testing.
+   - Current available functions are:
+     - `empty.empty_wasp_site_factors()`: Create empty site-factors dataset
+     - `empty.empty_bwc()`: Create empty binned wind climate dataset filled with meaningful random numbers, i.e. they are generated
+    with a weibull distribution, the sum of `wdfreq` is 1, etc.
+     - `empty.empty_wwc()`: Create empty weibull wind climate dataset filled with meaningful random numbers, e.g.the values from `A`
+    are generated from a uniform function between `5` and `10` and the values for `k` from a uniform function between `1.5` and `2.5`.
+     - `empty.empty_gwc()`: Create empty generalized wind climate dataset with the data variables are filled with meaningful random
+    numbers, e.g. the values from `A` are generated from a uniform function between `5` and `10` and the values for `k` from a uniform function between `1.5` and `2.5`.
+     - `empty.empty_met_fields()`: Create empty dataset filed with `met_fields`
+     - `empty.empty_z0meso()`: Create empty site_factors with only `z0meso` and `slfmeso`
+     - `empty.empty_pwc()`: Create empty predicted wind climate with optional variables
+* New plotting functions in module `windkit.plot`
+   - `windkit.plot.vertical_profile()` plots the vertical profile of the dataArray or dataArrays introduced.
+   - `windkit.plot.time_series()` creates a time series plot.
+   - `windkit.plot.single_curve()` and `windkit.plot.power_ct_curves()` add plots for electrical power output curve, thrust curve, RPM, Ct, Cp or any other curve the user wants to pass through.
+
+* New function `windkit.time_series.read_ts_windpro_txt`
+   -  Reads windpro format txt file into a `xarray.Dataset`
+### Removed
+* Removed `windkit.binned_wind_climate.dir_to_sec` in favor of `windkit.wd_to_sector`
+* Modules `WaspVectorMap` and `WaspRasterMap` were moved to `pywasp` and are no longer available.
+ - Removed legacy WTG objet `WindTurbine`
+
+### Changes
+
+* `windkit/windturbine.py`
+    - Added function `read_wtg` to create an `xarray.Dataset` representation of a turbine
+    - Added `wtg_power` and `wtg_ct` to return power and ct for given wind speed
+    - Added `wtg_to_pywake`, which makes a PyWake wind turbine class from the dataset
+ - renamed `vectormap` to `vector_map`
+ - renamed `rastermap` to `raster_map`
+ - renamed `elevationmap` to `elevation_map`
+ - renamed `roughnessmap` to `roughness_map`
+ - `windkit.read_elevation_map` and `windkit.read_roughness_map` now reads both raster and vector maps.
+ - Changed the default behavior of  `windkit.spatial.clip` for rasters/cuboids so small padding is added to the mask, ensuring that raster/cuboid pixels/cells on the edge of the mask are included in the clipping. The previous behavior can be achieved by setting `pad=False`
+ - Updated `windkit.spatial.clip` for `point` `xr.dataset` to always include points that are on the edges of the mask (previously only points inside were included)
+ - Updated BWC "header" and GWC "desc" to both be named "wasp_header" (only included for objects originating from WAsP files or methods).
+ - Updated to always use "crs" as spatial reference argument name
+ - Mirrored `spatial.to_raster` as `spatial.to_cuboid`
+
+### Improvements
+ - Updated `read_rsffile` and `read_wrgfile` to automatically infer the number of sectors in the data.
+ - Updated `windkit.spatial.clip` so it uses a faster clipping method (`_clip_to_bbox_raster`) when the masking is provided as a `BBox` object. This speeds up clipping speed significantly.
+ - Allowed `ws_bins` to be python `range` in `windkit.binned_wind_climate._freqs_to_dataset()`
+ - Functions `windkit.spatial.to_point` and `windkit.spatial.to_raster_point` now also work on scalar datasets.
+ - Fix bug in `metadata.update_var_attrs` to store "Object type" as the string "None" to properly write to netCDF
+ - Landcover GMLs use custom reader/write to ensure compatibility with WAsP & Map Editor
+ - Plotting documentation improved by rendering jupyter notebooks on it using `nbsphinx` extension for `sphinx`
+ - Plotting testing with jupyter notebooks