calculateEdgeCost.m

function [edgeResponse, commonSize] = calculateEdgeCost(edgeImg, labels,labelIndices, labOne, labTwo)
edgeResponse = 0;
commonSize = 0;
% Find common border..


[rC,cC] = size(labels);
temp = labelIndices{1,labOne};
rowArr = temp(:,1);
colArr = temp(:,2);
[realSize,~] = size(rowArr);

totalResp  = 0;
commonSize = 0;
for i=1:realSize
    xCor = rowArr(i,1);
    yCor = colArr(i,1);
    
    if xCor + 1 <= rC
        if labels(xCor+1,yCor) == labTwo % Common border!
            commonSize = commonSize + 1;
            if edgeImg(xCor,yCor) ~= edgeImg(xCor+1,yCor)
                totalResp = totalResp + 1;
            end
        end
    end
    
    if xCor - 1 > 0
         if labels(xCor-1,yCor) == labTwo % Common border!
            commonSize = commonSize + 1;
            if edgeImg(xCor,yCor) ~= edgeImg(xCor-1,yCor)
                totalResp = totalResp + 1;
            end
        end
    end
    
    if yCor + 1 <= cC
         if labels(xCor,yCor+1) == labTwo % Common border!
            commonSize = commonSize + 1;
            if edgeImg(xCor,yCor) ~= edgeImg(xCor,yCor+1)
                totalResp = totalResp + 1;
            end
        end
    end
    
    if yCor - 1 > 0
        if labels(xCor,yCor-1) == labTwo % Common border!
            commonSize = commonSize + 1;
            if edgeImg(xCor,yCor) ~= edgeImg(xCor,yCor-1)
                totalResp = totalResp + 1;
            end
        end
    end
end

if commonSize ~=0
    edgeResponse = totalResp / commonSize;
end

end