infil = '$MDI_CAL_TXT/distgrid617.gcross
if (n_elements(mat) eq 0) then mat = rd_tfile(infil, 11, 5, /conv)
;
imgarr = mat(0,*)
x = mat(3,*)
y = mat(4,*)
;
u_imgarr = imgarr(uniq(imgarr))
nimg = n_elements(u_imgarr)
;
plot, fltarr(2), xrange=[0,1000], yrange=[0,1000], xtit='X Pos', ytit='Y pos'
for i=0,nimg-1 do begin
    ss = where(imgarr eq u_imgarr(i))
    oplot, x(ss), y(ss), psym=1
end
pause
;
!p.multi = [0,2,2]
for i=0,nimg-1 do begin
    ss = where(imgarr eq u_imgarr(i))
    plot, x(ss), y(ss), xrange=[0,1000], yrange=[0,1000], xtit='X Pos', ytit='Y pos', tit='Image ' + string(u_imgarr(i)), psym=2
    if ((i mod 4) eq 3) then pause
end
;

;
end