## RDP: Ranked Differential Privacy for Facial Feature Protection in Multiscale Sparsified Subspace

Lu Ou Hunan University Changsha, China oulu9676@gmail.com

Shaolin Liao IIT

Chicago, USA sliao5@iit.edu

Shihui Gao Hunan University Changsha, China gaoqaq@hnu.edu.cn

Guandong Huang State Grid Hunan Materials Company Changsha, China huanggd@hnu.edu.cn

Zheng Qin Hunan University Changsha, Country zqin@hnu.edu.cn

Abstract -With the widespread sharing of personal face images in applications' public databases, face recognition systems faces real threat of being breached by potential adversaries who are able to access users' face images and use them to intrude the face recognition systems. Although large amounts of work on facial features privacy preserving methods have been proposed, it is difficult to provide accurate and lightweight protection without degrading the visualization quality of face images, or data utility, under the constraint of given privacy budget for effective facial features protection. In this paper, we propose a novel privacy protection method in the multiscale sparsified feature subspaces to protect sensitive facial features, by taking care of the influence or weight ranked feature coefficients on the privacy budget, named 'Ranked Differential Privacy (RDP)'. After the multiscale feature decomposition, the lightweight Laplacian noise is added to the dimension-reduced sparsified feature coefficients according to the geometric superposition method. Then, we rigorously prove that the RDP satisfies ε 0 -Differential Privacy. After that, the nonlinear Lagrange Multiplier (LM) method is formulated for the constraint optimization problem of maximizing the utility of the visualization quality protected face images with sanitizing noise, under a given facial features privacy budget ε 0 . Then, two methods are proposed to solve the nonlinear LM problem and obtain the optimal noise scale parameters: 1) the analytical Normalization Approximation (NA) method where all Laplacian noise scale parameters are normalized to the average noise scale parameter, which is good for real-time online applications; and 2) the LM optimization Gradient Descent (LMGD) numerical method to obtain the nonlinear solution through iterative updating, which is better for more accurate offline applications. Experimental results on two real-world datasets show that our proposed RDP outperforms other state-of-the-art methods: at a privacy budget of ε 0 = 0 2 . , the PSNR (Peak Signal-to-Noise Ratio) of the RDP is about ∼ 10 dB higher than (10 times as high as) the highest PSNR of all compared methods.

Index Terms -data privacy, multiscale feature space, differential privacy, optimization, face image databases.

## I. INTRODUCTION

In the current digital age, it has become a common phenomenon that individuals' static face images are shared in the public domain, such as Facebook, TikTok, Bing, etc. [1], [2]. Moreover, there are universal and unique facial features with rich biometric information in such images. Existing research shows that facial features provide a solid foundation for face recognition techniques, which plays a vital role in applications such as identity authentication and access control. For example, the face recognition payment function launched by Alipay allows a user to verify his/her identity through scanning his/her face [3]; and it is a common identity authentication application to filter permissions using face recognition techniques [2], [4]. Unfortunately, adversaries could breach face recognition systems by public face images [5]-[8]. Then they can use smart devices without authorization [9], access private systems [4], and even obtain personal sensitive information [2], [10], etc. . As shown in Figure 1, due to the sensitive facial features in face images, once adversaries have obtained original face images in public databases, they could breach face recognitionbased systems. From that, it means that facial features in public face images not only threaten personal privacy, but also may lead to a series of security threats, including identity theft, fraud, and even political and economic espionage, which has been brought up by [11]-[13]. What's more, the privacy leakage even could be caused by only one face image [4]. All of that said, it is necessary to develop an effective facial feature protection method for public face images.

In existing research, cryptography [14], [15], deep learning, [12], [16] and Differential Privacy (DP) [11], [17] are adopted to protect facial features effectively. Although cryptographybased facial feature privacy protection methods can provide strong security, they usually focus on defending against untrusted models or servers, and have expensive computation and communication overhead. Besides, deep learning [16] needs to be trained on large-scale data sets to avoid privacy risks, which is more used for preserving face recognition accuracy while protecting facial features. Different from the previous two kinds of methods, DP [11], [17], [18] is a solid mathematical method that can adds sanitizing noise for protecting facial features while preserving not only the recognition accuracy but also the visual quality of images with released demand. What's more important is that DP can be optimized to maximize the data utility of face images after the sanitizing noise is added, under the constraint of a given privacy budget. Existing methods based on DP [17], [19], [20] may be more likely to add perturbation during the training process, which comes with a large training cost.

In this paper, the data utility is the visualization quality of the noisy face images after the sanitizing noise is added. And the sensitive information is the facial features in face images that could be used to breach some face recognition-based systems. For hiding sensitive information of each individual's face image and conserving the face image's data utility as far

Fig. 1. The process of breaching face recognition systems: for each individual's face recognition, there are the corresponding standard facial features stored in face recognition-based systems, such as face recognition check-in systems. When adversaries download the particular individual's face images from public databases, corresponding to applications, e g . ., Facebook, Weibo, TikTok, or Bing, they may use the downloaded images to breach face recognition-based systems which are used by such individual. On the contrary, if all the face images are protected with a privacy protection method, adversaries will fail to achieve the threat aforementioned.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-216d1b77-2c77-40a8-a538-5061330814d0.jpg)

as possible, the perturbation on the facial features should be minimized while maximizing each noisy face image's visualization quality with a given privacy budget of this image's facial features. And then, we propose such an optimized facial features protection method based on Ranked Differential Privacy (RDP): first, we perform the dimension-reduced sparsification of each original face image for extracting the original facial details into the multiscale feature space and obtaning sparse feature coefficients; then, the influence or weight of each sparse feature coefficient on the total privacy budget ε 0 is evaluated; finally, sanitizing noise of the optimal scale parameters is added to the ranked significant sparse feature coefficients by taking care of the influence or weight of each coefficient. The main contributions of this paper are as follows,

- · First, after reducing the dimension of every face image and obtaining the sparse feature coefficients, the weights of each coefficient with respect to the facial feature vector and thus the total privacy budget ε 0 are evaluated and coefficients are ranked in accordance with such weights. Then, Lagrange Multiplier method (LM) constraint optimization problem is formulated: optimizing the Laplace scale parameters of all feature coefficients with significant influence or weights to maximize the visualization quality of the privacy-preserving noisy face image, under the constraint of the total privacy budget ε 0 of the face image's facial features.
- · What's more important, two methods are proposed to solve the nonlinear LM optimization problem: the Normalized Approximation (NA) analytical method by approximating all weighted Laplacian noise scale parameters to be identical, which is suitable for real-time online applications; and the iteratively updating LM optimization Gradient Descent (LMGD) numerical method, which is first proposed here to achieve more accurate offline applications.
- · Finally, experiments on two real-world datasets show that the facial features have been effectively protected, while the face image's visualization quality is only affected minimally, compared to other state-of-the-art DP meth-

ods.

## II. RELATED WORK

Due to facial features that existed in face images, the images could be used to breach face recognition-based systems [2], [21]. There is the related research on facial features protections. Then, we will review the research in recent years.

A lot of research focus on the problem of how to reserve the efficiency of face recognition while preventing facial features in face images from being leaking for individuals' sensitivity information or face images reconstruction: Chamikara et al . [17] proposed a face recognition privacy protection scheme based on Local Differential Privacy (LDP), which calculates by perturbing the eigenface and storing the perturbed data in a third-party server; Ding et al . [22] proposed a model with dual channels and multiple monitors to prevent DeepFake antiforensic attacks for protecting facial features; Meng et al . [20] proposed a face privacy framework, which extracts purified clusters locally through Differentially Private Local Clustering (DPLC) method and then encourages global consensus among clients to generate more discriminative features; Wang et al . [12] considered that facial features used by face recognition system for authentication may be attempted to reconstruct the original face image and proposed a privacy protection method based on disturbing the mapping from facial features to the image; Lin et al . [23] proposed a privacy-preserving facial features based on a new facial recognition paradigm with a feature compensation mechanism for the challenge between the public demand for powerful face recognition and the individuals' privacy concerns; Mi et al . [24] proposed the MinusFace facial features privacy protection framework for enhancing facial feature protection by random channel shuffling and achieving high accuracy in face recognition system. Others are more interested in preventing adversaries from attacking face recognition models or servers storing the relevant data to obtain facial features: Wang and Leng [15] proposed an image encryption algorithm that combines face recognition and bitonal sequence technologies, and designed a row-column scrambling algorithm to scramble facial features through the characteristics of bitonal sequences; Zhuo et al . [25] proposed an AdaBoost-based framework POR to protect the privacy of users' facial features and the service provider's key learning parameters in face recognition integrated deep learning algorithms. Additionally, there are some studies worrying about the facial features leaking of face images released to the public domain: Ji et al . [19] proposed a frequencydomain DP protection method DCT-DP for facial features, with an optimized privacy budget allocation method based on the loss of the face recognition network; Croft et al . [26] proposed a face obfuscation DP method for generating machine learning models, adding noise to the image pixel intensities for protecting facial features.

All above, such research effectively protects individuals' facial features in face images. However, most studies are more concerned with the trade-off between face recognition accuracy and faial features protection, or consider the unreliability of face recognition models and servers storing related data. In this paper we intend to solve how to release face images with high-quality visual effects on the premise of protecting facial features. Different from [19], we attempt to propose a lightweight perturbation method with adding less noises than [26] in the face images not in training model.

## III. PRELIMINARY WORK

In this section, we will formally define the notations involved in this paper and systematically introduce models and problem statement.

At beginning, notations are shown in Table I.

TABLE I NOTATIONS AND DEFINITIONS

