Linear Inversion Algorithm Research Articles

A Concrete Core Void Imaging Approach and Parameter Analysis of Concrete-Filled Steel Tube Members Using Travel Time Tomography: Multi-Physics Simulations and Experimental Studies.

Concrete-filled steel tube (CFST) members have been widely used in civil engineering due to their advanced mechanical properties. However, internal defects such as the concrete core voids and interface debonding in CFST structures are likely to weaken their load-carrying capacity and stiffness, which affects the safety and serviceability. Visualizing the inner defects of the concrete cores in CFST members is a critical requirement and a challenging task due to the obvious difference in the material mechanical parameters of the concrete core and steel tube in CFST members. In this study, a curved ray theory-based travel time tomography (TTT) with a least square iterative linear inversion algorithm is first introduced to quantitatively identify and visualize the sizes and positions of the concrete core voids in CFST members. Secondly, a numerical investigation of the influence of different parameters on the inversion algorithm for the defect imaging of CFST members, including the effects of the model weighting matrix, weighting factor and grid size on the void's imaging quality and accuracy, is carried out. Finally, an experimental study on six CFST specimens with mimicked concrete core void defects is performed in a laboratory and the mimicked defects are visualized. The results demonstrate that TTT can identify the sizes and positions of the concrete core void defects in CFST members efficiently with the use of optimal parameters.

Robust breast quantitative susceptibility mapping in the presence of silicone.

To (a) develop a preconditioned water-fat-silicone total field inversion (wfsTFI) algorithm that directly estimates the susceptibility map from complex multi-echo data in the breast in the presence of silicone and to (b) evaluate the performance of wfsTFI for breast quantitative susceptibility mapping (QSM) in silico and in vivo in comparison with formerly proposed methods. Numerical simulations and in vivo multi-echo gradient echo breast measurements were performed to compare wfsTFI to a previously proposed field map-based linear total field inversion algorithm (lTFI) with and without the consideration of the chemical shift of silicone in the field map estimation step. Specifically, a simulation based on an in vivo scan and data from five patients were included in the analysis. In the simulation, wfsTFI is able to significantly decrease the normalized root mean square error from lTFI without (4.46) and with (1.77) the consideration of the chemical shift of silicone to 0.68. Both the in silico and in vivo wfsTFI susceptibility maps show reduced shadowing artifacts in local tissue adjacent to silicone, reduced streaking artifacts and no erroneous single voxels of diamagnetic susceptibility in proximity to silicone. The proposed wfsTFI method can automatically distinguish between subjects with and without silicone. Furthermore wfsTFI accounts for the presence of silicone in the QSM dipole inversion and allows for the robust estimation of susceptibility in proximity to silicone breast implants and hence allows the visualization of structures that would otherwise be dominated by artifacts on susceptibility maps.

Optimization of Lai Estimation Method Based on Smartphones with Fisheye Lens

Leaf area index (LAI) is an important biological factor reflecting vegetation growth and forest ecosystem. LAI can be used to obtain plant health status, carbon cycle, and surrounding ecological environment effectively. In this study, the smartphone was equipped with a fisheye lens, and the optimization method was used to estimate LAI, which was compared with digital hemispherical photography (DHP) to investigate the possibility of the new method for LAI estimation. The hemispherical image was divided into blocks, and the optimized Otsu method was used for algorithm segmentation, which can effectively distinguish vegetation from the sky. Concurrently, when the gap fraction inversion LAI was performed, the linear inversion algorithm was improved based on single-angle inversion, and the LAI was obtained by inversion through the linear fitting of the mul-tiangle gap fraction. The experimental sample was located in Olympic National Forest Park in Beijing. Three coniferous mixed forests and three broadleaved forests were selected from the experimental sample. LAI measurements from smartphones were compared with those from DHP. In the samples for mixed coniferous forests, the values for coefficients of determination R^2 were 0.835, 0.802, and 0.809, and root mean square errors (REMS) were 0.137, 0.120, and 0.147. For the broadleaf forest samples, the values for R² were 0.629, 0.679, and 0.758, and REMS were 0.144, 0.135, and 0.137. The R^2 and RMES for the overall data was 0.810 and 0.134, respectively, and a good agreement between the LAI measurements from the proposed method and those from the DHP supports an accurate estimation. The results show that the use of a fisheye lens on a smartphone can effectively and accurately obtain tree canopy LAI. This provides a fast and effective new method to measure LAI of forest vegetation near the ground, which is of great significance for studying the interaction between plant growth status, ecological environment, and phenological changes.

