Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  488 1617 (N18E51 L= 17)   610    112        2 ( 460/1579- 517/1667) -40.6 E53-E46 N15-N21
      2  593 1631 (N20E43 L= 24)   332     44        2 ( 570/1607- 616/1658)  56.3 E44-E41 N18-N22
      3  769 1505 (N15E26 L= 41)   356    151        0 ( 736/1436- 808/1572)  90.0 E26-E25 N10-N20
      4  369 1370 (N01E50 L= 18)   730    170        2 ( 285/1347- 451/1389)   4.1 E57-E42 S01-N02
      5 1759 1319 (N29W43 L=111)   633    131        1 (1711/1264-1813/1353)  36.9 W38-W50 N23-N33
      6  764 1290 (N03E20 L= 47)   201     98        0 ( 740/1248- 797/1329)  78.7 E21-E19 N00-N05
      7  374 1199 (S09E45 L= 23)   673    197        2 ( 317/1119- 435/1239)   0.0 E49-E39 S14-S06
      8 1280 1062 (N01W14 L= 81)   223    120        0 (1225/1030-1320/1106)  -8.1 W09-W16 N00-N01
      9  333 1069 (S17E46 L= 21)   282     75        3 ( 297/1046- 360/1100)  31.0 E48-E43 S18-S15
     10 1657  994 (N06W39 L=106)   437     75        2 (1608/ 970-1701/1010)   8.1 W35-W42 N04-N07
     11  848  969 (S14E09 L= 58)   368    132        1 ( 784/ 947- 915/ 994)   7.1 E12-E04 S14-S12
     12  387  923 (S25E40 L= 28)   809    268        1 ( 301/ 832- 493/ 989)  38.7 E47-E31 S31-S19
     13 1531  774 (S09W36 L=103)   248     90        0 (1501/ 741-1574/ 806)  59.0 W34-W38 S11-S06
     14 1335  703 (S18W26 L= 93)   486    226        1 (1228/ 665-1441/ 750)  30.3 W20-W32 S22-S13
     15 1584  694 (S12W43 L=110)   361     88        1 (1544/ 675-1622/ 731)  53.1 W40-W46 S13-S09
     16 1680  544 (S17W58 L=125)   436     71        2 (1650/ 513-1713/ 569)   0.0 W54-W61 S18-S15
     17 1408  455 (S30W42 L=109)   464    105        2 (1349/ 438-1461/ 474)  33.7 W38-W45 S31-S27
     18 1575  389 (S28W60 L=128)   465     89        1 (1538/ 350-1608/ 419) -10.3 W54-W65 S29-S25
  

    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