3D Profile Research Articles

BRAS: Bidirectional Reflectance Adjustment Strategy for 3-D Reconstruction of Mirror-Like Surface

A mirror-like surface (MLS) reflects highlight, aggravating the image saturation in structured light three-dimensional (3D) reconstruction systems and precluding defect detection based on reconstructed 3D profiles. Previous studies have focused on a strategy of limiting the luminous flux entering the camera. However, the high intensity of the reflected highlight force the limitation to be strengthened, which heavily reduces the modulation in the images for the 3D reconstruction. Therefore, we propose a new strategy to adjust the source of the reflected highlight, i.e., the bidirectional reflectance (BR) of MLS, which fundamentally suppresses the highlight and removes the limitation on the luminous flux. To execute the proposed strategy, an unfixed view structured light system (UVSLS) is established. The UVSLS converts the reflection viewer from the real camera to the virtual camera, realizing the flexible adjustment of the MLS BR. Finally, an MLS 3D reconstruction framework is constructed to obtain the 3D profile of MLS. Experiments demonstrate that the proposed framework reduces the saturated pixels by 73.72 <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\%$</tex-math></inline-formula> compared with the conventional method. Compared with the previous methods, the saturated pixels are reduced by an average of 47.62 <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\%$</tex-math></inline-formula> .

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Analysis of the Effect of Tilling and Crop Type on Soil Structure Using 3D Laser Profilometry

Background and aim: Soil structure is an important indicator of the quality of soil, but detecting the early signs of soil degradation from soil structure is difficult. Developing precise instruments able to diagnose soil structure quickly is therefore critical to improve management practices. Here, the objective is to develop an instrument analysing the roughness of surfaces resulting from the fracture of soil cores, and to test the instrument’s ability to detect changes in soil structure cause by crop type and tillage. We have designed and constructed a 3D laser profilometer suitable for analysing standard soil cores. The 3D soil profiles were first assembled into a 3D surface using image analysis before roughness indicators could be computed. The method was tested by analysing how soil surface roughness was affected by crop varieties (barley and bean) and tillage (conventional tillage and no tillage). Results showed the method is precise and could reliably detect an influence of crop type and tillage on the roughness indicator. It was also observed that tillage reduced the difference in soil structure between the different cultures. Also, the soil in which barley grew had significantly lower roughness, irrespective of the tillage method. This could indicate that the roughness indicator is affected by biopores created by the root system. In conclusion, roughness indicators obtained from the fracture of soil cores can be easily obtained by laser profilometry and could offer a reliable method for assessing the effect of crop types and soil management on soil quality.

Cement sheath integrity at simulated reservoir conditions of pressure, temperature, and wellbore configuration in a laboratory setup

Cement is fundamental to maintaining integrity throughout the life of a well by preventing leakage of fluids. Cement is thus crucial in oil and natural gas producing wells as well as in wells used for storage and utilization of carbon dioxide and hydrogen. Even if the cement slurry has been properly placed, changes in downhole reservoir conditions can induce thermal and mechanical stresses to harm the integrity of the cement sheath.This work presents a novel apparatus to analyze the cement sheath performance at simulated reservoir conditions of pressure, temperature, and wellbore configuration. One of the major features demonstrated is three-dimensional (3D) monitoring of in-situ dynamic changes in cement using acoustic data. 3D density profiles allow visual observance of the evolution of possible settling over the cement curing process. A first arrival pick time decrease of 26.6% is observed over a 15-h curing time. Confinement improves the strength of cement as seen by decrease in travel time (10.8%) and increase of amplitude (32%) with 15 MPa confining pressure increase. Temperature increase of 66.4 °C has damaging effects on the strength of cement as seen by increase (11.9%) in travel time and decrease (20.5%) in amplitude. An increase in cement shear bond strength from 13% to 89.7% with longer curing time is measured using this device.

Exact solutions and bifurcations of the time-fractional coupled Boussinesq-Burgers equation