Video-Based Nearshore Bathymetric Inversion on a Geologically Constrained Mesotidal Beach during Storm Events

Although geologically constrained sandy beaches are ubiquitous along wave-exposed coasts, there is still a limited understanding of their morphological response, particularly under storm conditions, which is mainly due to a critical lack of nearshore bathymetry observations. This paper examines the potential to derive bathymetries from video imagery under challenging wave conditions in order to investigate headland control on morphological beach response. For this purpose, a video-based linear depth inversion algorithm is applied to three consecutive weeks of frames collected during daylight hours from a single fixed camera located at La Petite Chambre d’Amour beach (Anglet, SW France). Video-derived bathymetries are compared against in situ topo-bathymetric surveys carried out at the beginning and end of the field experiment in order to assess the performance of the bathymetric estimates. The results show that the rates of accretion/erosion within the surf zone are strongly influenced by the headland, whereas the beach morphological response can be classified into three main regimes depending on the angle of wave incidence θp: (1) under deflection configuration (θp>0°), the alongshore sediment transport was trapped at the updrift side of the headland, promoting sand accretion. (2) Under shadowed configuration (θp<0°), the interruption of the longshore current drove a deficit of sand supply at the downdrift side of the headland, leading to an overall erosion in the surf zone. (3) Under shore-normal configuration (θp=0°), rip channels developed, and up-state beach transition was observed. A comparison between video-derived bathymetries and surveys shows an overall root mean square error (RMSE) around 0.49 to 0.57 m with a bias ranging between −0.36 and −0.29 m. The results show that video-derived bathymetries can provide new insight into the morphological change driven by storm events. The combination of such inferred bathymetry with video-derived surface current data is discussed, showing great potential to address the coupled morphodynamics system under time-varying wave conditions.

Ultrasound Computed Tomography.

This chapter presents theoretical, numerical, and experimental frameworks for the use of Ultrasound Computed Tomography (USCT) for cortical bone tissue imaging. Most of the research conducted on this topic concerns adult bone, although some work presented in this chapter is specific to the study of child bone. USCT is recognized as a powerful method for soft tissue imaging. In bone imaging, the difficulties arise from the very high impedance contrast between tissues which alters the propagation of the ultrasonic waves and limits the linear inversion algorithms used. Solutions consist in optimally assessing non-linear effects in an iterative approach aiming at local linearization. When the problem can be reduced to the study of a fluid-like cavity buried in an elastic cylinder surrounded by water, the signal processing and/or compound algorithms can be added as an extension to the linear algorithms. The main limitation of these methods is the heavy experimental costs involved. We have then suggested the introduction of purely numerical non-linear full-waveform inversion algorithms. The performances and the limitations of these linear and non-linear methods applied to cortical bone tissue imaging problems are overviewed and discussed.

An Algorithm to Derive Temperature and Humidity Profile Changes Using Spatially and Temporally Averaged Spectral Radiance Differences

