Add online documentation for morphpy

diffpy · sbillinge · Jun 27, 2025 · Apr 29, 2025 · Apr 29, 2025 · Apr 29, 2025
commit 80b4a603ee23972e222bff4ed220d6f30cf10f55
diff --git a/doc/source/morphpy.rst b/doc/source/morphpy.rst
@@ -0,0 +1,173 @@
+.. _morphpy:
+
+Using diffpy.morph in Python
+############################
+
+On top of the command-line usage described in the `quickstart tutorial <quickstart.html>`__,
+``diffpy.morph`` also supports Python integration.
+This functionality is intended for those acquainted with the basic morphs
+described in the aforementioned quickstart tutorial who want to use ``diffpy.morph`` in their
+Python scripts.
+
+Python Morphing Functions
+=========================
+
+    1. In the quickstart tutorial, you were asked to try a combined scale, stretch, and smear
+       morph on the files `darkSub_rh20_C_01.gr` and `darkSub_rh20_C_44.gr` using the command-line
+       command ::
+
+          diffpy.morph --scale=0.8 --smear=-0.08 --stretch=0.5 --rmin=1.5 --rmax=30 darkSub_rh20_C_01.gr darkSub_rh20_C_44.gr
+
+    2. To do the same on Python, we must first create a new Python script in the same directory as the
+       data files `darkSub_rh20_C_01.gr` and `darkSub_rh20_C_44.gr`.
+    3. Then, in that script, import ::
+
+           from diffpy.morph.morphpy import morph
+
+    3. Finally, we run the ``morph`` function ::
+
+           morph_info, morph_table = morph("darkSub_rh20_C_01.gr", "darkSub_rh20_C_44.gr", scale=0.8, smear=-0.08, stretch=0.5, rmin=1.5, rmax=30)
+
+           * The ``morph`` function takes in two file names (or paths). You can also provide various parameters
+             for morphing (see the Full Parameter List below).
+           * If, let's say, the file `darkSub_rh20_C_01.gr` is in a subdirectory `subdir/darkSub_rh20_C_01.gr`,
+             you should replace ``"darkSub_rh20_C_01.gr"`` in the above example with ``"subdir/darkSub_rh20_C_01.gr"``.
+
+    4. The ``morph`` function returns a dictionary ``morph_info`` and a numpy array ``morph_table``.
+
+       * ``morph_info`` contains all morphs as keys (e.g. ``"scale"``, ``"stretch"``, ``"smear"``) with
+         the optimized morphing parameters found by ``diffpy.morph`` as values. ``morph_info`` also contains
+         the Rw and Pearson correlation coefficients found post-morphing.
+       * ``morph_table`` is a two-column array of the morphed function interpolated onto the grid of the
+         target function (e.g. in our example, it returns the contents of `darkSub_rh20_C_01.gr` after
+         the morphs are applied interpolated onto the grid of `darkSub_rh20_C_44.gr`).
+    5. Notice that most parameters you are able to use are the same as the options provided in the command-line
+       interface version of ``diffpy.morph``. For example, the ``--apply`` option becomes the ``apply=True`` parameter.
+    6. With that, you have already mastered the basics of using ``diffpy.morph`` on Python!
+    7. Note that instead of passing two files to ``diffpy.morph``, you might instead want to directly
+       pass arrays. For example, rather than passing `darkSub_rh20_C_01.gr`, I may want to pass
+       a two-column array named ``ds_rh20_c_01_array`` containing the data table contents of the file
+       `darkSub_rh20_C_01.gr`. In this case, we have a separate function ::
+
+           from diffpy.morph.morphpy import morph_arrays
+
+    8. Assuming we have loaded the data in `darkSub_rh20_C_01.gr` into ``ds_rh20_c_01_array`` and
+       `darkSub_rh20_C_44.gr` into ``ds_rh20_c_44_array``, we can apply the same morph as step 3
+       by running ::
+
+           morph_info, morph_table = morph_arrays(ds_rh20_c_01_array, ds_rh20_c_44_array, scale=0.8, smear=-0.08, stretch=0.5, rmin=1.5, rmax=30)
+
+    9. Notice that the two-column format of the input to ``morph_arrays`` is the same as the
+       output of ``morph`` and ``morph_arrays``. It is VERY IMPORTANT that the data is in two-column format
+       rather than the traditional two-row format. This is to reflect the file format used to store PDFs.
+    10. For a full list of parameters used by (both) ``morph`` and ``morph_arrays``, see the Full Parameter List
+        section below.
+
+Full Parameter List
+===================
+
+General Parameters
+------------------
+
+save: str or path
+    Save the morphed function to a the file passed to save. Use '-' for stdout.
+verbose: bool
+    Print additional header details to saved files. These include details about the morph
+    inputs and outputs.
+rmin: float
+    Minimum r-value (abscissa) to use for function comparisons.
+rmax: float
+    Maximum r-value (abscissa) to use for function comparisons.
+tolerance: float
+    Specify least squares refiner tolerance when optimizing for morph parameters. Default: 10e-8.
+pearson: bool
+    The refiner instead maximizes agreement in the Pearson function
+    (default behavior is to minimize the residual).
+    Note that this is insensitive to scale.
+addpearson: bool
+    Maximize agreement in the Pearson function as well as minimizing the residual.
+
+Manipulations
+-------------
+These parameters select the manipulations that are to be applied to the
+function. The passed values will be refined unless specifically
+excluded with the apply or exclude parameters.
+
+apply: bool
+    Apply morphs but do not refine.
+exclude: str
+    Exclude a manipulation from refinement by name.
+scale: float
+    Apply scale factor. This multiplies the function ordinate by scale.
+stretch: float
+    Stretch function grid by a fraction stretch. Specifically, this multiplies the function grid by 1+stretch.
+squeeze: list of float
+    Squeeze function grid given a polynomial
+    p(x) = squeeze[0]+squeeze[1]*x+...+squeeze[n]*x^n. n is dependent on the number
+    of values in the user-inputted comma-separated list.
+    The morph transforms the function grid from x to x+p(x).
+    When this parameter is given, hshift is disabled.
+    When n>1, stretch is disabled.
+smear: float
+    Smear the peaks with a Gaussian of width smear. This
+    is done by convolving the function with a Gaussian
+    with standard deviation smear. If both smear and
+    smear_pdf are used, only smear_pdf will be
+    applied.
+smear_pdf: float
+    Convert PDF to RDF. Then, smear peaks with a Gaussian
+    of width smear_pdf. Convert back to PDF. If both smear and
+    smear_pdf are used, only smear_pdf will be
+    applied.
+slope: float
+    Slope of the baseline used in converting from PDF to RDF.
+    This is used with the option smear_pdf. The slope will
+    be estimated if not provided.
+hshift: float
+    Shift the function horizontally by hshift to the right.
+vshift: float
+    Shift the function vertically by vshift upward.
+qdamp: float
+    Dampen PDF by a factor qdamp.
+radius: float
+    Apply characteristic function of sphere with radius
+    given by parameter radius. If pradius is also specified, instead apply
+    characteristic function of spheroid with equatorial
+    radius radius and polar radius pradius.
+pradius: float
+    Apply characteristic function of spheroid with
+    equatorial radius given by above parameter radius and polar radius pradius.
+    If only pradius is specified, instead apply
+    characteristic function of sphere with radius pradius.
+iradius: float
+    Apply inverse characteristic function of sphere with
+    radius iradius. If ipradius is also specified, instead
+    apply inverse characteristic function of spheroid with
+    equatorial radius iradius and polar radius ipradius.
+ipradius: float
+    Apply inverse characteristic function of spheroid with
+    equatorial radius iradius and polar radius ipradius.
+    If only ipradius is specified, instead apply inverse
+    characteristic function of sphere with radius ipradius.
+funcy: tuple (function, dict)
+    See Python-Specific Morphs below.
+
+Python-Specific Morphs
+======================
+
+Some morphs in diffpy.morph are supported only in Python. Here, we detail
+how they are used and how to call them.
+
+funcy: tuple (function, dict)
+    This morph applies the function funcy[0] with parameters given in funcy[1].
+    The function funcy[0] must be a function of both the abscissa and ordinate
+    (e.g. take in at least two inputs with as many additional parameters as needed).
+    For example, let's start with a two-column table with abscissa x and ordinate y.
+    let us say we want to apply the function ::
+
+        def linear(x, y, a, b, c):
+            return a * x + b * y + c
+
+    This function takes in both the abscissa and ordinate on top of three additional
+    parameters a, b, and c. To use the funcy parameter with initial guesses
+    a=1.0, b=2.0, c=3.0, we would pass ``funcy=(linear, {a: 1.0, b: 2.0, c: 3.0})``.
diff --git a/doc/source/quickstart.rst b/doc/source/quickstart.rst
@@ -4,7 +4,11 @@ diffpy.morph Tutorial
 #####################
 
 Welcome! This will be a quick tutorial to accquaint users with ``diffpy.morph``
