43results about How to "Improve spatial resolution" patented technology

The application relates to a mode filter, a topology optimization method and an ultrasonic guided wave damage imaging artifact suppression method, relates to the field of ultrasonic detection and imaging, and the topology optimization method comprises the following steps: acquiring initial design parameters, generating an initial population, each chromosome individual in the initial population is coded by a binary logic matrix representing material distribution and meeting design constraints; a population is optimized by using a genetic algorithm, an adaptability function in the genetic algorithm is constructed based on a mode purity index and an energy transmission index, the mode purity index is the ratio of in-plane displacement integral in a transmission area in finite element simulation with the mode filter to total displacement integral, and the energy transmission index is the ratio of relative energy at a transmission end in the finite element simulation with the mode filter to relative energy of an S0 mode; and final optimized material distribution of the mode filter is output. Compared with the prior art, the application has the advantages of effectively suppressing background artifacts, giving consideration to mode regulation performance and substrate structure integrity and the like.

This application discloses a crystal, a rare-earth-doped crystal, its preparation method, and its applications, belonging to the field of crystal materials. The chemical formula of the crystal is [Gd(L)]. 1.5 (DMA)] n Formula I; in Formula I, L is a ligand formed by deprotonated benzene polycarboxylic acid; benzene polycarboxylic acid is selected from any one of terephthalic acid, phthalic acid, and isophthalic acid; DMA is N , N -Dimethylacetamide; n ∞ represents continuous repetition and infinite extension. The crystal is rare-earth doped, and the chemical formula of the rare-earth doped crystal is [Gd]. 1‑x% (L) 1.5 (DMA)] n :x% Ln 3+ Equation II; in Equation II, x = 0.1~10; Ln 3+ Selected from Tb 3+ Or Eu 3+ Any one of the following. It has potential application value in fields such as fluorescent materials, multicolor display materials, long afterglow materials, X-ray detection and imaging materials, radiation detection dosimeters, and optical information encryption. The preparation method used in this application can be easily achieved through solution reaction under relatively mild conditions, with low cost and strong industrial feasibility.

This application provides a method, apparatus, equipment, and storage medium for assessing drought / flood disaster risks in tobacco cultivation. It relates to the field of hydrological technology. The method includes: analyzing the response relationship between drought / flood conditions and actual disasters in different months of the tobacco growing season based on multi-timescale SPI analysis; calculating the distribution sequence of drought / flood conditions with a high probability of causing drought / flood disasters using the Naive Bayes method; constructing a new paradigm for assessing the risk of drought and flood disasters during the tobacco growing season based on drought / flood type classification, Bayesian modeling, and disaster-causing effect quantification; and further establishing a three-dimensional simplified framework drought / flood risk assessment model. This application aims to establish a drought / flood risk assessment model by analyzing the response relationship between drought / flood conditions and actual disasters in different months of the tobacco growing season based on multi-timescale SPI analysis, calculating the distribution sequence of drought / flood conditions with a high probability of causing drought / flood disasters using the Naive Bayes method, and then establishing a drought / flood risk assessment model.

This invention relates to the field of pile foundation testing, and more particularly to a pile foundation testing device and method for civil engineering construction. The device includes a tamping apparatus for testing the pile body and testing components. In this invention, when the pile body tilts, the light spot undergoes vertical and / or lateral displacement on the receiving plane. If the tilt causes the light spot to move upward or downward, it may partially or completely enter the light aperture range of adjacent adjustment apertures. For example, when the top of the pile tilts to the right, due to the clockwise rotation of the transmitting end, the projection point of the 45° beam moves upward on the receiving end plane. The light spot originally illuminating the middle adjustment aperture may simultaneously cover the edge of the light aperture of one or two upper adjustment apertures. If the tilt angle is large, the center of the light spot may completely move out of the original aperture and into the adjacent aperture area, or even cross multiple apertures. Therefore, the tilt amount of the pile body can be determined by monitoring the number of infrared receivers that have sensed infrared radiation.

