Design of 3D Image Visual Communication System for Automatic Reconstruction of Digital Images

Abstract

In order to improve the visual communication effect of the 3D image visual communication system, this article combines the automatic reconstruction technology of digital images to design the 3D image visual communication system. Moreover, this article studies the shock stability of the shock capture scheme by combining entropy generation analysis, linear disturbance analysis, and numerical experiments. In addition to this, this article presents a general method that can be used to suppress numerical shock instabilities in various Godunov-type schemes. It can be seen from the experimental results that the proposed 3D image visual communication system for the automatic reconstruction of digital images has a good visual communication effect.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

There are three methods commonly used in 3D modeling technology. The first is to build a model on 3D production software, use AutoCAD to edit the spatial geometry of the survey area based on topographic map and vector data, and use 3DsMax texture mapping technology to perform texture mapping to generate a textured 3D model. This method is relatively simple, requires a lot of manual operations, has low efficiency, and is suitable for small-scale refined three-dimensional model production. The second method uses laser scanning technology to scan the object to be measured in all directions to obtain a large amount of point cloud data on the surface of the object, removes noise points, and splices the point clouds to obtain the 3D point cloud image of the object to achieve 3D model construction. The third method uses digital photogrammetry technology to obtain the original image of the ground object and uses the image information to restore the original ground object model to realize the reconstruction of the three-dimensional model [1]. The 3D modeling method based on oblique photography technology has become a mature and common method for urban 3D modeling and planning. It has the advantages of high degree of automation, high modeling efficiency, and no need for manual intervention. Moreover, it can quickly build digital products such as urban digital orthophoto (DOM), digital elevation model (DEM), digital orthophoto (DOM), and digital grid map (DRG), which is suitable for large-scale urban 3D model production [2].

When it comes to 3D technology, we have to mention computer graphics, because when we need to use 3D technology to create graphics, we often need to use 3D software to operate. Moreover, the development and progress of 3D software are developed along with the development of computer graphics [3]. Therefore, computer graphics technology is actually the basis for the development of 3D technology. Computer graphics is abbreviated as CG. Its content is very rich, and because it involves a wide range of fields and is closely related to many disciplines, it lacks a very clear definition. The most common explanation that can be found is generally that computer graphics is the science of using mathematical algorithms to convert two-dimensional or three-dimensional graphics into grid form on a computer display [4]. In simple terms, computer graphics mainly studies how to represent graphic elements in computers and how to use computers to perform graphics operations, and the development of algorithms and principles for graphics processing and display. At present, most CG is often considered to be the processing of three-dimensional images, and in fact, it also includes the production of two-dimensional images such as computer painting [5].

Through the development of computer graphics over the past few decades, the content it contains has exceeded the original content, and the current content is very extensive, including computer hardware, image algorithms, graphics standards, raster graphics generation algorithms, surface modeling, polygon modeling, virtual reality, computer animation, graphic interaction technology, augmented reality, mixed reality, natural scene simulation, and computer-aided design. In general, the content of computer graphics can be divided into four major parts: modeling, rendering, animation, and human-computer interaction [6].

When the created 3D models and scenes are available, these 3D models need to be visualized, that is, rendered. Only after rendering can a 3D model become a realistic image, rendering is also a very core existence in traditional computer graphics. In addition, with the different requirements of the final image quality, the rendering of the model is also divided into several levels. The more realistic it is, the higher the computer resources consumed [7]. Today’s computer hardware is becoming more and more powerful, and rendering has undergone great changes with the development. From the initial biased rendering, to the later physically-based rendering, to the recent popular GPU rendering and real-time rendering, computers can render various objects in reality, such as mountains, vegetation, fluids, hair, and skin, as real as the photos taken [8].

With the development of photogrammetry and computer vision, close-range photogrammetry has been used in reverse engineering, cultural relics protection, 3D reconstruction, and other fields. An image is the carrier of feature information in photogrammetry. In order to obtain high-precision measured object information, a high-precision camera is necessary. However, because professional cameras are relatively expensive, and ordinary digital cameras have relatively large distortion, they must be calibrated by the camera. After distortion correction, it can be directly applied to industrial production [9]. Camera calibration is a prerequisite for mapping from 2D to 3D. In recent years, many scholars in photogrammetry and computer vision have conducted in-depth research and proposed many mature methods. Generally, it can be divided into three categories: traditional camera calibration method, active vision camera calibration method, and camera self-calibration method. The traditional camera calibration method uses the known scene structure information, such as control field and calibration block, which is characterized by being suitable for any camera model and having high calibration accuracy, but the process is complicated and limited by the high-precision object information [10]. The active vision camera calibration method is a method of the camera self-calibration algorithm proposed in the field of computer vision using active vision technology. It can usually be solved linearly and has high robustness, but some motion information of the camera is required. Representative methods are based on planes. The orthogonal motion method of homography matrix and the orthogonal motion method based on outer pole are compared [11]; the camera self-calibration method was proposed by Faugeras, Luong, Maybank, etc., which does not require object control points but only needs to capture objects from different directions. It is more flexible to establish the relative relationship between the images, but because it is a nonlinear calibration, it is quite different from the actual results [12]. In order to make up for the shortcomings of the above methods, on this basis, subsequent scholars, combined with automatic coding recognition, use digital coding information to complete the automatic recognition of object square points and image square points and guide the establishment of relative models; representative algorithms are based on ring-shaped numbers. Code recognition, point-based digital code recognition, and digital color code recognition are different [13]. Regardless of the method, a set of information correspondence between the object side and the image side is established, which may be point information, that is, collinear equation, or other information, such as parallel, vanishing point, and other information. After the initial three-dimensional model is established, the calibration parameters are iteratively obtained by the beam method adjustment. The error equation established by the collinear equation is generally used. When there are enough redundant observations, the camera calibration parameters with high precision and reliability can be obtained [14].

