Mercurial > repos > imgteam > points2labelimage
changeset 5:4a49f74a3c14 draft default tip
planemo upload for repository https://github.com/BMCV/galaxy-image-analysis/tree/master/tools/points2labelimage/ commit 9a40a5d1e1008c26cc327c6d163df2a1af22a1a0
| author | imgteam |
|---|---|
| date | Mon, 12 May 2025 14:01:26 +0000 |
| parents | aef9911c0d5c |
| children | |
| files | creators.xml points2label.py points2label.xml test-data/output1_binary.tif test-data/output1_binary.tiff test-data/output2.tif test-data/output2.tiff test-data/output2_binary.tif test-data/output2_binary.tiff test-data/output3.tif test-data/output3.tiff test-data/output3_binary.tif test-data/output3_binary.tiff test-data/output4.tif test-data/output4.tiff test-data/output5.tiff test-data/output6.tiff test-data/rois-illegal1.geojson test-data/rois-illegal2.geojson test-data/rois-noname.geojson test-data/rois.geojson |
| diffstat | 21 files changed, 456 insertions(+), 27 deletions(-) [+] |
line wrap: on
line diff
--- a/creators.xml Wed Apr 23 14:37:42 2025 +0000 +++ b/creators.xml Mon May 12 14:01:26 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>
--- a/points2label.py Wed Apr 23 14:37:42 2025 +0000 +++ b/points2label.py Mon May 12 14:01:26 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, )
--- a/points2label.xml Wed Apr 23 14:37:42 2025 +0000 +++ b/points2label.xml Mon May 12 14:01:26 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>
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rois-illegal1.geojson Mon May 12 14:01:26 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" + } + } + ] +}
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rois-illegal2.geojson Mon May 12 14:01:26 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" + } + } + ] +}
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rois-noname.geojson Mon May 12 14:01:26 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" + } + } + ] +}
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/rois.geojson Mon May 12 14:01:26 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