55 lines
1.4 KiB
Python
Executable File
55 lines
1.4 KiB
Python
Executable File
#!/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()
|