# HG changeset patch
# User imgteam
# Date 1710240880 0
# Node ID 214f548481eba37ad3d824d3e370ae5862bb8b84

planemo upload for repository https://github.com/BMCV/galaxy-image-analysis/tree/master/tools/orientationpy commit fa4501fb81c17f0b56889f250cf92396804295d1

diff -r 000000000000 -r 214f548481eb orientationpy-cli.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/orientationpy-cli.py	Tue Mar 12 10:54:40 2024 +0000
@@ -0,0 +1,54 @@
+import argparse
+import csv
+
+import numpy as np
+import orientationpy
+import skimage.io
+import skimage.util
+
+
+if __name__ == '__main__':
+
+    parser = argparse.ArgumentParser()
+    parser.add_argument('input', type=str)
+    parser.add_argument('--mode', type=str, required=True)
+    parser.add_argument('--sigma', type=float, required=True)
+    parser.add_argument('--min_coherency', type=float, required=True)
+    parser.add_argument('--min_energy', type=float, required=True)
+    parser.add_argument('--max_precision', type=int, required=True)
+    parser.add_argument('--output_angle_tsv', type=str, default=None)
+    args = parser.parse_args()
+
+    im = skimage.io.imread(args.input)
+    im = skimage.util.img_as_float(im)
+    im = np.squeeze(im)
+    assert im.ndim == 2
+
+    Gy, Gx = orientationpy.computeGradient(im, mode=args.mode)
+    structureTensor = orientationpy.computeStructureTensor([Gy, Gx], sigma=args.sigma)
+    orientations = orientationpy.computeOrientation(structureTensor, computeEnergy=True, computeCoherency=True)
+
+    # Compute angle according to https://bigwww.epfl.ch/demo/orientationj/#dist:
+    mask = np.logical_and(
+        orientations['coherency'] >= args.min_coherency,
+        orientations['energy'] >= args.min_energy * orientations['energy'].max(),
+    )
+    angles = orientations['theta'][mask]
+    weights = orientations['coherency'][mask]
+    bin_size = 1 if args.max_precision == 0 else pow(10, -args.max_precision)
+    hist, bin_edges = np.histogram(
+        angles,
+        range=(-90, +90),
+        weights=weights,
+        bins=round(180 / bin_size),
+    )
+    hidx = np.argmax(hist)
+    angle = (bin_edges[hidx] + bin_edges[hidx + 1]) / 2
+    angle = round(angle, args.max_precision)
+
+    # Write results
+    if args.output_angle_tsv:
+        with open(args.output_angle_tsv, 'w') as fp:
+            writer = csv.writer(fp, delimiter='\t', lineterminator='\n')
+            writer.writerow(['Angle'])
+            writer.writerow([angle])
diff -r 000000000000 -r 214f548481eb orientationpy.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/orientationpy.xml	Tue Mar 12 10:54:40 2024 +0000
@@ -0,0 +1,97 @@
+<tool id="orientationpy" name="Compute image orientation" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="23.0">
+    <description>with OrientationPy</description>
+    <macros>
+        <token name="@TOOL_VERSION@">0.2.0.4</token>
+        <token name="@VERSION_SUFFIX@">0</token>
+    </macros>
+    <edam_operations>
+        <edam_operation>operation_3443</edam_operation>
+    </edam_operations>
+    <xrefs>
+        <xref type="biii">orientationj</xref>
+    </xrefs>
+    <requirements>
+        <requirement type="package" version="0.2.0.4">orientationpy</requirement>
+        <requirement type="package" version="0.22.0">scikit-image</requirement>
+    </requirements>
+    <command><![CDATA[
+    
+    ## Inputs
+    
+    python '$__tool_directory__/orientationpy-cli.py' '$image'
+    --mode $mode
+    --sigma $sigma
+    --min_energy $min_energy
+    --min_coherency $min_coherency
+    --max_precision $max_precision
+    
+    ## Outputs
+    
+    --output_angle_tsv '$output_angle_tsv'
+    
+    ]]>
+    </command>
+    <inputs>
+        <param name="image" type="data" format="png,tiff" label="Input image" />
+        <param argument="--mode" type="select" label="Gradient computation">
+            <option value="finite_difference" selected="True">Finite differences</option>
+            <option value="gaussian">Gaussian</option>
+            <option value="splines">Splines</option>
+        </param>
+        <param argument="--sigma" type="float" label="Spatial scale (in pixels)" value="2" />
+        <param argument="--min_energy" type="float" label="Minimum energy" value="0.1" min="0" max="1" />
+        <param argument="--min_coherency" type="float" label="Minimum coherency" value="0.2" min="0" max="1" />
+        <param argument="--max_precision" type="integer" label="Maximum possible precision (in decimals)" value="3" min="0" max="8" help="Governs the histogram bin size (increasing maximum possible precision leads to smaller bin sizes). Note that increasing the maximum possible precision too far may actually decrease the actual precision, because the bin size becomes too small." />
+    </inputs>
+    <outputs>
+        <data format="tabular" name="output_angle_tsv" from_work_dir="output_angle.tsv" />
+    </outputs>
+    <tests>
+        <!-- Test negative angle -->
+        <test>
+            <param name="image" value="input1.tif" />
+            <param name="mode" value="finite_difference" />
+            <param name="sigma" value="2" />
+            <param name="min_energy" value="0.1" />
+            <param name="min_coherency" value="0.2" />
+            <param name="max_precision" value="1" />
+            <output name="output_angle_tsv">
+                <assert_contents>
+                    <has_line line="Angle" />
+                    <has_line line="-45.0" />
+                    <has_n_lines n="2" />
+                    <has_n_columns n="1" />
+                </assert_contents>
+            </output>
+        </test>
+        <!-- Test positive angle -->
+        <test>
+            <param name="image" value="input2.tif" />
+            <param name="mode" value="finite_difference" />
+            <param name="sigma" value="2" />
+            <param name="min_energy" value="0.1" />
+            <param name="min_coherency" value="0.2" />
+            <param name="max_precision" value="1" />
+            <output name="output_angle_tsv">
+                <assert_contents>
+                    <has_line line="Angle" />
+                    <has_line line="45.0" />
+                    <has_n_lines n="2" />
+                    <has_n_columns n="1" />
+                </assert_contents>
+            </output>
+        </test>
+    </tests>
+    <help>
+        Compute the orientation of 2-D images in degrees using OrientationPy.
+
+        OrientationPy is the pythonic successor to the well-loved OrientationJ Fiji Plugin.
+        It is a library that takes in images and computes the orientation of the greylevels.
+        A key step is the computation of image gradients, for a number of different techniques is supported.
+
+        For more information on ``--min_energy`` and ``--min_coherency`` see: https://epfl-center-for-imaging.gitlab.io/orientationpy/orientationpy_examples/plot_fibres_2d.html#plot-hsv-composition
+    </help>
+    <citations>
+        <citation type="doi">10.1007/s10237-011-0325-z</citation>
+    </citations>
+</tool>
diff -r 000000000000 -r 214f548481eb test-data/input1.tif
Binary file test-data/input1.tif has changed
diff -r 000000000000 -r 214f548481eb test-data/input2.tif
Binary file test-data/input2.tif has changed