% 2.2.3 計算目標函數值
% calobjvalue.m函數的功能是實現目標函數的計算
%遺傳算法子程序
%Name: calobjvalue.m
%實現目標函數的計算
function [objvalue]=calobjvalue(pop)
[px,py]=size(pop);
q=imread('pic3.bmp');
q1=imread('pic.bmp');
[py1,px1]=size(q1);
ax=decodechrom(pop,1,8);%將pop每行轉化成十進制數
ay=decodechrom(pop,9,8);
for i=1:px;
    if ax(i,1)==0
       ax(i,1)=ax(i,1)+1;
   end
   if ay(i,1)==0
       ay(i,1)=ay(i,1)+1;
   end
    rect=[ax(i,1) ay(i,1) px1-1 py1-1];
    q2=imcrop(q,rect);
    r(i,1)=corr2(q1,q2);
    objvalue(i,1)=[(r(i,1)+1)/2]^2;%計算目標函數值
end