function [con_i,con_j]=reference(square);

ia=54;
ja=90;
ib=77;
jb=565;
ic=331;
jc=81;
id=348;
jd=553;

Im_out=zeros(475,275);
new_cor=[1 1 1 1;
         1 475 275 1;
         275 1 275 1;
         275 475 275*475 1];
ori_cor_i=[ia;ib;ic;id];
ori_cor_j=[ja;jb;jc;jd];

con_i=new_cor\ori_cor_i;
con_j=new_cor\ori_cor_j;