The basic information is the corresponding information between the object side and the image side, which can be point or line information. In order to obtain high-precision object-side information and image-side information, the traditional method is to use high-precision control measurement and manual-guided semi-automatic positioning technology [15]; with the development and high integration of computer vision and photogrammetry, and digitally encoded targets, the automatic identification algorithm has gradually replaced the tedious work of traditional manual measurement and puncture, and has overcome the shortcoming that high-precision matching cannot be effectively achieved due to the influence of ground object texture, photographic conditions, and noise. Due to the unique coding and high robustness characteristics of digital coding markers, they are widely used in close-range photogrammetry [16]. The most representative ones are circular coding, dot coding, and color coding. Ring-shaped coding has the characteristics of high precision and simplicity, but has limited coding capacity and poor robustness; point-shaped coding has large coding capacity and high robustness, but because the dots that make up point-shaped coding are too small, it is easily affected by the environment, resulting in a great boundary extraction error; while color coding is more common than ordinary coding, it not only has the characteristics of high accuracy of ring coding, but also has the characteristics of high robustness of point distribution [17].

The acquisition process of oblique images is a multiview photography process, so it is necessary to recover the three-dimensional information of real objects through these two-dimensional images of multiple observation values, while considering the accuracy and effect factors and practicability [18]. The accuracy and effect are mainly affected by various noise and texture environments. For example, under the condition of low texture or repeated texture and high noise, the matching effect is directly affected; practicability is based on the accuracy and effect, considering the complexity of time. It is an objective issue of degree, which is often considered comprehensively with the former [19]. The final result of PMVS to achieve 3D model reconstruction is to use dense patches with normal vectors and center coordinates to approximate the surface of the object; the essence of the algorithm is the regional diffusion algorithm, and the core idea is to perform sparse matching first and then perform sparse matching. The matching results are expanded and error filtered, and iteratively iterates to a certain threshold condition to generate a dense point cloud; unlike other dense matching algorithms, it does not require rough initial values, such as visible convex hulls or bounding boxes, but needs to be based on the grayscale consistency constraints and geometric consistency constraints of the image points [20].

This article combines the automatic reconstruction of digital image technology to design the 3D image visual communication system to improve the 3D film and visual communication effect.

2. Automated Reconstruction Algorithms for Digital Images

2.1. Analysis of Entropy Generation in Godunov-Type Format

The ultimate goal of analyzing the digital image reconstruction algorithm is to enable the Riemann solver to capture shock waves stably without introducing any additional numerical viscosity of linear or nonlinear waves. In order to further study the relationship between entropy generation and shock instability, this article firstly analyzes the numerical entropy generation of a typical HLL-type scheme. Then, by modifying the viscous term of the HLLEM scheme, this article constructs two novel shock capture schemes to control the entropy generation inside the shock structure. In the study, the modified equations in discrete format play an important role, and the method is still applicable when there is a shock solution.

Exact entropy distribution of viscous Navier–Stokes shocks is shown in Figure 1. The entropy generation in numerical format is closely related to the viscosity of its entropy wave. In order to analyze the entropy generation of the shock capture method, this section discusses the entropy wave viscosity of the Godunov-type format, especially the HLL-type format. First, this article considers the approximate solution of the HLLEM format to the Riemann problem, namely,$$\begin{array}{}\text{(1)}& \omega \frac{x}{t};U_L,U_R=\begin{array}{cc}{U}_L,& \frac{x}{t}<S_L,\\ U^{\ast }+x-\widehat{u}t\widehat{{\delta }_2}\widehat{{\alpha }_2}\widehat{R_2},& S_L<\frac{x}{t}<S_R,\\ U_R,& S_R<\frac{x}{t}.\end{array}\end{array}$$

[figure(s) omitted; refer to PDF]

