Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  558 1702 (N36E45 L=213)   385     50        2 ( 538/1677- 581/1731) -33.7 E49-E42 N34-N38
      2  434 1599 (N28E51 L=207)   794    138        4 ( 420/1524- 457/1663) -45.0 E56-E46 N22-N32
      3 1645 1555 (N36W44 L=303)   456     42        2 (1629/1523-1662/1588)  68.2 W42-W46 N33-N38
      4 1319 1505 (N28W16 L=274)   479    143        1 (1277/1456-1384/1532)  40.6 W12-W20 N23-N29
      5 1572 1208 (N11W33 L=291)   304     54        0 (1536/1186-1598/1226)   0.0 W30-W35 N09-N12
      6 1852 1075 (N07W59 L=317)   452     69        1 (1836/1037-1867/1115)  78.7 W56-W60 N04-N09
      7  697  973 (S10E20 L=239)   215    101        0 ( 679/ 930- 726/1021)  68.2 E20-E17 S12-S07
      8  204  847 (S20E61 L=197)   433    120        1 ( 158/ 811- 242/ 882) -21.8 E66-E56 S21-S17
      9 1522  745 (S17W35 L=293)   322     93        0 (1495/ 709-1547/ 792)  80.5 W33-W36 S19-S13
     10 1630  680 (S19W45 L=304)   293     77        0 (1598/ 651-1657/ 711)  53.1 W43-W47 S21-S17
     11 1706  647 (S20W54 L=312)   483     65        0 (1689/ 611-1729/ 676) -76.0 W51-W55 S22-S18
     12 1451  387 (S41W44 L=302)   788     97        2 (1417/ 349-1477/ 445) -60.3 W41-W48 S43-S36
     13  497  302 (S58E57 L=201)   734     77        1 ( 470/ 282- 538/ 327)   0.0 E66-E47 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