Curves of Flux, Circular Polarization, and Spectrum Index from a Partial Area

In order to obtain curves of flux, circular polarization degree, and spectrum index from a partial area of a series of images,
IDL> abox=[23,35,33,41] $<$CR$>$
IDL> fi17=norh_tb2flux(dataa,indexa,abox=abox) $<$CR$>$
IDL> fv17=norh_tb2flux(datas,indexs,abox=abox) $<$CR$>$
IDL> fi34=norh_tb2flux(dataz,indexz,abox=abox) $<$CR$>$
IDL> norh_polariz,indexa,fi17,indexs,fv17,indexp,polariz,mvdp $<$CR$>$
IDL> norh_alpha,indexa,fi17,indexz,fi34,indexal,alpha,mvda $<$CR$>$
Plot as, say,
IDL> utplot,indexp,polariz $<$CR$>$
IDL> utplot,indexal,alpha $<$CR$>$