IAJIT First Online Publication

Transcription

1 An Improved Feaure Exracion and Combinaion of Muliple Classifiers for Query-by- ming Naha Phiwma and Parinya Sanguansa 2 Deparmen of Compuer Science, Suan Dusi Rajabha Universiy, Thailand 2 Faculy of Engineering and Technology, Panyapiwa Insiue of Managemen, Thailand IAJIT Firs Online Publicaion Absrac: In his paper, we propose new mehods for feaure exracion and sof majoriy voing o adjus efficiency and accuracy of music rerieval. For our work, he inpu is humming sound which is sound wave and Musical Insrumen Digial Inerface (MIDI) is used as he reference song in daabase. A criical issue of humming sound are variaion such as duraion, sound, empo, key, and noise inerference from boh environmen and acquisiion insrumens. Besides all he problems of humming sound we have menioned earlier, wheher humming sound and MIDI in differen domain which will make he difficuly for wo domains o compare each oher. However, o make hese wo in he same domain, we conver hem ino he frequency domain. Our approach sars from pre-processing by using feaures for noe segmenaion by humming sound. The process consiss of four seps as follows: Firsly, he MIDI is already a sequence of pich while he pich in humming sound is needed o exrac by Subharmonic-o-Harmonic (SHR). Subsequenly, he exraced pich can be used o calculae all above aribues and hen muliple classifiers are applied o classify he muliple subses of hese feaures. Aferwards, he subse conain he muliple aribues, Muli-Dimensional Dynamic Time Warping (MD-DTW) is used for similariy measuremen. Finally, Neares Neighbours (NN) and sof majoriy voing are used o obain he rerieval resuls in case of equal scores. From he experimens, o achieve 00% accuracy rae a he early op-n rank in rerieving, he appropriae feaure se should consis of five classifiers. Keywords: Query-by-humming, feaure exracion, majoriy voing, muliple classifiers, MD-DTW, SHR. Received February 8, 202; acceped May 22, 202. Inroducion A presen, he music becomes par of our lives boh lisening and singing o enerain and relax ourselves. The prevalen of problem, mos users forge he name of he song, bu hey wan o find a song for lisening and singing. However, radiional approaches for rerieving music daa were based on he exual informaion such as iles, composers, file names or singers. Because of heir incompleeness, here are many difficulies in saisfying paricular requiremens of applicaions. Therefore, many researchers have proposed echniques o query a song base on humming which is called Query-by-ming (QBH) sysem [2, 3, 4, 20, 2, 26]. QBH sysem allows he user o rerieve an inended song based on humming some par of he song. The general framework of QBH sysem conains hree main componens, which are query processing module, melody daabase and maching engine [2]. Firsly, he sysem handles he Musical Insrumen Digial Inerface (MIDI) in daabase. Subsequenly, he sysem process he users inpu humming signal hen exrac signal fundamenal frequency, humming query is convered ino melody represenaion. Finally, when a search is iniiaed, melody represenaion is used o mach agains he melody in he feaure daabase, according o heir similariies and reurn a rank lis of songs. Normally, naural sounds are a composiion of a fundamenal frequency wih a se of harmonics. The frequency ha he human ear inerpres as he pich of a sound is his fundamenal frequency, even if i is absen in he sound. The pich of naural sounds is imporan in many conexs. Pich is he percepion of how high or low a musical noe sounds, which can be considered as a frequency which corresponds closely o he fundamenal frequency or main repeiion rae in he signal [5]. I is one of he mos imporan parameers in he voice signal analysis and can be deermined by he fundamenal frequency of he uni frame [6]. For QBH sysem, pich is he key feaure of melody. As he humming sound consiss of noise, pich needs o be exraced and in order o ge he mos significan informaion. 2. Relaed Works For early work, pich is used in QBH as a feaure [2], [8, 4]. There are many echniques o analyse and

