Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1742 (N56W08 L=218)   917     91        1 (1347/1716-1433/1780)  29.1 W03-W12 N51-N61
      2 1148 1640 (N38E09 L=201)   688    179        2 (1103/1569-1191/1712) -36.0 E15-E03 N34-N42
      3 1664 1495 (N44W38 L=248)   922    115        1 (1617/1452-1705/1529) -11.3 W32-W42 N41-N46
      4 1492  778 (S05W34 L=244)   244     27        1 (1475/ 756-1509/ 801) -33.7 W32-W35 S06-S04
      5  324  715 (S38E46 L=164)   325     63        0 ( 301/ 691- 357/ 733)  29.7 E49-E42 S40-S36
      6  377  553 (S49E43 L=167)   796     80        2 ( 344/ 516- 411/ 591)  33.7 E48-E37 S52-S45
      7  417  491 (S53E40 L=171)   423     68        0 ( 390/ 465- 454/ 510)  20.0 E44-E33 S55-S50
      8 1480  386 (S28W54 L=264)  2693    441        3 (1259/ 320-1636/ 455)  35.3 W39-W64 S37-S19
      9  507  404 (S58E29 L=181)   408     33        0 ( 482/ 383- 529/ 424)   0.0 E33-E25 S60-S55
  

    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