Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1209 1861 (N49W54 L=260)   743     73        1 (1187/1816-1226/1896)  17.4 W47-W63 N45-N52
      2  639 1618 (N51E11 L=194)   465     61        1 ( 611/1583- 666/1647)  68.2 E14-E08 N48-N53
      3  823 1608 (N44W04 L=209)   242     51        0 ( 796/1591- 848/1624) -53.1 W01-W05 N42-N46
      4  910 1475 (N33W05 L=211)   219     52        0 ( 895/1445- 926/1501)  45.0 W03-W07 N31-N34
      5  989 1240 (N17W03 L=209)   366    210        0 ( 920/1163-1042/1306) -74.7 W00-W04 N11-N22
      6 1908  932 (S26W65 L=270)  2772    198        5 (1876/ 849-1944/1010)  71.6 W59-W72 S31-S20
      7 1557  338 (S59W26 L=231)  1523    174        3 (1472/ 318-1635/ 368) -17.1 W14-W40 S64-S53
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y