def test_ATGP(data, wvl, mask, path):
print('Testing ATGP')
atgp = eea.ATGP()
pr = profile()
#U = atgp.extract(data, 8, normalize=True)
U = atgp.extract(data, 8, normalize=True, mask=mask)
stat(pr)
print(' End members indexes:', atgp.get_idx())
atgp.plot(path, axes=wvl, suffix='test1')
atgp.plot(path, suffix='test2')
U = U[[0, 1], :]
test_amap(data, U, 'ATGP', path, mask, amaps='UCLS')
def ATGP(data, wvl, mask, path):
print('Testing ATGP')
atgp = eea.ATGP()
U = atgp.extract(data, 4, normalize=True)
atgp.plot(path, suffix='test2')
return U