Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1329 1831 (N67W19 L= 70)  3511    317        3 (1193/1782-1490/1871)   2.2 E02-W51 N62-N71
      2 1477 1615 (N47W29 L= 80)  2051    285        3 (1333/1559-1576/1644)  -2.3 W14-W40 N43-N50
      3  414 1583 (N26E59 L=351)  1301    138        3 ( 366/1540- 465/1623)  26.6 E65-E54 N23-N29
      4  894 1539 (N31E17 L= 33)   489    135        1 ( 867/1484- 930/1603) -48.8 E20-E13 N27-N35
      5  353 1505 (N20E61 L=350)  1022     83        1 ( 317/1479- 401/1534)  14.0 E65-E54 N18-N22
      6  761 1402 (N20E23 L= 28)   557    169        1 ( 713/1360- 830/1465) -26.6 E27-E17 N17-N23
      7  797 1279 (N13E18 L= 33)   295    111        0 ( 755/1256- 842/1298)  21.8 E20-E15 N11-N14
      8 1877 1193 (N22W68 L=118)  1509    174        1 (1819/1127-1919/1245) -24.8 W59-W75 N18-N25
      9 1659  755 (S07W47 L= 98)   978    244        1 (1589/ 694-1718/ 827)  52.1 W42-W51 S11-S02
     10  523  704 (S27E30 L= 20)   925    268        0 ( 463/ 662- 576/ 786)  52.1 E35-E25 S29-S20
     11  987  723 (S18W02 L= 53)   288    103        0 ( 940/ 693-1034/ 749)  -9.5 E01-W05 S20-S17
     12  675  637 (S29E18 L= 32)   453    161        0 ( 610/ 595- 738/ 691)  41.6 E23-E13 S32-S24
     13 1710  627 (S15W57 L=108)  1007    138        4 (1652/ 589-1754/ 672)  32.0 W52-W61 S17-S11
     14  940  563 (S30W01 L= 52)  1728    583        2 ( 689/ 515-1198/ 608)  17.9 E16-W18 S36-S23
     15 1509  471 (S28W46 L= 97)   333    108        0 (1468/ 446-1552/ 511)  51.3 W43-W47 S29-S24
     16  931  160 (S68W16 L= 67)   779     56        1 ( 912/ 127- 951/ 193) -38.2 W10-W24 S74-S63
  

    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