Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1243 1831 (N57E16 L=354)  1459    143        2 (1170/1789-1320/1863)   0.0 E25-E06 N54-N60
      2  607 1709 (N26E52 L=318)  2701    376        6 ( 567/1559- 663/1829) -25.9 E73-E38 N16-N34
      3 1659 1629 (N57W39 L= 49)   839     78        1 (1637/1588-1678/1666)  38.2 W34-W50 N52-N63
      4  505  847 (S28E28 L=342)   442    114        0 ( 491/ 780- 523/ 899) -66.8 E30-E26 S31-S24
      5  649  737 (S31E16 L=354)   237     99        0 ( 596/ 723- 695/ 756)  29.7 E19-E12 S32-S28
      6  788  658 (S32E04 L=  6)   320     78        0 ( 752/ 637- 822/ 681)  45.0 E06-E02 S33-S29
      7 1340  600 (S19W32 L= 42)   296    103        0 (1283/ 586-1388/ 615)  33.7 W28-W34 S21-S17
      8  632  327 (S63E08 L=  2)   531     68        3 ( 594/ 302- 671/ 348)  26.6 E13-E03 S66-S59
      9 1399  291 (S35W58 L= 68)  1323    137        3 (1336/ 269-1460/ 324)  12.1 W49-W64 S37-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