| Symbol   | Description                                                                                                                                                                                                |
|----------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| P        | A pixel matrix of a face image.                                                                                                                                                                            |
| C        | A matrix of sparsed feature coefficients for the face image.                                                                                                                                               |
| ¯ P      | A one-dimensional pixel vector of the face image.                                                                                                                                                          |
| ¯ C      | A one-dimensional sparsed feature coefficient vector of the face image.                                                                                                                                    |
| ¯ C ∗    | A one-dimensional sparsed feature coefficient vector after sorting the elements of ¯ C .                                                                                                                   |
| M P      | The number of elements of ¯ C or ¯ P .                                                                                                                                                                     |
| ξ k      | The noise of the k -th element ¯ C k ∗′ of the noisy sparsed feature coefficients vector be ¯ C ∗′ , 1 ≤ k ≤ M P .                                                                                         |
| b k ¯ F  | The Laplacian scale parameter of ξ k , 1 ≤ k ≤ M P . A feature vector of the face image.                                                                                                                   |
| M F      | The number of elements of ¯ F .                                                                                                                                                                            |
| ∆ i      | The sensitivity of the i -th element ¯ F ′ i of the noisy feature vector ¯ F ′ , 1 ≤ i ≤ M F .                                                                                                             |
| w ik     | The weights between the i -th element ¯ F ′ i of the noisy feature vector ¯ F ′ and the k -th element ¯ C k ∗′ of the ranked noisy sparsed feature coefficient vector ¯ C ∗′ , 1 ≤ i ≤ M F , 1 ≤ k ≤ M P . |
| p        | The probability parameter of the geometric distribu- tion, 0 < p < 1 .                                                                                                                                     |
| K        | The geometric number, K ∼ Geo ( p )                                                                                                                                                                        |
| ε 0      | A constant, which is given as the total privacy budget of ¯ F ′ .                                                                                                                                          |

## A. Models and Problem Statement

Next, we will give the system and adversary models, and state the problem.

1) System Model: With a certain individual A , there exist a large amount of A 's static face images represented as a dataset Ω in public, which can be unrestricted observed and obtained. ¯( F ¯ A ) is the standard feature vector retrieved from A 's standard face ¯ A , used for authentication in some systems with face recognition techniques. For any face image A visually similar to the individual A , its feature vector is ¯ F in the same face recognition with ¯( F ¯ A ) . The system identifies that individual A is passing and can pass through A , if ¯ F and ¯( F ¯ A ) meet the following matching criterion

$$\left | \bar { F } _ { i } - \bar { F } _ { i } ^ { ( \bar { A } ) } \right | \leq \bar { R } _ { i } ; \ i = 1, 2, \cdots, M _ { F }, \quad \ \ ( 1 )$$

where ¯ F i and ¯ ( F i ¯ A ) are the i -th element of ¯ F and ¯( F ¯ A ) with M F elements respectively, and ¯ R i is the i -th element of the feature identification radius vector ¯ R for this authentication system.

2) Adversary and Problem: In this paper, adversaries are considered to be interest to and have the ability to download the face image dataset Ω of the certain individual A from the public domain, and they could obtain information about authentication devices equipped with the face recognition system, which are commonly used by individual A . Adversaries may attempt to breach these face recognition system by using such face image in Ω , and pass the authentication to obtain A 's privacy information.

For the face image dataset Ω of A , our goal is to develop an accurate and lightweight facial feature vector privacypreserving method f ′ that sanitizes these images by adding noise, i e . . A ′ = f ′ ( A ) for any face image A ∈ Ω . Such method can ensure that the noisy face image A ′ can prevent adversaries from breaching face recognition-based systems while having extremely high visual quality, meaning that we aim to optimize the utility of the sanitized images by maximizing images' visualization quality while providing privacy protection for facial features in face images.

## B. Basic Concepts

In this section, we give some preliminary definitions for DP. First, in the terminology of DP, we define a pair of neighboring face images, expressed in Definition 1 .

Definition 1 (Neighboring Face Images) . If a face image A belongs to the individual A and a face image B does not belong to the individual A , The face images A and B are a pair of neighboring face images.

Then, the sensitivity of facial features is defined as Definition 2 .

Definition 2 (The Sensitivity of Facial Features) . If M is a measurement of face images' facial features, its sensitivity is

$$\Delta _ { i } = \max _ { A, B } | | \mathcal { M } ( A ) - \mathcal { M } ( B ) | | _ { 1 },$$

where A and B are neighboring face images.

Finally, the definition of DP is as follows.

Definition 3 (Differential Privacy) . There is a measurement M ′ and a pair of neighboring face images A and B . M ′ is ε 0 -differentially private, if the following is satisfied

$$\left | \frac { \Pr \left [ \mathcal { M } ^ { \prime } ( A ) = \bar { F } ^ { \prime } \right ] } { \Pr \left [ \mathcal { M } ^ { \prime } ( B ) = \bar { F } ^ { \prime } \right ] } \right | \leq \exp \left ( \varepsilon _ { 0 } \right ).$$

## IV. RANKED DIFFERENTIAL PRIVACY

In this section, we introduce our proposed Ranked Differential Privacy (RDP) for facial features with weighted sparsified subspace. Overall, before releasing the original face images, the Laplacian noise is added to the face images' facial features by geometric superposition [27], [28], in order to hide the facial features which are used in the process of face recognition of the systems.

## A. Perturbation in Multiscale Feature Space

As well known, the face images' facial features play a key role in the implementation of face recognition. So that they should be protected to avoid adversaries' breaching face recognition-based systems. From this, before the perturbation, facial features should be extracted in the multiscale feature space. Because of its efficiency and rigorous mathematical expression, Multi-level Haar Wavelet Transform (HWT) [29] is adopted to reduce the dimensionality of the face image's pixels and obtain accurate information about the face image's facial features, i e . ., the wavelet coefficient vector constructed by the sparsified wavelet coefficients is obtained. After the perturbation on the wavelet coefficient vectors, the noisy face image could be generated by the noisy pixels matrix P ′ obtained by the noisy coefficients vector ¯ C ′ through Multilevel Inverse Wavelet Transform (MIWT).

During the extraction and dimensionality reduction of facial features from face images, we use the first-order Taylor expansion to describe the liner mapping relationship between the noise of the sparse wavelet space and the noisy reduced dimensionality facial features of the noisy face image. After adding Laplacian noise to the elements of the original wavelet coefficient vectors ¯ C , we obtain the noisy feature vector ¯ F ′ by the noisy face image, which comes from the noisy wavelet coefficient vectors ¯ C ′ with MIWT. According to the first-order Taylor formula for the i -th element ¯ F ′ i of the noisy facial feature vector ¯ F ′ , we have

$$\bar { F } _ { i } ^ { \prime } = \bar { F } _ { i } + \sum _ { j = 1 } ^ { M _ { P } } \frac { \partial \bar { F } _ { i } ^ { \prime } } { \partial \bar { C } _ { j } ^ { \prime } } \cdot \left ( \bar { C } _ { j } ^ { \prime } - \bar { C } _ { j } \right ) + o \left ( \bar { C } ^ { \prime } \right ), \quad \left ( 2 \right ) \quad \text{where}$$

where 1 ≤ i ≤ M F ; ¯ C ′ j is the j -th element of ¯ C ′ , where M P is the number of the elements of the face image's pixel matrix and 1 ≤ j ≤ M P ; and o ( ¯ C ′ ) is the minimal term, i e . ., o ( ¯ C ′ ) ∼ 0 . Thus, we think

$$\bar { F } _ { i } ^ { \prime } = \bar { F } _ { i } + \sum _ { j = 1 } ^ { M _ { P } } \frac { \partial \bar { F } _ { i } ^ { \prime } } { \partial \bar { C } _ { j } ^ { \prime } } \cdot ( \bar { C } _ { j } ^ { \prime } - \bar { C } _ { j } ).$$

Obviously, ¯ C ′ j -¯ C j represents the noise of ¯ C ′ j , and the noise of ¯ F ′ i can be considered as the weighted sum of elements' noises of ¯ C ′ , with that the first-order partial derivative of ¯ F ′ i on ¯ C ′ j is regarded as the weight for the noise of ¯ C ′ j .

k After this, with the perturbation provides by weighted sum, we can select the wavelet coefficients more significant in facial features from the sorted wavelet coefficients basis of such weights for adding Laplacian noises with appropriate scale to wavelet coefficients, thereby providing stricter privacy preservation. The definition of the noise mechanism based on such ranked wavelet coefficients is shown in Definition 4.

Then, we define the weight w ik for the i -th element ¯ F ′ i of the noisy feature vector ¯ F ′ and the k -th element ¯ C k ∗′ of the sorted noisy wavelet coefficient vector ¯ C ∗′ as w ik = ∂F ¯ ′ i ∂C ¯ ∗′ .

Definition 4 (Ranked Differential Privacy) . Given the geometric number K between 1 and M P with mean 1 p , 0 &lt; p &lt; 1 . For obtaining the noisy wavelet coefficients vector ¯ C ′ , K Laplacian noises with mean 0 are added to the top K of the wavelet coefficients vector ¯ C 's ranked elements of a face image, with ranking elements of ¯ C ′ by the corresponding weights between the noisy feature vector's elements and such elements. After the inverse transformation with MIWT of ¯ C ′ , the noisy pixel matrix P ′ is obtained from ¯ C ′ , and the noisy face image is generated by P ′ . Such method, providing perturbation to the elements of the feature vector ¯ F by the weighted sum of ranked wavelet coefficients' Laplacian noises for protecting facial features, is called Ranked Differential Privacy, referred as 'RDP'.

In RDP, after perturbation of wavelet coefficients, the noisy pixels matrix P ′ is obtained by ¯ C ∗′ through MIWT, and the noisy feature vector of the noisy face image generated by the noisy face image from P ′ is ¯ F ′ with M F elements. Let the noise of the i -th element ¯ F ′ i in ¯ F ′ be ζ i , ¯ F ′ i = ¯ + F i ζ i , with 1 ≤ i ≤ M F . Then we have Theorem 1 of RDP about ¯ F ′ .

Theorem 1. The noise ζ i of the i -th element ¯ F ′ i in ¯ F ′ ,which is obtained by RDP, satisfies Laplace distribution, i e . ., Pr[ ζ i ] = Lap ( 0 , ∆ i ε i ) , where ε i is the privacy budget for the i -th element ¯ F ′ i of ¯ F ′ , and 1 ≤ i ≤ M F .

Proof: Let ¯ C ∗′ and ¯ C ∗ be the noisy ranked wavelet coefficient vector and the ranked wavelet coefficient vector, separately. According to Eq. (2), we have the following firstorder Taylor expansion for ¯ F ′ i ,

$$\bar { F } _ { i } ^ { \prime } & = \bar { F } _ { i } + \sum _ { \substack { k = 1 \\ - } } ^ { M _ { P } } \frac { \partial \bar { F } _ { i } ^ { \prime } } { \partial \bar { C } _ { k } ^ { * \prime } } \cdot \left ( \bar { C } _ { k } ^ { * \prime } - \bar { C } _ { k } ^ { * } \right ) + o \left ( \bar { C } ^ { * \prime } \right ), \\ - \quad. \quad \, - \,. \quad \, -$$

