Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1402 1722 (N44W51 L= 13)   493     63        2 (1370/1699-1436/1739)   7.1 W45-W54 N42-N45
      2  639 1443 (N38E22 L=299)   533    166        0 ( 568/1417- 713/1468) -14.0 E28-E16 N35-N39
      3 1573 1321 (N15W42 L=  4)   748    152        1 (1543/1263-1616/1389)  83.7 W39-W43 N10-N19
      4  866 1194 (N19E07 L=314)   440    163        2 ( 845/1132- 882/1277)  71.6 E09-E05 N14-N23
      5  435 1170 (N22E39 L=283)   381     87        0 ( 397/1140- 474/1201) -39.8 E42-E35 N19-N24
      6 1712  707 (S25W45 L=  7)   909    193        1 (1656/ 623-1764/ 782)  45.0 W40-W50 S30-S21
      7  514  400 (S28E51 L=271)  2310    271        4 ( 470/ 284- 566/ 509)  27.9 E74-E40 S38-S20
      8  987  327 (S40E16 L=306)   283     65        0 ( 952/ 315-1022/ 344)   9.5 E18-E12 S40-S38
  

    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