Quantum Representations and Scaling Up Algorithms of Adaptive Sampled-Data in Log-Polar Coordinates

In log-polar coordinates, the conventional data sampling method is to sample uniformly in the log-polar radius and polar angle directions, which makes the sample at the fovea of the data denser than that of the peripheral. The central oversampling phenomenon of the conventional sampling method gives no more efficient information and results in computational waste. Fortunately, the adaptive sampling method is a powerful tool to solve this problem in practice, so the paper introduces it to quantum data processing. In the paper, the quantum representation model of adaptive sampled data is proposed first, in which the upper limit of the sampling number of the polar angles is related to the log-polar radius. Owing to this characteristic, its preparation process has become relatively complicated. Then, in order to demonstrate the practicality of the model given in the paper, the scaling up algorithm with an integer scaling ratio based on biarcuate interpolation and its circuit implementation of quantum adaptive sampled data is given. However, due to the special properties of the adaptive sampling method in log-polar coordinates, the interpolation process of adaptive sampled data becomes quite complicated as well. At the end of this paper, the feasibility of the algorithm is verified by a numerical example.


Introduction
Quantum computation [1], as the interdisciplinary subject of many disciplines, such as computer science, mathematics, physics and engineering technology, has generated much interest in the world. Quantum computation can solve many problems that classical computers cannot solve or are inefficient at, and has the ability for data processing beyond classical computers by utilizing the coherent superposition characteristics of quantum states. Quantum algorithms [2,3] use quantum interference, quantum entanglement and other unique quantum mechanical principles to calculate, which can achieve exponential acceleration or quadratic acceleration for some classical algorithms.
Although quantum data processing has made great progress, there are still many problems. Most of the existing quantum data representation models are carried out in Cartesian coordinates; only few studies on quantum data representation models are discussed in log-polar coordinates. Compared with Cartesian coordinates, log-polar coordinates have the advantages of rotation and scale invariance. Because any scale and rotation in Cartesian coordinates can be represented as shifting in the log-polar radius and the polar angular directions in log-polar coordinates, respectively, it is easier to rotate data in log-polar coordinates than that in Cartesian coordinates. The existing quantum data representation model in log-polar coordinates is based on sample uniformly in both the log-polar radius and polar angle directions. The main problem is the high computational cost in data processing, which comes from central oversampling. In other words, the sample density in the data center is higher than that in the periphery, and no more useful information is gained from oversampling. An effective way to solve this problem is to employ the adaptive sampling method [21], and the detailed information is given in the next section.
According to the characteristics of log-polar coordinates, the adaptive sampling method can effectively solve the problem of oversampling in the center because the upper limit of the sampling number of the polar angles is related to the log-polar radius. This greatly avoids the computational waste caused by central oversampling. Hence, the new quantum data representation model based on the adaptive sampling method in log-polar coordinates is proposed in the paper. In quantum representations of adaptive sampled data in log-polar coordinates (QRASD), three entangled qubit sequences are applied to deposit the data information, the log-polar radius information and the polar angle information. In order to demonstrate the practicality of the representation model, the quantum data scaling up algorithm with integer scaling ratio based on biarcuate interpolation is given. The feasibility of the algorithm is verified by an example at the end of the paper.
The rest of this paper is arranged as follows. In Section 2, the conventional sampling method and adaptive sampling method in log-polar coordinates are introduced. The novel representation model and its preparation procedure are described in Section 3. In Section 4, the quantum data scaling up algorithm with integer scaling ratio using biarcuate interpolation based on the QRASD model is given. Finally, a conclusion is presented in Section 5.

