Removing Random-Valued Impulse Noises by a Two-Staged Nonlinear Filtering Method

Digital images are frequently contaminated by impulse noise(IN)during acquisition and transmission.The removal of this noise from images is essential for their further processing.In this paper,a two-staged nonlinear filtering algorithm is proposed for removing random-valued impulse noise(RVIN)from digital images.Noisy pixels are identified and corrected in two cascaded stages.The statistics of two subsets of nearest neighbors are employed as the criterion for detecting noisy pixels in the first stage,while directional differences are adopted as the detector criterion in the second stage.The respective adaptive median values are taken as the replacement values for noisy pixels in each stage.The performance of the proposed method was compared with that of several existing methods.The experimental results show that the performance of the suggested algorithm is superior to those of the compared methods in terms of noise removal,edge preservation,and processing time.

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.The staggered DG scheme defines the discrete pressure on the primal triangular mesh,while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh.In this paper,a new pair of equal-order-interpolation velocity-pressure finite elements is proposed.On the primary triangular mesh(the pressure elements),the basis functions are piecewise polynomials of degree N and are allowed to jump on the boundaries of each triangle.On the dual mesh instead(the velocity elements),the basis functions consist in the union of piecewise polynomials of degree N on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries,while they are continuous inside.In other words,the basis functions on the dual mesh arc built by continuous finite elements on the subtriangles.This choice allows the construction of an efficient,quadrature-free and memory saving algorithm.In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations,the arbitrary high order of accuracy in time is achieved through tire use of time-dependent test and basis functions,in combination with simple and efficient Picard iterations.Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes,not only from a computer memory point of view,but also concerning the computational time.

Land Surface Temperature Change in Hangzhou Urban Area in Recent 20 Years Based on Landsat Image

Since the reform and opening-up in 1978, the urbanization level of our country has been continuously improved and the urban development has made great progress. However, with the rapid expansion of urban construction land, the population density and building density have been greatly increased, resulting in the urban heat island effect, which has negative impact on the urban thermal environment and restricts the high-quality development of urbanization. This paper focuses on how the urban surface thermal environment of Hangzhou changes in 20 years. In this paper, the characteristics of land surface temperature (LST) in Hangzhou urban area from 2000 to 2020 were studied by using Landsat images. The radiative transfer equation method is used to retrieve the land surface temperature, and the retrieval results are analyzed. The results show that: 1) the land surface temperature in Hangzhou city area has a slight upward trend in the past 20 years;2) the area of high temperature area is expanding;3) the land surface temperature in the city center area has decreased significantly in the past 20 years, while the ground temperature in other areas around the city center has increased significantly.

SEMI-IMPLICIT SPECTRAL DEFERRED CORRECTION METHODS BASED ON SECOND-ORDER TIME INTEGRATION SCHEMES FOR NONLINEAR PDES

In[20],a semi-implicit spectral deferred correction(SDC)method was proposed,which is efficient for highly nonlinear partial differential equations(PDEs).The semi-implicit SDC method in[20]is based on first-order time integration methods,which are corrected iteratively,with the order of accuracy increased by one for each additional iteration.In this paper,we will develop a class of semi-implicit SDC methods,which are based on second-order time integration methods and the order of accuracy are increased by two for each additional iteration.For spatial discretization,we employ the local discontinuous Galerkin(LDG)method to arrive at fully-discrete schemes,which are high-order accurate in both space and time.Numerical experiments are presented to demonstrate the accuracy,efficiency and robustness of the proposed semi-implicit SDC methods for solving complex nonlinear PDEs.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

Superconvergence of Direct Discontinuous Galerkin Methods:Eigen-structure Analysis Based on Fourier Approach

This paper investigates superconvergence properties of the direct discontinuous Galerkin(DDG)method with interface corrections and the symmetric DDG method for diffusion equations.We apply the Fourier analysis technique to symbolically compute eigenvalues and eigenvectors of the amplification matrices for both DDG methods with different coefficient settings in the numerical fluxes.Based on the eigen-structure analysis,we carry out error estimates of the DDG solutions,which can be decomposed into three parts:(i)dissipation errors of the physically relevant eigenvalue,which grow linearly with the time and are of order 2k for P^(k)(k=2,3)approximations;(ii)projection error from a special projection of the exact solution,which is decreasing over the time and is related to the eigenvector corresponding to the physically relevant eigenvalue;(iii)dissipative errors of non-physically relevant eigenvalues,which decay exponentially with respect to the spatial mesh sizeΔx.We observe that the errors are sensitive to the choice of the numerical flux coefficient for even degree P^(2)approximations,but are not for odd degree P^(3)approximations.Numerical experiments are provided to verify the theoretical results.

Influence of Humidity and Air Pressure on the Ion Mobility Based on Drift Tube Method

