72 results about "Main diagonal" patented technology

The invention relates to the technical field of electrical digital data processing, and discloses a lightweight-based multi-modal data enhancement method, which comprises the following steps: acquiring a multi-modal feature tensor; extracting a vector orthogonal projection scalar and determining the novelty; updating the global covariance matrix when the novelty is greater than a redundancy threshold, and maintaining the global covariance matrix in a register state when the novelty is not greater than the redundancy threshold; extracting a main diagonal variance component to determine a differential modulation coefficient; calculating a second-order moment manifold projection operator according to the global covariance matrix, and calibrating the operator by using a differential modulation coefficient; the calibrated operator is used for carrying out orthogonal projection on random noise to generate a structured disturbance vector, the structured disturbance vector is superposed to a multi-modal feature tensor to output enhanced features, a novelty judgment mechanism is used for restraining statistical deviation caused by steady-state redundant data, computing resource occupation is reduced, and it is ensured that semantic alignment between modals is maintained in enhanced feature distribution.

The invention discloses a document comparison method based on an improved LCS algorithm, and relates to the field of data processing, and the method comprises the steps: receiving N to-be-compared documents uploaded by a user; dividing the N documents into N groups, wherein each group comprises a main document and N-1 auxiliary documents; setting a minimum common subsequence length threshold value L0; for each group, constructing a virtual two-dimensional array, taking the character number of the main document as the row number, taking the character number of the auxiliary document as the column number, adopting a sparse matrix for storage, and only recording matched character positions; traversing the elements of the first row and the first column row by row from the elements of the last row and the last column of the virtual two-dimensional array along the main diagonal direction of the two-dimensional array, and recording the content of continuous matching characters and the starting positions and the ending positions of the continuous matching characters in the main document and the auxiliary document; in view of low comparison efficiency of tobacco purchasing files in the prior art, optimization strategies such as document grouping, multi-thread parallel comparison and clustering combination are introduced, so that the comparison efficiency is improved.

The invention discloses a method for generating an amplitude preparation circuit and a quantum state preparation method and device, and the method comprises the steps: setting a main diagonal element of a block matrix as a to-be-prepared amplitude, and setting other elements as 0, so that the block matrix is equivalent to a diagonal matrix, thereby facilitating the construction of a unitary matrix through employing the block matrix, and in the unitary matrix, the unitary matrix is not damaged; the to-be-prepared amplitude is arranged at the upper left corner of the unitary matrix, so that the amplitude preparation of the to-be-prepared quantum state can be carried out by utilizing a circuit corresponding to the unitary matrix, and the preparation of the to-be-prepared quantum state by utilizing a large-scale operator can also be understood, thereby being beneficial to more efficiently and accurately preparing the to-be-prepared quantum state.

The invention provides a confusion matrix scoring method and device based on entropy coding, and relates to the technical field of machine learning. The method comprises the following steps: acquiring an original confusion matrix; setting an element in a main diagonal line of the original confusion matrix to be zero to obtain an error matrix; binary entropy coding is carried out on terms of the error matrix in an entropy coding mode, and a coding matrix is obtained; calculating the code length of each code in the coding matrix according to the coding matrix to obtain a code length matrix; multiplying each item of the code length matrix by each item of the error matrix to obtain an average code length matrix; summing each item of the average code length matrix to obtain a total average code length; and the value of the total average code length is used as a scoring index. According to the invention, the error rate and the dispersity of error conditions can be reflected.

A method for compressing a plaintext matrix and a computing device, the plaintext matrix is obtained, the plaintext matrix is a Fourier transform matrix or an inverse Fourier transform matrix, the plaintext matrix is decomposed to obtain a plurality of sparse matrixes, the number of non-zero diagonals of each sparse matrix except a main diagonal is at most 2, and the product of the plurality of sparse matrixes is the plaintext matrix, dividing the plurality of sparse matrixes into n groups, aiming at any group of sparse matrixes, determining a product of the group of sparse matrixes as a combined matrix corresponding to the group of sparse matrixes, aiming at any combined matrix, determining an access period corresponding to the combined matrix, and respectively performing access on the combined matrix along each diagonal line according to the access period, the compression diagonals corresponding to the combined matrix are determined according to the access result, and the plaintext matrix can be compressed by utilizing the repetition rule of the diagonal elements of the combined matrix obtained after the plaintext matrix is decomposed and combined, so that the space required for storing the plaintext matrix can be reduced, and the efficiency of transmitting the plaintext matrix can be improved.

