PRO arpore,band,date,px1,px2,py1,py2,ardi ;Development of sizes and brightnesse of a pore (IR obs. 1997) ;INPUT: BAND = 1 (IR) or 2 (WL) ; DATE = 23 or 25 ; PX1,...,PY2 - limits of the pore's field ;OUTPUT: ARDI - a 2D array ([Area, Deff, Imin],[frames]), ; where the total area and absolute intensity minimum of all ; fragments of the observed pore is given for each frame. ;SUBROUTINES: AR2I.PRO, LABEL_REGION.PRO, HISTOGRAM.PRO ;AR2I.pro contains parameters: scale ("/pix), minimum size, QR subfield... if band eq 1 then fil1='/net/obus/scratch3/sob/obs97/irdata/i' $ else fil1='/net/obus/scratch3/sob/obs97/wldata/w' if date eq 23 then begin lista=indgen(104)+7 fil2='fi23' arr=intarr(408,320,118) ardi=fltarr(3,208) ;output array x1=10 x2=397 y1=10 y2=309 endif else begin lista=indgen(136)+7 fil2='fi25' arr=intarr(388,300,150) ardi=fltarr(3,272) ;output array x1=10 x2=377 y1=10 y2=289 endelse n=n_elements(lista) ;Processing 1st part openr,1,fil1+fil2+'a.bin' print,' reading ',fil1+fil2+'a.bin' readu,1,arr close,1 imgs=arr(x1:x2,y1:y2,lista) ;Measurement of size and brightness for i=0,n-1 do begin AR2I,imgs(*,*,i),date,ob,px1,px2,py1,py2 ardi(0,i)=total(ob(0,*),2) ;area ardi(1,i)=0.1658*sqrt(4*ardi(0,i)/!pi) ;Deff ardi(2,i)=min(ob(2,*)) ;Imin endfor ;Processing 2nd part openr,1,fil1+fil2+'b.bin' print,' reading ',fil1+fil2+'b.bin' readu,1,arr close,1 imgs=arr(x1:x2,y1:y2,lista) ;Measurement of size and brightness for i=0,n-1 do begin AR2I,imgs(*,*,i),date,ob,px1,px2,py1,py2 ardi(0,i+n)=total(ob(0,*),2) ;area ardi(1,i+n)=0.1658*sqrt(4*ardi(0,i+n)/!pi) ;Deff ardi(2,i+n)=min(ob(2,*)) ;Imin endfor END