This invention provides a method and system for direction-of-arrival (DOA) estimation of minimum redundancy linear arrays based on the quantum giant trevally mechanism, belonging to the field of array signal processing. To address the problem of drastic performance degradation of minimum redundancy linear arrays under harsh conditions of impulse noise, this invention proposes using the mean square error of DOA estimation and the probability of successful estimation as core optimization indicators. This method effectively improves the performance degradation problem of arrays under impulse noise by constructing a mapping relationship between the positions of minimum redundancy linear array elements and DOA estimation performance, significantly enhancing the accuracy and robustness of the direction-finding system. The quantum swarm intelligence optimization method designed in this invention provides a promising solution path for the optimization problem through a global search strategy and fitness function design. This invention not only overcomes the applicability and performance bottlenecks of traditional methods under impulse noise environments but also has significant theoretical and practical application value for promoting the engineering application of array technology in high-precision direction-finding.

A method and system for super-resolution reconstruction of infrared images based on convolutional neural networks, relating to the field of electronic digital data processing, is disclosed. The method includes: inputting a low-resolution infrared image into a preset convolutional neural network to obtain multi-layer feature maps of the low-resolution infrared image; scaling the feature maps of different layers in the multi-layer feature maps according to a preset ratio and then stitching and fusing them to obtain a fused feature map; generating feature vectors from the fused feature map using global average pooling, and transforming the feature vectors using a fully connected layer to obtain attention weights; weighting the attention weights with the fused feature map to obtain an enhanced feature map; and performing upsampling reconstruction processing on the enhanced feature map using an upsampling structure, introducing a residual connection structure during the upsampling reconstruction process to generate a high-resolution infrared image. Implementing this method generates high-resolution infrared images with more detail.

This invention discloses a high spatiotemporal resolution millimeter-wave radar imaging system, comprising: a U-shaped antenna array module, wherein the transmitting and receiving antennas are symmetrically distributed in a U-shape, with an element spacing greater than or equal to half a wavelength, adapted to the SAR imaging mechanism, and using an electronic scanning scheme to achieve beam scanning; a radio frequency (RF) module, employing a cyclic one-transmit-multiple-receive mode, completing data acquisition for all transmitting and receiving channels by cyclically switching the transmitting elements to ensure high system refresh rate and spatiotemporal resolution; a signal processing module, receiving multi-channel echo signals, rearranging the echo signals, and using a hardware-parallel accelerated time-domain backpropagation (BP) algorithm for imaging; and a control module, using a microcontroller to control the RF chip of the RF module. This invention verifies the applicability of antenna spacing greater than or equal to half a wavelength in near-field SAR imaging mechanisms, overcomes the dependence of traditional MIMO systems on antenna half-wavelength spacing, achieves high spatiotemporal resolution near-field imaging, and reduces system complexity and cost.

This invention discloses a method for encoding structured light patterns, a structured light projection module, and a three-dimensional measurement system, belonging to the field of optical measurement. Addressing the problem of fixed and inaccurate structured light encoding units in existing methods, this invention provides a method for encoding structured light patterns. The method divides the structured light pattern into multiple types of grid units, forming a grid encoding group with a predetermined number of consecutively arranged grid units. The grid encoding group is encoded based on the type of each grid unit, ensuring that each grid encoding group has a unique and non-repeating code. The code of the grid encoding group maps to its spatial position within the structured light pattern. Decoding includes acquiring the structured light pattern, determining the type of each grid unit, and matching the grid encoding group with the encoding set to obtain its position. This invention achieves flexible, high-density, and highly occlusion-resistant absolute position encoding by uniquely and non-repeating grid encoding groups composed of multiple consecutive grid units, making it suitable for high-precision three-dimensional measurement.

This invention discloses a fine motor imagery decoding system based on PINN and multimodal fusion. The system employs a fine motor imagery decoding method based on PINN and multimodal fusion, and includes a multimodal physiological signal synchronous acquisition module, a data preprocessing and correlation analysis module, a dual-stream decoding network module, a physical manifold constraint layer module, an adaptive gradient balance optimization module, and a classification decision module. By deeply embedding biophysical mechanisms into a deep learning framework, this invention addresses the ill-posedness of the source localization inverse problem, overcomes overfitting in small-sample scenarios, and significantly improves the decoding accuracy, robustness, and physiological interpretability of fine hand movements.

This invention provides a highly efficient photon scanning method and ranging technique, primarily combining rotational scanning and push-broom techniques, and employing a multi-beam emission method. When the aircraft is moving in a straight line, rotational scanning and push-brooming can be simultaneously achieved simply by rotating the light source. A single-photon array detector is used to collect the echo photons reflected from the detected surface by the probe laser. This invention significantly advances the rapid measurement of depth features of large-area targets. Rotational scanning makes the sampled target points denser and improves the spatial resolution of the detection; push-brooming allows the system to measure more target points in the scanning direction, increasing the detection area. This invention has good development potential in fields such as terrain mapping and single-photon 3D imaging.

