Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1156 1859 (N59W50 L=213)   504     63        0 (1134/1833-1185/1876)   4.4 W43-W56 N57-N60
      2 1416 1400 (N21W34 L=196)   636    177        0 (1342/1380-1497/1421) -21.8 W29-W39 N18-N23
      3 1741 1083 (S05W49 L=211)   720    159        1 (1717/1011-1760/1155)  83.7 W46-W51 S09-N00
      4 1638 1015 (S06W39 L=201)   715    222        0 (1594/ 931-1684/1121) -71.6 W36-W41 S12-N00
      5 1841  947 (S16W57 L=220)   605     90        2 (1807/ 905-1865/ 985) -78.7 W53-W60 S18-S13
      6 1517  792 (S17W27 L=189)   608    217        0 (1474/ 715-1538/ 892)  59.0 W22-W29 S21-S11
      7 1728  702 (S29W45 L=207)   546    119        0 (1691/ 659-1770/ 751)  29.7 W40-W48 S30-S26
      8 1279  697 (S18W09 L=172)   323    147        0 (1220/ 673-1341/ 732)   0.0 W05-W13 S18-S16
      9  383  554 (S12E55 L=107)  1227    255        2 ( 343/ 459- 425/ 637)  69.0 E61-E50 S19-S06
     10 1695  534 (S41W44 L=207)   762    101        0 (1671/ 485-1726/ 578)  81.9 W41-W47 S44-S37
     11  633  558 (S15E34 L=128)   502    156        0 ( 557/ 542- 707/ 572)  -5.7 E39-E29 S16-S13
     12 1192  396 (S37E02 L=161)  3495    591        7 ( 989/ 232-1385/ 551)  22.8 E24-W14 S47-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