2 exrac pich conour, pich inerval and duraion from voice humming query [2]. In general, radiional mehods for deecing piches have been proposed in he pas, i can be divided roughly ino wo domains o idenify he pich: ime-domain based, frequency domain based [, 5]. Pich and fundamenal frequency are imporan feaures, herefore i mus be exraced pich. A pich deerminaion algorihm (PDA) based on Subharmonico-Harmonic Raio (SHR) is developed in he frequency domain and describes he ampliude raio beween subharmonics and harmonics [23], [24]. For our sysem, we have implemened pich racking using SHR. The Mel Frequency Cepsral Coefficiens (MFCC) was adoped in many speech analysis applicaions. This ype of feaure exracion is being widely used in robus speech recogniion sysems inspired by human audiory percepion and focusing on effecive signal processing in he ear using cochlear filerbanks []. MFCCs were also used as feaures [7], [2]. From hese experimens i shows ha using MFCC wih he dimension 3 and audio recogniion will give beer resuls han oher dimensions. MFCC is used in our pre-processing. Generally, o gaher all aribues o use all a once migh no give good resul. Some feaures are appropriae bu some are no. However, we need o find many classifiers o help wih he resul. To improve he performance of he sysem, here are many researches used a lo of informaion, such as pich, duraion, rhyhm, iner-onse inerval, sar and end ime o be muually considered and make feaure in [0], [4]. The mos ofen used classifiers combinaion approaches in Muliple Classifiers Sysem (MCS) include classifier selecion, he majoriy voing, he weighed combinaion (weighed averaging), he probabilisic schemes, various rank-ordered rules and ec. [4]. Besides, MCS and majoriy voe is applied for off-line Arabic handwriing recogniion, he accuracy is higher han individual classifier [9]. Therefore, MCS will be used o find he resul and he easies way o do is majoriy voing. However, he feaure sill has variable lengh in he form of melody conour, hence he radiional Dynamic Time Warping (DTW) canno be used for his feaure. Muli-dimensional ime series (MD-DTW) algorihm was proposed for DTW on muli-dimensional ime series, which he algorihm uilises all dimensions o find he bes synchronizaion [3]. Muli-dimensional (ime) series are series in which muliple measuremens are made simulaneously. MD-DTW is applied wih image exure [9], gesure recogniion [3], ime series [25], hus we have an idea o apply his o he QBH. The segmen of a noe in he humming waveform is model by a Hidden Markov model (HMM) while he pich of he noe is model by a pich model using a Gaussian mixure model. The frame based analysis is performed on a noe segmen which usually has several frames. Muliple frames of a segmened noe are used for pich model analysis. Afer applying auocorrelaion o hose frames, pich feaures are exraced. The firs sage of he proposed algorihm is noe segmenaion, where he process of segmening noes of a humming piece is conduced. Firs, a feaure se which can characerize a noe is chosen. Nex, he HMM definiion is chosen before raining. During he raining phase, noes phone level HMMs are rained using he seleced feaure se. The rained noe models are hen used by he noe decoder for noe segmenaion. Finally, he duraion of a segmened noe is label according o is relaive duraion change [5, 8, 9, 22]. IAJIT Firs Online Publicaion 3. Maerials and Mehods 3.. Melody Conour Exracion Algorihm The following algorihm describes how o exrac pich from humming sound o obain he melody conour. Melody Conour Exracion, we have proposed in [6]. Le m represens melody conour and le p be he pich. The variables of algorihm are described as follows: s is he size of he window for filering, g is he gap of pich difference, T is he hreshold of sandard deviaion, and v is he variance of pich inerval. The Algorihm proceeds as follows: Require: p, g, T, s Ensure: m Sep : smoohing p by median filer. Sep 2: iniial m p Sep 3: N lengh of p Sep 4: j Sep 5: while N do d = p p Y { p v, p v+,..., p+ v, p+ v} SY Sandard deviaion of Y if d > g and SY < T hen mj p end if + s j j+ end while Sep 6: reurn m The firs sep of his echnique is o ake a pich o pass hrough he noise filering process which uses he median filer in order o make he signal go smoohly. Then, find he differen value of p by comparing wih he defined g value by selecing only he exceed value. The value of s is deermined in order o apply o find he range of signal ha changes a lile a ha period of ime. In oher words, i discards he signal ha changes rapidly in a shor ime comparing wih his

