Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1446 1740 (N32W55 L=101)   807    203        2 (1390/1702-1503/1765)  -4.8 W49-W61 N30-N34
      2 1029 1249 (N18W06 L= 52)   695    215        3 ( 933/1225-1155/1275) -30.3 W00-W13 N14-N21
      3 1524 1221 (N03W34 L= 80)   410     71        2 (1487/1189-1554/1247) -36.9 W32-W36 N00-N04
      4 1894 1150 (S14W64 L=110)  2204    345        4 (1841/1103-1951/1192)  -3.4 W56-W73 S17-S10
      5 1159  848 (S07W03 L= 49)  1000    381        1 (1037/ 782-1321/ 958) -32.7 E02-W12 S11-N01
      6 1678  770 (S28W34 L= 80)  1813    371        1 (1562/ 681-1740/ 934)  40.2 W24-W41 S32-S19
      7  807  898 (N05E15 L= 31)   212     69        0 ( 769/ 883- 847/ 914) -36.9 E17-E13 N03-N06
      8 1266  677 (S20W04 L= 50)   685    229        2 (1240/ 585-1312/ 805)  47.7 W00-W10 S25-S13
      9  289  667 (N03E60 L=346)  1275    304        1 ( 210/ 625- 356/ 717)   0.0 E68-E50 N00-N06
  

    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