Among them, the viscous term corresponding to shear waves is omitted in the one-dimensional case. The flux function in HLLEM format can be expressed as$$\begin{array}{}\text{(2)}& F_{\text{HLLEM}}=\begin{array}{cc}F{U}_L,& S_L>0,\\ \frac{S_{R}F{U}_L-S_{L}F{U}_R}{S_R-S_L}+\frac{S_R{S}_L}{S_R-S_L}{U}_R-U_L-\widehat{{\delta }_2}\widehat{{\alpha }_2}\widehat{R_2},& S_L\le 0\le <S_R,\\ F{U}_R,& S_R<0.\end{array}\end{array}$$

In order to facilitate further analysis, the following wave velocity estimation is used in formula (2):$$\begin{array}{c}\text{(3)}& S_L=\widehat{u}-\widehat{\alpha },S_R=\widehat{u}-\widehat{\alpha }.\end{array}$$

For the one-dimensional stationary shock problem, the velocity estimation is reasonable. It should be noted that HLLEM format (2) including the definition of wave velocity estimation (3) is equivalent to Roe format (1). Therefore, the conclusions of this analysis also apply to the Roe format. In the subsonic case, the above flux function can be expressed as$$\begin{array}{}\text{(4)}& F_{\text{HLLEM}}=\frac{1}{2}{F}_L-F_R-\frac{1}{2}{Q}_{i+1/2}{U}_R-U_L,\end{array}$$where the numerical viscosity matrix $Q_{i+1/2}$ can be expressed as$$\begin{array}{}\text{(5)}& Q_{i+1/2}=\frac{S_L+S_R}{S_R-S_L}\widehat{A}-\frac{2S_L{S}_R}{S_R-S_L}I+\frac{2S_L{S}_R}{S_R-S_L}\widehat{{\delta }_2}\widehat{R}T{\widehat{R}}^{-1},\end{array}$$and$$\begin{array}{}\text{(6)}& T=\begin{array}{ccc}0& 0& 0\\ 0& 1& 0\\ 0& 0& 0\end{array},\end{array}$$where $\widehat{A}$ represents the Jacobian matrix calculated with the Roe mean variable, and $R={\widehat{R}}_1,{\widehat{R}}_2,{\widehat{R}}_3$ represents the right eigenvector matrix of the Jacobian matrix. The last term on the right-hand side of (5) is negative, which has an effect on reducing the numerical viscosity of the entropy wave. A possible way to increase entropy production is to increase the coefficient, so less viscosity will be imposed on the entropy wave and the corresponding dissipation will be smaller. In order to further clarify this issue, this article considers the modified equation of the conservation scheme of the Euler equation calculated by flux function (4).

The HLLEM format of the one-dimensional Euler equation can be expressed as the following semi-discrete conservation form:$$\begin{array}{}\text{(7)}& \frac{\text{d}}{\text{d}t}{U}_i=-\frac{1}{\Delta x}{F}_{i+1/2}-F_{i-1/2},\end{array}$$where $F_{i+1/2}$ represents the numerical flux, which is used to approximate the flux at interface $x_{i+1/2}$. $\mathrm{\Delta }x$ represents the spatial grid scale and $x_i=\mathrm{\Delta }x,i=1,\mathrm{...},N$. For the Godunov-type forms defined by formulas (4), (5), and (6), (7) can be expressed as the following equivalent viscous form:$$\begin{array}{}\text{(8)}& \frac{\text{d}}{\text{d}t}{U}_i=-\frac{1}{2\Delta x}{F}_{i+1/2}-F_{i-1/2}+\frac{1}{2\Delta x}{Q}_{i+1/2}{U}_{i+1}-U_i-Q_{i-1/2}{U}_i-U_{i-1},\end{array}$$where $Q_{i+1/2}$ is the numerical viscosity coefficient matrix. The modified equation of formula (8) can be expressed as$$\begin{array}{}\text{(9)}& {\tilde{U}}_t+F{\tilde{U}}_x=\frac{1}{2}Q\tilde{U}{\tilde{U}}_{xx}\Delta x,\end{array}$$where $\tilde{U}$ represents the polynomial obtained by interpolating the conserved variable U around the grid point x. Under the modified equation, the second- and higher-order error terms in the equation are omitted.

In order to obtain the modified equation in the HLLEM format, the numerical viscosity coefficient matrix in formula (9) needs to be determined. By making U and U approach U, it can be expressed as$$\begin{array}{}\text{(10)}& Q_{\text{HLLEM}}U=\frac{S_L+S_R}{S_R-S_L}AU-\frac{2S_L{S}_R}{S_R-S_L}I+\frac{2S_L{S}_R}{S_R-S_L}{\delta }_{2}RUT{R}^{-1}U,\end{array}$$where A (U) is the Jacobian matrix of the flux function F(U), $AU=\partial FU/\partial U$, I represents the identity matrix, sL = SL (U, U), and SR = SR(U, U).