3 inerval. There is he spread around he signal and i only needs he group of significan signals. Hence, i finds he range of signal which has a small value of he spread when comparing wih he hreshold of sandard deviaion (T) as shown in Figure. The oupu of he algorihm melody conour conains significan pich. Finally, when his echnique is applied o rerieval ask, i o do rerieval process, he resul will be more correc han he radiional mehod. Figure. Example of pich exracion by Melody conour exracion Noe Segmenaion by ming Sound For his paper, we have proposed he mehod of noe segmenaion by humming sound o differeniae he sounds par from he silence pars in order o choose he mos imporan par, which is he sound par, o use in he nex process [7]. From he sound wave in Figure 2, he silence inerval is removed manually as preprocessing before being fed o he HMM. Figure 2. Sound wave from humming sandard noe in C major scale (do, re, me,...,do). Figure 3. Noe Model and Silence Model. As shown in Figure 3, he HMM conain 3 saes wih lef-o-righ opology using 2 Gaussian mixure disribuions. Boh he noe and he silence are used o rain hese HMMs Muli-Dimensional Dynamic Time Warping (MD-DTW) IAJIT Firs Online Publicaion Muli-dimensional series consis of a number of measuremens made a each insance. The number of measuremens is he dimensionaliy of he series, he number of ime insances is lengh. Noe ha mulidimensional series need no be ime signals, any siuaion in which several measuremens are made simulaneously depending on one variable ha gives a mulidimensional series. They assume ha measuremens are sored in a marix, in which columns are feaures and rows are ime insances. Muli-Dimensional Dynamic Time Warping (MD- DTW) was proposed [3], an approach o calculae he DTW by synchronizing muli-dimensional series, which is basically an exension of he original DTW, where he marix D is creaed by compuing he disance beween k-dimensional poins (where, differenly from he original approach, k can be larger han ). This approach pre-processes he mulidimensional series, which mus have he same number of dimensions. The las sep of his algorihm is he execuion of he radiional DTW. However, in many cases, all dimensions will conain informaion needed for synchronisaion herefore propose MD-DTW for synchronising such series. The MD-DTW algorihm runs in 4 seps: Le A, B be wo series of dimension K and lengh M, N respecively. Sep : Normalize each dimension of A and B separaely o a zero mean and uni variance Sep 2: If desired, smooh each dimension wih a Gaussian filer Sep 3: Fill he M by N disance marix K D according o: D( i, j) = Ai( k) Bj( k). k = Sep 4: Use his disance marix o find he bes synchronizaion wih he regular DTW algorihm. Take wo series A and B. DTW involves he creaion of a marix in which he disance beween every possible combinaion of ime insances A( i) B( j) is sored. This disance is calculaed in erms of he feaure values of he poins. Various norms are possible. In D-DTW, he disance is usually calculaed by aking he absolue or he squared disance beween he feaure values of each combinaion of poins. For MD-DTW, a disance measure for wo K- dimensional poins mus be calculaed. This disance can be any p-norm. They use he -norm, i.e. he sum

