Integrated vision-based attitude estimation

For many years, computer vision has lured researchers with promises of a low-cost, passive, lightweight and information-rich sensor suitable for navigation purposes. The prime difficulty in vision-based navigation is that the navigation solution will continually drift with time unless external information is available, whether it be cues from the appearance of the scene, a map of features (whether built online or known a priori), or from an externally-referenced sensor. It is not merely position that is of interest in the navigation problem. Attitude (i.e. the angular orientation of a body with respect to a reference frame) is integral to a visionbased navigation solution and is often of interest in its own right (e.g. flight control). This thesis examines vision-based attitude estimation in an aerospace environment, and two methods are proposed for constraining drift in the attitude solution; one through a novel integration of optical flow and the detection of the sky horizon, and the other through a loosely-coupled integration of Visual Odometry and GPS position measurements. In the first method, roll angle, pitch angle and the three aircraft body rates are recovered though a novel method of tracking the horizon over time and integrating the horizonderived attitude information with optical flow. An image processing front-end is used to select several candidate lines in a image that may or may not correspond to the true horizon, and the optical flow is calculated for each candidate line. Using an Extended Kalman Filter (EKF), the previously estimated aircraft state is propagated using a motion model and a candidate horizon line is associated using a statistical test based on the optical flow measurements and location of the horizon in the image. Once associated, the selected horizon line, along with the associated optical flow, is used as a measurement to the EKF. To evaluate the accuracy of the algorithm, two flights were conducted, one using a highly dynamic Uninhabited Airborne Vehicle (UAV) in clear flight conditions and the other in a human-piloted Cessna 172 in conditions where the horizon was partially obscured by terrain, haze and smoke. The UAV flight resulted in pitch and roll error standard deviations of 0.42° and 0.71° respectively when compared with a truth attitude source. The Cessna 172 flight resulted in pitch and roll error standard deviations of 1.79° and 1.75° respectively. In the second method for estimating attitude, a novel integrated GPS/Visual Odometry (GPS/VO) navigation filter is proposed, using a structure similar to a classic looselycoupled GPS/INS error-state navigation filter. Under such an arrangement, the error dynamics of the system are derived and a Kalman Filter is developed for estimating the errors in position and attitude. Through similar analysis to the GPS/INS problem, it is shown that the proposed filter is capable of recovering the complete attitude (i.e. pitch, roll and yaw) of the platform when subjected to acceleration not parallel to velocity for both the monocular and stereo variants of the filter. Furthermore, it is shown that under general straight line motion (e.g. constant velocity), only the component of attitude in the direction of motion is unobservable. Numerical simulations are performed to demonstrate the observability properties of the GPS/VO filter in both the monocular and stereo camera configurations. Furthermore, the proposed filter is tested on imagery collected using a Cessna 172 to demonstrate the observability properties on real-world data. The proposed GPS/VO filter does not require additional restrictions or assumptions such as platform-specific dynamics, map-matching, feature-tracking, visual loop-closing, gravity vector or additional sensors such as an IMU or magnetic compass. Since no platformspecific dynamics are required, the proposed filter is not limited to the aerospace domain and has the potential to be deployed in other platforms such as ground robots or mobile phones.

Frank-starling control of a left ventricular assist device

A physiological control system was developed for a rotary left ventricular assist device (LVAD) in which the target pump flow rate (LVADQ) was set as a function of left atrial pressure (LAP), mimicking the Frank-Starling mechanism. The control strategy was implemented using linear PID control and was evaluated in a pulsatile mock circulation loop using a prototyped centrifugal pump by varying pulmonary vascular resistance to alter venous return. The control strategy automatically varied pump speed (2460 to 1740 to 2700 RPM) in response to a decrease and subsequent increase in venous return. In contrast, a fixed-speed pump caused a simulated ventricular suction event during low venous return and higher ventricular volumes during high venous return. The preload sensitivity was increased from 0.011 L/min/mmHg in fixed speed mode to 0.47L/min/mmHg, a value similar to that of the native healthy heart. The sensitivity varied automatically to maintain the LAP and LVADQ within a predefined zone. This control strategy requires the implantation of a pressure sensor in the left atrium and a flow sensor around the outflow cannula of the LVAD. However, appropriate pressure sensor technology is not yet commercially available and so an alternative measure of preload such as pulsatility of pump signals should be investigated.

Prandtl number scaling of natural convection of the flow on an inclined at plate due to uniform surface heat flux