-and some of what it can do. To see more details and definitions about
+and some of what it can do on the command-line.
+For those wishing to integrate ``diffpy.morph`` into their Python scripts,
+see the `morphpy tutorial <morphpy.html>`__.
+
+To see more details and definitions about
 the morphs please see the publication describing ``diffpy.morph``.
 
 To be published:

diff --git a/src/diffpy/morph/morphpy.py b/src/diffpy/morph/morphpy.py
@@ -67,6 +67,15 @@ def morph(
             pmorph_value = kwargs.pop(pmorph)
             pymorphs.update({pmorph: pmorph_value})
 
+    # Special handling of parameters with dashes
+    kwargs_copy = kwargs.copy()
+    kwargs = {}
+    for key in kwargs_copy.keys():
+        new_key = key
+        if "_" in key:
+            new_key = key.replace("_", "-")
+        kwargs.update({new_key: kwargs_copy[key]})
+
     # Wrap the CLI
     parser = create_option_parser()
     params = {
@@ -92,7 +101,7 @@ def morph(
 
 
 # Take in array-like objects as input.
-def morphpy(
+def morph_arrays(
     morph_table,
     target_table,
     scale=None,

diff --git a/tests/test_morphio.py b/tests/test_morphio.py
@@ -10,7 +10,7 @@
     multiple_targets,
     single_morph,
 )
-from diffpy.morph.morphpy import morphpy
+from diffpy.morph.morphpy import morph_arrays
 
 # Support Python 2
 try:
@@ -201,7 +201,7 @@ def quadratic(x, y, a0, a1, a2):
         r = np.linspace(0, 10, 101)
         gr = np.linspace(0, 10, 101)
 
-        morphpy(
+        morph_arrays(
             np.array([r, gr]).T,
             np.array([r, quadratic(r, gr, 1, 2, 3)]).T,
             squeeze=[0, 0, 0],

diff --git a/tests/test_morphpy.py b/tests/test_morphpy.py
@@ -6,7 +6,7 @@
 import pytest
 
 from diffpy.morph.morphapp import create_option_parser, single_morph
-from diffpy.morph.morphpy import morph, morphpy
+from diffpy.morph.morphpy import morph, morph_arrays
 from diffpy.morph.tools import getRw
 
 thisfile = locals().get("__file__", "file.py")
@@ -99,7 +99,7 @@
         for target_file in self.testfiles[1:]:
             _, grm0 = morph(morph_file, morph_file)
             _, grt = morph(target_file, target_file)
-            mr, grm = morphpy(
+            mr, grm = morph_arrays(
                 grm0, grt, scale=1, stretch=0, sort_by="temperature"
             )
             morph_results.update({target_file.name: mr})
@@ -130,12 +130,12 @@
         target_r = np.linspace(0, 100, 1001)
         target_gr = 0.5 * gaussian(target_r, 50, 5) + 0.05
 
-        morph_info, _ = morphpy(
+        morph_info, _ = morph_arrays(
             np.array([morph_r, morph_gr]).T,
             np.array([target_r, target_gr]).T,
             scale=1,
-            vshift=0.01,
             smear=3.75,
+            vshift=0.01,
             funcy=(gaussian_like_function, {"mu": 47.5}),
             tolerance=1e-12,
         )
@@ -152,7 +152,7 @@
         def linear(x, y, s):
             return s * (x + y)
 
-        morph_info, _ = morphpy(
+        morph_info, _ = morph_arrays(
             np.array([r, gr]).T,
             np.array([r, gr]).T,
             squeeze=[1, 2, 3, 4, 5],
@@ -167,4 +167,4 @@


 if __name__ == "__main__":
    TestMorphpy()