Allowing TPS transformation

preprocessor/preprocessor/steps/georeference.py

@@ -30,7 +30,8 @@ def georeference_step(source_dir: os.PathLike, target_dir: os.PathLike, type: st
         georef_func(filename, target_filename, **options)
 
 
-def gcp_georef(input_filename: os.PathLike, target_filename: os.PathLike, order: int=1, projection: str='EPSG:4326'):
+def gcp_georef(input_filename: os.PathLike, target_filename: os.PathLike, order: int=1, projection: str='EPSG:4326',
+               tps: bool=False):
     succeded = False
 
     # simple case: get the geotransform from some GCPs
@@ -57,6 +58,7 @@ def gcp_georef(input_filename: os.PathLike, target_filename: os.PathLike, order:
             options=gdal.WarpOptions(
                 polynomialOrder=order,
                 dstSRS=projection,
+                tps=tps,
             )
         )