where ¯ C k ∗′ and ¯ C ∗ k are the k -th element of ¯ C ∗′ and ¯ C ∗ , separately, and o ( ¯ C ∗′ ) are the minimal term, i e . ., o ( ¯ C ∗′ ) ∼ 0 . Thus, we think

$$\bar { F } _ { i } ^ { \prime } = \bar { F } _ { i } + \sum _ { k = 1 } ^ { M _ { P } } \frac { \partial \bar { F } _ { i } ^ { \prime } } { \partial \bar { C } _ { k } ^ { * ^ { \prime } } } \cdot ( \bar { C } _ { k } ^ { * ^ { \prime } } - \bar { C } _ { k } ^ { * } ).$$

The weight for the i -th element ¯ F ′ i of the noisy feature vector ¯ F ′ and the k -th element ¯ C k ∗′ of the sorted noisy wavelet coefficient vector ¯ C ∗′ is w ik , w ik = ∂F ¯ ′ i ∂C ¯ k ∗′ . K is the given geometric number between 1 and M P with mean 1 p , 0 &lt; p &lt; 1 . After we add Laplacian noises to the top K elements of the ranked wavelet coefficient vector ¯ C ∗ according to RDP, the first-order Taylor expansion of ¯ F ′ i is further

$$\bar { F } _ { i } ^ { \prime } = \bar { F } _ { i } + \sum _ { k = 1 } ^ { M _ { P } } \frac { \partial F _ { i } ^ { \prime } } { \partial \bar { C } _ { k } ^ { * ^ { \prime } } } \left ( \bar { C } _ { k } ^ { * ^ { \prime } } - \bar { C } _ { k } ^ { * } \right ) = \bar { F } _ { i } + \sum _ { k = 1 } ^ { M _ { P } } w _ { i k } \xi _ { k }.$$

Then due to there are K noises, ξ k = { ξ k , 1 ≤ k ≤ K 0 , k ≤ k ≤ M P , and ¯ F ′ i = ¯ + F i K ∑ k =1 w ik ξ k . From that, we have ζ i = K ∑ k =1 w ik ξ k , and can get the probability of ζ i , Pr[ ζ i ] = Pr [ K ∑ k =1 w ik ξ k ] . And based on the geometric superposition, ζ i satisfies Laplace distribution, i e . ., ζ i ∼ Lap ( 0 , ∆ i ε i ) with sensitivity ∆ i and privacy budget ε i of ¯ F ′ i . Especially, the scale parameter b k of the noise ξ k is related to the privacy budget ε i of ¯ F ′ i , 1 ≤ k ≤ M P .

This theorem has been proved.

□

Next, the distribution of the noisy feature vector ¯ F ′ can be derived from elements of ¯ F ′ . Based on Theorem 1 , we can obtain the relationship between the total privacy budget of ¯ F ′ and the scale parameters of Laplacian noises in ¯ C ∗′ .

## B. Privacy Budget Analysis

In face recognition, the facial feature vector is obtained by projecting the face image onto the eigenfaces, which are constructed by a set of orthogonal eigenvectors through SVD decomposition and dimensionality reduction from the original face images dataset [30]. Such eigenvectors capture the most significant variations of features and represent the main change directions in the feature space, and their corresponding singular values represent the importance of eigenfaces. Due to their orthogonality, these eigenvectors are independent and can serve as the basis vectors of the feature space, enabling the representation of the original face image. Consequently, the feature vector's elements, which are projections of the original image onto these eigenvectors, exhibit independence.

According to Theorem 1 , ζ i = K ∑ k =1 w ik ξ k , where K is the given geometric number with mean 1 p , 0 &lt; p &lt; 1 , and w ik is the weight between ¯ F ′ i and ¯ C k ∗′ , w ik = ∂F ¯ ′ i ∂C ¯ k ∗′ , 1 ≤ k ≤ M P . Let the sensitivity ∆ i and the privacy budget ε i of ¯ F ′ i , we have Pr[ ζ i ] = Lap ( 0 , ∆ i ε i ) . According to the variance operation, the variance of ζ i is

$$V a r \left [ \zeta _ { i } \right ] = 2 \left ( \frac { \Delta _ { i } } { \varepsilon _ { i } } \right ) ^ { 2 }.$$

$$\text{is} \quad \text{and the variance of } \sum _ { k = 1 } ^ { K } w _ { i k } \xi _ { k } \text{ is}$$

$$\varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varp \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \vp \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi \varphi
\vdots \leq K \\ \vdots \sum _ { k = 1 } ^ { M _ { P } } w _ { i k } \xi _ { k }, \\ w _ { i k } \xi _ { k } \Big. \quad \text{Because of } V a r [ \zeta _ { i } ] \, = \, V a r \Big [ \sum _ { k = 1 } ^ { K } w _ { i k } \xi _ { k } ] \, \text{ from } \Pr [ \zeta _ { i } ] \, =$$

Because of V ar [ ζ i ] = V ar [ K ∑ k =1 w ik ξ k ] from Pr[ ζ i ] =

Pr [ K ∑ k =1 w ik ξ k ] , we get the expression of ε i with respect to b k as

$$\stackrel { \dots } { \colon } \quad \varepsilon _ { i } \left ( b _ { k } \right ) = \frac { \Delta _ { i } } { \sqrt { \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { M _ { P } + 1 } \right ] } }.$$

And based on that the elements of ¯ F ′ are independent, the total privacy budget of ¯ F ′ is M F ∑ i =1 ε i ( b k ) , then we have the total privacy budget's expression ε ( b k ) of ¯ F ′ with respect to b k as

$$\int _ { \substack { \int _ { i = 1 } \\ \text{and} \\ \text{h are} \\ \text{SVD} } } ^ { \int _ { \substack { \int _ { M _ { P } } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 } \\ \int _ { k = 1 }$$

## C. Optimal Data Utility under Privacy Budget

Our optimization goal is to maximize the noisy face image's visualization quality under the constraint of the privacy budget, i e . ., optimizing the noise scale parameters of all wavelet coefficients, which is a constraint optimization problem.

Since the variance of the noise quantifies the negative impact of the noise on the visualization quality and a constant does not have influence on the variance, the cost function for the visualization quality of the face image can be defined by deriving the variance of the noisy pixel matrix, of which value is smaller with the smaller loss of visualization quality.

Further, the constrained optimization problem can be converted into minimizing the cost function under the constraint of privacy budget. Then, we define the cost function as Definition 5.

Definition 5 (Cost Function) . After multi-level HWT for the original face image, ¯ C is the one-dimensional vector that the wavelet coefficient matrix expansion to. ¯ C ∗ is a onedimensional vector after sorting the elements of ¯ C . According to RDP (Definition 4), the Laplacian noises are added to ¯ C ∗

and get the ranked noisy one-dimensional wavelet coefficient vector ¯ C ∗′ . Define that b k is the scale parameter of the Laplacian noise of the k -th element of ¯ C ∗′ , 1 ≤ k ≤ M P , then the expression of the cost function f ( b k ) of the noisy face image with respect to b k is

$$f \left ( b _ { k } \right ) = \sum _ { k = 1 } ^ { M _ { P } } \left ( b _ { k } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { 1 + M _ { P } } \right ]. \quad \left ( 4 \right ) \text{We} \quad \text{grams} \quad \text{$row$ can} \quad.$$

Proof: Given the pixel matrix P of a face image, and it undergoes a multi-level wavelet transform to obtain a matrix of wavelet coefficients C . Assume that the transformation matrix of wavelet transform is ˜ H [31]-[33], which is a √ M P × √ M P orthogonal matrix. After N -level wavelet transform, let H =

$$\tilde { H } ^ { N }, \, \text{and} \, \mathbf H ^ { T } = \left ( \tilde { \mathbf H } ^ { T } \right ) ^ { N }. \, \text{Then}$$

$$C = \tilde { H } ^ { N } P = { \mathrm H } P, \text{ and } P = \left ( { \mathrm H } ^ { T } \right ) ^ { N } C = { \mathrm H } ^ { T } C,$$

where H H T = I , and I is an identity matrix.

Define the Laplacian noise matrix of wavelet coefficient matrix C is N , and the noisy wavelet coefficient matrix is C ′ = C + N . Then for the noisy pixel matrix P ′ ,

$$\mathbf P ^ { \prime } = \left ( \tilde { \mathbf H } ^ { T } \right ) ^ { N } \mathbf C ^ { \prime } = \mathbf H ^ { T } \left ( \mathbf C + \mathbf N \right ) = \mathbf P + \mathbf H ^ { T } \mathbf N.$$

Let B be the Laplacian scale parameter matrix of N N , mn ∼ Lap (0 , B mn ) , 1 ≤ m ≤ M P , 1 ≤ n ≤ M P . Thus,

$$L a p & \left ( 0, \mathbf B _ { m n } \right ), \, 1 \leq m \leq M _ { P }, \, 1 \leq n \leq M _ { P }. \, \text{ Thus}, \\ V a r [ \mathbf P ^ { \prime } ] & \ = 2 \sum _ { m = 1 } ^ { \sqrt { M _ { P } } \sqrt { M _ { P } } \sqrt { M _ { P } } } \sum _ { n = 1 } ^ { \sqrt { M _ { P } } } \left [ \left ( \mathbf H ^ { T } \right ) _ { m t } \right ] ^ { 2 } \cdot \left ( \mathbf B _ { t n } \right ) ^ { 2 } } \\ & \ = 2 \sum _ { n = 1 } ^ { \sqrt { M _ { P } } \sqrt { M _ { P } } } \sum _ { t = 1 } ^ { M _ { P } } \left [ \left ( \mathbf H ^ { T } \right ) _ { m t } \right ] ^ { 2 } } \\ & \ = 2 \sum _ { n = 1 } ^ { \sqrt { M _ { P } } \sqrt { M _ { P } } } \sum _ { t = 1 } ^ { \sqrt { M _ { P } } } \left ( \mathbf B _ { t n } \right ) ^ { 2 } } \\ & \ = V a r [ \mathbf C ^ { \prime } ]. & \quad \ \text{c} \\ \text{For the convenience of calculation, we flatten the pixel}$$

For the convenience of calculation, we flatten the pixel matrix P and the wavelet coefficient matrix C into the onedimensional pixel vector ¯ P and the wavelet coefficient vector ¯ C C . ¯ ∗ is a one-dimensional vector after sorting the elements of ¯ C . Then the variance of the noisy pixel vector ¯ P ′ and the noisy wavelet coefficient vector ¯ C ′ is equal, i e . .

