Kaiser filter for antialiasing in digital photogrammetry

Image aliasing is a problem appearing as artefacts in digitally resampled images, which degrades the quality of the image. In digital rectification and texture mapping, pixels from an input image are transformed to pixels of an output image. The discrete nature of a digital image causes aliasing in the transformed image. In this paper the source of aliasing and the theory of antialiasing are described. The necessity of a precise filter design in antialiasing is discussed and a filter based on a Kaiser adjustable window is designed. Different practical antialiasing methods are described as well as interpolation methods, which are conventional in photogrammetry. Selected antialiasing methods are implemented and applied to a close range image. An objective analysis is carried out by applying inverse transformations to rectified images and deriving some measures to estimate the information loss for each method by comparing original and reconstructed images. Results indicate that interpolation methods are not capable of removing or reducing aliasing in highly decimating transformations. The output images of interpolation methods therefore suffer from edge corruption and interfusion of small features. Applying a Kaiser filter with a precise antialiasing method results in the least information loss and considerably reduces aliasing at the expense of higher computation load.     

Lithospheric structure of the Tornquist Zone resolved by nonlinear P and S teleseismic tomography along the TOR array

The main aim of the TOR project is to study the lithospheric–asthenospheric boundary structure under the Sorgenfrei–Tornquist Zone, across northern Germany, Denmark and southern Sweden. Relative arrival-time residuals of teleseismic P and S phases from 51 earthquakes, recorded by 150 seismic stations along the TOR array, were used to delineate the transition zone in the studied area. The effects of crustal structures were investigated by correcting the teleseismic residuals for travel-time variations in the crust based on a 3D crustal model derived from other data. The inversion was carried out for S phases. The results were then compared with the corresponding P-wave models. As expected, the derived models show that the relatively old and cold Baltic Shield has higher velocity at depth than the younger lithosphere farther South. The models show two sharp and distinct increases in depth to velocities which are low compared to our reference model, as we move from South to North. The location and sharpness of these boundaries suggests that the features resolved are, at least partially, compositional in origin, presumably related to mantle depletion. A sharp and steep subcrustal boundary is found roughly coincident with the southern edge of Sweden. This is below where the edge of the Baltic Shield is usually placed, based on surface geological evidence (the Sorgenfrei–Tornquist Zone). Another less significant transition is recognised more or less beneath the Elbe-lineament. Relatively high d(V_p / V_s) ratios under the central part of the profile (Denmark) indicate relatively low S-velocity in an area where a gravity high supports the hypothesis of extensive mafic intrusions.     

Transient scattering of elastic waves by dipping layers of arbitrary shape. Part 2: Plane strain model

Scattering of elastic waves by dipping layers of arbitrary shape embedded within an elastic half-space is investigated for a plane strain model by using a boundary method. Unknown scattered waves are expressed in the frequency domain in terms of wave functions which satisfy the equations of motion and appropriate radiation conditions at infinity. The steady state displacement field is evaluated throughout the elastic medium for different incident waves so that the continuity conditions along the interfaces between the layers and the traction-free conditions along the surface of the half-space are satisfied in the least-squares sense. Transient response is constructed from the steady state one through the Fourier synthesis. The results presented show that scattering of waves by dipping layers may cause locally very large amplification of surface ground motion. This amplification depends upon the type and frequency of the incident wave, impedance contrast between the layers, component of displacement which is being observed, location of the observation station and the geometry of the subsurface irregularity. These results are in agreement with recent experimental observations.     

Automatic Surface Reconstruction Based on the Fusion of Fuzzy Logic and Robust Estimation Techniques

Automatic surface reconstruction by means of digital image matching essentially involves computation in both image and object spaces. In image space, image matching is performed according to certain radiometric and geometric similarity criteria for conjugate point determination. The matching operations may therefore be regarded as decision–making processes. In practice, however, the similarity criteria for conjugate point determination demonstrate non–deterministic behaviour. Thus, the nature of these decision–making processes is more compatible with fuzzy reasoning concepts. In object space, on the other hand, surface reconstruction procedures, which operate on the basis of previously determined conjugate points, are essentially deterministic in nature and hence demand rigorous geometric modelling. In this paper, a new approach for automatic surface modelling is proposed, based on two different mathematical treatments: a fuzzy logic reasoning method for the conjugate point determination in image space; and a robust finite elements approach for the surface modelling in object space. Tests carried out on real data demonstrate the high potential of the strategy proposed.     

Groundwater Origin in Qanats,Chemo-Isotopic,and Hydrogeological Evidence

