کد تشخیص لبه و چین و چروک پنجره حداقل

مبنای این کد به این شکل است که با حرکت در طول عکس، و با انتخاب یک بازه در محور y ، کد عددی کل پنجره  در صورتی که دارای مقدار حداقل نباشد، برابر مقدار حداقل می کند.

 

کد متلب:

% Run Demo
[FileName,PathName] = uigetfile('*.jpg','Select the MATLAB code file');
I=imread([PathName,FileName]);
% I=imrotate(I,35,'bilinear');
I = im2double(imread([PathName,FileName]));
img=imrotate(I,35,'bilinear');
imshow(I)
%   I=Icrop ;
im2=I(:,:,1)*0+1;

figure;imshow( I);hold on;


facBoundryX2=[1:1:size(I,2)];


facBoundryY2=[1:1:size(I,1) ]*0+ (size(I,1) );

% plot( facBoundryX2,facBoundryY2,'.')
flagout=0;


%     imshow(I);hold on;
for ind2=1:numel(facBoundryX2)
    for ind1=1:numel(facBoundryY2)
        xI=facBoundryX2(ind2);
        yI=max(facBoundryY2(ind1),1);
        
        
        yVect=[yI-4:yI+4];
        %         yVect=[yI-10:yI+10 ];
        
        yVect=yVect(yVect>0 & yVect<=size(I,1));
        if isempty(yVect)
            flagout=1;
            break
        end
        
        minboc=rgb2gray(I(round( yVect) ,round( xI), 1:3));
        
        [v,in]=min( minboc  );
        facBoundryY2=[1:1:size(I,1) ]*0+ (size(I,1) )-ind1*1;
        
        %         hold on;plot(xI,yI,'ro');
        try
            if v== I( yI,xI)
                
            else
                
                
                im2(yVect(in),xI)=.5;double(v);   %abs(double(v)- double(I1( yI,xI)));
                ss=abs(double(v)- double(I( yI,xI)));
            end
            
        catch
            
        end
        
        
    end
    
    imshow(im2 )
    
    
    
    if flagout==1
        
        break
    end
    
    pause(0.000005)
end


%  mesh((double((im2))))

% end


%%%%%%%%%%%%%%%%%%%%%%%%%%-----------software details-----------%%%%%%%%%%%%%%%%%%%%%%

%devloped in Iran(ardebil)                                       %
%Copyright (C) 2017(1396) by moussa hasanzadeh as the head devloper %%
%All rights reserved.                                             % 
%programed in mathlab 2016                                 
%gmail:[e-mail]
%phone [phone]

 

پیام عذرخواهی از دوستان

من شرمندم بخدا، اینقدر به خاطر کرونا(البته کرونا نگرفتما! ) مشکلات روی سرم ریخت که اصلا نتونستم پست جدیدی بفرستم انشاالله که بتونیم در خدمتتون باشیم....

دوستانی که پیام دادند و نتونستم جواب بدم از اونا هم عذرخواهی می کنم
پست ها را  دوباره  پی می‌گیرم...

/

نحوه نوشتن یک تابع فیتنس(making a fitness function )

برداشت کلی این است که نوشتن یک تابع فیتنس عموما با الگوریتم‌های فراابتکاری معنی پیدا می‌کند، با این وجود یک تابع فیتنس چیزی نیست جز انتخاب جواب های بهینه از بین مجموعه‌هایی که به شکل تصادفی، حریصانه... انتخاب می‌شوند. تابع فیتنس می‌تواند از نوع مینیمم، ماکزیمم یا همگرا به سمت مقادیر یا بردارها و... باشد.

به عنوان یک مثال ساده برای تابع فیتنس برای ماکزیمم سازی  را  می‌توان به شکل زیر اورد.

 فقط حواستون باشه، برای نوشتن تابع فیتنس برای مینیمم سازی باید مقدار اولیه بردار پاسخ را برابر بی نهایت یا inf قرار بدید و برای ماکزیمم سازی این مقدار را برابر با 0 یا حداقل مقدار قرار بدید.

 

clc
clear;
close all;


[X,Y,Z] = peaks(100);
[val,ind]=min(Z(:));
v.val=val;
v.ind=ind;
cn=0;
hs=figure('name' ,  ['period: '  num2str(cn) ' of '  num2str((length(X(:))) )   ' curen  value:  '  num2str(v.val(end)),  ' real max value:  '  num2str( max(Z(:)))] );
 subplot(1,2,1);surf(X,Y,Z);hold on
 subplot(1,2,1);plot3(X(ind),Y(ind),Z(ind),'r*');hold on

while(true)
     ind= randi(length(X(:)),1);
     if  v.val(end)<Z(ind)
      v.val(end+1)=Z(ind);

          v.ind(end+1)= ind ;
 subplot(1,2,1);plot3(X(v.ind),Y(v.ind),Z(v.ind),'r*');hold on ; 
  subplot(1,2,2);plot3(X(v.ind),Y(v.ind),Z(v.ind),'r-*');hold on
  

     else
    disp('');
     end
      cn=cn+1;
 

     
 
     pause(0.002)

    if cn>=length(X(:)) || v.val(end)==max(Z(:))
           break;
    end
     hs.Name=['period: '  num2str(cn) ' of '  num2str((length(X(:))) )   ' curen  value:  '  num2str(v.val(end)),  ' real max value:  '  num2str( max(Z(:)))];
end

 

خروجی کد بالا به شکل زیر خواهد بود:

 

 

 

توزیع نودها در فضاهای نامتقارن(Distribution of nodes in asymmetric spaces)

