%\documentstyle[twocolumn]{article}
\documentclass[twocolumn]{article}
\usepackage{amsmath,epsfig}
%\usepackage{mdwlist}
%\usepackage[small,compact]{titlesec}
%\usepackage[small,it]{caption}
\usepackage{picins}


\pagestyle{empty}
\setlength{\textwidth}{7in}
\setlength{\textheight}{9.125in}
\setlength{\columnsep}{0.5in}
\setlength{\topmargin}{-0.8in}
\setlength{\oddsidemargin}{-0.25in}
\setlength{\parindent}{5 ex}

\makeatletter 
\def\@normalsize{\@setsize\normalsize{12pt}\xpt\@xpt \abovedisplayskip 11pt 
plus2pt minus5pt\belowdisplayskip \abovedisplayskip \abovedisplayshortskip \z@ 
plus3pt\belowdisplayshortskip 6pt plus3pt 
minus3pt\let\@listi\@listI} 
%the following line was changed 
\def\subsize{\@setsize\subsize{12pt}\xipt\@xipt} 
\def\section{\@startsection{section}{1}{\z@}{24pt plus 2 pt 
minus 2 pt} {12pt plus 2pt minus 2pt}{\large\bf}} 
\def\subsection{\@startsection {subsection}{2}{\z@}{12pt 
plus 2pt minus 2pt}{12pt plus 2pt minus 2pt}{\subsize\bf}} 
%\makeatother 

%\makeatletter
%\newcommand\std@startsection{}
%\let\std@startsection\@startsection
%\def\@startsection{\renewcommand*\baselinestretch{1}\std@startsection}
%\makeatother

\begin{document} 
\date{} 
%\title{\Large\bf Voice-Based DRM System For Accces Control } 
\title{\Large\bf Watermarking of Space Curves using Wavelet Decomposition} 
\author{ R. C. Motwani, M. C. Motwani, and F. C. Harris, Jr. \\ 
 Computer Science and Engineering Department \\ 
 University of Nevada, Reno \\ 
 Reno, NV ~~89557  
} 
\maketitle 
\thispagestyle{empty} 
\subsection*{\centering Abstract} 
\vspace*{-3mm} 
This paper describes an imperceptible, non-blind, fragile watermarking technique for space curves. The proposed technique employs a wavelet-based approach, and computes a multi-resolution representation of the space curve to embed a watermark so that it has widespread presence in the curve. Watermarks that are widely spread within the host data can not be easily damaged by cropping and replacement attacks that result in localized alteration of the host data.   A variety of wavelet families are exploited and experimental results provide a comparison of the performance of different wavelets in terms of the watermark's imperceptibility and tolerance to attacks. To quantify space curve distortion, a signal-to-noise ratio is used, and a linear correlation measure is employed to determine the resistance of the watermark to modifications. Results indicate that watermark insertion using wavelet packet decomposition outperforms orthogonal wavelet decomposition in terms of imperceptibility and diversified presence of the embedded watermark.


\section{Introduction}
Motion capture (MoCap) technology yields appealing computer graphics animations but entails high investments in terms of cost, time and effort. The digital nature of MoCap data makes it vulnerable to piracy and plagiarism, thereby discouraging MoCap studios and labs from publishing such data. This paper focuses on tamper detection in trajectories (space curves) derived from motion capture data, to assist in detecting modifications that violate copyrights of motion data extracted from published MoCap datasets. 

\begin{figure}[t!]
	\centering
		\includegraphics[width=0.35\textwidth]{mocap.png}
	\caption{Trajectory Plot(red) of Left Hand Joint of Human Skeleton(blue). Motion Sequence from \textit{climb.bvh} }
	\label{fig:3dcurve}
\end{figure}

Watermarking techniques have been used for copyright protection, ownership authentication, and tamper proofing of digital data. Watermarking schemes insert information in the digital content in such a way that the embedded information is imperceptible to the the human eye. Robust watermarking techniques strive to embed information in such a way that it is difficult to remove without causing perceivable distortions to the original data. However, this is a challenging research problem therefore such schemes are only tolerant to a limited set of attacks. Fragile watermarking schemes, on the other hand, embed watermarks that have low resistance to modifications and are destroyed at the slightest variation to the host content. Therefore, fragile schemes find applications in tamper proofing digital data, since a damaged watermark is an indication of a malicious modification attempt to the data.  

Research related to watermarking of 3D data is still in its infancy, and finds applications to 3D meshes and motion data streams. The work presented here explores the use of various wavelet decomposition techniques for watermark insertion and extraction. The contribution of this paper lies in analyzing the impact of embedding random noise as watermark in different wavelet subbands and wavelet coefficients obtained using wavelet packets. Elaborate experimentation has yielded results that indicate improved performance of wavelet packet technique over the rest. 

