How to create a function to define a square

1 次查看（过去 30 天）

显示更早的评论

TOKEN 2019-4-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/457612-how-to-create-a-function-to-define-a-square

评论： Walter Roberson 2019-4-21

How to create a function to define a square

3 个评论
显示 1更早的评论隐藏 1更早的评论

TOKEN 2019-4-21

在 MATLAB Online 中打开

I am actually doing the shapes recoginition using optimization tools, I have to define a square first.

here is the code

clc
clear
M0=zeros(1024,1024); % Background of zeros
x=1:1024; 
y=1:1024;
[X,Y]=meshgrid(x,y); % grid of pixel positions
for in=11
    A=double(imread([num2str(in,'%4.3d'),'.jpg']))/255; % Load image and scale to 1
    M=@(p) M0+double((X-p(1)).^2+(Y-p(2)).^2<(p(3))^2); % define your teplate 
                                                        % You may need a separate funtion file
    
    fun=@(p) sum(sum(abs(A-M(p)).^2)); % The objective funtion to minimize
     
    % find the avergae center position and ditribution of the nonzero
    % pixels in image
    pos=find(A>0);
    x0=mean(X(pos));
    y0=mean(Y(pos));
    d=(range(X(pos))+range(Y(pos)))/2;
    
    trl=[x0 y0 d/2]; % The trial solution % number of parameters depend on the shape
    
    lb=[0 0 0];% lower bounds of the prameters
    ub=[1024 1024 512];% upper bounds of the prameters
    
    
    % Optimization process
    %You dont have to change the below lines
    opts=optimoptions('ga'); 
    opts.InitialPopulationMatrix=trl;
    opts.Display='iter';     
    [sol,fval]=ga(fun,length(trl),[],[],[],[],lb,ub,[],opts);
    
    if fval<100 
        'circle'
    else
        'unknown'
    end
    
    figure(1)
    clf
    surf(M(sol),'linestyle','none')
    view(2)
    daspect([1 1 1])
    
    figure(2)
    clf
    surf(A,'linestyle','none')
    view(2)
    daspect([1 1 1])     
end