The present disclosure relates to the technical field of seismic processing, and particularly relates to an interpolation method and device based on prior data, equipment, a storage medium and a computer program. The method comprises the following steps: acquiring seismic data, performing τ-P positive conversion on the seismic data in a preset data processing domain to obtain converted data; selecting unaliasing and preset numerical aliasing data in the converted data, and collecting the selected unaliasing and preset numerical aliasing data into to-be-processed data; performing five-dimensional interpolation on the to-be-processed data to obtain an interpolation result; and performing five-dimensional regularization processing on the seismic data according to the interpolation result to obtain a processing result. The present application can adapt to interpolation processing based on prior data.

This invention discloses a snapshot-type hyperspectral three-dimensional skin imaging method and device, belonging to the field of biomedical imaging technology. The method includes: constructing a snapshot-type spectral-light field coupled imaging system containing a light source and an acquisition model; achieving full-domain skin illumination with broadband illumination; synchronously coupling and acquiring spectral and three-dimensional spatial information through the synergistic action of a lens group and a hyperspectral sensor to generate spectral coupled light field data; establishing a correlation dataset between the spectral coupled light field image and the hyperspectral light field image; training a spectral reconstruction network and solving the data; reconstructing the hyperspectral light field image and obtaining multi-view images through macro-pixel sampling; using a light field super-resolution network to improve image resolution; selecting the optimal skin penetration band image for stereo matching to obtain a target band depth map; and fusing the high-resolution center-view images from each channel with the target band depth map to complete the generation of a hyperspectral three-dimensional image. This invention achieves single-exposure hyperspectral three-dimensional imaging of the superficial skin layer.

The invention discloses a soil body heat flux profile determination method based on an AH-DTS technology, and the method comprises the steps: determining a reference depth at which the soil body temperature is stable along with time through measuring the soil body temperature at different depths and different moments and the corresponding volume water content; calculating the volume heat capacity of the soil body, layering the soil body above the reference depth, and obtaining the average temperature and heat capacity of each layer; and calculating the single-layer heat storage amount and accumulating to obtain the total heat storage amount, and finally obtaining the heat flux profile of the soil body at any target depth. The method comprises the steps of reference depth determination, soil volume heat capacity calculation, soil layer equal division and layered heat storage capacity calculation, heat storage capacity accumulation and heat flux profile solving. The method is easy and convenient to operate and high in spatial resolution, the influence of soil moisture content change on the heat flux can be effectively represented, in-situ continuous measurement of the soil heat flux profile can be achieved, and a reliable technical means is provided for research on the soil moisture and heat transfer law.

This invention relates to the field of power chip testing technology and discloses a power chip testing data processing system and method based on big data. The method performs spatial and temporal correction on raw time-series observation data synchronously acquired by multiple testing devices to form a three-dimensional observation data volume. A virtual sensing grid is constructed within this volume, and a continuously spatially distributed gridded parameter field is generated through interpolation and iterative resampling. The local geometric features of this field are identified and adaptively divided into feature sub-regions, mapped to a graph structure, and a hierarchical feature map is generated through topological agglomeration. This map drives a parameter generation network to synthesize a high-resolution parameter field. The synthesized field is then fused with the original gridded parameter field across scales to obtain an enhanced observation data volume. Physical decoupling and signal reconstruction are performed based on a physical constraint model, and the test results are output. This invention achieves intelligent identification and enhanced reconstruction of the internal features of the test data, improving the resolution and physical reliability of the results.

This invention provides a method for assessing prestress loss in prestressed beams based on digital simulation, relating to the field of concrete structure performance assessment technology. The specific steps include: first, establishing a finite element model and dividing the beam into several sub-regions, calculating the regional sensitivity coefficient by combining the material properties and environmental parameters of each region; second, collecting measured data to invert and calibrate the model, optimizing duct friction and anchor deformation parameters; subsequently, inputting historical environmental and load sequences, calculating the time-varying loss of each sub-region through time-step simulation, and summarizing to generate a comprehensive loss index for the entire region; further, quantifying the loss components caused by concrete shrinkage, creep, and prestressing tendon relaxation through independent simulation; and finally, predicting the loss development trend based on future working conditions, providing a reliable basis for long-term structural performance assessment and maintenance decisions.

