Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1522 1692 (N28W57 L= 79)  2688    318        7 (1398/1644-1682/1736) -19.8 W45-W72 N22-N31
      2 1611 1533 (N17W53 L= 75)   774    171        1 (1534/1484-1688/1574) -48.4 W48-W57 N12-N21
      3  245 1223 (N32E52 L=330)  1419    256        2 ( 156/1121- 325/1296) -31.3 E67-E43 N23-N37
      4 1836 1173 (S08W58 L= 80)  2265    335        7 (1769/1100-1907/1243) -39.3 W52-W66 S12-S02
      5 1694 1158 (S05W44 L= 66)   471    140        3 (1669/1081-1724/1201)  35.5 W40-W47 S08-S02
      6 1407 1104 (S01W23 L= 45)   626    229        1 (1360/1038-1482/1157)   0.0 W19-W28 S04-N01
      7  256 1057 (N21E51 L=331)   355     60        1 ( 223/1038- 280/1081)   9.5 E54-E48 N19-N22
      8 1141  918 (S05W04 L= 26)  2006    496        1 ( 994/ 750-1277/1064) -66.8 E01-W08 S17-N04
      9 1546 1002 (S10W30 L= 52)   284     64        1 (1513/ 972-1568/1034)  26.6 W27-W31 S12-S08
     10 1444 1000 (S08W23 L= 45)   419     99        3 (1401/ 961-1483/1027)  -9.5 W20-W26 S10-S05
     11  565  968 (N10E28 L=354)   518    127        0 ( 491/ 949- 633/ 988) -24.0 E32-E23 N07-N12
     12 1603  743 (S27W30 L= 52)  1406    246        2 (1520/ 692-1704/ 819) -26.6 W24-W38 S30-S20
     13 1404  637 (S28W14 L= 36)  1453    289        3 (1335/ 513-1456/ 776)  86.8 W11-W16 S36-S18
     14  765  748 (S06E21 L=  1)   198     36        2 ( 740/ 726- 784/ 766) -63.4 E22-E19 S07-S05
     15  664  468 (S21E37 L=345)  3740    718        4 ( 441/ 281- 840/ 698)  18.4 E72-E18 S28-S10
     16  692  635 (S11E29 L=353)   873    194        3 ( 611/ 579- 784/ 696)  -4.8 E35-E23 S13-S07
     17 1394  459 (S40W10 L= 32)   842    144        2 (1335/ 396-1439/ 508)  26.6 W04-W14 S43-S37
     18 1252  317 (S47E05 L= 17)  1630    152        4 (1176/ 276-1335/ 358)  -8.7 E11-W02 S49-S43
  

    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