Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1145 1744 (N46W37 L=158)   301     44        1 (1129/1715-1164/1773)  38.7 W34-W39 N44-N48
      2 1395 1582 (N28W42 L=163)  1041    197        1 (1319/1523-1484/1627)   0.0 W33-W50 N26-N29
      3 1518 1431 (N16W42 L=164)   515    177        1 (1456/1389-1596/1468)  -4.4 W36-W49 N15-N17
      4 1077 1359 (N25W12 L=133)   751    169        4 (1017/1311-1156/1390)  -6.3 W07-W16 N22-N26
      5  788 1238 (N25E09 L=113)   192     64        0 ( 762/1214- 815/1259)   0.0 E10-E06 N23-N25
      6 1787 1032 (S13W50 L=171)  4697    807       14 (1599/ 829-1886/1245)  18.9 W30-W68 S20-S03
      7 1056 1170 (N15W05 L=127)   237    122        0 (1016/1121-1083/1214) -81.9 W04-W06 N11-N18
      8 1183 1109 (N08W11 L=132)   244    111        0 (1151/1060-1218/1159)  26.6 W08-W14 N06-N10
      9  854 1038 (N12E09 L=112)   253    100        0 ( 805/1022- 889/1062) -26.6 E11-E07 N10-N13
     10  563  974 (N14E29 L= 93)   422     83        1 ( 527/ 951- 609/1008) -31.0 E31-E26 N12-N16
     11  785  899 (N06E16 L=105)   316    114        2 ( 766/ 837- 805/ 951)  71.6 E18-E15 N02-N08
     12  852  758 (S04E16 L=105)   516    187        2 ( 782/ 709- 942/ 801) -53.1 E19-E12 S08-N00
     13 1811  740 (S32W50 L=171)   616    114        0 (1759/ 715-1851/ 774)   4.8 W43-W55 S33-S31
     14 1527  686 (S26W23 L=144)   574    105        2 (1488/ 653-1561/ 723) -71.6 W20-W25 S28-S22
     15 1515  551 (S35W20 L=141)  1385    360        0 (1404/ 437-1611/ 701)  31.0 W09-W29 S39-S27
     16 1180  566 (S23E02 L=119)   465    177        0 (1113/ 518-1233/ 631)  24.0 E07-W02 S24-S20
     17 1669  486 (S46W36 L=157)   971    114        1 (1621/ 447-1712/ 531)   5.7 W29-W40 S49-S42
     18  992  425 (S28E18 L=103)   372    115        0 ( 963/ 401-1038/ 474) -80.5 E19-E15 S29-S23
     19  434  398 (S18E61 L= 61)   443     49        0 ( 420/ 368- 446/ 428)  45.0 E63-E58 S20-S15
  

    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