The lecture notes describe the Delaunay Tessellation Field Estimator for Cosmic Web analysis. The high sensitivity of Voronoi/Delaunay tessellations to the local point distribution is used to obtain estimates of density and related quantities. The adaptive and minimum triangulation properties of Delaunay tessellations are utilized for interpolation to continuous volume-covering density and velocity fields. The adaptivity to the local density and geometry preserves three key characteristics of the cosmic matter distribution: 1) its hierarchical nature,2) its anisotropic weblike morphology, 3) the presence and structure of voids. The lecture notes discuss the mathematical context and background of the DTFE method, in essence a first order version of Natural Neighbour methods, smooth and higher order spatial interpolation techniques. This is followed by an outline of the DTFE procedure. A series of tests on spatially complex point distributions provides quantitative information on its performance and noise characteristics. We conclude with a few cosmological applications. The reconstruction of the spatial weblike patterns in N-body simulations and in the 2dFGRS provide a visual impression of DTFE's performance. The analysis of the velocity divergence pdf and maps obtained of the velocity flow in the Local Universe (implied by PSCz) illustrate its promise for probing the dynamics of cosmic structure. The Watershed Void Finder (Platen et al. 2007), for the study of voids, and Multiscale Morphology Filter (Aragon-Calvo et al. 2007), for identifying and analyzing weblike networks of filaments, sheets and clusters, are based upon DTFE density fields. We also mention SimpleX, a related random lattice radiative transfer formalism (Ritzerveld & Icke 2006).

We present the results of a study confronting density maps reconstructed by the Delaunay Tessellation Field Estimator (DTFE) and by regular SPH kernel-based techniques. The comparison between the two methods clearly demonstrates the superior performance of the DTFE with respect to conventional SPH methods, in particular at locations where SPH appears to fail. The DTFE is a fully self-adaptive technique for reconstructing continuous density fields from discrete particle distributions, and is based upon the corresponding Delaunay tessellation. Its principal asset is its complete independence of arbitrary smoothing functions and parameters specifying the properties of these. As a result it manages to faithfully reproduce the anisotropies of the local particle distribution and through its adaptive and local nature proves to be optimally suited for uncovering the full structural richness in the density distribution. Through the improvement in local density estimates, calculations invoking the DTFE will yield a much better representation of physical processes which depend on density. The presented results form an encouraging step towards the application and insertion of the DTFE in astrophysical hydrocodes. We describe an outline for the construction of a particle hydrodynamics code in which the DTFE replaces kernel-based methods. Further discussion addresses the issue and possibilities for a moving grid based hydrocode invoking the DTFE, and Delaunay tessellations, in an attempt to combine the virtues of the Eulerian and Lagrangian approaches.

We use a simple analytic model to deproject 2-d luminosity functions (LF) of galaxies in the Coma cluster measured by Beijersbergen et al. 2002. We demonstrate that the shapes of the LFs change after deprojection. It is therefore essential to correct LFs for projection effects. The deprojected LFs of the central area have best-fitting Schechter parameters of M^{*}_U=-18.31^{+0.08}_{-0.08} and \alpha_U=-1.27^{+0.018}_{-0.018}, M^{*}_B=-19.79^{+0.14}_{-0.15} and \alpha_B=-1.44^{+0.016}_{-0.016} and M^{*}_r=-21.77^{+0.20}_{-0.28} and \alpha_r=-1.27^{+0.012}_{-0.012}. The corrections are not significant enough to change the previously observed trend of increasing faint end slopes with increasing distance to the cluster center. The weighted U, B, and r band slopes of the deprojected LFs show a slightly weaker steepening with increasing projected cluster radius.

Optimally resolved one-dimensional density and velocity profiles through cosmological N-body simulations are constructed by means of the Voronoi-Delaunay tessellation reconstruction technique. In a fully self-adaptive fashion a strikingly detailed view of the density features and the corresponding cosmic motions is recovered.

Astronomical observations, physical experiments as well as computer simulations often involve discrete data sets supposed to represent a fair sample of an underlying smooth and continuous field. Reconstructing the underlying fields from a set of irregularly sampled data is therefore a recurring key issue in operations on astronomical data sets. Conventional methods involve artificial filtering through a grid or a smoothing kernel and fail to achieve an optimal result. Here we describe a fully self-adaptive geometric method which does not make use of artificial filtering, and which makes optimal use of the available information.

The application of Voronoi and Delaunay tessellation based methods for reconstructing continuous fields from discretely sampled data sets is discussed. The succesfull operation as ``multidimensional interpolation'' method is corroborated through their ability to reproduce even intricate statistical aspects of the analytical predictions of perturbation theory of cosmic velocity field evolution. The newly developed and fully self-adaptive technique for density field estimation by means of Delaunay tessellations, basically exploiting their ``minimum triangulation'' quality, is shown to succesfully reproduce the morphology of the foamlike structure in N-body simulations of cosmic structure formation. The full hierarchy of structure is implicitly and directly reproduced at every spatial resolution scale present in the particle distribution, at the same time automatically and sharply rendering its characteristic anisotropic filamentary and wall-like features.

Here we introduce the Delaunay Density Estimator Method. Its purpose is rendering a fully volume-covering reconstruction of a density field from a set of discrete data points sampling this field. Reconstructing density or intensity fields from a set of irregularly sampled data is a recurring key issue in operations on astronomical data sets, both in an observational context as well as in the context of numerical simulations. Our technique is based upon the stochastic geometric concept of the Delaunay tessellation generated by the point set. We shortly describe the method, and illustrate its virtues by means of an application to an N-body simulation of cosmic structure formation. The presented technique is a fully adaptive method: automatically it probes high density regions at maximum possible resolution, while low density regions are recovered as moderately varying regions devoid of the often irritating shot-noise effects. Of equal importance is its capability to sharply and undilutedly recover anisotropic density features like filaments and walls. The prominence of such features at a range of resolution levels within a hierarchical clustering scenario as the example of the standard CDM scenario is shown to be impressively recovered by our scheme.

The neutral hydrogen line profiles along the major axis of the nearby spiral galaxy NGC 2403 show a wing towards the sytemic velocity. This asymmetry can be explained with the presence of an abnormally thick HI disk (FWHM ~ 5 kpc) or with a two-component structure: a thin disk and a slowly rotating, thicker (1-3 kpc) HI layer. The latter model gives a better representation of the observations. These results throw a new light on the disk-halo connection. In particular, the decrease of rotational velocity with height above the plane may be the result of a galactic fountain flow. A vertically extended, slowly rotating HI layer may be common among spiral galaxies with high levels of star formation.