Python MarSystemManager Examples

Example #1

File: pvoc.py Project: Amos-zq/marsyas

def classic(infile):
	print "Classic phasevocoding for " + infile 
	# Build the processing network
	mng = MarSystemManager()
	a = ("Series", "pvseries")
	pvseries = mng.create("Series", "pvseries")
	pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
	pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))
	
	filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
	inSamples = pvseries.getControl("mrs_natural/inSamples")
	winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
	fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
	analysisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Decimation")
	synthesisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Interpolation")
	convertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/convertMode")
	unconvertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/unconvertMode")
	outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
	notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
	phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
	magHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes"); 
	peakHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks");
	sinusoids = pvseries.getControl("PhaseVocoder/pvoc/PvConvert/sinusoids");
	
	pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos");
	
	
	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = name + "_classic" + ext
	
    #classic

	filename.setValue_string(infile);
	inSamples.setValue_natural(512);
	analysisHop.setValue_natural(512);
	synthesisHop.setValue_natural(1024);
	winSize.setValue_natural(4096);
	fftSize.setValue_natural(4096);
	sinusoids.setValue_natural(40);
	convertMode.setValue_string("sorted");
	unconvertMode.setValue_string("classic");
	outfname.setValue_string(outsfname);
	
	i = 0 
	while notempty.to_bool():
		if (i == 0):
			phaselock.setValue_bool(True)
			
		i = i + 1
		pvseries.tick()

		
	print "Classic phasevocoding done. Output in "  + outsfname 

Example #2

File: pvoc.py Project: Amos-zq/marsyas

def oscbank(infile):
	print "Osc Bank phasevocoding for " + infile 
	# Build the processing network
	mng = MarSystemManager()
	a = ("Series", "pvseries")
	pvseries = mng.create("Series", "pvseries")
	pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseries.addMarSystem(mng.create("PhaseVocoderOscBank", "pvoc"))
	pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))
	
	filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
	inSamples = pvseries.getControl("mrs_natural/inSamples")
	winSize = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_natural/winSize")
	fftSize = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_natural/FFTSize")
	analysisHop = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_natural/Decimation")
	synthesisHop = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_natural/Interpolation")
	convertMode = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_string/convertMode")
	outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
	notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
	phaselock = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_bool/phaselock")
	pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos");
	sinusoids = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_natural/Sinusoids");
	pitchShift = pvseries.getControl("PhaseVocoderOscBank/pvoc/mrs_real/PitchShift");
	
	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = name + "_oscbank" + ext

	winSize.setValue_natural(2048);
	fftSize.setValue_natural(2048);
	synthesisHop.setValue_natural(512);
	analysisHop.setValue_natural(512);
	sinusoids.setValue_natural(120);
	pitchShift.setValue_real(0.66);
	
	convertMode.setValue_string("sorted");
	inSamples.setValue_natural(512);
	filename.setValue_string(infile);
	


	outfname.setValue_string(outsfname);
	


	i = 0 
	while notempty.to_bool():
		i = i + 1
		pvseries.tick()

Example #3

def classic(infile):
    print "Classic phasevocoding for " + infile
    # Build the processing network
    mng = MarSystemManager()
    a = ("Series", "pvseries")
    pvseries = mng.create("Series", "pvseries")
    pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
    pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
    inSamples = pvseries.getControl("mrs_natural/inSamples")
    winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
    fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
    analysisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Decimation")
    synthesisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Interpolation")
    convertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/convertMode")
    unconvertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/unconvertMode")
    outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
    notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
    phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
    magHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes")
    peakHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks")
    sinusoids = pvseries.getControl("PhaseVocoder/pvoc/PvConvert/sinusoids")

    pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos")

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = name + "_classic" + ext

    #classic

    filename.setValue_string(infile)
    inSamples.setValue_natural(512)
    analysisHop.setValue_natural(512)
    synthesisHop.setValue_natural(1024)
    winSize.setValue_natural(4096)
    fftSize.setValue_natural(4096)
    sinusoids.setValue_natural(40)
    convertMode.setValue_string("sorted")
    unconvertMode.setValue_string("classic")
    outfname.setValue_string(outsfname)

    i = 0
    while notempty.to_bool():
        if (i == 0):
            phaselock.setValue_bool(True)

        i = i + 1
        pvseries.tick()

    print "Classic phasevocoding done. Output in " + outsfname

Example #4