The remainder of this paper is organized as follows. Section 2 presents the related work in this relatively immature field. Section 3 describes the proposed watermarking approach. Section 4 provides the results of experiments. Conclusions with future work are summed up in Section 5.
\begin{figure*}[t!]
	\centering
		\includegraphics[width=0.92\textwidth]{wav_x.png}
	\caption{Original Signal in $x$-dimension and its Multiresolution Wavelet Decomposition at Levels 1 through 5  }
	\label{fig:wav_x}
\end{figure*}

\section{Related Work}
Related work on curve watermarking has been investigated for planar curves (in the 2D contex) for copyright protection of digitally distributed maps (~\cite{wm2dcurve},\cite{curvewm},\cite{wmline}), vector fonts \cite{wmkfont}, hand drawn curves and topographic maps \cite{fingprintmaps}. However, limited work has been done on curve/trajectory/motion-data watermarking in the 3D domain.  



The authors in \cite{13} propose a progressive watermarking scheme for 3D motion capture data that uses frame decimation.  A robust, blind 3D motion capture data watermarking algorithm for human motion animation is proposed in \cite{15}, that is cluster-based and uses quantization techniques. The authors in \cite{14}, describe a spatial domain technique to watermark 3D motion capture data. Pu \textit{et al.} \cite{wmocap1}, adopt singular value decomposition to consider both the time varying relations among the motion frames and the spatial correlations  among the different joints in motion. The motion data matrix is decomposed into two eigen vector matrices and a singular values matrix. The watermark is added to the singular values matrix. Agarwal and Prabhakaran \cite{wmocap2} propose a tamper-proofing mechanism for MoCap data that applies hash functions to the data matrix and embed identifiers as watermarks to detect attacks such as row/column shuffling and element shuffling.

Most watermarking techniques \cite{4} adopt a certain level of randomness in the algorithm to battle attacks on watermark removal by brute force approach. However, this is the simplest approach and has its drawbacks. Embedding the watermark directly in the spatial domain makes it vulnerable to removal or replacement attacks. It is preferred to transform the motion data into frequency domain. This assures that the watermark is spread across the 3D curve such that removal or replacement of parts of the curve does not destroy the watermark completely. 
In \cite{6}, Yamazaki proposes segmentation of the motion data followed by a discrete cosine transform operation on each segment to embed the watermark in the spread spectrum domain. In \cite{yamazaki}, Yamazaki employs wavelet-based spectral analysis for watermark insertion. \begin{figure*}[t!]
	\centering
		\includegraphics[width=1.0\textwidth]{blockdiagram.png}
	\caption{Watermark Insertion and Extraction Process using Wavelet Decomposition}
	\label{fig:system_block}
\end{figure*} 

The watermarking approach presented in this paper also employs wavelets but differs from Yamazaki's approach as it utilizes a multiresolution wavelet representation of the 3D curve for watermark insertion. Moreover, the proposed approach isolates the trajectories of the human skeletal joints and applies to the space curve generated by each joint. In addition, a variety of wavelet families are experimented with to determine the best performer. 



\section{Methodology}
For MoCap data, a space curve is a sequence of coordinates in 3D space. This space curve is derived from the motion of one joint (denoted by dot marker) of the human skeleton, as shown in Fig. \ref{fig:3dcurve}. The MoCap dataset used for this figure is obtained from \textit{BeyondMotion Studio} \cite{modata} and represents the \textit{climb.bvh} sequence. The space curve is the trajectory represented by red markers in the plot shown in  Fig. \ref{fig:3dcurve}. 
 
\begin{figure}[h!]
\centering
\begin{tabular}{cc}
\includegraphics[width=0.23\textwidth]{spacecurve.png} &
\includegraphics[width=0.23\textwidth]{reconstruct_level2.png} \\
\includegraphics[width=0.23\textwidth]{reconstruct_level3.png} &
\includegraphics[width=0.23\textwidth]{reconstruct_level4.png} \\
\end{tabular}
\caption{Multiresolution Analysis of the Space Curve- The space curve is represented at decreasing scales C[2],C[3] and C[4] (Level-2 wavelet transform yields a higher scale, Level-4 results in a lower scale). At lower resolutions the finer details are lost during reconstruction.}
	\label{fig:mra}
\end{figure}  


The proposed approach transforms the spatial representation of the 3D curve to the spread spectrum domain using wavelets \cite{11}. Wavelet transform is preferred over Fourier or Discrete Cosine transforms because it captures both the global pattern (i.e. averages or approximations) and the local variations (i.e. fluctuations or details) in the curve. Wavelet functions decompose a space curve into multiple resolutions thereby facilitating examination of the gross and finer details of the curve at different scales or resolutions (see Fig. \ref{fig:wav_x}). The following subsections describe the  proposed multiresolution wavelet-based, non-blind watermarking approach.