In this paper, the time-fractional coupled Boussinesq-Burgers equations are studied to construct exact solutions by utilizing the modified extended tanh-function method, the improved method, the two variables -expansion method and bifurcation analysis. Firstly, based on the modified Riemann-Liouville derivative and traveling wave transformation, the equations are rewritten as an ordinary differential equation. Some exact solutions, including singular, periodic, kink, anti-kink, soliton and periodic-singular solutions, are derived. Then, we graphically describe the 3D, density plots and solution profiles of some solutions by setting suitable parameters. Finally, the effectiveness of the methods applied is verified via bifurcation theory.

Costal Cartilage Graft Repair Osteochondral Defect in a Mouse Model.

Osteochondral defects develop into osteoarthritis without intervention. Costal cartilage can be utilized as an alternative source for repairing osteochondral defect. Our previous clinical study has shown the successful osteochondral repair by costal cartilage graft with integration into host bone bed. In this study, we investigate how cartilaginous graft adapt to osteochondral environment and the mechanism of bone-cartilage interface formation. Costal cartilage grafting was performed in C57BL/6J mice and full-thickness osteochondral defect was made as control. 3D optical profiles and micro-CT were applied to evaluate the reconstruction of articular cartilage surface and subchondral bone as well as gait analysis to evaluate articular function. Histological staining was performed at 2, 4, and 8 weeks after surgery. Moreover, costal cartilage from transgenic mice with fluorescent markers were transplanted into wild-type mice to observe the in vivo changes of costal chondrocytes. At 8 weeks after surgery, 3D optical profiles and micro-CT showed that in the graft group, the articular surface and subchondral bone were well preserved. Gait analysis and International Cartilage Repair Society (ICRS) score evaluation showed a good recovery of joint function and histological repair in the graft group. Safranin O staining showed the gradual integration of graft and host tissue. Costal cartilage from transgenic mice with fluorescent markers showed that donor-derived costal chondrocytes turned into osteocytes in the subchondral area of host femur. Costal cartilage grafting shows both functional and histological repair of osteochondral defect in mice. Graft-derived costal chondrocytes differentiate into osteocytes and contribute to endochondral ossification.

Estimation of 2D profile dynamics of electrostatic potential fluctuations using multi-scale deep learning

Abstract Dynamics in magnetically confined plasmas are dominated by turbulence driven by spatial inhomogeneities in density and temperature. Simultaneous measurement of velocity field and density fluctuations is necessary to observe the particle transport, but the measurement of the velocity field fluctuations is often challenging. Here, we propose a method to estimation velocity field fluctuations from density fluctuations by using plasma turbulence simulations and a deep technique learning. In order to take multi-scale characteristics into account, the several number of spatial filters are used in the convolutional neural network. The velocity field fluctuations are successfully predicted, and the particle transport estimated from the predicted velocity field fluctuations is within 93.1% accuracy. The deep learning could be used for the prediction of physical variables which are difficult to be measured.

Wear Behaviour of Additive Manufactured Aluminium Alloy ER 5356

In the automotive industry, parts are mostly made from aluminium alloy due to its lightweight properties and high corrosion resistance. However, the drawback is that the aluminium alloy is easily worn due to wear and friction and will end up in the scrap yard. In order to salvage the aluminium component, the worn part can be repaired. Currently, wire arc additive manufacturing (WAAM) offers flexible remanufacturing of the worn part. However, the wear behaviour of the additively manufactured part needs to be studied first to improve the wear performance of the material. In this study, the gas metal arc welding (GMAW) or MIG-based WAAM machine was utilised to produce a 3D profile from the available aluminium alloy wire grade ER 5356. The wear test was carried out in accordance with ASTM G-99, using a pin-on disc in both dry and wet sliding conditions. It was found that on dry sliding, the specific wear rates are decreasing from 5.3632 x 10-11 mm3/Nm to 4.3496 x 10-11 mm3/Nm and 4.1513 x 10-11 mm3/Nm as the speed increases from 200 to 400 RPM at the constant 20 N load. Meanwhile, for wet sliding, it has been observed that the specific wear rate increases as similar speed values are used in dry sliding conditions, which are 6.8122 x 10-12 mm3/Nm, 1.1931 x 10-11 mm3/Nm and 3.7561 x 10-11 mm3/Nm with a similar constant 20 N load. Next, the coefficient of friction for dry sliding shows that as the speed decreases. In contrast, for wet sliding, it is observed that the coefficient of friction increases.