def oscbank(infile):
    print "Osc Bank phasevocoding for " + infile
    # Build the processing network
    mng = MarSystemManager()
    a = ("Series", "pvseries")
    pvseries = mng.create("Series", "pvseries")
    pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseries.addMarSystem(mng.create("PhaseVocoderOscBank", "pvoc"))
    pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
    inSamples = pvseries.getControl("mrs_natural/inSamples")
    winSize = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_natural/winSize")
    fftSize = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_natural/FFTSize")
    analysisHop = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_natural/Decimation")
    synthesisHop = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_natural/Interpolation")
    convertMode = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_string/convertMode")
    outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
    notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
    phaselock = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_bool/phaselock")
    pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos")
    sinusoids = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_natural/Sinusoids")
    pitchShift = pvseries.getControl(
        "PhaseVocoderOscBank/pvoc/mrs_real/PitchShift")

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = name + "_oscbank" + ext

    winSize.setValue_natural(2048)
    fftSize.setValue_natural(2048)
    synthesisHop.setValue_natural(512)
    analysisHop.setValue_natural(512)
    sinusoids.setValue_natural(120)
    pitchShift.setValue_real(0.66)

    convertMode.setValue_string("sorted")
    inSamples.setValue_natural(512)
    filename.setValue_string(infile)

    outfname.setValue_string(outsfname)

    i = 0
    while notempty.to_bool():
        i = i + 1
        pvseries.tick()

Example #5

File: testPylab.py Project: xcode2010/Extracting-Features-from-audio

# EXAMPLE OF USING Marsyas together with PyLab
# for MATLAB-like plotting
# run using ipython -pylab testPylab.py foo.wav

#!/usr/bin/env python

from pylab import *
import sys
from marsyas import MarSystemManager, MarControlPtr

msm = MarSystemManager()

file = msm.create("SoundFileSource", "file")
sl = msm.create("StretchLinear", "sl")
gain = msm.create("Gain", "g")
pipe = msm.create("Series", "pipe")

pipe.addMarSystem(file)
pipe.addMarSystem(sl)
pipe.addMarSystem(gain)

filename = pipe.getControl("SoundFileSource/file/mrs_string/filename")
notempty = pipe.getControl("SoundFileSource/file/mrs_bool/hasData")
outData1 = pipe.getControl("SoundFileSource/file/mrs_realvec/processedData")
outData2 = pipe.getControl("StretchLinear/sl/mrs_realvec/processedData")
alpha = pipe.getControl("StretchLinear/sl/mrs_real/stretch")

alpha.setValue_real(2.0)

try:
    filename.setValue_string(sys.argv[1])

Example #6

File: wekatest.py Project: Amos-zq/marsyas

#!/usr/bin/env python

from pylab import *
from marsyas import MarSystemManager,MarControlPtr

msm = MarSystemManager()

weka= msm.create("WekaSource", "weka")
classifier= msm.create("GaussianClassifier", "classifier")
summary= msm.create("Summary", "summary")
pipe = msm.create("Series","pipe")

pipe.addMarSystem(weka)
pipe.addMarSystem(classifier)
pipe.addMarSystem(summary)

filename = pipe.getControl("WekaSource/weka/mrs_string/filename")
validationmode = pipe.getControl("WekaSource/weka/mrs_string/validationMode")


validationmode.setValue_string("PercentageSplit,50%")
filename.setValue_string("ms.arff")

insamples= pipe.getControl("mrs_natural/inSamples")
insamples.setValue_natural(1)

wekanclasses= pipe.getControl("WekaSource/weka/mrs_natural/nClasses")
summarynclasses= pipe.getControl("Summary/summary/mrs_natural/nClasses")
classifiernclasses= pipe.getControl("GaussianClassifier/classifier/mrs_natural/nClasses")
summarynclasses.setValue_natural(wekanclasses.to_natural());
classifiernclasses.setValue_natural(wekanclasses.to_natural())

Example #7

File: aim.py Project: puentev/nupic.audio

#!/usr/bin/env python


# Plotting auditory image model stuff

import pickle
from pylab import *
import sys
from matplotlib import pyplot
from marsyas import MarSystemManager,MarControlPtr

# create a global MarSystemManager
msm = MarSystemManager()







# helper function to create a list of MarSystems from a list of names
def create(msys_list_names):
	return map(msm.create,msys_list_names)

# create a Series MarSystem composed of a list of MarSystems
def series(name, msys_list):
	net = msm.create("Series/"+name)
	map(net.addMarSystem, msys_list)
	return net

# create a Fanout MarSystem composed of a list of MarSystems

Example #8

File: testPylab.py Project: Amos-zq/marsyas

# EXAMPLE OF USING Marsyas together with PyLab
# for MATLAB-like plotting
# run using ipython -pylab testPylab.py foo.wav 

#!/usr/bin/env python

from pylab import *
import sys
from marsyas import MarSystemManager,MarControlPtr

msm = MarSystemManager()

file = msm.create("SoundFileSource","file")
sl = msm.create("StretchLinear", "sl")
gain = msm.create("Gain", "g") 
pipe = msm.create("Series","pipe")

