python 计算骨架中轴的距离

计算骨架中轴的距离可以使用scikit-image库中的函数`skimage.morphology.medial_axis`和`scipy.spatial.distance.cdist`。

首先，使用`skimage.morphology.medial_axis`函数计算骨架的中轴：

from skimage.morphology import medial_axis
import numpy as np

# 生成一个随机的二值图像骨架
skeleton = np.zeros((100,100))
skeleton[20:80, 20:80] = 1
skeleton[30:70, 30:70] = 0
skeleton = medial_axis(skeleton)

接下来，使用`scipy.spatial.distance.cdist`函数计算中轴上每个点之间的距离：

from scipy.spatial.distance import cdist

# 提取中轴上的点坐标
coords = np.column_stack(np.where(skeleton))

# 计算中轴上每个点之间的距离
distances = cdist(coords, coords)

这将生成一个`distances`数组，其中第`i`行和第`j`列的元素表示中轴上第`i`个点和第`j`个点之间的距离。可以使用这个数组来计算中轴的总长度：

# 计算中轴的总长度
total_distance = np.sum(distances) / 2

注意，由于`distances`数组是对称的，所以需要将总距离除以2。