# HG changeset patch
# User imgteam
# Date 1747058469 0
# Node ID 64c155acb864432f6a5d2d3a52936ae2bdc1d318
# Parent  2ae122d5d85a68c754e5a34b0c507f6a52ad0085
planemo upload for repository https://github.com/BMCV/galaxy-image-analysis/tree/master/tools/points2labelimage/ commit 9a40a5d1e1008c26cc327c6d163df2a1af22a1a0

diff -r 2ae122d5d85a -r 64c155acb864 creators.xml
--- a/creators.xml	Wed Apr 23 14:37:31 2025 +0000
+++ b/creators.xml	Mon May 12 14:01:09 2025 +0000
@@ -24,5 +24,10 @@
         <person givenName="Till" familyName="Korten"/>
         <yield/>
     </xml>
-    
+
+    <xml name="creators/pavanvidem">
+        <person givenName="Pavan" familyName="Videm"/>
+        <yield/>
+    </xml>
+
 </macros>
diff -r 2ae122d5d85a -r 64c155acb864 points2label.py
--- a/points2label.py	Wed Apr 23 14:37:31 2025 +0000
+++ b/points2label.py	Mon May 12 14:01:09 2025 +0000
@@ -1,15 +1,82 @@
 import argparse
+import json
 import os
 import warnings
+from typing import (
+    Dict,
+    List,
+    Tuple,
+    Union,
+)
 
 import giatools.pandas
 import numpy as np
+import numpy.typing as npt
 import pandas as pd
 import scipy.ndimage as ndi
 import skimage.io
 import skimage.segmentation
 
 
+def is_rectangular(points: Union[List[Tuple[float, float]], npt.NDArray]) -> bool:
+    points = np.asarray(points)
+
+    # Rectangle must have 5 points, where first and last are identical
+    if len(points) != 5 or not (points[0] == points[-1]).all():
+        return False
+
+    # Check that all edges align with the axes
+    edges = points[1:] - points[:-1]
+    if any((edge == 0).sum() != 1 for edge in edges):
+        return False
+
+    # All checks have passed, the geometry is rectangular
+    return True
+
+
+def geojson_to_tabular(geojson: Dict):
+    rows = []
+    labels = []
+    for feature in geojson['features']:
+        assert feature['geometry']['type'].lower() == 'polygon', (
+            f'Unsupported geometry type: "{feature["geometry"]["type"]}"'
+        )
+        coords = feature['geometry']['coordinates'][0]
+
+        # Properties and name (label) are optional
+        try:
+            label = feature['properties']['name']
+        except KeyError:
+            label = max(labels, default=0) + 1
+        labels.append(label)
+
+        # Read geometry
+        xs = [pt[0] for pt in coords]
+        ys = [pt[1] for pt in coords]
+
+        x = min(xs)
+        y = min(ys)
+
+        width = max(xs) + 1 - x
+        height = max(ys) + 1 - y
+
+        # Validate geometry (must be rectangular)
+        assert is_rectangular(list(zip(xs, ys)))
+
+        # Append the rectangle
+        rows.append({
+            'pos_x': x,
+            'pos_y': y,
+            'width': width,
+            'height': height,
+            'label': label,
+        })
+    df = pd.DataFrame(rows)
+    point_file = './point_file.tabular'
+    df.to_csv(point_file, sep='\t', index=False)
+    return point_file
+
+
 def rasterize(point_file, out_file, shape, has_header=False, swap_xy=False, bg_value=0, fg_value=None):
 
     img = np.full(shape, dtype=np.uint16, fill_value=bg_value)
@@ -122,16 +189,16 @@
                 img[y, x] = label
 
     else:
-        raise Exception("{} is empty or does not exist.".format(point_file))  # appropriate built-in error?
+        raise Exception('{} is empty or does not exist.'.format(point_file))  # appropriate built-in error?
 
     with warnings.catch_warnings():
         warnings.simplefilter("ignore")
         skimage.io.imsave(out_file, img, plugin='tifffile')  # otherwise we get problems with the .dat extension
 
 
-if __name__ == "__main__":
+if __name__ == '__main__':
     parser = argparse.ArgumentParser()