In order to analyze the numerical viscosity of specific eigenwaves, especially the artificial viscosity of entropy waves, it is necessary to obtain a modified equation of the characteristic form. The right eigenvectors of the Jacobian matrix A(U) are R₁, R₂, and R₃, and the corresponding eigenvalues can be expressed as$$\begin{array}{c}\text{(11)}& {\lambda }_1=u-a,{\lambda }_2=u,\\ {\lambda }_3=u+a.\end{array}$$

The transformation matrix R and the corresponding inverse matrix R⁻¹ are applied to modified formula (9), and in the case where $\tilde{U}$ approaches U, it can be known that$$\begin{array}{}\text{(12)}& R^{-1}{U}_t+R^{-1}AR{R}^{-1}{U}_x=\frac{1}{2}{R}^{-1}QUR{R}^{-1}{U}_{xx}\Delta x.\end{array}$$

The above formula can further be simplified by the relation $A=R\wedge R^{-1}$ and $dW=R^{-1}dU$, and the modified equation in characteristic form can be expressed as$$\begin{array}{}\text{(13)}& W_t+\wedge W_x=\frac{1}{2}v{W}_{xx}\Delta x,\end{array}$$where $\wedge$  = diag (u − a,, u + a) is the eigenvalue of the matrix, and the eigenvariable can be defined as$$\begin{array}{}\text{(14)}& dW=\begin{array}{c}d{w}_1,\\ d{w}_2,\\ d{w}_3,\end{array}=\begin{array}{c}\frac{1}{2a^2}dp-\frac{\rho }{2a}du\\ d\rho -\frac{dp}{a^2}\\ \frac{1}{2a^2}dp+\frac{\rho }{2a}du\end{array}\end{array}.$$

V is the viscosity coefficient matrix of characteristic form, which can be expressed as the following form:$$\begin{array}{}\text{(15)}& V=\frac{S_L+S_R}{S_R-S_L}\wedge U-\frac{2S_L{S}_R}{S_R-S_L}I+\frac{2S_L{S}_R}{S_R-S_L}{\delta }_{2}T.\end{array}$$

Formula (13) can be re-expressed in the form of explicit representation of eigenwaves, that is,$$\begin{array}{c}\text{(16)}& \frac{\partial w_i}{\partial t}={\lambda }_i\frac{\partial w_i}{\partial x},=\frac{1}{2}{u}_i\frac{{\partial }^2}{\partial x^2}{w}_i\Delta x,\end{array}$$where u can be interpreted as a numerical viscosity coefficient corresponding to a specific eigenwave. To evaluate the numerical entropy generation, this article only considers the artificial viscosity coefficient corresponding to the entropy wave dw₂, that is,$$\begin{array}{}\text{(17)}& {\mu }_2=\frac{S_L+S_R}{S_R-S_L}{\lambda }_2-\frac{2S_L{S}_R}{S_R-S_L}+\frac{2S_L{S}_R}{S_R-S_L}{\delta }_2.\end{array}$$

It can be observed that the last term on the right-hand side of (17) is nonpositive, which acts to reduce the numerical viscosity of the entropy wave. For the HLL format, the coefficient ${\delta }_2$ is equal to zero. That’s why it contains too much entropy wave viscosity. One way to increase entropy production is to increase the coefficient ${\delta }_2$, so that less dissipation will be imposed on the entropy wave. However, this scheme is infeasible because the coefficients in formula (5) should satisfy the TVD condition, so that the viscous matrix is valid, that is,$$\begin{array}{c}\text{(18)}& {\lambda }_k^{\ast }\ge {\widehat{\lambda }}_k,k=2.\end{array}$$where ${\lambda }_k^{\ast }$ represents the eigenvalue of $Q_{i+1/2}$. From formula (18), it can be known that the inverse dissipation coefficient $\widehat{{\delta }_2}$ must satisfy the following conditions:$$\begin{array}{}\text{(19)}& \widehat{{\delta }_2}\le \frac{\widehat{a}}{\widehat{a}+\widehat{u}}.\end{array}$$

Considering the definition of the coefficient $\widehat{{\delta }_2}$ in formula (19), it is unreasonable to increase the coefficient $\widehat{{\delta }_2}$. Therefore, this article adopts the method of modifying the viscosity term in the wave intensity $\widehat{{\alpha }_2}$ to control the dissipation of the entropy wave, so as to realize the control of the entropy generation.

The average state quantity Un12 in formula (1) can be corrected as$$\begin{array}{}\text{(20)}& {\omega }_{\text{HLLEM}-\rho }\frac{x}{t};U_L,U_R=\begin{array}{cc}{U}_L,& \frac{x}{t}<S_L,\\ U^{\ast }+\widehat{{\delta }_2}x-\widehat{u}t{\widehat{{\alpha }_2}}_2^{\ast }\widehat{R_2},& S_L<\frac{x}{t}<S_R,\\ U_R,& S_R<\frac{x}{t}.\end{array}\end{array}$$