\subsection{MultiResolution Analysis}
The space curve is represented by a three-dimensional discrete signal $C$ of length $n$. The wavelet transform is applied to the $x$, $y$ and $z$ co-ordinates of $C$ separately. As depicted by Eq. \ref{eq:1}, a discrete wavelet transform (DWT) applied to $C$ decomposes the signal into two sub-signals, $S_{i}$ and $W_{i}$, each of half its length 
\begin{figure*}[t!]
	\centering
		\includegraphics[width=0.92\textwidth]{wavpckt_x.png}
	\caption{Original Signal in $x$-dimension and its Wavelet Packet Decomposition at Levels 1 through 3 }
	\label{fig:wavpckt_x}
\end{figure*} 
($m=\frac{n}{2}$ where $n$ is an integral power of 2 with zero padding), where $i$ represents the multiresolution level of wavelet transform. 



\begin{equation}
	DWT(C[n]) = S_{i}[m]+W_{i}[m]
	\label{eq:1}
\end{equation}


The first sub-signal constitutes the scalar co-efficients that represent the approximation of the original signal and is computed by the following equation:             
\begin{equation}
	S_{i} = \sum_{k}{C(k)\phi_{i}(k)}  
	\label{eq:2}
\end{equation}
where $\phi(k)$ represents the scaling function of the chosen wavelet family.
   

                     
The second sub-signal represents the wavelet co-efficients that constitute the differences  between the subsequent components of the original signal and is denoted by:

\begin{equation}
 W_{i} =  \sum_{k}{C(k)\varphi_{i}(k)}
	\label{eq:3}
\end{equation}
where $\varphi(k)$ represents the wavelet function of the chosen wavelet family.  



The functions $\phi(k)$ and $\varphi(k)$ are defined by the chosen wavelet. \textit{Haar, Daubechies, Biorthogonal, Meyer,} and \textit{Mexican Hat} are different  families of wavelets. Readers are advised to refer to \cite{mallat}, \cite{multiresolutioncurves}, and \cite{dwt} for further details on wavelet transform.   

At level $i=1$, the Level-1 resolution space curve $C[1] = S_{1}+W_{1}$. The Level-2 resolution space curve is obtained by applying DWT only on the approximation coefficients $S_{1}$ and Level-n wavelet transform of $C$ is obtained by $DWT(S_{n-1})$. 

A multi-resolution representation of the space curve decomposed at levels with decreasing resolution, is demonstrated in Fig. \ref{fig:mra}. A visual representation of multiresolution wavelet decomposition of $C$ into approximation and detail co-efficients, in the $x$-dimension, is shown in Fig. \ref{fig:wav_x}. In this figure, the scalar and wavelet coefficients at Level 2,3,4, and 5 are obtained by taking the DWT on the scalar coefficients of the previous level. 

\begin{figure}[h!]
	\centering
		\includegraphics[width=0.48\textwidth]{wpdec.png}
	\caption{Original Signal (0,0) and its Tree representation of Wavelet Packet Decomposition at Levels 1 through 3. (1,0) represents the scalar coefficients at Level-1. (1,1) represents the wavelet coefficients at Level-1. For Level-$N$, scalar-coefficients are denoted by ($N$, even index) and wavelet coefficients are denoted by ($N$,odd index).}
	\label{fig:tree-wp}
\end{figure}


\subsection{Wavelet Packets Analysis}
The wavelet packet decomposition provides a richer signal analysis. In the orthogonal wavelet decomposition which is described in the previous sub-section, the approximation coefficients are decomposed into two parts, i.e. scalar and wavelet coefficients, at each level of DWT. The wavelet coefficients are not further decomposed at any levels. In the wavelet packets technique, the wavelet coefficients are also decomposed \begin{figure*}[t!]
	\centering
		\includegraphics[width=1.0\textwidth]{blockdiagram-wp.png}
	\caption{Watermark Insertion and Extraction Process using Wavelet Packets}
	\label{fig:system_block-wp}
\end{figure*} at each level into scalar and wavelet coefficients by applying the DWT. The complete binary tree produced by the wavelet packet transform is illustrated in Fig. \ref{fig:tree-wp}. This results in a much detailed analysis of the space curve, as demonstrated by Fig. \ref{fig:wavpckt_x}. 

                                 
\subsection{Watermark Embedding}
The steps underlying the process of watermark insertion and extraction are demonstrated by 
Fig.~\ref{fig:system_block} and \ref{fig:system_block-wp}. The watermark insertion process adds a random watermark $R_{j}$ to the multiresolution wavelet coefficients $W_{j}$ or Level-$N$ wavelet packet coefficients (represented by ($N$,1) through ($N$,$2^{N}-1$), as indicated by the wavelet packet tree in Fig. \ref{fig:tree-wp}.) selected by a key $K_{j}$, which is derived from a pseudo-random number generator function,  where $j$ represents the $x,y,z$ dimension. The watermark $R$ is a sequence of pseudo-random numbers. The watermark is multiplied by a scaling factor $M$, which determines the embedding strength. Experimental values for $M$ lie in the range  $10^{-4}$ to $10^{-5}$.


