细胞图像处理大作业
更新时间:2024-01-17 20:05:01 阅读量: 教育文库 文档下载
实验课程:数字图像处理实验名称:实验地点:计算机中心学生姓名:学 号:指导教师:实验时间:
实验报告
骨髓细胞计数 405
年 月
一、实验要求:
将被染成紫色的白细胞从较暗的红细胞和背景中分割出来,去除细胞粘连,用数字图像处理技术实现图像中紫色细胞的自动计数。 二、实验过程
1、对原图像进行预处理
(1)由于源图像是彩色图像,为了便于进行仿真实验,在进行处理之前需要将彩色图像转换为灰度图像。用到的函数为rgb2gray
(2)获取图像的全局阈值并将图像二值化,以方便后面实验中对图像进行形态学修饰的进行
(3)在图像获取过程中,由于热传输等影响,会产生高斯噪声,较适合采用中值滤波的方法进
行去噪处理。因为,中值滤波的方法能够较好的保护细胞的边缘信息,满足本实验的需求;
(4)删除小面积图像
(5)对图像进行开闭运算。开和闭这两种运算可以除去比结构元素小的特定图像细节,同时保证不产生全局几何失真。开运算可以把比结构元素小的突刺滤掉,切断细长搭接而起到分离作用;闭运算可以把比结构元素小的缺口或孔填充上,搭接短的间隔而起到连接作用。 (6)对图像进行孔洞填充, (7)开闭运算会留下一些偏离的孤立像素,使用bwareaopen来移除少于特定像素个数的斑点。BW2= bwareaopen(BW,P)从二值图像中移除所以少于P像素值的连通块,得到另外的二值图像BW3。 实验图片显示如下:
1
2、分割计数
本次实验采用分水岭法对进行过预处理的图像进行分割和计数,最后计数得紫色细胞个数为35
三、实验心得体会
通过这学期的学习,我获益匪浅,学到了很多有趣的知识,知道了什么是灰度图,二值化,阈值,分水岭等很多有趣的图像处理方法,知道了怎么使用matlab帮助查找自己不会的函数,也成功的自己独立完成了细胞的计数。这门课也激起了我对图像处理的兴趣,课程结束后也会继续学习更多的函数知识,不会停止这方面的学习。
2
四、源代码
Image = imread('b1.bmp'); subplot(321); imshow(Image); title('原图'); Image=rgb2gray(Image); subplot(322); imshow(Image); title('灰度图');
Theshold = graythresh(Image); BW = im2bw(Image,1.3*Theshold);
subplot(323),imshow(BW); title('二值化图像'); BW1= medfilt2(BW,[13 13]); BW1=bwareaopen(BW,50,4) subplot(3,2,4); imshow(BW1);
title('删除小面积后图像'); se=strel('disk',5); fo=imopen(BW1,se); foc=imclose(fo,se) subplot(3,2,5); imshow(foc);
title('开闭运算后图像'); BW2=imfill(foc,'holes'); subplot(3,2,6); imshow(BW2);
title('孔洞填充后图像'); BW3=bwareaopen(BW2,20); D=bwdist(~BW3); D=-D;
D(~BW3)=-Inf; L=watershed(D);
rgb=label2rgb(L,'jet');
3
figure; imshow(rgb,[]);
title('分水岭分割计算后图像'); number=max(max(L))-1; number
4
方法一源程序:
A=imread('y1.bmp'); B=rgb2gray(A);
h=[0 1 0,1 -4 1,0 1 0];
J1=conv2(B,h,'same'); K1=medfilt2(J1,[3 3]);¨ BW=edge(K1,'canny');
SE= strel('line',11,90); BW1=imdilate(BW,SE);
BW2=bwmorph(BW1,'thin',Inf);
[H, theta, rho] = hough(BW2, 'ThetaResolution', 0.2);
subplot(221),imshow(H, [], 'XData', theta, 'YData', rho, 'InitialMagnification', 'fit')
axis on, axis normal
xlabel('\\theta'), ylabel('\\rho') , peaks = houghpeaks(H, 5); hold on
plot(theta(peaks(:, 2)), rho(peaks(:, 1)), ... 'linestyle', 'none', 'marker', 's', 'color', 'w'); lines = houghlines(BW2, theta, rho, peaks); subplot(222),imshow(BW2), hold on for k = 1:length(lines)
5
xy = [lines(k).point1 ; lines(k).point2];
plot(xy(:,1), xy(:,2), 'LineWidth', 4, 'Color', [.8 .8 .8]); end
6
正在阅读:
细胞图像处理大作业01-17
2002年江西省南昌市中考数学试卷11-27
请求协助函范本_关于请求协助函范文03-26
2022央视中秋晚会观后感(6篇)03-26
Xpath工具使用教程01-04
我并不孤单作文800字07-14
部编版三年级上册语文教学计划及进度表03-26
学校开展学生课后服务工作计划(精选三篇)03-26
中秋加班通知范文03-26
速卖通创业策划书03-09
- exercise2
- 铅锌矿详查地质设计 - 图文
- 厨余垃圾、餐厨垃圾堆肥系统设计方案
- 陈明珠开题报告
- 化工原理精选例题
- 政府形象宣传册营销案例
- 小学一至三年级语文阅读专项练习题
- 2014.民诉 期末考试 复习题
- 巅峰智业 - 做好顶层设计对建设城市的重要意义
- (三起)冀教版三年级英语上册Unit4 Lesson24练习题及答案
- 2017年实心轮胎现状及发展趋势分析(目录)
- 基于GIS的农用地定级技术研究定稿
- 2017-2022年中国医疗保健市场调查与市场前景预测报告(目录) - 图文
- 作业
- OFDM技术仿真(MATLAB代码) - 图文
- Android工程师笔试题及答案
- 生命密码联合密码
- 空间地上权若干法律问题探究
- 江苏学业水平测试《机械基础》模拟试题
- 选课走班实施方案
- 图像处理
- 细胞
- 作业
- 2017-2018学年度第一学期北师大版五年级语文上册期末试卷及答案
- 上海临港新城
- 列方程求比一个数的几倍多几或少几的练习题
- 财务绩效评价文献综述
- 赔偿法习题
- 年产57000万片高效硅片(1200MW光伏电池)项目可行性研究报告 - 图文
- 景观工程技术标 2 - 图文
- 通风技术管理制度
- PEP小学英语三年级上册英语试卷分析
- 含7套高考模拟卷高考英语语法复习专题讲练-非谓语动词 - 图文
- 全国2013年10月自学考试00182《公共关系学》试卷及答案 - 图文
- 中层干部廉政谈话和诫勉谈话制度
- 2019届高考语文三年高考试题分类汇编解析版 名句默写
- 校园文化建设实施方案(深河中学 朱红伟)
- 2015届高三数学一轮复习练习:达标练习11算法、统计、概率 必修三
- 建筑工程专业二级建造师继续教育选修课考试题库及参考答案精华版
- 第四章社会主义建设道路初步探索的理论成果
- 小儿穴位帖敷处方
- 二手车鉴定评估师(三级)操作A卷参考答案
- 最新人教版四年级语文下册单元测试题及答案全套