Preliminaries
In this section, the conventional sampling method and adaptive sampling method in log-polar coordinates are introduced in detail first, and then the QUALPI model based on the conventional sampling method is also discussed.
The conventional mathematical expression for a mapping procedure in which logpolar coordinates are interchangeable with Cartesian coordinates is shown below: where (x c , y c ) represents the center datum of the transformation in Cartesian coordinates; (x, y) denotes the sampling datum in Cartesian coordinates; (ρ, θ) denotes the log-polar radius and the polar angular position in log-polar coordinates. For the sake of simplicity, the base is the natural logarithm. Figure 1 shows a conventional comparison between the two sampling methods. Figure 1a represents the conventional sampling method in log-polar coordinates, and Figure 1b represents the conventional sampling method in Cartesian coordinates.
(a) (b)  Figure 1a indicates that the sample at the fovea of the data is denser than that of the peripheral using the conventional sampling method in log-polar coordinates, so the data processing focuses more on the fovea or central area while the peripheral area is given less consideration. Fortunately, the adaptive sampling method in log-polar coordinates [21] avoids this problem, which is introduced below.
Let n ρ and n θ be the numbers of samples in the log-polar radius and the polar angular directions, respectively; R i , i = 0, 1, . . . , n ρ , be the radius size which is sampled equally to cover 360 • in the angular direction for θ = 0, 1, . . . , n θ . Furthermore, to prevent undersampling in the log-polar radius direction, R n ρ − R n ρ −1 ≤ 1 must be satisfied. It is easy to see that the following equation should be true: In the polar angular direction, the maximum sampling circumference, R max , has approximately 2πR max data. Hence, the number of samples in the polar angular direction is equal or greater than that amount, which is necessary to prevent data being missed in the sampling procedure, that is: To make the computation simpler, n ρ and n θ are selected as a multiple power of two. Then, Equations (4) and (5) can be expressed as: where A is the nearset integer equal to or greater than A. It is found that the most suitable sampling method is to set the number of polar angles to eight times the number of log-polar radii. Therefore, the adaptive sampling method uses an eight-fold sampling method for each circumference. The adaptive sampling method is shown in Figure 2a. Figure 2b describes the resulting sample distribution in the polar angular and log-polar radius directions, a series of sample points are arranged in the step-like manner, that is, as ρ increases, the length of the sample points (number of samples represented by θ) increases. The total number of samples is less than what is needed for the conventional sampling method to have the same sampling resolution. The QUALPI model [9] based on conventional sampling method stores and processes data in log-polar coordinates. It utilizes three entangled qubit sequences to store the gray scale information, the log-polar radius position and the polar angular position information for each pixel in a log-polar image. The model is shown below: The QUALPI model realizes the quantum representation of log-polar images for the first time, but the problem of central oversampling in log-polar coordinates still exists. Therefore, it is of great significance for the development of quantum data processing in the future to improve the existing quantum data representation models and find a better method to represent quantum data in log-polar coordinates.

Quantum Data Representations and Preparation Based on Adaptive Sampling Method
The premise of the application of data in log-polar coordinates is the existence of its representation and preparation. In this section, details of representations and preparation procedures are given in Sections 3.1 and 3.2, respectively. The time complexity of the preparation procedure is also discussed in Section 3.3.

Quantum Data Representations Based on Adaptive Sampling Method
Motivated by the QUALPI model, a novel quantum representation model of adaptive sampled-data in log-polar coordinates (QRASD) to store and process data is proposed as follows: where ρ encodes the information of log-polar radius and θ encodes the information of the polar angle. For each fixed ρ, i.e., in the circumference with a log-polar radius ρ, the circumference is sampled equidistantly with sampling number of 8ρ and sampling interval of 2π/8(ρ + 1). It is worth noting that both ρ and θ are coded from 0, which have a certain impact on subsequent calculations. In order to improve the flexibility and breadth of the model, three different data information representation methods are presented.
Equation (11) is referred to as the FRQI model based on the adaptive sampling method, shorthand for the AS-FRQI sub-model.

2.
Representing the data information as in the NEQR model is taken the form of the data information in the NEQR model [8], where , C t ρθ ∈ {0, 1}, t = 0, 1, . . . , q − 1, then Equation (9) can be rewritten as: The new representation method is called the NEQR model based on an adaptive sampling method, or the AS-NEQR sub-model for short.
3. Representing the data information as in the QR2-DD model The IEEE (Institute of Electrical and Electronics Engineers) is an international organization that has designed specific binary formats for storing floating-point numbers. IEEE-754 is the most widely used standard for floating point arithmetic operations. Figure 3 shows the format of floating-point numbers in IEEE-754. The QR2-DD model [10] uses the IEEE-754 standard to represent quantum data information, which lays a foundation for quantum data processing.

