Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1078 1861 (N55W49 L=272)   841     90        0 (1046/1822-1116/1895)   6.3 W40-W58 N53-N56
      2 1701 1582 (N15W69 L=292)  1814    192        1 (1626/1509-1765/1637) -70.0 W64-W74 N09-N20
      3 1526 1583 (N22W52 L=275)   356    138        0 (1471/1533-1566/1612) -74.1 W49-W53 N17-N24
      4  909 1463 (N34W06 L=229)  1515    428        3 ( 724/1421-1130/1499) -30.3 E06-W18 N25-N39
      5 1637 1415 (N09W51 L=274)   567    123        1 (1599/1368-1663/1477)  32.0 W46-W54 N07-N12
      6 1104 1260 (N17W11 L=234)   618    271        0 ( 968/1234-1215/1303) -18.4 W03-W18 N14-N20
      7 1312 1219 (N09W21 L=245)   629    203        0 (1233/1144-1366/1278) -74.7 W18-W23 N03-N14
      8  749 1177 (N22E12 L=211)   683    206        2 ( 657/1143- 841/1209) -18.4 E18-E06 N19-N23
      9  272 1132 (N30E51 L=172)   344     60        0 ( 254/1097- 287/1159) -63.4 E52-E49 N28-N32
     10 1806 1112 (S12W53 L=277)   517    100        1 (1767/1081-1850/1145) -26.6 W49-W58 S14-S09
     11  756 1059 (N15E14 L=209)   536    200        0 ( 671/1020- 840/1084) -20.0 E20-E09 N11-N16
     12  247  893 (N17E57 L=166)   898    216        0 ( 167/ 840- 334/ 938) -26.6 E66-E49 N12-N20
     13 1571  817 (S21W27 L=250)   378     92        0 (1537/ 775-1609/ 853)  18.4 W24-W30 S22-S19
     14 1750  588 (S44W43 L=266)  2723    315        3 (1646/ 507-1837/ 723)  20.0 W30-W54 S47-S37
     15  760  436 (S21E34 L=189)   668    168        1 ( 737/ 350- 777/ 504)  45.0 E38-E29 S26-S16
     16 1471  393 (S48W12 L=235)   724    145        1 (1414/ 349-1540/ 428)   8.7 W06-W19 S49-S45
     17 1007  239 (S43E31 L=193)   424     61        0 ( 978/ 214-1036/ 260)   8.1 E34-E27 S43-S41
     18 1171  215 (S51E22 L=202)  1483    255        2 (1055/ 173-1284/ 248)  -7.1 E34-E10 S54-S48
  

    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