; File: image_interp
; Author: Erik Brisson
a = bytarr(192,192)
openr, 3, 'dat/people.dat'
readu, 3, a
close, 3

window, 0, retain=2, xsize=192, ysize=192
tv, a

p = float(a)
x = findgen(400) * (192.0 / 400.0) 
y = findgen(400) * (192.0 / 400.0)
r = interpolate(p, x, y, /grid)
ar = byte(r)
print, "min(r), max(r) = ", min(r), max(r)
print, "min(ar), max(ar) = ", min(ar), max(ar)

window, 1, retain=2, xsize=400, ysize=400
tv, ar
loadct

end