2296 results about "Geodesy" patented technology

Disclosed is a three-dimensional deformation field modeling method fusing global navigation satellite system (GNSS) / synthetic aperture radar interferometry (InSAR) observation data. The three-dimensional deformation field modeling method includes: five parts of solving an InSAR line-of-sight velocity field, unifying a spatial reference datum of the GNSS data and the InSAR data, verifying precision of a result of the InSAR data, solving a three-dimensional deformation field having high precision based on a least square method, and analyzing and verifying a result. According to the method, the InSAR data is processed according to a high-precision time series method, a long-wavelength error of the InSAR data is corrected on the basis of a difference of two data spatial reference datums, and the accuracy of the result is ensured through verification of internal coincidence precision.

The invention provides a road slope data monitoring method under a complex geological condition, and relates to the technical field of road slope monitoring, and the method comprises the steps: carrying out the temperature compensation of a change rate, analyzing the thermal expansion characteristic of a material, calculating the deformation caused by the temperature change, correcting the monitored side length, and generating a temperature compensation side length; according to the temperature compensation side length and the surrounding rock rheological equation, correcting the two adjacent side lengths through the strain difference of the adjacent sides so as to obtain the final side length; combining the final side length with the reference interior angle to generate a topological structure; based on the topological structure, constructing a three-dimensional deformation index matrix containing a diagonal displacement difference, an adjacent side curvature ratio and an interior angle variation coefficient; and inputting the index matrix, the support structure strain space-time matrix and the rock mass fragmentation tensor into a space-time attention graph neural network together, and outputting a deformation evolution trend including a diagonal displacement difference, an adjacent side curvature ratio and an interior angle variation coefficient. According to the invention, high-precision real-time monitoring can be realized.

The invention discloses a Doppler and pseudo-range fusion positioning method and system for low-orbit satellite conduction signals, and relates to satellite positioning. The method comprises the steps that a Doppler and pseudo-range observation equation is established; performing least square residual calculation on the calculated Doppler and pseudo-range theoretical values and the observed Doppler and pseudo-range measured values, and taking a grid center point coordinate corresponding to a two-norm minimum value as a receiver position initial value; performing Taylor series expansion on the Doppler and pseudo-range observation equation, and representing a state updating matrix by using the established state transition matrix and receiver position, clock error and frequency offset correction; the inverse matrix of the observation noise covariance is used as a weight matrix to be distributed to the state transition matrix; and iterating the initial position value of the receiver, updating the state updating matrix, and stopping iteration until the two-norm of the residual solution is smaller than a set threshold value, so as to obtain the position, clock error and frequency offset results of the receiver. According to the method, the low-orbit satellite positioning algorithm model can be optimized, the convergence time is shortened, and the positioning precision is improved.

The invention discloses a single-epoch GNSS-oriented PAR ambiguity precision factor determination method, and relates to the technical field of satellite navigation and positioning. Comprising the following steps: dividing pseudo-range and carrier wave observed quantities in a GNSS double-difference mathematical model into two groups, and constructing a function relationship between a floating point solution # imgabs1 # and a floating point solution # imgabs2 # variance-covariance matrix of an ambiguity vector # imgabs0 #; constructing a single-frequency single-epoch single-system P-ADOP expression and establishing a shrinkage factor k de-inversion transformation analytic expression, and simultaneously obtaining a single-frequency single-epoch single-system GNSS non-inverse transformation P-ADOP approximation formula; and performing analogy derivation on the single-frequency single-epoch single-system non-inverse transformation P-ADOP approximate formula to obtain a single-frequency multi-system or multi-frequency single-epoch GNSS non-inverse transformation P-ADOP general approximate formula. The method has important theoretical value and practical guiding significance for GNSS optimal ambiguity subset rapid selection and positioning and partial ambiguity fixing success rate influence mechanism analysis.

