# HG changeset patch
# User thomaswollmann
# Date 1546857567 18000
# Node ID 470fd1c1c10cc0e4ec920492b77a649b41cc222f

planemo upload for repository https://github.com/BMCV/galaxy-image-analysis/tools/projective_transformation_points/ commit 787ebcc8daa1834214bc92c201c921c704ef2d1f

diff -r 000000000000 -r 470fd1c1c10c projective_transformation_points.py
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/projective_transformation_points.py	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,35 @@
+from skimage.transform import ProjectiveTransform
+import numpy as np
+import pandas as pd
+import argparse
+
+
+def warp_coords_batch(coord_map, coords, dtype=np.float64, batch_size=1000000):
+    tf_coords = coords.astype(np.float32)
+
+    for i in range(0, (tf_coords.shape[0]//batch_size+1)):
+        tf_coords[batch_size*i:batch_size*(i+1)] = coord_map(tf_coords[batch_size*i:batch_size*(i+1)])
+
+    return np.unique(np.round(tf_coords).astype(coords.dtype),axis=0)
+
+
+def transform(coords, warp_matrix, out):
+    indices = np.array(pd.read_csv(coords, delimiter="\t"))
+    a_matrix = np.array(pd.read_csv(warp_matrix, delimiter=",", header=None))
+    
+    trans = ProjectiveTransform(matrix=a_matrix)
+    warped_coords = warp_coords_batch(trans, indices)
+
+    df = pd.DataFrame()
+    df['x'] = warped_coords[:,0]
+    df['y'] = warped_coords[:,1]
+    df.to_csv(out, index = False, sep="\t")
+
+
+if __name__ == "__main__":
+    parser = argparse.ArgumentParser(description="Transform coordinates")
+    parser.add_argument("coords", help="Paste path to .csv with coordinates to transform (tab separated)")
+    parser.add_argument("warp_matrix", help="Paste path to .csv that should be used for transformation (, separated)")
+    parser.add_argument("out", help="Paste path to file in which transformed coords should be saved (tab separated)")
+    args = parser.parse_args()
+    transform(args.coords, args.warp_matrix, args.out)
diff -r 000000000000 -r 470fd1c1c10c projective_transformation_points.xml
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/projective_transformation_points.xml	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,37 @@
+<?xml version="1.0"?>
+<tool id="ip_projective_transformation_points" name="Projective Transformation Points" version="0.0.2">
+    <description>Projective Transformation for Points</description>
+    <requirements>
+        <requirement type="package" version="0.14.0">scikit-image</requirement>
+        <requirement type="package" version="0.23.4">pandas</requirement>
+        <requirement type="package" version="1.15.2">numpy</requirement><!--conflict of 1.15.4 with mkl_fft!?-->
+    </requirements>
+    <command>
+    <![CDATA[
+        python '$__tool_directory__/projective_transformation_points.py'
+            '$moving_points'
+            '$warp_matrix'
+            '$out'
+    ]]>
+    </command>
+    <inputs>
+        <param name="moving_points" type="data" format="csv" label="Moving Points" />
+        <param name="warp_matrix" type= "data" format="csv" label="Warp Matrix" />
+    </inputs>
+    <outputs>
+       <data format="csv" name="out" />
+    </outputs>
+    <tests>
+      <test>
+        <param name="moving_points" value="points_tsv.csv" />
+        <param name="warp_matrix" value="warp_matrix.csv" />
+        <output name="out" value="out.csv" ftype="csv" />
+      </test>
+    </tests>
+    <help>
+This tool performs a projective transformation of the input (moving) points.
+    </help>
+    <citations>
+        <citation type="doi">10.1016/j.jbiotec.2017.07.019</citation>
+    </citations>
+</tool>
diff -r 000000000000 -r 470fd1c1c10c test-data/out.csv
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/out.csv	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,7 @@
+x	y
+18	24
+18	25
+18	26
+18	27
+20	22
+25	35
diff -r 000000000000 -r 470fd1c1c10c test-data/out2.csv
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/out2.csv	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,7 @@
+x	y
+18	24
+18	25
+18	26
+18	27
+20	22
+25	35
diff -r 000000000000 -r 470fd1c1c10c test-data/points_tsv.csv
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/points_tsv.csv	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,7 @@
+5	3
+10	2
+8	4
+15	15
+8	5
+8	6
+8	7
diff -r 000000000000 -r 470fd1c1c10c test-data/warp_matrix.csv
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/warp_matrix.csv	Mon Jan 07 05:39:27 2019 -0500
@@ -0,0 +1,3 @@
+1,0,10
+0,1,20
+0,0,1