pipe.addMarSystem(file)
pipe.addMarSystem(sl)
pipe.addMarSystem(gain)

filename = pipe.getControl("SoundFileSource/file/mrs_string/filename")
notempty = pipe.getControl("SoundFileSource/file/mrs_bool/hasData")
outData1  = pipe.getControl("SoundFileSource/file/mrs_realvec/processedData")
outData2 = pipe.getControl("StretchLinear/sl/mrs_realvec/processedData")
alpha = pipe.getControl("StretchLinear/sl/mrs_real/stretch")

alpha.setValue_real(2.0)

try:
	filename.setValue_string(sys.argv[1])

Example #9

def multires(infile):
    print "MultiResolution phasevocoding for " + infile

    # Build the processing network
    mng = MarSystemManager()
    pvseries = mng.create("Series", "pvseries")
    pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseries.addMarSystem(mng.create("ShiftInput", "si"))
    pvseries.addMarSystem(mng.create("Windowing", "fo"))
    pvseries.addMarSystem(mng.create("Spectrum", "spk"))
    pvseries.addMarSystem(mng.create("PvConvert", "conv"))
    pvseries.addMarSystem(mng.create("PvOscBank", "pob"))
    pvseries.addMarSystem(mng.create("ShiftOutput", "so"))
    pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))

    pvseriesLong = mng.create("Series", "pvseriesLong")
    pvseriesLong.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseriesLong.addMarSystem(mng.create("ShiftInput", "si"))
    pvseriesLong.addMarSystem(mng.create("Windowing", "fo"))
    pvseriesLong.addMarSystem(mng.create("Spectrum", "spk"))
    pvseriesLong.addMarSystem(mng.create("PvConvert", "conv"))
    pvseriesLong.addMarSystem(mng.create("PvOscBank", "pob"))
    pvseriesLong.addMarSystem(mng.create("ShiftOutput", "so"))
    pvseriesLong.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
    inSamples = pvseries.getControl("mrs_natural/inSamples")
    awinSize = pvseries.getControl("ShiftInput/si/mrs_natural/winSize")
    swinSize = pvseries.getControl("PvOscBank/pob/mrs_natural/winSize")
    analysisHop = pvseries.getControl("PvConvert/conv/mrs_natural/Decimation")
    synthesisHop = pvseries.getControl(
        "PvOscBank/pob/mrs_natural/Interpolation")
    synthesisHop1 = pvseries.getControl(
        "ShiftOutput/so/mrs_natural/Interpolation")
    convertMode = pvseries.getControl("PvConvert/conv/mrs_string/mode")
    outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
    notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
    sinusoids = pvseries.getControl("PvConvert/conv/mrs_natural/Sinusoids")
    pitchShift = pvseries.getControl("PvOscBank/pob/mrs_real/PitchShift")
    zeroPhasing = pvseries.getControl("Windowing/fo/mrs_bool/zeroPhasing")
    wintype = pvseries.getControl("Windowing/fo/mrs_string/type")
    zeroPadding = pvseries.getControl("Windowing/fo/mrs_natural/zeroPadding")
    outData1 = pvseries.getControl("Spectrum/spk/mrs_realvec/processedData")

    filenameLong = pvseriesLong.getControl(
        "SoundFileSource/src/mrs_string/filename")
    inSamplesLong = pvseriesLong.getControl("mrs_natural/inSamples")
    outData = pvseriesLong.getControl(
        "ShiftInput/si/mrs_realvec/processedData")
    awinSizeLong = pvseriesLong.getControl("ShiftInput/si/mrs_natural/winSize")
    swinSizeLong = pvseriesLong.getControl("PvOscBank/pob/mrs_natural/winSize")
    analysisHopLong = pvseriesLong.getControl(
        "PvConvert/conv/mrs_natural/Decimation")
    synthesisHopLong = pvseriesLong.getControl(
        "PvOscBank/pob/mrs_natural/Interpolation")
    synthesisHop1Long = pvseriesLong.getControl(
        "ShiftOutput/so/mrs_natural/Interpolation")
    convertModeLong = pvseriesLong.getControl("PvConvert/conv/mrs_string/mode")
    outfnameLong = pvseriesLong.getControl(
        "SoundFileSink/dest/mrs_string/filename")
    notemptyLong = pvseriesLong.getControl(
        "SoundFileSource/src/mrs_bool/hasData")
    sinusoidsLong = pvseriesLong.getControl(
        "PvConvert/conv/mrs_natural/Sinusoids")
    pitchShiftLong = pvseriesLong.getControl(
        "PvOscBank/pob/mrs_real/PitchShift")
    zeroPhasingLong = pvseriesLong.getControl(
        "Windowing/fo/mrs_bool/zeroPhasing")
    wintypeLong = pvseriesLong.getControl("Windowing/fo/mrs_string/type")
    zeroPaddingLong = pvseriesLong.getControl(
        "Windowing/fo/mrs_natural/zeroPadding")
    outData1Long = pvseriesLong.getControl(
        "Spectrum/spk/mrs_realvec/processedData")

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = name + "_multires" + ext
    outsfnameLong = name + "_multiresLong" + ext

    zeroPhasing.setValue_bool(True)
    wintype.setValue_string("Hanning")
    zeroPadding.setValue_natural(3 * 512)
    awinSize.setValue_natural(512)
    swinSize.setValue_natural(2048)
    synthesisHop.setValue_natural(128)
    synthesisHop1.setValue_natural(128)
    analysisHop.setValue_natural(128)
    sinusoids.setValue_natural(120)
    pitchShift.setValue_real(0.66)
    convertMode.setValue_string("sorted")
    inSamples.setValue_natural(128)
    filename.setValue_string(infile)
    outfname.setValue_string(outsfname)

    zeroPhasingLong.setValue_bool(True)
    wintypeLong.setValue_string("Hanning")
    zeroPaddingLong.setValue_natural(0)
    awinSizeLong.setValue_natural(2048)
    swinSizeLong.setValue_natural(2048)
    synthesisHopLong.setValue_natural(128)
    synthesisHop1Long.setValue_natural(128)
    analysisHopLong.setValue_natural(128)
    sinusoidsLong.setValue_natural(120)
    pitchShiftLong.setValue_real(0.66)
    convertModeLong.setValue_string("sorted")
    inSamplesLong.setValue_natural(128)
    filenameLong.setValue_string(infile)
    outfnameLong.setValue_string(outsfnameLong)

    l
    xLong = zeros(2048, float)
    xcorr = zeros(2048, float)
    mpx = zeros(2048, float)
    correlations = zeros(13525, float)
    times = zeros(13525, float)
    iterations = 0
    time_in_seconds = 0.0
    while notempty.to_bool():
        iterations = iterations + 1
        if (iterations > 0):
            #			raw_input("Press enter to continue")
            timedata = outData.to_realvec()
            plotdata = outData1.to_realvec()
            plotdataLong = outData1Long.to_realvec()

            for i in range(0, len(plotdata)):
                x[i] = plotdata[i]
                xLong[i] = plotdataLong[i]
                xcorr[i] = x[i] * xLong[i]

            mx = abs(x)
            mxLong = abs(xLong)
            correlations[iterations] = mean(xcorr)
            times[iterations] = time_in_seconds
            df = mx - mpx
            df2 = mx - mxLong
            # subplot(321)
            # 			plot(mx[0:500], 'r-')
            # 			subplot(322)
            # 			plot(mxLong[0:500], 'b-')
            # 			subplot(323)
            # 			plot(timedata, 'b-');
            # 			subplot(324)
            # 			plot(df[0:500] * 1000, 'y-')
            # 			subplot(325)
            # 			plot(df2[0:500] * 1000, 'b-')
            #plot(mx[0:500], 'r-')
            #hold(True)
            #plot(mpx[0:500], 'g-')
            #hold(True)
            #plot(mxLong[0:500], 'b-')
            #			hold(False)
            #			show()
            mpx = mx

        pvseries.tick()
        pvseriesLong.tick()
        time_in_seconds = time_in_seconds + (128.0 / 44100.0)
        print time_in_seconds
        print iterations
    plot(times, correlations * 100000000.0, 'b-')
    show()
    raw_input("Hello")

