请用Python语言写一个代码，可以用来计算裂缝图像中裂缝的长度和宽度。步骤如下： 1.首先对输入的图片进行二值化 2.然后找出并标记所有连通域 3.利用medial_axis()函数对每个连通域提取骨架，然后计算骨架长度，绘制并保存提取后的骨架图像 4.分别获得每个骨架上的所有点，然后计算每一点到两边轮廓边缘的法向量的长度相加得到裂缝宽度。 5.分别输出每个裂缝的长度和最大宽度

# 以下是Python代码
import cv2
import numpy as np# 读取图像
img = cv2.imread('crack.jpg')# 二值化处理
ret, thresh = cv2.threshold(img,127,255,cv2.THRESH_BINARY)# 找出并标记所有连通域
contours, hierarchy = cv2.findContours(thresh, cv2.RETR_TREE, cv2.CHAIN_APPROX_SIMPLE)# 计算每个连通域的骨架长度
lengths = []
for cnt in contours:
# 提取骨架
skeleton = cv2.medial_axis(cnt, return_distance=False)
# 计算骨架长度
length = 0
for i in range(len(skeleton)-1):
length += np.sqrt(np.sum((skeleton[i]-skeleton[i+1])**2))
lengths.append(length)# 绘制并保存提取后的骨架图像
skel_img = np.zeros(img.shape, np.uint8)
for i in range(len(contours)):
skel = cv2.medial_axis(contours[i], return_distance=False)
cv2.drawContours(skel_img, [skel], 0, 255, 2)
cv2.imwrite('skel.jpg', skel_img)# 计算裂缝宽度
widths = []
for cnt in contours:
# 获得每个骨架上的所有点
points = cv2.findNonZero(cv2.medial_axis(cnt, return_distance=False))
# 计算每一点到两边轮廓边缘的法向量的长度相加得到裂缝宽度
width = 0
for pt in points:
dist1 = cv2.pointPolygonTest(cnt, pt, True)
dist2 = cv2.pointPolygonTest(cnt, pt, False)
width += dist1 + dist2
widths.append(width)# 输出每个裂缝的长度和最大宽度
for i in range(len(contours)):
print("Length of crack %d is %.2f, width of crack %d is %.2f" % (i+1, lengths[i], i+1, widths[i]))