Machine learning-directed electrical impedance tomography to predict metabolically vulnerable plaques.

The characterization of atherosclerotic plaques to predict their vulnerability to rupture remains a diagnostic challenge. Despite existing imaging modalities, none have proven their abilities to identify metabolically active oxidized low-density lipoprotein (oxLDL), a marker of plaque vulnerability. To this end, we developed a machine learning-directed electrochemical impedance spectroscopy (EIS) platform to analyze oxLDL-rich plaques, with immunohistology serving as the ground truth. We fabricated the EIS sensor by affixing a six-point microelectrode configuration onto a silicone balloon catheter and electroplating the surface with platinum black (PtB) to improve the charge transfer efficiency at the electrochemical interface. To demonstrate clinical translation, we deployed the EIS sensor to the coronary arteries of an explanted human heart from a patient undergoing heart transplant and interrogated the atherosclerotic lesions to reconstruct the 3D EIS profiles of oxLDL-rich atherosclerotic plaques in both right coronary and left descending coronary arteries. To establish effective generalization of our methods, we repeated the reconstruction and training process on the common carotid arteries of an unembalmed human cadaver specimen. Our findings indicated that our DenseNet model achieves the most reliable predictions for metabolically vulnerable plaque, yielding an accuracy of 92.59% after 100 epochs of training.

Design Optimization of a Cavitating 2D Profile in Proximity of the Free Surface

The appearance of AC72 foiling catamarans in the scenario of sailing yacht competitions in 2013 raised attention to this ship design concept, although not brand new in the yacht design history. The drastic drag reduction connected with the elevation of the ship hull outside the water is obtained by the use of a foil, or a system of foils, acting as the wings of a plane, providing a lift force balancing the weight of the ship. Since this lift is proportional (non-linearly) to the ship hull speed, the take-off speed of the hull cannot be low. As a result, since we are travelling in water at high speeds, the occurrence of the phenomenon of cavitation cannot be completely avoided, and the performance of the ship undergoes deterioration. Shaping of the foil profile must consider this peculiar situation, so the design tools commonly adopted for the aero-hydrodynamic hull design optimization are no longer adequate. In this paper, we are considering the optimization of the 2D profile of a foil in three different physical conditions: single fluid, two fluids and two fluids with cavitation. The first is typical of aeronautic wing design, the second of the appendages of a displacement ship, and the third of a foiling ship. Results give evidence of the different requirements for the three different conditions.

Femtosecond laser micro-machining of three-dimensional surface profiles on flat single crystal sapphire

In this study, a regression model is proposed using the machine learning (ML) method for femtosecond laser micro-machining of three-dimensional (3D) surface profiles on flat single-crystal sapphire. Based on the ML regression model, a quantitative relationship is established between the femtosecond laser processing parameters and the processed 3D functional curved surface topography. The combination of laser processing parameters (i.e. laser fluence F0, hatch dh, and scanning speed v) is also optimized according to the predicted parameters (i.e. surface quality Ra, ablation depth Z, and material removal rate MRR). Finally, the optimized combination of laser processing parameters was employed to process 3D surface profiles on flat single-crystal sapphire sheets. The maximum predicted profile error of the 3D structures was approximately 12.6% without any prior optimization experiments. Both the numerical simulation and experimental results demonstrate the feasibility of the proposed ML regression model in optimizing laser processing parameter combinations to achieve the desired 3D surface geometries on single-crystal sapphire by femtosecond laser micro-machining without trial-and-error.

Optical soliton for (2+1)-dimensional coupled integrable NLSE using Sardar-subequation method

The [Formula: see text]-dimensional coupled nonlinear Schrödinger equation has versatile applications for modeling nonlinear waves in different areas such as optics, atmospheric science, fluid dynamics, and plasma physics. This study focuses on the propagation of optical solitons and their interaction within various mediums such as multi-mode fiber and fiber arrays. The Sardar-subequation method is applied to achieve the bright, dark, combined bright–dark, periodic, combined dark-singular and singular optical solitons solutions. The 3D, contour and 2D profiles for some of the assimilation solutions are also plotted to show the physical behavior of the solutions. The obtained results demonstrate the effectiveness of the adopted methodology in obtaining traveling wave solutions and have many applications in applied mathematics and engineering. The novelty of our work lies in the successful application of Sardar-subequation method in obtaining diverse soliton solutions and exploring their physical behavior. This study goes beyond previous efforts in the literature by presenting a comprehensive analysis of soliton dynamics. This ability to describe and predict the behavior of optical solitons in diverse mediums opens up new opportunities for investigating the existing systems.