Example #10

def identity(infile):
    print "Identity-phaselocking phasevocoding " + infile

    # Build the processing network
    mng = MarSystemManager()
    a = ("Series", "pvseries")
    pvseries = mng.create("Series", "pvseries")
    pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
    pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
    inSamples = pvseries.getControl("mrs_natural/inSamples")
    winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
    fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
    analysisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Decimation")
    synthesisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Interpolation")
    convertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/convertMode")
    unconvertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/unconvertMode")
    outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
    notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
    phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
    magHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes")
    peakHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks")
    pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos")

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = name + "_identity" + ext

    #identity phaselocking

    filename.setValue_string(infile)
    inSamples.setValue_natural(512)
    analysisHop.setValue_natural(512)
    synthesisHop.setValue_natural(1024)
    winSize.setValue_natural(4096)
    fftSize.setValue_natural(4096)
    convertMode.setValue_string("full")
    unconvertMode.setValue_string("identity_phaselock")
    outfname.setValue_string(outsfname)

    i = 0
    while notempty.to_bool():
        if (i == 0):
            phaselock.setValue_bool(True)

        i = i + 1
        pvseries.tick()


# 		if (i==100):
# 			magdata = magHandle.to_realvec()
# 			peakdata = peakHandle.to_realvec()
# 			hold(False)
# 			plot(magdata)
# 			hold(True)
# 			plot(peakdata)
# 			hold(False)
# 			show()

    print "Identity-phaselocking phasevocoding done. Output in " + outsfname

