Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1765 (N30E54 L=312)  1436    141        2 ( 530/1709- 637/1815) -16.3 E69-E45 N26-N33
      2  838 1572 (N21E25 L=342)   743    255        1 ( 711/1550- 948/1602)  10.6 E33-E17 N18-N24
      3 1524 1378 (N26W23 L= 30)  2321    683        2 (1253/1181-1717/1556) -14.4 W01-W40 N18-N30
      4  144  964 (S26E69 L=298)   836    141        2 ( 107/ 908- 177/1019)  30.3 E77-E63 S29-S22
      5 1385  982 (S01W22 L= 28)   197     48        0 (1357/ 962-1410/1001)  45.0 W20-W23 S02-N00
      6  370  902 (S27E42 L=325)   886    194        2 ( 302/ 848- 425/ 990)  42.3 E48-E37 S30-S20
      7  203  925 (S28E60 L=307)   350     69        0 ( 179/ 909- 238/ 953)  -8.1 E63-E55 S28-S26
      8  480  740 (S35E31 L=335)  1153    306        1 ( 372/ 657- 600/ 805)  35.0 E42-E22 S42-S28
  

    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