The correction wave intensity used in it can be expressed as$$\begin{array}{}\text{(21)}& \widehat{{\alpha }_2^{\ast }}=f_{\rho }\cdot \Delta p-\frac{\Delta p}{\widehat{{\alpha }_2}}.\end{array}$$

To identify such corrections, the correction format is denoted by $\text{HLLEM}-\rho$. The modification of the intermediate state quantity $\text{HLLEM}-\rho$ does not change the integral relation (24), so Riemann solver (20) can be expressed in the following conservative form:$$\begin{array}{}\text{(22)}& F_{i+1.2}^{\text{HLLEM}-\rho }=\frac{S_{R}F{U}_L-S_{L}F{U}_R}{S_R-S_L}+\frac{S_L{S}_R}{S_R-S_L}{U}_R-U_L-{\widehat{\delta }}_2\widehat{{\alpha }_2^{\ast }}\widehat{R_2}.\end{array}$$

The corrected wave intensity $\widehat{{\alpha }_2^{\ast }}$ in the formula is shown in formula (21).

In addition to this, the differential pressure term in $\widehat{{\alpha }_2^{\ast }}$ can also be modified. Similarly, the new format can be identified as $\text{HLLEM}-\rho$, which is defined as follows:$$\begin{array}{}\text{(23)}& F_{i+1.2}^{\text{HLLEM}-\rho }=\frac{S_{R}F{U}_L-S_{L}F{U}_R}{S_R-S_L}+\frac{S_L{S}_R}{S_R-S_L}{U}_R-U_L-\widehat{{\delta }_2}\widehat{{\alpha }_2^{\ast }}\widehat{R_2}.\end{array}$$

The corrected wave intensity in the formula is defined as$$\begin{array}{}\text{(24)}& \widehat{{\alpha }_2^{\ast }}=\Delta p-f_P\frac{\Delta p}{\widehat{{\alpha }_2}}.\end{array}$$

The coefficients $f_{\rho }$ and $f_P$ are used to control the numerical viscosity on the entropy wave. In order to elucidate their mechanism of action, the dissipative mechanisms of these Godunov-type schemes are analyzed [21].

If the numerical viscosity coefficient matrix in HLLEM-p format is determined, its modified equation can be obtained. The flux functions defined in formulas (23) and (24) can be expressed as$$\begin{array}{}\text{(25)}& F_{\text{HLLEM}-\rho }=\frac{1}{2}{F}_L+F_R-\frac{1}{2}{Q}_{i+1/2}{U}_R-U_L.\end{array}$$

The numerical viscosity coefficient matrix $Q_{i+1/2}$ in formula is defined as follows:$$\begin{array}{}\text{(26)}& Q_{i+1/2}=\frac{S_L+S_R}{S_R-S_L}\widehat{A}-\frac{2S_L{S}_R}{S_R-S_L}I+\frac{2S_L{S}_R}{S_R-S_L}{\delta }_2\widehat{R}T\widehat{R^{-1}}-\frac{2S_L{S}_R}{S_R-S_L}{\delta }_2\widehat{R}T\widehat{R^{\ast ,-1}},\end{array}$$and contains the following relations:$$\begin{array}{}\text{(27)}& \widehat{{\alpha }_2^{\ast }}\widehat{R_2}=\widehat{R}T\widehat{R^{\ast ,-1}}-U_R-U_L.\end{array}$$

The modified equation using the conservation scheme (7) of the $\text{HLLEM}-\rho$ flux can be expressed as$$\begin{array}{}\text{(28)}& U_t+F{U}_x=\frac{1}{2}{Q}_{\text{HLLEM}-\rho }U{U}_{xx}\Delta x.\end{array}$$

The viscosity coefficient matrix can be expressed as$$\begin{array}{c}\text{(29)}& Q_{\text{HLLEM}}U=\frac{S_L+S_R}{S_R-S_L}AU-\frac{2S_L{S}_R}{S_R-S_L}I+\frac{2S_L{S}_R}{S_R-S_L}{\delta }_{2}RUT{R}^{-1}U-\frac{2S_L{S}_R}{S_R-S_L}{\delta }_{2}RUT{R}^{\ast ,-1}U.\end{array}$$

Following the steps in (12), the modified equation of the characteristic form can be expressed as$$\begin{array}{}\text{(30)}& W_t+\wedge W_x=\frac{1}{2}V{W}_{xx}\Delta x+f_p-1\frac{2S_L{S}_R}{S_R-S_L}{\delta }_{2}T{W}_{xx}^2\Delta x.\end{array}$$

The following relationship is used in the formula:$$\begin{array}{}\text{(31)}& dW=d{W}_1-d{W}^2,\end{array}$$and$$\begin{array}{c}\text{(32)}& d{W}^1=\begin{array}{c}d{w}_1\\ d{w}_2\\ d{w}_3\end{array}=\begin{array}{c}\frac{1}{2a^2}dp-\frac{\rho }{2a}du\\ d\rho \\ \frac{1}{2a^2}dp+\frac{\rho }{2a}du\end{array}d{W}^2\\ =\begin{array}{c}0\\ d{w}_2^2\\ 0\end{array}=\begin{array}{c}0\\ \frac{dp}{a^2}\\ 0\end{array}.\end{array}$$

