<tool id="ip_landmark_registration" name="Landmark Registration" version="0.0.3" profile="20.05">
    <description>estimates the affine transformation matrix</description>
    <requirements>
        <requirement type="package" version="0.18.1">scikit-image</requirement>
        <requirement type="package" version="1.6.2">scipy</requirement>
        <requirement type="package" version="1.2.4">pandas</requirement>
        <requirement type="package" version="1.20.2">numpy</requirement>
    </requirements>
    <command detect_errors="aggressive">
    <![CDATA[ 
        python '$__tool_directory__/landmark_registration.py'
        '$fn_pts1'
        '$fn_pts2'
        '$fn_tmat'
        #if $algo_option.algo == 'ransac'
        --res_th $algo_option.res_thr
        --max_ite $algo_option.max_iter
        #end if
]]></command>
    <inputs>
        <param name="fn_pts1" type="data" format="tabular" label="Coordinates of SRC landmarks (tsv file)" />
        <param name="fn_pts2" type="data" format="tabular" label="Coordinates of DST landmarks (tsv file)" />
        <conditional name="algo_option">
            <param name="algo" type="select" label="Select the algorithm">
                <option value="ransac" selected="True">RANSAC</option>
                <option value="ls">Least Squares</option>
            </param>
            <when value="ransac">
                <param name="res_thr" type="float" value="2.0" label="Maximum distance for a data point to be classified as an inlier." />
                <param name="max_iter" type="integer" value="100" label="Maximum number of iterations for random sample selection." />
            </when>
            <when value="ls"></when>
        </conditional>
    </inputs>
    <outputs>
        <data format="tabular" name="fn_tmat" />
    </outputs>
    <tests>
        <test>
            <param name="fn_pts1" value="points_moving.tsv"/>
            <param name="fn_pts2" value="points_fixed.tsv"/>
            <param name="algo_option.algo" value="ls"/>
            <output name="fn_tmat" value="tmat.tsv" ftype="tabular" compare="diff" lines_diff="6"/>
        </test>
        <test>
            <param name="fn_pts1" value="points_moving.tsv"/>
            <param name="fn_pts2" value="points_fixed.tsv"/>
            <param name="algo_option.algo" value="ransac"/>
            <param name="algo_option.res_thr" value="2"/>
            <param name="algo_option.max_iter" value="100"/>
            <output name="fn_tmat" value="tmat.tsv" ftype="tabular" compare="diff" lines_diff="6"/>
        </test>
    </tests>
    <help>
    **What it does**

    This tool estimates the transformation matrix between two sets of 2d points.

    About the format of landmark coordinates in the input TSV table: Columns with header "x" and "y" are for x- and y-coordinate, respectively.
    </help>
    <citations>
        <citation type="doi">10.1016/j.jbiotec.2017.07.019</citation>
    </citations>
</tool>