A computer-implemented method is for reducing computing time when determining an inverse matrix from a symmetrical input matrix comprising n rows and n columns. The method includes forming i blocks with mi rows and mi columns, with mi less than n and i being greater than or equal to 2. The i blocks are each formed from entries of main diagonals which are most strongly correlated and mi entries of secondary diagonals in the same rows and columns as the entries of the main diagonals. The method further includes generating a block diagonal matrix from the i blocks, and generating the inverse matrix to the input matrix by inverting the blocks in the block diagonal matrix.

Matrix multiplication operations can be implemented, at least in part, on one or more tensor cores of a parallel processing unit. An efficiency of the matrix multiplication operations can be improved in cases where one of the input operands or the output operand of the matrix multiplication operation is a square matrix having a triangular data pattern. In such cases, the number of computations performed by the tensor cores of the parallel processing unit can be reduced by dropping computations and / or masking out elements of the square matrix input operand on one side of the main diagonal of the square matrix. In other cases where the output operand exhibits the triangular data pattern, computations can be dropped or masked out for the invalid side of the main diagonal of the square matrix. In an embodiment, a library implementing the matrix multiplication operations is provided.

The application discloses a multiple wave prediction method and device, a computing device and a storage medium. The method comprises the following steps: acquiring frequency domain data of any shot point at any seismic trace; generating a seismic data matrix according to the frequency domain data, the seismic data matrix being an upper triangular matrix with a main diagonal line of 0, and any non-first column of the seismic data matrix containing frequency domain data of at least one seismic trace of the same shot point; constructing a multiple wave matrix, the multiple wave matrix being an upper triangular matrix with a main diagonal line of 0, and any non-first column of the multiple wave matrix containing a multiple wave model of at least one seismic trace of the same shot point; establishing a functional relationship between the multiple wave matrix and the seismic data matrix; and obtaining model frequency domain data corresponding to any multiple wave model in the multiple wave matrix according to the functional relationship. According to the scheme, the prediction efficiency of the multiple wave can be improved under the premise of ensuring the prediction accuracy of the multiple wave.

Disclosed is an operation method for executing an addition operation on quantized matrices. This operation method comprises the steps in which: an NPU matrix multiplication unit executes a matrix multiplication operation on a first matrix and a second matrix, wherein the first matrix has a size of R0*2C0 and is obtained by concatenating a first input tensor and a second input tensor, each of which is composed of quantized values and has a size of R0*C0, and the second matrix is obtained by concatenating a first diagonal matrix having a size of C0*C0 and having all main diagonal element values equal to integer N1 and a second diagonal matrix having a size of C0*C0 and having all main diagonal element values equal to integer N2; and an NPU vector processing unit calculates elements of a non-quantized tensor obtainable by adding the second input tensor to the first input tensor, by adding real number B to a product of real number R and a third matrix calculated by the matrix multiplication operation.

The present application discloses a document comparison method based on an improved LCS algorithm, which relates to the field of data processing, including: receiving N documents to be compared uploaded by a user; dividing the N documents into N groups, each group comprising a main document and N-1 sub-documents; setting a minimum common subsequence length threshold L0; for each group, constructing a virtual two-dimensional array, with the number of characters in the main document as the number of rows, the number of characters in the sub-document as the number of columns, using sparse matrix storage, and only recording the matching character positions; starting from the elements in the last row and last column of the virtual two-dimensional array, traversing row by row to the elements in the first row and first column along the main diagonal direction of the two-dimensional array, recording the content of continuous matching characters, the starting position and the ending position in the main document and the sub-document; in view of the low efficiency of tobacco procurement document comparison in the prior art, the present application introduces optimization strategies such as document grouping, multi-threaded parallel comparison, cluster merging, etc., to improve the comparison efficiency.

