COMPUTATION OF DISTANCES OF LONG ARCS FOR RADIO PURPOSES

Abstract

In the E.S.R., No. 59, Vol. VIII, of January, 1946, I gave a formula which I had worked out to give a rapid and easy means of computing the lengths of long arcs, up to 1000 kilometres, between two points whose latitude and longitudes are known on a definite figure of the earth.     

THE ACCURAOY OF ASTRONOMIOAL OBSERVATIONS FOR POSITION

Abstract

The following paper describes some of the results obtained on it, tour of the Katsina and Kano Provinces of Northern Nigeria between November 1953 and February 1954. The purpose of the tour was to obtain latitudes and longitudes by astronomical observations for air photo mapping control in the northern parts of these provinces where triangulation is not available. The area is mostly flat and sandy.     

MARGINAL SCALES OF LATITUDE AND LONGITUDE ON TRANSVERSE MERCATOR MAPS

Abstract

The problem of computing marginal scales of latitude and longitude on a rectangular map on the Transverse Mercator projection, where the sheet boundaries are projection co-ordinate lines, may be solved in various ways. A simple method is to compute the latitudes and longitudes of the four corners of the sheet, and then, assuming a constant scale, to interpolate the parallels and meridians between these corner values. Although it is probably sufficiently accurate for practical purposes, this method is not precise. It is not difficult to adapt the fundamental formulce of the projection to give a direct solution of the problem.     

LONG LINES ON THE EARTH: THE DIRECT PROBLEM

Abstract

19. Formulae.—In Nos. 6, vol. i, and 9, vol. ii, pp. 259 and 156, there has been described a new method for dealing with long geodesics on the earth's surface. There the so-called “inverse” problem has claimed first attention: given the latitudes and longitudes of the extremities of a geodesic, to find its length and terminal azimuths. It remains to discuss the “direct” problem : a geodesic of given length starts on a given azimuth from a station of known latitude and longitude; to find the latitude and longitude of its extremity and the azimuth thereat. The solution of this direct problem demands a certain recasting of the formulae previously given. In order of working the several expressions now assume the forms below.     

GEOGRAPHICAL POSITIONS IN MAURITIUS AND RODRIGUEZ

Abstract

The desirability of determining as accurately as possible the latitudes and longitudes of a number of points in Mauritius appears to have first appealed to the French astronomer M. L'Abbé de la Caille, who was sent to the island by the French East India Company in 1753. His observatory was situated 4,730 feet east and 2,610 feet north of Port Louis Time Ball, and in addition to determining the geographical position of this. 0bservatory (the house, later demolished, of a Mr Mabile, where Mr D'Après had made observations the previous year)—

Latitude 20° 09′ 42″ S.; Longitude 55° 08′ 15″ E. of Paris—,

he succeeded in effecting a triangulation of the island. Four bases were measured with wooden scales which previously had been compared with an iron toise (6.394 feet—the French fathom) approved by the Academy of Sciences, Paris; a I4-inch quadrant fitted with micrometer was used to measure the angles.     

User Desires and Graphics Capability in the Academic Environment

Abstract

In this paper the cartographic grids engraved on two antique instruments from Iran for finding the sacred direction and distance to Mecca are discussed. It appears that these grids can be well explained in terms of the Mecca-centred retro-azimuthal orthographic projection described in 1968 by J. E. Jackson. In this projection the lines of constant latitude reduce to a set of ellipses with their major axes parallel to the equator and the lines of the constant longitudes reduce to a set of non-equidistant straight lines parallel to the north-south direction. It is shown that the curves actually engraved on the instrument conform to this projection and can be fairly easily constructed. This interpretation of the grid on the Iranian instruments stands in contrast with another explanation, recently proposed by King (1999), which is based on medieval Arabic concepts such as the so-called 'methods of the zijes'. Insufficiently accurate workmanship makes it impossible to distinguish between the two explanations through the study of the instruments themselves. The newly gained insight into the projection itself, however, shows that a direct relation between the Iranian maps and Islamic mappings insight knows from the ninth century, as suggested by King, does not exit. Thus, it is concluded that it is as yet completely unknown when and where the very idea behind the Iranian cartographic grid was first conceived, and that the quest for their historical background is still open.     

A FIFTY-YEARS RETROSPECT

Abstract