AbstractA linear inversion algorithm to derive changes of surface skin temperature and atmospheric temperature and specific humidity vertical profiles using spatially and temporally averaged spectral radiance differences is developed. The algorithm uses spectral radiative kernels, which is the top-of-atmosphere spectral radiance change caused by perturbations of skin temperature and air temperature and specific humidity in the atmosphere, and is an improved version of the algorithm used in earlier studies. Two improvements are the inclusion of the residual and cloud spectral kernels in the form of eigenvectors of principal components. Three and six eigenvectors are used for, respectively, the residual and cloud spectral kernels. An underlying assumption is that the spectral shape of the principal components is constant and their magnitude varies temporally and spatially. The algorithm is tested using synthetic spectral radiances with the spectral range of the Atmospheric Infrared Sounder averaged over 16 days and over a 10° × 10° grid box. Changes of skin temperature, air temperature, and specific humidity vertical profiles are derived from the difference of nadir-view all-sky spectral radiances. The root-mean-square difference of retrieved and true skin temperature differences is 0.59 K. The median of absolute errors in the air temperature change is less than 0.5 K above 925 hPa. The median of absolute errors in the relative specific humidity changes is less than 10% above 825 hPa.

Seismogenic fault and tectonic significance of 1996 Karakoram Pass earthquake (Ms 7.1) based on InSAR

Due to hard observation condition of the western Tibet region, the slip behaviors of the Ms7.1 Karakoram Pass earthquake occurred in Hetian, Xinjiang on November 19, 1996 remains unclear. Using ERS 1/2 SAR data and InSAR technique, we obtain the co-seismic deformation of the earthquake. The north and south deformation areas show asymmetric pattern, with the maximum LOS displacement of the southern part approximately 24.6 cm, and the maximum LOS displacement in the northern part approximately − 18.5 cm. Nonlinear and linear inversion algorithms are used to determine the geometric parameters and slip distribution of the earthquake fault. Our results show that the co-seismic displacement is dominated by deformation fields are clearly visible sinistral strike-slip accompanied by a small amount of normal slip component. The co-seismic slip occurred between 0 and 18 km at depth. The maximum slip is ~ 81 cm, occurring at a depth of 8.5 ± 0.5 km at (35.36°N 78.03°E), indicating a shallow event with a moment magnitude of Mw 6.5. The seismogenic fault is a secondary fault in the Karakoram fault zone with strike 96°, dip 84°, and rake – 24°. This earthquake shows that the Karakoram fault zone undergoes a complex tectonic deformation process, with central part of the fault zone showing minor tensional deformation behaviors.

Global activity shaping strategies for a retinal implant

Objective. Retinal prostheses provide visual perception via electrical stimulation of the retina using an implanted array of electrodes. The retinal activation resulting from each electrode is not point-like; instead each electrode introduces a spread of retinal activation that may overlap with activations from other electrodes. With most conventional stimulation strategies this overlap leads to image blur. Here we propose a ‘shaping’ algorithm that uses multiple electrodes to manipulate the current between electrodes in a desired way. Approach. We assume a forward model for the conversion of electrode strengths to retinal activation. Three alternative global shaping algorithms are developed by calculating reverse models under different assumptions: linear inversion using singular value decomposition to produce the pseudoinverse, a linearly constrained quadratic program, and a binary quadratic program to partition the target pattern. The algorithms were assessed using both the mean squared error between the resulting images and desired images, as well as their adherence to the maximum allowed electrode currents. Main results. Under wide activation spreads the linear inversion algorithm gave improved solutions but faced two limitations: under low-noise conditions the electrode amplitudes exceeded their set limit; the set of solutions did not include the possibility of using negative local currents to induce retinal activation. The linearly constrained quadratic program and binary quadratic program respectively addressed these problems, but required much greater computation time. Significance. This provides a framework for improving the resolution of future retinal implants, especially those with high density electrode arrays.

Untangling the nonlinearity in inverse scattering with data-driven reduced order models