The watermark is inserted into the multiresolution wavelet or wavelet packet coefficients according to the following equation: 
\begin{equation}
	W^{'}_{j}(k)= W_{j}(k) + R_{j}M
\end{equation}
where $W^{'}$ denotes the watermarked wavelet coefficient, $k$ denotes the wavelet coefficient's index selected by key $K_{j}$, and $j$ represents the $x,y$ and $z$ coordinates of the space curve. 

Inverse transform applied to the unmodified scalar coefficients and the modified wavelet coefficients yields the watermarked space curve as shown in Fig. \ref{fig:wm}. The space curve in this figure represents the trajectory generated by red markers plotted in Fig. \ref{fig:3dcurve}, but it looks different since it has been plotted independently of the skeleton with the x,y,z axes swapped and does not incorporate the scaling of the coordinate axes in the plot. 

\begin{figure}[h!]
	\centering
		\includegraphics[width=0.5\textwidth]{bior2.png}
		\caption{Original(blue) and	Watermarked(green) Space Curves}
	\label{fig:wm}
\end{figure}


\subsection{Watermark Detection}
To detect if a space curve has been modified, wavelet domain representation of the original 3D curve is subtracted from the wavelet domain representation of the test space curve. The extraction process requires the key $K$, hence the watermarking technique is non-blind.   Correlation of the subtraction result with the original watermark determines if the curve has been tampered with or not.

A linear correlation coefficient $corr$ is used as the metric for similarity between the original and 
\begin{figure*}[t!]
	\centering
	\begin{tabular}{cccc}
		\includegraphics[width=0.23\textwidth]{haar1.png} &
		\includegraphics[width=0.23\textwidth]{db3.png} &
		\includegraphics[width=0.23\textwidth]{bior4.png} &
		\includegraphics[width=0.23\textwidth]{rbio5.png} \\
		\includegraphics[width=0.23\textwidth]{db1.png} &
		\includegraphics[width=0.23\textwidth]{coif2.png} &
		\includegraphics[width=0.23\textwidth]{dmey3.png} &
		\includegraphics[width=0.23\textwidth]{sym4.png} \\
	\end{tabular}
	\caption{Distortion Analysis-Original Space Curve(blue) and Watermarked Space Curve(green) at Different Levels of Transform for Different Wavelet Families }
	\label{fig:distortion}
\end{figure*} extracted watermark. Given pairs of quantities (i.e. two sets of data $A$ and $B$) $(A_{j}, B_{j})$, where $j = 1,\ldots, N$ and $\bar{A}$ is the mean of all $A_{j}$'s and $\bar{B}$ is the mean of all $B_{j}$'s, $corr$ is given by the formula:
\begin{equation}
  corr =  \frac{\sum_{j}  \left(A_{j} - \bar{A} \right) \left(B_{j} - \bar{B} \right)}{\sqrt{\sum_{j}  \left(A_{j} - \bar{A} \right)^{2}} \sqrt{\sum_{j}  \left(B_{j} - \bar{B} \right)^{2}}}
    \label{eq:5}
\end{equation}

When $corr=1$, the extracted watermark is identical to the original watermark, which implies that the test curve has not been tampered with.  



\section{Experiments}
The experiments were done in Matlab using the Wavelet toolbox and Motion Capture toolbox \cite{mocaptoolbox}. The data used in this project is obtained from \cite{cmumodata}. Distortion analysis of the original and watermarked space curves is based on the signal-to-noise ratio (SNR) metric which is given by the following equation:

\begin{equation}
	SNR(C, C^{'}) = 20Log_{10} \left( \frac{RMS(C)}{RMS(C-C^{'})}\right)
