Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1749 (N52E23 L= 69)   880    111        1 ( 616/1711- 735/1784)  29.1 E28-E17 N48-N54
      2  351 1435 (N31E45 L= 47)  1456    314        4 ( 254/1371- 414/1471) -27.9 E56-E39 N25-N34
      3  823 1333 (N19E10 L= 83)   583    162        0 ( 767/1292- 881/1369) -23.2 E13-E06 N16-N21
      4 1527 1270 (N09W34 L=126)   450     95        0 (1474/1229-1562/1321) -56.3 W31-W36 N06-N12
      5  614  814 (S10E27 L= 65)   245     60        0 ( 582/ 798- 647/ 832) -14.0 E29-E25 S11-S08
      6  501  497 (S28E43 L= 49)  2927    641        5 ( 380/ 345- 584/ 700)  35.9 E61-E32 S37-S16
      7 1730  521 (S39W52 L=145)   573     96        1 (1694/ 496-1754/ 540)   9.5 W48-W54 S41-S37
      8 1639  471 (S42W45 L=137)   406     94        0 (1607/ 458-1673/ 487) -15.9 W41-W48 S43-S40
  

    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