亚像素边缘检测的matlab代码，比一般的边缘检测代码精确

亚像素边缘检测的matlab代码，比一般的边缘检测代码精确

%   感谢亲亲使用此代码，此代码解决您的问题了吗~（@^_^@）~
%   没解决的话告诉亲亲一个好消息，我这里可以1分钱帮助代码改错，还提供1分钱成品代码（′▽‘〃）哦~
%   登录淘宝店铺“大成软件工作室”便可领取
%   是的，亲亲真的没有看错，挠破头皮的问题真的1分钱就可以解决了\（^o^）/YES!
%   小的这就把传送门给您，记得要收藏好哦（づ￣3￣）づ╭～
%   传送门：https://item.taobao.com/item.htm?spm=a1z10.1-c.w4004-15151018122.5.uwGoq5&id=538759553146
%   如果传送门失效，亲亲可以来店铺讨要，客服MM等亲亲来骚扰哦~（*/ω╲*）
function [oimgXY]=baseZernike1（img）
load M2;     %获取模板7*7的
%卷积运算
Z00=conv2（M00img）;
Z11I=conv2（M11Iimg）;
Z11R=conv2（M11Rimg）;
Z20=conv2（M20img）;
Z31I=conv2（M31Iimg）;
Z31R=conv2（M31Rimg）;
Z40=conv2（M40img）;
%截掉多余部分
Z00=Z00（4:end-34:end-3）;
Z11I=Z11I（4:end-34:end-3）;
Z11R=Z11R（4:end-34:end-3）;
Z20=Z20（4:end-34:end-3）;
Z31I=Z31I（4:end-34:end-3）;
Z31R=Z31R（4:end-34:end-3）;
Z40=Z40（4:end-34:end-3）;
%设置k，oimg，theta，Zz11，l等的初始矩阵，以防后来矩阵运算大小不匹配。
k=zeros（size（img1）size（img2））;
oimg=k;
theta=k;
theta31=k;
Zz11=k;
Zz31=k;
z11=k;
z31=k;
Zz40=k;
Zz20=k;
h=k;
bounder=k;
l=100*ones（size（img1）size（img2））;
l1=100*ones（size（img1）size（img2））;
l2=100*ones（size（img1）size（img2））;
%         防止分子、分母均为0而产生的NAN
a=abs（Z11R）<0.0001;
b=abs（Z11I）<0.0001;
c=a&b;Z11R（c）=1;Z11I（c）=1;
theta=atan（Z11I./Z11R）;
g=abs（Z31R）<0.0001;
j=abs（Z31I）<0.0001;
f=g&j;Z31R（f）=1;Z31I（f）=1;
theta31=atan（Z31I./Z31R）;

Zz11=Z11R.*cos（theta）+Z11I.*sin（theta）;
Zz31=Z31R.*cos（theta31）+Z31I.*sin（theta31）; 
Zz20=Z20;
Zz40=Z40;
l1=sqrt（（5*Zz40+3*Zz20）./（8*Zz20））;
l2=sqrt（（5*Zz31+Zz11）./（6*Zz11））;
l=（l2+l1）/2;
 e=abs（l）>1/（3.5*sqrt（2））;
 k=1.5.*Zz11./（（1-l2.^2）.^1.5）;
 k（e）=0;
 h=（Z00-（k*pi）/2+k.*asin（l2）+k.*l2.*sqrt（1-l2.^2））./pi;
a=abs（l2-l1）<1;
 b=abs（k）>=1.7;
%b=abs（k）>max（img（:））/10;
% amp=sqrt（Z11R.^2+Z11I.^2）;
%  t=amp>15;
% ab分别为距离和边缘强度判断结果c为之前需要去掉的NAN部分
 %d=a&b&t&~c&~f;
d=a&b&~c&~f;
oimg（d）=1;
bounder（3:end-23:end-2）=1;
oimg（~bounder）=0;

figureimshow（oimg）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2016-09-23 05:36  baseZernike1\
     文件        2116  2016-09-23 10:11  baseZernike1\baseZernike1.m
     文件        1255  2009-03-13 09:40  baseZernike1\M2.mat
     文件      786872  2009-03-11 16:01  baseZernike1\基于Zernike正交矩的图像亚像素边缘检测算法改进.pdf