An Advanced, Three-Dimensional Plotting Library for Astronomy

We present a new, three-dimensional (3D) plotting library with advanced features, and support for standard and enhanced display devices. The library - S2PLOT - is written in C and can be used by C, C++ and FORTRAN programs on GNU/Linux and Apple/OSX systems. S2PLOT draws objects in a 3D (x,y,z) Cartesian space and the user interactively controls how this space is rendered at run time. With a PGPLOT inspired interface, S2PLOT provides astronomers with elegant techniques for displaying and exploring 3D data sets directly from their program code, and the potential to use stereoscopic and dome display devices. The S2PLOT architecture supports dynamic geometry and can be used to plot time-evolving data sets, such as might be produced by simulation codes. In this paper, we introduce S2PLOT to the astronomical community, describe its potential applications, and present some example uses of the library.Comment: 12 pages, 10 eps figures (higher resolution versions available from http://astronomy.swin.edu.au/s2plot/paperfigures). The S2PLOT library is available for download from http://astronomy.swin.edu.au/s2plo

Halo Coronal Mass Ejections during Solar Cycle 24: reconstruction of the global scenario and geoeffectiveness

In this study we present a statistical analysis of 53 fast Earth-directed halo CMEs observed by the SOHO/LASCO instrument during the period Jan. 2009-Sep. 2015, and we use this CME sample to test the capabilities of a Sun-to-Earth prediction scheme for CME geoeffectiveness. First, we investigate the CME association with other solar activity features by means of multi-instrument observations of the solar magnetic and plasma properties. Second, using coronagraphic images to derive the CME kinematical properties at 0.1 AU, we propagate the events to 1 AU by means of the WSA-ENLIL+Cone model. Simulation results at Earth are compared with in-situ observations at L1. By applying the pressure balance condition at the magnetopause and a solar wind-Kp index coupling function, we estimate the expected magnetospheric compression and geomagnetic activity level, and compare them with global data records. The analysis indicates that 82% of the CMEs arrived at Earth in the next 4 days. Almost the totality of them compressed the magnetopause below geosynchronous orbits and triggered a geomagnetic storm. Complex sunspot-rich active regions associated with energetic flares result the most favourable configurations from which geoeffective CMEs originate. The analysis of related SEP events shows that 74% of the CMEs associated with major SEPs were geoeffective. Moreover, the SEP production is enhanced in the case of fast and interacting CMEs. In this work we present a first attempt at applying a Sun-to-Earth geoeffectiveness prediction scheme - based on 3D simulations and solar wind-geomagnetic activity coupling functions - to a statistical set of potentially geoeffective halo CMEs. The results of the prediction scheme are in good agreement with geomagnetic activity data records, although further studies performing a fine-tuning of such scheme are needed.Comment: Accepted for publication in the Journal of Space Weather and Space Climate (SWSC

Contributions to virtual reality

153 p.The thesis contributes in three Virtual Reality areas: ¿ Visual perception: a calibration algorithm is proposed to estimate stereo projection parameters in head-mounted displays, so that correct shapes and distances can be perceived, and calibration and control procedures are proposed to obtain desired accommodation stimuli at different virtual distances.¿ Immersive scenarios: the thesis analyzes several use cases demanding varying degrees of immersion and special, innovative visualization solutions are proposed to fulfil their requirements. Contributions focus on machinery simulators, weather radar volumetric visualization and manual arc welding simulation.¿ Ubiquitous visualization: contributions are presented to scenarios where users access interactive 3D applications remotely. The thesis follows the evolution of Web3D standards and technologies to propose original visualization solutions for volume rendering of weather radar data, e-learning on energy efficiency, virtual e-commerce and visual product configurators

Cosmic Rays: The Second Knee and Beyond

We conduct a review of experimental results on Ultra-High Energy Cosmic Rays (UHECR's) including measurements of the features of the spectrum, the composition of the primary particle flux and the search for anisotropy in event arrival direction. We find that while there is a general consensus on the features in the spectrum -- the Second Knee, the Ankle, and (to a lesser extent) the GZK Cutoff -- there is little consensus on the composition of the primaries that accompany these features. This lack of consensus on the composition makes interpretation of the agreed upon features problematic. There is also little direct evidence about potential sources of UHECRs, as early reports of arrival direction anisotropies have not been confirmed in independent measurements.Comment: 46 pages, 30 figures. Topical Review to appear in J. Physics

3D Reconstruction of Small Solar System Bodies using Rendered and Compressed Images

Synthetic image generation and reconstruction of Small Solar System Bodies and the influence of compression is becoming an important study topic because of the advent of small spacecraft in deep space missions. Most of these missions are fly-by scenarios, for example in the Comet Interceptor mission. Due to limited data budgets of small satellite missions, maximising scientific return requires investigating effects of lossy compression. A preliminary simulation pipeline had been developed that uses physics-based rendering in combination with procedural terrain generation to overcome limitations of currently used methods for image rendering like the Hapke model. The rendered Small Solar System Body images are combined with a star background and photometrically calibrated to represent realistic imagery. Subsequently, a Structure-from-Motion pipeline reconstructs three-dimensional models from the rendered images. In this work, the preliminary simulation pipeline was developed further into the Space Imaging Simulator for Proximity Operations software package and a compression package was added. The compression package was used to investigate effects of lossy compression on reconstructed models and the possible amount of data reduction of lossy compression to lossless compression. Several scenarios with varying fly-by distances ranging from 50 km to 400 km and body sizes of 1 km and 10 km were simulated and compressed with lossless and several quality levels of lossy compression using PNG and JPEG 2000 respectively. It was found that low compression ratios introduce artefacts resembling random noise while high compression ratios remove surface features. The random noise artefacts introduced by low compression ratios frequently increased the number of vertices and faces of the reconstructed three-dimensional model

Field measurements of horizontal forward motion velocities of terrestrial dust devils: towards a proxy for ambient winds on Mars and Earth

Dust devils – convective vortices made visible by the dust and debris they entrain – are common in arid environments and have been observed on Earth and Mars. Martian dust devils have been identified both in images taken at the surface and in remote sensing observations from orbiting spacecraft. Observations from landing craft and orbiting instruments have allowed the dust devil translational forward motion (ground velocity) to be calculated, but it is unclear how these velocities relate to the local ambient wind conditions, for (i) only model wind speeds are generally available for Mars, and (ii) on Earth only anecdotal evidence exists that compares dust devil ground velocity with ambient wind velocity. If dust devil ground velocity can be reliably correlated to the ambient wind regime, observations of dust devils could provide a proxy for wind speed and direction measurements on Mars. Hence, dust devil ground velocities could be used to probe the circulation of the martian boundary layer and help constrain climate models or assess the safety of future landing sites. We present results from a field study of terrestrial dust devils performed in the southwest USA in which we measured dust devil horizontal velocity as a function of ambient wind velocity. We acquired stereo images of more than a hundred active dust devils and recorded multiple size and position measurements for each dust devil. We used these data to calculate dust devil translational velocity. The dust devils were within a study area bounded by 10 m high meteorology towers such that dust devil speed and direction could be correlated with the local ambient wind speed and direction measurements. Daily (10:00 to 16:00 local time) and two-hour averaged dust devil ground speeds correlate well with ambient wind speeds averaged over the same period. Unsurprisingly, individual measurements of dust devil ground speed match instantaneous measurements of ambient wind speed more poorly; a 20-minute smoothing window applied to the ambient wind speed data improves the correlation. In general, dust devils travel 10-20% faster than ambient wind speed measured at 10 m height, suggesting that their ground speeds are representative of the boundary layer winds a few tens of meters above ground level. Dust devil ground motion direction closely matches the measured ambient wind direction. The link between ambient winds and dust devil ground velocity demonstrated here suggests that a similar one should apply on Mars. Determining the details of the martian relationship between dust devil ground velocity and ambient wind velocity might require new in-situ or modelling studies but, if completed successfully, would provide a quantitative means of measuring wind velocities on Mars that would otherwise be impossible to obtain

Extraction of spatial information from sterioscopic SAR images

Synthetic Aperture Radar (SAR) is now widely used for generating Digital Elevation Models (DEMs) and has advantages over optical data in terms of availability as it allows all-day and all-weather operations. The stereoscopic SAR method, which allows direct extraction of spatial information in three-dimensional space, has been established for decades. However, the traditional stereoscopic methods developed for SAR data depend on many human operations and need ground control points (GCPs), to set up geometric models. The aims of the thesis are not only to propose a refined rigorous stereoscopic SAR method and a new error model to predict theoretic errors, but also to achieve a higher level of automation and accuracy. By using a weighting matrix, which is derived by considering different observations in the space intersection algorithm, the minimal number of the GCPs required for the refined algorithm is only two. To achieve a high degree of automation, an optimized strategy of parameter selection for the pyramidal image correlation scheme employing a region-growing technique has been proposed. This avoids a trial-and-error approach to produce digital parallax data from the same-side SAR image pairs. A new method to derive GCPs automatically has been developed using a SAR image simulation technique, under the condition that a known DEM chip is available, to minimize human interventions and operator error. The proposed method for providing GCPs and the DEMs generated from space intersection have been incorporated into the procedures for geocoding SAR images to validate the proposed algorithms. The results derived show that the stereoscopic SAR data can be applied to geometric rectification in flat-to-moderate areas, and other applications of extraction of spatial information are promising