Example #11

File: pvoc.py Project: Amos-zq/marsyas

def multires(infile):
	print "MultiResolution phasevocoding for " + infile
	
	# Build the processing network
	mng = MarSystemManager()
	pvseries = mng.create("Series", "pvseries")
	pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseries.addMarSystem(mng.create("ShiftInput", "si"))
	pvseries.addMarSystem(mng.create("Windowing", "fo"));
	pvseries.addMarSystem(mng.create("Spectrum", "spk"));
	pvseries.addMarSystem(mng.create("PvConvert", "conv"));
	pvseries.addMarSystem(mng.create("PvOscBank", "pob"));
	pvseries.addMarSystem(mng.create("ShiftOutput", "so"));
	pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))


	pvseriesLong = mng.create("Series", "pvseriesLong")
	pvseriesLong.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseriesLong.addMarSystem(mng.create("ShiftInput", "si"))
	pvseriesLong.addMarSystem(mng.create("Windowing", "fo"));
	pvseriesLong.addMarSystem(mng.create("Spectrum", "spk"));
	pvseriesLong.addMarSystem(mng.create("PvConvert", "conv"));
	pvseriesLong.addMarSystem(mng.create("PvOscBank", "pob"));
	pvseriesLong.addMarSystem(mng.create("ShiftOutput", "so"));
	pvseriesLong.addMarSystem(mng.create("SoundFileSink", "dest"))


	
	filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
	inSamples = pvseries.getControl("mrs_natural/inSamples")
	awinSize = pvseries.getControl("ShiftInput/si/mrs_natural/winSize")
	swinSize = pvseries.getControl("PvOscBank/pob/mrs_natural/winSize")
	analysisHop = pvseries.getControl("PvConvert/conv/mrs_natural/Decimation")
	synthesisHop = pvseries.getControl("PvOscBank/pob/mrs_natural/Interpolation")
	synthesisHop1 = pvseries.getControl("ShiftOutput/so/mrs_natural/Interpolation")
	convertMode = pvseries.getControl("PvConvert/conv/mrs_string/mode")
	outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
	notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
	sinusoids = pvseries.getControl("PvConvert/conv/mrs_natural/Sinusoids");
	pitchShift = pvseries.getControl("PvOscBank/pob/mrs_real/PitchShift");
	zeroPhasing = pvseries.getControl("Windowing/fo/mrs_bool/zeroPhasing");
	wintype = pvseries.getControl("Windowing/fo/mrs_string/type");
	zeroPadding = pvseries.getControl("Windowing/fo/mrs_natural/zeroPadding");
	outData1  = pvseries.getControl("Spectrum/spk/mrs_realvec/processedData")



	filenameLong = pvseriesLong.getControl("SoundFileSource/src/mrs_string/filename")
	inSamplesLong = pvseriesLong.getControl("mrs_natural/inSamples")
	outData = pvseriesLong.getControl("ShiftInput/si/mrs_realvec/processedData");
	awinSizeLong = pvseriesLong.getControl("ShiftInput/si/mrs_natural/winSize")
	swinSizeLong = pvseriesLong.getControl("PvOscBank/pob/mrs_natural/winSize")
	analysisHopLong = pvseriesLong.getControl("PvConvert/conv/mrs_natural/Decimation")
	synthesisHopLong = pvseriesLong.getControl("PvOscBank/pob/mrs_natural/Interpolation")
	synthesisHop1Long = pvseriesLong.getControl("ShiftOutput/so/mrs_natural/Interpolation")
	convertModeLong = pvseriesLong.getControl("PvConvert/conv/mrs_string/mode")
	outfnameLong = pvseriesLong.getControl("SoundFileSink/dest/mrs_string/filename")
	notemptyLong = pvseriesLong.getControl("SoundFileSource/src/mrs_bool/hasData")
	sinusoidsLong = pvseriesLong.getControl("PvConvert/conv/mrs_natural/Sinusoids");
	pitchShiftLong = pvseriesLong.getControl("PvOscBank/pob/mrs_real/PitchShift");
	zeroPhasingLong = pvseriesLong.getControl("Windowing/fo/mrs_bool/zeroPhasing");
	wintypeLong = pvseriesLong.getControl("Windowing/fo/mrs_string/type");
	zeroPaddingLong = pvseriesLong.getControl("Windowing/fo/mrs_natural/zeroPadding");
	outData1Long  = pvseriesLong.getControl("Spectrum/spk/mrs_realvec/processedData")

	
	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = name + "_multires" + ext
	outsfnameLong = name + "_multiresLong" + ext


	zeroPhasing.setValue_bool(True);
	wintype.setValue_string("Hanning");
	zeroPadding.setValue_natural(3 * 512);
	awinSize.setValue_natural(512);
	swinSize.setValue_natural(2048);
	synthesisHop.setValue_natural(128);
	synthesisHop1.setValue_natural(128);
	analysisHop.setValue_natural(128);
	sinusoids.setValue_natural(120);
	pitchShift.setValue_real(0.66);
	convertMode.setValue_string("sorted");
	inSamples.setValue_natural(128);
	filename.setValue_string(infile);
	outfname.setValue_string(outsfname);


	zeroPhasingLong.setValue_bool(True);
	wintypeLong.setValue_string("Hanning");
	zeroPaddingLong.setValue_natural(0);
	awinSizeLong.setValue_natural(2048);
	swinSizeLong.setValue_natural(2048);
	synthesisHopLong.setValue_natural(128);
	synthesisHop1Long.setValue_natural(128);
	analysisHopLong.setValue_natural(128);
	sinusoidsLong.setValue_natural(120);
	pitchShiftLong.setValue_real(0.66);
	convertModeLong.setValue_string("sorted");
	inSamplesLong.setValue_natural(128);
	filenameLong.setValue_string(infile);
	outfnameLong.setValue_string(outsfnameLong);
	
	l
	xLong = zeros(2048, float)
	xcorr = zeros(2048, float)
	mpx = zeros(2048, float)
	correlations = zeros(13525, float)
	times=zeros(13525, float)
	iterations = 0
	time_in_seconds = 0.0;
	while notempty.to_bool():
		iterations = iterations + 1
		if (iterations>0):