-    parser.add_argument('point_file', type=argparse.FileType('r'), help='point file')
+    parser.add_argument('in_file', type=argparse.FileType('r'), help='Input point file or GeoJSON file')
     parser.add_argument('out_file', type=str, help='out file (TIFF)')
     parser.add_argument('shapex', type=int, help='shapex')
     parser.add_argument('shapey', type=int, help='shapey')
@@ -141,11 +208,25 @@
 
     args = parser.parse_args()
 
+    point_file = args.in_file.name
+    has_header = args.has_header
+
+    try:
+        with open(args.in_file.name, 'r') as f:
+            content = json.load(f)
+            if isinstance(content, dict) and content.get('type') == 'FeatureCollection' and isinstance(content.get('features'), list):
+                point_file = geojson_to_tabular(content)
+                has_header = True  # header included in the converted file
+            else:
+                raise ValueError('Input is a JSON file but not a valid GeoJSON file')
+    except json.JSONDecodeError:
+        print('Input is not a valid JSON file. Assuming it a tabular file.')
+
     rasterize(
-        args.point_file.name,
+        point_file,
         args.out_file,
         (args.shapey, args.shapex),
-        has_header=args.has_header,
+        has_header=has_header,
         swap_xy=args.swap_xy,
         fg_value=0xffff if args.binary else None,
     )
diff -r 2ae122d5d85a -r 64c155acb864 points2label.xml
--- a/points2label.xml	Wed Apr 23 14:37:31 2025 +0000
+++ b/points2label.xml	Mon May 12 14:01:09 2025 +0000
@@ -1,10 +1,10 @@
-<tool id="ip_points_to_label" name="Convert point coordinates to label map" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="20.05">
+<tool id="ip_points_to_label" name="Convert coordinates to label map" version="@TOOL_VERSION@+galaxy@VERSION_SUFFIX@" profile="20.05">
     <description></description>
     <macros>
         <import>creators.xml</import>
         <import>tests.xml</import>
         <token name="@TOOL_VERSION@">0.4.1</token>
-        <token name="@VERSION_SUFFIX@">0</token>
+        <token name="@VERSION_SUFFIX@">1</token>
     </macros>
     <creator>
         <expand macro="creators/bmcv" />
@@ -15,12 +15,12 @@
     <xrefs>
         <xref type="bio.tools">galaxy_image_analysis</xref>
     </xrefs>
-    <requirements> 
-        <requirement type="package" version="0.21">scikit-image</requirement> 
+    <requirements>
+        <requirement type="package" version="0.21">scikit-image</requirement>
         <requirement type="package" version="1.26.4">numpy</requirement>
         <requirement type="package" version="1.2.4">pandas</requirement>
         <requirement type="package" version="2024.6.18">tifffile</requirement>
-        <requirement type="package" version="0.3.1">giatools</requirement> 
+        <requirement type="package" version="0.3.1">giatools</requirement>
     </requirements>
     <command detect_errors="aggressive"><![CDATA[
 
@@ -34,13 +34,13 @@
         $binary
 
     ]]></command>
-    <inputs> 
-        <param name="input" type="data" format="tabular" label="Tabular list of points"/> 
+    <inputs>
+        <param name="input" type="data" format="tabular,geojson" label="List of points in tabular or geojson format"/>
         <param name="shapex" type="integer" value="500" min="1" label="Width of output image" />
         <param name="shapey" type="integer" value="500" min="1" label="Height of output image" />
-        <param name="has_header" type="boolean" checked="true" truevalue="--has_header True" falsevalue="" optional="true" label="Tabular list of points has header" help="Turning this off will ignore the first row and assume that the X and Y coordinates correspond to the first and second column, respectively." /> 
+        <param name="has_header" type="boolean" checked="true" truevalue="--has_header True" falsevalue="" optional="true" label="Tabular list of points has header" help="Turning this off will ignore the first row and assume that the X and Y coordinates correspond to the first and second column, respectively. Ignored, if GeoJSON is used for input." />
         <param name="swap_xy" type="boolean" checked="false" falsevalue="" truevalue="--swap_xy True" optional="true" label="Swap X and Y coordinates" help="Swap the X and Y coordinates, regardless of whether the tabular list has a header or not." />