The measurement of ion mobility is important for calculating the corona loss of the transmission lines,and for establishing corona loss and ion current models that take humidity and air pressure into consideration.This paper describes a needle-ring corona discharge experiment setup based on a prior study,which simulates different humidity and air pressure conditions;it is observed that the execution of the ion gate is improved in this experiment from off-on-off to off-on to achieve higher ion current amplitude.Additionally,positive and negative ion mobility under different humidity and air pressure is obtained,with the positive ion mobility measurements in the range of 1.1126 to 1.9167 cm2/V·s,and the negative ion mobility measurements in the range of 1.3574 to 2.5643 cm2/V·s.The experimental results indicate that ion mobility decreases nonlinearly with increase in humidity,but trends towards saturation in the 30%–70%relative humidity range.Finally,a parameter correction method for calibrating the relationship among the ion mobility,humidity,and air pressure is proposed.

Improved inclination correction method applied to the guidance system of agricultural vehicles

This research introduces a new inclination correction method with increased accuracy applied to the guidance system of an agricultural vehicle.The method considers the geometry of a robot tractor and uses an Inertial Measurement Unit(IMU)to correct the lateral error of the RTK-GPS antenna measurements raised by the tractor's inclinations.A parameters optimization experiment and an automatic guidance experiment under real working conditions were used to compare the accuracy of a traditional correction method with the new correction method,by calculating the RMSE of the lateral error and the error reduction percentage.An additional tuned correction method was found by using a simple analytical method to find the optimal variables that reduced the lateral error to a minimum.The results indicate that the new correction method and the tuned correction method display a significant error reduction percentage compared to the traditional correction method.The methods could correct the GPS lateral error caused by the roll inclinations in real-time.The resulting lateral deviation caused by the tractor's inclinations could be reduced up to 23%for typical travelling speeds.

An Improved Method for Defining Short-Term Climate Anomalies

A conventional method to define short-term climate anomalies for atmospheric and oceanic variables,recommended by the World Meteorological Organization(WMO),is the departure from a 30-yr climatological mean in the preceding three decades.Such a method,however,introduces spurious errors such as sudden jumps and artificial trends.A new method,named a trend correctional method,is introduced to eliminate the errors.To demonstrate the capability of this new method,we examine a set of idealized cases first by superposing a"true"interannual or interdecadal signal onto a linear or a nonlinear trend.Comparing to the conventional method,the trend correctional method is able to retain,to a large extent,the"true"anomaly signals.Next,we examined real-time indices.The anomaly time series derived based on the trend correctional method show a better agreement with the observed anomaly series.The rootmean-square error is greatly improved,comparing to that calculated based on the conventional method.Therefore,the results from both the idealized and real cases demonstrate that the new method has a clear advantage to the conventional method in deriving true climate anomalies,in particular under the ongoing global warming circumstance.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research,and lithology classification,mineral analysis and porosity inversion are important links in reservoir evaluation.However,affected by the diverse lithology,complicated mineral and widespread alteration,conventional logging lithology classification and mineral inversion become considerably difficult.At the same time,owing to the limitation of the wireline log response equation,the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations.To overcome those issues,this study takes the South China Sea as an example,and combines multi-scale data such as micro rock slices,petrophysical experiments,wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks.Specifically,we define the lithology and mineral characteristics through core slices and mineral data,and establish an igneous multi-mineral volumetric model.Then we determine element cutting log correction method based on core element data,and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity.However,it is always difficult to determine the elemental eigenvalues of different minerals in inversion.This paper uses multiple linear regression methods to solve this problem.Finally,an integrated inversion technique for altered igneous formations was developed.The results show that the corrected element cutting log are in good agreement with the core element data,and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction.The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

Semi-Implicit Spectral Deferred Correction Method Based on the Invariant Energy Quadratization Approach for Phase Field Problems

This paper presents a high order time discretization method by combining the semi-implicit spectral deferred correction method with energy stable linear schemes to simulate a series of phase field problems.We start with the linear scheme,which is based on the invariant energy quadratization approach and is proved to be linear unconditionally energy stable.The scheme also takes advantage of avoiding nonlinear iteration and the restriction of time step to guarantee the nonlinear system uniquely solvable.Moreover,the scheme leads to linear algebraic system to solve at each iteration,and we employ the multigrid solver to solve it efficiently.Numerical re-sults are given to illustrate that the combination of local discontinuous Galerkin(LDG)spatial discretization and the high order temporal scheme is a practical,accurate and efficient simulation tool when solving phase field problems.Namely,we can obtain high order accuracy in both time and space by solving some simple linear algebraic equations.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Fixing an illusion–an empirical assessment of correction methods for the terrain reversal effect in satellite images

Identifying land forms and land cover classes are important tasks in image interpretation.Sometimes,a phenomenon called terrain reversal effect(TRE)causes an inverted perception of 3D forms.When this inversion occurs,valleys appear as ridges and vice versa.While the TRE can severely impair the ability to identify 3D land forms,‘correcting’for the TRE in imagery can introduce new problems.Importantly,one of most commonly-proposed methods–shaded relief map(SRM)overlay–appears to impair the ability to identify land cover classes.In this paper,we report a comparative empirical evaluation of an SRM overlay solution,and its‘enhanced’versions supported by various other cues(stereopsis,motion,labels).In response to the different solutions,we measure the effectiveness,efficiency,confidence and preferences of our participants in land form and land cover identification tasks.All examined methods significantly improve the ability to detect land forms accurately,but they also impair the ability to identify the land cover classes to different degrees.Additionally,participants’visualization preferences contradict their performance with them,calling for reflection on the visual effects of the applied correction methods.Based on the study,recommendations concerning the correction of the TRE are drawn,and gaps are identified.

