function [img, featurelist] = mark_features(imname, datname)
% [img, size]: image matrix and size vector(width,height, max_gray)
img = readp5(imname);
cla;
image(img);
colormap(gray(256));
axis('image');
grid on;
featurelist = read_motion_dat(datname)+1;

hold on;
for i=1:size(featurelist,1);
  plot(featurelist(i,2), featurelist(i,1),'r+');
end;
zoom on;
return