The invention provides an internal force calculation method and system for a door-shaped tower multi-nominal-height tower leg component for realizing a digital power grid, and relates to the technical field of digital power grids. According to the method, basic parameter information of the door-shaped tower is firstly obtained, then load influence factors under different non-uniform icing working conditions are considered, and the load design value at the line hanging point of the door-shaped tower is calculated. The load design value is used for solving the main material internal force and the inclined material internal force of the tower leg under each door-shaped tower subtower nominal height combination, and the most unfavorable internal force is determined; based on the most unfavorable internal force of the shortest and longest tower leg main diagonal members, the method further utilizes an interpolation method to calculate the main diagonal member internal force of each tower leg. Based on the method, the requirements of tower legs with different heights can be met, and the design flexibility and adaptability are improved. In addition, the method is beneficial to reducing potential safety hazards caused by improper design, reducing operation risks of a power grid system, promoting digital application of power grid design and promoting the power grid to develop in a more intelligent and safer direction.

The embodiment of the invention provides a multidirectional seismic attribute tensor fusion method and device, equipment and a medium, and relates to the technical field of oil-gas exploration and development. The method comprises the following steps: extracting seismic attribute data volumes in a plurality of observation orientations from five-dimensional seismic data, and respectively calculating a seismic attribute tensor corresponding to the seismic attribute data volume in each observation orientation, each seismic attribute tensor represents the change rate of the gradient of the seismic attribute in the corresponding observation orientation along three orthogonal directions of a three-dimensional space; adding the values of the elements with the same element positions in each seismic attribute tensor to obtain a multi-azimuth seismic attribute tensor fusion operator; and calculating the sum of main diagonal elements in the multidirectional seismic attribute tensor fusion operator to obtain a multidirectional seismic attribute tensor fusion result. Thus, according to the obtained fusion result, the influence of the included angle between the observation coordinate system and the constitutive coordinate system is eliminated, and effective fusion of multi-azimuth seismic attribute information is achieved.

A gridless coherent signal direction-of-arrival estimation method based on alternating projection according to the present invention is characterized in that: the patented method uses a hydrophone array installed on an underwater vehicle to receive signal data to construct a positive semi-definite matrix S0, and uses the projection method to alternately project S0 onto a rank-constrained set and a Toeplitz matrix with averaged main diagonal elements. Secondly, the convergence degree of S0 is judged, and then T i (u) is subjected to reduced-dimensional space smoothing-singular value decomposition. Finally, z is obtained by using the polynomial root-finding method, and the incident angle of the coherent signal is solved by using the corresponding mapping relationship. The patented method effectively realizes the angle estimation of coherent signals. Secondly, the polynomial root-finding method is used instead of spectral peak search, which effectively reduces the computational amount, improves the robustness and real-time performance of the algorithm, and has high estimation accuracy, having important engineering application value.

The present invention relates to the technical field of radar signal, and specifically provides a radar interference identification method, system, terminal and storage medium based on bispectrum phase angle, including: acquiring observation data received by a radar; calculating a bispectrum phase angle matrix of the observation data; converting the bispectrum phase angle matrix into a gray-level co-occurrence matrix, and extracting image texture features based on the gray-level co-occurrence matrix; extracting the numerical values of the main diagonal from the bispectrum phase angle matrix to obtain a diagonal sequence, and performing spectrum analysis and sequence trend feature analysis based on the diagonal sequence; performing principal component analysis based on the image texture features and sequence trend features, and inputting the features after principal component analysis into a preset classification model to obtain the result of identifying a true target or a false target. The present invention improves the reliability of identifying a true target and an interference signal, and solves the problem that it is difficult to distinguish a towed interference signal by traditional methods.

The present application discloses a sparse matrix storage access method, device, equipment and medium for fluid computing, which relates to the field of fluid mechanics, including: based on the grid units in the current fluid mechanics scene, discretely solving the flow control equation to obtain a coefficient matrix; the coefficient matrix is ​​a sparse matrix; based on the mapping relationship between the connection surface between the grid units and the matrix position storage information of the non-diagonal non-zero elements in the coefficient matrix, a target mapping table is determined, and a matrix storage operation is triggered based on a preset matrix storage rule, a target mapping table and the main diagonal elements and non-diagonal non-zero elements in the coefficient matrix; after completing the matrix storage, the matrix position index corresponding to the connection surface to be accessed is determined based on the corresponding target surface number in the target mapping table, and the target non-diagonal non-zero element value is obtained using the matrix position index to complete the connection surface data access operation. The present application solves the problem of large waste of memory space in sparse matrix storage and poor random data access performance.

