CN105374360A - Interleaved additive spread spectrum audio watermark embedding method and detection method and system - Google Patents
Interleaved additive spread spectrum audio watermark embedding method and detection method and system Download PDFInfo
- Publication number
- CN105374360A CN105374360A CN201510833111.5A CN201510833111A CN105374360A CN 105374360 A CN105374360 A CN 105374360A CN 201510833111 A CN201510833111 A CN 201510833111A CN 105374360 A CN105374360 A CN 105374360A
- Authority
- CN
- China
- Prior art keywords
- signal
- signal frame
- watermark
- frequency
- domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/018—Audio watermarking, i.e. embedding inaudible data in the audio signal
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Signal Processing For Digital Recording And Reproducing (AREA)
Abstract
The invention discloses an interleaved additive spread spectrum audio watermark embedding method and detection method and system. The embedding method comprises the steps that time domain audio signals are sub-framed and then converted into frequency domain signals; the frequency domain signals are divided into odd-even interleaved frames; positive watermark bits and negative watermark bits are embedded in the odd number of signal frames and the even number of signal frames respectively by adopting an additive spread spectrum watermark embedding method; and the odd number of signal frames and the even number of signal frames with embedding of watermarks are combined into a frame of signals. The detection method comprises the steps that received signals are divided into odd-even interleaved frames; related detection is performed on the odd number of signal frames and the even number of signal frames respectively so that odd decision statistic quantity and even decision statistic quantity are obtained; difference of the odd decision statistic quantity and the even decision statistic quantity acts as decision statistic quantity; and the watermarks are extracted according to the decision statistic quantity. Variance of the decision statistic quantity can be greatly reduced so that error rate of watermark detection can be reduced.
Description
Technical field
The present invention relates to Digital Audio Watermarking Techniques field, particularly relate to a kind of additivity spread spectrum audio frequency watermark embedding grammar, detection method and system of intersecting.
Background technology
Digital audio frequency watermark in sound signal, adds some digital watermark information to reach the object signal processing operations such as file True-false distinguish, copyright protection, Information hiding.Intersect and refer to that odd indexed point frame in a frame audio sample point and even number sequence number point frame are cross arrangement distributions.The short-term stationarity of sound signal determines odd indexed point frame and even number sequence number point frame has stronger correlativity.In traditional additivity spread spectrum audio frequency watermark system, owing to not utilizing the high correlation of intersecting between parity frame, watermark signal after spread spectrum is directly added with main signal by it, during the employing coherent detection of test side, the decision statistic obtained has larger variance, causes the error rate of watermark detection high.
Summary of the invention
For prior art Problems existing, the invention provides a kind of the intersection additivity spread spectrum audio frequency watermark embedding grammar, detection method and the system that reduce watermark detection error rate.
For solving the problems of the technologies described above, the present invention adopts following technical scheme:
One, intersection additivity spread spectrum audio frequency watermark embedding grammar, comprises step:
S1 reads audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, is transformed to frequency-region signal xf by after time-domain audio signal xt framing;
S2 is according to the audio frequency range of sample rate f s1, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range;
It is the bipolarity pseudorandom spreading sequence u of (fmax1-fmin1+1) that S3 utilizes random number seed seed to generate length by rand function, fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum value;
Signal frame in frequency domain signal X F is numbered by S4 in turn, the signal frame of odd indexed is combined into odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven;
S5 is according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed;
Odd indexed signal frame YFodd and even number sequence number signal frame YFeven is combined into signal frame YF by sequence number by S6;
Signal frame YF is transformed to time-domain signal by S7, generates audio file, the i.e. audio file of embed watermark.
S3 comprises further:
RandStream function is utilized to carry out initialization to rand function and generate random number seed seed;
Random number seed seed and rand function is utilized to generate the random number sequence that length is (fmax1-fmin1+1);
Random number each in random number sequence is rounded up respectively, the binary sequence that must be made up of 0 and 1;
Binary number each in binary sequence is all taken advantage of to 2 and subtracted 1, obtains bipolarity pseudorandom spreading sequence u.
Two, intersection additivity spread spectrum audio frequency watermark detection method, comprises step:
The audio file of S1 reading embed watermark obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, is transformed into frequency-region signal sf by after time-domain audio signal rt framing;
S2 is according to the audio frequency range of sample rate f s2, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range;
It is the bipolarity pseudorandom spreading sequence u of (fmax1-fmin1+1) that S3 utilizes random number seed seed to generate length by rand function, fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum;
Signal frame in frequency-region signal SF is numbered by S4 in turn, the signal frame of odd indexed is combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven;
S5, based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively to odd indexed signal frame SFodd and even number sequence number signal frame SFeven, obtains strange decision statistic rodd and even decision statistic reven;
S6 obtains decision statistic rall=rodd-reven according to strange decision statistic rodd and even decision statistic reven, by judging the symbol extraction watermark of decision statistic rall.
S3 comprises further:
RandStream function is utilized to carry out initialization to rand function and generate random number seed seed;
Random number seed seed and rand function is utilized to generate the random number sequence that length is (fmax1-fmin1+1);
Random number each in random number sequence is rounded up respectively, the binary sequence that must be made up of 0 and 1;
Binary number each in binary sequence is all taken advantage of to 2 and subtracted 1, obtains bipolarity pseudorandom spreading sequence u.
Three, intersection additivity spread spectrum audio frequency watermark embedded system, comprising:
Time-frequency domain conversation module, is used for reading audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, be transformed to frequency-region signal xf by after time-domain audio signal xt framing;
Frequency range selects module, is used for according to the audio frequency range of sample rate f s1, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range;
Frequency expansion sequence generation module, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum value;
Divide odd-even interleaving frame module, be used for signal frame in frequency domain signal X F to number in turn, the signal frame of odd indexed be combined into odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven;
Watermark intersection merge module, be used for according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed;
Odd-even interleaving frame composite module, is used for odd indexed signal frame YFodd and even number sequence number signal frame YFeven to be combined into signal frame YF by sequence number;
Time-frequency domain inverse transform module, is used for signal frame YF to be transformed to time-domain signal, generates audio file, the i.e. audio file of embed watermark.
Four, intersection additivity spread spectrum audio frequency watermark detection system, comprising:
Time-frequency domain conversation module, the audio file being used for reading embed watermark obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, is transformed into frequency-region signal sf by after time-domain audio signal rt framing;
Frequency range selects module, is used for according to the audio frequency range of sample rate f s2, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range;
Frequency expansion sequence generation module, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum;
Divide odd-even interleaving frame module, be used for signal frame in frequency-region signal SF to number in turn, the signal frame of odd indexed be combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven;
Coherent detection module, is used for based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively, obtain strange decision statistic rodd and even decision statistic reven to odd indexed signal frame SFodd and even number sequence number signal frame SFeven;
Intersection judging module, is used for obtaining decision statistic rall=rodd-reven according to strange decision statistic rodd and even decision statistic reven, by judging the symbol extraction watermark of decision statistic rall.
Compared to the prior art, tool of the present invention has the following advantages and beneficial effect:
The present invention can reduce decision statistic variance greatly, thus reduces the error rate of watermark detection.
Accompanying drawing explanation
Fig. 1 is that the present invention intersects the structured flowchart of additivity spread spectrum audio frequency watermark embedded system;
Fig. 2 is that the present invention intersects the structured flowchart of additivity spread spectrum audio frequency watermark detection system;
Fig. 3 is that the present invention intersects the particular flow sheet of additivity spread spectrum audio frequency watermark embedding grammar;
Fig. 4 is that the present invention intersects the particular flow sheet of additivity spread spectrum audio frequency watermark detection method.
Embodiment
Below in conjunction with accompanying drawing, specific embodiments of the present invention is described in detail.
See Fig. 1, the intersection additivity spread spectrum audio frequency watermark embedded system that the present embodiment provides, comprise time-frequency domain conversation module 1, frequency range selects module 2, frequency expansion sequence generation module 3, divide odd-even interleaving frame module 4, watermark intersects merge module 5, odd-even interleaving frame composite module 6 and time-frequency domain inverse transform module 7, software firming bechnology can be adopted to realize each module during concrete enforcement.
Described time-frequency domain conversation module 1, is used for reading audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, be transformed to frequency-region signal xf by after time-domain audio signal xt framing.Frequency-region signal xf is exported to frequency range and select module 2.
Described frequency range selects module 2, is used for according to the audio frequency range of sample rate f s1, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range.The maximal value of domain samples point range of embed watermark and minimum value can export to frequency expansion sequence generation module 3, frequency domain signal X F be exported to a point odd-even interleaving frame module 4.
Described frequency expansion sequence generation module 3, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum value.Bipolarity pseudorandom spreading sequence u exports to watermark intersection merge module 5.
Described point odd-even interleaving frame module 4, is used for signal frame in frequency domain signal X F to number in turn, by the signal frame of odd indexed combination odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven.Odd indexed signal frame XFodd and even number sequence number signal frame XFeven is exported to watermark intersection merge module 5.
Described watermark intersection merge module 5, be used for according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed.Odd indexed signal frame YFodd and even number sequence number frame frequency-region signal YFeven is exported to odd-even interleaving frame composite module 6.
Described odd-even interleaving frame composite module 6, is used for odd indexed signal frame YFodd and even number sequence number signal frame YFeven to be combined into signal frame YF by sequence number.Signal frame YF is exported to time-frequency domain inverse transform module 7.
Described time-frequency domain inverse transform module 7, is used for signal frame YF to be transformed to time-domain signal, generates audio file, the i.e. audio file of embed watermark.
See Fig. 2, the intersection additivity spread spectrum audio frequency watermark embedded system that the present embodiment provides, comprise time-frequency domain conversation module 8, frequency range selection module 9, frequency expansion sequence generation module 10, point odd-even interleaving frame module 11, coherent detection module 12 and judging module 13 of intersecting, during concrete enforcement, software firming bechnology can be adopted to realize each module.
Time-frequency domain conversation module 8, frequency range select module 9, frequency expansion sequence generation module 10, the function of point odd-even interleaving frame module 11 selects module 2, frequency expansion sequence generation module 3, point odd-even interleaving frame module 4 with time-frequency domain conversion module 1, frequency range respectively.
Described time-frequency domain conversation module 8, the audio file being used for getting embed watermark obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, is transformed into frequency-region signal sf by after time-domain audio signal rt framing.Frequency-region signal xf is exported to frequency range and select module 9.
Described frequency range selects module 9, is used for according to the audio frequency range of sample rate f s2, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range.The maximal value of domain samples point range of embed watermark and minimum value can export to frequency expansion sequence generation module 10, frequency-region signal SF will be exported to a point odd-even interleaving frame module 11.
Described frequency expansion sequence generation module 10, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum.Bipolarity pseudorandom spreading sequence u exports to coherent detection module 12.
Described point odd-even interleaving frame module 11, is used for signal frame in frequency-region signal SF to number in turn, the signal frame of odd indexed is combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven.Odd indexed signal frame SFodd and even number sequence number signal frame SFeven is exported to coherent detection module 12.
Described coherent detection module 12, is used for based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively, obtain strange decision statistic rodd and even decision statistic reven to odd indexed signal frame SFodd and even number sequence number signal frame SFeven.Strange decision statistic rodd and even decision statistic reven is exported to intersection judging module 13.
Described intersection judging module 13, is used for obtaining decision statistic rall=rodd-reven according to strange decision statistic rodd and even decision statistic reven, by judging the symbol extraction watermark of decision statistic rall.
See Fig. 3, the intersection additivity spread spectrum audio frequency watermark embedding grammar that the present embodiment provides, can adopt computer software automatically to carry out, specifically comprise the following steps:
Step 1, reads audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, carry out time-frequency domain conversation, obtain frequency-region signal xf by after time-domain audio signal xt framing.
To the first framing of time-domain audio signal xt, frame length is frame.Time-frequency domain conversation can adopt Fast Fourier Transform (FFT) method (FFT), discrete cosine transform method (DCT), multiple modulation lapped transform method (MCLT) etc.
Step 2, according to the audio frequency range of sample rate f s1, frame length frame and people's ear sensitivity, obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range;
Can the routine techniques being calculated as watermark embedding technique field of domain samples point range of embed watermark, for ease of being beneficial to, provide detailed computing method below.
Adopting formula (1) and (2) to calculate respectively can the maximal value fmax1 of domain samples point range of embed watermark and minimum value fmin1:
fmax1=fix((FMAX*2.0/fs1)*frame)(1)
fmin1=fix((FMIN*2.0/fs1)*frame)(2)
In formula (1) ~ (2), FMAX, FMIN represent highest frequency and the low-limit frequency of the audio frequency range of people's ear sensitivity respectively, and fix () is bracket function.
Step 3, utilizes random number seed seed to generate bipolarity pseudorandom spreading sequence u by rand function.
This step adopts function in MATLAB to realize, and detailed process is as follows:
First, RandStream function is utilized to carry out initialization to rand function and generate random number seed seed.Then, random number seed seed and rand function is utilized to generate the random number sequence that length is (fmax1-fmin1+1).Random number due to the generation of rand function is the number between 0 ~ 1, to random number each in random number sequence through rounding up, obtain the binary sequence be made up of 0 and 1, binary number each in binary sequence is all taken advantage of to 2 and deducted 1, namely obtain bipolarity pseudorandom spreading sequence u.
Step 4, numbers in turn by signal frame in frequency domain signal X F, the signal frame of odd indexed is combined into odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven.
The detailed process of this step is as follows:
If frequency domain signal X F has 2N frequency samples point, frequency domain signal X F can be write as following form:
XF=[xf
1odd,xf
1even,xf
2odd,xf
2even,...,xf
Nodd,xf
Neven](3)
XFodd is as follows for odd indexed signal frame:
XFodd=[xf
1odd,xf
2odd,...,xf
Nodd](4)
XFeven is as follows for even number sequence number signal frame:
XFeven=[xf
1even,xf
2even,...,xf
Neven](5)
Xf
iodd, xf
ievenrepresent 2i-1 and 2i signal frame in frequency domain signal X F respectively, i=1,2 ... N.
Step 5, according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed.
The specific implementation process of this step is as follows:
Adopt formula (6) ~ (7) that watermark bit symbol and-1*symbol is embedded odd indexed signal frame XFodd and even number sequence number signal frame XFeven respectively:
YFodd=XFodd+smybol*u(6)
YFeven=XFeven-smybol*u(7)
Step 6, is combined into signal frame YF by odd indexed signal frame YFodd and even number sequence number signal frame YFeven by sequence number, sees formula (8).
YF=[yf
1odd,yf
1even,yf
2odd,yf
2even,...,yf
Nodd,yf
Neven](8)
Yf
ioddrepresent i-th signal frame in odd indexed signal frame YFodd, yf
ievenrepresent i-th signal frame in even number sequence number signal frame YFeven, i=1,2 ... N.
Step 7, carries out time-frequency domain inverse transformation to signal frame YF, obtains time-domain signal, and generate audio file, i.e. the audio file of embed watermark.
See Fig. 4, the intersection additivity spread spectrum audio frequency watermark detection method that the present embodiment provides, can adopt computer software automatically to carry out, specifically comprise the following steps:
Step 8, reads the audio file of embed watermark, obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, carry out time-frequency domain conversation by after time-domain audio signal rt framing, obtain frequency-region signal sf.
To the first framing of time-domain audio signal rt, frame length is with frame length frame in step 1.Time-frequency domain conversation method is with step 1.
Step 9, according to the audio frequency range of sample rate f s2, frame length frame and people's ear sensitivity, obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range.
Can formula (1) and (2) be adopted to calculate acquisition by the maximal value fmax1 of domain samples point range of embed watermark and minimum value fmin1.
Step 10, utilizes random number seed seed to generate bipolarity pseudorandom spreading sequence u by rand function.
This step adopts function in MATLAB to realize, and detailed process is with step 3, and gained bipolarity pseudorandom spreading sequence is identical with step 3 gained bipolarity pseudorandom spreading sequence.
Step 11, numbers in turn by signal frame in frequency-region signal SF, the signal frame of odd indexed is combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven.
This step is with step 4, and its specific implementation process repeats no more.
Step 12, based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively to odd indexed signal frame SFodd and even number sequence number signal frame SFeven, obtains corresponding strange decision statistic rodd and even decision statistic reven.
The computing formula of strange decision statistic rodd and even decision statistic reven is as follows:
In formula (9) ~ (10), < > represents coherent detection; SFodd
i, SFeven
i, u
ibe respectively i-th data point in SFodd, SFeven, u.
Step 12, utilizes strange decision statistic rodd and even decision statistic reven to obtain decision statistic rall, by judging the watermark of decision statistic rall symbol extraction.
Formula (11) is adopted to obtain decision statistic rall:
rall=rodd-reven(11)
If rall >=0, the watermark received_bit=1 so extracted, otherwise, received_bit=0.
Claims (6)
1. intersection additivity spread spectrum audio frequency watermark embedding grammar, is characterized in that, comprise step:
S1 reads audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, is transformed to frequency-region signal xf by after time-domain audio signal xt framing;
S2 is according to the audio frequency range of sample rate f s1, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range;
It is the bipolarity pseudorandom spreading sequence u of (fmax1-fmin1+1) that S3 utilizes random number seed seed to generate length by rand function, fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum value;
Signal frame in frequency domain signal X F is numbered by S4 in turn, the signal frame of odd indexed is combined into odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven;
S5 is according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed;
Odd indexed signal frame YFodd and even number sequence number signal frame YFeven is combined into signal frame YF by sequence number by S6;
Signal frame YF is transformed to time-domain signal by S7, generates audio file, the i.e. audio file of embed watermark.
2. the additivity spread spectrum audio frequency watermark embedding grammar that intersects as claimed in claim 1, is characterized in that:
S3 comprises further:
RandStream function is utilized to carry out initialization to rand function and generate random number seed seed;
Random number seed seed and rand function is utilized to generate the random number sequence that length is (fmax1-fmin1+1);
Random number each in random number sequence is rounded up respectively, the binary sequence that must be made up of 0 and 1;
Binary number each in binary sequence is all taken advantage of to 2 and subtracted 1, obtains bipolarity pseudorandom spreading sequence u.
3. intersection additivity spread spectrum audio frequency watermark detection method, is characterized in that, comprise step:
The audio file of S1 reading embed watermark obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, is transformed into frequency-region signal sf by after time-domain audio signal rt framing;
S2 is according to the audio frequency range of sample rate f s2, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range;
It is the bipolarity pseudorandom spreading sequence u of (fmax1-fmin1+1) that S3 utilizes random number seed seed to generate length by rand function, fmax1 and fmin1 is can the minimum and maximum value of domain samples point range of embed watermark;
Signal frame in frequency-region signal SF is numbered by S4 in turn, the signal frame of odd indexed is combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven;
S5, based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively to odd indexed signal frame SFodd and even number sequence number signal frame SFeven, obtains strange decision statistic rodd and even decision statistic reven;
S6 obtains decision statistic rall=rodd-reven according to strange decision statistic rodd and even decision statistic reven, by judging the symbol extraction watermark of decision statistic rall.
4. additivity spread spectrum audio frequency watermark detection method of intersecting as claimed in claim 2, is characterized in that:
S3 comprises further:
RandStream function is utilized to carry out initialization to rand function and generate random number seed seed;
Random number seed seed and rand function is utilized to generate the random number sequence that length is (fmax1-fmin1+1);
Random number each in random number sequence is rounded up respectively, the binary sequence that must be made up of 0 and 1;
Binary number each in binary sequence is all taken advantage of to 2 and subtracted 1, obtains bipolarity pseudorandom spreading sequence u.
5. intersection additivity spread spectrum audio frequency watermark embedded system, is characterized in that, comprising:
Time-frequency domain conversation module, is used for reading audio file and obtains time-domain audio signal xt and sample rate f s1 thereof, be transformed to frequency-region signal xf by after time-domain audio signal xt framing;
Frequency range selects module, is used for according to the audio frequency range of sample rate f s1, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency domain signal X F of this domain samples point range;
Frequency expansion sequence generation module, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum value;
Divide odd-even interleaving frame module, be used for signal frame in frequency domain signal X F to number in turn, the signal frame of odd indexed be combined into odd indexed signal frame XFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame XFeven;
Watermark intersection merge module, be used for according to bipolarity pseudorandom spreading sequence u, adopt additivity spread-spectrum watermark embedding inlay technique, to XFodd and XFeven embed watermark bit symbol and-1*symbol respectively, obtain the odd indexed signal frame YFodd after embed watermark and even number sequence number signal frame YFeven, symbol is the watermark bit that should embed;
Odd-even interleaving frame composite module, is used for odd indexed signal frame YFodd and even number sequence number signal frame YFeven to be combined into signal frame YF by sequence number;
Time-frequency domain inverse transform module, is used for signal frame YF to be transformed to time-domain signal, generates audio file, the i.e. audio file of embed watermark.
6. intersection additivity spread spectrum audio frequency watermark detection system, is characterized in that, comprising:
Time-frequency domain conversation module, the audio file being used for reading embed watermark obtains time-domain audio signal rt and the sample rate f s2 thereof of embed watermark, is transformed into frequency-region signal sf by after time-domain audio signal rt framing;
Frequency range selects module, is used for according to the audio frequency range of sample rate f s2, frame length and people's ear sensitivity, and obtaining can the domain samples point range of embed watermark, thus obtains the frequency-region signal SF of this domain samples point range;
Frequency expansion sequence generation module, be used for utilizing random number seed seed to generate by rand function the bipolarity pseudorandom spreading sequence u that length is (fmax1-fmin1+1), fmax1 and fmin1 is can the maximal value of domain samples point range of embed watermark and minimum;
Divide odd-even interleaving frame module, be used for signal frame in frequency-region signal SF to number in turn, the signal frame of odd indexed be combined into odd indexed signal frame SFodd, the signal frame of even number sequence number is combined into even number sequence number signal frame SFeven;
Coherent detection module, is used for based on bipolarity pseudorandom spreading sequence u, carries out coherent detection operation respectively, obtain strange decision statistic rodd and even decision statistic reven to odd indexed signal frame SFodd and even number sequence number signal frame SFeven;
Intersection judging module, is used for obtaining decision statistic rall=rodd-reven according to strange decision statistic rodd and even decision statistic reven, by judging the symbol extraction watermark of decision statistic rall.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510833111.5A CN105374360B (en) | 2015-11-25 | 2015-11-25 | Intersect additivity spread spectrum audio frequency watermark embedding grammar, detection method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510833111.5A CN105374360B (en) | 2015-11-25 | 2015-11-25 | Intersect additivity spread spectrum audio frequency watermark embedding grammar, detection method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105374360A true CN105374360A (en) | 2016-03-02 |
CN105374360B CN105374360B (en) | 2018-12-14 |
Family
ID=55376490
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510833111.5A Active CN105374360B (en) | 2015-11-25 | 2015-11-25 | Intersect additivity spread spectrum audio frequency watermark embedding grammar, detection method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105374360B (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106409302A (en) * | 2016-06-22 | 2017-02-15 | 华中师范大学 | Audio frequency watermark method and system based on embedding area selection |
CN106599630A (en) * | 2016-12-12 | 2017-04-26 | 中国矿业大学 | Improved additive spread spectrum watermarking method |
CN109036442A (en) * | 2018-07-27 | 2018-12-18 | 北京感动无限科技有限公司 | Watermark recognition methods and device for broadcasting content |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1898396A1 (en) * | 2006-09-07 | 2008-03-12 | Deutsche Thomson-Brandt Gmbh | Method and apparatus for encoding/decoding symbols carrying payload data for watermarking of an audio or video signal |
CN101185121A (en) * | 2005-06-02 | 2008-05-21 | 汤姆森许可贸易公司 | Method and apparatus for watermarking an audio or video signal with watermark data using a spread spectrum |
CN101290773A (en) * | 2008-06-13 | 2008-10-22 | 清华大学 | Adaptive MP3 digital watermark embedding method |
CN103957099A (en) * | 2014-05-16 | 2014-07-30 | 武汉大学 | Media encryption and decryption method combined with watermark |
CN104217725A (en) * | 2014-09-29 | 2014-12-17 | 北京理工大学 | Audio watermarking method based on multi-echo core |
CN104320719A (en) * | 2014-11-14 | 2015-01-28 | 武汉大学 | Television program interaction participating method and system based on audio watermarking |
CN104658542A (en) * | 2015-03-16 | 2015-05-27 | 武汉大学 | Additive spread spectrum audio watermarking embedding method, additive spread spectrum audio watermarking detection method and additive spread spectrum audio watermarking embedding system based on orthogonality |
-
2015
- 2015-11-25 CN CN201510833111.5A patent/CN105374360B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101185121A (en) * | 2005-06-02 | 2008-05-21 | 汤姆森许可贸易公司 | Method and apparatus for watermarking an audio or video signal with watermark data using a spread spectrum |
EP1898396A1 (en) * | 2006-09-07 | 2008-03-12 | Deutsche Thomson-Brandt Gmbh | Method and apparatus for encoding/decoding symbols carrying payload data for watermarking of an audio or video signal |
CN101290773A (en) * | 2008-06-13 | 2008-10-22 | 清华大学 | Adaptive MP3 digital watermark embedding method |
CN103957099A (en) * | 2014-05-16 | 2014-07-30 | 武汉大学 | Media encryption and decryption method combined with watermark |
CN104217725A (en) * | 2014-09-29 | 2014-12-17 | 北京理工大学 | Audio watermarking method based on multi-echo core |
CN104320719A (en) * | 2014-11-14 | 2015-01-28 | 武汉大学 | Television program interaction participating method and system based on audio watermarking |
CN104658542A (en) * | 2015-03-16 | 2015-05-27 | 武汉大学 | Additive spread spectrum audio watermarking embedding method, additive spread spectrum audio watermarking detection method and additive spread spectrum audio watermarking embedding system based on orthogonality |
Non-Patent Citations (1)
Title |
---|
A VALIZADEH,ZJ WANG: ""Improved multiplicative spread spectrum embedding for image data hiding"", 《IEEE INTERNATIONAL CONFERENCE ON IMAGE PROCESSING》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106409302A (en) * | 2016-06-22 | 2017-02-15 | 华中师范大学 | Audio frequency watermark method and system based on embedding area selection |
CN106409302B (en) * | 2016-06-22 | 2019-07-09 | 华中师范大学 | Audio-frequency water mark method and system based on insertion regional choice |
CN106599630A (en) * | 2016-12-12 | 2017-04-26 | 中国矿业大学 | Improved additive spread spectrum watermarking method |
CN109036442A (en) * | 2018-07-27 | 2018-12-18 | 北京感动无限科技有限公司 | Watermark recognition methods and device for broadcasting content |
Also Published As
Publication number | Publication date |
---|---|
CN105374360B (en) | 2018-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104658542A (en) | Additive spread spectrum audio watermarking embedding method, additive spread spectrum audio watermarking detection method and additive spread spectrum audio watermarking embedding system based on orthogonality | |
CN105976823A (en) | Adaptive audio watermarking method based on phase coding and system | |
CN100559466C (en) | A kind of audio-frequency watermark processing method of anti-DA/AD conversion | |
CN105374360A (en) | Interleaved additive spread spectrum audio watermark embedding method and detection method and system | |
CN110163787B (en) | Audio digital robust blind watermark embedding method based on dual-tree complex wavelet transform | |
CN103455972A (en) | Multiple-bit spread spectrum watermark method with robustness enhanced by utilization of carrier | |
Hua et al. | Cepstral analysis for the application of echo-based audio watermark detection | |
CN103413552A (en) | Audio watermark embedding and extracting method and device | |
CN104795071A (en) | Blind audio watermark embedding and watermark extraction processing method | |
CN113782041B (en) | Method for embedding and positioning watermark based on audio variable frequency domain | |
Maha et al. | DCT based blind audio watermarking scheme | |
CN101494053A (en) | Method and apparatus for determining whether or not a reference pattern is present in a received and possibly water-marked signal | |
CN103050120A (en) | High-capacity digital audio reversible watermark processing method | |
CN111641422A (en) | System and method for determining self-adaptive anti-interference detection threshold of high-dynamic digital receiver | |
CN102867513B (en) | Pseudo-Zernike moment based voice content authentication method | |
CN102509551A (en) | Audio steganalysis algorithm based on variance of difference of sum cepstrum | |
Li et al. | Design of a general hardware in the loop underwater communication emulation system | |
Wu et al. | Audio watermarking algorithm with a synchronization mechanism based on spectrum distribution | |
CN107689815A (en) | A kind of MSK band spread receivers Adaptive Anti-jamming detection threshold determines method | |
CN106409302B (en) | Audio-frequency water mark method and system based on insertion regional choice | |
CN104217725A (en) | Audio watermarking method based on multi-echo core | |
CN106205627B (en) | Digital audio reversible water mark algorithm based on side information prediction and histogram translation | |
Guo et al. | An audio digital watermarking algorithm against A/D and D/A conversions based on DCT domain | |
US9742554B2 (en) | Systems and methods for detecting a synchronization code word | |
CN102710564A (en) | Channel time domain impulse response filter method and device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |