function [ T ] = getTemplate ( I ) imagesc(I); zoom on; fprintf (1, 'Zoom until desired region is displayed, then hit any key.\n'); pause; ax = round(axis); T = I(ax(3) : ax(4), ax(1) : ax(2));