The invention discloses a matrix simplification method applied to electromagnetic transient simulation and a related device, and the method comprises the steps: converting a quantified component admittance matrix of a power system into a graph structure, obtaining a graph structure matrix, enabling nodes of the graph structure matrix to represent matrix rows or columns, and enabling edges of the graph structure matrix to represent a non-zero element connection relation; performing augmented path traversal search in the graph structure matrix based on a DFS algorithm, and recording a non-zero line sequence number and a non-zero column sequence number; reordering the graph structure matrix according to the non-zero row sequence number and the non-zero column sequence number to obtain a non-zero diagonal matrix; and performing path traversal search on the non-zero diagonal matrix along the main diagonal non-zero element nodes based on a DFS algorithm, and performing non-zero row and column reordering operation to obtain a triangular block structure matrix. The technical problem of low electromagnetic transient simulation efficiency caused by the problems of high calculation complexity, lack of stability of numerical values, increase of matrix understanding difficulty, time consumption of calculation and the like in the prior art can be solved.

This invention belongs to the field of sparse matrix storage and distributed parallel computing technology. It discloses an adaptive storage and parallel computing method for sparse matrices based on vertex cutting and two-stage compression. The method includes: initializing a cross-linked list in memory for dynamic assembly of sparse matrix elements; reducing the dimensionality of the dynamically assembled cross-linked list data structure to a one-dimensional continuous compressed sparse row format using a two-stage traversal algorithm; constructing a heuristic diagnostic model based on row non-zero variance and diagonal occupancy to achieve adaptive matrix storage; extracting the pure graph topology after removing self-loop elements on the main diagonal of the sparse matrix, and performing multi-way recursive graph segmentation based on vertex cutting on the pure graph topology to generate a distributed mapping array; and constructing a communication scheduling table based on the distributed mapping array to achieve deep overlap between non-blocking data exchange and local matrix computation. This invention significantly improves the storage compression rate and parallel computing scalability of large sparse matrices.

This invention relates to the field of digital signal processing technology, and discloses a DSP multi-core parallel QR decomposition method and apparatus based on tile algorithm. The method includes: dividing the matrix to be decomposed into blocks in the horizontal direction according to the number of cores of the DSP processor and setting tile blocks; performing QR decomposition on the main diagonal tile blocks; and updating all tile blocks to the right of the main diagonal tile blocks in the horizontal direction; pairwise coupling of the main diagonal tile blocks with each tile block in the vertical direction of the same column; performing QR decomposition on the coupling matrix formed by each pair of couplings; and updating all coupling matrices to the right of the pair of coupling matrices in the horizontal direction; and cyclically executing the decomposition and horizontal update of the main diagonal tile blocks, the coupling decomposition of the tile blocks in the same column, and the update of the right coupling matrix until all main diagonal tile blocks in the matrix to be decomposed have been processed, resulting in an upper triangular matrix. This invention can fully utilize the parallel processing advantages of multi-core processors to achieve inter-core parallel data processing.

The invention relates to the technical field of radar signals, and particularly provides a radar interference identification method and system based on a bispectrum phase angle, a terminal and a storage medium, and the method comprises the steps: obtaining observation data received by a radar; calculating a bispectrum phase angle matrix of the observation data; converting the bispectrum phase angle matrix into a gray level co-occurrence matrix, and extracting image texture features based on the gray level co-occurrence matrix; extracting numerical values of main diagonals from the bispectrum phase angle matrix to obtain a diagonal sequence, and performing spectrum analysis and sequence trend characteristic analysis based on the diagonal sequence; and performing principal component analysis based on the image texture features and the sequence trend features, and inputting the features after the principal component analysis into a preset classification model to obtain a result of the identified real target or false target. According to the method, the reliability of identifying the real target and the interference signal is improved, and the problem that the dragging interference signal is difficult to distinguish by a traditional method is solved.

The present invention belongs to the field of image transmission technology, and specifically relates to an image transmission method for low-power monitoring equipment, the method comprising: determining a standard order according to the definition requirement of the monitoring video collected by the monitoring equipment, adjusting the standard order according to the ratio of the actual power of the monitoring equipment to the estimated power, obtaining the adjusted order, evenly dividing all kinds of gray values ​​into kinds of gray levels, constructing an order gray co-occurrence matrix of the image according to the gray level pairs of all pixel pairs in any frame image of the monitoring video, selecting multiple elements from the main diagonal of the gray co-occurrence matrix and its vicinity as the first element, and recording the remaining elements as the second element, respectively, compressing and encoding all first pixel pairs and second pixel pairs according to the code words of the first element and the second element, obtaining the compressed encoding result of the image and transmitting it. The present invention reduces the energy consumption of the monitoring equipment when storing and transmitting images.