The invention provides a deformation rate difference-based time sequence InSAR landslide deformation monitoring fusion method, which comprises the following steps of: acquiring and preprocessing a synthetic aperture radar image of a monitoring area to obtain processed SAR image set data, then acquiring a linear deformation rate by using Stacking-InSAR, extracting PS point high-precision deformation speed by using PS-InSAR, and acquiring a nonlinear deformation rate by using SBAS-InSAR; coordinates are generated according to the deformation rate of each technical data point in combination with the orbit and geometric information, and creep, acceleration and temporary sliding deformation areas are divided by using ArcGIS and K-Means clustering; then, SBAS-InSAR processing is carried out on the newly-added image in the regional boundary, and dynamic deformation information is obtained; and finally, carrying out weighted averaging on the deformation rates of the three monitoring technologies. The landslide deformation monitoring precision under complex terrains and environments is effectively improved, errors caused by vector monitoring point density and deformation rate problems in single interference measurement are overcome, and powerful scientific support is provided by accurately monitoring and analyzing landslide evolution, determining the range, predicting the trend and early warning.

The invention discloses an InSAR (Interferometric Synthetic Aperture Radar) and GNSS (Global Navigation Satellite System) robust adaptive fusion method for high-precision three-dimensional surface deformation monitoring, which belongs to the technical field of space geodetic survey and remote sensing, and comprises the following steps of: constructing a time sequence deformation field based on a small baseline set time sequence InSAR processing technology, and resampling GNSS observation data to a grid consistent with the InSAR by using Kriging spatial interpolation; a Helmert variance component estimation method is adopted to adaptively estimate variance components of InSAR and GNSS data, and reasonable weight fixing of a multi-source observation value is achieved; iterative reweighted least square estimation is introduced, outlier influences in various observation data are dynamically restrained through a Tukey double-weight function, robust optimization of a fusion model is achieved, and therefore a high-precision three-dimensional deformation field is reconstructed. The method comprises the steps of Helmert weight fixing, IRLS robust, variance component updating and a closed-loop feedback mechanism of weight matrix optimization.

The invention relates to the technical field of pole and tower inclination monitoring, and provides a pole and tower inclination monitoring method and system based on point cloud rapid reconstruction, and the method comprises the following steps: obtaining original point cloud data of a pole and tower through a hierarchical scanning and adaptive sampling density method; performing down-sampling, noise filtering and preliminary classification on the original point cloud data to obtain tower structure point cloud data; extracting tower geometric features based on the tower structure point cloud data, and constructing a tower skeleton three-dimensional model through the tower geometric features; calculating the inclination angle of each part of the tower based on the central axis of the tower skeleton three-dimensional model; and carrying out deformation analysis and deformation prediction on the tower based on the inclination angle of each part of the tower, and generating a tower inclination monitoring report. According to the method, the whole-process deformation monitoring from data acquisition, processing, modeling to analysis and prediction is realized, and the tower inclination monitoring efficiency is improved.

The invention discloses a PPP-RTK time delay processing method and device in conversion from SSR to OSR, and relates to the technical field of global navigation satellite system positioning, and the method comprises the steps: converting orbit, clock error, ionosphere and troposphere delay SSR parameters provided by a server into OSR virtual observation values through geometric distance calculation and projection function mapping; based on a constant velocity model and a static random process, predicting a correction amount of a future moment, and correcting an OSR virtual observation value through state transition matrix compensation time delay; combining server side parameters and user side observation data, recursively generating an error covariance matrix, dynamically adjusting a state prediction weight, and obtaining covariance information through adaptive filtering optimization positioning calculation; and realizing ambiguity rapid fixation and centimeter-level positioning by using the corrected OSR observation value and covariance information. According to the method, real-time centimeter-level positioning in a high-delay scene is realized through a dynamic correction value prediction method and a recursive covariance estimation method.

The invention provides a synthetic aperture radar sea surface flow velocity inversion method and device based on machine learning. The method comprises the following steps: acquiring an original synthetic aperture radar image, and OCN data, mixed coordinate ocean model data and NDBC buoy data corresponding to the original synthetic aperture radar image; calculating a sea wave spectrum of the synthetic aperture radar image by using a PFSM algorithm, and exploring a relationship between the flow velocity and synthetic aperture radar image parameters; a synthetic aperture radar sea surface flow velocity inversion algorithm based on machine learning is constructed; and verifying the accuracy of the inversion algorithm. The invention provides an effective method for SAR sea surface flow velocity inversion.

