#!/usr/bin/python
# -*- coding: utf-8 -*-
# \file dfc_sdaps.py
# \brief TODO
# \author Florent Guiotte <florent.guiotte@gmail.com>
# \version 0.1
# \date 27 août 2018
#
# TODO details

import numpy as np

import sys
sys.path.append('..')
import ld2dap

def run(rasters, treshold=1e4, areas=None, sd=None, moi=None):
    # Parse parameters type
    treshold = float(treshold)
    areas = None if areas is None else np.array(areas).astype(np.float).astype(np.int)
    sd = None if sd is None else np.array(sd).astype(np.float)
    moi = None if moi is None else np.array(moi).astype(np.float)

    # Pipelines
    loader = ld2dap.LoadTIFF(rasters)
    dfc_filter = ld2dap.Treshold(treshold)
    dfc_filter.input = loader
    sdaps = ld2dap.SelfDualAttributeProfiles(area=areas, sd=sd, moi=moi)
    sdaps.input = dfc_filter
    out_vectors = ld2dap.RawOutput()
    out_vectors.input = sdaps

    # Compute vectors
    out_vectors.run()

    return out_vectors.data

def version():
    return 'v0.0'