In order to evaluate the numerical viscosity of the entropy wave, the following modified equation corresponding to the entropy wave $d{w}_2$ needs to be considered, which can be obtained from formula (30):$$\begin{array}{}\text{(33)}& \frac{\partial w_2}{\partial t}+{\lambda }_2\frac{\partial w_2}{\partial t}=\frac{1}{2}{\mu }_2\frac{{\partial }^2}{\partial x^2}{w}_2\Delta x+f_p-1\frac{S_L{S}_R}{S_R-S_L}{\delta }_2\frac{{\partial }^2}{\partial x^2}{w}_2^2\Delta x.\end{array}$$

The definition of the numerical viscosity coefficient pz is shown in formula (17). By combining formulas (31)–(33), it can be split into two subformulas as follows:$$\begin{array}{}\text{(34)}& \frac{\partial w_1^1}{\partial t}+{\lambda }_2\frac{\partial w_2^1}{\partial x}=\frac{1}{2}{\mu }_2\frac{{\partial }^2}{\partial x^2}{w}_2^1\Delta x,\text{(35)}& \frac{\partial w_2^1}{\partial t}+{\lambda }_2\frac{\partial w_2^2}{\partial x}=\frac{1}{2}{\mu }_2\frac{{\partial }^2}{\partial x^2}{w}_2^2\Delta x+f_p-1\frac{S_L{S}_R}{S_R-S_L}{\delta }_2\frac{{\partial }^2}{\partial x^2}{w}_2^2\Delta x.\end{array}$$

The last term on the right side of (35) is the artificial viscosity term, which acts as a dissipative wave $d{w}_2^2$. Observing the last term on the right side of formula (35), we can see that if $f_p$ decreases, then $d{w}_2^2$ also decreases. Combining formulas (31) and (32), it can be seen that the entropy generation $d{w}_2^2$ will increase. Similarly, if $f_p$ increases, then the entropy production $d{w}_2^2$ also increases. However, it should be noted that when $f_p$ is greater than unity, the stability condition (18) for contact discontinuity no longer holds. Therefore, the HLLEM-p format cannot be used for multidimensional flow field simulations because of the possibility of contact discontinuities in the flow field. However, for one-dimensional stationary shock problems, the scheme still works.

2.2. Quantitative Analysis of the Relationship between Entropy Generation and Shock Stability

The research shows that the dissipation corresponding to the entropy wave can be effectively weakened by increasing the density difference term Δ$p$ or decreasing the pressure difference term Δ$p$ in the wave intensity term $\widehat{{\alpha }_2}$. If these correction schemes are applied at the shock, the entropy generation inside the numerical shock structure will become larger. The coefficients $f_{\rho }$ and $f_p$ play an important role in controlling the entropy generation. Since it is difficult to determine the values of the coefficients $f_{\rho }$ and $f_p$ analytically, this article resorts to numerical experiments to determine these coefficients.

Considering the one-dimensional steady normal shock problem, the numerical calculation adopts the first-order HLLEM format. The three-step Runge–Kutta method with first-order accuracy is applied, the CFL number was set to 0.5, and the calculation process is carried out until 40,000 time steps. At the same time, the Mach number is set to M₀ = 6.0, and the shock position is set to $\epsilon$  = 0.3. Numerical fluxes at the interface of grid M are calculated using the modified formats HLLEM-p and HLLEM-p to control entropy generation within grid M₀. The HLLEM format will have a one-dimensional numerical phenomenon. Therefore, in order to obtain a stable and convergent solution, this article gradually increases the magnitude of the coefficient $f_{\rho }$ or decreases the magnitude of the coefficient $f_{\rho }$ to increase the entropy generation in the intermediate grid M. In Figure 2, under different values of $f_{\rho }$ and $f_p$ coefficients, the entropy distribution curve of the stable solution is obtained. The corresponding residual value convergence curve is shown in Figure 3. It can be seen from the figure that as the coefficient $f_{\rho }$ becomes larger or the coefficient $f_{\rho }$ becomes smaller, the entropy $S=\log p/{\rho }^{\gamma }-\log p_L/{\rho }_L$ in the numerical shock structure (the grid M) will increase. This result is consistent with the entropy generation analysis presented in subsection 1. These results suggest that if sufficient entropy can be generated inside the shock structure, the instability can be successfully removed. The convergence curves presented in Figure 3 show that more entropy production will allow the computation to converge to a stable solution at a faster rate. It should be noted that similar results can be obtained for other Mach numbers M⁰ and shock position parameter $\epsilon$. In Figure 4, this article gives the stability thresholds of $f_{\rho ,\text{threshold}}$ and $f_{p,\text{threshold}}$ under different Mach numbers M₀, and the corresponding entropy thresholds are shown in Figure 5.