The motivation of this work is an inverse problem for the acoustic wave equation, where an array of sensors probes an unknown medium with pulses and measures the scattered waves. The goal of the inversion is to determine from these measurements the structure of the scattering medium, modeled by a spatially varying acoustic impedance function. Many inversion algorithms assume that the mapping from the unknown impedance to the scattered waves is approximately linear. The linearization, known as the Born approximation, is not accurate in strongly scattering media, where the waves undergo multiple reflections before they reach the sensors in the array. Thus, the reconstructions of the impedance have numerous artifacts. The main result of the paper is a novel, linear-algebraic algorithm that uses a reduced order model (ROM) to map the data to those corresponding to the single scattering (Born) model. The ROM construction is based only on the measurements at the sensors in the array. The ROM is a proxy for the wave propagator operator, that propagates the wave in the unknown medium over the duration of the time sampling interval. The output of the algorithm can be input into any off-the-shelf inversion software that incorporates state of the art linear inversion algorithms to reconstruct the unknown acoustic impedance.

Regional structural mapping using a combined geological and geophysical approach – A preliminary study at Cairo-Suez district, Egypt

Contrary to the usual application of individual geological and/or geophysical surveys, this paper represented a sequential use of surface structural mapping and DCR method to afford a sensible structural model in a cost effective way. Understanding the geological setting of a non-mapped area at the central part of Cairo-Suez District (CSD), Egypt, is an immediate engineering concern due to its complex deformation history. Field work is carried out to measure stratigraphic sections and establish structural characteristics of the rock units. Surface geophysical DCR sounding and electrical resistivity tomography (ERT) surveys are performed. To overcome the non-uniqueness problem of DCR interpretation, the available borehole data, detailed surface geological measurements and a sequential use of linear and non-linear inversion algorithms are applied. The practical benefits of detailed surface structural mapping are demonstrating in considering all geological parameters to a conceptual subsurface model, which would have not been identified or misinterpreted otherwise. In an integrative approach, such detailed surface data allow us to confirm the uncertain inferred regional subsurface structures based on DCR data inversion. For detailed subsurface structural characterization, the results of ERTs using coarse triangular mesh and iteratively reweighted algorithm show a good agreement with prior known surface and underground structural features. Additionally, ERTs significantly resolve the complex geological features, which constitute uncertain structures in the surface structural model. The main impact of presented integrative approach is to (i) constraint the DCR data interpretation and (ii) inform a preliminary study plan of regional structural investigation, which can be used as a platform for engineering decision–making and/or further detailed geophysical and geotechnical investigations.

Evaluation of video-based linear depth inversion performance and applications using altimeters and hydrographic surveys in a wide range of environmental conditions

The performance of a linear depth inversion algorithm, cBathy, applied to coastal video imagery was assessed using observations of water depth from vessel-based hydrographic surveys and in-situ altimeters for a wide range of wave conditions (0.3 < significant wave height < 4.3 m) on a sandy Atlantic Ocean beach near Duck, North Carolina. Comparisons of video-based cBathy bathymetry with surveyed bathymetry were similar to previous studies (root mean square error (RMSE) = 0.75 m, bias = −0.26 m). However, the cross-shore locations of the surfzone sandbar in video-derived bathymetry were biased onshore 18–40 m relative to the survey when offshore wave heights exceeded 1.2 m or were greater than half of the bar crest depth, and broke over the sandbar. The onshore bias was 3–4 m when wave heights were less than 0.8 m and were not breaking over the sandbar. Comparisons of video-derived seafloor elevations with in-situ altimeter data at three locations onshore of, near, and offshore of the surfzone sandbar over ∼1 year provide the first assessment of the cBathy technique over a wide range of wave conditions. In the outer surf zone, video-derived results were consistent with long-term patterns of bathymetric change (r2 = 0.64, RMSE = 0.26 m, bias = −0.01 m), particularly when wave heights were less than 1.2 m (r2 = 0.83). However, during storms when wave heights exceeded 3 m, video-based cBathy over-estimated the depth by up to 2 m. Near the sandbar, the sign of depth errors depended on the location relative to wave breaking, with video-based depths overestimated (underestimated) offshore (onshore) of wave breaking in the surfzone. Wave speeds estimated by video-based cBathy at the initiation of wave breaking often were twice the speeds predicted by linear theory, and up to three times faster than linear theory during storms. Estimated wave speeds were half as fast as linear theory predictions at the termination of wave breaking shoreward of the sandbar. These results suggest that video-based cBathy should not be used to track the migration of the surfzone sandbar using data when waves are breaking over the bar nor to quantify morphological evolution during storms. However, these results show that during low energy conditions, cBathy estimates could be used to quantify seasonal patterns of seafloor evolution.