$$\text{my wavelet connection vector} & \colon = \text{equal, } i e. \\ \ V a r \left [ \bar { P } ^ { \prime } \right ] & \quad = \ V a r \left [ \bar { C } ^ { * \prime } \right ] & \quad \text{th} \\ & \quad = 2 \sum _ { K = 1 } ^ { M _ { P } } p ( 1 - p ) ^ { K } \sum _ { k = 1 } ^ { K } \left ( b _ { k } \right ) ^ { 2 } & M \\ & \quad = 2 \sum _ { k = 1 } ^ { M _ { P } } \left ( b _ { k } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { M _ { P } + 1 } \right ], \\ \text{where} & \bar { C } ^ { * \prime } \text{ is the banding noise} \text{ in a minimal natural equal}$$

where ¯ C ∗′ is the ranked noisy one-dimensional wavelet coefficient vector.

From that, the expression of the cost function of ¯ P ′ with respect to b k is as follows

$$f \left ( b _ { k } \right ) = \sum _ { k = 1 } ^ { M _ { P } } \left ( b _ { k } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { M _ { P } + 1 } \right ],$$

which is the cost function of the noisy face image.

□

From that, a constant ε 0 is given as the total privacy budget of the noisy feature vector ¯ F ′ , and the constraint is ε ( b k ) = ε 0 . Thus, the optimization problem is as follows,

$$b _ { k } = \underset { b _ { k } } { \arg \min } \left \{ f \left ( b _ { k } \right ) \left | \varepsilon \left ( b _ { k } \right ) = \varepsilon _ { 0 } \right \}.$$

We solve this optimization problem according to the Lagrangian Multiplier method (LM) to better balance the need to sanitize the facial features and preserve the visualization quality of the face image. According to LM, the Lagrangian function with respect to b k is

$$L \left ( b _ { k }, \lambda \right ) = f \left ( b _ { k } \right ) + \lambda \left ( \varepsilon \left ( b _ { k } \right ) - \varepsilon _ { 0 } \right ),$$

which can be solved by setting the first derivatives with respect to λ and b k to zeros,

$$\begin{cases} \ \frac { \partial L ( b _ { k } ) } { \partial \lambda } = 0 ; \\ \ \frac { \partial L ( b _ { k } ) } { \partial b _ { k } } = 0, \end{cases}$$

where the first equation reduces to the following,

$$\int _ { m n } ^ { M _ { F } } \frac { \lambda \Delta _ { i } ( w _ { i k } ) ^ { 2 } } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 3 } { 2 } } } = 2. \ ( 5 )$$

It is clear that Eq. (5) is nonlinear. Here we propose two ways to solve the nonlinear LM equations for two different applications, i e . ., 1) the Normalization Approximation (NA) for the real-time online applications; and 2) the LM optimization Gradient Descent (LMGD) for the offline applications that need more accurate result.

- 1) Optimization with Normalization Approximation: It is clear that Eq. (5) is nonlinear, which makes it difficult to find the analytical solution. Therefore, we propose the optimization with Normalization Approximation (NA), which assumes that all wavelet coefficients of interest contribute equally to the total privacy budget ε 0 .

To achieve such goal, it would be intuitive to choose the weighted Laplacian noise scale parameters b k of Eq. (5) to be the average value, i e . ., w ik b k ≈ ¯ = b MF ∑ i =1 w ik b k M F , where 1 ≤ i ≤ M F and 1 ≤ k ≤ M P . And substituting ¯ b into Eq. (5), we get

$$\underset { \text{oe} } { \sum } ^ { M _ { F } } \frac { \lambda \Delta _ { i } ( w _ { i k } ) ^ { 2 } } { \bar { b } ^ { 3 } \left \{ \sum _ { k = 1 } ^ { M _ { P } } ( w _ { i k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 3 } { 2 } } } = 2.$$

It can be substituted by w ik b k to replace ¯ b for obtaining the equations about every b k and λ , which can be used to solve ∂L b ( k ) ∂b k = 0 . And then, Eq. (5) can be solved as follows,

$$b _ { k } = \frac { \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } } { w _ { i k } } } { \varepsilon _ { 0 } \left [ \frac { 1 - p - ( 1 - p ) ^ { 1 + M _ { P } } } { p } - M _ { P } ( 1 - p ) ^ { 1 + M _ { P } } \right ] ^ { \frac { 1 } { 2 } } }.$$

2) LM Optimization via Gradient Descent: Although NA can solve the system of nonlinear equations in this optimization problem in real time, the assumption as its basis is are very idealistic, which may cause errors far from the best result. For the case allowing solving this problem offline, we propose solving the LM optimization problem via Gradient Descent (LMGD), which continuously tunes the initial value of b k until an optimal solution is obtained through iterative gradient descent updating.

It aims to find the best solution for b k and λ by solving following equations,

$$\begin{cases} \ \frac { \partial L ( b _ { k } ) } { \partial \lambda } = 0 ; \\ \ \frac { \partial L ( b _ { k } ) } { \partial b _ { k } } = 0, \end{cases}$$

which can be solved with two alternating gradient descent updates:

- 1) gradient descent update with respect to b k ;
- 2) gradient descent update with respect to λ .

However, simple alternating gradient descent updates have the problem of finding the optimal step length or learning rate, which is difficult for non-convex optimization problems. Also, the nonlinear operation is required to achieve the minimum b k , which motivates us to develop the LMGD to solve the nonlinear LM problem.

According to LM, the purpose of LMGD is to obtain the optimal solution of b k and λ , through solving the following equations,

$$\begin{cases} \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } } { \sqrt { \sum _ { k = 1 } ^ { M _ { P } } ( w _ { i k } b _ { k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] } } - \varepsilon _ { 0 } = 0, \\ \sum _ { i = 1 } ^ { M _ { F } } \frac { \lambda \Delta _ { i } ( w _ { i k } ) ^ { 2 } } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } ( w _ { i k } b _ { k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 3 } { 2 } } } - 2 = 0. \\ \intertext { W a d i n a t h o r f i n t o n e o f t w o n d i t i o n c h t a c l e i n } \end{cases}$$

We define the error functions of two updating sub-tasks in LMGD as follows,

$$& L M G D \text{ as follows, } & \frac { \partial \mathcal { F } _ { 2 } } { F _ { 1 } } \left ( b _ { k } ^ { ( h ) } \right ) & \frac { \frac { \partial \mathcal { F } _ { 2 } } { \longrightarrow } } \\ & = \left | \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } } { \sqrt { \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { 1 + M _ { P } } \right ] } } - \varepsilon _ { 0 } \right | ; \quad = \left | 3 \frac { \Lambda } { 2 } \\ & \quad \partial F _ { 2 }$$

$$& F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) \\ = & \left | \sum _ { i = 1 } ^ { M _ { F } } \frac { \lambda ^ { ( h ) } \Delta _ { i } ( w _ { i k } ) ^ { 2 } } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 3 } { 2 } } } - 2 \right |, \quad = \left | \sum _ { i = 1 } ^ { M _ { F } } \\ \text{where $h$ represents the $h$-th parameter update. From that, the} \quad \text{$c$} \cdots \varepsilon$}.$$

where h represents the h -th parameter update. From that, the total loss function is

$$F \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) = F _ { 1 } \left ( b _ { k } ^ { ( h ) } \right ) + F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ).$$

,

When F 1 ( b ( h ) k ) → 0 with F 2 ( b ( h ) k , λ ( h ) ) → 0 , the closer b k and λ are to the optimal solution.

̸

The LMGD takes an initial solution b (0) k of b k and λ (0) of λ as input, where ε ( b (0) k ) ≈ ε 0 , however ε ( b (0) k ) = ε 0 . After the initial values are inputted, the optimal solution of b k and λ is obtained through the iterative updating process updated by two alternating gradient descent until the loss function F ( b ( h ) k , λ ( h ) ) is less than a minimum value δ , which is very close to 0 .

Given the step length or learning rate of η . Now the iterative updating process takes two alternating updates as follows,

(1) Step 1: firstly, b ( h ) k is updated iteratively by minimizing the loss function F 1 ( b ( h ) k ) using the gradient descent, while b ( h ) k and λ ( h ) are updated iteratively by minimizing the loss function F 2 ( b ( h ) k , λ ( h ) ) using gradient descent, i e . .,

$$b _ { k } ^ { ( h ) } & = b _ { k } ^ { ( h ) } - \eta \left [ \frac { \partial F _ { 1 } \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } } + \frac { \partial F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } } \right ] ; \\ \lambda ^ { ( h ) } & = \lambda ^ { ( h ) } - \eta \frac { \partial F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) } { \partial \lambda ^ { ( h ) } }, \\ \cdots$$

where the gradient ∂F 1 ( b ( h ) k ) ∂b ( h ) k can be obtained as follows

$$\underset { \left ( \sum } { \frac { \partial F _ { 1 } \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } } } = \left | \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } b _ { k } ^ { ( h ) } ( w _ { i k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 4 } { 2 } } } \right |.$$

(2) Step 2: then, the gradients of ∂F 2 ( b ( h ) k ,λ ( h ) ) ∂b ( h ) k and ∂F 2 ( b ( h ) k ,λ ( h ) ) ∂λ ( h ) are calculated by

$$\frac { \partial F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } }$$

$$\Big. \, \varepsilon _ { 0 } \Big | \, ; \quad = \Big | 3 \sum _ { i = 1 } ^ { M _ { F } } \frac { \lambda ^ { ( h ) } b _ { k } ^ { ( h ) } \Delta _ { i } ( w _ { i k } ) ^ { 4 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 5 } { 2 } } \Big | } \, ; \\ \quad \frac { \partial F _ { 2 } \left ( b _ { k } ^ { ( h ) }, \lambda ^ { ( h ) } \right ) } { \partial \lambda ^ { ( h ) } } \\ \Big | \,. \, \Big | \,.$$

$$\bar { \psi } ^ { - } \, & \, = \left | \sum _ { i = 1 } ^ { M _ { F } } \frac { \partial \lambda ^ { ( h ) } } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { 1 + M _ { P } } \right ] \right \} ^ { \frac { 3 } { 2 } } } \right | \,. \\ & \, \cap \, \epsilon. \quad \epsilon. \quad \epsilon. \quad \epsilon. \quad \cdot. \quad \cdot. \quad \cdot. \quad \cdot.$$

