Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  551 1748 (N29E57 L=142)  1440    195        4 ( 510/1712- 602/1775)  33.7 E62-E52 N25-N31
      2  356 1558 (N12E63 L=136)   994    125        3 ( 326/1506- 392/1585)  78.7 E67-E59 N08-N13
      3 1857  865 (N10W63 L=262)   681    117        1 (1836/ 835-1882/ 903) -14.0 W59-W67 N08-N11
      4 1158  831 (S11W13 L=212)  1052    287        1 (1053/ 755-1261/ 886)  50.7 W08-W19 S18-S06
      5  741  829 (S22E12 L=187)   354    198        0 ( 694/ 786- 836/ 850)  36.9 E14-E06 S24-S18
      6 1474  797 (S05W32 L=232)   424     84        1 (1445/ 752-1492/ 846) -78.7 W29-W33 S07-S02
      7  283  656 (S45E49 L=151)   452     95        0 ( 253/ 638- 323/ 681)  -8.1 E51-E44 S45-S43
      8  605  621 (S38E17 L=182)  1077    324        0 ( 477/ 594- 757/ 643)  17.7 E28-E06 S41-S33
      9 1373  458 (S27W40 L=239)  2372    679        2 (1057/ 386-1597/ 511)  12.8 W17-W61 S34-S22
  

    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