Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1615 1718 (N24W73 L= 91)  1073     85        3 (1569/1700-1666/1734) -19.7 W65-W82 N21-N26
      2 1891 1281 (S07W69 L= 87)   599     54        3 (1867/1247-1913/1310)  26.6 W65-W74 S08-S05
      3  323 1229 (N35E44 L=334)   793    176        1 ( 243/1172- 390/1264) -26.6 E53-E37 N30-N39
      4  751 1082 (N16E14 L=  4)   415    140        0 ( 699/1039- 783/1127) -63.4 E16-E11 N13-N19
      5  128 1051 (N28E70 L=308)   995    103        1 (  75/1020- 170/1073)  -5.0 E85-E62 N25-N29
      6  718 1022 (N13E17 L=  1)   222     52        0 ( 690/1001- 742/1042) -45.0 E19-E15 N11-N14
      7 1588  996 (S11W32 L= 50)   300     70        0 (1555/ 963-1608/1032)  31.0 W28-W33 S12-S09
      8 1633  906 (S18W33 L= 51)   585    320        0 (1538/ 833-1711/1018)  13.2 W24-W41 S20-S14
      9 1789  953 (S21W47 L= 65)   349     49        0 (1765/ 933-1817/ 970)   0.0 W44-W50 S21-S19
     10 1401  934 (S09W18 L= 36)   263    107        0 (1362/ 917-1454/ 963) -45.0 W16-W21 S11-S06
     11 1330  810 (S14W11 L= 29)   725    244        1 (1236/ 761-1376/ 899)   6.3 W05-W14 S16-S09
     12 1491  815 (S19W21 L= 39)   393     90        0 (1470/ 769-1519/ 863) -90.0 W20-W22 S22-S15
     13 1784  833 (S28W45 L= 63)   522    128        0 (1744/ 787-1835/ 861)  11.3 W40-W50 S29-S26
     14  882  718 (S07E16 L=  2)   919    246        0 ( 849/ 603- 922/ 822)  42.3 E21-E10 S13-S02
     15  616  383 (S21E46 L=332)  3307    457        5 ( 504/ 218- 780/ 526)  22.1 E65-E28 S32-S11
     16  381  482 (S10E60 L=318)   596     92        2 ( 357/ 445- 412/ 511)  20.6 E64-E56 S11-S07
     17 1418  459 (S39W08 L= 26)   506     86        0 (1401/ 409-1438/ 501)  82.9 W06-W09 S43-S35
     18 1304  325 (S46E05 L= 13)  1621    248        3 (1227/ 248-1391/ 389)  18.4 E14-W04 S50-S42
  

    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