#			raw_input("Press enter to continue")
			timedata = outData.to_realvec()
			plotdata = outData1.to_realvec()
			plotdataLong = outData1Long.to_realvec()

			
			for i in range(0, len(plotdata)):
				x[i] = plotdata[i]
				xLong[i] = plotdataLong[i]
				xcorr[i] = x[i] * xLong[i] 

			mx = abs(x)
			mxLong = abs(xLong)
			correlations[iterations] = mean(xcorr)
			times[iterations] = time_in_seconds
			df = mx - mpx;
			df2 = mx - mxLong;
			# subplot(321)
# 			plot(mx[0:500], 'r-')
# 			subplot(322)
# 			plot(mxLong[0:500], 'b-')
# 			subplot(323)
# 			plot(timedata, 'b-');
# 			subplot(324)
# 			plot(df[0:500] * 1000, 'y-')
# 			subplot(325)
# 			plot(df2[0:500] * 1000, 'b-')
			#plot(mx[0:500], 'r-')
			#hold(True)
			#plot(mpx[0:500], 'g-')
			#hold(True)			
			#plot(mxLong[0:500], 'b-')
#			hold(False)
#			show()
			mpx = mx
			
		pvseries.tick()
		pvseriesLong.tick()
		time_in_seconds = time_in_seconds + (128.0 / 44100.0);
		print time_in_seconds
		print iterations
	plot(times, correlations * 100000000.0, 'b-')
	show()
	raw_input("Hello")

Example #12

File: pvoc.py Project: Amos-zq/marsyas

def identity(infile):
	print "Identity-phaselocking phasevocoding " + infile 

	# Build the processing network
	mng = MarSystemManager()
	a = ("Series", "pvseries")
	pvseries = mng.create("Series", "pvseries")
	pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
	pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))
	
	filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
	inSamples = pvseries.getControl("mrs_natural/inSamples")
	winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
	fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
	analysisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Decimation")
	synthesisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Interpolation")
	convertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/convertMode")
	unconvertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/unconvertMode")
	outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
	notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
	phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
	magHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes"); 
	peakHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks");
	pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos");
	
	
	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = name + "_identity" + ext
	
    #identity phaselocking 

	filename.setValue_string(infile);
	inSamples.setValue_natural(512);
	analysisHop.setValue_natural(512);
	synthesisHop.setValue_natural(1024);
	winSize.setValue_natural(4096);
	fftSize.setValue_natural(4096);
	convertMode.setValue_string("full");
	unconvertMode.setValue_string("identity_phaselock");
	outfname.setValue_string(outsfname);
	
	i = 0 
	while notempty.to_bool():
		if (i == 0):
			phaselock.setValue_bool(True)
			
		i = i + 1
		pvseries.tick()
# 		if (i==100):
# 			magdata = magHandle.to_realvec()
# 			peakdata = peakHandle.to_realvec()
# 			hold(False)
# 			plot(magdata)
# 			hold(True)
# 			plot(peakdata)
# 			hold(False)
# 			show()


		
	print "Identity-phaselocking phasevocoding done. Output in "  + outsfname 

Example #13

File: pvoc.py Project: Amos-zq/marsyas

