Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  571 1522 (N47E26 L= 96)   266     44        0 ( 554/1498- 589/1542)  26.6 E27-E23 N45-N49
      2  397 1518 (N50E47 L= 74)   361     53        1 ( 377/1489- 417/1540)  71.6 E50-E44 N48-N52
      3  681 1491 (N43E15 L=106)   388     94        1 ( 636/1463- 720/1526) -45.0 E18-E11 N40-N46
      4 1728 1410 (N12W60 L=181)   681    132        2 (1691/1355-1767/1447)  29.7 W55-W65 N09-N14
      5  199 1134 (N26E63 L= 58)  1498    271        1 ( 133/1033- 252/1207) -32.9 E75-E55 N19-N30
      6 1869 1150 (S06W65 L=186)   859     93        3 (1827/1116-1904/1184) -31.0 W60-W70 S08-S02
      7  605  783 (N00E30 L= 91)   705    149        1 ( 572/ 710- 659/ 840) -74.1 E32-E28 S04-N03
      8 1586  653 (S28W30 L=152)   578     98        1 (1561/ 607-1609/ 698)  71.6 W28-W32 S31-S25
      9  407  497 (S15E56 L= 65)   569    147        1 ( 367/ 456- 436/ 562)  23.2 E64-E50 S17-S11
     10 1010  331 (S38E18 L=103)   336    105        0 ( 979/ 310-1057/ 369) -51.3 E20-E15 S40-S34
  

    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