About this time an excellent Instructor in Surveying was appointed to the School of Military Engineering in the person of Major A. C. MacDonnell. He had served in India,—though not on the Survey of India,—and, being well acquainted with the excellent Indian frontier survey methods, resolved to introduce them into the course at Chatham. So he started using the system of computing latitudes and longitudes from trigonometrical data by Puissant's formulæ, in the form used by the Survey of India. But he had reckoned without his host, the higher authorities. His dreadful deed became known, and the matter was referred to three eminent officers for their opinion. The three officers were Sir Charles Wilson, Director of Military Education, Sir John Ardagh, Director of Military Intelligence, and Sir John Farquharson, Director-General of the Ordnance Survey; none of the three had had any personal acquaintance with the method in question, although two of them had directed the Ordnance Survey, and Sir Charles Wilson in the sixties had carried out some very interesting surveys in Palestine and Sinai. Well, these three distinguished officers solemnly condemned the Indian method as being unsuitable for use at Chatham, and MacDonnell had to revert to more primitive ways, which later on would have made impossible the conduct of a properly managed boundary commission or such surveys as that of the Orange Free State, Uganda, or Northern Sinai, or much of the technical work on the Western Front during the War. And that was that.     

The terrestrial and lunar reference frame in lunar laser ranging

The initial value problem and the stability of solution in the determination of the coordinates of three observing stations and four retro-reflectors by lunar laser ranging are discussed. Practical iterative computations show that the station coordinates can be converged to about 1 cm, but there will be a slight discrepancy of the longitudinal components computed by various analysis centers or in different years. There are several factors, one of which is the shift of the right ascension of the Moon, caused by the orientation deviation of the adopted lunar ephemeris, which can make the longitudinal components of all observing stations rotate together along the longitudinal direction with same angle. Additionally, the frame of selenocentric coordinates is stable, but a variation or adjustment of lunar third-degree gravitational coefficients will cause a simultaneous shift along the reflectors' longitudes or rotation around the Y axis. Received: 21 August 1996 / Accepted: 17 November 1998     

A statistical filter for geodetic observations

A method for filtering of geodetic observationwhich leaves the final result normally distributed, is presented. Furthermore, it is shown that if you sacrifice100.a% of all the observations you may be (1−β).100% sure that a gross error of the size Δ is rejected. Another and, may be intuitively, more appealing method is presented; the two methods are compared and it is shown why Method 1 should be preferred to Method 2 for geodetic purposes. Finally the two methods are demonstrated in some numerical examples.     

On non-linear low–low SST observation equations for the determination of the geopotential based on an analytical solution

In satellite data analysis, one big advantage of analytical orbit integration, which cannot be overestimated, is missed in the numerical integration approach: spectral analysis or the lumped coefficient concept may be used not only to design efficient algorithms but overall for much better insight into the force-field determination problem. The lumped coefficient concept, considered from a practical point of view, consists of the separation of the observation equation matrix A=BT into the product of two matrices. The matrix T is a very sparse matrix separating into small block-diagonal matrices connecting the harmonic coefficients with the lumped coefficients. The lumped coefficients are nothing other than the amplitudes of trigonometric functions depending on three angular orbital variables; therefore, the matrix N=B ^T B will become for a sufficient length of a data set a diagonal dominant matrix, in the case of an unlimited data string length a strictly diagonal one. Using an analytical solution of high order, the non-linear observation equations for low–low SST range data can be transformed into a form to allow the application of the lumped concept. They are presented here for a second-order solution together with an outline of how to proceed with data analysis in the spectral domain in such a case. The dynamic model presented here provides not only a practical algorithm for the parameter determination but also a simple method for an investigation of some fundamental questions, such as the determination of the range of the subset of geopotential coefficients which can be properly determined by means of SST techniques or the definition of an optimal orbital configuration for particular SST missions. Numerical results have already been obtained and will be published elsewhere. Received: 15 January 1999 / Accepted: 30 November 1999     

Higher-order ionospheric effects in GPS time and frequency transfer

