Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1493 1507 (N27W34 L=227)   281     97        0 (1445/1474-1529/1530) -31.0 W30-W35 N25-N28
      2  323 1403 (N21E52 L=142)   401     83        1 ( 302/1378- 347/1447) -39.8 E55-E49 N19-N24
      3 1735 1414 (N22W53 L=247)   438     48        0 (1706/1387-1766/1440) -36.9 W50-W56 N20-N23
      4  712 1313 (N15E20 L=174)   724    226        2 ( 625/1276- 828/1342) -13.0 E25-E11 N12-N16
      5  270 1243 (N11E53 L=141)  1336    289        2 ( 199/1156- 331/1334)  45.0 E59-E47 N06-N17
      6  354 1175 (N07E45 L=149)   852    173        3 ( 310/1108- 400/1251)  53.1 E48-E41 N02-N11
      7  916  838 (S14E07 L=187)  1002    330        0 ( 844/ 701- 957/ 990)  67.6 E11-E04 S22-S05
      8 1468  805 (S16W29 L=222)   457    180        1 (1412/ 776-1520/ 860) -35.5 W24-W32 S17-S12
      9  908  472 (S38E09 L=185)   711    216        1 ( 803/ 452-1003/ 503)  10.6 E17-E01 S39-S36
     10 1694  444 (S39W64 L=257)   556     67        0 (1660/ 426-1729/ 462)  -3.8 W57-W72 S40-S37
     11  882  304 (S52E14 L=179)   514    114        0 ( 821/ 287- 939/ 319)  -5.2 E20-E08 S53-S50
  

    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