The invention provides a satellite-borne SAR ocean image geometric positioning precision improving method and system, and the method comprises the steps: collecting AIS data at the same imaging moment as an SAR satellite and in an imaging region range, correcting the orbit position of the satellite-borne SAR satellite according to the actual imaging moment of a ship, and fusing the AIS data with the SAR satellite data, thereby improving the geometric positioning precision of the satellite-borne SAR ocean image. The method comprises the following steps: optimizing a ship to obtain a matching point pair of an SAR image position and an AIS position of the ship, correcting an RD model error by taking the AIS position as a control point according to an optimized matching point pair result and an orbit position of a satellite-borne SAR satellite corresponding to an actual imaging moment of each ship, and performing iterative solution again to obtain new geographic coordinate information of each image pixel value, finally, a system geometric correction image product with improved positioning precision is obtained, and the problem that in the prior art, ship positioning precision is not high only depending on SAR images obtained by SAR satellites is solved.

The invention belongs to the field of geophysical monitoring, and discloses an earth crust deformation monitoring method and system based on satellite radar interference measurement, electronic equipment and a storage medium, and the method comprises the steps: obtaining a measurement image of a satellite radar for the earth surface, and carrying out the image registration of the measurement image, and obtaining a registered image; pre-processing the registered image to obtain a pre-processed image; performing interference processing on the preprocessed image to generate a differential interferogram containing a terrain phase and a deformation phase; phase unwrapping and terrain phase correction are carried out on the differential interferogram, a deformation phase is separated out, and earth surface deformation information is obtained; and performing time sequence analysis on the earth surface deformation information of the continuous time to obtain a change trend of earth crust deformation. According to the invention, based on satellite radar interference measurement, a large area can be covered in a short time, and effective monitoring can be carried out on wide plains, complex mountainous regions and remote regions.

The invention belongs to the technical field of marine surveying and mapping, and discloses an unmanned aerial vehicle laser / shipborne multi-beam data land and water topographic mapping method, which comprises the following steps of: coarsely extracting a data overlapping region through a virtual grid, finely extracting the overlapping region by using an alpha-shape algorithm, and obtaining boundary points; removing edge triangles in an irregular triangulation network constructed by the unmanned aerial vehicle laser sounding data and the shipborne multi-beam sounding data; calculating a vertex normal vector, taking a vertex tangent plane of the multi-beam triangular irregular network as a registration element, carrying out weight distribution according to the distance from a point to the vertex, and carrying out coarse registration through a point-to-plane ICP registration algorithm; and constructing a precise registration model based on a point-surface distance, and eliminating the elevation difference of the overlapping region by minimizing the normal distance from the unmanned aerial vehicle laser sounding point to the multi-beam sounding triangular surface element. By adopting the point cloud fusion model for processing, the phenomenon that different source point cloud elevations in the overlapping area are inconsistent is effectively improved, and underwater topographic data fusion is realized.

The invention discloses an InSAR (Interferometric Synthetic Aperture Radar) deformation inversion method for lifting rail offset correction, and mainly relates to the technical field of surface deformation monitoring. In order to solve the problem that the deformation monitoring quantity of the existing multi-dimensional small baseline set formed aperture radar interference is obviously lower than the actual deformation quantity, a time sequence deformation observation equation is established by using multi-orbit synthetic aperture radar data, an ascending-orbit small baseline set InSAR settlement monitoring point cloud and a descending-orbit small baseline set InSAR settlement monitoring point cloud are calculated respectively, point cloud registration is carried out, and the settlement monitoring point cloud of the ascending-orbit small baseline set InSAR settlement monitoring point cloud is calculated to obtain the settlement monitoring point cloud of the descending-orbit small baseline set InSAR settlement monitoring point cloud. Calculating the offset of the monitoring result of the lifting rail in the vertical, east-west and south-north directions; calculating to obtain an ascending orbit incidence angle correction number and a descending orbit incidence angle correction number of the synthetic aperture radar image by using the offsets of the ascending orbit monitoring result in the vertical direction, the east-west direction and the south-north direction; and finally, constructing a multi-dimensional small baseline set InSAR deformation inversion model containing lifting rail offset correction, solving to obtain a corrected deformation monitoring quantity, and improving the deformation monitoring precision.

