import numpy as np
import yaml

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

def run(rasters, treshold=1e4, areas=None, sd=None, moi=None):
    # Parse attribute 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)

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

    # Compute vectors
    out_vectors.run()
    
    return out_vectors.data

def version():
    return 'v0.0'