In high-precision geodetic time and frequency transfer, which requires precise modeling of code and carrier phase GPS data, the ionosphere-free combinations P ₃ and L ₃ of the codes and carrier phases, made on the two GPS frequencies, are used to remove the first-order ionospheric effect. We quantify the impact of the residual second- and third-order ionospheric effects on geodetic time and frequency transfer solutions for continental and intercontinental baselines. All time transfer computations are done using the ATOMIUM software, developed at the Royal Observatory of Belgium. In order to avoid contamination by some imperfect modeling of the second- and third-order ionospheric effects in the satellite clock products, only single-difference, common-view processing is used, based on code and carrier phase measurements. The results are shown for weak and strong solar activity, as well as for particular epochs of ionospheric storms. Second-order ionospheric delays can lead to corrections up to about 10 ps in the common-view clock solution of intercontinental baselines with very different longitudes. However, realistic values of the geomagnetic field in the ionosphere are required to assess the amplitude of second-order ionospheric effects in time and frequency transfer during an ionospheric storm.     

THE CONFORMAL TRANSFORMATION OF A NETWORK OF TRIANGULATION

Abstract

In a recent issue of this Review, an example is given of the conformal transformation of a network of triangulation using Newton's interpolation formula with divided differences. While the application of the method appears to be new, attention should be drawn to the fact that Kruger employed Lagrange's interpolation formula in a discussion and extension of the Schols method in a paper which was published in the Zeitschrift für Vermessungswesen in 1896. A reference to this paper was given at the end of the paper, “Adjustment of the Secondary Triangulation of South Africa”, published in a previous issue of the E.S.R. (iv, 30, 480).     

Gravimetric and astro-geodetic geoids and mean free-air anomalies in India

A preliminary gravimetric geoid with respect to the International Spheroid and the latest astro-geodetic geoid computed on the Everest and International Spheroids are given in the form of undulation maps over the Indian Sub—continent. 1⁰x1⁰ mean free-air anomalies (modified) on the Geodetic Reference System, 1967 (GRS-67) are also given for the whole country in the form of a chart. For the purpose of computing the gravimetric geoid, 5⁰x5⁰ mean free-air anomalies were used outside the area bounded by latitudes 0⁰ to 40⁰ N and longitudes 60⁰ to 100⁰ E and 1⁰x1⁰ mean free-air anomalies within these limits. The anomalies partly computed by Survey of India and mostly collected from other sources (such as B.G.I.) were utilised for this purpose. The astro-geodetic geoid is based on the astronomical data observed in India up to 1978.     

Statistical analysis of discrepancies in high precision levelling

An investigation was made of the behaviour of the variable (where ρ_ij are the discrepancies between the direct and reverse measurements of the height of consecutive bench marks and theR _ij are their distance apart) in a partial net of the Italian high precision levelling of a total length of about1.400 km. The methods of analysis employed were in general non-parametric individual and cumulative tests; in particular randomness, normality and asymmetry tests were carried out. The computers employed wereIBM/7094/7040. From the results evidence was obtained of the existence of an asymmetry in respect to zero of thex _ij confirming the well-known results given firstly by Lallemand. A new result was obtained from the tests of randomness which put in evidence trends of the mean values of thex _ij and explained some anomalous behaviours of the cumulative discrepancy curves. The extension of this investigation to a broader net possibly covering other national nets would be very useful to get a deeper insight into the behaviour of the errors in high precision levelling. Ad hoc programs for electronic computers are available to accomplish this job quickly. Presented at the 14^th International Assembly of Geodesy (Lucerne, 1967).     

COMPUTATION OF DISTANCES OF LONG ARCS FOR RADIO PURPOSES

Abstract

When we require the distance between two points on a spheroid, there are at least six different lines between the points which might be taken.     

Mixed Integer-Real Valued Adjustment (IRA) Problems: GPS Initial Cycle Ambiguity Resolution by Means of the LLL Algorithm

