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

.. only:: html

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

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

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

.. _sphx_glr_gallery_custom_materials_plot_measured_material.py:


Create and export a measured material
=====================================

This example builds a complex refractive-index model from measured arrays,
evaluates the interpolation, and verifies a YAML export/reload round trip.

.. GENERATED FROM PYTHON SOURCE LINES 10-33

.. code-block:: Python

    from tempfile import TemporaryDirectory

    import matplotlib.pyplot as plt
    import numpy
    from TypedUnit import ureg

    from PyOptik import TabulatedMaterial


    wavelength_samples = numpy.array([400, 475, 550, 625, 700]) * ureg.nanometer
    measured_n = [1.58, 1.56, 1.545, 1.535, 1.528]
    measured_k = [0.030, 0.018, 0.010, 0.006, 0.004]

    material = TabulatedMaterial.from_arrays(
        "measured-film",
        wavelength_samples,
        n=measured_n,
        k=measured_k,
        interpolation="pchip",
        reference="Illustrative ellipsometry dataset",
        conditions={"temperature": "293 K"},
    )








.. GENERATED FROM PYTHON SOURCE LINES 34-35

Exported files use the same validated schema as catalog materials.

.. GENERATED FROM PYTHON SOURCE LINES 35-40

.. code-block:: Python

    with TemporaryDirectory() as directory:
        path = material.to_yaml(f"{directory}/measured-film.yml")
        restored = TabulatedMaterial("measured-film", file_path=path, interpolation="pchip")
        assert restored.reference == material.reference








.. GENERATED FROM PYTHON SOURCE LINES 41-55

.. code-block:: Python

    wavelengths = numpy.linspace(400, 700, 301) * ureg.nanometer
    index = material.compute_refractive_index(wavelengths, out_of_range="raise")

    figure, (axis_n, axis_k) = plt.subplots(2, 1, sharex=True, layout="constrained")
    axis_n.plot(wavelengths.magnitude, index.real, label="PCHIP interpolation")
    axis_n.scatter(wavelength_samples.magnitude, measured_n, label="Measurements")
    axis_n.set(ylabel="Refractive index n", title="User-defined optical constants")
    axis_n.legend()

    axis_k.plot(wavelengths.magnitude, index.imag, color="tab:red")
    axis_k.scatter(wavelength_samples.magnitude, measured_k, color="tab:red")
    axis_k.set(xlabel="Wavelength [nm]", ylabel="Extinction coefficient k")

    plt.show()



.. image-sg:: /gallery/custom_materials/images/sphx_glr_plot_measured_material_001.png
   :alt: User-defined optical constants
   :srcset: /gallery/custom_materials/images/sphx_glr_plot_measured_material_001.png
   :class: sphx-glr-single-img






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

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


.. _sphx_glr_download_gallery_custom_materials_plot_measured_material.py:

.. only:: html

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

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

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

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

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


.. only:: html

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

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