Quantum Data Preparation Based on Adaptive Sampling Method
Since the preparation process for the three quantum data representation sub-models given above is similar, only one detailed preparation process is given in this subsection. Compared with the AS-FRQI sub-model and the AS-NEQR sub-model, floating-point numbers can represent decimal numbers, which is more suitable for data processing in actuality. Therefore, the AS-QR2-DD sub-model is taken as an example for quantum preparation.
Due to the inherent characteristics of the binary operation process of quantum computer, redundancy is an inevitable result in the process of binary computation in quantum computers. In Figure 5, the gray part contains data information, while the white part is redundant without data information. The empty quantum state with (p + q + 2m + 3) qubits is constructed as follows: The preparation process of adaptive sampled-data is shown in Figure 6. With the empty quantum state, the position information is first prepared in Figure 6a, and then the data information, which is the shaded part in Figure 6b, is prepared by using the comparator. Step 1: In this step, the single-qubit identity gate I and the Hadamard gate H are utilized to transform |ψ 0 to the intermediate state |ψ 1 . The whole quantum operation U 1 in Step 1 is defined as follows: which transforms the initial state |ψ 0 to the intermediate state |ψ 1 : Step 2: Compared with the preparation process of the existing quantum data representation models, the upper limit of the sampling number of θ in the QRASD model is related to ρ, that is 0 ≤ θ ≤ (8ρ + 7), so a comparison between (8ρ + 7) and θ is necessary. In order to achieve this goal, some quantum circuit modules are introduced first.

Plus one and minus one modules
This part is designed to perform plus one and minus one [17] operations. The quantum circuit of plus one and minus one modules are shown in Figure 7a,b, respectively.

Multiply by eight module
This module implements the multiplication by eight, i.e., |X → |8X , and the quantum circuit is shown in Figure 8. The essence of multiplying by eight is that the qubit shifts three places to the left. It can be seen that, the effect of the module is |000x n−1 x n−2 . . . x 0 → |x n−1 x n−2 . . . x 0 000 .

Comparator module
Comparator [22] plays an important role in quantum data processing. The function of this module is to compare |x and |y , where |x = |x n−1 x n−2 . . . x 0 , |y = |y n−1 y n−2 . . . y 0 , x t , y t ∈ {0, 1}, t = 0, 1, . . . , n − 1. The qubit |z is the output. If |z = |1 , then x ≥ y; if |z = |0 , then x < y. The quantum circuit of the comparator module is shown in Figure 9. These are all the modules required for the process of comparison. As shown in Figure 10, the specific operation of this step is as follows: The qubit representing ρ is executed first by the operation of adding one to obtain (ρ + 1). Then, (ρ + 1) is multiplied by eight and subtracted by one to obtain (8ρ + 7). Finally, (8ρ + 7) is compared with θ. When the output |z is equal to |1 , the next step is executed. Step 3: Since there are 8(ρ + 1) polar angles on each circumference of log-polar radius ρ, this step is divided into 8(1 + 2 + . . . + M) = 4M(M + 1) sub-operations to set the data information at the corresponding position. The sub-operation U ρθ is expressed as below: where the value-setting operation Ω ρθ can be shown as: It is obviously that Ω ρθ only affect data at the specified position. The process for Ω ρθ to obtain data is as follows: The following formula shows the result of the intermediate state |ψ 1 under the influence of U ρθ .
To set the 4M(M + 1) log-polar data, the whole operation U 2 which transforms the intermediate state |ψ 1 to the final state |ψ 2 is defined below: Hitherto, the initialized state |0 (p+q+2m+3) is transformed to the desired state.

Time Complexity Analysis of the Preparation Procedure
Time complexity refers to the time needed to complete an algorithm in the field of computer science and engineering, and is an important parameter to evaluate the merits of an algorithm. Therefore, it is necessary to discuss the time complexity of the preparation procedure. To calculate the time complexity, the number of operation units of the preparation procedure must be estimated, and each operation unit runs for the same amount of time.
Proof. Firstly, in Step 1, U 1 is consisted of (2m + 3) Hadamard gates, and: where A is the nearset integers equal to or greater than A. Then it's obvious that the time complexity of this step is O(m). Secondly, U 2 can be decomposed into four operations: plus one, multiply by eight, minus one and comparator operations.
The plus one module consists of one NOT gate, one C-NOT gate, one C 2 -NOT gate, . . . , and one C m−1 -NOT gate. While a Toffoli gate (C 2 -NOT gate) can be decomposed into six C-NOT gates and some 1-bit gates, such as two Hadamard gates H, one phase gate S, three π 8 gates T and four conjugate transpose of the π 8 gate T + [23]. Summed up, a Toffoli gate has a time complexity of 16. Figure 11 depicts the process of implementing the Toffoli gate with elementary quantum gates. When t = 3 and 4, one C t -NOT gate can be decomposed into (3 · 2 t − 4) C-NOT gates and 2 · 2 t 1-bit gates. When t = 5, 6, . . . , m − 1, each controlled quantum gate C t -NOT can be simulated by 8(t − 3) Toffoli gates [24]. Therefore, the time complexity of the plus one module is O(m 2 ). Similar to the plus one module, the minus one module requires one NOT gate, one C-NOT gate, one C 2 -NOT gate, . . . , and one C m+2 -NOT gate to complete. Then, the time complexity of the minus one module is also O(m 2 ).
The comparator operation consists of one C 2(m+3) -NOT gate, one C 2m+5 -NOT gate, three C 2 -NOT gates, three C 4 -NOT gates, . . . and three C 2(m+2) -NOT gates. The further decomposition has the time complexity of the comparator module O(m 2 ).
To sum up, the time complexity of Step 2 is O(m 2 ). Last but not least, the operation U 2 in Step 3 can be divided into 4M(M + 1) suboperations U ρθ to set all the data at the corresponding position. Each U ρθ is composed of (p + q) C 2m+3 -NOT gates. It is known that sub-operation U ρθ costs no more than O((p + q)m). Hence, the time complexity of Step 3 will cost no more than O((p + q)mM 2 ).
As the discussion above, the time complexity of the whole preparation procedure is about O(m 2 + (p + q)mM 2 ).
The quantum circuit of preparation for the example data shown in Figure 4 is described in Figure 12.