(3) Step 3: finally, after alternatively repeating Step 1 and Step 2 above, the optimal solution of b k and λ are obtained until F ( b ( h ) k , λ ( h ) ) &lt; δ .

The process by Eq. (6) is proved as follow:

Proof: If the noise's scale parameter of the i -th element F i ′ of the noisy feature vector ¯ F ′ is b i ( ¯ F ) , then the expression of b i ( ¯ F ) with respect to b ( h ) k is

$$b _ { i } ^ { \left ( \bar { F } \right ) } \left ( b _ { k } ^ { \left ( h \right ) } \right ) = \sqrt { \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { \left ( h \right ) } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { M _ { P } + 1 } \right ] },$$

and the expression of the total privacy budget ε with respect to b ( h ) k is

$$\varepsilon \left ( b _ { k } ^ { ( h ) } \right ) = \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } } { \sqrt { \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ \left ( 1 - p \right ) ^ { k } - \left ( 1 - p \right ) ^ { M _ { P } + 1 } \right ] } }, \begin{array} { c } B. \\ \text{Let } \bar { F } \\ \text{visual} \mathbb { A } \\ i \text{-th elem} \end{array}$$

where ∆ i is the sensitivity of the i -th element ¯ F ′ i of ¯ F ′ , 1 ≤ i ≤ M F , w ik is the weights between ¯ F ′ i and ¯ C k ∗′ , 1 ≤ k ≤ M P , and p is the probability parameter greater than 0 and less than 1 .

According to the chain rule,

$$\frac { \partial \varepsilon \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } } = \frac { \partial \varepsilon \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { i } ^ { \left ( \bar { F } \right ) } \left ( b _ { k } ^ { ( h ) } \right ) } \cdot \frac { \partial b _ { i } ^ { \left ( \bar { F } \right ) } \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } }.$$

Due to the independence of the elements of the feature vector,

$$\hat { \Phi } _ { \quad } \frac { \partial \varepsilon \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { i } ^ { \left ( \bar { F } \right ) } } = - \frac { \Delta _ { i } } { \left ( b _ { i } ^ { \left ( \bar { F } \right ) } \right ) ^ { 2 } },$$

and

$$\frac { \partial b _ { i } ^ { ( F ) } } { \partial b _ { k } ^ { ( h ) } } & = \frac { b _ { k } ^ { ( h ) } ( w _ { i k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { M _ { P } + 1 } \right ] } { \sqrt { \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { M _ { P } + 1 } \right ] } }. \\ \quad \_ \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

Next, we substitute the above expression into the derivative ∂ε ( b ( h ) k ) ∂b ( h ) k , and obtain it as follows

$$\frac { \partial F _ { 1 } \left ( b _ { k } ^ { ( h ) } \right ) } { \partial b _ { k } ^ { ( h ) } } & & \quad \partial F _ { 1 } \left ( b _ { k } ^ { ( h ) } \right ) & & \quad \text{Now} \\ = & \left | \sum _ { i = 1 } ^ { M _ { F } } \frac { \Delta _ { i } b _ { k } ^ { ( h ) } ( w _ { i k } ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { M _ { P } + 1 } \right ] } { \left \{ \sum _ { k = 1 } ^ { M _ { P } } \left ( w _ { i k } b _ { k } ^ { ( h ) } \right ) ^ { 2 } \left [ ( 1 - p ) ^ { k } - ( 1 - p ) ^ { M _ { P } + 1 } \right ] \right \} ^ { \frac { 3 } { 2 } } \right | } { \left | }. & & \quad \text{from wl} \\ & & \quad \text{Data} \\ & & \quad \text{(Peak S}. \\ & & \quad \Box \quad ( i. e. \ \text{vis} \quad \text{measure} \\ & & \quad \text{V. ANALYSIS OF PRIVACY AND UTILITY}$$

## V. ANALYSIS OF PRIVACY AND UTILITY

Privacy. Now let's look at the privacy settings of the RDP for the feature vectors ¯ F ′ of the noisy face image according to Definition 4. In this paper, the elements of the feature vector are independent of each other with the given total privacy budget ε 0 .

Theorem 2. RDP satisfies ε 0 -differential privacy.

Proof: The face images A and B are a pair of neighboring face images. According to Definition 3, proving that RDP satisfies ε 0 -differentially privacy is equivalent to prove

$$\left | \frac { \Pr \left [ \bar { F } ^ { ( A ) \prime } \right ] } { \Pr \left [ \bar { F } ^ { ( B ) \prime } \right ] } \right | \leq \exp \left ( \varepsilon _ { 0 } \right ),$$

where ¯ F ( A ) ′ and ¯ F ( B ) ′ is the noisy feature vectors of A and B .

Let ¯( F ¯ A ) be the standard feature vector of a certain individual A 's standard face ¯ A . Then, ∆ i is the sensitivity fot the i -th element of the noisy feature vector, and the Laplace scale parameter of that is b ( ¯ F ) i . Define the privacy budget forthe i -th element of the noisy feature vector as ε i , ε i = ∆ i b ( ¯ F ) i ,

$$\stackrel { - 1 } { \colon s s } 1 \leq i \leq M _ { F }. \text{ Let } \varepsilon _ { 0 } = \sum _ { i = 1 } ^ { M _ { F } } \varepsilon _ { i },$$

$$1 \leq k \leq \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ \text{$1\leq i\leq M_{F}$. Let $\varepsilon_{0}=\sum_{i=1}$} \\ \left | \Pr \left [ \bar { F } ^ { ( A ) ^ { \prime } } \right ] \right | = \prod _ { i = 1 } ^ { M _ { F } } \left | \Pr \left [ \bar { F } ^ { ( A ) ^ { \prime } } \right ] \right | \\ \left | \Pr \left [ \bar { F } ^ { ( B ) ^ { \prime } } \right ] \right | \\ \text{$e\text{ vector,} } = \prod _ { i = 1 } ^ { M _ { F } } \left | \frac { \exp \left ( - \left | \bar { F } ^ { ( A ) ^ { \prime } } - \bar { F } ^ { ( A ) } \right | } { b _ { i } ^ { ( \bar { F } ) } \right ) \right | } \\ \exp \left ( - \frac { \left | \bar { F } ^ { ( B ) ^ { \prime } } - \bar { F } ^ { ( i } } { b _ { i } ^ { ( \bar { F } ) } \right | } \right | } { b _ { i } ^ { ( \bar { F } ) } \right ) \right | \\ = \prod _ { i = 1 } ^ { M _ { F } } \exp \left ( \frac { \left | \bar { F } ^ { ( B ) ^ { \prime } } - \bar { F } ^ { ( i } } { \right | } - \left | \bar { F } ^ { ( A ) ^ { \prime } } - \bar { F } ^ { ( i } } ^ { ( \bar { A } ) } \right | } { b _ { i } ^ { ( \bar { F } ) } \right ) \right | \\ \left | \frac { b _ { i } ^ { ( \bar { F } ) } } { b _ { i } ^ { ( \bar { F } ) } \right | } \right | \\ \right | \text{$b_{i}^{\prime}$} } \right | \\ \text{$e\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{ $\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{	$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text{$\text{$\text{$\text{$\text{$\text$$\text$$\text$$\text$$\text$$\text$$\text$$\text$$\text$\text$$\text$$\text$\text$$\text$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$$\text$\text$\text$\text$$\text$\text$\text$\text$\text$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$	\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\attrib$}$}$}$} \right | } \\ \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$.} } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } ) } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } \ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } _ } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } } }$$

Now we have

$$\left | \frac { \Pr \left [ \bar { F } ^ { ( A ) \prime } \right ] } { \Pr \left [ \bar { F } ^ { ( B ) \prime } \right ] } \right | \leq \exp \left ( \varepsilon _ { 0 } \right ),$$

from which Theorem

2 is proved. □

Data utility. In this paper, we use the face image's PSNR (Peak Signal-to-Noise Ratio) value to represent the data utility ( . i e . visualization quality) of the face image. PSNR is a measure to evaluate image quality. The higher the value, the less the face image distortion. PSNR of the noisy face image is

$$P S N R = 1 0 \cdot \log _ { 1 0 } \left ( \frac { \left ( \max \left \{ \bar { P } ^ { \prime } \right \} \right ) ^ { 2 } } { \sigma _ { r } ^ { 2 } } \right ), \quad \ \ ( 7 )$$

where ¯ P ′ is the one-dimensional pixel vector after the noisy pixel matrix P ′ flattened, max { ¯ P ′ } is the maximum element of ¯ P ′ , and σ r 2 is the real variance of ¯ P ′ . The higher PSNR is, the more image quality of the noisy face image retains, i e . ., the smaller the difference with the original face image.

## VI. THE OPTIMIZED RDP ALGORITHM

Algorithm 1 summarizes the optimized RDP: the inputs include the original face image pixel matrix P , the total privacy budget ε 0 , and the geometric superposition parameter p ∈ [0 , 1] ; and the outputs are the noise sanitized face image pixel matrix P ′ and the noisy feature vector ¯ F ′ .

First, we obtain the original sparse feature coefficient matrix C by multiscale feature dimensionality reduction of the original face image, which can be cast into the one-dimensional feature coefficient vector ¯ C . Then, the influence or weight w ik the k -th element ¯ C k ∗′ of the sorted noisy wavelet coefficient vector ¯ C ∗′ , with respect to the i -th element ¯ F ′ i of the noisy feature vector ¯ F ′ , is calculated, with 1 ≤ i ≤ M F , 1 ≤ k ≤ M P . After that, we obtain the optimized scale parameter b k of the Laplacian noise for the k -th element ¯ C k ∗′ of the sorted noisy feature coefficient vector ¯ C ∗′ according to the NA or the LMGD, with 1 ≤ k ≤ M P . Finally, the sorted noisy feature coefficient vector ¯ C ∗′ is obtained through the Definition 4, from which we obtain the noisy feature coefficient matrix C ′ and the noisy feature vector ¯ F ′ of the noisy face image generated from C ′ .

Now let's look at the time complexity and storage complexity of Algorithm 1, which is mainly determined by two steps: 1) the step of calculating the weights between the noisy feature vector ¯ F ′ and the sorted noisy feature coefficient vector ¯ C ∗′ ; and 2) the step of obtaining the optimized Laplacian noise scale parameters of ¯ C ∗′ .

Step 1) calculating the weights between the noisy feature vector ¯ F ′ and the sorted noisy feature coefficient vector ¯ C ∗′ (lines 4-8): because ¯ F ′ has M F element and ¯ C ∗′ has M P elements, the number of weights is M F × M P . Then the time complexity is O ( M M F P ) and the storage complexity of this stage is O ( M M F P ) .

Step 2) obtaining the optimized Laplacian noise scale parameters of ¯ C ∗′ (lines 9-11): the number of elements in ¯ C ∗′ is M P . When the RDP is optimized by NA, the temporal complexity is O (1) and the spatial complexity is O ( M P ) . When the optimization is LMGD, the time complexity is O ( M P ) and the spatial complexity is O ( M P ) .

Since these two steps are sequentially composed, the time complexity is O ( M M F P ) and the spatial complexity is O ( M M F P + M P ) if the scale parameters are optimized by NA; and the time complexity is O ( M M F P + M P ) and the spatial complexity is O ( M M F P + M P ) if the scale parameters are optimized by the LMGD.

## VII. EXPERIMENTAL EVALUATIONS

We implement our simulations with Python 3.10 on a laptop with Intel Core i7-13700KF, 3.40GHz and Windows 10 system equipped with 64GB main memory.

## Algorithm 1 The Optimized RDP Algorithm

Input: The original face image pixel matrix P , the total privacy budget ε 0 , the probability parameters p

Output: The noisy face image pixel matrix P ′ , the noisy feature vector ¯ F ′

- 1: Obtain C from P through multiscale feature dimensionality reduction;
- 2: C is cast into ¯ C ;
- 3: sort ¯ C to get ¯ C ∗ according to the amplitudes of its elements;
- 4: for i = 1 → M F do
- 5: for k = 1 → M P do ¯ ′

6:

w

∂F

ik

=

¯

i

∂C

k

∗′

;

- 7: end for
- 8: end for
- 9: for k = 1 → M P do

10:

Calculate b k according to the NA or the LMGD;

- 11: end for
- 12: Sample the geometric number K ∼ Geo p ( ) ;
- ¯
- 13: Initialize C ;
- ∗′
- 14: for k = 1 → M P do
- 15: if k ≤ K then
- 16: Sample the Laplacian noise ξ k ∼ Lap (0 , b k ) ; 17: ¯ C k ∗′ = ¯ C ∗ k + ξ k ;
- 18: else

$$1 9 \colon \quad \tilde { C } _ { k } ^ { \ast \prime } = \tilde { C } _ { k } ^ { \ast } ;$$

- 20: end if
- 21: end for
- 22: ¯ C ∗′ is cast into C ′ ;
- 23: Obtain P ′ from C ′ through MIWT;
- 24: Obtain ¯ F ′ according to P ′ ;
- 25: return P ′ , ¯ F ′ ;

## A. Experimental Setup

1) Datasets: We use the following datasets to conduct experiments:

Labeled Faces in the Wild (LFW) [34] contains more than 13000 images of 5749 identities and provides a standard benchmark for face verification, composed of 6000 face pairs with 3000 matched and the other 3000 non-matched.

PubFig83 [35] contains 8300 cropped face images, made up of 100 images for each of 83 public figures.

2) Evaluation Metrics: We use the following metrics in the experiment for performance comparisons:

Variance . This metric measures the visualization quality of the noise sanitized face images. Because noise will have a negative impact on the quality of the face image, the smaller the variance of the face image's noise, the higher the retained quality of the face image. In this paper, we compare the theoretical variance and the true variance of the noisy image of five methods over the given range from ε 0 = 0 2 . to ε 0 = 1 0 . of the total privacy budget with probability parameter p = 0 02 . . Since the noise follows the Laplace distribution, the theoretical variance σ t 2 = 2 K ∑ k =1 ( b k ) 2 , where b k is the scale parameter of

Laplacian noise of the k -th element ¯ C k ∗′ of the sorted feature coefficient vector ¯ C ∗′ and K is the geometric number with mean 1 p .

The real or actual variance during the numerical experiment, denoted as σ r 2 , is given by σ r 2 = M P ∑ ( k =1 ¯ P ′ k -¯ P k ) 2 , where ¯ P ′ k is the k -th element of the noisy one-dimensional pixel vector ¯ P ′ of the face image, and ¯ P k is the k -th element of the original one-dimensional pixel vector ¯ P of the face image.

Structural Similarity Index (SSIM) [36]. We conduct SSIM between the noise sanitized face images and the original face images. It imitates the human visual system and quantifies the attributes of images from brightness, contrast and structure. It uses the mean to estimate brightness, the variance to estimate contrast, and the covariance to estimate structural similarity.

Peak Signal-to-Noise Ratio (PSNR) . This metric evaluates the visualization quality of the sanitized face images. It is calculated by Eq. (7). The higher the PSNR value, the better the visualization quality of the sanitized face images or the less the face image distortion.

False Negative Rate (FNR) . This metric measures the effect of protection against the bypassing attack. It refers to the ratio of face images that should be identified as this individual but not. The larger the FNR is, the smaller the probability that the real facial features in this face image can be obtained, and the stronger protection for the method provides.

In particular, for the accuracy of face recognition after adding noise, Ω ( A ) represents the face image dataset of the individual A , and F A is a function that outputs the face recognition results, representing the feature vector of any face image A and the standard of A whether the feature vector of the face satisfies Eq. (1). If it does, the output is 1 , i e . ., F A ( A ) = 1 , indicating that the result of face recognition is correct, otherwise the output is 0 . f ′ represents the differential privacy mechanism with noise. Then, before adding noise, the accuracy is 1, that is Pr [ F A ( A ) = 1 | A ∈ Ω ( A ) ] = 1 . Let Pr [ F A ( f ′ ( A )) = 0 | A ∈ Ω ( A ) ] = ˆ p , then

$$\cdots \omega _ { 1 } \vee \cdots \nu _ { 1 } \vee \cdots \cdots \cdots \cdots \cdots \cdots \\ \hat { p } = \prod _ { i = 1 } ^ { M _ { F } } \frac { \varepsilon _ { i } } { 2 \Delta _ { i } } \cdot \exp \left ( - \sum _ { i = 1 } ^ { M _ { F } } \varepsilon _ { i } \right ) = \prod _ { i = 1 } ^ { M _ { F } } \frac { \varepsilon _ { i } } { 2 \Delta _ { i } } \cdot \exp \left ( - \varepsilon _ { 0 } \right ) \quad \underset { \sigma _ { r } ^ { 2 } } { \text{dev} } \\ \leq & \left ( \frac { \varepsilon _ { 0 } } { 2 } \right ) ^ { M _ { F } } \frac { \exp \left ( - \varepsilon _ { 0 } \right ) } { M _ { F } }, \quad \underset { \prod _ { i = 1 } ^ { 2 } \Delta _ { i } } { \text{shc} } \\ \text{flux} \\. \quad & ( \hat { P } ) \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \$$

where ε i ( ¯ F ) is the privacy budget of the i -th element of the noisy feature vector. Thus, we can get

$$M _ { F } \ln \left ( \varepsilon _ { 0 } \right ) - \varepsilon _ { 0 } \geq \ln \left ( \hat { p } \right ) + \sum _ { i = 1 } ^ { M _ { F } } \ln \left ( \Delta _ { i } \right ) + M _ { F } \ln 2. \quad \begin{array} { c } \mathbb { I } \\ a \\ c \end{array}$$

From that, the total privacy budget ε 0 of the noisy feature vector has the ability to represent the accuracy of face recognition, i e . ., equal total privacy budget ε 0 can be considered as equal face recognition accuracy. Therefore, we will no longer separately show the comparison of the face recognition accuracy of the five methods on the two data in this paper.

3) Evaluated Methods: In experiments, we compare the proposed method with the existing advanced methods. And the methods are depicted as follows.

DCT-DP. [19] proposed a privacy-preserving face recognition method based on frequency-domain Discrete Cosine Transform (DCT) differential privacy, i e . ., DCT-DP.

Pixel-DP. This method adds noises to all pixels of the face images. According to the geometric superposition, the Laplacian noises with the same scale parameter ¯ b P are added to the selected geometric number of pixels of the face image, and the sum of the privacy budgets of the facial feature vector elements of the noisy image is the given total privacy budget.

RDP. This method is to add Laplacian noises with the same scale parameter ¯ b W , as shown in Definition 4.

RDP-NA. This method is the optimized RDP with NA proposed in Sec. IV-C1). And its Laplacian noise with the optimized scale parameter b k is calculated by Eq. (6).

RDP-LMGD. The method is the optimized RDP through the LMGD proposed in Sec. IV-C2). It calculates the optimized Laplacian noise scale parameter b k by LMGD, with the step length or learning rate set to η ∈ [0 01 . , 0 1] . .

B. Evaluation on Privacy Budget and Utility

Fig. 2. Normalized difference between real and theoretical variances, i e . ., | σ r 2 -σ 2 t | σ 2 t .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-8691949f-a792-49c6-a0a9-76a945d3bdd7.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-234b94e1-4908-4904-a8dd-93fbb450251e.jpg)

1) Variance: Now, let's look at the normalized difference or deviation between the real variance of numerical experiment σ 2 r and the theoretical variance σ t 2 , which are shown in Figure 2(a) and Figure 2(b). As the privacy budget ε 0 increased, for LFW data, the differences of the methods other than RDP (as shown in Definition 4) do not undergo drastic changes, but fluctuate within a narrow range; while the differences of the five methods consistently fluctuate within a small value range for PubFig83 data. Specifically, RDP-NA (RDP Optimized with Normalization Approximation, as seen in Sec.IV-C1) and RDP-LMGD fluctuate slightly near the minimum value of 0 and do not exceed 0 05 . for LFW data and PubFig83 data. This demonstrates the effectiveness of our defined cost function for the privacy budget allocation optimization problem. In addition, among the normalized differences between the other three methods, RDP fluctuates greatly, but is never higher than 0 15 . , while DCT-DP and Pixel-DP will never exceed (0 05 . , 0 15) . range. Obviously, our experiment is reasonable and valid

Fig. 3. Logarithm of real variance log ( σ r 2 ) vs. ε 0 .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-6f406158-26a3-4c7d-b185-2a3375a2977e.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-246d9710-bd97-42de-b613-d453b952e5ad.jpg)

Finally, let's look at the real variance vs. the privacy budget ε 0 . Figure 3(a) and Figure 3(b) show the fitted curves of the logarithm log ( σ r 2 ) of the real variance σ r 2 of the five methods for the above two datasets under the given privacy budget. When the privacy budget ε 0 = 0 2 . , the log ( σ r 2 ) of RDP-NA and RDP-LMGD on both data is less than 4 0 . . Meanwhile, log ( σ r 2 ) of DCT-DP is the smallest of the other three methods, around 4 75 . , yet it is significantly higher than that of RDP-NA and RDP-LMGD, which shows that the real variance σ r 2 of RDP-NA and RDP-LMGD on the two data is significantly smaller than that of other methods when ε 0 = 0 2 . . Furthermore, as the privacy budget ε 0 increases, although the variance of all methods decreases significantly, no significant changes are found in this case, and log ( σ r 2 ) of RDP-LMGD is always slightly less than RDP-NA for LFW data, meaning that the real variance of RDP-LMGD sometimes is slightly less than RDP-NA. These findings suggest that RDP-NA and RDPLMGD have less visual damage to the original face images than other methods.

2) SSIM: In this paper, we conduct simulation experiments of SSIM on two real datasets for five methods within a given privacy budget range, as shown in Table II. It can be seen that with the increase of privacy budget ε 0 , RDP-NA and RDP-LMGD always have higher SSIM than other methods. This indicates that the scale parameters of Laplacian noise optimized by NA or LMGD can make the noisy image more similar to the original face image. In particular, the SSIM of RDP-LMGD is always slightly higher than that of RDP-NA in most cases, except that SSIM of RDP-NA in LFW data is slightly higher than that of RDP-LMGD when ε 0 = 0 2 . , and SSIM of RDP-NA and RDP-LMGD are equal when ε 0 = 1 0 . on both data.

3) PSNR: Finally, let's look at the visualization quality of the sanitized face images after the sanitizing noise is added.

In this paper, we compared the PSNR of the noise sanitized face images produced by the tested methods within a given range of total privacy budgets from ε 0 = 0 2 . to ε 0 = 1 0 . with the geometric superposition of parameter p = 0 02 . . We obtained face images from two datasets within the given privacy budgets. The fitting curves of PSNR vs. face images are shown in Figure 4, from which we can observe that as the privacy budget increases, the PSNR becomes higher. As shown in Figure 4(a), the PSNR of RDP-NA and RDP-LMGD

for LFW data higher than all compared methods. In particular, they are very close. Also, as the privacy budget increases, the advantages of RDP-NA and RDP-LMGD become more pronounced. Similarly, Figure 4(b) shows the result For PubFig83 data, from which it is clear that the PSNRs of both RDP-NA and RDP-LMGD are much higher than those of compared methods. For both data at high level of noise of ε 0 = 0 2 . , the PSNRs of both RDP-NA and RDP-LMGD significantly exceeded 50 dB, while the best PSNR of the compared methods is only around 42 dB, which is ∼ 10 dB less than that of RDPNA. In sum, both RDP-NA and RDP-LMGD achieve better visualization quality of the sanitized face images.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-fd7e3dbb-f5e5-4e6e-a015-80fe85ba2fc8.jpg)

Fig. 4. Privacy Budget ε 0 vs. PSNR.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-a5f43c0d-d2d7-4413-8467-7004be568146.jpg)

Now, let's look at the visualization effect. As shown in Figure 5 and Figure 6, we show the sanitized face images of five methods obtained with a privacy budget ε 0 = 0 2 . and the geometric superposition parameter of p = 0 02 . for LFW data and FubFig83 data. It can be shown that the Pixel-DP obtains more local mosaics on images than other wavelet transform based methods, which means that adding noise to selective wavelet coefficients of large significance on the privacy budget is better than adding noise directly to image pixels. Also, compared to the other four methods, the DCT-DP method contains deviation cross the whole face image, which is due to the global effect of the cosine basis used. On the contrary, no significant distortion nor visable mosaic effect can be seen in the sanitized face images of both RDP-NA and RDP-LMGD, which means better data utility of visualization quality of the sanitized face images has been achieved.

## C. Evaluation of Bypass Attacks Against

We use the open-source facial recognition library facerecognition to calculate the Euclidean distance of the encoding vectors of the noisy facial image and the original face image. If the Euclidean distance is less than the set threshold 0 975 . , it is considered that the noisy facial image belongs to the certain individual to whom the original facial image corresponds. Otherwise, we consider it not the same face as the original facial image.

Table III shows the false negative rate of the face recognition corresponding to the five methods under the same utility settings on LFW data.

(a) Original face image

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-96dba0fe-2f33-4426-8761-a2ae06b5b5bb.jpg)

TABLE II

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-182bf09f-32a2-40c0-8e6e-8a4d849cd37c.jpg)

## PRIVACY BUDGET ε 0 VS. SSIM

| Datasets   | Methods             | ε 0 = 0 . 2                 | ε 0 = 0 . 4                         | ε 0 = 0 . 6          | ε 0 = 0 . 8                 | ε 0 = 1 . 0             |
|------------|---------------------|-----------------------------|-------------------------------------|----------------------|-----------------------------|-------------------------|
| LFW        | DCT-DP Pixel-DP RDP | 0 . 9767 0 . 0 . 0 . 0.9955 | 0 . 9908 0 . 9918 0 . 9926 0 . 9980 | 0 . 9937 0 . 0 . 0 . | 0 . 9956 0 . 0 . 0 . 0.9993 | 0 . 9965 0 . 0 . 0.9994 |
| LFW        |                     | 9826                        |                                     | 9953                 | 9962                        | 9972                    |
| LFW        |                     | 9844                        |                                     | 9953                 | 9963                        | 9973                    |
| LFW        | RDP-NA              | 9954                        |                                     | 9988                 | 9992                        |                         |
| LFW        | RDP-LMGD            |                             | 0.9983                              | 0.9991               |                             | 0.9994                  |
| PubFig83   | DCT-DP              | 0 . 9728                    | 0 . 9882                            | 0 . 9931             | 0 . 9947                    | 0 . 9958                |
| PubFig83   | Pixel-DP            | 0 . 9836                    | 0 . 9923                            | 0 . 9951             | 0 . 9964                    | 0 . 9971                |
| PubFig83   | RDP                 | 0 . 9815                    | 0 . 9905                            | 0 . 9940             | 0 . 9960                    | 0 . 9969                |
| PubFig83   | RDP-NA              | 0.9968                      | 0 . 9985                            | 0 . 9989             | 0 . 9992                    | 0.9994                  |
| PubFig83   | RDP-LMGD            | 0 . 9967                    | 0.9986                              | 0.9991               | 0.9993                      | 0.9994                  |

TABLE III FALSE NEGATIVE RATE UNDER THE SAME UTILTY SETTING

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-3b605993-843d-427f-8a09-578e429efc86.jpg)

(b) Dct-DP

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-48654c56-e526-4715-84a7-3e6a7e8329fb.jpg)

(c) Pixel-DP

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-09c4ec10-5dd2-4ec9-8be2-616c202e788e.jpg)

(d) RDP

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-6c45737f-d0e8-40ed-9909-683056dc5f0f.jpg)

(e) RDP-NA

(f) RDP-LMGD

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-5e6c4484-1681-47f6-af24-94dfbf0e8236.jpg)

Fig. 5. Visualization of LFW Data.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-b1ef16ad-435f-4097-81c3-af6d767fff2d.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-1aba5d3e-2201-4b95-b36d-3803212fd2c5.jpg)

(a) Original face image

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-67481a07-c023-4c1e-bd32-756648e1d41b.jpg)

(b) Dct-DP

(c) Pixel-DP

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-6c4f45f0-c6be-421d-9452-e47fbb5d8ea7.jpg)

(d) RDP

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00294v1-figure-517ff5c4-a8ab-479c-a960-61350c7d495d.jpg)

(e) RDP-NA

(f) RDP-LMGD

Fig. 6. Visualization of PubFig83 Data.

| Dataset   | Method                     | False Negative Rate                        |
|-----------|----------------------------|--------------------------------------------|
| LFW       | DCT-DP Pixel-DP RDP RDP-NA | 37 . 75% 39 . 44% 40 . 05% 77 . 18% 81.38% |
| LFW       | RDP-LMGD                   |                                            |

## VIII. CONCLUSION

In this paper, a novel and efficient RDP is proposed to protect the facial features privacy for face image data, which can be formulated as the constraint optimization problem of maximizing the data utility or visualization quality of the noise sanitized face images, for a given privacy budget ε 0 . First, the RDP calculates the influence or weights of the noisy feature coefficients in the multiscale transform subspaces with respect to the total privacy budget ε 0 ; then, the RDP adds Laplacian noises to the weight-ranked feature coefficients by geometric superposition, followed by the inverse transform to obtain the noise sanitized face image. It is rigorously proved that the noisy facial feature vector obtained with RDP satisfies ε 0 -differential privacy. After that, the constraint optimization problem is formulated by the LM method by making the first-order Taylor expansion approximation on the influence or weight of the feature coefficients on the facial features. Furthermore, to solve the nonlinear LM, two methods are proposed, one for the real-time online applications and the other for the accurate offline applications: 1) the analytical NA method for the real-time online applications: by assuming that all Laplacian noise scale parameters to be identical to their average value; and 2) the LMGD method for the accurate offline applications: by computing the optimal scale parameters of Laplacian noise through iterative updating. Experiments on two real-world datasets are performed, which show that the RDP does achieve better data utility for a given facial feature privacy budget ε 0 , compared to other state-of-the-art DP methods.

## REFERENCES

- [1] Z. Kou, L. Shang, Y. Zhang, S. Duan, and D. Wang, 'Can i only share my eyes? a web crowdsourcing based face partition approach towards privacy-aware face recognition,' in Proceedings of the ACM Web Conference 2022 , ser. WWW '22. New York, NY, USA: Association for Computing Machinery, 2022, p. 3611-3622. [Online]. Available: https://doi.org/10.1145/3485447.3512256
- [2] B. Durmaz and E. Ayday, 'Entering watch dogs: Evaluating privacy risks against large-scale facial search and data collection,' in Proceedings of IEEE INFOCOM 2021 - IEEE Conference on Computer Communications Workshops , 2021, pp. 1-6.
- [3] H. Wei, 'Alipay upgrades facial-recognition system,' China Daily, 12 2018, accessed: 2024-06-03. [Online]. Available: https://www.chinadaily.com.cn/a/201812/14/ WS5c12f272a310eff303290f11.html
- [4] X. Zhang, H. Ye, Z. Huang, X. Ye, Y. Cao, Y. Zhang, and M. Yang, 'Understanding the (in)security of cross-side face verification systems in mobile apps: A system perspective,' in Proceedings of 2023 IEEE Symposium on Security and Privacy (SP) , 2023, pp. 934-950.
- [5] Y. Ma, J. Shen, Z. Zhao, H. Liang, Y. Tan, Z. Liu, K. Qian, M. Yang, and B. Hu, 'What can facial movements reveal? depression recognition and analysis based on optical flow using bayesian networks,' IEEE Transactions on Neural Systems and Rehabilitation Engineering , vol. 31, pp. 3459-3468, 2023.
- [6] I. Avcibas, 'Morphed face detection with wavelet-based co-occurrence matrices,' IEEE Signal Processing Letters , vol. 31, pp. 1344-1348, 2024.
- [7] H. Dastmalchi and H. Aghaeinia, 'Super-resolution of very low-resolution face images with a wavelet integrated, identity preserving, adversarial network,' Signal Processing: Image Communication , vol. 107, p. 116755, 2022. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0923596522000753
- [8] S. M. Eragi, F. Bensaid, and A. M. Alimi, 'Efficient human face recognition in real-life applications using the discrete wavelet transformation (hfrdwt),' Multimedia Tools And Applications , 2023.
- [9] W. Xu, J. Liu, S. Zhang, Y. Zheng, F. Lin, J. Han, F. Xiao, and K. Ren, 'Rface: Anti-spoofing facial authentication using cots rfid,' in Proceedings of IEEE INFOCOM 2021 - IEEE Conference on Computer Communications , 2021, pp. 1-10.
- [10] G. Brown, J. Martinez-del Rincon, and P. Miller, 'Least privilege learning for attribute obfuscation,' in Pattern Recognition . Cham: Springer International Publishing, 2022, pp. 142-156.
- [11] L. Ou, Y. He, S. Liao, Z. Qin, Y. Hong, D. Zhang, and X. Jia, 'Faceidp: Face identification differential privacy via dictionary learning neural networks,' IEEE Access , vol. 11, pp. 31 829-31 841, 2023.
- [12] Z. Wang, H. Wang, S. Jin, W. Zhang, J. Hut, Y. Wang, P. Sun, W. Yuan, K. Liu, and K. Rent, 'Privacy-preserving adversarial facial features,' in Proceedings of 2023 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , 2023, pp. 8212-8221.
- [13] V. Ravindra and A. Grama, 'De-anonymization attacks on neuroimaging datasets,' in Proceedings of the 2021 International Conference on Management of Data , ser. SIGMOD '21. New York, NY, USA: Association for Computing Machinery, 2021, p. 2394-2398. [Online]. Available: https://doi.org/10.1145/3448016.3457234
- [14] X. Kou, Z. Zhang, Y. Zhang, and L. Li, 'Efficient and privacypreserving distributed face recognition scheme via facenet,' in ACM Turing Award Celebration Conference-China , ser. ACM TURC 2021. New York, NY, USA: Association for Computing Machinery, 2021, p. 110-115. [Online]. Available: https://doi.org/10.1145/3472634.3472661
- [15] X. Wang and Z. Leng, 'Image encryption algorithm based on face recognition, facial features recognition and bitonic sequence,' Multimedia Tools and Applications , vol. 83, no. 11, pp. 31 603-31 627, 2024. [Online]. Available: https://doi.org/10.1007/s11042-023-16787-8
- [16] J. Yang, J. Liu, R. Han, and J. Wu, 'Transferable face image privacy protection based on federated learning and ensemble models,' Complex &amp; Intelligent Systems , vol. 7, pp. 2299-2315, 2021.
- [17] M. Chamikara, P. Bertok, I. Khalil, D. Liu, and S. Camtepe, 'Privacy preserving face recognition utilizing differential privacy,' Computers &amp; Security , vol. 97, p. 101951, 2020. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S0167404820302273
- [18] C. Dwork and A. Roth, 'The algorithmic foundations of differential privacy,' Foundations and Trends in Theoretical Computer Sci-
- ence , vol. 9, no. 3-4, pp. 1-27,29-117,119-141,143-191,193-219,221235,237-259,261-267,269-277,279-286,a1-a2, 2013.
- [19] J. Ji, H. Wang, Y. Huang, J. Wu, X. Xu, S. Ding, S. Zhang, L. Cao, and R. Ji, 'Privacy-preserving face recognition with learnable privacy budgets in frequency domain,' in Proceedings of Computer Vision, ECCV 2022 , vol. 13672, 2022, pp. 475-491.
- [20] Q. Meng, F. Zhou, H. Ren, T. Feng, G. Liu, and Y. Lin, 'Improving federated learning face recognition via privacy-agnostic clusters,' in IProceedings of nternational Conference on Learning Representations , 2022. [Online]. Available: https://openreview.net/forum?id=7l1IjZVddDW
- [21] M. Xue, C. He, J. Wang, and W. Liu, 'Backdoors hidden in facial features: a novel invisible backdoor attack against face recognition systems,' Peer-to-Peer Networking and Applications , vol. 14, no. 1, 2021.
- [22] F. Ding, B. Fan, Z. Shen, K. Yu, G. Srivastava, K. Dev, and S. Wan, 'Securing facial bioinformation by eliminating adversarial perturbations,' IEEE Transactions on Industrial Informatics , vol. 19, no. 5, pp. 66826691, 2023.
- [23] L. Yuan, W. Chen, X. Pu, Y. Zhang, H. Li, Y. Zhang, X. Gao, and T. Ebrahimi, 'Pro-face c: Privacy-preserving recognition of obfuscated face via feature compensation,' IEEE Transactions on Information Forensics and Security , vol. 19, pp. 4930-4944, 2024.
- [24] Y. Mi, Z. Zhong, Y. Huang, J. Ji, J. Xu, J. Wang, S. Wang, S. Ding, and S. Zhou, 'Privacy-preserving face recognition using trainable feature subtraction,' in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) , June 2024, pp. 297-307.
- [25] Z. Ma, Y. Liu, X. Liu, J. Ma, and K. Ren, 'Lightweight privacypreserving ensemble classification for face recognition,' IEEE Internet of Things Journal , vol. 6, no. 3, pp. 5778-5790, 2019.
- [26] W. L. Croft, J.-R. Sack, and W. Shi, 'Obfuscation of images via differential privacy: From facial images to general images,' Peer-to-Peer Networking and Applications , vol. 14, pp. 1705-1733, 2021.
- [27] L. Ou, Z. Qin, S. Liao, J. Weng, and X. Jia, 'An optimal noise mechanism for cross-correlated iot data releasing,' IEEE Transactions on Dependable and Secure Computing , vol. 18, no. 4, pp. 1528-1540, 2021.
- [28] A. A. Toda, 'Weak limit of the geometric sum of independent but not identically distributed random variables,' arXiv: Probability , 2011. [Online]. Available: https://api.semanticscholar.org/CorpusID:88512037
- [29] A. N. Akansu and R. A. Haddad, 'Multiresolution signal decomposition, transforms, subbands and wavelets,' Academic press , 2000.
- [30] M. Turk and A. Pentland, 'Eigenfaces for recognition,' Journal of Cognitive Neuroscience , vol. 3, no. 1, pp. 71-86, 1991.
- [31] C. K. Chui, An Introduction to Wavelets . Academic Press, 1992, vol. 1.
- [32] S. Mt Aznam and M. Chowdhury, 'Generalized haar wavelet operational matrix method for solving hyperbolic heat conduction in thin surface layers,' Results in Physics , vol. 11, pp. 243-252, 2018. [Online]. Available: https://www.sciencedirect.com/science/article/pii/S2211379718314682
- [33] D. R. Bull, 'Transforms for image and video coding.' Elsevier, 2014, pp. 133-169.
- [34] G. B. Huang, M. Mattar, T. Berg, and E. Learned-Miller, 'Labeled faces in the wild: A database for studying face recognition in unconstrained environments,' in Proceedings of Workshop on Faces in 'Real-Life' Images: Detection, Alignment, and Recognition . Marseille, France: Erik Learned-Miller and Andras Ferencz and Fr´ ed´ eric Jurie, Oct. 2008. [Online]. Available: https://inria.hal.science/inria-00321923
- [35] N. Pinto, Z. Stone, T. Zickler, and D. Cox, 'Scaling up biologicallyinspired computer vision: A case study in unconstrained face recognition on facebook,' in Proceedings of CVPR 2011 Workshops , 2011, pp. 3542.
- [36] Z. Wang, A. Bovik, H. Sheikh, and E. Simoncelli, 'Image quality assessment: From error visibility to structural similarity,' IEEE Transactions on Image Processing , vol. 13, no. 4, pp. 600-612, 2004.
Lu Ou Shaolin Liao Shihui Gao Guandong Huang Zheng Qi
RDP: Ranked Differential Privacy for Facial Feature Protection in Multiscale Sparsified Subspace
With the widespread sharing of personal face images in applications' public
databases, face recognition systems faces real threat of being breached by
potential adversaries who are able to access users' face images and use them to
intrude the face recognition systems. In this paper, we propose a novel privacy
protection method in the multiscale sparsified feature subspaces to protect
sensitive facial features, by taking care of the influence or weight ranked
feature coefficients on the privacy budget, named "Ranked Differential Privacy
(RDP)". After the multiscale feature decomposition, the lightweight Laplacian
noise is added to the dimension-reduced sparsified feature coefficients
according to the geometric superposition method. Then, we rigorously prove that
the RDP satisfies Differential Privacy. After that, the nonlinear Lagrange
Multiplier (LM) method is formulated for the constraint optimization problem of
maximizing the utility of the visualization quality protected face images with
sanitizing noise, under a given facial features privacy budget. Then, two
methods are proposed to solve the nonlinear LM problem and obtain the optimal
noise scale parameters: 1) the analytical Normalization Approximation (NA)
method with identical average noise scale parameter for real-time online
applications; and 2) the LM optimization Gradient Descent (LMGD) numerical
method to obtain the nonlinear solution through iterative updating for more
accurate offline applications. Experimental results on two real-world datasets
show that our proposed RDP outperforms other state-of-the-art methods: at a
privacy budget of 0.2, the PSNR (Peak Signal-to-Noise Ratio) of the RDP is
about ~10 dB higher than (10 times as high as) the highest PSNR of all compared
methods.