The invention provides a global magnetic disturbance compensation method in aviation magnetic survey. The method comprises the following steps: S1, performing high-altitude calibration flight through a motion platform, and acquiring total magnetic field data based on detection data of an optical pump magnetometer; s2, obtaining geomagnetic background gradient information and obtaining geomagnetic field data based on the geomagnetic background gradient information; s3, obtaining motion platform magnetic interference data according to the total magnetic field data and the geomagnetic field data; s4, motion platform attitude information, motion platform current information and fluxgate magnetometer detection data are acquired; s5, obtaining a motion platform magnetic interference coefficient according to the motion platform magnetic interference data, the motion platform attitude information, the motion platform current information and the fluxgate magnetometer detection data; and S6, according to the geomagnetic background gradient information and the motion platform magnetic interference coefficient, performing data compensation on the total magnetic field data in real time during actual aviation magnetic measurement.

The invention relates to the technical field of coal mine goaf exploration, in particular to a goaf combined detection and identification method fusing gravity, earthquake and transient electromagnetism, which comprises the following steps of: (1) arranging an earthquake, transient electromagnetism and gravity method combined observation system according to a detection geological task and a field condition of a coal goaf; (2) seismic data, transient electromagnetic data and gravity data are collected through the combined observation system; (3) performing independent inversion on the gravity data, the transient electromagnetic data and the seismic data to obtain density distribution, resistivity distribution and wave velocity distribution of the goaf; (4) carrying out a joint inversion method based on improved cross gradient seismic data, transient electromagnetic data and gravity data; (5) obtaining a joint inversion resistivity profile map, a joint inversion wave velocity profile map and a joint inversion density profile map of the goaf through joint inversion of the seismic data, the transient electromagnetic data and the gravity data of the goaf; and identifying development characteristics of the goaf.

The invention discloses a data processing method and system for micro-motion detection array station distribution, and relates to the technical field of micro-motion detection, and the method comprises the steps: according to the actually measured starting point coordinates and the azimuth angle of a measuring line, setting the measuring point distance, and the radius and number of nested triangle circumcircles, generating field station coordinates in batches; collecting and preprocessing station waveform data; calculating a spatial autocorrelation coefficient of any center station and a non-center station; dBSCAN clustering is carried out based on the station distance; performing spatial averaging on the spatial autocorrelation coefficient according to a clustering result; fitting the averaged spatial autocorrelation coefficient and a zero-order Bessel function to obtain a frequency dispersion spectrum, picking up a surface wave phase velocity in the frequency dispersion spectrum, and calculating an apparent shear wave velocity; inverting the shear wave velocity of the waveform data by adopting a self-adaptive gradient inversion strategy; and generating a velocity comprehensive profile of a depth domain and a frequency domain according to the surface wave phase velocity, the apparent shear wave velocity and the inverted shear wave velocity. The invention provides a new method for station distribution to data processing for micro-motion detection.

The invention relates to the technical field of reservoir storage capacity measurement, in particular to a reservoir storage capacity measurement method and system based on remote sensing and point cloud data, and the method comprises the following steps: extracting a terrain surface point set based on input remote sensing image data and point cloud data, and calculating a three-dimensional coordinate value, a curvature value and a point density value of each point; and comparing and classifying curvature and density data in combination with the change of the reservoir boundary. According to the method, through comprehensive application of remote sensing image data and point cloud data, three-dimensional coordinates, curvature values and point density values are adopted to accurately extract topographic features, scientific division and feature description of sub-regions are achieved, refined three-dimensional grid parameters are generated through fitting and interpolation calculation, and the accuracy of the three-dimensional grid parameters is improved. The terrain modeling precision and the space geometry description capability of a complex area are improved, boundary node coordinates and normal vectors are dynamically adjusted, the grid continuity is optimized, the integrity of a boundary transition area is ensured, and the boundary complex volume is accurately described by adopting unit volume subdivision and boundary correction.

