我有两张图像说 gt 地面实况图像和 a 分割图像。我想计算这两者之间的MSD。我写了一个代码。只是想知道下面的代码是否正常。
grd=np.array(measure.find_contours(gt, 0.99))
grd_x=np.array(np.transpose(grd[:,:,0]))
grd_y=np.array(np.transpose(grd[:,:,1]))
seg=np.array(measure.find_contours(a,0.99))
seg_x=np.array(np.transpose(seg[:,:,0]))
seg_y=np.array(np.transpose(seg[:,:,1]))
distance=[]
for i in range(0,len(seg_x)):
temp_vector_x=seg_x[i]*np.ones(len(grd_x))
temp_vector_y=seg_y[i]*np.ones(len(grd_y))
tt=np.array(np.sqrt((np.square(temp_vector_x-grd_x))+(np.square(temp_vector_y-grd_y))))
ttmin=tt.min()
distance.append(ttmin)
results=(np.sum(distance))/len(seg_x)