Quantum Data Scaling Up Algorithm with Integer Scaling Ratio Based on Biarcuate Interpolation
Data scaling up as a kind of data geometric transformation has been widely studied and applied in the classical data processing; however, the quantum version of which, in logpolar coordinates does not exist. Here, by employing biarcuate interpolation, a quantum algorithm to scale up adaptive sampled-data in log-polar coordinates is given. To make it easier to understand, the AS-NEQR sub-model is taken as an example.

Quantum Circuit Modules for Arithmetic Operations
Designing a quantum scaling up circuit of biarcuate interpolation is a complex task. A series of quantum circuit modules to realize special functions are introduced first, including multiplying the C-NOT operation module, adder module [25], subtractor module [25], multiplier module [26] and divider module [27].

Multiply C-NOT Operation Module
The multiply C-NOT operation is to copy the determined quantum states. The quantum circuit of the multiply C-NOT operation is shown in Figure 13. It utilizes n C-NOT gates to copy the n-qubit information of |X = |x n−1 x n−2 . . . x 0 into the n ancillary qubits |0 ⊗n , where |x n−1 , |x n−2 , . . . , |x 1 and |x 0 are the control qubits and the n ancillary qubits |0 ⊗n are the target qubits. The time complexity of the multiply C-NOT operation module is O(n).

Adder Module
The adder module [25] is not only the core of the multiplier, but also one of the core modules of constructing the biarcuate interpolation method. The adder (A) adding a n-qubit unsigned integer |Y to a n-qubit unsigned integer |X is designed by one half adder (HA) and (n − 1) full adders (FA) as shown in Figure 14. Here, the input |X comes out as the sequence |s n s n−1 . . . s 0 , which represents the sum of |X + Y , and the input |Y stays the same. Other unremarked qubits are the garbage outputs and the input qubit |0 is the ancillary constant input. For convenience, the block diagram of A omits the ancillary inputs and the garbage outputs. The time complexity of the adder module is O(n).

Subtractor Module
The function of the subtractor (S) is to implement the subtraction of binary unsigned integers |X and |Y [25]. The subtractor is designed by one half subtractor (HS) and (n − 1) full subtractors (FS). The input |X come out as the sequence |d n d n−1 . . . d 0 , which represents the difference of |X − Y , and the input |Y stays the same. If X ≥ Y, then |d n = |0 , if X < Y, then |d n = |1 . The quantum circuit of the subtractor module is shown in Figure 15. The time complexity of the subtractor module is O(n).

Multiplier Module
The purpose of the multiplier module [26] is to take the product of two numbers, which can be done by addition. Let |X = |x m−1 x m−2 . . . x 0 and |Y = |y n−1 y n−2 . . . y 0 , then: First, the auxiliary qubit |0 ⊗n is given as the initial value of the partial product. When |x 0 = |1 , copy |Y to the partial product, add '0' to the left of the partial product, and add '0' to the right of |Y . Then, if |x t = |1 , t = 1, 2, . . . , m − 1, add the padded |Y to the partial product. If |x t = |0 , the partial product remains unchanged. Since the adder is a carry addition, only after the proccess do we need to add '0' to the right of |Y . The block diagram of multiplier (M) is shown in Figure 16, where it consists of multiple adders. The time complexity of the multiplier module is O(mn).

