Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1683 1552 (N26W65 L= 55)   628     63        3 (1660/1519-1711/1592)  68.2 W62-W67 N23-N28
      2 1450 1496 (N28W39 L= 29)   628    185        0 (1418/1431-1503/1577) -90.0 W37-W42 N23-N33
      3 1618 1538 (N27W56 L= 46)   445     76        0 (1580/1517-1656/1555) -26.6 W53-W59 N24-N28
      4 1778 1482 (N20W72 L= 63)  2245    157        4 (1741/1415-1823/1545)  60.9 W68-W80 N14-N23
      5  693 1384 (N33E18 L=332)  1044    234        2 ( 647/1286- 747/1480) -63.4 E21-E14 N25-N40
      6 1057 1380 (N27W08 L=358)   279     98        0 (1022/1340-1100/1417) -56.3 W05-W09 N24-N30
      7 1674 1145 (N03W45 L= 35)   883    105        3 (1628/1105-1715/1199) -71.6 W41-W48 N00-N06
      8 1807 1069 (S05W56 L= 46)   480     85        1 (1788/1017-1819/1105) -90.0 W54-W57 S07-S02
      9  354  685 (S08E52 L=298)  1156    206        0 ( 308/ 603- 390/ 773)  42.3 E59-E46 S12-S02
     10 1097  653 (S18E01 L=349)   278    138        0 (1034/ 628-1152/ 687)   8.1 E04-W03 S18-S16
     11 1344  576 (S27W15 L=  5)   291     97        1 (1294/ 556-1391/ 601)   0.0 W11-W18 S28-S25
     12  883  425 (S31E20 L=330)   302     92        0 ( 871/ 379- 900/ 474)  76.0 E21-E17 S34-S26
  

    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