Subsurface mapping of petroleum product (diesel) spillage in Apapa, Lagos, Nigeria from geophysical, lithological, and hydrogeological studies

A multi-disciplinary method has been applied for the subsurface mapping of a diesel (refined petroleum product) spillage in the Apapa area of Lagos, Nigeria. The methods include geological, geophysical, hydrogeological, and geoinformatics approaches. Apapa is located within the eastern sector of the Dahomey Basin, Nigeria. Twelve profiles were probed, adopting a 2D Wenner configuration. Three 3D resistivity images were produced from the 2D profiles using ZONDRES 3D and SURFER 16. Additionally, six boreholes were drilled for lithological information and the determination of groundwater flow directions of the area. The study area is characterized by loose, fine to coarse grained sand underlain by peat. The sands were observed to be partly contaminated by diesel (spilled diesel). The 3D resistivity images revealed that the area is composed of 3 geoelectric layers; an organic debris layer (<15 Ωm), overlain by a clean sandy layer (15–54 Ωm), as well as a contaminated sand layer (>55 Ωm). The results showed that higher resistivity values were obtained at shallow depths ≤3 m which depicts the diesel contaminated sand. The 3D model also showed that the magnitude of contamination decreases from the southern part (close to the source of contamination) towards the north of the area. These results are in tandem with the principal flow direction of water which is from south to north, confirming that the pollution source is from the southern part of the area. The combination of the resistivity, groundwater flow direction, and lithological information, therefore, unraveled the source, direction, and extent of contamination in the study area for a successful remediation plan.

Tunable ultrafast thermionic emission from femtosecond-laser hot spot on a metal surface by control of laser polarization and angle of incidence: A numerical investigation

Ultrafast laser induced thermionic emission from metal surfaces has several applications. Here, we investigate the role of laser polarization and angle of incidence on the ultrafast thermionic emission process from laser driven gold coated glass surface. The spatio-temporal evolution of electron and lattice temperatures are obtained using an improved three-dimensional (3D) two-temperature model (TTM) which takes into account the 3D laser pulse profile focused obliquely onto the surface. The associated thermionic emission features are described through the modified Richardson-Dushman equation, including dynamic space-charge effects and are included self-consistently in our numerical approach. We show that temperature-dependent reflectivity influences laser energy absorption. The resulting peak electron temperature on the metal surface monotonically increases with the angle of incidence for the P polarization, while for the S polarization it shows the opposite trend. We observe that thermionic emission duration shows a strong dependence on the angle of incidence and contrasting polarization dependent behavior. The duration of the thermionic current shows strong correlation to the intrinsic electron-lattice thermalization time, in a fluence regime well below the damage threshold of gold. The observations and insights have important consequences in designing ultrafast thermionic emitters using a metal based architecture.

Initializing and accelerating Stereo-DIC computation using semi-global matching with geometric constraints

Stereo Digital Image Correlation (Stereo-DIC) has become a mainstream optical metrology technique for quantitatively analyzing full-field 3D shape, displacement, or deformation of materials and structures. Whether it is to measure 3D profile or deformation, stereo matching is essential for Stereo-DIC to reconstruct 3D point clouds from stereo images. Traditional feature-based (e.g., SIFT) methods provide initial 2D displacements for stereo matching with the aid of extracting massive features, but at the cost of expensive computational overhead. In addition, these methods preclude precise measurement of objects with steep and ridged surfaces or undergoing large rotation and/or deformation due to low feature matching accuracy of complex regions caused by perspective differences. In this paper, we propose a fast and robust stereo matching method using semi-global matching with geometric constraints (GC-SGM) for initializing and accelerating Stereo-DIC computation. For GC-SGM, an optimized semi-global matching (SGM) algorithm based on GPU acceleration is first utilized to quickly estimate dense and reliable disparity maps between the rectified stereo images. The global pixel-wise 2D correspondence between raw stereo images can be established inversely using epipolar constraints and 1D disparity information, and then converted to accurate and initial second-order deformation parameters for 2D-DIC-based sub-pixel refinement by least-squares-based surface fitting. Experimental results prove that the proposed GC-SGM enhances the matching correctness and robustness for complex objects while improving the processing speed on GPU by 3∼10 times compared with SIFT-based methods, enabling high-precision and computationally efficient 3D shape and deformation measurement.