4 of he absolue differences in all dimensions. To combine differen dimensions in his way, i is necessary o normalize each dimension o a zero mean and uni variance. For his, he dimensions mus be comparable. If for insance one dimension conains real valued measuremens and one is binary, comparing hem direcly is no possible and a more sophisicaed disance measure mus be found [3, 9]. 4. Our Approach In our approach, i consiss of wo seps which are preprocessing and processing. Pre-processing is process of noe segmenaion by humming sound. Processing, i consiss of feaure exracion and sof majoriy voing. F2. Normalized Duraion of ime: N ime ( T ) = (2) T T where represens noe duraions in seconds and T is he summaion of duraion ime. F3. Melody conour exracion (Melslope) M ( p ) = melslope( p ) (3) IAJIT Firs Online Publicaion F4. Sring numeric relaive (UDR) 0, p p < ε udr ( p) =, p p > ε 2, p p < ε 4.2. Feaure Exracion of MIDI (4) Conducing feaure exracion of MIDI has four approaches as following: Normalized Pich MIDI P ( p) = N( p) (5) Normalized Duraion of ime MIDI 2 (6) P T N T ime Figure 4. Block diagram of our approach. We propose wo mehods in our framework as shown in Figure 4 which are feaure exracion (Secion 4., 4.2 and 4.3) and majoriy voing exracion (Secion 4.4). Our framework sars from preprocessing by using a feaure o faciliae noe segmenaion by a humming sound. The process consiss of four seps as follows: Firsly, he MIDI is already a sequence of pich while he pich in humming sound is needed o exrac by Subharmonic-o- Harmonic (SHR) [23, 24]. Consequenly, he pich is exraced by our new feaure exracion mehod and hen muliple classifiers are applied o classify he muliple subses of hese feaures. Aferwards, MD- DTW is used for similariy measuremen. Finally, Neares Neighbors (NN) and sof majoriy voing are used o obain he rerieval resuls in case of equal scores. 4. Feaure Exracion In his process, he principle funcion used for making feaure exracion of inpu (humming sound) and reference songs in daabase (MIDI). F. Normalized Pich N ( p ) log p log p = logσ p () Normalized Duraion of pich MIDI 3 (7) P p N p ime Sring numeric relaive (UDR) MIDI 4 (8) P p udr p 4.3. Feaure Exracion of Inpu Conducing feaure exracion of hum has six approaches as following: Normalized Pich (9) P p N p Normalized Duraion of ime 2 (0) P p N T ime Normalized Duraion of pich 3 () P p N p ime Sring numeric relaive (UDR) 4 (2) P p udr p Melody Conour Exracion (Melslope) 7 (3) P p M p Melslope of pich pass noe segmenaion 8 (4) P pseg M pseg

5 Where pseg represens pich passed noe segmenaion. The MIDI is already a sequence of pich while he pich in humming sound is needed o exrac by SHR [23, 24]. Due o he difference characerisic of MIDI and humming sound, he feaure exracion P o P 6 are performed o boh of hem as P 7 o P 0 are only used for humming sound. While P, P8,, P 0 MIDI MIDI MIDI MIDI are compared wih P, P, P5, P 6 respecively. Aferwards, he exraced pich can be used o calculae all above aribues and hen muliple classifiers are applied o classify he muliple subses of hese feaures. In case of he subse conain he muliple aribues, MD-DTW is used insead of DTW for similariy measuremen. Finally, Neares Neighbours (NN) and sof majoriy voing are used o obain he rerieval resuls Sof Majoriy Voing Majoriy voing mehod is widely used in many asks classificaion. Voing is a mehod for a group o make a decision. By he principle of voing, in general, he final decision is based on highes score. Neverheless, in erms of equal voe, here are many ways of making decision, depending on paricular siuaion. Thus, we propose o make imporan decisions if he voe is equal, based on he principle of minimum disance, which i is called sof majoriy voing. daabase, feaure exracion, op-n rank and combinaion of feaure. 5.. Daase In his daase, here are 500 MIDI forma songs and hey are divided ino hree subses which are 00, 300 and 500. We used 00 ess humming sound o query songs in daabase. The es query is a humming sound which consiss of hummed unes wih Da Da Da. We used 00 humming sounds from differen people o es our sysem. The recording was done a 8 khz sampling rae, mono and ime duraion 0 seconds, sar a he beginning of song. The resuls are showed ha when he number of MIDI in daabase is smaller, he accuracy rae is higher. IAJIT Firs Online Publicaion 5.2. Variaion of Feaure Ses In his paper, some single aribue are used for creaing muliple aribues, such as P, P2, P4 and P7, as described in Tables -2. Experimens have shown he effeciveness of he sysem and according o he various condiions such as he variaions of number of songs in daabase, feaure exracion, op-n rank and combinaion of feaure. Name P P2 P3 P4 P7 P8 Table. Lis of single aribue. Single Aribue Descripion Normalized Pich Normalized Duraion Time Normalized Duraion of Pich Sring Numeric Relaive (UDR) Melody Conour Exracion (Melslope) Melslope of pich pass noe segmenaion Table 2. Lis of muliple aribues. Name P5 P6 Muliple Aribues Descripion P,P2 P,P4 P2,P7 P4,P7 Figure 5. Sof Majoriy Voing mehod. From Figure 5, if here is he only one highes score, i will reurn he class ha has he highes score as he resul. Bu if here are muliple highes scores, all members will be reconsidered minimum disance by sof majoriy voing mehod. The principle of sof majoriy voing makes all members reconsidered hrough disance. 5. Resuls We have conduced exensive experimens o measure rerieval performance in erms of accuracy. Experimens have shown he effeciveness of he sysem and according o he various condiions. For effeciveness of his sysem, he measures were seup o explore such as he variaion of number of songs in In his experimen, he number of classifiers was varied from wo o en. The feaure ses for each classifier are defined in Table 4. I is fixed as P7 in every classifiers. Since, our experimen, we have found ha combining aribue P7 wih oher feaures can achieve 00 % accuracy rae, which faser han using random mehod. The deail of classifier is used in experimenal as shown in Tables 3-4. Table 3 shows he aribue of each feaure. Table 4 conains a feaure se ha is used as each classifier.