def scaled(infile):
	print "Scaled-phaselocking phasevocoding " + infile 

	# Build the processing network
	mng = MarSystemManager()
	pvseries = mng.create("Series", "pvseries")
	pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
	pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
	pvseries.addMarSystem(mng.create("StretchLinear", "sl"))
	pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))
	
	filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
	inSamples = pvseries.getControl("mrs_natural/inSamples")
	winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
	fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
	analysisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Decimation")
	synthesisHop = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/Interpolation")
	convertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/convertMode")
	unconvertMode = pvseries.getControl("PhaseVocoder/pvoc/mrs_string/unconvertMode")
	outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
	notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
	phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
	magHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes"); 
	peakHandle = pvseries.getControl("PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks");
	pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos");

	sfactor = pvseries.getControl("StretchLinear/sl/mrs_real/stretch");
	
	
	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = "temp.wav"
	
    #scaled phaselocking 

	filename.setValue_string(infile);
	inSamples.setValue_natural(512);
	analysisHop.setValue_natural(512);
	synthesisHop.setValue_natural(341);
	winSize.setValue_natural(2048);
	fftSize.setValue_natural(2048);
	convertMode.setValue_string("analysis_scaled_phaselock");
	unconvertMode.setValue_string("scaled_phaselock");
	sfactor.setValue_real(1.0);
	outfname.setValue_string(outsfname);
	
	i = 0 
	while notempty.to_bool():
		if (i == 0):
			phaselock.setValue_bool(True)
			
		i = i + 1
		pvseries.tick()


	[name, ext] = os.path.splitext(os.path.basename(infile))
	outsfname = name + "_scaled" + ext

	pnet = mng.create("Series", "pnet")
	pnet.addMarSystem(mng.create("SoundFileSource", "src"))
	pnet.addMarSystem(mng.create("StretchLinear", "sl"))
	pnet.addMarSystem(mng.create("SoundFileSink", "dest"))

	filename = pnet.getControl("SoundFileSource/src/mrs_string/filename")
	notempty = pnet.getControl("SoundFileSource/src/mrs_bool/hasData")
	sfactor = pnet.getControl("StretchLinear/sl/mrs_real/stretch")
	outfname = pnet.getControl("SoundFileSink/dest/mrs_string/filename")
	inSamples = pnet.getControl("mrs_natural/inSamples")
	
	sfactor.setValue_real(1.5);
	filename.setValue_string("temp.wav");
	outfname.setValue_string(outsfname);
	size = pnet.getControl("SoundFileSource/src/mrs_natural/size");
	inSamples.setValue_natural(size.to_natural())
	
	while notempty.to_bool():
		pnet.tick()

	print "Scaled-phaselocking phasevocoding done. Output in "  + outsfname 

Example #14

File: test.py Project: Amos-zq/marsyas

#!/usr/bin/env python


from marsyas import MarSystemManager,MarControlPtr

msm = MarSystemManager()

file = msm.create("SoundFileSource","file")
sink = msm.create("AudioSink","sink")
gain = msm.create("Gain", "gain")
pipe = msm.create("Series","pipe")
pipe.addMarSystem(file)
pipe.addMarSystem(gain)
pipe.addMarSystem(sink)

filename = pipe.getControl("SoundFileSource/file/mrs_string/filename")
notempty = pipe.getControl("SoundFileSource/file/mrs_bool/hasData")
iniAudio = pipe.getControl("AudioSink/sink/mrs_bool/initAudio")
outData  = pipe.getControl("Gain/gain/mrs_realvec/processedData")
import sys

for arg in sys.argv[1:] :
	filename.setValue_string(arg)

	if notempty.to_bool(): # If file has data
		iniAudio.setValue_bool(1)
		
		while notempty.to_bool(): # Play it
			pipe.tick()
			plotdata = outData.to_realvec()

Example #15

