.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "auto_examples/resource_grid.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_auto_examples_resource_grid.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_auto_examples_resource_grid.py:


Calculate Resource grid
===========================

Example of using PyWAsP to calculate a resource grid

.. GENERATED FROM PYTHON SOURCE LINES 9-13

Prepare TopographyMap
------------------------
First, we need to prepare the topography map. This is done by reading the
elevation and roughness maps and creating a TopographyMap object.

.. GENERATED FROM PYTHON SOURCE LINES 13-39

.. code-block:: Python



    import numpy as np
    import windkit as wk
    import pywasp as pw

    bwc = wk.read_bwc(
        "../../../modules/examples/tutorial_1/data/SerraSantaLuzia.omwc", crs="EPSG:4326"
    )
    bwc = wk.spatial.reproject(bwc, to_crs="EPSG:32629")

    elev_map = wk.read_vector_map(
        "../../../modules/examples/tutorial_1/data/SerraSantaLuzia.map",
        map_type="elevation",
        crs="EPSG:32629",
    )

    lc_map, lc_tbl = wk.read_vector_map(
        "../../../modules/examples/tutorial_1/data/SerraSantaLuzia.map",
        map_type="roughness",
        crs="EPSG:32629",
    )

    topo_map = pw.wasp.TopographyMap(elev_map, lc_map, lc_tbl)









.. GENERATED FROM PYTHON SOURCE LINES 40-44

Define output locations
------------------------
Second, we need to define the output locations. This is done by creating a
dataset with the coordinates of the output locations.

.. GENERATED FROM PYTHON SOURCE LINES 44-63

.. code-block:: Python


    dx, dy = 300, 300
    nx, ny = 51, 51
    x0, y0 = 510500, 4613500
    steps = np.arange(nx) * dx
    x = x0 + steps
    y = y0 + steps
    xmin, xmax = x.min(), x.max()
    ymin, ymax = y.min(), y.max()

    output_locs = wk.create_dataset(
        west_east=x,
        south_north=y,
        height=[50],
        crs="EPSG:32629",
        struct="cuboid",
    )









.. GENERATED FROM PYTHON SOURCE LINES 64-70

Calculate resource grid
------------------------
Finally, we can calculate the resource grid. This is done by calling the
predict_wwc function. This function takes the output locations,
the boundary conditions, and the topography map as input.
The output is a weibull wind climate dataset at the output locations.

.. GENERATED FROM PYTHON SOURCE LINES 70-74

.. code-block:: Python


    wwc = pw.wasp.predict_wwc(bwc, topo_map, output_locs)






.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    /opt/conda/lib/python3.11/site-packages/windkit/spatial/_bbox.py:263: UserWarning: 'use_bounds' currently defaults to False in BBox.reproject, in the future this will change to True
      warnings.warn(
    /opt/conda/lib/python3.11/site-packages/windkit/spatial/_bbox.py:263: UserWarning: 'use_bounds' currently defaults to False in BBox.reproject, in the future this will change to True
      warnings.warn(
    /opt/conda/lib/python3.11/site-packages/windkit/spatial/_bbox.py:263: UserWarning: 'use_bounds' currently defaults to False in BBox.reproject, in the future this will change to True
      warnings.warn(
    /opt/conda/lib/python3.11/site-packages/windkit/spatial/_bbox.py:263: UserWarning: 'use_bounds' currently defaults to False in BBox.reproject, in the future this will change to True
      warnings.warn(




.. GENERATED FROM PYTHON SOURCE LINES 75-78

Plot the mean wind speed
------------------------
We can plot the mean wind speed to see the result.

.. GENERATED FROM PYTHON SOURCE LINES 78-80

.. code-block:: Python


    wwc["wspd"].plot()



.. image-sg:: /auto_examples/images/sphx_glr_resource_grid_001.png
   :alt: height = 50 [m], crs = 0
   :srcset: /auto_examples/images/sphx_glr_resource_grid_001.png
   :class: sphx-glr-single-img






.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 11.437 seconds)


.. _sphx_glr_download_auto_examples_resource_grid.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: resource_grid.ipynb <resource_grid.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: resource_grid.py <resource_grid.py>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_