Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  915 1478 (N37W04 L=257)  1013    211        2 ( 818/1403- 998/1536) -50.7 E02-W08 N30-N42
      2 1674 1396 (N12W53 L=306)   334     75        0 (1633/1372-1705/1416) -38.7 W50-W55 N10-N14
      3 1013 1310 (N24W06 L=259)   218     22        0 ( 993/1287-1033/1331)  71.6 W04-W07 N22-N25
      4 1378 1258 (N12W27 L=280)   240     63        1 (1346/1232-1400/1287) -90.0 W24-W27 N10-N14
      5  383 1138 (N26E41 L=213)  2816    588        6 ( 216/1029- 536/1252) -29.1 E56-E29 N17-N34
      6 1724  887 (S18W43 L=296)  2062    345        2 (1636/ 748-1771/1072)  46.8 W33-W49 S24-S08
      7  712  987 (N11E19 L=234)   303     62        2 ( 677/ 961- 741/1016)  14.0 E21-E17 N09-N12
      8  840  937 (N06E13 L=241)   463     79        2 ( 797/ 906- 881/ 975) -36.9 E15-E10 N03-N08
      9  822  743 (S05E18 L=236)   534    157        0 ( 742/ 710- 897/ 774)   0.0 E23-E12 S06-S04
     10 1213  600 (S22W02 L=255)   264     81        0 (1187/ 558-1230/ 631)  78.7 W00-W03 S25-S20
     11 1558  504 (S39W26 L=279)   427     76        0 (1518/ 473-1597/ 530)   0.0 W22-W29 S40-S37
     12  565  393 (S23E47 L=206)   374     77        1 ( 531/ 357- 608/ 421) -71.6 E49-E45 S26-S20
  

    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