Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1365 1438 (N31W23 L=179)  1086    325        3 (1257/1356-1464/1552) -42.9 W15-W29 N25-N38
      2 1607 1517 (N37W47 L=203)  1439    178        3 (1528/1476-1677/1546)  -7.6 W37-W54 N34-N39
      3  385 1444 (N25E53 L=104)   388     58        1 ( 363/1420- 414/1464)  36.9 E55-E50 N23-N26
      4 1813 1328 (N24W64 L=220)  1426    212        1 (1727/1256-1874/1390) -31.0 W55-W72 N19-N28
      5 1748  855 (S05W53 L=209)   494     25        1 (1726/ 833-1773/ 882)  71.6 W50-W55 S06-S03
      6  298  735 (S20E54 L=103)   655     94        1 ( 257/ 708- 341/ 763)  -6.3 E57-E48 S21-S18
      7 1672  683 (S17W49 L=205)   365     88        0 (1654/ 644-1691/ 725)  80.5 W47-W50 S19-S13
      8 1785  613 (S21W66 L=222)  1310    182        2 (1751/ 543-1813/ 698) -61.4 W61-W72 S26-S15
      9  356  604 (S29E51 L=105)  2051    310        1 ( 248/ 493- 426/ 695) -38.2 E63-E46 S36-S22
     10  621  506 (S34E28 L=128)   457     79        0 ( 580/ 490- 664/ 527)  -8.1 E31-E24 S34-S32
  

    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