-        <param name="binary" type="boolean" checked="false" truevalue="--binary True" falsevalue="" optional="true" label="Produce binary image" help="Use the same label for all points (65535)." /> 
+        <param name="binary" type="boolean" checked="false" truevalue="--binary True" falsevalue="" optional="true" label="Produce binary image" help="Use the same label for all points (65535)." />
     </inputs>
     <outputs>
         <data name="output" format="tiff" />
@@ -54,7 +54,7 @@
             <param name="has_header" value="false" />
             <param name="swap_xy" value="true" />
             <param name="binary" value="true" />
-            <expand macro="tests/binary_image_diff" name="output" value="output1_binary.tif" ftype="tiff" />
+            <expand macro="tests/binary_image_diff" name="output" value="output1_binary.tiff" ftype="tiff" />
         </test>
         <!-- Binary / TSV with header / Circles -->
         <test>
@@ -64,7 +64,7 @@
             <param name="has_header" value="true" />
             <param name="swap_xy" value="false" />
             <param name="binary" value="true" />
-            <expand macro="tests/binary_image_diff" name="output" value="output2_binary.tif" ftype="tiff" />
+            <expand macro="tests/binary_image_diff" name="output" value="output2_binary.tiff" ftype="tiff" />
         </test>
         <!-- Labeled / TSV with header / Circles -->
         <test>
@@ -74,7 +74,7 @@
             <param name="has_header" value="true" />
             <param name="swap_xy" value="false" />
             <param name="binary" value="false" />
-            <expand macro="tests/label_image_diff" name="output" value="output2.tif" ftype="tiff" />
+            <expand macro="tests/label_image_diff" name="output" value="output2.tiff" ftype="tiff" />
         </test>
         <!-- Binary / TSV with header / TSV with labels / Circles -->
         <test>
@@ -84,7 +84,7 @@
             <param name="has_header" value="true" />
             <param name="swap_xy" value="false" />
             <param name="binary" value="true" />
-            <expand macro="tests/binary_image_diff" name="output" value="output3_binary.tif" ftype="tiff" />
+            <expand macro="tests/binary_image_diff" name="output" value="output3_binary.tiff" ftype="tiff" />
         </test>
         <!-- Labeled / TSV with header / TSV with labels / Circles -->
         <test>
@@ -94,7 +94,7 @@
             <param name="has_header" value="true" />
             <param name="swap_xy" value="false" />
             <param name="binary" value="false" />
-            <expand macro="tests/label_image_diff" name="output" value="output3.tif" ftype="tiff" />
+            <expand macro="tests/label_image_diff" name="output" value="output3.tiff" ftype="tiff" />
         </test>
         <!-- Labeled / TSV with header / TSV with labels / Rectangles -->
         <test>
@@ -104,20 +104,66 @@
             <param name="has_header" value="true" />
             <param name="swap_xy" value="false" />
             <param name="binary" value="false" />
-            <expand macro="tests/label_image_diff" name="output" value="output4.tif" ftype="tiff" />
+            <expand macro="tests/label_image_diff" name="output" value="output4.tiff" ftype="tiff" />
+        </test>
+        <!-- Labeled / ROI / Rectangles -->
+        <test>
+            <param name="input" value="rois.geojson" />
+            <param name="shapex" value="300" />
+            <param name="shapey" value="300" />
+            <param name="has_header" value="false" />
+            <param name="swap_xy" value="false" />
+            <param name="binary" value="false" />
+            <expand macro="tests/label_image_diff" name="output" value="output5.tiff" ftype="tiff" />
+        </test>
+        <test>
+            <param name="input" value="rois-noname.geojson" />
+            <param name="shapex" value="300" />
+            <param name="shapey" value="300" />
+            <param name="has_header" value="false" />
+            <param name="swap_xy" value="false" />
+            <param name="binary" value="false" />
+            <expand macro="tests/label_image_diff" name="output" value="output5.tiff" ftype="tiff" />
+        </test>
+        <test>
+            <param name="input" value="rois.geojson" />
+            <param name="shapex" value="300" />
+            <param name="shapey" value="300" />
+            <param name="has_header" value="false" />
+            <param name="swap_xy" value="true" />
+            <param name="binary" value="false" />
+            <expand macro="tests/label_image_diff" name="output" value="output6.tiff" ftype="tiff" />
+        </test>
+        <test expect_failure="true">
+            <!-- The test should fail because the GeoJSON is not rectangular -->
+            <param name="input" value="rois-illegal1.geojson" />
+            <param name="shapex" value="300" />
+            <param name="shapey" value="300" />
+            <param name="has_header" value="false" />
+            <param name="swap_xy" value="false" />
+            <param name="binary" value="false" />
+        </test>
+        <test expect_failure="true">
+            <!-- The test should fail because the GeoJSON is not of `Polygon` type -->
+            <param name="input" value="rois-illegal2.geojson" />
+            <param name="shapex" value="300" />
+            <param name="shapey" value="300" />
+            <param name="has_header" value="false" />
+            <param name="swap_xy" value="false" />
+            <param name="binary" value="false" />
         </test>
     </tests>
     <help>
 
