Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  594 1756 (N71E17 L=108)  1236    128        2 ( 549/1709- 640/1797)  28.3 E27-E07 N66-N75
      2 1204 1716 (N42W40 L=165)   377     54        1 (1181/1695-1229/1736)   0.0 W36-W42 N39-N43
      3 1590 1611 (N22W60 L=186)   546    122        0 (1548/1561-1619/1654)  74.1 W57-W63 N19-N26
      4 1281 1621 (N34W38 L=163)   332    123        0 (1219/1605-1337/1643) -39.8 W34-W40 N30-N36
      5 1427 1463 (N20W38 L=164)   552    263        0 (1318/1412-1534/1494) -42.3 W32-W43 N14-N25
      6  347 1301 (N39E44 L= 82)   443     83        2 ( 319/1261- 372/1334) -45.0 E47-E41 N36-N42
      7  807 1253 (N25E07 L=119)   608    197        1 ( 730/1222- 906/1298) -38.7 E11-E01 N21-N30
      8  729 1020 (N14E17 L=108)   214    106        0 ( 702/ 982- 758/1063)  31.0 E19-E14 N12-N15
      9  530  941 (N14E32 L= 94)   639    147        1 ( 464/ 921- 601/ 968) -24.0 E36-E26 N11-N16
     10  364  909 (N15E46 L= 80)   315     88        1 ( 338/ 878- 387/ 946)  53.1 E48-E43 N13-N17
     11  768  908 (N07E17 L=108)   293    101        3 ( 740/ 858- 798/ 941) -90.0 E18-E15 N03-N09
     12  208  897 (N17E62 L= 64)   530     58        1 ( 185/ 864- 231/ 929) -63.4 E65-E59 N14-N18
     13 1301  782 (S14W10 L=135)   281     74        0 (1284/ 742-1314/ 822)  78.7 W08-W10 S16-S11
     14  280  755 (N07E58 L= 68)   818    122        1 ( 229/ 729- 338/ 783) -16.7 E62-E52 N05-N09
     15  724  608 (S10E29 L= 97)   257    104        0 ( 699/ 561- 737/ 654)  80.5 E29-E27 S12-S06
     16 1586  620 (S33W26 L=152)   425     76        1 (1558/ 581-1604/ 652)  76.0 W23-W27 S36-S30
     17  524  522 (S11E46 L= 80)   438     68        1 ( 502/ 490- 550/ 559)  63.4 E47-E44 S12-S08
     18  528  450 (S15E49 L= 76)   552     83        2 ( 496/ 410- 555/ 489)  71.6 E51-E46 S18-S12
     19  508  365 (S21E57 L= 68)   640    112        0 ( 454/ 348- 568/ 391)   3.8 E65-E50 S21-S19
  

    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