6 Table 3. Lis of feaures. Feaure Aribues Feaure Aribues Feaure Aribues P 4 P2 P3 P8 27 P P3 P6 P8 2 P2 5 P2 P4 28 P P4 P6 P8 3 P3 6 P2 P4 P6 29 P P6 P8 4 P4 7 P2 P6 30 P2 P3 P4 P6 5 P5 8 P2 P8 3 P2 P5 P6 P8 6 P6 9 P3 P6 P8 32 P2 P4 P8 7 P7 20 P4 P6 P8 33 P3 P4 P5 P6 8 P8 2 P P2 P4 P4 P5 P6 P8 34 P P P2 P3 P3 P4 P8 35 P P P4 P6 P2 P4 P5 P6 36 P8 P8 P2 P4 24 P2 P4 P8 P3 P4 P5 P P2 P6 25 P2 P7 P8 P P2 P3 P4 38 P5 P6 P8 3 P P6 26 P3 P5 P8 P2 P3 P4 P6 39 P8 Table 4. Lis of classifiers. # Classifier Feaures se The performance evaluaions vary op-n from op- o op-60. The experimens are shown in Tables 5-7 for each daase. The resuls of using five classifiers give he bes performance in all daases. Tha is, i can achieve 00% of op-n which are op-0 in case of 00 MIDI songs, op-20 in case of 300 MIDI songs, and op-35 in case of 500 MIDI songs. The feaure se of five classifiers consiss of feaures, , as shown in Table 4. Tha is he ses of aribues {P7}, {P2, P3, P8}, {P2, P4, }, {P2, P4, P6} and {P, P4, P6, P8, } as shown in Table 3. From his resul, we found ha P5 is no included in his se. While P8, which noe segmenaion was processed, is employed. Feaure ses ha use wo classifiers can achieve 00% accuracy rae a op 60 or feaure se of en classifiers can achieve 00 % accuracy rae a op-30 for MIDI 00 songs in daabase, as shown in Table 5. Meanwhile, if MIDI songs in daabase increase feaure se of 2 classifiers can only achieve 97% accuracy rae while feaure se of 0 classifiers can achieve 00 % accuracy rae a op-50 and op-60 for MIDI 300 and 500 songs in daabase, as shown in Tables 6-7 and Figures 6-8. In addiion, query ime is used o measure he complexiy of our proposed echnique as shown in Table 8. We performed all he ess on a noebook wih a CPU of Inel Core 2 Duo processor 2.26 GHz, 2 GB of RAM. Normally, MCS wih more classifiers ake more query imes. 6. Conclusions In his paper, we propose new mehod for feaure exracion and sof majoriy voing o make imporan decision if he voe is equal in applicaion of QBH. Our approach consiss of wo processes which make humming sound go hrough noe segmenaion and hen exrac he feaure o creae many feaure ses by using six approaches for inpu and four approaches for MIDI. The main feaure we use in each se, i obains from melody conour exracion algorihm, which we have proposed earlier. Nex, sof majoriy voing will be used for making a decision o choose he bes resul. The advanage of our approach is o increase efficiency and accuracy for rerieving daa. From he use of muliple classifiers sysem by using sof majoriy voing as we have proposed, if he score is equal, all he members will ge o reconsider by finding minimum disance, which we can look a his as an advanage. Moreover, using more han one feaure can achieve beer accuracy rae han one feaure because of including more informaion and obaining muliple aspecs of ha. From he experimens, using feaure se which consiss of 5 classifiers will ge 00% accuracy a he early op-n rank in rerieving. Neverheless, using a greaer number of classifiers makes he sysem higher complexiy and longer query ime. IAJIT Firs Online Publicaion Table 5. Tes resuls of experimen wih 00 MIDI songs wih variaions of feaure ses. Top-n Classifier Rae (%)