اغلب نواحی که برای شبیه سازی شبکه‌‌های بی‌سیم مورد استفاده قرار می‌گیرند، نواحی متقارن و عموما مربعی شکل می‌باشند ، حال آنکه به شکل طبیعی امکان توزیع متقارن نودها وجود نداره، در این پست سعی کردیم که کد توزیع نودها در فضاهای غیرمتقارن را منتشر کنیم، امیدواریم که براتون مفید واقع شده باشه.

در توزیع نودها از روش انتشار نقاط منحنی استفاده شده است.

کد متلب برای توضیحات بالا

%%%%%%%%%%%%%%%%%%%%%%%%%%-----------software details-----------%%%%%%%%%%%%%%%%%%%%%%
%devloped in Iran(ardebil)                                       
%Copyright (C) 2020(1398) by moussa hasanzadeh as the head devloper %%
%All rights reserved.                                            
%programed in mathlab 2016                                 
%gmail:[e-mail]
%phone [phone]

%% clear commands
delete(allchild(0));% delete all figures and graphical objects
clear;% Remove items from workspace, freeing up system memory
clc;% clear the Console

%% path commands
file=matlab.desktop.editor.getActive;% get current script address
try % if dont occur any error this block will run
      fileDetail=dir(file.Filename);% get current script details
      folderName=fileDetail.folder;% get current script folder name
catch% else if an  error occurs then
      index=strfind( file.Filename,'\');% find index of back space(\) in path text
      folderName=file.Filename(1:index(end)-1);% select all path text from 1 to last Occurrence of \ as script folder name
end
paths=genpath(folderName);% make current path and all sub paths
addpath(paths);% add all paths in known paths of matlab
cd(folderName);% go to path that current script runned from that


%% net building  commands
net=[];% create an empty var as net object
net.para.width=1000;% width param of net
net.para.length=1000;% width param of net
net.para.height=5;% width param of net
net.para.side=6;% side of Polygon
net.para.nodNum=1000 ;% number of nodes of net work
net.para.simPer=1000;% number of simulation period
net.para.sInd=0;%period counter



 
ph=pi/net.para.side;% used for shifting of area to be Attractive!
t=linspace(0,1,net.para.side+1);% Generate linearly spaced vector( used for making polygon)
x1=sin(2*pi*t+ph); x1=(((x1-min(x1))/(max(x1)-min(x1)))-0.5)* net.para.length;% x coordinate of area(after normalising)
y1=-cos(2*pi*t+ph);y1=(((y1-min(y1))/(max(y1)-min(y1)))-0.5)*net.para.width;% y coordinate of area(after normalising)


  [x1,y1]=divider(x1,y1,10,1);% adding points between each pair of vertices

rs1=randsample(normalVect(rand(1,length(x1)),1,.5),length(x1));
figure;
x1=x1.*rs1;
y1=y1.*rs1;
[x1,y1]=divider(x1,y1,500);% adding points between each pair of vertices
x1(end+1)=x1(1);
y1(end+1)=y1(1);

plot(x1,y1,'r');hold on
rs=sqrt(randsample(normalVect(rand(1,length(x1)),.0,1),length(x1)));

plot(x1.*rs,y1.*rs,'.')






% give a vector or matrix for normalizing between tow number or 0 and 1
% vect=normalVect([1,-2,4,1,11,4,5])
% vect=normalVect([1,-2,4,1,11,4,5],-2,1)
% vect=normalVect([1,-2,4;11,4,5])

function vect=normalVect(varargin)
if nargin==0
    msgbox('give a 1 argoment atleast');
    vect=[];
elseif ~isnumeric(varargin{1})
    msgbox('abnormal argoment');
    vect=[];
else
   vect=varargin{1}; 
if nargin==1
    vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
elseif nargin==3 ||  nargin>3
    if varargin{2}~=varargin{3}
        vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
        vMin=min([varargin{2},varargin{3}]);
        vMax=max([varargin{2},varargin{3}]);
        vDif=vMax-vMin;
        vect=(vect*vDif)+vMin;
        if nargin>3
            msgbox('we used 3 argoment only');
        end
    else
        vect=(vect-min(vect(:)))./(max(vect(:))-min(vect(:)));
        msgbox('abnormal limits: normaled 0: 1');
    end
end

end
end





%version 1.1
%isFair=1 adding points is based on distance isFair=2 adding of points based on edges
function   [X,Y]=divider(X,Y,n,isFair)



try
    if isFair==1
    else
        isFair=0 ;
    end
catch
    
isFair=0;
end
switch isFair
    case 0
T1=0;
T2=0;

for i=1:length(X)-1
    temp1=0;
    temp1=linspace(X(i),X(i+1),n+2)  ;
    
    T1=[T1,temp1(1:end-1)];
end
T1=[T1,temp1(end)];
X=T1(2:end);


for i=1:length(Y)-1
    temp2=0;
    temp2=linspace(Y(i),Y(i+1),n+2)  ;
    
    T2=[T2,temp2(1:end-1)];
end
T2=[T2,temp2(end)];
Y=T2(2:end);
    case 1

dists= sqrt( (X(2:end)-X(1:end-1)).^2 + (Y(2:end)-Y(1:end-1)).^2) ;


n1=round((dists./min(dists))*n);
T1=0;
T2=0;

for i=1:length(X)-1
  
    temp1=0;
    temp1=linspace(X(i),X(i+1),n1(i)+2)  ;
    
    T1=[T1,temp1(1:end-1)];
end
T1=[T1,temp1(end)];
X=T1(2:end);


for i=1:length(Y)-1
    temp2=0;
    temp2=linspace(Y(i),Y(i+1),n1(i)+2)  ;
    
    T2=[T2,temp2(1:end-1)];
end
T2=[T2,temp2(end)];
Y=T2(2:end);
end


end

نمونه خروجی برای کد بالا(Output sample for the above code):