In order to achieve to GPS solutions of first-order accuracy and integrity, carrier phase observations as well as pseudorange observations have to be adjusted with respect to a linear/linearized model. Here the problem of mixed integer-real valued parameter adjustment (IRA) is met. Indeed, integer cycle ambiguity unknowns have to be estimated and tested. At first we review the three concepts to deal with IRA: (i) DDD or triple difference observations are produced by a properly chosen difference operator and choice of basis, namely being free of integer-valued unknowns (ii) The real-valued unknown parameters are eliminated by a Gauss elimination step while the remaining integer-valued unknown parameters (initial cycle ambiguities) are determined by Quadratic Programming and (iii) a RA substitute model is firstly implemented (real-valued estimates of initial cycle ambiguities) and secondly a minimum distance map is designed which operates on the real-valued approximation of integers with respect to the integer data in a lattice. This is the place where the integer Gram-Schmidt orthogonalization by means of the LLL algorithm (modified LLL algorithm) is applied being illustrated by four examples. In particular, we prove that in general it is impossible to transform an oblique base of a lattice to an orthogonal base by Gram-Schmidt orthogonalization where its matrix enties are integer. The volume preserving Gram-Schmidt orthogonalization operator constraint to integer entries produces “almost orthogonal” bases which, in turn, can be used to produce the integer-valued unknown parameters (initial cycle ambiguities) from the LLL algorithm (modified LLL algorithm). Systematic errors generated by “almost orthogonal” lattice bases are quantified by A. K. Lenstra et al. (1982) as well as M. Pohst (1987). The solution point of Integer Least Squares generated by the LLL algorithm is = (L')⁻¹[L'◯] ∈ ℤ ^m where L is the lower triangular Gram-Schmidt matrix rounded to nearest integers, [L], and = [L'◯] are the nearest integers of L'◯, ◯ being the real valued approximation of z ∈ ℤ ^m , the m-dimensional lattice space Λ. Indeed due to “almost orthogonality” of the integer Gram-Schmidt procedure, the solution point is only suboptimal, only close to “least squares.” ? 2000 John Wiley & Sons, Inc.     

The Earliest Maps of Southampton

Abstract

In editing the early maps of Southampton for reproduction, the City Archivist has made a detailed study of their origins. This account demonstrates the precautions which must be taken in dealing with old maps and plans.     

On the night errors in astronomical determinations

Summary The discrepancy between precision and accuracy in astronomical determinations is usually explained in two ways: on the one hand by ostensible large refraction anomalies and on the other hand by variable instrumental errors which are systematic over a certain interval of time and which are mainly influenced by temperature.In view of the research of several other persons and the author’s own investigations, the authors are of the opinion that the large night-errors of astronomical determinations are caused by variable, systematic instrumental errors dependent on temperature. The influence of refraction anomalies is estimated to be smaller than 0″.1 for most of the field stations. The possibility of determining the anomalous refraction from the observations by the programme given by Prof. Pavlov and Anderson has also been investigated. The precision of the determination of the anomalous refraction is good as long as no other systematic error working in a similar way is present.The results, which are interpreted as an effect of the anomalous refraction by Pavlov and Sergijenko, could also be interpreted as a systematic instrumental error. It is furthermore maintained thatthe latitude and longitude of a field station can be determined in a few hours of one night if the premisses given in [3, p.68]are kept. It has been deplored that the determination of the azimuth has not been given the necessary attention. It is therefore proposed to intensify the research on this problem. The profession has been called upon to acquaint itself better with the valuable possibilities of astronomical determinations and to apply them in a useful and appropriate manner. At the same time, attention has been called to the possibility of improving astronomical determinations with regard to accuracy as well as effectiveness.     

Analysis on the global morphology of stratospheric gravity wave activity deduced from the COSMIC GPS occultation profiles

Monthly mean global morphologies of potential energy density E _p from stratospheric gravity waves are revealed by observations of COSMIC GPS radio occultation. The E _p is obtained from vertical wavelengths ranging from 2 to 10 km over cells of 1° × 2° in latitude and longitude. The computed values confirm previous results and obtain new ones. The large gravity wave E _p values found in the tropics between 25°N and 25°S could be mainly due to the strong tropical cumulus convection; July values are larger than those for January (2007). In mid and high latitudes, the most prominent features of the northern winter hemisphere are the enhanced densities above the Eurasian continent and the North Atlantic and the depressed E _p values above the North Pacific and North America for which topography, wind sources and wind filtering may be responsible. In southern winter hemisphere, large E _p values are found around 180° and 300° longitudes that are likely due to the topography of the Antarctic plateau, the Antarctic Peninsula and South America. Enhanced E _p values are found over Scandinavia. However, there is no clear evidence to show that gravity waves are localized over the Rocky Mountains, the Himalayas and the Andes. Topography and planetary wave modulations are proposed to interpret the large-scale longitudinal variations and inter-hemisphere asymmetry of the GW activity.