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

.. only:: html

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

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

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

.. _sphx_glr_examples_morpho_example_maxtree.py:


Maxtree example
===================================

This example shows how to use the component trees with Pylena

.. GENERATED FROM PYTHON SOURCE LINES 9-18

.. code-block:: Python


    import matplotlib.pyplot as plt
    from skimage.data import camera
    import pylena as pln
    import pandas as pd
    import numpy as np

    img = camera()








.. GENERATED FROM PYTHON SOURCE LINES 19-20

Compute the max-tree and the area of component and a custom attribute

.. GENERATED FROM PYTHON SOURCE LINES 20-41

.. code-block:: Python


    maxtree = pln.morpho.maxtree(img, 4)
    area = maxtree.compute_area()


    mass_centers = maxtree.compute_attribute(lambda init, p: (init[0] + p[0], init[1] + p[1]), (0, 0))
    mass_centers = np.array(list(mass_centers), ndmin=2)

    df = pd.DataFrame(
        {
            "parent": maxtree.parent,
            "gray": maxtree.values,
            "area": area,
            "mc_x": mass_centers[:, 0] / area,
            "mc_y": mass_centers[:, 1] / area,
        }
    )

    df







.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>parent</th>
          <th>gray</th>
          <th>area</th>
          <th>mc_x</th>
          <th>mc_y</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>0</th>
          <td>-1</td>
          <td>0</td>
          <td>262144</td>
          <td>255.500000</td>
          <td>255.500000</td>
        </tr>
        <tr>
          <th>1</th>
          <td>0</td>
          <td>1</td>
          <td>262143</td>
          <td>255.499498</td>
          <td>255.500525</td>
        </tr>
        <tr>
          <th>2</th>
          <td>1</td>
          <td>2</td>
          <td>262142</td>
          <td>255.498993</td>
          <td>255.501049</td>
        </tr>
        <tr>
          <th>3</th>
          <td>2</td>
          <td>3</td>
          <td>262122</td>
          <td>255.492103</td>
          <td>255.510003</td>
        </tr>
        <tr>
          <th>4</th>
          <td>3</td>
          <td>4</td>
          <td>261495</td>
          <td>255.318052</td>
          <td>255.834716</td>
        </tr>
        <tr>
          <th>...</th>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
          <td>...</td>
        </tr>
        <tr>
          <th>48994</th>
          <td>3</td>
          <td>4</td>
          <td>1</td>
          <td>328.000000</td>
          <td>119.000000</td>
        </tr>
        <tr>
          <th>48995</th>
          <td>3</td>
          <td>4</td>
          <td>1</td>
          <td>314.000000</td>
          <td>152.000000</td>
        </tr>
        <tr>
          <th>48996</th>
          <td>3</td>
          <td>4</td>
          <td>1</td>
          <td>308.000000</td>
          <td>152.000000</td>
        </tr>
        <tr>
          <th>48997</th>
          <td>3</td>
          <td>4</td>
          <td>3</td>
          <td>327.666667</td>
          <td>136.666667</td>
        </tr>
        <tr>
          <th>48998</th>
          <td>3</td>
          <td>4</td>
          <td>1</td>
          <td>329.000000</td>
          <td>135.000000</td>
        </tr>
      </tbody>
    </table>
    <p>48999 rows × 5 columns</p>
    </div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 42-43

Filtering and reconstruction

.. GENERATED FROM PYTHON SOURCE LINES 43-59

.. code-block:: Python


    maxtree.filter((df["area"] > 200).values, inplace=True)
    filtered = maxtree.reconstruct()

    fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(10, 5))
    axes[0].imshow(img, cmap="gray")
    axes[1].imshow(filtered, cmap="gray")
    axes[0].set_title("Original")
    axes[1].set_title("Filtered")
    for i in range(2):
        axes[i].set_xticks([])
        axes[i].set_yticks([])
    plt.tight_layout()

    plt.show()




.. image-sg:: /examples/morpho/images/sphx_glr_example_maxtree_001.png
   :alt: Original, Filtered
   :srcset: /examples/morpho/images/sphx_glr_example_maxtree_001.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_examples_morpho_example_maxtree.py:

.. only:: html

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

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

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

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

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

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: example_maxtree.zip <example_maxtree.zip>`


.. only:: html

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

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