Zernike 亚像素边缘检测的matlab代码

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

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;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     786872  2009-03-11 16:01  基于Zernike正交矩的图像亚像素边缘检测算法改进.pdf

     文件       1255  2009-03-13 09:40  M2.mat

     文件       1611  2009-03-18 17:13  baseZernike1.m

----------- ---------  ---------- -----  ----

               789738                    3