Enhancing the mechanical stability of restored teeth with interfacial cracks: Finite element analysis

ObjectivesThis study aims to enhance the mechanical stability of restored molar teeth with class II occlusal-distal (OD) cavities. We seek to achieve this goal through a comprehensive investigation of three primary factors: (1) the choice of restoration material properties, (2) internal cavity geometries, and (3) the impact of double-layered restoration configurations. MethodsTo achieve our objectives, we initiated by creating two-dimensional (2D) models of restored teeth featuring class II OD cavities, utilizing scanned and segmented images of maxillary molar teeth. We drew 2D profiles of dentine and enamel, which were then imported into finite element analysis (FEA) software. To explore various cavity geometries, we implemented a total of thirteen different designs, encompassing straight, oblique, grooved, curved, and double-layered configurations. We utilized a semi-circular stone to simulate the application of contact load on the restored tooth. We applied identical boundary conditions and contact loading across all models. To assign material properties, we developed a Python code, enabling the automatic assignment of seven elastic moduli ranging from 2 GPa to 26 GPa to the restoration materials. Meanwhile, constant material properties were assigned to the enamel and dentine. In total, we conducted 133 FEA simulations to comprehensively analyse the effects of the aforementioned factors on the strength and performance of restored molar teeth. ResultsOur analysis revealed two key factors significantly influencing the mechanical resistance of treated teeth, particularly in the presence of a crack or debonding: (1) the marginal geometry of the OD cavity and (2) the elastic modulus of the restorative material. However, altering the internal cavity angle and implementing a double-layered restoration did not significantly influence the restored tooth's overall strength and performance in the face of crack or debonding situations. SignificanceThe findings of this study have substantial implications for designing and restoring class II OD cavities to enhance resistance to cracks or debonding. The use of curved marginal geometries in restoration design can significantly improve fracture resistance, with double-curved geometries reducing stress concentrations by approximately 43% compared to straight cavities. These results offer valuable guidance for strengthening the structural integrity of restored teeth, calling for further experimental investigations to explore practical applications and benefits.

Laser Profilometry Application in Welding Geometrical Characterization

Welding processes are commonly used in the industry for the manufacture of large parts or due to their complex geometry that does not allow the part to be manufactured as a whole. Nevertheless, the union can show worst mechanical properties than the rest of the piece, affecting negatively its service behavior, so it is necessary evaluate weald seam to ensure the correct process application. Electrical welding operations are commonly used due to the reduced equipment size or their possibilities application in numerous metallic materials. Notwithstanding, different variables have to be taken into account during the metal deposition, as intensity or speed deposition, among others. Weald seam geometrical evaluation is usually utilized to validate the union surface conditions. Despite this, surface irregularities caused during the process make its difficult to measure correctly with conventional equipment. However, laser profilometry is a non-contact technique that can be used to generate 3D profile of weald seam, facilitating its measurement with high accuracy. Therefore, in this work an initial analysis of the influence of material deposition speed and arc welding intensity on the weald seam geometry will be carried out using a laser profilometry equipment. In addition, to ensure a correct information acquisition, the laser profilometer requires a constant speed movement in relation with the weald seam analyzed, so new equipment has been manufactured, using additive manufacturing techniques, to support the profilometer throughout the information acquisition process.

Pseudo-3D cubes from densely spaced subbottom profiles via projection onto convex sets interpolation: An open-source workflow applied to a pockmark field

