IDL> fits_read, 'polar_45d_1s_14_redbias.fits', p1, header IDL> fits_read, 'polar_45d_1s_15_redbias.fits', p2, header IDL> fits_read, 'polar_45d_1s_16_redbias.fits', p3, header ----------------------------------------------------------------------------- IDL> shift_div_graph, p2, p1, pn, .95, 1.05, xrange=[0,900], yrange=[400,1299] enter x shift: .65 enter y shift: 0 Mean, Max, Min, StdDev 0.995505 2.07908 0.316830 0.0270855 IDL> shift_div_graph, p2, p1, pn, .95, 1.05, xrange=[0,1299], yrange=[0,1299], scale =.5 enter x shift: .6 enter y shift: 0 Mean, Max, Min, StdDev 0.998256 206.343 -181.895 0.332756 ----------------------------------------------------------------------------- IDL> shift_div_graph, p3, p1, pn, .9, 1.1, xrange=[0,1299], yrange=[1,1299], scale =1 Mean, Max, Min, StdDev 0.996595 84.7382 -9.22679 0.207563 enter x shift: -.45 enter y shift: .1 ----------------------------------- IDL> p2_f = shift_im_cbc(p2, .6, .1) IDL> p3_f = shift_im_cbc(p3, -.45, .1) IDL> fits_write, 'polar_45d_1s_15_shifted.fits', p2 IDL> fits_write, 'polar_45d_1s_16_shifted.fits', p3 ------------------------------------------------------------- IDL> fits_read, 'polar_45d_1s_14_redbias.fits', p45_0 IDL> fits_read, 'polar_45d_1s_15_shifted.fits', p45_1 IDL> fits_read, 'polar_45d_1s_16_shifted.fits', p45_2 IDL> p45 = (float(p45_0) + float(p45_1) + float(p45_2)) / 3.0 IDL> fits_write, 'polar_45d_1s_mean_aligned.fits', p45