Characteristics in hypocenters of microseismic events due to hydraulic fracturing and natural faults: a case study in the Horn River Basin, Canada

For two to three decades, microseismic monitoring has been popular in the development of unconventional resources, because the fracture network generated by hydraulic fracturing mainly controls the productivity, and microseismic monitoring enables direct measurements for imaging the fracture network. Nevertheless, some refinements are required to make this method more practical. One challenge is to quantify the effects of pre-existing natural fractures for generating microseismic events. We determine the hypocenters of microseismic events occurring in a shale gas play in the Horn River Basin, Canada, and report several interesting spatial and temporal features of the hypocenter distributions. Automatic phase-picking is applied to waveform data recorded at 98 shallow buried three-component geophones, and phases thought to be from the same event are associated. The initial hypocenters of events are determined by iterative linear inversion algorithm then relocated using a double-difference algorithm, where relative travel time measurements are obtained with the waveform cross-correlation. We group events into many clusters based on fracking stages and their hypocenters, and then define the best-fitting plane of hypocenters for each cluster. Most strikes of the best-fitting planes are consistent with the direction of local horizontal stress maximum, indicating that hydraulic fracturing induces most microseismic events. However, the best-fitting planes of several clusters have strikes similar to those of pre-existing faults or fractures, indicating that pre-existing natural faults or fractures can affect the generation of microseismic events. In addition, some observations suggest that natural fractures can affect the temporal evolution of the spatial occurrence pattern of microseismic events. We observed specific migration patterns of microseismic events around known faults in the study area. Although further work is required for complete understanding of this phenomenon, our observations help elucidate the nature of microseismic generation.

Three-dimensional S-velocity structure of the crust in the southeast margin of the Tibetan plateau and geodynamic implications

The lower crustal flow model is one of several competing models considered to interpret the growth and expansion of southeastern Tibet. However, the dynamic processes involved in the evolution and deformation of the Tibetan plateau remain poorly understood due to the absence of reliable geophysical observations. In this study, 159 earthquakes of magnitude Ms≥6.2 were selected, which were recorded by 50 permanent broadband stations deployed in southwest China, and 1873 pairs of P-wave receiver functions with high signal-to-noise ratio were isolated. On the basis of a linearized inversion algorithm, a two-step inversion procedure was implemented that not only reduced the dependence of the inversion results on the initial models, but also provided a statistical estimation of the solution. Thus, we obtained an accurate 3D image of the S-wave velocity structure of the crust and uppermost mantle in the southeast margin of the Tibetan plateau. Our results reveal that an extensive intracrustal low-velocity zone spreads beneath southwest China, and further suggest that a lower crustal flow coming from eastern Tibet is blocked by the Jinshajiang-Red River fault to the west and the Xiaojiang fault to the east and extends largely by the Sichuan-Yunnan diamond-shaped block. This crustal flow along the southeastern margin of Tibet does not appear to be restricted to two narrow low-velocity channels as was advanced, but rather reflects the southeastward extrusion of crustal material as a way of tectonic escape.

Seismic inversion for geologic fractures and fractured media

