Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  892 1826 (N44E31 L=162)   565     81        0 ( 849/1811- 938/1843)  15.9 E34-E27 N42-N46
      2  689 1713 (N31E39 L=154)   927    155        1 ( 625/1653- 736/1737) -17.1 E45-E32 N27-N33
      3  578 1658 (N26E45 L=148)   661    100        2 ( 550/1606- 603/1706) -63.4 E47-E41 N21-N29
      4 1673 1563 (N42W42 L=235)   908     73        2 (1642/1527-1708/1602) -26.6 W37-W45 N39-N45
      5  264 1519 (N13E69 L=124)  2425    130        4 ( 218/1460- 299/1575)  90.0 E76-E63 N09-N17
      6  327 1505 (N13E60 L=133)   445     62        1 ( 300/1474- 353/1535)  45.0 E61-E57 N10-N14
      7 1513 1497 (N32W24 L=218)   307     95        1 (1479/1465-1556/1520)  -9.5 W21-W28 N30-N33
      8 1561 1380 (N25W29 L=222)   603    146        2 (1507/1325-1617/1419) -20.6 W24-W33 N22-N27
      9 1771 1272 (N24W48 L=241)   790    174        0 (1719/1182-1824/1326) -28.6 W42-W53 N19-N26
     10  182 1073 (S15E63 L=130)  5945    311       14 ( 108/ 922- 277/1184)  58.0 E74-E52 S23-S07
     11 1121  976 (S08W06 L=199)   221     33        1 (1093/ 955-1146/ 998)  45.0 W04-W07 S08-S06
     12 1503  769 (S12W35 L=228)   632    218        0 (1448/ 691-1579/ 859)  71.6 W33-W37 S17-S05
     13  508  784 (S29E31 L=162)  1029    270        1 ( 405/ 745- 644/ 819)  -3.0 E39-E20 S30-S27
     14 1712  647 (S14W56 L=249)   892    185        2 (1679/ 575-1769/ 703)  76.0 W53-W60 S17-S09
     15  207  608 (S43E74 L=119)   788     48        0 ( 181/ 584- 229/ 631)   7.8 E88-E66 S44-S41
     16 1267  601 (S27W24 L=217)   203     60        0 (1231/ 589-1297/ 617)  33.7 W21-W25 S27-S25
     17  386  452 (S54E53 L=140)  1241    150        1 ( 328/ 397- 451/ 493)  12.2 E77-E40 S58-S50
     18 1542  445 (S29W52 L=245)  2640    357        2 (1384/ 397-1702/ 487)  28.4 W41-W65 S35-S22
  

    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