Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1378 1441 (N31W24 L=179)  1049    314        3 (1258/1355-1478/1555) -40.9 W15-W30 N25-N38
      2 1616 1517 (N37W48 L=204)  1365    179        3 (1539/1482-1688/1545)   3.2 W38-W56 N34-N38
      3  388 1438 (N25E52 L=103)   489     65        2 ( 355/1414- 419/1460)  18.4 E55-E48 N23-N26
      4 1852 1305 (N22W69 L=225)   903    120        0 (1817/1258-1878/1353) -56.3 W65-W73 N19-N25
      5  354  929 (S07E46 L=110)   313     64        2 ( 331/ 902- 379/ 955)  45.0 E47-E43 S08-S05
      6 1752  856 (S05W53 L=208)   540     27        1 (1728/ 831-1776/ 882)  90.0 W50-W55 S06-S03
      7  310  733 (S20E52 L=103)   555    100        0 ( 268/ 710- 351/ 762)  -6.3 E56-E47 S21-S18
      8 1680  680 (S17W50 L=205)   365     59        0 (1663/ 647-1696/ 718)  90.0 W48-W50 S18-S14
      9  286  675 (S24E57 L= 98)   511     53        2 ( 262/ 646- 317/ 698)  12.5 E61-E52 S26-S22
     10  391  581 (S30E48 L=107)  1391    235        1 ( 321/ 490- 436/ 684) -69.0 E52-E45 S36-S23
     11 1791  609 (S22W67 L=222)  1230    160        2 (1756/ 544-1820/ 683)  90.0 W62-W74 S26-S16
     12  630  508 (S33E27 L=128)   537     92        1 ( 584/ 492- 681/ 529)   0.0 E31-E23 S34-S31
  

    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