def scaled(infile):
    print "Scaled-phaselocking phasevocoding " + infile

    # Build the processing network
    mng = MarSystemManager()
    pvseries = mng.create("Series", "pvseries")
    pvseries.addMarSystem(mng.create("SoundFileSource", "src"))
    pvseries.addMarSystem(mng.create("PhaseVocoder", "pvoc"))
    pvseries.addMarSystem(mng.create("StretchLinear", "sl"))
    pvseries.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pvseries.getControl("SoundFileSource/src/mrs_string/filename")
    inSamples = pvseries.getControl("mrs_natural/inSamples")
    winSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/winSize")
    fftSize = pvseries.getControl("PhaseVocoder/pvoc/mrs_natural/FFTSize")
    analysisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Decimation")
    synthesisHop = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_natural/Interpolation")
    convertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/convertMode")
    unconvertMode = pvseries.getControl(
        "PhaseVocoder/pvoc/mrs_string/unconvertMode")
    outfname = pvseries.getControl("SoundFileSink/dest/mrs_string/filename")
    notempty = pvseries.getControl("SoundFileSource/src/mrs_bool/hasData")
    phaselock = pvseries.getControl("PhaseVocoder/pvoc/mrs_bool/phaselock")
    magHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/magnitudes")
    peakHandle = pvseries.getControl(
        "PhaseVocoder/pvoc/PvUnconvert/uconv/mrs_realvec/peaks")
    pos = pvseries.getControl("SoundFileSource/src/mrs_natural/pos")

    sfactor = pvseries.getControl("StretchLinear/sl/mrs_real/stretch")

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = "temp.wav"

    #scaled phaselocking

    filename.setValue_string(infile)
    inSamples.setValue_natural(512)
    analysisHop.setValue_natural(512)
    synthesisHop.setValue_natural(341)
    winSize.setValue_natural(2048)
    fftSize.setValue_natural(2048)
    convertMode.setValue_string("analysis_scaled_phaselock")
    unconvertMode.setValue_string("scaled_phaselock")
    sfactor.setValue_real(1.0)
    outfname.setValue_string(outsfname)

    i = 0
    while notempty.to_bool():
        if (i == 0):
            phaselock.setValue_bool(True)

        i = i + 1
        pvseries.tick()

    [name, ext] = os.path.splitext(os.path.basename(infile))
    outsfname = name + "_scaled" + ext

    pnet = mng.create("Series", "pnet")
    pnet.addMarSystem(mng.create("SoundFileSource", "src"))
    pnet.addMarSystem(mng.create("StretchLinear", "sl"))
    pnet.addMarSystem(mng.create("SoundFileSink", "dest"))

    filename = pnet.getControl("SoundFileSource/src/mrs_string/filename")
    notempty = pnet.getControl("SoundFileSource/src/mrs_bool/hasData")
    sfactor = pnet.getControl("StretchLinear/sl/mrs_real/stretch")
    outfname = pnet.getControl("SoundFileSink/dest/mrs_string/filename")
    inSamples = pnet.getControl("mrs_natural/inSamples")

    sfactor.setValue_real(1.5)
    filename.setValue_string("temp.wav")
    outfname.setValue_string(outsfname)
    size = pnet.getControl("SoundFileSource/src/mrs_natural/size")
    inSamples.setValue_natural(size.to_natural())

    while notempty.to_bool():
        pnet.tick()

    print "Scaled-phaselocking phasevocoding done. Output in " + outsfname

Example #16

#!/usr/bin/env python

from pylab import *
from marsyas import MarSystemManager, MarControlPtr

msm = MarSystemManager()

weka = msm.create("WekaSource", "weka")
classifier = msm.create("GaussianClassifier", "classifier")
summary = msm.create("Summary", "summary")
pipe = msm.create("Series", "pipe")

pipe.addMarSystem(weka)
pipe.addMarSystem(classifier)
pipe.addMarSystem(summary)

filename = pipe.getControl("WekaSource/weka/mrs_string/filename")
validationmode = pipe.getControl("WekaSource/weka/mrs_string/validationMode")

validationmode.setValue_string("PercentageSplit,50%")
filename.setValue_string("ms.arff")

insamples = pipe.getControl("mrs_natural/inSamples")
insamples.setValue_natural(1)

wekanclasses = pipe.getControl("WekaSource/weka/mrs_natural/nClasses")
summarynclasses = pipe.getControl("Summary/summary/mrs_natural/nClasses")
classifiernclasses = pipe.getControl(
    "GaussianClassifier/classifier/mrs_natural/nClasses")
summarynclasses.setValue_natural(wekanclasses.to_natural())
classifiernclasses.setValue_natural(wekanclasses.to_natural())

Example #17

#!/usr/bin/env python

from marsyas import MarSystemManager, MarControlPtr

msm = MarSystemManager()

file = msm.create("SoundFileSource", "file")
sink = msm.create("AudioSink", "sink")
gain = msm.create("Gain", "gain")
pipe = msm.create("Series", "pipe")
pipe.addMarSystem(file)
pipe.addMarSystem(gain)
pipe.addMarSystem(sink)

filename = pipe.getControl("SoundFileSource/file/mrs_string/filename")
notempty = pipe.getControl("SoundFileSource/file/mrs_bool/hasData")
iniAudio = pipe.getControl("AudioSink/sink/mrs_bool/initAudio")
outData = pipe.getControl("Gain/gain/mrs_realvec/processedData")
import sys

for arg in sys.argv[1:]:
    filename.setValue_string(arg)

    if notempty.to_bool():  # If file has data
        iniAudio.setValue_bool(1)

        while notempty.to_bool():  # Play it
            pipe.tick()
            plotdata = outData.to_realvec()