The invention provides a building deformation monitoring method and system based on a multi-temporal InSAR technology, and relates to the technical field of building monitoring, and the method comprises the steps: obtaining an initial difference interferogram according to a plurality of SAR images and DEM data; resolving the initial differential interferogram by adopting a multi-temporal InSAR technology to obtain an initial terrain error and an initial time sequence deformation quantity; according to the newly-added SAR image and the DEM data, obtaining a newly-added differential interferogram; according to the previous topographic error, the previous time sequence deformation quantity and the newly added differential interferogram, a multi-temporal InSAR technology and a sequential adjustment method are adopted for calculation, and a current topographic error and a current time sequence deformation quantity are obtained; according to the current time sequence deformation quantity, the horizontal linear distance between the multiple deformation points and the radar wave incidence angle are combined, and the inclination angle of the building is obtained; and monitoring the deformation condition of the building according to the inclination angle of the building and the inclination angle threshold. The analysis complexity is reduced, the analysis efficiency is relatively high, and rapid and continuous deformation monitoring of the building can be realized.

The invention provides an ocean buoy data real-time acquisition method based on low earth orbit satellite communication, and belongs to the technical field of electrical digital data processing.The ocean buoy data real-time acquisition method comprises the steps that firstly, a low earth orbit satellite network communication link is built, data characteristics are analyzed through a compression conversion gain index model, and a low earth orbit satellite network communication link is established; and establishing a multi-stage compression mechanism. The system adopts an optimization equation set to dynamically adjust compression parameters, the optimization equation set comprises a resource consumption equation, a compression quality equation and a transmission efficiency equation, and self-adaptive compression of data is achieved. And meanwhile, through a hierarchical algorithm selection mechanism and a parameter optimization strategy and in combination with real-time feedback control, continuous optimization of the compression efficiency is ensured. According to the compression scheme based on multi-level optimization, the technical problem of low compression efficiency caused by poor algorithm adaptability in the ocean buoy data compression process in the prior art is effectively solved.

The invention belongs to the field of oil-gas exploration seismic data processing, and relates to a first-motion wave travel time chromatography velocity inversion method, system and device and a medium. According to the method, the first arrival travel time and the grid unit travel time are acquired, and the minimum grid unit travel time is selected; an initial velocity field closer to a real structure is established through micro-logging constraint so as to replace a conventionally used gradient velocity field. The initial velocity field and the first arrival travel time are combined to construct the target function, the slowness information and the regularization information are added in the target function, the multi-information constrained target function is obtained, and the inversion stability and the calculation precision are improved. And adopting a Gaussian-Newton method to minimize the target function of the multi-information constraint, solving the minimized target function of the multi-information constraint, and carrying out iterative updating on a solving result to obtain an inversion velocity field. Explicit calculation and sensitive kernel function storage are not needed, the method is more suitable for large-scale data processing, and the calculation time of speed inversion and the requirement for storage space are reduced.

The invention relates to the technical field of seismology, in particular to a seismic oscillation field stochastic simulation method based on spatial correlation, which is composed of a single-field-point seismic oscillation model and a spatial correlation model, and specifically comprises the following steps: step 1, building a single-field-point seismic oscillation time-frequency power spectrum model; 2, building a prediction equation of time-frequency power spectrum parameters; step 3, converting a time-frequency power spectrum to a time history; 4, building a spatial correlation model of the time-frequency power spectrum parameters; according to the method, the simulation speed of the seismic motion field is greatly improved, meanwhile, spatial correlation of model parameters is introduced, correlation of time-frequency characteristics between field points in space is well represented, and the seismic motion field simulated through the method can conform to physical cognition.