The application provides a compressor rotor outlet unsteady flow field distribution measurement and analysis method, comprising: arranging a movable slot extending in the circumferential direction on a compressor casing corresponding to a rotor blade outlet, so that a probe in the movable slot moves in the circumferential direction and the radial direction to perform spatial scanning on the rotor blade outlet; performing lock-in measurement, and controlling the probe to collect dynamic pressure signal data according to set acquisition parameters; and establishing a time axis and a space axis respectively according to dynamic pressure signal data obtained by a probe measuring point in a certain blade height direction of the rotor blade and corresponding circumferential space positions and measurement time sequences, and then establishing a time-space correlation axis; performing interpolation on dynamic pressure signal data on the time-space correlation axis to obtain a flow field reconstruction virtual channel, so as to obtain a change rule of dynamic pressure signal distribution of the rotor blade blade height position with time; and repeating the above steps for measuring points at different radial positions of the rotor blade to obtain a transient two-dimensional distribution of the rotor outlet flow field parameters.

This invention provides a non-invasive stimulation device based on acoustic field modulation using a concave fiber bundle. The device comprises a pulsed laser module, a laser modulation control module, a multi-channel concave fiber bundle, and a photoacoustic transducer coating. The pulsed laser, after processing by the modulation module, is input into the fiber bundle. Its output end is concave, and the photoacoustic transducer coating on the surface converts the laser into ultrasonic waves. The concave structure achieves initial geometric focusing of the ultrasonic waves, forming the geometric focus of the acoustic field. The modulation module has control channels corresponding one-to-one with the fiber transmission channels, allowing for dual-dimensional control of the emission time delay and pulse intensity of each channel's pulsed laser. This enables hardware-free axial and lateral scanning of the focus; intensity control maximizes the focus sound pressure and minimizes the sidelobe sound pressure. This invention significantly improves sound pressure gain and spatial resolution through the combination of concave fiber bundle geometric focusing and a high-efficiency photoacoustic transducer coating.

The application discloses a wide-field magnetic imaging system and a fault detection method based on a diamond nitrogen-vacancy color center. The wide-field magnetic imaging system comprises a laser excitation unit for generating and forming a uniform excitation light field to irradiate a diamond NV color center sensor; a fluorescence collection and imaging unit for collecting fluorescence emitted by the NV color center under excitation and forming a wide-field fluorescence image, comprising an objective lens, a filter and a camera; a signal processing and control unit electrically connected with the laser excitation unit and the fluorescence collection and imaging unit, for controlling the system timing, and processing the intensity information of the wide-field fluorescence image to obtain a two-dimensional distribution map of the measured magnetic field.

The application is suitable for the technical field of agricultural monitoring, and particularly relates to a smart agricultural monitoring method and system.The method comprises the following steps: delimiting a target area requiring disease and pest monitoring, and dividing the target area into a plurality of sections, selecting monitoring points, irradiating the sections by using light source devices previously arranged at the monitoring points, and clustering the light source devices into a fixed group and an adjustable group; starting the light source devices of the fixed group, collecting spectral reflectance data at each monitoring point, wherein each wave band corresponds to a reflectance, drawing a plurality of spectral curves, and identifying the change trend of each spectral curve.The application can adjust the irradiation angle, optimize the light coverage, reduce the monitoring dead angle, meet the monitoring requirements of different leaf orientations and plant structures, greatly improve the monitoring accuracy, and improve the disease and pest control efficiency and intelligent level of crops.

This invention provides a method and system for analyzing the surface shape changes of coatings on ultrathin crystal devices based on finite element analysis, belonging to the field of data processing technology. The method includes: calculating the displacement gradient tensor of subdomain elements to obtain entropy correction coefficients; correcting the amplitude of a uniformly distributed mechanical load based on the entropy correction coefficients to obtain a final uniformly distributed mechanical load; reloading the final uniformly distributed mechanical load into a numerically corrected configuration, repeating the solution and load correction process until the absolute difference of the entropy correction coefficients between two adjacent iterations is less than a preset tolerance, obtaining a final convergent potential field; and extracting the nodal displacement distribution on the substrate surface based on the final convergent potential field to obtain the surface shape change result. This invention can more accurately solve the displacement field and surface shape changes of ultrathin crystal devices after coating, effectively improving the convergence stability and numerical reliability of surface shape analysis.

