%-- re-map small mask to big mask

function out_mask = remap_mask(in_mask,m_cnt,I)

    [y_max, x_max]=size(I);

    M=size(in_mask,1);

    % center

    m_cnt_x=m_cnt(1);    

    m_cnt_y=m_cnt(2);

    % top left corner

    x1=m_cnt_x-M/2;

    y1=m_cnt_y-M/2;

    % bottom right corner

    x4=m_cnt_x+M/2-1;

    y4=m_cnt_y+M/2-1;

    out_mask=zeros(x_max,y_max);

    out_mask(y1:y4,x1:x4)=in_mask;

end