Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  526 1716 (N47E58 L= 61)   492     51        1 ( 512/1684- 547/1738) -26.6 E63-E53 N44-N49
      2 1704 1486 (N34W58 L=177)  2390    361        1 (1578/1356-1831/1582) -38.7 W47-W70 N24-N41
      3 1394 1417 (N29W25 L=144)   422    105        1 (1366/1364-1418/1463) -74.1 W23-W26 N25-N32
      4  824 1385 (N25E15 L=104)   462    139        0 ( 764/1347- 881/1421)  35.5 E18-E11 N22-N27
      5 1639 1367 (N26W45 L=164)   440     91        2 (1599/1338-1670/1401)   8.1 W41-W48 N23-N28
      6 1656 1041 (N05W43 L=161)   446    141        1 (1637/ 984-1677/1108)  82.9 W41-W44 N01-N09
      7  209  846 (S12E62 L= 56)  1964    210        5 ( 179/ 754- 243/ 944)  59.7 E68-E57 S18-S06
      8 1749  856 (S06W52 L=171)   546      9        1 (1724/ 832-1775/ 878)  14.0 W49-W55 S07-S04
      9  335  758 (S17E49 L= 69)   573     63        3 ( 312/ 720- 365/ 788)  76.0 E51-E46 S19-S15
     10  492  690 (S20E36 L= 82)  1229    216        2 ( 415/ 617- 559/ 768)  35.5 E44-E30 S25-S15
     11  775  654 (S21E15 L=103)   816    230        0 ( 698/ 602- 874/ 714) -30.3 E20-E08 S24-S17
     12 1666  625 (S21W50 L=168)   500     82        1 (1629/ 595-1708/ 656)  26.6 W46-W54 S23-S19
     13  299  578 (S30E60 L= 58)  1132    132        0 ( 242/ 547- 341/ 626)  13.4 E73-E52 S32-S26
     14 1043  555 (S27W03 L=122)   401    105        0 (1016/ 510-1084/ 597)  56.3 W01-W05 S30-S24
     15 1747  569 (S25W62 L=181)   533     77        0 (1709/ 551-1783/ 583)   5.7 W58-W69 S26-S24
     16 1379  479 (S32W29 L=147)   288     84        0 (1346/ 468-1427/ 498)  -7.1 W25-W33 S32-S30
     17 1126  465 (S34W10 L=128)   604    114        1 (1065/ 436-1169/ 493) -14.0 W05-W13 S35-S31
  

    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