\end{equation}
where $C$ is the original space curve and $C^{'}$ represents the watermarked space curve. $RMS$ denotes the root-mean-square value. The imperceptibility of the watermarking algorithm is measured by this SNR value. 

Results for distortion analysis for the space curve, 
shown in Fig. \ref{fig:distortion} (defined by 352 points in 3D), are listed in Table \ref{tab:distortion}. Experiments are conducted on a seven families of wavelets to determine the best performers. The payload value in Table \ref{tab:distortion} represents the length of the watermark (i.e. the number of wavelet coefficients that are modified to accommodate the watermark). The payload capacity increases as the level of wavelet transform increases since the watermark is inserted into the wavelet coefficients from all levels 1 through $N$, where $N$ is the level of applied wavelet transform. For example, SNR at Level-3 for \textit{Haar} wavelet indicates presence of watermark in all Levels 1,2 and 3. Thus, SNR in Table \ref{tab:distortion} decreases as number of levels of the wavelet transform increases, since noise(watermark) is added at more levels. As depicted by Fig. \ref{fig:distortion}, a visual distortion is observed in the watermarked space curve for SNR values lower than 50. 

\begin{table}[b!]
	\centering
		\begin{tabular}{l|c|c|c|c}
		\hline
		\hline
    Payload              &113  &160 &194 & 218   \\ 
		\hline
		\hline
		Wavelet &Level-1	 &Level-2	 &Level-3  &Level-4  \\
		Family  & SNR      &SNR  &SNR  &SNR   \\
		\hline  
		\hline 
     \textit{Haar}                &50.12 	&42.02 &35.71 & 30.25   \\
    \hline
    \textit{Daubechies}	&70.24 		&59.82 &49.08 &38.52  \\
    \hline
    \textit{Biorthogonal}	   &66.03   &56.88  &48.99 &36.61  \\
    \hline
    \textit{Rev. Bior.} &59.67   &48.91 &38.63 &29.88   \\
    \hline
		\textit{Coiflets}	 &63.85 		&53.98 &45.37 &34.25   \\
    \hline
    \textit{Symlets}	&64.24 		&53.07 &44.67 &33.99   \\
    \hline
    \textit{Meyer}	 &69.48 		&61.89 &51.25 &39.61  \\
    \hline
    \hline
		\end{tabular}
		\caption{Imperceptibility Measure and Payload Capacity of the watermarking algorithm at different levels of Multiresolution Wavelet Transform for a space curve comprised of 352 points}
	\label{tab:distortion}
\end{table}



Results for various attacks on the watermarked space curve using multiresolution wavelet analysis are shown in Fig. \ref{fig:attacks} and outlined in Table \ref{tab:attacks}. The correlation measure $corr$ determines the performance of the algorithm under the following attacks:  i) cropping - in this attack parts of the space curve are removed by an adversary, ii) replacement - this attack involves modification of sections of the space curve by different data, and iii) concatenation - this attack appends data from different space curves to yield a new space curve. 

Since the proposed watermarking scheme is fragile, the watermark is destroyed at the slightest variation to the space curve caused by attacks. When $corr$ is not equal to 1, it signals a modification to the watermarked space curve thereby indicating violation of copyrights. A $corr$ value of 1 indicates proof of ownership. 

\begin{table}[h!]
	\centering
		\begin{tabular}{l|c|c|c}
		\hline
		Wavelet              &Crop	 &Replace	 &Concat  \\
		Family   &      & &  \\
		\hline
		\hline
    \textit{Haar}	        &0.3781 &0.3129 &1.0000   \\ 
    \hline
    \textit{Daubechies}	  &0.4273 &0.6741 &1.0000   \\
    \hline
    \textit{Biorthogonal}	 &0.3618 &0.5392 &1.0000   \\
    \hline
    \textit{Reverse} 	     &0.2346 &0.4075 &1.0000   \\
    \textit{Biorthogonal}  &&&\\
    \hline
		\textit{Coiflets}	   &0.3351 &0.5813 & 1.0000  \\
    \hline
    \textit{Symlets}	    &0.3468 &0.4927 & 1.0000 \\
    \hline
    \textit{Meyer}	      &0.3687 &0.5360 & 1.0000   \\
    \hline
    \hline
		\end{tabular}
	\caption{Correlation Measure for Attacks on Watermarked Space Curves obtained by MultiResolution Wavelet Analysis}
	\label{tab:attacks}
\end{table}

\begin{figure*}[t!]
	\centering
	\begin{tabular}{cccc}
		\includegraphics[width=0.23\textwidth]{wp5haar.png} &
		\includegraphics[width=0.23\textwidth]{wp4db.png} &
		\includegraphics[width=0.23\textwidth]{wp5dmey.png} &
		\includegraphics[width=0.23\textwidth]{wp2rbior.png} \\
		\includegraphics[width=0.23\textwidth]{wp3sym.png} &
		\includegraphics[width=0.23\textwidth]{wp3coif.png} &
		\includegraphics[width=0.23\textwidth]{wp1dmey.png} &
		\includegraphics[width=0.23\textwidth]{wp1sym2.png} \\
	\end{tabular}
	\caption{Distortion Analysis-Original Space Curve(blue) and Watermarked Space Curve(green) at Different Levels of Wavelet Packet Transform for Different Wavelet Families }
	\label{fig:wpdistortion}
\end{figure*} 

Results for distortion analysis for the space curve obatined by applying wavelet packets technique, shown in Fig. \ref{fig:wpdistortion} (defined by 352 points in 3D), are listed in Table \ref{tab:wpdistortion}. Experiments are conducted on a seven families of wavelets to determine the best performers. The payload value in Table \ref{tab:wpdistortion} represents the length of the watermark (i.e. the number of wavelet coefficients that are modified to accommodate the watermark). The payload capacity increases as the level of wavelet transform increases since the watermark is inserted into higher number of wavelet coefficients. SNR values in Table \ref{tab:wpdistortion} decreases as number of levels of the wavelet transform increases, since noise(watermark) is added at more levels. As depicted by Fig. \ref{fig:wpdistortion}, no visual distortion is observed even at Level-1 of wavelet transform. It has been observed that the $Meyer$ wavelet family seems to perform poorly as compared to the rest.  