High-resolution subbottom profiler (SBP) data are commonly recorded during academic or industrial surveys and analyzed as 2D profiles. However, there is growing interest in 3D imaging of the shallow seafloor subsurface, e.g., site surveys prior to the installation of offshore wind farms. Despite ongoing advancements in 3D systems, SBPs continue to be used mainly for 2D profiling. In other seismic imaging applications, “pseudo-3D” seismic cubes have been successfully generated from dense 2D seismic surveys. We developed a novel open-source Python-based workflow to interpolate densely spaced 2D SBP profiles acquired with a hull-mounted parametric SBP into pseudo-3D cubes. This workflow is applied to a study area on the Chatham Rise east of New Zealand’s South Island, comprising numerous seafloor pockmarks and subsurface paleo-pockmarks. Our workflow consists of two stages: (1) processing of individual 2D profiles to compensate for existing vertical offsets and (2) sparse data interpolation by applying the iterative projection onto convex sets algorithm. We correct the residual static effect caused by the sea state during acquisition, compensate for the varying tidal elevations, and apply vertical time shifts to individual profiles to correct misties at profile intersections. The interpolation of the binned sparse 3D cube is performed in the frequency domain using the fast Fourier transform and executed in parallel processes, which significantly increases the computational efficiency. We validate our workflow by testing different bin sizes, sparse transforms, thresholding approaches, and varying total iterations, which noticeably impact the resulting interpolation quality. Our developed workflow generates pseudo-3D subsurface images from densely spaced SBP profiles, with numerous potential applications for academic and industrial surveys.

The effects of uplift and erosion on the petroleum systems in the southwestern Barents Sea: Insights from seismic data and 2D petroleum systems modelling

In the past three decades, the hydrocarbon prospectivity has been enigmatic in the southwestern Barents Sea as gas discoveries have dominated over oil. However, more recent discoveries indicate oil potential, especially in the Polheim Platform and Hoop Complex areas. Cenozoic uplift and erosion episodes are critical because they control the duration and depth of maximum burial of source rocks and the reduction of porosity of reservoir rocks due to diagenetic processes at larger depth and increased temperatures. Modelling of interpreted horizons from the southwestern Barents Sea provide a best-fit realization of the basin-scale sedimentary filling from the post-rifting in the Jurassic until the Last Glacial Maximum. 2D basin modelling was performed across one regional seismic section crossing the southwestern Barents Sea. For calibration of the 2D model, 1D modelling was performed for a selected number of wells in the region by integrating different estimates of net apparent erosion (velocity inversion of wells, interpreted profiles and structure maps as well as from temperature and vitrinite reflectance data) in order to assess the burial, thermal and maturity history of the source rocks at various well locations as well as along the 2D profile. Simulation results show that the Upper Jurassic source rock is immature to marginal mature in the central and eastern parts of the Norwegian Barents Sea. The western part is characterised by high chargeability of traps associated with rich source kitchens and porous reservoir rocks. Although not proven by commercial discoveries, there are evidences that the Permian Ørret and the Lower-Middle Triassic basal Klappmyss and basal Kobbe shales may be important gas-prone source rocks in the eastern Norwegian Barents Sea. Consequently, the exploration must rely on the Triassic and/or Palaeozoic source rocks in the eastern area.

3D wind field profiles from hyperspectral sounders: revisiting optic-flow from a meteorological perspective

In this work, we present an efficient optic flow algorithm for the extraction of vertically resolved 3D atmospheric motion vector (AMV) fields from incomplete hyperspectral image data measures by infrared sounders. The model at the heart of the energy to be minimized is consistent with atmospheric dynamics, incorporating ingredients of thermodynamics, hydrostatic equilibrium and statistical turbulence. Modern optimization techniques are deployed to design a low-complexity solver for the energy minimization problem, which is non-convex, non-differentiable, high-dimensional and subject to physical constraints. In particular, taking advantage of the alternate direction of multipliers methods (ADMM), we show how to split the original high-dimensional problem into a recursion involving a set of standard and tractable optic-flow sub-problems. By comparing with the ground truth provided by the operational numerical simulation of the European Centre for Medium-Range Weather Forecasts (ECMWF), we show that the performance of the proposed method is superior to state-of-the-art optical flow algorithms in the context of real infrared atmospheric sounding interferometer (IASI) observations.