The application discloses a frequency coupling grid-connected inverter stability improvement method, which introduces the q-axis voltage and the q-axis current after the suppression link into the q-axis control in the form of positive feedback as the reference value of the q-axis current through the q-axis voltage loop compensation mode. The strategy can ingeniously utilize the sign difference of the frequency domain convolution itself in the positive sequence and the negative sequence, realize the increase of the amplitude of the main diagonal line subsystem in the admittance matrix of the grid-connected inverter, and greatly reduce the amplitude of the coupling subsystem, so that the coupling subsystem is negligible compared with the main diagonal line subsystem, and the frequency coupling oscillation is essentially suppressed and the stability of the system is improved.

The invention discloses a sparse matrix storage access method and device for fluid calculation, equipment and a medium, and relates to the field of fluid mechanics, and the method comprises the steps: carrying out the discrete solving of a flow control equation based on a grid unit in a current fluid mechanics scene, and obtaining a coefficient matrix; the coefficient matrix is a sparse matrix; determining a target mapping table based on a mapping relation between a connection surface between the grid units and matrix position storage information of the off-diagonal non-zero elements in the coefficient matrix, and triggering a matrix storage operation based on a preset matrix storage rule, the target mapping table and main diagonal elements and the off-diagonal non-zero elements in the coefficient matrix; and after matrix storage is completed, determining a matrix position index corresponding to the to-be-accessed connection surface based on a corresponding target surface number in the target mapping table, and obtaining a target off-diagonal non-zero element value by using the matrix position index so as to complete connection surface data access operation. According to the method, the problems of large memory space waste and poor random data access performance of sparse matrix storage are solved.

The application discloses a running-in state identification method based on a friction signal threshold-free recursive analysis result, extracts nonlinear components of a friction signal time sequence output in a friction and wear process by collecting the friction signal time sequence; constructs a threshold-free recursive matrix by combining a sliding window and an avoidance embedding phase space reconstruction method and realizes visualization of the threshold-free recursive matrix in a two-dimensional plane; divides a lower triangular matrix of the threshold-free recursive matrix into several regions by a group of equidistant line segments parallel to a main diagonal line of the recursive matrix, extracts mean values of recursive elements in the regions, analyzes a change rate of element values of the threshold-free recursive matrix when the threshold-free recursive matrix transitions along the main diagonal line to edges and corners as a recursive measure, analyzes variation rules of the recursive measure along with a running-in and wear process, and further realizes running-in state online identification driven by the friction signal threshold-free recursive analysis result. The application optimizes a recursive characteristic representation and a measure extraction method of the friction signal time sequence, and realizes running-in state online identification from a perspective of dynamics-friction coupling.

The present invention provides a method for constructing a discontinuous finite element monomial basis function applicable to any grid, belongs to the field of computer-aided engineering, and is used for hybrid grid calculations. The present invention divides the calculation area into a plurality of unstructured grids to maximize the filling of the calculation area. The expressions of the basis functions are consistent for different units and can be directly used for the discontinuous finite element method on the unstructured hybrid grid. The diagonal part of the mass matrix of the basis function space is 1, and the absolute value of the non-diagonal part is less than 1, which satisfies the main diagonal dominance and makes the calculation more stable; the basis function is defined in global coordinates, and when calculating the gradient, it is not necessary to use the Jacobian matrix of the unit isoparametric transformation, which is more direct and efficient.

