Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1087 1911 (N58W59 L=315)   733     62        2 (1050/1896-1116/1929)  -5.7 W54-W66 N55-N60
      2 1011 1464 (N31W12 L=268)   201     52        0 ( 981/1451-1045/1476) -26.6 W10-W14 N29-N31
      3 1707  844 (S24W39 L=294)   595    149        0 (1644/ 798-1769/ 888)   0.0 W32-W45 S25-S23
      4 1499  664 (S28W19 L=275)   302     98        0 (1484/ 622-1512/ 716)  71.6 W17-W21 S31-S25
      5 1432  626 (S29W14 L=270)   514    127        0 (1392/ 570-1465/ 677)  39.8 W10-W16 S31-S26
      6 1566  626 (S33W24 L=279)   477     72        0 (1521/ 611-1610/ 646) -15.9 W20-W27 S34-S31
      7 1242  574 (S26W00 L=256)   506    151        0 (1216/ 517-1273/ 641)  49.4 E02-W04 S29-S22
      8  595  375 (S22E49 L=207)   474     96        1 ( 566/ 346- 625/ 414)  53.1 E51-E46 S23-S19
  

    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