The invention discloses a millimeter wave radar and range finder golf sports monitoring system and method and a medium, and relates to the field of sports monitoring, the system comprises a handheld integrated terminal, a laser ranging module, a millimeter wave radar module, a processing unit and a display unit are integrated in the handheld integrated terminal, and a first communication module is connected with the processing unit; the data interaction module is used for performing data interaction with an external mobile terminal; the mobile terminal comprises a second communication module, a trajectory prediction module and a recommendation decision module, and is used for comparing the predicted flight trajectory with a target trajectory preset by a user and generating a swing suggestion in a training mode; and in the motion mode, the central processing unit is used for planning and recommending a ball hitting strategy based on the expected drop point distance input by the user and the real-time environment data measured by the handheld integrated terminal in combination with historical training data of the user. According to the invention, accurate monitoring and dynamic analysis of the golf ball flight path are realized.

The invention discloses a high-precision three-dimensional first-motion wave travel time chromatography velocity inversion method and device, and relates to the field of seismic exploration. A multi-template windward difference algorithm and an undulating surface narrow-band continuation strategy are adopted to determine three-dimensional global forward first-motion wave propagation travel time; ray tracing is carried out according to the three-dimensional global forward first-motion wave propagation travel time, and a first-motion wave propagation path and a Frechet differential quotient matrix are obtained; establishing a three-dimensional initial inversion velocity field constrained by micro-logging information; constructing a tomographic inversion equation based on the apparent slowness constraint and the regularization constraint according to the travel time residual error and the Frechet differential quotient matrix; and solving the tomographic inversion equation to obtain a tomographic inversion slowness solution, and converting the slowness solution into a velocity solution. The problem of non-uniform travel time delay of seismic wave propagation of a complex near-surface structure can be solved, a more accurate superficial layer velocity structure can be obtained, and the quality of seismic data processing and imaging inversion can be improved.

To eliminate the influence of station position fluctuations in a satellite navigation system. [Solution] In satellite navigation systems, including the US GPS and Japan's Quasi-Zenith Satellite System, when multiple reference stations are installed to form a wide-area differential correction system (WADGPS), previously, position errors caused by fluctuations in station position would remain in the position information calculated by the user receiver.However, by adding a correction amount to the correction information provided by WADGPS in advance to cancel this, the position errors caused by fluctuations in station position at the user station can be reduced without changing the calculation processing content at the WADGPS user station, thereby improving the accuracy of the position information calculated by the user station.

The invention is suitable for the technical field of aeromagnetic data processing, and provides an aeromagnetic data compensation method and system for a multi-rotor unmanned aerial vehicle, and the method comprises the following steps: obtaining the aeromagnetic data of the multi-rotor unmanned aerial vehicle; the aeromagnetic data comprises three-component fluxgate data; determining total geomagnetic field data according to the three-component fluxgate data; according to the three-component fluxgate data and the total geomagnetic field data, calculating to obtain a direction cosine matrix; based on a preset PINN neural network, taking the direction cosine matrix and the total geomagnetic field data as input values, and outputting a magnetic interference prediction value; and compensating the aeromagnetic data according to the magnetic interference predicted value. According to the method, the PINN neural network is combined with physical information contained in the data, so that the fitting compensation capability is relatively high, meanwhile, the method does not need to depend on a large amount of data support, and the compensation precision of the aeromagnetic data is relatively high.

The invention relates to a ship berthing position detection method and system based on shore end laser radars, and the method comprises the steps: deploying a plurality of laser radars at a shore end, and collecting original point clouds in parallel through an asynchronous thread model; constructing a rigid transformation matrix based on the radar installation attitude parameters, and further generating fusion point cloud data under a unified world coordinate system; then extracting 3D geometric edge points in a preset range of the front edge of the wharf as near-shore point cloud data, and further calculating the minimum ship-shore distance; based on the ship bow direction and a preset segmentation threshold value, the coastal point cloud data are segmented into ship bow area point cloud and ship stern area point cloud in the x-axis direction, the coordinate mean value of the ship bow area point cloud and the ship stern area point cloud is calculated and compared with an expected position in a wharf operation system, and berthing accuracy is judged; and the drift distance is calculated in real time, the drift distance, the attitude angle and the change rate of the attitude angle are compared with corresponding preset threshold values respectively for graded early warning and drift alarm, and high-precision ship berthing position dynamic detection is achieved.

