function Iout = cropCenter(Iin)
	h=(size(Iin,1)+1)/4; w=(size(Iin,2)+1)/4;
	Iout = Iin(h:3*h-1, w:3*w-1); %cut the center rigion
return
