Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  557 1731 (N36E59 L=247)   736    110        0 ( 534/1690- 587/1783) -18.4 E68-E50 N32-N39
      2  721 1606 (N30E34 L=272)   298     73        1 ( 688/1575- 750/1629) -15.9 E37-E30 N28-N31
      3 1208 1494 (N32W03 L=309)   809    264        1 (1126/1408-1295/1563) -33.7 E03-W10 N27-N35
      4  895 1414 (N21E15 L=290)   341     93        0 ( 850/1378- 933/1453) -23.2 E19-E12 N18-N22
      5 1654  684 (S10W50 L=355)   645    147        1 (1596/ 656-1724/ 713)  26.6 W45-W55 S12-S07
      6  764  617 (S30E10 L=296)   775    266        2 ( 651/ 581- 891/ 649)  19.4 E18-E01 S34-S26
      7 1232  463 (S32W27 L=332)  1547    262        3 (1184/ 371-1300/ 567) -30.5 W19-W36 S37-S25
      8 1434  380 (S33W47 L=353)   988    117        1 (1373/ 357-1496/ 411)  16.7 W42-W53 S34-S31
  

    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