This invention belongs to the field of spectroscopic instrument technology, and particularly relates to a snapshot-type FP interferometric array spectroscopic imager. The snapshot-type FP interferometric array spectroscopic imager includes: a front imaging mirror group for imaging target light onto a primary image plane and collimating the light emitted from the primary image plane into parallel light; an FP interferometric array for multi-beam interference modulation of the parallel light, outputting modulated light carrying interference information; a microlens array for aperture-splitting imaging of the modulated light to form an interferometric image array; and a relay imaging mirror group for secondary imaging of the interferometric image array and coupling the imaging results to a detector. This invention achieves a good balance between spatial resolution, spectral bandwidth, and spectral resolution, while enabling real-time snapshot detection of rapidly changing targets. This invention eliminates the need for moving mirror scanning, resulting in higher structural stability; and by eliminating the beam splitter, it effectively reduces the difficulty of system assembly and alignment.

The application relates to a space non-continuous deformation monitoring method based on a time-series InSAR technology, and specifically comprises the following steps: firstly, an arc network is constructed based on selected coherent points, and an arc segment solving model without a deformation mode constraint is used to perform time-series solving on all arc segments in the network; then, based on a phase jump mechanism at a non-continuous deformation position, an arc segment mean square error threshold method is used to detect the displacement discontinuous position and divide the coherent network into multiple independent subnetworks; based on the mechanical deformation characteristics of the structure, a center point in each subnetwork is selected as a reference point, and then multi-reference-point and adaptive unwrapping of all subnetworks are completed, and finally, relatively reliable deformation monitoring results are obtained. The method provided by the application can automatically detect the displacement discontinuous position and perform adaptive segmented unwrapping, can overcome the inapplicability of existing time-series InSAR algorithms in space non-continuous deformation monitoring, and improves the reliability of the deformation monitoring results.

The utility model discloses a kind of force-sensitive trigger probe of annular microstructure optical fiber sensing, it is related to precision measurement technical field, to solve the integration of existing probe structure and FBG sensor generally not high, the problem of limited measurement accuracy improvement;A kind of force-sensitive trigger probe of annular microstructure optical fiber sensing, including support unit, support unit is fixedly connected with annular sensitive unit, the central position of annular sensitive unit is provided with trigger probe unit, annular sensitive unit includes several rings, and sensitive beam is arranged between adjacent rings and between trigger probe unit and nearest ring, and optical fiber grating sensing module is arranged on sensitive beam.The utility model technical scheme considers sensitivity and stiffness, through the combination of multi-ring structure and sensitive beam, both ensure overall stiffness, and realize high sensitivity strain detection by sensitive beam.

This invention proposes a resistive tactile imaging method and system based on graph networks and attention mechanisms. The method includes: acquiring measured spatial boundary voltage sequences and measured boundary voltage sequences under pressure using a resistive tactile sensor to be imaged, and calculating the difference to obtain the conductivity distribution, which is then mapped to the graph structure features of the corresponding graph nodes in the finite element simulation model of the sensor; inputting these graph structure features into a preset hybrid graph neural network model to obtain the predicted conductivity values ​​of the corresponding graph nodes, ultimately achieving imaging; wherein each layer of the hybrid graph neural network model employs a hybrid graph convolutional module composed of sequentially connected graph convolutional networks and graph attention networks. This invention combines graph neural networks to extract spatial topological features and utilizes attention mechanisms to enhance local responses, suppress imaging noise and artifacts, achieve robust transfer from simulation to actual measurement, and improve the reconstruction accuracy and spatial resolution of the tactile contact area.

This invention provides a metasurface structure and design method for realizing broadband measurement substrates, relating to the field of micro-spectral imaging technology. It enables spectral imaging devices to achieve higher spatial resolution while maintaining spectral resolution. The metasurface structure includes a metasurface modulation layer, a CMOS image sensor wafer layer, and a signal processing circuit layer. The metasurface modulation layer comprises multiple metasurface modulation units, each independently performing broadband modulation of incident light. The transmission spectrum corresponding to each metasurface modulation unit includes multiple absorption peaks. The CMOS image sensor wafer layer receives the optical signal encoded by the multiple metasurface modulation units of the metasurface modulation layer and converts the optical signal into an electrical signal. The signal processing circuit layer processes the electrical signal from the CMOS image sensor wafer layer to obtain the transmission spectrum corresponding to each metasurface modulation unit.