Results for various attacks on the watermarked space curve using multiresolution wavelet analysis are shown in Fig. \ref{fig:attacks} and outlined in Table \ref{tab:wpattacks}. The correlation measure $corr$ determines the performance of the algorithm under the following attacks:  i) cropping - in this attack parts of the space curve are removed by an adversary, ii) replacement - this attack involves modification of sections of the space curve by different data, and iii) concatenation - this attack appends data from different space curves to yield a new space curve. 

\begin{table}[h!]
	\centering
		\begin{tabular}{l|c|c|c|c}
		\hline
		\hline
    Payload              & 139  &160 &200 & 218  \\ 
		\hline
		\hline
		Wavelet  &Level-1	 &Level-2	 &Level-3  &Level-4   \\
	 Family  & SNR     &SNR &SNR &SNR  \\
		\hline
		\hline
    \textit{Haar}	 &209.86 	&207.34 &206.90 & 207.38   \\
    \hline
    \textit{Daubechies}	&210.58		&207.86 &207.30 &207.02   \\
    \hline
    \textit{Biorthogonal} &209.64   &205.79  &205.38 &206.40  \\
    \hline
    \textit{Rev. Bior.} &210.50   & 209.03 &208.58 &207.90  \\
    \hline
		\textit{Coiflets}	&209.98 		&208.38 &207.87 &207.72  \\
    \hline
    \textit{Symlets} &210.95 		&208.30 &206.82 &207.02  \\
    \hline
    \textit{Meyer}	  &130.01 		&122.48 & 117.09 &113.95  \\
    \hline
    \hline
		\end{tabular}
		\caption{Imperceptibility Measure and Payload Capacity of the watermarking algorithm at different levels of Wavelet Packet Transform for a space curve comprised of 352 points}
	\label{tab:wpdistortion}
\end{table}
\begin{figure}[h!]
\centering
	\begin{tabular}{cc}
		\includegraphics[width=0.23\textwidth]{crop.png} &
		\includegraphics[width=0.23\textwidth]{replace.png} \\
		\multicolumn{2}{c}{\includegraphics[width=0.23\textwidth]{concat.png}}\\
	\end{tabular}
\caption{Attacks(magenta) on Watermarked Space Curve(green): Cropping(Top Left), Replacement(Top Right), and Concatenation(Bottom)}
	\label{fig:attacks}
\end{figure}

Since the proposed watermarking scheme is fragile, the watermark is destroyed at the slightest variation to the space curve caused by attacks. When $corr$ is not equal to 1, it signals a modification to the watermarked space curve thereby indicating violation of copyrights. A $corr$ value of 1 indicates proof of ownership. 

\begin{table}[h!]
	\centering
		\begin{tabular}{l|c|c|c}
		\hline
		Wavelet              &Crop	 &Replace	 &Concat  \\
		Family   &      & &  \\
		\hline
		\hline
    \textit{Haar}	        &0.6134 &0.6254 &1.0000   \\ 
    \hline
    \textit{Daubechies}	  &0.7932 &0.7216 &1.0000   \\
    \hline
    \textit{Biorthogonal}	 &0.6418 &0.7872 &1.0000   \\
    \hline
    \textit{Reverse} 	     &0.6652 &0.6523 &1.0000   \\
    \textit{Biorthogonal}  &&&\\
    \hline
		\textit{Coiflets}	   &0.7539 &0.8047 & 1.0000  \\
    \hline
    \textit{Symlets}	    &0.8317 &0.7891 & 1.0000 \\
    \hline
    \textit{Meyer}	      &0.3814 &0.4024 & 1.0000   \\
    \hline
    \hline
		\end{tabular}
	\caption{Correlation Measure for Attacks on Watermarked Space Curves obtained by Wavelet Packets }
	\label{tab:wpattacks}
\end{table}


\section{Conclusions And Future work}
This paper presents an imperceptible, fragile, non-blind watermarking technique
for space curves derived from motion capture data. The proposed watermarking algorithm is based on multiresolution wavelet analysis of the space curve. The implementation embeds information into the wavelet coefficients to minimize perceivable distortion to the space curve since the human eye can not perceive changes in the higher frequencies.  The algorithm maximizes the presence of the watermark across the entire space curve by modifying the wavelet coefficients at multiple resolution levels. The performance of various wavelet families at different levels of transform has been evaluated and experimental results indicate that the \textit{Daubechies}, \textit{Biorthogonal}, and \textit{Meyer}  wavelets yield better SNR and provide optimal performance at Level-3. Space curves with sharp discontinuities can be efficiently represented with \textit{Haar} wavelet. Motion curves do not exhibit such abruptness and therefore the experiments have demonstrated improved performance using smoother wavelets. Future work entails varying the scaling factor $M$ in accordance with the transform level of the wavelet coefficients.