Amplitude variation with offset (AVO) inversion attempts to use the available surface seismic data to estimate the density, P-wave velocity, and S-wave velocity of the earth model. Under linear slip interface theory, synthetic seismograms for models with fractures prove that fractures are also reflection generators. Consequently, observed reflections are not necessarily due to lithologic variations only, but they could be due in part to the effect of fractures. To obtain approximate equations for AVO inversion for fractured media, denoted by AVO with fracture (AVOF), we derived new equations for PP-wave reflection and transmission coefficients that are based on nonwelded contact boundary conditions. In particular, along with the fracture compliances, azimuth has also been taken into account in the equations because the fractures can have any orientation. The new approximate AVOF equations for a horizontally fractured medium with impedance contrast are developed by simplifying the equations for the new PP-wave reflection and transmission coefficients. In the new approximate AVOF equations, the reflection coefficients are divided into a welded contact part (a conventional impedance contrast part) and a nonwelded contact part (a fracture part). This makes the equations flexible enough to separately invert for the rock properties of the fracture and the background medium in the case of a fractured medium with impedance contrast. The new approximate AVOF equations state that fractures could cause the seismic reflectivity to be frequency dependent, and that the fractures not only influence the wave amplitude but also change the wave phase. The linear least-squares and nonlinear conjugate gradient inversion algorithms are applied to estimate the elastic reflectivity using the new approximate AVOF equations. The inverted results for seismic data for a horizontally fractured medium with impedance contrast are evaluated to find a more accurate delineation of the subsurface rock properties.

A study on the seismogenic structure of the 2016 Zaduo, Qinghai Ms6.2 earthquake using InSAR technology

On October 17th, 2016, a Ms6.2 earthquake occurred in Zaduo County of Qinghai Province, China. The aim of this study is to use synthetic aperture radar (SAR) technology aboard the Sentinel-1A satellite to obtain high-resolution co-seismic surface displacement data and then to confirm the geometric parameters of the fault and slip distribution model. To this end, linear and non-linear inversion algorithms based on an elastic half-space dislocation model were used. The results showed that a distributed slip model can explain the surface deformation field measured by InSAR very well. The surface deformation field caused by the earthquake was an oval-shaped region of subsidence with a maximum displacement of 5 cm along the line of sight of the radar waves. This earthquake was mainly the result of a normal-slip fault process with 72°N strike and 65° dip. The slip was mainly concentrated at depths of 9–15 km. The maximum slip was 0.17 m, located at a depth of 12 km. The moment magnitude given by inversion was Mw5.9. This was basically in agreement with the moment magnitudes and surface magnitudes measured by USGS and CENC.

Satellite Observations of Phytoplankton Functional Type Spatial Distributions, Phenology, Diversity, and Ecotones

Phytoplankton functional diversity plays a key role in structuring the ocean carbon cycle and can be estimated using measurements of phytoplankton functional type (PFT) groupings. Concentrations of 18 phytoplankton pigments were calculated using a linear matrix inversion algorithm, with an average r2 value of 0.70 for all pigments with p-values below the statistical threshold of 0.05. The inversion algorithm was then used with a chlorophyll-based absorption spectra model and Moderate Resolution Imaging Spectroradiometer (MODIS-Aqua) chlorophyll observations to calculate phytoplankton pigment concentrations in an area of the Atlantic Ocean off the United States east coast. Pigment distributions were analyzed to assess the distribution of PFTs. Five unique PFTs were found and delineated into three distinct offshore, transition, and open ocean groups. Group 1 (Diatoms) had highest abundance along the coast. Group 2 (prymnesiophytes, prokaryotes, and green algae) was a year-round stable offshore community that extended at reduced levels into the coast. Group 3 (dinoflagellates) dominated offshore between the Groups 1 and 2. Phytoplankton communities were delineated into coastal and offshore populations, with Group 2 having a dampened seasonal cycle, relative to the coastal populations. Shannon Diversity Indices (H) for the PFTs showed both spatial and temporal variability and had a clear non-linear relationship with chlorophyll. Diversity levels varied seasonally with changes in chlorophyll a levels. Peak PFT H was observed on the shelf where frontal features dominate, with diversity levels declining nearshore and offshore. This region marks an ecotone for phytoplankton in the study domain, and is associated with the coastal-side boundary of dinoflagellate dominance. Highest levels of diversity were observed in the tidally well-mixed regions of the Gulf of Maine and along a band that ran along the shelf region of the study area that was narrowest in the summer periods and broadened during the winter. These peak diversity zones were associated with moderate levels (~0.8 mg m-3) of chlorophyll a. While the sign in the linear trends in chlorophyll between 2002-2016 varied depending on the region, the trends in the PFT H values were nearly all negative due to the non-linear relationship between chlorophyll levels and H.

