Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  752 1753 (N30E42 L=140)  2003    367        1 ( 683/1637- 862/1818) -86.2 E47-E34 N21-N36
      2 1190 1697 (N38E10 L=172)   840    168        2 (1131/1636-1233/1760) -26.6 E16-E04 N34-N42
      3  367 1508 (N06E58 L=124)   864    137        0 ( 323/1476- 413/1543)  45.0 E63-E54 N04-N09
      4 1728 1411 (N37W39 L=221)  1147    256        0 (1673/1287-1770/1515) -60.3 W33-W43 N29-N43
      5  408 1476 (N05E52 L=129)   354     54        1 ( 386/1452- 431/1498) -11.3 E54-E49 N03-N06
      6  429 1275 (S06E43 L=139)   330     64        1 ( 400/1247- 448/1302)  63.4 E44-E41 S08-S04
      7  234  974 (S27E56 L=126)  1836    428        4 ( 166/ 802- 275/1149)  63.4 E68-E52 S38-S16
      8 1863 1026 (N18W57 L=239)  1326    231        3 (1833/ 929-1911/1110) -59.0 W52-W65 N12-N22
      9  356  865 (S32E42 L=139)  1257    421        1 ( 312/ 722- 414/1040) -68.7 E46-E35 S40-S21
     10  479  956 (S24E33 L=149)   573    140        2 ( 448/ 921- 531/1016) -33.7 E35-E28 S26-S20
     11  691  925 (S21E17 L=164)   595    188        0 ( 613/ 853- 736/ 978) -21.8 E23-E13 S24-S18
     12 1460  810 (S07W31 L=212)   373     87        1 (1431/ 767-1487/ 848)  80.5 W28-W32 S10-S04
     13 1627  640 (S11W49 L=231)  1407    249        0 (1521/ 570-1724/ 687)  45.0 W44-W56 S17-S06
     14 1328  501 (S28W36 L=217)   650    151        0 (1255/ 478-1392/ 534)  32.0 W30-W38 S30-S24
     15 1513  379 (S28W58 L=239)  1765    244        1 (1394/ 351-1626/ 407)  15.6 W47-W72 S32-S24
     16  612  335 (S64E13 L=169)  1744    181        1 ( 524/ 295- 693/ 381)   2.3 E26-E01 S67-S60
     17 1194  328 (S42W39 L=220)  1548    228        5 (1066/ 296-1285/ 363)  28.1 W28-W46 S46-S37
  

    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