Benefits of using wavelet packets is higher diversified presence and higher payload insertion capacity of the watermark as evident by the higher SNR. Results clearly indicate that wavelet packet decomposition outperforms the orthogonal wavelet decomposition technique. The improved performance is owed to the difference in the signal decomposition outlined in Fig. \ref{fig:wav_x} and Fig. \ref{fig:wavpckt_x}. Wavelet packet decomposition yields higher number of wavelet coefficients that serve as excellent hosts for the watermark signal due to their high frequency content. The higher the number of modified coefficients, the wider is the spread of the watermark. 

Watermarking of space curves can only provide copyright protection for MoCap data. Protecting copyright ownership of the skinned mesh animations generated from MoCap data is a different area of research all together, since skinned mesh animations are generated by interpolation of keyframes. Authors in \cite{last} have suggested a technique to watermark skinned mesh animations by randomly inserting watermark in mesh skin weights.

The work presented here is preliminary and focuses only on the space curve generated by one joint of the human skeleton used for MoCap animation. Further work is required to incorporate the motion constraints (temporal and spatial) while modifying the space curves of all joints of the skeleton. Future work also involves refining the algorithm such that it is resistant to various motion editing tasks \cite{motionediting} such as motion enhancement/attenuation, blending, stitching, shuffling, and noise removal.

\section*{Acknowledgment}
The data used in this project was obtained from mocap.cs.cmu.edu.
The database was created with funding from NSF EIA-0196217.