Estimate of carbonyl sulfide tropical oceanic surface fluxes using Aura Tropospheric Emission Spectrometer observations

Author(s): Kuai, L; Worden, JR; Campbell, JE; Kulawik, SS; Li, KF; Lee, M; Weidner, RJ; Montzka, SA; Moore, FL; Berry, JA; Baker, I; Denning, AS; Bian, H; Bowman, KW; Liu, J; Yung, YL | Abstract: Quantifying the carbonyl sulfide (OCS) land/ocean fluxes contributes to the understanding of both the sulfur and carbon cycles. The primary sources and sinks of OCS are very likely in a steady state because there is no significant observed trend or interannual variability in atmospheric OCS measurements. However, the magnitude and spatial distribution of the dominant ocean source are highly uncertain due to the lack of observations. In particular, estimates of the oceanic fluxes range from approximately 280 Gg S yr-1 to greater than 800 Gg S yr-1, with the larger flux needed to balance a similarly sized terrestrial sink that is inferred from NOAA continental sites. Here we estimate summer tropical oceanic fluxes of OCS in 2006 using a linear flux inversion algorithm and new OCS data acquired by the Aura Tropospheric Emissions Spectrometer (TES). Modeled OCS concentrations based on these updated fluxes are consistent with HIAPER Pole-to-Pole Observations during 4th airborne campaign and improve significantly over the a priori model concentrations. The TES tropical ocean estimate of 70±16GgS in June,when extrapolated over the whole year (about 840 ± 192 Gg S yr-1), supports the hypothesis proposed by Berry et al. (2013) that the ocean flux is in the higher range of approximately 800 Gg S yr-1.

Single and Double Sequence of Investigation Domains in 2-D Linear Inversions Applied to GPR Data

We focus on the choice of the investigation domains for a 2-D linear inversion algorithm based on the Born approximation. The application of an inversion algorithm to a wide B-scan requires the subdivision of the data into several observation domains and the partition of the buried scenario into different investigation domains, and this can be done in many ways, which raises the problem of examining the well advisedness of the adopted choice. In particular, here, we will propose the possibility of a double sequence of investigation domains. The results will also be compared with a migration algorithm.

A Linear Inverse Scattering Algorithm for Radar Imaging in Multipath Environments

This letter deals with the electromagnetic imaging in the presence of multipath propagation of interest for through-wall and urban sensing scenarios. The 2-D tomographic approach here presented combines a linear inverse scattering model, based on the Kirchhoff approximation, with the finite-difference time-domain (FDTD) technique. In particular, FDTD is exploited to evaluate the incident field and Green's function in noncanonical scenarios, so that the kernel of the linear integral equation is completely built. After, an inversion scheme based on the truncated singular value decomposition is applied to obtain a regularized solution of the problem. Numerical results demonstrate that the proposed approach yields well-focused images free of multipath ghosts, thus allowing to discriminate the actual target position. Moreover, it permits to highlight the capabilities offered by multipath exploitation such as improved crossrange resolution and detection of targets in the non-line-of-sight region of the radar.

Using Linear and Nonlinear Inversion Algorithm Combined with Simple Dislocation Model Inversion of Coal Mine Subsidence Mechanism

The long-term exploiting activity of coal mine threats the local ecological environment. Combined with surface level monitoring data of a mine in Jiangxi province, discussed the simplified dislocation model suited for coal mine, through the least squares method and Monte Carlo method to invert the Tensile component of mining area faulting. Comparing forward calculation data with level data, the accuracy of the inversion method has been analysis. The research results have realistic meanings for reducing the surface deformation mechanism due to the excavation in coal mine. It can provided theoretical basis for early warning of our country coal mining surface deformation.