Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  993 1747 (N43E08 L= 63)   483     85        0 ( 960/1714-1032/1770)  38.7 E10-E05 N40-N45
      2 1368 1600 (N35W21 L= 92)  3435    770       11 (1098/1465-1673/1713) -12.8 W00-W44 N27-N41
      3  739 1625 (N32E24 L= 46)   692    177        5 ( 697/1563- 777/1677)  -3.3 E28-E20 N27-N36
      4 1247 1279 (N12W12 L= 83)   438    174        1 (1225/1187-1277/1363) -70.0 W10-W14 N06-N17
      5  610 1306 (N10E28 L= 43)   360    131        1 ( 558/1259- 652/1329) -29.7 E31-E24 N07-N11
      6 1823 1219 (N14W55 L=126)   475    129        2 (1789/1190-1866/1252) -20.6 W51-W60 N12-N15
      7  290 1196 (N03E51 L= 20)   279     76        0 ( 263/1165- 315/1234)  53.1 E53-E48 N01-N05
      8  679  707 (S26E21 L= 50)   278     92        0 ( 646/ 673- 715/ 746)  45.0 E23-E18 S27-S22
      9  299  571 (S35E55 L= 16)  1005    257        3 ( 243/ 535- 393/ 622) -13.0 E60-E46 S37-S31
     10 1635  510 (S31W50 L=121)   821    260        2 (1541/ 467-1742/ 563)  23.6 W42-W58 S34-S26
  

    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