The application discloses a method for direction finding of unknown coherent sources based on double relaxation solution, and belongs to the technical field of array signal processing. The application calculates the covariance matrix of array receiving signals, and performs eigenvalue decomposition on the covariance matrix, then sorts the eigenvalues to obtain the eigenvectors corresponding to the maximum eigenvalue and the minimum eigenvalue, constructs a dictionary matrix of beam domain steering vectors, and realizes convex optimization solution of the linear relationship between the dictionary matrix and the maximum and minimum eigenvectors by introducing a relaxation factor, and completes accurate calculation of the signal wave direction. The application can provide a high-precision direction finding method for the passive detection field under the condition that there are related sources and the number of sources is unknown.

The present application relates to the technical field of online monitoring, and more particularly to an online spectrum monitoring and early warning device for UV lining curing quality, which comprises a multi-channel UVLED array arranged in a circumferential partition and an axial partition in a lining pipeline, a variable light spot optical mechanism arranged at an out-light end of the multi-channel UVLED array, a reflector cover matched with the variable light spot optical mechanism, a reflector cover angle adjusting mechanism for driving the reflector cover to change the angle, a spectrum acquisition unit for acquiring the spectrum of the inner wall of the lining pipeline, a human-computer interaction unit for displaying information to a construction worker, and a control unit electrically connected with the multi-channel UVLED array in the circumferential partition and the axial partition, the variable light spot optical mechanism, the reflector cover, the reflector cover angle adjusting mechanism, the spectrum acquisition unit and the human-computer interaction unit respectively. The present application can significantly improve the accuracy, uniformity and construction reliability of the lining curing quality.

This disclosure provides a transistor active amplification sensor device and its formation method. The transistor active amplification sensor device includes an optical sensing amplification circuit, which comprises a first complementary field-effect transistor (CFPT), a second complementary field-effect transistor (CFPT), and a sensing unit. The output terminal of the first CFPT is connected to the input terminal of the second CFPT. The first and second CFPTs are monolithically integrated three-dimensional structures. The optical sensing amplification circuit is used for forward signal amplification. The N-type semiconductor layer of the second CFPT and colloidal quantum dots form the sensing unit, which is used to sense the external environment and generate a potential input signal. This disclosure improves the integration density of the transistor active amplification sensor, reduces its size, achieves higher spatial resolution, enhances the sensor's response to weak signals, and improves the performance of the semiconductor structure.

The application belongs to the field of magnetic particle imaging, and particularly relates to an MPI acceleration calibration method based on a system matrix super-resolution network, aiming to solve the problems of long calibration time and poor accuracy of reconstructed MPI images of the existing magnetic particle imaging system matrix reconstruction method. The method comprises the following steps: placing a unit MNP sample in the imaging field of view of a magnetic particle imaging device for sparse calibration scanning, and constructing a system matrix as a low-resolution system matrix; inputting the low-resolution system matrix into a system matrix super-resolution network model after pretreatment, to obtain a super-resolved system matrix row set; performing post-processing on the super-resolved system matrix row set to obtain a high-resolution system matrix; and combining the high-resolution system matrix to solve the magnetic particle concentration distribution of a target object to be imaged, and reconstruct an MPI image of the target object to be imaged. The application can greatly reduce the calibration time while accurately reconstructing the MPI image.

The present application belongs to the field of sewage treatment, and provides a chemical industry park groundwater pollution tracing method based on geophysical exploration cooperation, which uses geophysical methods to detect target areas in the chemical industry park, and obtains underground electrical parameter data; based on the spatial form characteristics and physical parameter gradient change law of the resistivity anomaly area, a layout scheme of groundwater chemical sampling target points is determined; groundwater samples are collected according to the layout scheme, and chemical detection analysis is performed; the pollutant data obtained by the chemical detection analysis is subjected to spatial position corresponding and correlation analysis with the resistivity three-dimensional model; the pollutant concentration spatial distribution map, the park enterprise layout and the regional hydrogeological conditions are comprehensively considered to identify the pollutant migration path, lock the pollution source and output the pollution tracing result. The method improves the accuracy of pollution tracing.