%pylab inline
from __future__ import print_function
from __future__ import division
import IPython
Populating the interactive namespace from numpy and matplotlib
img=imread('swirl.jpg')
imshow(img)
<matplotlib.image.AxesImage at 0x5bf49d0>
img.shape
(720, 1280, 3)
subplot(131)
imshow(img[:,:,0]) #Red
colorbar()
subplot(132)
imshow(img[:,:,1]) #Blue
colorbar()
subplot(133)
imshow(img[:,:,2]) #Green
colorbar()
gcf().set_figwidth(18)
img.dtype
dtype('uint8')
sig = img[350,:, 1]
plot(sig)
sig.dtype
dtype('uint8')
import wave
from scipy.io import wavfile
sigout = list(sig.astype(uint8))*50
plot(sigout)
#wavfile.write('swirl1line2.wav', 8000, array(sigout, dtype=uint8))
gcf().set_figwidth(18)
xlim(0,5000)
IPython.display.Audio(data=sigout, filename=None, url=None, embed=True, rate=8000, autoplay=False)
len(sigout)
64000
wt = linspace(0, 6 * pi, 64000)
oscillation = 125*cos(wt)
plot(oscillation);
gcf().set_figwidth(20)
plot(oscillation)
audioout=sigout*oscillation
plot(audioout);
gcf().set_figwidth(20)
#wavfile.write('swirl3final.wav', 8000, array(audioout, dtype=int16))
IPython.display.Audio (data=audioout, filename=None, url=None, embed=True, rate=8000, autoplay=False)
from scipy.special import jn
x = linspace(0,64000, 64000)
plot( x,jn(48,x))
osc=1000*jn(48,x)
xlim(0,1000)
gcf().set_figwidth(20)
plot(osc)
audioout=sigout*osc
plot(audioout);
gcf().set_figwidth(20)
xlim(0,25000)
(0, 25000)
#wavfile.write('swirl3final1.wav', 8000, array(audioout, dtype=int16))
IPython.display.Audio (data=audioout, filename=None, url=None, embed=True, rate=8000, autoplay=False)