def find_facade_cuts(facade_sig, dilation_amount=60):
edges = facade_sig > facade_sig.mean() + facade_sig.std()
edges = binary_closing(edges, structure=np.ones(dilation_amount))
run, start, val = run_length_encode(edges)
result = []
for s, e in zip(start[val], start[val] + run[val]):
result.append(s + facade_sig[s:e].argmax())
result = [0] + result + [len(facade_sig) - 1]
result = np.unique(result)
return np.array(result)
评论列表
文章目录