8000 How to get extended_persistent_barcode? · Issue #45 · GUDHI/TDA-tutorial · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content
How to get extended_persistent_barcode? #45
@tanjia123456

Description

@tanjia123456

Hello, sorry to bother you again.
Q1: I have use brain cortical curv of each vertex as filteration value, but when i use the following code:
for i in range(len(data_vertices)):
st.assign_filtration([i], filtration = filteration_value[i])
#st.assign_filtration([0], filtration = filteration_value[0])#success!
#st.assign_filtration([0], filtration = filteration_value[i])#success!
error:Process finished with exit code -1073741819 (0xC0000005)
I have tried some solution, but it's doesn't work! I'm curious why filter values can only be fed one node at a time.

Q2: I add filter value manually
st.assign_filtration([0], filtration = filteration_value[0])
st.assign_filtration([1], filtration = filteration_value[1])
st.assign_filtration([2], filtration = filteration_value[2])
st.assign_filtration([3], filtration = filteration_value[3])

But you said use extended_persistence, I have use it, but :
st2 = gd.SimplexTree()
st2.extend_filtration()
dgms2 = st2.extended_persistence(min_persistence=1e-5)
gd.plot_persistence_barcode(dgms2)#IndexError: list index out of range
plt.show()
the error message:
IndexError: list index out of range
but when I use gd.plot_persistence_barcode(dgm),use ordinary persistent barcode, there's no problem. What code should I use to get the extended continuous bar code directly?

Q3: I've added filter values to each node. What pair of distances do I need to use to weight the edges?

Q4: Afetr I get the extended_persistent_barcode, How to get betti???

My datasets and codes are in:https://github.com/tanjia123456/Brain_area_PH/blob/main/compute_betti_gudhi3.py

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions

      0