[figure(s) omitted; refer to PDF]

As shown in Figure 6, the conserved quantities inside the numerical shock structure suffer perturbation errors when crossing the shock. Due to the difference of physical quantities in the lateral direction of the shock wave, a new Riemann problem will arise in the lateral direction, as shown by the red solid line in Figure 6.

[figure(s) omitted; refer to PDF]

3. 3D Image Visual Communication System for Automatic Reconstruction of Digital Images

The structure of the depth estimation network constructed in this article is shown in Figure 7. For a reference image whose depth is to be estimated and several adjacent images, image features are extracted based on the pyramid structure, and then, the initial depth map is estimated. After that, based on the multiscale optimization strategy, residual depth estimation is performed on the initial results, and the optimized depth map is finally output.

[figure(s) omitted; refer to PDF]

The functional block diagram of the system is shown in Figure 8.

(1) User interaction. It refers to the interaction between the user and the volume data, such as dragging, rotating, zooming in, and zooming out.

(2) Data loading. It includes the loading of volume data and the loading of transfer functions. Volume data loading loads volume data according to the size, format, and bit depth of the volume data, and encapsulates the volume data into a 3D texture for use by the fragment shader program. Transfer function loading generates a transfer function lookup table based on the transfer function curve for use by fragment shader programs.

(3) The shader program is loaded. Load different shader programs according to the different effects to be achieved (pre-interpolation shading classification, postinterpolation shading classification, pre-integration shading classification, etc.). Shader programs include vertex shader programs and fragment shader programs.

(4) Matrix transformation. Before rendering, the volume data model is subjected to necessary model-view transformation and projection transformation to prepare for rendering processing. Finally, the transformation matrix is sent to the GPU for use by the vertex shader program.

(5) Image output. According to the user's choice, the output mode of the controlled drawing result is two-dimensional image or three-dimensional stereoscopic image, and the three-dimensional stereoscopic image is output for autostereoscopic display.

[figure(s) omitted; refer to PDF]

On the basis of the above research, this article verifies the effect of the model proposed in this article, and quantitatively processes the algorithm proposed in this article and some other 3D reconstruction methods based on multiview images, statistical accuracy error, integrity error, and integrity error. The statistical error results are shown in Table 1.

Table 1

Quantitative test results of 3D reconstruction quality.

Num	Accuracy error	Integrity error	Global error
1	0.81092	0.53835	0.67512
2	0.59364	0.91083	0.75272
3	0.33271	1.1543	0.74399
4	0.26578	1.15721	0.71198
5	0.4365	1.01171	0.72362
6	0.44232	0.62662	0.53447
7	0.37151	0.43844	0.40449
8	0.39382	0.42098	0.4074
9	0.35017	0.40837	0.37927
10	0.34241	0.39576	0.3686

From the above test results, it can be seen that the 3D image visual communication system for automatic reconstruction of digital images proposed in this article has a good 3D image reconstruction effect. On this basis, the visual communication effect of the 3D image visual communication system for the automatic reconstruction of digital images is evaluated, as shown in Figure 9.

[figure(s) omitted; refer to PDF]

From the above test results, it can be seen that the 3D image visual communication system for the automatic reconstruction of digital images proposed in this article has a good visual communication effect.

4. Conclusion

As a classic problem in computer vision, 3D reconstruction has experienced a long process of development, produced rich results, and played an important role in people's life and production. In some traditional fields, such as in industrial manufacturing, reverse engineering based on 3D reconstruction technology can effectively improve the visual communication effect of 3D film and television. The 3D reconstruction technology using multiview image depth map has become one of the most mainstream research directions at present because of its high flexibility and good reconstruction accuracy. This article combines the automatic reconstruction of digital image technology to design the 3D image visual communication system to improve the 3D film and visual communication effect. It can be seen from the test results that the 3D image visual communication system for the automatic reconstruction of digital images proposed in this article has a good visual communication effect.

Acknowledgments

This study was sponsored by Xinyang Vocational and Technical College.

References

[1] A. Abbasi, H. Chen, A. Salem, "Sentiment analysis in multiple languages," ACM Transactions on Information Systems, vol. 26 no. 3,DOI: 10.1145/1361684.1361685, 2008.

[2] S. Agrawal, R. K. Singh, Q. Murtaza, "Prioritizing critical success factors for reverse logistics implementation using fuzzy-TOPSIS methodology," Journal of Industrial Engineering International, vol. 12 no. 1, pp. 15-27, DOI: 10.1007/s40092-015-0124-8, 2016.

[3] P. Bosc, D. Dubois, H. Prade, "Fuzzy functional dependencies and redundancy elimination," Journal of the American Society for Information Science, vol. 49 no. 3, pp. 217-235, DOI: 10.1002/(sici)1097-4571(199803)49:3<217::aid-asi4>3.0.co;2-d, 1998.

