Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  612 1707 (N61E12 L=126)   444     79        0 ( 577/1679- 652/1742)  10.3 E17-E05 N59-N63
      2  528 1498 (N47E23 L=115)   301     52        0 ( 510/1467- 547/1521)  76.0 E25-E21 N44-N48
      3 1177 1405 (N22W20 L=158)   297    112        0 (1135/1378-1226/1428)   0.0 W16-W23 N20-N23
      4 1802  974 (S20W51 L=189)   707    127        2 (1769/ 924-1857/1008)  10.3 W46-W57 S21-S17
      5 1723  732 (S32W40 L=178)   795    197        1 (1647/ 677-1796/ 797)   7.1 W32-W48 S33-S30
      6  909  668 (S12E16 L=122)   194     48        1 ( 891/ 638- 923/ 699)  63.4 E17-E14 S14-S10
      7  564  591 (S09E40 L= 98)   655    127        0 ( 534/ 527- 598/ 648)  45.0 E43-E36 S12-S05
      8 1370  466 (S39W07 L=145)   261     49        2 (1342/ 451-1398/ 484) -14.0 W04-W08 S39-S37
      9  632  340 (S26E49 L= 89)  1955    191        3 ( 583/ 282- 719/ 417) -71.6 E54-E43 S30-S20
     10 1021  363 (S35E21 L=117)   428    102        0 ( 981/ 331-1057/ 405)  14.0 E24-E16 S36-S32
     11  521  322 (S25E61 L= 77)   535     66        0 ( 489/ 309- 558/ 335)  -5.7 E65-E55 S25-S23
  

    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