#!/usr/bin/ipython
# -*- coding: utf-8 -*-
# \file test.py
# \brief TODO
# \author Florent Guiotte <florent.guiotte@gmail.com>
# \version 0.1
# \date 03 avril 2018
#
# TODO details

import numpy as np
import ld2dap
import logger

def diff_test():
    layers_files = [
        '../Data/phase1_rasters/DEM+B_C123/UH17_GEM051_TR.tif',
        '../Data/phase1_rasters/DEM_C123_3msr/UH17_GEG051_TR.tif',
        '../Data/phase1_rasters/DEM_C123_TLI/UH17_GEG05_TR.tif',
        '../Data/phase1_rasters/DSM_C12/UH17c_GEF051_TR.tif',
        '../Data/phase1_rasters/Intensity_C1/UH17_GI1F051_TR.tif',
        '../Data/phase1_rasters/Intensity_C2/UH17_GI2F051_TR.tif',
        '../Data/phase1_rasters/Intensity_C3/UH17_GI3F051_TR.tif'
    ]

    l = ld2dap.LoadTIFF(layers_files[0:3])
    t = ld2dap.Treshold(1e4)
    p = ld2dap.SelfDualAttributeProfiles(area=[100,1e4])
    f = ld2dap.Differential()
    #p = ld2dap.AttributeProfiles(area=[100,1e4,1e6])
    m = ld2dap.Merger()
    a = ld2dap.LocalFeatures([np.std], 7)#, sd=[.4,.6,.8], moi=[.5,.9])
    d = ld2dap.ShowFig(stack_id='all', symb=True)
    o = ld2dap.RawOutput()

    d.input = m
    m.input = f
    m.second.input = a
    a.input = f
    f.input = p
    p.input = t
    t.input = l

    d.run()


#from core import Input, Output, Filter

def main():
    logger.setup_logging()
    diff_test()

if __name__ == '__main__':
    main()