Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  863 1762 (N60W10 L=271)   362     91        0 ( 816/1735- 895/1788)   9.5 W03-W15 N58-N62
      2 1438 1748 (N40W60 L=321)   518     90        0 (1399/1715-1466/1784)   7.1 W52-W68 N39-N41
      3 1164 1633 (N41W27 L=289)   477    107        2 (1138/1583-1195/1679)  40.6 W23-W32 N37-N43
      4 1300 1642 (N38W38 L=299)   641    188        0 (1212/1618-1378/1663) -31.0 W31-W43 N34-N41
      5  500 1637 (N59E39 L=223)   422     53        0 ( 471/1617- 526/1661)  15.9 E42-E34 N57-N60
      6 1572  740 (S22W30 L=291)   263     74        1 (1544/ 714-1595/ 775)  26.6 W27-W31 S23-S20
      7  334  547 (S11E61 L=200)   754    119        0 ( 276/ 518- 381/ 573)   6.7 E71-E54 S13-S09
      8 1429  489 (S36W16 L=278)   346     98        1 (1394/ 445-1464/ 532)  45.0 W13-W19 S38-S33
      9  433  411 (S22E60 L=201)   998    118        4 ( 381/ 376- 475/ 443) -50.2 E64-E56 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