% Generated by IEEEtran.bst, version: 1.13 (2008/09/30)
\begin{thebibliography}{10}
\providecommand{\url}[1]{#1}
\csname url@samestyle\endcsname
\providecommand{\newblock}{\relax}
\providecommand{\bibinfo}[2]{#2}
\providecommand{\BIBentrySTDinterwordspacing}{\spaceskip=0pt\relax}
\providecommand{\BIBentryALTinterwordstretchfactor}{4}
\providecommand{\BIBentryALTinterwordspacing}{\spaceskip=\fontdimen2\font plus
\BIBentryALTinterwordstretchfactor\fontdimen3\font minus
  \fontdimen4\font\relax}
\providecommand{\BIBforeignlanguage}[2]{{%
\expandafter\ifx\csname l@#1\endcsname\relax
\typeout{** WARNING: IEEEtran.bst: No hyphenation pattern has been}%
\typeout{** loaded for the language `#1'. Using the pattern for}%
\typeout{** the default language instead.}%
\else
\language=\csname l@#1\endcsname
\fi
#2}}
\providecommand{\BIBdecl}{\relax}
\BIBdecl

\bibitem{wm2dcurve}
J.~Kim, D.~Im, H.~Lee, and H.~Lee, ``Watermarking curves using 2d mesh spectral
  transform,'' \emph{IEEE International Symposium on Circuits and Systems}, pp.
  2969 --2972, 2008.

\bibitem{curvewm}
H.~Yongjian, L.~Heung-Kyu, and Z.~Huafei, ``Curve watermarking technique for
  fingerprinting digital maps,'' in \emph{Proceedings of the International
  Conference on Intelligent Information Hiding and Multimedia Signal
  Processing}, 2008, pp. 223--226.

\bibitem{wmline}
V.~Solachidis and I.~Pitas, ``Watermarking polygonal lines using fourier
  descriptors,'' \emph{IEEE Computer Graphics Applications}, vol.~24, no.~3,
  pp. 44--51, 2004.

\bibitem{wmkfont}
S.~Thiemert, M.~Steinebach, and P.~Wolf, ``A digital watermark for vector-based
  fonts,'' in \emph{ACM Proceedings of the 8th workshop on Multimedia and
  security}, 2006, pp. 120--123.

\bibitem{fingprintmaps}
H.~Gou and M.~Wu, ``Data hiding in curves with application to fingerprinting
  maps,'' \emph{IEEE Transactions on Signal Processing}, vol.~53, no.~10, pp.
  3988 -- 4005, 2005.

\bibitem{13}
S.~Li and M.~Okuda, ``Iterative frame decimation and watermarking for human
  motion animation,'' \emph{ICGST International Journal on Graphics, Vision and
  Image Processing,Special Issue on Watermarking}, vol.~07, pp. 41--48, 2010.

\bibitem{15}
P.~Agarwal and B.~Prabhakaran, ``Blind robust watermarking of 3d motion data,''
  \emph{ACM Transactions on Multimedia Computing and Communication
  Applications}, vol.~6, no.~1, pp. 1--32, 2010.

\bibitem{14}
P.~Agarwal, K.~Adi, and B.~Prabhakaran, ``Robust blind watermarking mechanism
  for motion data streams,'' in \emph{Proceedings of the 8th workshop on
  Multimedia and security}.\hskip 1em plus 0.5em minus 0.4em\relax ACM, 2006,
  pp. 230--235.

\bibitem{wmocap1}
Y.~Pu, J.~Lin, M.~Hung, B.~Chen, and I.~Jou, ``Digital watermarking of mocap
  data based on singular value decomposition,'' \emph{IEEE International
  Symposium on Knowledge Acquisition and Modeling Workshop}, pp. 932 --935,
  2008.

\bibitem{wmocap2}
P.~Agarwal and B.~Prabhakaran, ``Tamper proofing mechanisms for motion capture
  data,'' in \emph{Proceedings of the 10th ACM workshop on Multimedia and
  security}, 2008, pp. 91--100.

\bibitem{4}
S.~Li and M.~Okuda, ``Watermarking for progressive human motion animation,'' in
  \emph{IEEE International Conference on Multimedia and Expo}, 2007, pp. 1259
  -- 1262.

\bibitem{6}
S.~Yamazaki, ``Watermarking motion data,'' in \emph{Proceedings of Pacific Rim
  Workshop on Digital Steganography}, 2004, pp. 177--185.

\bibitem{yamazaki}
S.~Yamazaki, M.~Mochimaru, and T.~Kanade, ``Watermarking motion clips,'' in
  \emph{Proceedings of the International Conference on Computer Animation and
  Social Agents}, 2005, pp. 171--176,.

\bibitem{modata}
B.~Studio, ``Free motion capture data samples,''
  http://www.beyondmotion.com.au/free\_motion.html, 2009.

\bibitem{11}
R.~Polikar, ``Wavelets in signal processing,'' 
http://users.rowan.edu/\~polikar/WAVELETS/ WTtutorial.html, last Accessed August 5, 2010.

\bibitem{mallat}
S.~Mallat, ``A theory for multiresolution signal decomposition: the wavelet
  representation,'' \emph{IEEE Transactions on Pattern Analysis and Machine
  Intelligence,}, vol.~11, no.~7, pp. 674 --693, 1989.

\bibitem{multiresolutioncurves}
A.~Finkelstein and D.~Salesin, ``Multiresolution curves,'' in \emph{SIGGRAPH:
  Proceedings of the 21st annual conference on Computer graphics and
  interactive techniques}, 1994, pp. 261--268.

\bibitem{dwt}
D.~Philips, ``Filter banks and discrete  wavelet transform,''
http://www.engmath.dal.ca/ courses/engm6610/notes/node6.html, 2003.

\bibitem{mocaptoolbox}
N.~Lawrence, ``The university of manchester: Matlab motion capture toolbox,''
  http://www.cs.man.ac.uk/~neill/mocap/.

\bibitem{cmumodata}
CMU, ``Carnegie mellon university graphics lab motion capture database,''
  http://mocap.cs.cmu.edu/.

\bibitem{last}
R.~Motwani, A.~Ambardekar, M.~Motwani, and F.~Harris, ``Robust watermarking of
  3d skinning mesh animations,'' 2008, pp. 1752 --1756.

\bibitem{motionediting}
J.~Lee and S.~Y. Shin, ``Multiresolution motion analysis and synthesis,''
  Computer Science Department, Korea Advanced Institute of Science and
  Technology, Tech. Rep., 2000.

\end{thebibliography}

\parpic{\includegraphics[width=1in,clip,keepaspectratio]{rakhi.png}}
\noindent {\bf R. C. Motwani} received her B.E. degree in Computer Science from the University of Pune, India, in 2000. She received her M.S. degree from the University of Nevada, Reno (UNR), in 2002 and Ph.D. degree in Computer Science and
Engineering from UNR, in 2010. She is presently an Adjunct Faculty at UNR and provides consulting services to the IT industry. Her research interests lie in the areas of information security, applied artificial intelligence and service oriented architecture.


\parpic{\includegraphics[width=1in,clip,keepaspectratio]{mukesh2.png}}
\noindent {\bf M.C. Motwani} received his B.E. degree in electronics engineering from the University of Pune, India, in 1999, and an M.S. degree from UNR, in 2002, in Computer Ccience. He is presently pursuing his Ph.D. degree in Computer Science and Engineering at UNR and works as a solutions architect to provide consulting services to the IT industry. His research interests lie in Service Oriented Architecture, Digital Rights Management systems, Watermarking, and applied Computational Intelligence

\parpic{\includegraphics[width=1in,clip,keepaspectratio]{fred.png}}
\noindent {\bf F. C. Harris, Jr.} is currently a Professor in the Department of Computer Science and Engineering and the Director of the High Performance Computation and Visualization Lab at the University of Nevada, Reno, USA. He received his B.S. and M.S. in Mathematics and Educational Administration from Bob Jones University in 1986 and 1988 respectively. He got his second M.S. and Ph.D. in Computer Science from Clemson University in 1991 and 1994 respectively. He is a member of ACM, IEEE, and ISCA. His research interests are in Parallel Computation, Graphics and Virtual Reality, and Bioinformatics.

%\begin{biography}[{\includegraphics[width=1in,height=1.25in,clip,keepaspectratio]{rakhi.png}}]{R.C. Motwani}

\end{document}