7 Table 6. Tes resuls of experimen wih 300 MIDI songs wih variaions of feaure ses. Top-n Classifier Rae (%) Figure 7. The performance of feaure ses wih 300 MIDI songs. IAJIT Firs Online Publicaion Table 7. Tes resuls of experimen wih 500 MIDI songs wih variaions of feaure ses. Top-n Classifier Rae (%) Table 8. Tes resuls of query ime. # Classifier Query Time (Second) 00 MIDI songs 300 MIDI songs 500 MIDI songs Figure 6. The performance of feaure ses wih 00 MIDI songs. Figure 8. The performance of feaure ses wih 500 MIDI songs. Acknowledgemens This work was assised by Suan Dusi Rajabha Universiy hrough suppor wih a scholarship and Rangsi Universiy by providing he laboraory room for daa processing. We would like o hanks all people who fain hummed a lo of unes for us. Addiionally, he invaluable recommendaion and supervision from he anonymous reviewers are much appreciaed. References [] Behroozmand R., and Almasganj F., Comparison of neural neworks and suppor vecor machines applied o oimized feaures exraced from paiens' speech signal for classificaion of vocal fold inflammaion, Proc. Inernaional Symposium on Signal Processing and Informaion Technology, pp , [2] Ghias A., e al., Query by humming: musical informaion rerieval in an audio daabase, Proc. 3rd ACM inernaional conference on Mulimedia, New York, NY, USA, pp , 995. [3] Hol en G.A., Reinders M.J.T., and Hendriks E.A., Muli-dimensional dynamic ime warping for gesure recogniion, Proc. Thireenh annual conference of he Advanced School for