-        **Converts a tabular list of points to a label map by rasterizing the point coordinates.**
+        **Converts a list of points to a label map by rasterizing the coordinates.**
 
         The created image is a single-channel image with 16 bits per pixel (unsigned integer). The points are
         rasterized with unique labels, or the value 65535 (white) for binary image output. Pixels not corresponding to
         any points in the tabular file are assigned the value 0 (black).
 
-        The tabular list of points can either be header-less. In this case, the first and second columns are expected
-        to be the X and Y coordinates, respectively. Otherwise, if a header is present, it is searched for the
-        following column names:
+        **Using a tabular input file:** The tabular list of points can either be header-less. In this case, the first
+        and second columns are expected to be the X and Y coordinates, respectively. Otherwise, if a header is present,
+        it is searched for the following column names:
 
         - ``pos_x`` or ``POS_X``: This column corresponds to the X coordinates.
         - ``pos_y`` or ``POS_Y``: This column corresponds to the Y coordinates.
@@ -126,10 +172,13 @@
         - If ``width`` or ``WIDTH`` and ``height`` or ``HEIGHT`` columns are present, then the points will be rasterized
           as rectangles of the corresponding size.
         - If a ``label`` or ``LABEL`` column is present, then the corresponding labels will be used for rasterization
-          (unless "Produce binary image" is activated). Different points are allowed to use the same label.
+          (unless "Produce binary image" is activated). Different points are allowed to use the same label. If used, the
+          label must be numeric and integer.
+
+        **Using a GeoJSON input file:** Only rectangular specifications of `Polygon` type geometry is supported.
 
     </help>
     <citations>
