Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1541 1748 (N32W65 L=295)   636     46        0 (1516/1722-1569/1767)   0.0 W60-W70 N30-N33
      2 1313 1703 (N36W41 L=271)   504    146        1 (1264/1640-1362/1754) -90.0 W39-W42 N30-N40
      3  671 1617 (N48E10 L=220)   270     59        0 ( 653/1587- 697/1636)  63.4 E12-E07 N45-N49
      4 1594 1560 (N20W53 L=283)   925    238        4 (1534/1493-1647/1625)  83.7 W47-W56 N15-N24
      5  516 1578 (N49E26 L=204)   435     98        1 ( 493/1541- 539/1618) -81.9 E28-E23 N45-N52
      6 1077 1551 (N32W17 L=247)  1312    303        5 ( 914/1519-1209/1580) -20.2 W06-W25 N29-N36
      7 1513 1294 (N07W35 L=266)  1034    314        4 (1375/1250-1608/1343) -34.7 W27-W42 N02-N12
      8  462 1110 (N20E33 L=197)  1209    443        3 ( 293/1006- 618/1204) -38.0 E46-E23 N10-N28
      9  111 1077 (N24E72 L=159)   598     62        1 (  74/1057- 141/1092)   0.0 E82-E65 N22-N24
     10  214  912 (N13E58 L=172)   704    188        2 ( 188/ 840- 237/1014) -90.0 E60-E55 N08-N19
     11  856  954 (N03E11 L=219)   177     46        1 ( 831/ 927- 881/ 978)  45.0 E12-E09 N01-N03
     12 1268  877 (S10W11 L=241)   666    231        5 (1202/ 835-1347/ 943) -53.1 W06-W14 S13-S05
     13  782  903 (N01E16 L=214)   192     69        0 ( 762/ 876- 797/ 941)  56.3 E17-E15 N00-N03
     14 1799  716 (S34W50 L=280)   792    122        0 (1738/ 690-1839/ 758)   4.8 W42-W55 S35-S32
     15  843  687 (S12E18 L=212)   283    104        0 ( 792/ 668- 887/ 716)   0.0 E21-E14 S13-S11
     16  717  651 (S11E26 L=204)   245    106        0 ( 700/ 604- 737/ 694)  68.2 E27-E23 S14-S09
     17  719  394 (S27E37 L=194)  1794    411        4 ( 573/ 320- 911/ 493)   5.4 E51-E19 S29-S24
     18 1129  299 (S45E14 L=217)  1628    269        2 (1020/ 242-1249/ 387)   9.9 E24-E01 S48-S40
  

    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