import numpy as np

from astrodendro import Dendrogram
from astrodendro.analysis import PPStatistic

npixel = 3

data = np.loadtxt("3x3.txt")

data = np.reshape(data,(npixel,npixel))

print(data)

d = Dendrogram.compute(data, min_value=0.5, min_delta=0., min_npix=1)

# there is only one structure in total
print(d.trunk[0])
print(d.trunk[0].indices())

stat = PPStatistic(d.trunk[0])

xcen = stat.x_cen.value
ycen = stat.y_cen.value

print(xcen,ycen)

# Gives "None"
print(d.structure_at((round(xcen),round(ycen))))

#Gives the right result, but x- and y-coordinate are switched
print(d.structure_at((round(ycen),round(xcen))))