-        <citation type="doi">10.1016/j.jbiotec.2017.07.019</citation> 
+        <citation type="doi">10.1016/j.jbiotec.2017.07.019</citation>
     </citations>
 </tool>
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output1_binary.tif
Binary file test-data/output1_binary.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output1_binary.tiff
Binary file test-data/output1_binary.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output2.tif
Binary file test-data/output2.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output2.tiff
Binary file test-data/output2.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output2_binary.tif
Binary file test-data/output2_binary.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output2_binary.tiff
Binary file test-data/output2_binary.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output3.tif
Binary file test-data/output3.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output3.tiff
Binary file test-data/output3.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output3_binary.tif
Binary file test-data/output3_binary.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output3_binary.tiff
Binary file test-data/output3_binary.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output4.tif
Binary file test-data/output4.tif has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output4.tiff
Binary file test-data/output4.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output5.tiff
Binary file test-data/output5.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/output6.tiff
Binary file test-data/output6.tiff has changed
diff -r 2ae122d5d85a -r 64c155acb864 test-data/rois-illegal1.geojson
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rois-illegal1.geojson	Mon May 12 14:01:09 2025 +0000
@@ -0,0 +1,40 @@
+{
+    "type": "FeatureCollection",
+    "features": [
+        {
+            "type": "Feature",
+            "id": "a5d9de43-1a4a-4194-b06d-a6c6d0f81f91",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            201,
+                            48
+                        ],
+                        [
+                            292,
+                            48
+                        ],
+                        [
+                            292,
+                            184
+                        ],
+                        [
+                            201,
+                            48
+                        ],
+                        [
+                            201,
+                            48
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation",
+                "name": "1"
+            }
+        }
+    ]
+}
diff -r 2ae122d5d85a -r 64c155acb864 test-data/rois-illegal2.geojson
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rois-illegal2.geojson	Mon May 12 14:01:09 2025 +0000
@@ -0,0 +1,40 @@
+{
+    "type": "FeatureCollection",
+    "features": [
+        {
+            "type": "Feature",
+            "id": "a5d9de43-1a4a-4194-b06d-a6c6d0f81f91",
+            "geometry": {
+                "type": "Point",
+                "coordinates": [
+                    [
+                        [
+                            201,
+                            48
+                        ],
+                        [
+                            292,
+                            48
+                        ],
+                        [
+                            292,
+                            184
+                        ],
+                        [
+                            201,
+                            184
+                        ],
+                        [
+                            201,
+                            48
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation",
+                "name": "1"
+            }
+        }
+    ]
+}
diff -r 2ae122d5d85a -r 64c155acb864 test-data/rois-noname.geojson
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rois-noname.geojson	Mon May 12 14:01:09 2025 +0000
@@ -0,0 +1,104 @@
+{
+    "type": "FeatureCollection",
+    "features": [
+        {
+            "type": "Feature",
+            "id": "9ef0487b-29a5-4e53-8eca-5c9dbf7bfc80",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            124,
+                            149
+                        ],
+                        [
+                            183,
+                            149
+                        ],
+                        [
+                            183,
+                            275
+                        ],
+                        [
+                            124,
+                            275
+                        ],
+                        [
+                            124,
+                            149
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation"
+            }
+        },
+        {
+            "type": "Feature",
+            "id": "a5d9de43-1a4a-4194-b06d-a6c6d0f81f91",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            201,
+                            48
+                        ],
+                        [
+                            292,
+                            48
+                        ],
+                        [
+                            292,
+                            184
+                        ],
+                        [
+                            201,
+                            184
+                        ],
+                        [
+                            201,
+                            48
+                        ]
+                    ]
+                ]
+            }
+        },
+        {
+            "type": "Feature",
+            "id": "b7b348f7-5438-47e3-a4ae-dbc923aa4e3b",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            151,
+                            95
+                        ],
+                        [
+                            260,
+                            95
+                        ],
+                        [
+                            260,
+                            162
+                        ],
+                        [
+                            151,
+                            162
+                        ],
+                        [
+                            151,
+                            95
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation"
+            }
+        }
+    ]
+}
diff -r 2ae122d5d85a -r 64c155acb864 test-data/rois.geojson
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/rois.geojson	Mon May 12 14:01:09 2025 +0000
@@ -0,0 +1,110 @@
+{
+    "type": "FeatureCollection",
+    "features": [
+        {
+            "type": "Feature",
+            "id": "9ef0487b-29a5-4e53-8eca-5c9dbf7bfc80",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            124,
+                            149
+                        ],
+                        [
+                            183,
+                            149
+                        ],
+                        [
+                            183,
+                            275
+                        ],
+                        [
+                            124,
+                            275
+                        ],
+                        [
+                            124,
+                            149
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation",
+                "name": "2"
+            }
+        },
+        {
+            "type": "Feature",
+            "id": "a5d9de43-1a4a-4194-b06d-a6c6d0f81f91",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            201,
+                            48
+                        ],
+                        [
+                            292,
+                            48
+                        ],
+                        [
+                            292,
+                            184
+                        ],
+                        [
+                            201,
+                            184
+                        ],
+                        [
+                            201,
+                            48
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation",
+                "name": "1"
+            }
+        },
+        {
+            "type": "Feature",
+            "id": "b7b348f7-5438-47e3-a4ae-dbc923aa4e3b",
+            "geometry": {
+                "type": "Polygon",
+                "coordinates": [
+                    [
+                        [
+                            151,
+                            95
+                        ],
+                        [
+                            260,
+                            95
+                        ],
+                        [
+                            260,
+                            162
+                        ],
+                        [
+                            151,
+                            162
+                        ],
+                        [
+                            151,
+                            95
+                        ]
+                    ]
+                ]
+            },
+            "properties": {
+                "objectType": "annotation",
+                "name": "3"
+            }
+        }
+    ]
+}
\ No newline at end of file