[4] R. Carrasco, M. G. Ordaz, F. J. M. López, "Science fiction and bodies of the future: alternative gender realities in Hollywood Cinema," Journal of Futures Studies, vol. 20 no. 2, pp. 67-80, 2015.

[5] E. Castano, "Art films foster theory of mind," Humanities and Social Sciences Communications, vol. 8 no. 1,DOI: 10.1057/s41599-021-00793-y, 2021.

[6] S. C. Chang, "Market size matters? an approach to illustrate the market preference of Hong Kong-mainland China co-production cinema," Journal of International Communication, vol. 26 no. 1, pp. 125-149, DOI: 10.1080/13216597.2020.1728358, 2020.

[7] R. D. Cohen, "The delinquents: censorship and youth culture in recent U. S. History," History of Education Quarterly, vol. 37 no. 3,DOI: 10.2307/369445, 1997.

[8] A. Dhall, R. Goecke, S. Lucey, T. Gedeon, "Collecting large, richly annotated facial-expression databases from movies," IEEE multimedia, vol. 19 no. 3, pp. 34-41, DOI: 10.1109/mmul.2012.26, 2012.

[9] M. Gulistan, I. Beg, N. Yaqoob, "A new approach in decision making problems under the environment of neutrosophic cubic soft matrices," Journal of Intelligent and Fuzzy Systems, vol. 36 no. 1, pp. 295-307, DOI: 10.3233/jifs-181296, 2019.

[10] E. Taiebi JavidTaiebi Javid, M. Nazari, M. R. GhaeliGhaeli, "Social media and e-commerce: a scientometrics analysis," International Journal of Data and Network Science, vol. 3 no. 3, pp. 269-290, DOI: 10.5267/j.ijdns.2019.2.001, 2019.

[11] L. Jayyusi, "Hollywood's transnational imaginaries: colonial agency and vision fromIndiana JonestoWorld War Z," Continuum, vol. 32 no. 3, pp. 355-369, DOI: 10.1080/10304312.2018.1449806, 2018.

[12] T. Jojola, "Absurd reality: hollywood goes to the Indians," Film and History: An Interdisciplinary Journal of Film and Television Studies, vol. 23 no. 1, 1993.

[13] A. Kaplan, M. Haenlein, "Rulers of the world, unite! the challenges and opportunities of artificial intelligence," Business Horizons, vol. 63 no. 1, pp. 37-50, DOI: 10.1016/j.bushor.2019.09.003, 2020.

[14] B. Kuklick, "Fascism comes to America," International Journal for History, Culture and Modernity, vol. 6 no. 1,DOI: 10.18352/hcm.547, 2018.

[15] S. Kumar, K. De, P. P. Roy, "Movie recommendation system using sentiment analysis from microblogging data," IEEE Transactions on Computational Social Systems, vol. 7 no. 4, pp. 915-923, DOI: 10.1109/tcss.2020.2993585, 2020.

[16] R. A. M. NoorNoor, Z. Ahmad, M. M. Don, M. H. Uzir, "Modelling and control of different types of polymerization processes using neural networks technique: a review," Canadian Journal of Chemical Engineering, vol. 88 no. 6, pp. 1065-1084, DOI: 10.1002/cjce.20364, 2010.

[17] C. Porcel, A. Tejeda-Lorente, M. A. Martínez, E. Herrera-Viedma, "A hybrid recommender system for the selective dissemination of research resources in a technology transfer office," Information Sciences, vol. 184 no. 1,DOI: 10.1016/j.ins.2011.08.026, 2012.

[18] S. S. Sundar, "Rise of machine agency: a framework for studying the psychology of human-AI interaction (haii)," Journal of Computer-Mediated Communication, vol. 25 no. 1, pp. 74-88, DOI: 10.1093/jcmc/zmz026, 2020.

[19] W. C. V. Wu, R. J. Wang, N. S. Chen, "Instructional design using an in-house built teaching assistant robot to enhance elementary school English-as-a-foreign-language learning," Interactive Learning Environments, vol. 23 no. 6, pp. 696-714, DOI: 10.1080/10494820.2013.792844, 2015.

[20] B. Zou, M. NurudeenNurudeen, C. Zhu, Z. Zhang, R. Zhao, L. Wang, "A nfuzzy crime prediction model based on video analysis," Chinese Journal of Electronics, vol. 27 no. 5, pp. 968-975, DOI: 10.1049/cje.2018.02.019, 2018.

[21] W. Xie, Z. Tian, Ye Zhang, H. Yu, W. Ren, "Further studies on numerical instabilities of Godunov-type schemes for strong shocks," Computers & Mathematics with Applications, vol. 102, pp. 65-86, DOI: 10.1016/j.camwa.2021.10.008, 2021.