Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  870 1765 (N60W11 L=272)   351     82        0 ( 824/1738- 900/1789)  10.3 W04-W15 N59-N62
      2 1444 1752 (N40W62 L=322)   519     95        0 (1404/1719-1472/1785)   6.7 W53-W70 N39-N41
      3 1262 1644 (N39W35 L=296)   977    329        2 (1150/1600-1397/1670) -25.2 W24-W44 N33-N43
      4  498 1637 (N59E39 L=222)   513     67        0 ( 462/1616- 527/1662)  20.6 E44-E34 N57-N60
      5 1581  744 (S22W30 L=291)   304     84        1 (1553/ 708-1606/ 784)  53.1 W28-W32 S24-S20
      6  333  544 (S11E62 L=199)   724    110        0 ( 283/ 521- 386/ 571)   7.1 E70-E54 S12-S09
      7 1434  490 (S36W17 L=277)   371     94        1 (1401/ 448-1464/ 539)  45.0 W13-W19 S38-S33
      8  435  413 (S21E60 L=201)  1127    101        2 ( 385/ 378- 478/ 447) -45.0 E63-E55 S24-S18
  

    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