Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1620 1541 (N21W58 L=187)   608    125        0 (1569/1521-1684/1560) -29.1 W53-W62 N17-N22
      2 1245 1485 (N28W27 L=156)   363     96        1 (1213/1457-1274/1518)  18.4 W23-W29 N26-N29
      3 1622 1462 (N16W53 L=182)   353     51        0 (1598/1435-1646/1482) -56.3 W51-W54 N14-N17
      4  189 1050 (N24E64 L= 66)  1741    202        7 ( 144/ 964- 225/1117) -45.0 E71-E58 N18-N27
      5 1848  942 (S20W57 L=186)   727     98        2 (1813/ 916-1882/ 973)   0.0 W52-W61 S21-S17
      6 1620  857 (S17W34 L=163)   682    166        2 (1573/ 790-1680/ 914)  26.6 W29-W39 S20-S14
      7  835  511 (S19E24 L=105)  1998    414        2 ( 692/ 365- 945/ 650)  26.6 E39-E13 S26-S12
      8  591  363 (S25E49 L= 80)  1846    353        2 ( 522/ 297- 696/ 495)  82.4 E54-E41 S31-S15
      9  930  458 (S25E20 L=109)   519    171        1 ( 875/ 412- 991/ 493) -39.8 E24-E17 S29-S21
     10 1070  435 (S30E12 L=117)   305     97        0 (1026/ 414-1113/ 465)   8.1 E15-E08 S30-S28
  

    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