想要利用matlab code来计算影像中软骨厚度,
文献中提到是取得软骨表面与骨头接触面之最短欧几里德距离(Euclidean distance),
目前的code如下:
I=dicomread('图档路径');imagesc(I)
colormap(gray)
[x,y]=ginput(5) 可手动选取5个点
p=polyfit(x,y,2);newy=polyval(p,x) 用二次多项式fitting得到新的y值
p=polyfit(x,y,2) 可得到fitting后方程式的系数
hold on
plot(x,y,'r',x,newy) 可画出fitting的那条线
应该可以得到软骨表面以及软骨与骨头接触面两条曲线,
但之后就不知道要怎么继续写才能得到两条曲线间的欧几里德距离。
因为还是新手,想要请教各位高手关于这方面的意见,非常感谢!!