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>