Merging Slope
=============

Main functions
--------------

.. autofunction:: SignalCreation.Treatment.MergingSlope.merge

Specific merging functions
--------------------------

.. autofunction:: SignalCreation.Treatment.MergingSlope.linear_merging

.. autofunction:: SignalCreation.Treatment.MergingSlope.cossin_merging

Tools functions
---------------

.. autofunction:: SignalCreation.Treatment.MergingSlope.channel_shift

Test process
------------

.. autofunction:: SignalCreation.Treatment.MergingSlope.test_merge

Examples
--------

.. code-block:: python

    # Input files
    lidar_file = ["/home/lotto/Documents/osu_reunion/python/data/s2161614.044649",
                  "/home/lotto/Documents/osu_reunion/python/data/s2161614.054759"]
    atmosphere_file = "/home/lotto/Documents/osu_reunion/data/2021/06/arletty_210615.txt"
    cross_section_file = "/home/lotto/Documents/osu_reunion/python/data/SerdyuchenkoGorshelev5digits.dat"

    # Reading Methods
    lidar = ReadLidar(lidar_file, LICEL)
    atmosphere = Atmosphere(atmosphere_file, ARLETTY)
    ozone_cross_section = CrossSection(cross_section_file, SERDYUNCHENKO_VALUES_FILE)

    # Geoscene creation
    geoscene = GeoSceneLidar(lidar, atmosphere, ozone_cross_section)

    # Add Altitude Vector
    altitude_step = geoscene.output_dataset["Raw_Data_Range_Resolution"].data[0]
    altitude_vector = np.array(range(0, geoscene.output_dataset.dims["points"] * altitude_step, altitude_step))

    geoscene.output_dataset["Altitude"] = xr.DataArray(altitude_vector, dims=["points"])

    # Merge variables
    s1 = geoscene.output_dataset["Raw_Lidar_Data"][0]
    s2 = geoscene.output_dataset["Raw_Lidar_Data"][2]
    altitude = geoscene.output_dataset["Altitude"].data
    alt_min = 8000  # meter
    alt_max = 12000  # meter
    method = "linear"
    shift = 1

    # Merging...
    merge(s1, s2, altitude, alt_min, alt_max, method, shifted_signal=shift)