Divider Module
Based on restoring the division algorithm, Khosropour et al. [27] proposed the quantum division circuit. In order to make full use of the operation modules given in this paper, appropriate modifications have been made in this subsection. In Figure 17, a pseudo code for this algorithm is given. According to the algorithm, to get the quotient (|Q ) and the remainder (|R ), first of all, the lowest n-qubit of |A are equal to the dividend |X , and the highest n-qubit are equal to zero. Furthermore, register |B is initiated such that its highest n-qubit contains the divisor |Y and its lowest n-qubit are set to zero. Then, |A is multiplied by two and compared to |B . If the result is nonnegative, then the qubit |q t in the quotient is set to one and |A should be subtracted |B ; otherwise, |q t = |0 . This process is repeated n times, each time one qubit of the quotient is determined. A concrete calculation example is given to clarify this procedure. Let |X = |1001 and |Y = |0010 , then the initial state is: The first iteration: The second iteration: The third iteration: The forth iteration: Therefore, it turns out that: Obviously, the calculation process of the divider can be achieved through the comparators, subtractors and multiply by two operations. The corresponding quantum circuit is shown in Figure 18. In order to facilitate the subsequent operations, the upper output is marked as Y, the middle output is marked as R and the lower output is marked as Q in the quantum simple circuit of the divider. The time complexity of the divider module is O(n 3 ).

Quantum Biarcuate Interpolation Method in Log-Polar Coordinates
The principle of the biarcuate interpolation method in log-polar coordinates is similar to that of the bilinear interpolation method in Cartesian coordinates. The difference is that bilinear interpolation is based on line segments, but biarcuate interpolation is based on angular segments.

Arcuate Interpolation Method
The datum at position (ρ, θ) of the interpolated data can be reconstructed from the datum at positions (ρ, θ 1 ) and (ρ, θ 1 ) of the original data. As shown in Figure 19, the value of the desired datum f (ρ, θ) can be obtained by the following equation: where f (ρ, θ 1 ) and f (ρ, θ 1 ) are two known data. That is to say, (28) Figure 19. Arcuate interpolation method.

