想要利用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的那條線
應該可以得到軟骨表面以及軟骨與骨頭接觸面兩條曲線,
但之後就不知道要怎麼繼續寫才能得到兩條曲線間的歐幾里德距離。
因為還是新手,想要請教各位高手關於這方面的意見,非常感謝!!