A new scaling analysis has been performed for the unsteady natural convection boundary layer under a downward facing inclined plate with uniform heat flux. The development of the thermal or viscous boundary layers may be classified into three distinct stages including an early stage, a transitional stage and a steady stage, which can be clearly identified in the analytical as well as numerical results. Earlier scaling shows that the existing scaling laws of the boundary layer thickness, velocity and steady state time scales for the natural convection flow on a heated plate of uniform heat flux provide a very poor prediction of the Prandtl number dependency. However, those scalings performed very well with Rayleigh number and aspect ratio dependency. In this study, a modifed Prandtl number scaling has been developed using a triple-layer integral approach for Pr > 1. It is seen that in comparison to the direct numerical simulations, the new scaling performs considerably better than the previous scaling.

Mastering the canonical loop of serine protease inhibitors : enhancing potency by optimising the internal hydrogen bond network

Background Canonical serine protease inhibitors commonly bind to their targets through a rigid loop stabilised by an internal hydrogen bond network and disulfide bond(s). The smallest of these is sunflower trypsin inhibitor (SFTI-1), a potent and broad-range protease inhibitor. Recently, we re-engineered the contact β-sheet of SFTI-1 to produce a selective inhibitor of kallikrein-related peptidase 4 (KLK4), a protease associated with prostate cancer progression. However, modifications in the binding loop to achieve specificity may compromise structural rigidity and prevent re-engineered inhibitors from reaching optimal binding affinity. Methodology/Principal Findings In this study, the effect of amino acid substitutions on the internal hydrogen bonding network of SFTI were investigated using an in silico screen of inhibitor variants in complex with KLK4 or trypsin. Substitutions favouring internal hydrogen bond formation directly correlated with increased potency of inhibition in vitro. This produced a second generation inhibitor (SFTI-FCQR Asn14) which displayed both a 125-fold increased capacity to inhibit KLK4 (Ki = 0.0386±0.0060 nM) and enhanced selectivity over off-target serine proteases. Further, SFTI-FCQR Asn14 was stable in cell culture and bioavailable in mice when administered by intraperitoneal perfusion. Conclusion/Significance These findings highlight the importance of conserving structural rigidity of the binding loop in addition to optimising protease/inhibitor contacts when re-engineering canonical serine protease inhibitors.

Modelling passenger flow at airport terminals : individual agent decision model for stochastic passenger behaviour

Airport system is complex. Passenger dynamics within it appear to be complicate as well. Passenger behaviours outside standard processes are regarded more significant in terms of public hazard and service rate issues. In this paper, we devised an individual agent decision model to simulate stochastic passenger behaviour in airport departure terminal. Bayesian networks are implemented into the decision making model to infer the probabilities that passengers choose to use any in-airport facilities. We aim to understand dynamics of the discretionary activities of passengers.

On the oscillatory behavior of transient Rayleigh Benard convection of air for 2D channel flow at moderate Rayleigh numbers

Unsteady numerical simulation of Rayleigh Benard convection heat transfer from a 2D channel is performed. The oscillatory behavior is attributed to recirculation of ascending and descending flows towards the core of the channel producing organized rolled motions. Variation of the parameters such as Reynolds number, channel outlet flow area and inclination of the channel are considered. Increasing Reynolds number (for a fixed Rayleigh number), delays the generation of vortices. The reduction in the outflow area leads to the later and the less vortex generation. As the time progresses, more vortices are generated, but the reinforced mean velocity does not let the eddies to enter the core of the channel. Therefore, they attach to the wall and reduce the heat transfer area. The inclination of the channel (both positive and negative) induces the generated vortices to get closer to each other and make an enlarged vortex.

Evidence for carbonate platform failure during rapid sea-level rise; ca 14 000 year old bioclastic flow deposits in the Lesser Antilles

