Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1839 (N34E59 L=176)   839     90        4 ( 655/1802- 711/1878) -20.6 E69-E53 N31-N37
      2  403 1643 (N14E65 L=170)   900     83        3 ( 378/1607- 432/1680) -26.6 E72-E61 N11-N16
      3  488 1635 (N15E55 L=179)  1325    123        3 ( 436/1587- 535/1670)  45.0 E59-E51 N12-N18
      4  603 1643 (N18E46 L=188)   349     41        1 ( 579/1625- 626/1667)   0.0 E49-E44 N16-N19
      5  355 1601 (N11E67 L=167)   691     73        2 ( 328/1577- 389/1631)  -8.7 E74-E61 N09-N12
      6  536 1540 (N10E46 L=189)   370     54        2 ( 518/1516- 560/1568) -18.4 E48-E42 N08-N12
      7  659 1523 (N12E36 L=198)   399     92        2 ( 637/1482- 684/1566) -33.7 E39-E33 N09-N14
      8 1681 1451 (N40W33 L=268)   624    122        0 (1634/1396-1716/1497) -24.0 W28-W37 N36-N41
      9 1627 1334 (N29W29 L=264)   361    100        0 (1614/1288-1648/1382) -51.3 W27-W31 N26-N32
     10  464 1147 (S14E37 L=197)   261     61        1 ( 449/1119- 483/1177) -45.0 E38-E35 S15-S12
     11  672  979 (S19E19 L=215)   253     90        1 ( 630/ 961- 718/ 991)  31.0 E21-E16 S19-S16
     12 1844  966 (N15W56 L=291)   371     40        0 (1821/ 946-1868/ 986)   0.0 W53-W59 N14-N16
     13 1900  950 (N17W65 L=299)   608     88        0 (1864/ 922-1927/ 986)  36.9 W59-W68 N14-N20
     14  504  904 (S27E29 L=205)   364     72        0 ( 460/ 883- 547/ 922)   8.1 E32-E25 S28-S25
     15 1688  595 (S09W57 L=292)   688    153        0 (1626/ 550-1734/ 654)  66.0 W53-W59 S13-S04
     16 1115  582 (S29W19 L=254)   305     76        0 (1080/ 565-1159/ 601)   0.0 W16-W22 S29-S27
     17 1601  437 (S20W61 L=295)   657    117        0 (1559/ 396-1642/ 481)  74.1 W58-W63 S23-S16
  

    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