Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  587 1684 (N20E49 L= 23)   421     80        1 ( 565/1647- 623/1710) -78.7 E51-E47 N17-N22
      2  418 1669 (N16E64 L=  8)   950    109        4 ( 381/1635- 467/1709)  51.3 E69-E59 N13-N19
      3  516 1649 (N17E53 L= 19)   465    110        0 ( 490/1621- 563/1686) -14.0 E55-E47 N15-N19
      4  588 1615 (N16E45 L= 27)   274     41        1 ( 566/1592- 608/1634)  90.0 E46-E43 N14-N17
      5 1223 1508 (N26E02 L= 70)  1017    300        1 (1109/1398-1304/1564) -15.9 E09-W05 N21-N29
      6 1617 1459 (N37W26 L= 98)   550    111        0 (1601/1400-1629/1516)  90.0 W25-W27 N32-N41
      7  392 1451 (N02E52 L= 20)  1632    237        4 ( 292/1397- 502/1488)  20.6 E60-E43 S02-N05
      8  387 1273 (S08E46 L= 26)   257     57        1 ( 363/1251- 406/1301)  71.6 E47-E43 S09-S06
      9 1652 1146 (N18W34 L=106)   243     53        0 (1619/1127-1683/1161)  45.0 W31-W36 N15-N19
     10  234 1068 (S23E56 L= 16)  2560    428        7 (  77/1001- 400/1152)  14.3 E84-E41 S28-S15
     11 1812 1118 (N22W49 L=121)   610    158        0 (1746/1086-1855/1147)   9.5 W41-W53 N20-N24
     12 1742 1072 (N17W43 L=115)   395    155        0 (1700/1004-1763/1141) -36.9 W38-W46 N13-N19
     13  575 1020 (S18E26 L= 46)  1191    382        0 ( 387/ 986- 707/1072)  24.4 E40-E18 S22-S12
     14 1683  939 (N07W42 L=114)   382     57        1 (1655/ 909-1716/ 962)  53.1 W39-W43 N05-N09
     15 1860  868 (N10W60 L=132)   422     72        2 (1832/ 838-1889/ 898)  -6.3 W56-W65 N09-N11
     16 1775  796 (N03W54 L=126)   419     63        1 (1752/ 760-1789/ 830) -63.4 W51-W56 N01-N05
     17 1252  694 (S18W23 L= 95)   354     93        1 (1198/ 680-1295/ 711)  18.4 W19-W25 S20-S16
     18 1617  475 (S18W58 L=130)  1136    172        5 (1574/ 413-1663/ 555) -66.0 W54-W61 S22-S13
  

    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