矫直从x射线断层扫描图像

1视图(30天)
redmef
redmef 2018年7月4日
评论道: 乔纳斯 2018年7月6日

亲爱的所有,

我扫描一些没有刚性材料x射线断层扫描。我有一堆图片,不幸的是不可能连续样本像下面你可以看到。

我已经与ImageJ试图做一个图像,结果并不是那么糟糕。

问题是:我有一堆图片1500张图片和它与ImageJ是很长的。你有一个想法的直接材料堆栈和Matlab的图像吗?

最好的

接受的答案

乔纳斯
乔纳斯 2018年7月4日
编辑:乔纳斯 2018年7月4日
这里是我缝合,似乎工作得很好。
% %加载图片,旋转和转换为二进制
灰色= imread (“scan_00123.tif”);
灰色= rot90(灰色);
阈值= graythresh(灰色);
BW = im2bw(灰色、阈值);
% %找到图像的坐标和光滑的曲线
[y、x] =找到(BW = = 1);
y = y-mean (y);
图;
次要情节(1、2、1);
情节(x, y)
赢得= 1500;
y2 =光滑(x, y,赢得);
情节(x, y2,“r”)
% %修正向量的大小
[x,指数]=独特的(x);
xq = 1:大小(灰色,2);
yq = interp1 (x, y2(指数),xq,“线性”,“extrap”);
情节(yq xq,“——k”)
% %调整图像列
NewImage = uint16 (nan(大小(灰色,1)* 3,大小(灰色,2)));
ymid =大小(灰色,1)* 3/2;
i = 1:尺寸(灰色,2)
s = yq(我);
NewImage (ymid-s: ymid-1-s +大小(灰色,1),i) =灰色(:,我);
结束
图;
NewImage (isnan (NewImage)) = 0;
imshow (NewImage)
10评论
乔纳斯
乔纳斯 2018年7月6日
“具体的这种材料吗?”
一般来说,没有。大多数的脚本应该为任何工作薄材料。你可能略有调整分割取决于材料的颜色/亮度。
具体方向?
脚本,,适用于材料垂直对齐。简单地删除线
灰色= rot90(灰色);
材料水平拉伸。

登录置评。

更多的答案(0)

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!