A systematic study of the chemo-isotopic characteristics and origin of the groundwater was carried out at six major qanats in the hyper-arid Gonabad area, eastern Iran. These qanats as a sustainable groundwater extraction technology have a long history, supporting human life for more than a thousand years in this region. The Gonabad qanats are characterized by outlet electrical conductivity (EC) values of 750 to 3900 µS/cm and HCO₃-Na-Mg and Cl-Na water types. The Gonabad meteoric water line (G_nMWL) was drawn at the local scale as δ²H = 6.32×δ¹⁸O + 8.35 (with R² = 0.90). It has a lower slope and intercept than the global meteoric water line due to different water vapor sources and isotope kinetic fractionation effects during precipitation in this arid region. The altitude effects on isotopic content of precipitation data were derived as δ¹⁸O = (−0.0031 × H_(m.a.s.l))−1.3). The δ²H and δ¹⁸O isotopes signatures demonstrate a meteoric origin of the groundwater of these qanats. The shift of the qanat's water samples from the local meteoric water line (LMWL) in a dry period with higher temperatures is most probably due to evaporation during the infiltration process and water movement in qanat gallery. Based on the isotopic results and mass balance calculations, the qanats are locally recharged from an area between 2000 to 2400 m.a.s.l of nearby carbonate formations and coarse alluvial sediments. The dissolution of evaporate interlayers in Neogene deposits deteriorates the groundwater quality, especially in Baidokht qanat.     

Investigation of Headcut Erosion in Cohesive Soils

Headcut, known as knickpoint migration too, is developed due to sudden change in channel bed followed by bed scour and erosion which progressing upstream. The results are the downstream morphological change and transporting massive sediment to the downstream reservoir. Most of the past studies focus on non-cohesive soils, although many problems occur because of cohesive soils. In this study, 10 different samples of cohesive soils in long term consolidation with different composition of silt and clay were tested under different circumstances of waterfall height and flow velocity to investigate the neck migration rate and the sediment yield. Tests were continued to reach a constant migration rate. One of the effective phenomena in all tests was tensional cracks on soil surface. The size and number of these cracks have inverse relation with percent of clay. Because of these cracks, massive erosion occurs at the beginning of all tests. By reducing percent of clay, headcut, waterfall height and sediment yield were increased and by reducing waterfall height and flow velocity these parameters were reduced. In lower percent of clay, headcut erosion will occur quickly with more slants. Caving phenomenon was not observed in any tests and massive erosion rate was more quickly.     

Experimental and Numerical Investigation of the Flow Hydraulic in Gradual Transition Open Channels

Channel expansions are common in both natural and artificial open channels. With increasing cross-sectional dimensions in an expansion, the flow decelerates. Due to separation of flow and subsequent eddy formation, a significant head loss is occurred along the transition. This study presents the results of experimental investigations on subcritical flow along the expansive transition of rectangular to trapezoidal channels. Also, a numerical simulation was developed using the finite volume method with Reynolds Stress turbulent model. Water surface profiles and velocity distributions of flow through the transition were measured experimentally and compared with the numerical results. Also, hydraulic efficiency of the transition and coefficient of energy head loss were calculated. The results show that with increasing the upstream Froude number, hydraulic efficiency of the transition and coefficient of energy head loss are decreased and increased, respectively. The results also showed the ability of numerical simulation for simulating the flow separation zones and bed shear stress along the transition for different inlet discharges and inflow Froude numbers.     

Effects of non-submerged boulder on flow characteristics – A field investigation

The effect of fully submerged boulders on the flow structure in channels has been studied by some researchers. However, many natural streams have bed material with boulders that are not fully submerged under water. In many natural streams, boulders cover between 1% and 10% of the area of the stream reach. The effect of non-submerged boulders on the velocity profile and flow characteristics is very important for assessing riverbed deformation. The objectives of this paper are to find the pattern of velocity distribution around a non-submerged boulder and to compare it with the classical studies on flow resistance and Reynolds stress distribution in open channels. Also, by considering the variation in the Reynolds stress distribution at different locations around a non-submerged boulder, the effect of a non-submerged boulder on the estimation of shear velocity and resistance to flow has been investigated. Results indicates that inside the scour hole caused by a non-submerged boulder in a river velocity distributions are irregular. However, velocity distributions are regular outside the scour hole. The presence of the boulder causes a considerable deviation of the Reynolds shear stress from the classic distribution, showing a non-specific distribution with negative values. The classical methods for calculating shear velocity are not suitable because these methods do not give detailed velocity and Reynolds stress distributions in natural rivers with a lot of boulders. Thus, the effect of a non-submerged boulder on the estimation of the resistance to flow by considering the variations in velocity and Reynolds stress distributions at different locations around a non-submerged boulder is important and needs to be studied in a natural river instead of just in laboratory flumes. The negative values in Reynolds stress distribution around a boulder indicate that the classical methods are unable to predict resistance to flow, and also show strong turbulence inside the scour hole where the complex flow conditions present ambiguous Reynolds stress distributions. In the current study, to obtain a reasonable estimation of parameters in natural rivers, the classical method has been modified by considering velocity and Reynolds stress distributions through the boundary layer method.