The invention provides a measurement mark-free section deformation resolving method for a high-speed rail tunnel, and relates to the technical field of tunnel detection. The method comprises the following steps: acquiring point cloud data of a high-speed rail tunnel, and carrying out mileage alignment on the point cloud data and datum point cloud data; performing section registration on the point cloud data and the reference point cloud data; establishing a corresponding relation of two periods of point cloud data through a nearest neighbor point of any point position in the datum point cloud data, and setting a displacement difference threshold value and an intensity grading threshold value of a point pair based on the corresponding relation; performing clustering analysis on the point clouds exceeding the displacement difference threshold value or the intensity grading threshold value based on a target clustering algorithm, and generating a potential deformation area continuous mask; the potential deformation area continuous mask represents the area where the high-speed rail tunnel has diseases. According to the method, the detection efficiency and the detection accuracy of high-speed rail tunnel deformation detection can be improved.

The invention discloses a radar distributed scatterer interference measurement method and device for a high and steep area. The method comprises the following steps: selecting a main image and a non-main image; calculating differential interference data; screening PS candidate pixels of which the amplitude deviation indexes are smaller than an amplitude deviation threshold value; high-coherence DS candidate pixels are screened out; constructing a distance-constrained Delaunay triangulation network for fusing the PS candidate pixels and the high-coherence DS candidate pixels; resolving a differential residual elevation, a differential deformation rate, a set coherence coefficient of a connected edge and a spatial differential residual phase; calculating residual terrain, linear deformation, deformation rate and residual phase after unwrapping; calculating an atmospheric orbit phase, nonlinear deformation and a noise phase; and obtaining time sequence deformation based on the linear deformation and the nonlinear deformation. According to the method, the number of monitoring points of the orbit descending InSAR deformation monitoring result is large, the coverage range of the monitoring result is widest, and MAE and RMSE are small.

The invention relates to the technical field of global satellite navigation systems, in particular to a network RTK (Real Time Kinematic) method for generating a virtual reference station based on a PPP (Point-to-Point Protocol) SSR (Simple Sequence Repeat) correction, which comprises the following steps of: calculating a moving station approximate coordinate of a target moving station; an observation value error of each target satellite is corrected through the target server by using the rough coordinate of the moving station and the PPP SSR information of the plurality of target satellites, and a position coordinate and an observation value of a virtual reference station are generated; on the basis of an RTCM protocol, the target server side encodes the position coordinates and the observation values to obtain OSR information corresponding to the position coordinates and the observation values, and broadcasts the OSR information to the target moving station; and based on the OSR information and a double-difference relative positioning strategy, enabling the target moving station to obtain relative position information between the target moving station and the virtual reference station so as to determine position coordinates of the target moving station. Therefore, the problems of slow convergence, limited coverage range, poor receiver compatibility and the like of the existing navigation positioning technology are solved.

The invention relates to the technical field of ocean engineering digital twinning and three-dimensional visualization, and provides a method and system for generating dredging terrain twinning based on a right-angle triangular mesh, and the method comprises the steps: obtaining regular grid elevation unified to an engineering plane coordinate system or preprocessed sampling height data, determining a covered space range, reading the space range and density parameters, and calculating the total row number and the total column number of the grids based on the space range and the density parameters; generating regularly arranged vertexes in the covered space, and constructing a topological structure of an isosceles right triangle grid; dividing the global grid into a plurality of Chunk blocks according to a fixed size by taking a vertex as a unit, establishing a triangular index mapping relationship between global rows and columns and the Chunk blocks, filling three vertexes of an isosceles right triangle and elevation / height data in each Chunk block in parallel, and generating local VBO, IBO and AABB; and performing batch rendering and coloring according to a height color table on all Chunk blocks to realize visualization of the terrain.