The application discloses a multifunctional radar wave element extraction method and system based on a dot matrix model, and the method comprises the following steps: for a multifunctional radar (MFR) pulse sequence, a dot matrix (D) based on a pulse repetition interval (PRI) value, a pulse width (PW) value and a radio frequency (RF) value is respectively constructed according to the threshold values of the corresponding PRI value, PW value and RF value of the MFR pulse sequence; PRI ,D RF ,D PW ; the constructed dot matrix (D PRI ,D RF ,D PW is respectively subjected to image preprocessing, image hole filling and image closing operation; an L-shaped corner point on a main diagonal line of the preprocessed dot matrix is detected by using a Harris corner point detection operator, that is, a wave element starting point; and a union set of the detection results is obtained, that is, a wave element starting pulse set. The application does not depend on cluster parameter setting, can fully utilize the overall modulation law of the wave element, reduces the false alarm probability while improving the wave element detection probability, and still has good performance in the case that the pulse is lost or the parameter measurement error is large.

An electronic device according to an embodiment disclosed herein includes: a battery module including a plurality of cells; the sensor circuit is used for obtaining a state value of each battery cell in the plurality of battery cells; and a processor, where the processor is capable of: obtaining, via the sensor circuitry, a value related to a state of each of the plurality of cells as input data; obtaining feature data based on a feature matrix derived from the input data by using a specified matrix decomposition algorithm; and identifying whether each of the plurality of cells is abnormal based on the feature data. The feature matrix may be a matrix in which values of a specified number of main diagonal components among a plurality of diagonal components of a diagonal matrix for input data are replaced with 0.

A remote sensing image encryption method based on two-dimensional discrete memristor hyperchaotic mapping includes (1) chaotic sequence generation; (2) main diagonal row-column diffusion; (3) pixel position scrambling operation; and (4) sub-diagonal row-column diffusion. This method employs a diffusion-scrambling-diffusion structure, using the obtained chaotic matrix to perform main diagonal row-column diffusion, pixel position scrambling, and sub-diagonal row-column diffusion operations on the plaintext image to obtain the ciphertext image. Experimental results show that this invention can resist differential attacks and brute-force attacks, making it difficult for attackers to obtain effective information from the ciphertext image. It not only has fast and efficient encryption capabilities but also high security performance, successfully ensuring the secure transmission of remote sensing images in the network.

The invention discloses an SRP‑PHAT sound source localization method and device based on FPGA, and relates to the technical field of sound source localization. The method comprises: receiving multi-channel audio data of a single-frame microphone array, grouping and cache rearrangement of input data according to channel sequence numbers, performing FFT operations and bandpass filtering operations in sequence from a first group of cached data streams to an nth group of cached data streams, each group of cached data streams being used as a subarray for cross-power spectrum operation, and outputting filtered data forming an n×n array from the subarrays; calculating the cross-power spectrum results of each subarray above the main diagonal of the n×n array, obtaining the cross-power spectrum of the entire array by taking a conjugate, and reordering the cross-power spectrum results; calculating characteristic frequency points of parallel input data, calculating a steering vector, and performing a plane scanning operation with the cross-power spectrum results. The invention realizes the SRP‑PHAT sound source localization method on the FPGA by optimizing the operation structure and method inside the FPGA.

The present invention discloses a method and system for analyzing ground inclination of a pole tower position based on inverse distance weight. The method comprises: obtaining terrain data in a three-dimensional scene, extracting slope data from the terrain data; obtaining an original slope value sampling matrix according to the slope data; calculating the distance from each sampling point of each pole tower in the transmission line to the corresponding pole tower point to obtain an inverse second-order distance matrix; performing inverse distance weight calculation according to the inverse second-order distance matrix to obtain an inverse distance weight matrix; transposing the inverse distance weight matrix to obtain a transposed matrix; multiplying the original slope value sampling matrix with the transposed matrix to obtain a result matrix; obtaining the main diagonal elements of the result matrix to obtain the slope value of each pole tower; and calculating the ground inclination of each pole tower according to the slope value of each pole tower. The present invention fully considers the terrain near the pole tower, and is a comprehensive value of ground inclination with high accuracy.

Taking advantage of the strong robustness of the frequency-domain digital watermarking algorithm, the present invention discloses a blind watermarking method for color digital images based on Contourlet transform and Schur decomposition. According to the characteristics of multi-scale, multi-directionality, translational invariance of the Contourlet transform and low complexity of Schur decomposition, the present invention first performs hierarchical processing on the host image and conducts Contourlet transform on the hierarchical image. Then, the selected coefficient blocks are subjected to Schur decomposition to obtain an upper triangular matrix, and the embedding and blind extraction of the color digital watermark are completed by quantifying the first main diagonal element of the upper triangular matrix. The present invention belongs to the technical field of cyberspace security, not only has good watermark concealment, but also has the characteristics of strong robustness and high security, and is suitable for occasions where color image copyright protection is carried out quickly and efficiently.