Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  563 1700 (N36E45 L=213)   353     46        1 ( 546/1677- 585/1728) -38.7 E48-E42 N34-N38
      2  438 1604 (N28E51 L=207)   686    127        3 ( 420/1533- 457/1661) -45.0 E56-E46 N23-N32
      3 1647 1558 (N36W45 L=303)   464     57        1 (1634/1520-1662/1592)  59.0 W42-W47 N33-N38
      4 1302 1497 (N27W14 L=272)   286     88        1 (1283/1455-1337/1529)  63.4 W13-W16 N23-N29
      5 1576 1207 (N11W34 L=291)   267     57        0 (1539/1186-1601/1221)   0.0 W31-W35 N10-N12
      6 1855 1074 (N07W59 L=317)   411     70        0 (1838/1035-1869/1113)  78.7 W57-W60 N04-N09
      7  702  976 (S10E19 L=239)   232     91        0 ( 687/ 930- 727/1019)  78.7 E20-E17 S12-S07
      8  227  832 (S21E58 L=199)   480     91        1 ( 186/ 800- 261/ 862) -29.7 E63-E54 S22-S18
      9 1532  757 (S16W36 L=293)   303    102        0 (1506/ 715-1556/ 810)  71.6 W34-W36 S18-S12
     10 1632  679 (S19W46 L=303)   266     80        0 (1598/ 650-1658/ 709)  53.1 W43-W47 S21-S17
     11 1707  644 (S20W54 L=312)   511     54        0 (1686/ 607-1729/ 676)  76.0 W51-W55 S22-S18
     12 1453  385 (S41W44 L=302)   772     99        2 (1417/ 349-1481/ 445) -66.8 W41-W49 S43-S36
     13  498  301 (S58E57 L=201)   652     50        1 ( 478/ 281- 530/ 325) -18.4 E65-E50 S59-S55
  

    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