Bioclastic flow deposits offshore from the Soufrie`re Hills volcano on Montserrat in the Lesser Antilles were deposited by the largest volume sediment flows near this active volcano in the last 26 kyr. The volume of these deposits exceeds that of the largest historic volcanic dome collapse in the world, which occurred on Montserrat in 2003. These flows were most probably generated by a large submarine slope failure of the carbonate shelf comprising the south west flank of Antigua or the east flank of Redonda; adjacent islands that are not volcanically active. The bioclastic flow deposits are relatively coarse-grained and either ungraded or poorly graded, and were deposited by non cohesive debris flow and high density turbidity currents. The bioclastic deposit often comprises multiple sub-units that cannot be correlated between core sites; some located just 2 km apart. Multiple sub-units in the bioclastic deposit result from either flow reflection, stacking of multiple debris flow lobes, and/or multi-stage collapse of the initial landslide. This study provides unusually precise constraints on the age of this mass flow event that occurred at ca 14 ka. Few large submarine landslides have been well dated, but the slope failures that have been dated are commonly associated with periods of rapid sea-level change.

Anatomy of a submarine pyroclastic flow and associated turbidity current : July 2003 dome collapse, Soufrière Hills volcano, Montserrat, West Indies

The 12 to 13 July 2003 andesite lava dome collapse at the Soufrière Hills volcano, Montserrat, provides the first opportunity to document comprehensively both the sub-aerial and submarine sequence of events for an eruption. Numerous pyroclastic flows entered the ocean during the collapse, depositing approximately 90% of the total material into the submarine environment. During peak collapse conditions, as the main flow penetrated the air–ocean interface, phreatic explosions were observed and a surge cloud decoupled from the main flow body to travel 2 to 3 km over the ocean surface before settling. The bulk of the flow was submerged and rapidly mixed with sea water forming a water-saturated mass flow. Efficient sorting and physical differentiation occurred within the flow before initial deposition at 500 m water depth. The coarsest components (∼60% of the total volume) were deposited proximally from a dense granular flow, while the finer components (∼40%) were efficiently elutriated into the overlying part of the flow, which evolved into a far-reaching turbidity current.

Fusing loop detector and probe vehicle data to estimate travel time statistics on signalized urban networks

This article presents a methodology that integrates cumulative plots with probe vehicle data for estimation of travel time statistics (average, quartile) on urban networks. The integration reduces relative deviation among the cumulative plots so that the classical analytical procedure of defining the area between the plots as the total travel time can be applied. For quartile estimation, a slicing technique is proposed. The methodology is validated with real data from Lucerne, Switzerland and it is concluded that the travel time estimates from the proposed methodology are statistically equivalent to the observed values.

Modelling the flow of cane constituents in the milling process

A model has been developed to track the flow of cane constituents through the milling process. While previous models have tracked the flow of fibre, brix and water through the process, this model tracks the soluble and insoluble solid cane components using modelling theory and experiment data, assisting in further understanding the flow of constituents into mixed juice and final bagasse. The work provided an opportunity to understand the factors which affect the distribution of the cane constituents in juice and bagasse. Application of the model should lead to improvements in the overall performance of the milling train.

In vitro and in vivo studies on protease-activated receptor-2

Protease-activated receptor-2 (PAR2) is a G protein coupled receptor (GPCR) that is activated by proteolytic cleavage of its amino terminal domain by trypsin-like serine proteases. Cleavage of this receptor exposes a neoepitope, termed the tethered ligand (TL), which binds intramolecularly within the receptor to stimulate signal transduction via coupled G proteins. PAR2-mediated signal transduction is also experimentally stimulated by hexapeptides (agonist peptides; APs) that are homologous to the TL sequence. Due to the irreversible nature of PAR2 proteolysis, downstream signal transduction is tightly regulated. Following activation, PAR2 is rapidly uncoupled from downstream signalling by the post-translational modifications phosphorylation and ubiquination which facilitate interactions with â- arrestin. This scaffolding protein couples PAR2 to the internalisation machinery initiating its desensitisation and trafficking through the early and late endosomes followed by receptor degradation. PAR2 is widely expressed in mammalian tissues with key roles for this receptor in cardiovascular, respiratory, nervous and musculoskeletal systems. This receptor has also been linked to pathological states with aberrant expression and signalling noted in several cancers. In prostate cancer, PAR2 signalling induces migration and proliferation of tumour derived cell lines, while elevated receptor expression has been noted in malignant tissues. Importantly, a role for this receptor has also been suggested in prostate cancer bone metastasis as coexpression of PAR2 and a proteolytic activator has been demonstrated by immunohistochemical analysis. Based on these data, the primary focus of this project has been on two aspects of PAR2 biology. The first is characterisation of cellular mechanisms that regulate PAR2 signalling and trafficking. The second aspect is the role of this receptor in prostate cancer bone metastasis. In addition, to permit these studies, it was first necessary to evaluate the specificity of the commercially available anti-PAR2 antibodies SAM11, C17, N19 and H99. The evaluation of the four commercially available antibodies was assessed using four techniques: immunoprecipitation; Western blot analysis; immunofluorescence; and flow cytometry. These approaches demonstrated that three of the antibodies efficiently detect ectopically expressed PAR2 by each of these techniques. A significant finding from this study was that N19 was the only antibody able to specifically detect N-glycosylated endogenous PAR2 by Western blot analysis. This analysis was performed on lysates from prostate cancer derived cell lines and tissue derived from wildtype and PAR2 knockout mice. Importantly, further evaluation demonstrated that this antibody also efficiently detects endogenous PAR2 at the cell surface by flow cytometry. The anti-PAR2 antibody N19 was used to explore the in vitro role of palmitoylation, the post-translational addition of palmitate, in PAR2 signalling, trafficking, cell surface expression and desensitization. Significantly, use of the palmitoylation inhibitor 2-bromopalmitate indicated that palmitate addition is important in trafficking of PAR2 endogenously expressed by prostate cancer cell lines. This was supported by palmitate labelling experiments using two approaches which showed that PAR2 stably expressed by CHO cells is palmitoylated and that palmitoylation occurs on cysteine 361. Another key finding from this study is that palmitoylation is required for optimal PAR2 signalling as Ca2+ flux assays indicated that in response to trypsin agonism, palmitoylation deficient PAR2 is ~9 fold less potent than wildtype receptor with a reduction of about 33% in the maximum signal induced via the mutant receptor. Confocal microscopy, flow cytometry and cell surface biotinylation analyses demonstrated that palmitoylation is required for efficient cell surface expression of PAR2. Importantly, this study also identified that palmitoylation of this receptor within the Golgi apparatus is required for efficient agonist-induced rab11amediated trafficking of PAR2 to the cell surface. Interestingly, palmitoylation is also required for receptor desensitization, as agonist-induced â-arrestin recruitment and receptor degradation were markedly reduced in CHO-PAR2-C361A cells compared with CHO-PAR2 cells. Collectively, these data provide new insights on the life cycle of PAR2 and demonstrate that palmitoylation is critical for efficient signalling, trafficking, cell surface localization and degradation of this receptor. This project also evaluated PAR2 residues involved in ligand docking. Although the extracellular loop (ECL)2 of PAR2 is known to be required for agonist-induced signal transduction, the binding pocket for receptor agonists remains to be determined. In silico homology modelling, based on a crystal structure for the prototypical GPCR rhodopsin, and ligand docking were performed to identify PAR2 transmembrane (TM) amino acids potentially involved in agonist binding. These methods identified 12 candidate residues that were mutated to examine the binding site of the PAR2 TL, revealed by trypsin cleavage, as well as of the soluble ligands 2f-LIGRLO-NH2 and GB110, which are both structurally based on the AP SLIGRLNH2. Ligand binding was evaluated from the impact of the mutated residues on PAR2-mediated calcium mobilisation. An important finding from these experiments was that mutation of residues Y156 and Y326 significantly reduced 2f-LIGRLO-NH2 and GB110 agonist activity. L307 was also important for GB110 activity. Intriguingly, mutation of PAR2 residues did not alter trypsin-induced signalling to the same extent as for the soluble agonists. The reason for this difference remains to be further examined by in silico and in vitro experimentation and, potentially, crystal structure studies. However, these findings identified the importance of TM domains in PAR2 ligand docking and will enhance the design of both PAR2 agonists and potentially agents to inhibit signalling (antagonists). The potential importance of PAR2 in prostate cancer bone metastasis was examined using a mouse model. In patients, prostate cancer bone metastases cause bone growth by disrupting bone homeostasis. In an attempt to mimic prostate cancer growth in bone, PAR2 responsive 22Rv1 prostate cancer cells, which form mixed osteoblastic and osteolytic lesions, were injected into the proximal aspect of mouse tibiae. A role for PAR2 was assessed by treating these mice with the recently developed PAR2 antagonist GB88. As controls, animals bearing intra-tibial tumours were also treated with vehicle (olive oil) or the prostate cancer chemotherapeutic docetaxel. The effect of these treatments on bone was examined radiographically and by micro-CT. Consistent with previous studies, 22Rv1 tumours caused osteoblastic periosteal spicule formation and concurrent osteolytic bone loss. Significantly, blockade of PAR2 signalling reduced the osteoblastic and osteolytic phenotype of 22Rv1 tumours in bone. No bone defects were detected in mice treated with docetaxel. These qualitative data will be followed in the future by quantitative micro-CT analysis as well as histology and histomorphometry analysis of already collected tissues. Nonetheless, these preliminary experiments highlight a potential role for PAR2 in prostate cancer growth in bone. In summary, in vitro studies have defined mechanisms regulating PAR2 activation, downstream signalling and trafficking and in vivo studies point to a potential role for this receptor in prostate cancer bone metastasis. The outcomes of this project are that a greater understanding of the biology of PAR2 may lead to the development of strategies to modulate the function of this receptor in disease.