亚像素边缘检测matlab代码

执行快速准确的亚像素边缘检测，基于改进的zernike矩方法，有需要的可以下载~~

function img = CreateCircleEdgeImage（nrowsncolsucvcradius）
%
% Support function for Examples 3 and 4. This draws a bright circle on a 
% dark bakground in a 128x128 image.
%
% 
% Version History:
%   v1.0 - Initial release 06 Apr 2018
%
% Developed by:
% Dr. John A. Christian
% Rensselaer Polytechnic Institute （RPI）
%
% Copyright 2018 John Christian
%


%Upsampling of image
upsample = 51;

%Init image with all zeros
nrows_up = upsample * nrows;
ncols_up = upsample * ncols;
img_up = zeros（nrows_upncols_up）;

%Draw circle in higher resolution image
A = （upsample-1）/2 + 1 - upsample;
vmat_up = [1:1:nrows_up]‘*ones（1ncols_up）;
umat_up = ones（nrows_up1）*[1:1:ncols_up];
d_up = （（umat_up-A）/upsample - uc）.^2 + （（vmat_up-A）/upsample-vc）.^2;
img_up（d_up<=radius^2） = 1;

%Init loop parameters
delta = （upsample-1）/2;
BoxSize = upsample^2;

%Init 128x128 image
img = zeros（nrowsncols）;

%Sum-down large inmage into desired image size
for row_ii = 1:nrows
    
    row_ii_up = A + upsample*row_ii;
    row_range = [（row_ii_up-delta）:（row_ii_up+delta）];
    
    for col_ii = 1:ncols
        
        col_ii_up = A + upsample*col_ii;
        
        col_range = [（col_ii_up-delta）:（col_ii_up+delta）];
        
        box = img_up（row_rangecol_range）;
        
        img（row_iicol_ii） = sum（ box（:） ） / BoxSize ;
        
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1347  2018-04-06 22:24  CreateCircleEdgeImage.m
     文件         210  2018-04-06 22:24  __MACOSX\._CreateCircleEdgeImage.m
     文件        1315  2018-04-06 22:24  CreateStraightEdgeImage.m
     文件         210  2018-04-06 22:24  __MACOSX\._CreateStraightEdgeImage.m
     文件        1478  2018-04-06 22:24  Example1.m
     文件         210  2018-04-06 22:24  __MACOSX\._Example1.m
     文件        1823  2018-04-06 22:24  Example2.m
     文件         210  2018-04-06 22:24  __MACOSX\._Example2.m
     文件        1696  2018-04-06 22:24  Example3.m
     文件         210  2018-04-06 22:24  __MACOSX\._Example3.m
     文件        2054  2018-04-06 22:24  Example4.m
     文件         210  2018-04-06 22:24  __MACOSX\._Example4.m
     文件        7150  2018-04-06 22:24  SubpixelEdge.m
     文件         210  2018-04-06 22:24  __MACOSX\._SubpixelEdge.m
     文件        1314  2018-04-06 22:24  license.txt