Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  695 1739 (N60E01 L=122)   318     50        1 ( 680/1710- 713/1765)  45.0 E03-W02 N57-N62
      2 1136 1573 (N33W24 L=147)   630    122        2 (1072/1537-1187/1599) -50.2 W20-W26 N29-N35
      3  678 1563 (N46E08 L=115)   277     76        2 ( 642/1540- 716/1582)   0.0 E11-E04 N44-N47
      4 1357 1505 (N22W34 L=158)   400    136        1 (1309/1461-1412/1539)   5.7 W29-W39 N21-N23
      5  617 1178 (N23E21 L=103)   369    119        1 ( 594/1122- 645/1237)  82.9 E21-E19 N19-N27
      6  211 1127 (N29E57 L= 66)  1049     86        2 ( 169/1088- 254/1169)  20.6 E63-E51 N26-N31
      7 1840  957 (S22W54 L=177)   815     98        1 (1820/ 912-1864/1000) -66.8 W51-W57 S25-S18
      8 1681  646 (S36W35 L=159)  2808    633        2 (1468/ 508-1896/ 858)   8.5 W15-W62 S40-S30
      9  756  669 (S09E24 L= 99)   487     94        1 ( 728/ 622- 794/ 712)  26.6 E27-E21 S11-S06
     10  812  393 (S27E32 L= 91)   829    144        2 ( 759/ 372- 892/ 424) -40.6 E35-E27 S30-S23
     11  618  342 (S26E49 L= 74)   580    129        2 ( 575/ 302- 672/ 367)  10.3 E55-E43 S27-S24
  

    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