Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  847 1806 (N36E40 L=166)   631    210        1 ( 819/1744- 870/1856) -30.3 E46-E34 N32-N40
      2  534 1075 (S16E30 L=177)   638    165        1 ( 464/1008- 593/1124)  59.0 E33-E27 S21-S11
      3 1833 1011 (N18W53 L=260)   555     93        0 (1812/ 966-1850/1060) -38.7 W49-W55 N15-N20
      4 1001  966 (S10W00 L=207)   226     87        0 ( 963/ 936-1038/ 988)  14.0 E01-W03 S11-S08
      5  661  916 (S22E18 L=188)   253     96        0 ( 639/ 883- 698/ 968) -31.0 E20-E15 S23-S19
      6 1752  876 (N08W50 L=257)  1232    246        0 (1710/ 783-1815/ 959) -21.8 W43-W59 N04-N11
      7 1274  855 (S08W19 L=225)   341     94        0 (1251/ 807-1287/ 909) -59.0 W17-W20 S10-S05
      8  886  826 (S21E02 L=205)   194     55        0 ( 856/ 812- 921/ 842)  14.0 E04-W00 S21-S19
      9  446  685 (S42E30 L=176)   898    285        2 ( 345/ 664- 566/ 707)  19.3 E40-E20 S45-S37
     10 1268  526 (S26W31 L=238)   526    180        0 (1180/ 493-1367/ 552)  38.7 W26-W36 S30-S21
     11 1532  490 (S19W50 L=257)   652    246        1 (1457/ 468-1618/ 522)  14.9 W42-W57 S20-S16
  

    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