Quantum Realization of Biarcuate Interpolation Method
Because the number of samples of polar angles on each circumference is different, the interpolation process of adaptive sampled-data becomes quite complicated. Assume that a ρ × 8ρ (ρ = 0, 1, . . . , M − 1) quantum data |D is scaled up to a ρ × 8ρ (ρ = 0, 1, . . . , M − 1) quantum data |D based on biarcuate interpolation. The scale ratio is r, i.e., M = rM, where r ∈ N and represented by k qubits. For the sake of simplicity, the interpolation process of adaptive sampled-data is divided into two steps according to the log-polar radius. The first step is interpolation based on the arcuate interpolation method. The second step is interpolation based on the biarcuate interpolation method using the data interpolated in the previous step. The specific process is discussed below.
Step 1: Circuit realization of arcuate interpolation method. For each circumference with log-polar radius ρ in the scaled up data, if (ρ + 1) is an integer multiple of r, i.e., ρ + 1 = r(ρ 0 + 1), then the interpolated data is the corresponding circumference. As shown in Figure 21, ρ 0 and θ 0 in (a) are the log-polar radius and polar angle corresponding to ρ and θ in (b), respectively. That is, The data after arcuate interpolation is shown in Figure 21c as the solid line. In this step, the data on the dotted line has not been interpolated. It can be seen from Figure 21b, there are (r − 1) points between θ and (θ + r) that need to be interpolated. In order to get the data of these points, the value of ρ 0 and θ 0 in the original data must be found first. Obviously, Thus,  (36), the value f (ρ, θ) in position (ρ, θ) of the interpolated data in Figure 21b can be calculated to arcuate interpolation according to the following equation: In order to prepare the data information f (ρ, θ) in position (ρ, θ) of the resulting data, the data information f (ρ 0 , θ 0 ) and f (ρ 0 , θ 0 + 1) in positions (ρ 0 , θ 0 ) and (ρ 0 , θ 0 + 1) of the original data need to be prepared first. To achieve the data information for the corresponding position, two quantum operators Ω ρ 0 ,θ 0 and Ω ρ 0 ,θ 0 +1 are used to achieve the original data information, respectively. The quantum operator Ω ρ,θ can realize the aim of assigning the data information as follows: Due to the quantum superposition, all data are processed simultaneously, so all original data must be assigned. Let the assignment operation θ=0 Ω ρ,θ , then the time complexity of the assignment operation is O(qmM 2 ).
The quantum circuit of arcuate interpolation is shown in Figure 22, where m =  Step 2: Circuit realization of biarcuate interpolation method. As shown in Figure 23, this step interpolates the part which is not interpolated in Step 1, that is, the data in the dotted part of Figure 21c or Figure 23a. Figure 23b shows the specific schematics of biarcuate interpolation for adaptive sampled-data, where points (ρ 1 , θ 1 ), (ρ 1 , θ 1 + 1), (ρ 2 , θ 2 ) and (ρ 2 , θ 2 + 1) in Figure 23a are the four closest points to (ρ, θ). Figure 23c describes the data after biarcuate interpolation, which is the scaled up data. It is obvious that the four interpolated data f (ρ 1 , θ 1 ), f (ρ 1 , θ 1 + 1), f (ρ 2 , θ 2 ) and f (ρ 2 , θ 2 + 1) from the previous step are needed to complete the interpolation for this step. The mathematical expression is as follows: From Step 1, (ρ 1 + 1) and (ρ 2 + 1) are both integer multiples of r, assuming that: where k = Q[ρ 1 + 1, r]. That is: Oobviously, Because there are 8(ρ + 1) polar angles on each circumference with log-polar radius ρ, and the common multiple of 8(ρ 1 + 1) and 8(ρ + 1) is 8(ρ + 1)(ρ 1 + 1). The principle of common multiples is used, there are: In the same manner, Thus, the Equation (39) can be rewritten as: As shown in Figure 24, the quantum circuit of biarcuate interpolation is given, where: m = log 2 M , M > 1; 1, M = 1, and the +(r − 1) operation can be realized by (r − 1) plus one operation. The dotted line serves no purpose but as an indicator. R = R[ρ + 1, r] represents the remainder of (ρ + 1) divided by r. When R = 0, the output f (ρ, θ) is the required data. Similarly, f (ρ, θ) has an error of less than one. However, due to the characteristics of the adaptive sampling method, there may be additional errors, which will be explained in detail in the next subsection. The time complexity of biarcuate interpolation is O(m 3 + rqm M 2 ). Therefore, the time complexity of the whole interpolation process is In classical computers, the algorithm complexity is mainly evaluated by two parameters, one is the number of rows of pre-cached data required to complete interpolation, and the other is the number of operation steps such as addition, subtraction, multiplication, division and displacement required to complete interpolation. Here, the time complexity of the proposed algorithm is estimated by 1-bit addition and subtraction as the unit of operation on classical computers. The m × n multiplication consists of O(mn) additions, and the n × n division can be decomposed into O(n 2 ) subtractions, assuming that all of the original data is known. By Equations (37) and (49), the time complexity of arcuate interpolation and biarcuate interpolation are O((m 2 + k 2 q 2 )rM 2 ) and O((k 2 + q 2 )rq 2 m 2 M 2 ), respectively. Therefore, the time complexity of the whole interpolation process in classical computers is O((k 2 + q 2 )rq 2 m 2 M 2 ). Obviously, the time complexity is much higher than that of quantum computers.
where f (2, 24) = f (2, 0) = 0. However, when ρ < r = 3, the process becomes a little more complicated. When ρ = 0, the Equation (49) turns into: where f (−1, 0) should be the data of the pole in the log-polar coordinates. However, the adaptive sampling method does not sample at the poles; thus, the characteristics of the minus one module are f (−1, 0) = f (7, 0). In this example, f (7, 0) is a redundant datum that has not been prepared, so f (7, 0) = 0. If f (7, 0) = 0, there is a certain error. The effect of all the data after biarcuate interpolation is shown in Figure 26.  Figure 4 is magnified three times.

Conclusions
In the existing quantum data representation models, the two variables representing positions are independent of each other. However, in the quantum data representation model proposed in this paper, the upper limit of the sampling number of the polar angles is related to the log-polar radius. This is an inspiration for the future research work about the representation of irregularly shaped data.
The quantum scaling up algorithm with integer scaling ratio based on biarcuate interpolation and its circuit of adaptive sampled-data in log-polar coordinates are also given in this paper. In the quantum scaling up circuit, the operation modules required are constructed on the basis of the existing methods in the references. Meanwhile, the feasibility of the quantum algorithm is also verified by an example. However, the quantum scaling up algorithm still has certain shortcomings, and the quantum scaling down circuit has not been given. That's one of the most important things we'll be looking at in the future.
In addition, how to change the scaling ratio from integer multiples to real multiples and how to optimize the circuit are also problems worth studying in future work.