Adaptive local discontinuous Galerkin methods with semi-implicit time discretizations for the Navier-Stokes equations

In this paper,we present a mesh adaptation algorithm for the unsteady compressible Navier-Stokes equations under the framework of local discontinuous Galerkin methods coupled with implicit-explicit Runge-Kutta or spectral deferred correction time discretization methods.In both of the two high order semi-implicit time integration methods,the convective flux is treated explicitly and the viscous and heat fluxes are treated implicitly.The remarkable benefits of such semi-implicit temporal discretizations are that they can not only overcome the stringent time step restriction compared with time explicit methods,but also avoid the construction of the large Jacobian matrix as is done for fully implicit methods,thus are relatively easy to implement.To save computing time as well as capture the flow structures of interest accurately,a local mesh refinement(h-adaptive)technique,in which we present detailed criteria for selecting candidate elements and complete strategies to refine and coarsen them,is also applied for the Navier-Stokes equations.Numerical experiments are provided to illustrate the high order accuracy,efficiency and capabilities of the semi-implicit schemes in combination with adaptive local discontinuous Galerkin methods for the Navier-Stokes equations.

A local discontinuous Galerkin method for the pattern formation dynamical model in polymerizing action flocks

In this paper,we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks.Optimal error estimates for the density and filament polarization in different norms are established.We use a semi-implicit spectral deferred correction time method for time discretization,which allows a relative large time step and avoids computation of a Jacobian matrix.Numerical experiments are presented to verify the theoretical analysis and to show the capability for simulations of action wave formation.

Optimal reduction of anthropogenic emissions for air pollution control and the retrieval of emission source from observed pollutantsⅢ:Emission source inversion using a double correction iterative method

Using the incomplete adjoint operator method in part I of this series of papers,the total emission source S can be retrieved from the pollutant concentrationsρob obtained from the air pollution monitoring network.This paper studies the problem of retrieving anthropogenic emission sources from S.Assuming that the natural source Sn is known,and as the internal source Sc due to chemical reactions is a function of pollutant concentrations,if the chemical reaction equations are complete and the parameters are accurate,Sc can be calculated directly fromρob,and then Sa can be obtained from S.However,if the chemical reaction parameters(denoted asγ)are insufficiently accurate,bothγand Sc should be corrected.This article proposes a"double correction iterative method"to retrieve Sc and correctγand proves that this iterative method converges.

Review of Real-time Simulation of Power Electronics

Real-time simulation of power electronics has been recognized by the industry as an effective tool for developing power electronic devices and systems.Since there is no energy transfer during the course of the usage,real-time simulation has a lot of advantages in the process of development and experimentation.From the perspective of real-time simulation,this paper focuses on the main problems in modeling accuracy,system bandwidth and stability,limitations on communication interface and energy interface,and the cost of platform construction.Finally,we provide further research directions.

Binary glowworm swarm optimization for unit commitment

This paper proposes a new algorithm—binary glowworm swarm optimization(BGSO)to solve the unit commitment(UC)problem.After a certain quantity of initial feasible solutions is obtained by using the priority list and the decommitment of redundant unit,BGSO is applied to optimize the on/off state of the unit,and the Lambda-iteration method is adopted to solve the economic dispatch problem.In the iterative process,the solutions that do not satisfy all the constraints are adjusted by the correction method.Furthermore,different adjustment techniques such as conversion from cold start to hot start,decommitment of redundant unit,are adopted to avoid falling into local optimal solution and to keep the diversity of the feasible solutions.The proposed BGSO is tested on the power system in the range of 10–140 generating units for a 24-h scheduling period and compared to quantuminspired evolutionary algorithm(QEA),improved binary particle swarm optimization(IBPSO)and mixed integer programming(MIP).Simulated results distinctly show that BGSO is very competent in solving the UC problem in comparison to the previously reported algorithms.

FINITE ELEMENT ANALYSIS OF LAMINAR FLOW AND HEAT TRANSFER IN A BUNDLE OF CYLINDERS

In this paper two dimensional Navier Stokes equations and energy equation governing incompressible laminar flow past a bundle of cylinders were numerically solved by using the finite element method. The velocity correction method was used for time advancement, and spatial discretization was carried out with the Galerkin weighted residual method. Viscous flows past the cylinder banks arranged in in line cylinder bundles and staggered cylinder bundles, coupled with heat transfer, were investigated for pitch diameter ratios of 1.5 and 2.0 and the Reynolds numbers from 50 to 500. Flow structures and heat transfer behavior were discussed. The results obtained agree well with available numerical data.