(sys.argv[0]))
try:
opts, args = getopt.getopt(sys.argv[1:], "hl:", ["help", "lda="])
except getopt.GetoptError:
usage()
sys.exit(2)
if len(args) < 2:
usage()
sys.exit(2)
ldafn = None
for o, a in opts:
if o in ('-h', '--help'):
usage()
sys.exit()
if o in ('-l', '--lda'):
ldafn = a
mlltfn = args[0]
accumdirs = args[1:]
gauden = s3gaucnt.accumdirs_full(accumdirs)
m = MLLTModel(gauden)
mllt = m.train()
if ldafn != None:
# Compose this with the LDA transform if given
lda = s3lda.open(ldafn)[0]
ldadim = mllt.shape[1]
ldamllt = dot(mllt, lda[0:ldadim])
s3lda.open(mlltfn, 'w').writeall(ldamllt[newaxis, :])
else:
s3lda.open(mlltfn, 'w').writeall(mllt[newaxis, :])
sb += d * numpy.outer(diff, diff)
print "Sw:\n", sw
print "Sb:\n", sb
BinvA = numpy.dot(numpy.linalg.inv(sw), sb)
u, v = numpy.linalg.eig(BinvA)
top = list(u.argsort())
top.reverse()
u = u.take(top)
# Remember, the eigenvalues are in the columns, but Sphinx expects
# them to be in the rows.
v = v.T.take(top, 0)
print "u:\n", u
print "v:\n", v
return v
if __name__ == '__main__':
if len(sys.argv) < 3:
sys.stderr.write("Usage: %s OUTFILE ACCUMDIRS...\n" % (sys.argv[0]))
sys.exit(1)
ldafn = sys.argv[1]
accumdirs = sys.argv[2:]
gauden = s3gaucnt.accumdirs_full(accumdirs)
lda = makelda(gauden)
s3lda.open(ldafn, 'w').writeall(lda[numpy.newaxis, :])
sys.stderr.write("Usage: %s [-l INFILE] OUTFILE ACCUMDIRS...\n" % (sys.argv[0]))
try:
opts, args = getopt.getopt(sys.argv[1:], "hl:", ["help", "lda="])
except getopt.GetoptError:
usage()
sys.exit(2)
if len(args) < 2:
usage()
sys.exit(2)
ldafn = None
for o, a in opts:
if o in ('-h', '--help'):
usage()
sys.exit()
if o in ('-l', '--lda'):
ldafn = a
mlltfn = args[0]
accumdirs = args[1:]
gauden = s3gaucnt.accumdirs_full(accumdirs)
m = MLLTModel(gauden)
mllt = m.train()
if ldafn != None:
# Compose this with the LDA transform if given
lda = s3lda.open(ldafn)[0]
ldadim = mllt.shape[1]
ldamllt = dot(mllt, lda[0:ldadim])
s3lda.open(mlltfn, 'w').writeall(ldamllt[newaxis,:])
else:
s3lda.open(mlltfn, 'w').writeall(mllt[newaxis,:])
diff = m / d - globalmean
sb += d * numpy.outer(diff, diff)
print "Sw:\n", sw
print "Sb:\n", sb
BinvA = numpy.dot(numpy.linalg.inv(sw), sb)
u, v = numpy.linalg.eig(BinvA)
top = list(u.argsort())
top.reverse()
u = u.take(top)
# Remember, the eigenvalues are in the columns, but Sphinx expects
# them to be in the rows.
v = v.T.take(top, 0)
print "u:\n", u
print "v:\n", v
return v
if __name__ == '__main__':
if len(sys.argv) < 3:
sys.stderr.write("Usage: %s OUTFILE ACCUMDIRS...\n" % (sys.argv[0]))
sys.exit(1)
ldafn = sys.argv[1]
accumdirs = sys.argv[2:]
gauden = s3gaucnt.accumdirs_full(accumdirs)
lda = makelda(gauden)
s3lda.open(ldafn, 'w').writeall(lda[numpy.newaxis,:])
