1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
###############################################################################
# pytest_composite.py: composite min/max
# Author(s): Pieter.Kempeneers@ec.europa.eu
# Copyright (c) 2016-2019 European Union (Joint Research Centre)
# License EUPLv1.2
#
# This file is part of jiplib
###############################################################################
# History
# 2017/10/24 - Created by Pieter Kempeneers (pieter.kempeneers@ec.europa.eu)
# Change log
import argparse
import os
import jiplib as jl
parser=argparse.ArgumentParser()
parser.add_argument("-input","--input",help="Path of the raster dataset",dest="input",required=True,type=str)
args = parser.parse_args()
try:
jim0=jl.createJim(args.input)
stats0=jim0.getStats({'function':['min','max'], 'band':0})
stats5=jim0.getStats({'function':['min','max'], 'band':5})
stats10=jim0.getStats({'function':['min','max'], 'band':10})
jl0=jl.JimList([jim0.cropBand({'band':0}),jim0.cropBand({'band':5}),jim0.cropBand({'band':10})])
print("stack bands")
# jimstack=jl0.stackBand({'verbose':2})
jimstack=jl0.stackBand()
print("creating min composite")
jim1=jl0.composite({'crule':'minband'}).clearNoData()
print("get min statistics")
statmin=jim1.getStats({'function':['min']})
if statmin['min']>min(stats0['min'],stats5['min'],stats10['min']):
print("Failed: min",statmin['min'],stats0['min'],stats5['min'],stats10['min'])
throw()
print("creating statistics profile")
jimsp=jimstack.statProfile({'function':'min'})
if not jim1.isEqual(jimsp):
print("Failed: statProfile min")
throw()
print("creating max composite")
jim1=jl0.composite({'crule':'maxband'}).clearNoData()
jl0.close()
print("get max statistics")
statmax=jim1.getStats({'function':['max']})
if statmax['max']<max(stats0['min'],stats5['min'],stats10['min']):
print("Failed: max",statsc['max'],stats0['max'],stats5['max'],stats10['max'])
throw()
print("creating statistics profile")
jimsp=jimstack.statProfile({'function':'max'})
if not jim1.isEqual(jimsp):
print("Failed: statProfile max")
throw()
jl0.close()
jim1.close()
jimsp.close()
print("Success: composite min/max")
except:
print("Failed: composite min/max")
|
