i=imread('H:\enabie\IrisEye.jpg');
i=rgb2gray(i);
%imshow(i);
[row,col]=size(i);
j=floor(row/2);
x=0;
for k=1:col
if(i(j,k)>200)
x=x+1;
else
break;
end
end
y=x;
z=0;
for k=y+1:col
if(i(j,k)<200)
z=z+1;
else
break;
end
end
y2=y+z;
w=0;
for k=y2+1:col
if(i(j,k)>200)
w=w+1;
else
break;
end
end
mohitinternal=floor(((y2+w)-y2)*3.14);
y3=y2+w;
q=0;
for k=y3+1:col
if(i(j,k)<200)
q=q+1;
else
break;
end
end
y4=y3+q;
mohitexternal=floor((y4-y)*3.14);
redius=(y4-y)/2;
xx=0;
for deg=0:365
x1=j+redius*sin((deg*3.14)/180);
y1=(y4-redius)+redius*cos((deg*3.14)/180);
disp(i(floor(x1),floor(y1)))
xx=xx+1;
end
rec=[y2-y,2000];
for kk=1:(y2-y)
col2=1;
for deg=0:365
x1=j+(redius)*sin((deg*3.14)/180);
y1=(y4-redius)+(redius)*cos((deg*3.14)/180);
rec(kk,col2)=i(mod(floor(x1),row)+1,mod(floor(y1),col)+1);
col2=col2+1;
end
end
rec=uint8(rec);
[rr,cc]=size(rec);
for i=1:(y2-y)
for j=1:cc-1
rec(i,j)=floor((rec(i,j)+rec(i,j+1))/2);
end
end
imshow(rec);
matlab.rar_in_matlap code
版权申诉
177 浏览量
2022-09-20
23:23:16
上传
评论
收藏 1KB RAR 举报
![avatar](https://profile-avatar.csdnimg.cn/a64d28507aff44a1a968cd823e7c3cbe_weixin_42665255.jpg!1)
寒泊
- 粉丝: 78
- 资源: 1万+
最新资源
- 一个简单的MATLAB仿真示例,展示如何使用MATLAB进行基本的信号处理仿真:生成一个正弦波信号,并对其进行低通滤波处理
- 串口空闲中断 cubemax 任意长度数据
- GDAL-3.7.3-cp311-cp311-win-amd64.whl
- IE8或IE9浏览器下提升:请升级浏览器版本
- mat2pngmat2pngmat2png
- 一个简单的Python爬虫示例,用于爬取时光网电影排行榜上的电影名称和评分信息 这个示例仅用于教育目的,展示如何使用Python
- React全家桶框架.rar
- (PC+WAP)货物运输快递物流网站pbootcms模板 汽车贸易网站源码下载
- svg动画-源文件压缩包
- Python 中文数据结构和算法教程.rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback](https://img-home.csdnimg.cn/images/20220527035711.png)
![feedback-tip](https://img-home.csdnimg.cn/images/20220527035111.png)