8 Compuing and Imaging, The Neherlands, pp , [4] Jiangao H. and Minghui W., Dynamic combinaion of muliple classifiers based on normalizing decision space, Proc. Inernaional Conference on Informaion Engineering (ICIE), vol., pp , 200. [5] Jing Q., e al., A novel mir approach based on dynamic hresholds segmenaion and weighed synhesis maching, Proc., IET Conference on Wireless, Mobile and Sensor Neworks, pp , [6] Jun Han B., Rho S., and Hwang E., An efficien voice ranscripion scheme for music rerieval, Proc. Inernaional Conference on Mulimedia and Ubiquious Engineering, pp , [7] Kim H. G., and Sikora T, Audio specrum projecion based on several basis decomposiion algorihms applied o general sound recogniion and audio segmenaion, in XII European Signal Processing Conference, Vienna, Ausria, pp , [8] Kosugi N., e al., A pracical query-by-humming sysem for a large music daabase, Proc. 8h ACM inernaional conference on Mulimedia, New York, NY, USA, pp , [9] Leila C., Maamai K., and Salim C., Combining Neural Neworks for Arabic Handwriing Recogniion, in The Inernaional Arab Journal of Informaion Technology, IAJIT. vol. 9, 20. [0] Lemsrom K., Laine P., and Peru S., Using relaive inerval slope in music informaion rerieval, Proc. ICMC 999, pp , 999. [] Li P., e al., A novel mir sysem based on improved melody conour definiion, Proc. Inernaional Conference on Mulimedia and Informaion Technology, pp , [2] Liu Y., e al. The sudy of he classificaion of Chinese folk songs by regional syle, Proc. Inernaional Conference on Semanic Compuing (ICSC), Washingon, DC, USA, pp , [3] Lu L., You H., and Zhang H.J., A new approach o query by humming in music rerieval, Proc. Inernaional Conference on Mulimedia and Expo, pp , 200. [4] McNab R. J., e al., Towards he digial music library: une rerieval from acousic inpu, Proc.s ACM inernaional conference on Digial libraries, New York, NY, USA, pp. 8, 996. [5] Mcnab R.J., Smih L. A., and Wien I.H, Signal processing for melody ranscripion, Proc. 9h Ausralasian Compuer Science Conf., 30307, pp , 996. [6] Phiwma N. and Sanguansa P., A music informaion sysem based on improved melody conour exracion, Proc. Inernaional Conference on Signal Acquisiion and Processing, pp.85 89, 200. [7] Phiwma N. and Sanguansa P., An improved noe segmenaion and normalizaion for Queryby-ming, in RANGSIT JOURNAL OF ARTS AND SCIENCES. vol., no.2, pp.39-48, 20. [8] Raphael C., Auomaic segmenaion of acousic musical signals using hidden markov models, IEEE Transacions on Paern Analysis and Machine Inelligence, vol. 2, pp , 998. [9] De Mello R.F. and Gondra I., Mulidimensional dynamic ime warping for image exure similariy, SBIA 2008, pp , [20] Ryyndnen M., and Klapuri A., Query by humming of midi and audio using localiy sensiive hashing, Proc. Inernaional Conference on Acousics, Speech and Signal Processing, pp , [2] Shih H. H., Narayanan S. S., and Kuo C. C. J, A saisical mulidimensional humming ranscripion using phone level hidden markov models for query by humming sysems, Proc. Inernaional Conference on Mulimedia and Expo, vol., pp. I 6 4, [22] Shih H. H., Narayanan S. S., and Kuo C. C. J., Mulidimensional humming ranscripion using a saisical approach for query by humming sysems, Proc. Inernaional Conference on Acousics, Speech,and Signal Processing, pp , [23] Sun X, A pich deerminaion algorihm based on subharmonico-harmonic raio, Proc. 6h Inernaional Conference of Spoken Language Processing, pp , [24] Sun X., Pich deerminaion and voice qualiy analysis using subharmonic-o-harmonic raio, Proc. Inernaional Conference on Acousics, Speech, and Signal, Orlondo, Florida, USA, pp , [25] Vlachos M., e al., Indexing muli-dimensional ime-series wih suppor for muliple disance measures, pp , [26] Zhu Y., Xu C., and Kankanhalli M., Melody curve processing for music rerieval, Proc. Inernaional Conference on Mulimedia and Expo, pp , 200. IAJIT Firs Online Publicaion

9 Naha Phiwma received he Ph.D. degrees in Informaion Technology from Rangsi Universiy, Thailand in 20. She is an assisan professor in he Deparmen of Compuer Science a Suan Dusi Rajabha Universiy, Thailand. Her research areas are music informaion rerieval and digial signal processing. IAJIT Firs Online Publicaion Parinya Sanguansa received he B.Eng., M.Eng. and Ph.D. degrees in Elecrical Engineering from he Chulalongkorn Universiy, Thailand. He is an assisan professor in he Faculy of Engineering and Technology, Panyapiwa Insiue of Managemen, Thailand in 200, 2004 and 2007 respecively. His research areas are digial signal processing in paern recogniion including on-line handwrien recogniion, face and auomaic arge recogniion.