MACHINE METHOD OF INTERSECTION

Abstract

The following method will be found better and quicker than the usual logarithmic process in computing the co-ordinates of intersected points in minor triangulation and traverse work. Let A and B be two stations whose co-ordinates (x ₁ y ₁), (x ₂ y ₂) are known. Let P be an intersected point whose co-ordinates (x, y) we wish to determine. Let α and β be the observed angles at A and B respectively.     

MINOR TRIANGULATION ON THE TRANSVERSE MERCATOR PROJECTION

Abstract

In a previous article in this Review, the writer endeavoured to show that chains of minor triangulation could be adjusted by plane rectangular co-ordinates ignoring the spherical form of the earth with little loss of accuracy, provided that the two ends were held fixed in position. It was demonstrated that the plane co-ordinates produced by the rigorous adjustment between the fixed starting and closing sides, differ by only a comparatively small amount from the projection co-ordinates produced by a rigorous adjustment on the Transverse Mercator projection. The saving in time when computing by plane co-ordinates as opposed to rigorous computation on the projection by any method will be apparent to any computer with experience of both methods.     

THE CLARKE FORMULÆ FOR LATITUDE,LONGITUDE AND REVERSE AZIMUTH,WITH CORRECTIVE TERMS FOR USE ON VERY LONG LINES

Abstract

This extremely simple and elegant method of computing geographical co-ordinates, given the initial azimuth and length of line from the standpoint, was published by Col. A. R. Clarke in 1880. There is no other known method giving the same degree of accuracy with the use of only three tabulated spheroidal factors. Clarke himself regarded this as an approximate formula (vide his remark in section 5, p. 109, “Geodesy”); but as this article demonstrates, it is capable of a high degree of precision in all occupied lati tudes when certain corrections are applied to the various terms. These corrections are comparatively easy to compute, require no further spheroidal factors, and some of them may be tabulated directly once and for all.     

INTERVISIBLE FEATURES

Abstract

It is frequently required to find whether a feature A of height h ₀ will interrupt the view between two other features A₁ and A₂, of heights h ₁ and h ₂ respectively. Suppose that the right line from A₁ to A₂, whose zenith distance is ζ at A₁, has a height h at A; it is then obvious that no more is necessary than to compute h and compare it with the known height h ₀ of the feature A.     

TRANSVERSE MERCATOR FORMULAE

Abstract

The Transverse Mercator Projection, now in use for the new O.S. triangulation and mapping of Great Britain, has been the subject of several recent articles in the “Empire Surpey Review. The formulae of the projection itself have been given by various writers, from Gauss, Schreiber and Jordan to Hristow, Tardi, Lee, Hotine and others—not, it is to be regretted, with complete agreement, in all cases. For the purpose for which these formulae have hitherto been employed, in zones of restricted width and in relatively low latitudes, the completeness with which they were given was adequate, and the omission of certain smaller terms, in the fourth and higher powers of the eccentricity, was of no practical importance. In the case of the British grid, however, we have to cover a zone which must be considered as having a total width of some ten to twelve degrees of longitude at least, and extending to latitude 61 °north. This means, firstly, that terms which have as their initial co-efficients the fourth and sixth powers of the longitude ω (or of y) will be of greater magnitude than usual, and secondly that tan² ? and tan⁴ ? are likewise greatly increased. Lastly, an inspection of the formulae (as hitherto available) shows a definite tendency for the numerical co-efficients of terms to increase as the terms themselves decrease—e.g. terms in η⁴, η⁶, etc.     

Intuitive derivation of loop inverses and array algebra

Array algebra forms the general base of fast transforms and multilinear algebra making rigorous solutions of a large number (millions) of parameters computationally feasible. Loop inverses are operators solving the problem of general matrix inverses. Their derivation starts from the inconsistent linear equations by a parameter exchangeX→L ₀, where X is a set of unknown observables,A ₀ forming a basis of the so called “problem space”. The resulting full rank design matrix of parameters L₀ and its ℓ-inverse reveal properties speeding the computational least squares solution expressed in observed values . The loop inverses are found by the back substitution expressing ∧X in terms ofL through . Ifp=rank (A) ≤n, this chain operator creates the pseudoinverseA ⁺. The idea of loop inverses and array algebra started in the late60's from the further specialized case,p=n=rank (A), where the loop inverse A ₀ ⁻¹ (AA ₀ ⁻¹ )^ℓ reduces into the ℓ-inverse A^ℓ=(A^TA)⁻¹A^T. The physical interpretation of the design matrixA A ₀ ⁻¹ as an interpolator, associated with the parametersL ₀, and the consideration of its multidimensional version has resulted in extended rules of matrix and tensor calculus and mathematical statistics called array algebra.     

TABULAR AZIMUTH

Abstract

In a known latitude Φ the azimuth A of a heavenly body, of known declination δ, is determinable if its altitude h is measured by a theodolite or other instrument. The azimuth of course is that for the instant of observation.     

GRID BEARINGS AND DISTANCES ON THE CONICAL ORTHOMORPHIC PROJECTION

Abstract

When computing and adjusting traverses or secondary and tertiary triangulation in countries to which the Transverse Mercator projection has been applied, it is often more convenient to work directly in terms of rectangular co-ordinates on the projection system than it is to work in terms of geographical coordinates and then convert these later on into rectangulars. The Transverse Mercator projection is designed in the first place to cover a country whose principal extent is in latitude and hence work on it is generally confined to a belt, or helts, in which the extent of longitude on either side of the central meridian is so limited as seldom to exceed a width of much more than about 200 miles.     

CORRECTION TO BEARING FOR HEIGHT OF SIGNAL

Abstract

When a beacon B _h stands on a mountain of height h, the bearing of B _h as seen from another station A is in general affected by its elevation. The correction never exceeds one second of arc, but in primary triangulation it is not always negligible.     

A MACHINE METHOD OF RESECTION

Abstract

With the modern calculating machine in easy reach of every computer, the problem of determining the position of an occupied point from which direction observations have been made to three or more known points has become quite simple. The method outlined below is quite elegant in form and exceedingly simple on the machine. Let A, B, C be the three points whose co-ordinates (X₁Y₁), (X₂Y₂), (X₃Y₃) are known, and let (XY) be the co-ordinates of the point P which we wish to fix.     

THE ADJUSTMENT OFA BLOCK OF AERIAL TRIANGULATION EVALUATED WITH THE WILD A5

Abstract

The procedure for aerotriangulation on the Wild A5 and similar plotting instruments is well known. The first overlap is set up in absolute orientation on well spaced plan and height control and successive overlaps are set up relatively, each to the previous overlap, by eliminating want of correspondence and preserving the height agreement of points falling in the common portion of successive overlaps. When each overlap is correctly set, the co-ordinates of selected points are measured on the instrument (machine co-ordinates). These co-ordinates differ from true ground co-ordinates only in origin, azimuth and scale, provided the settings and measurements are precisely done on error-free models, precisely connected together. However, such ideal conditions are never obtained, and the errors in azimp.th, scale and height datum increase with the number of overlaps added along a strip.     

Harmonic maps

Harmonic maps are generated as a certain class of optimal map projections. For instance, if the distortion energy over a meridian strip of the International Reference Ellipsoid is minimized, we are led to the Laplace–Beltrami vector-valued partial differential equation. Harmonic functions x(L,B), y(L,B) given as functions of ellipsoidal surface parameters of Gauss ellipsoidal longitude L and Gauss ellipsoidal latitude B, as well as x(,q), y(,q) given as functions of relative isometric longitude =L–L₀ and relative isometric latitude q=Q–Q₀ gauged to a vector-valued boundary condition of special symmetry are constructed. The easting and northing {x(b,),y(b,)} of the new harmonic map is then given. Distortion energy analysis of the new harmonic map is presented, as well as case studies for (1) B[–40°,+40°], L[–31°,+49°], B₀= ±30°, L₀=9° and (2) B[46°,56°], L{[4.5°, 7.5°]; [7.5°, 10.5°]; [10.5°,13.5°]; [13.5°,16.5°]}, B₀= 51°, L₀ {6°,9°,12°,15°}.     

NORMAL EQUATIONS RESOLVED BY APPROXIMATION

Abstract

We are indebted to Professor R. V. Southwell for the approximate method of computation known as the systematic relaxation of constraints. In an article to the Empire Survey Review, 1938, Mr A. N. Black showed how Southwell's ideas could be applied to the adjustment of the co-ordinates of a point.     

MINI-MACtm—A new generation dual—Band surveyor

Summary Litton Divisions presently produce both high accuracy GPS surveyors and low—cost GPS navigation sets. Aero Service'sMACROMETER ^R Interferometric Surveyors, have become the standard against which GPS surveying equipment is measured. Litton Aero Products has developed a highly digitized, low costL ₁,C/A code GPS card set. The integration of these technologies had led to the development of a low-cost, high-precision, GPS survey system which can be configured with or without a codelessL ₂ capability. TheMINI-MAC surveying system is the first member of the new generation of GPS survey systems resulting from this joint development. The system design is described in this paper, and initial survey test results using a prototypeMINI-MAC surveying system are presented.     

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.     

LENGTH AND AZIMUTH OF LONG LINES ON THE EARTH

Abstract

In the E.S.R., viii, 59, 191–194 (January 1946), J.H. Cole gives a very simple formula for finding the length of long lines on the spheroid (normal section arcs), given the coordinates of the end points. In the course of the computation the approximate azimuth of one end of the line is found, the error over a 500-mile line being of the order of 3″ or 4″. If the formula is amended so that the azimuth at the other end of the line is used in computing the length of the arc, the error is then less than 0″·1 over such a distance. An extra term is now given which makes this azimuth virtually correct over any distance. Numerical tests show that Cole's formula for length and the new one for azimuth are very accurate and convenient in all azimuths and latitudes.     

Vector methods to compute azimuth, elevation, ellipsoidal normal, and the Cartesian ( X , Y , Z ) to geodetic ( φ , λ , h ) transformation

Vector-based algorithms for the computation of azimuth, elevation and the ellipsoidal normal unit vector from 3D Cartesian coordinates are presented. As a by-product, the formulae for the ellipsoidal normal vector can also be used to iteratively transform rectangular Cartesian coordinates (X, Y, Z) into geodetic coordinates (φ, λ, h) for a height range from −5600 km to 10⁸ km. Comparisons with existing methods indicate that the new transformation can compete with them.     

SOME SUGGESTIONS FOR SPEEDING UP WORK ON TOPOGRAPHICAL SURVEYS

Abstract

A Few notes will now be given on the subject of triangulation on which practically all the methods already outlined depend. If we have a triangulation ready for us on which to base our work, so much the better; but, if not, we must make every effort to carry one through either from our own measured base or from any existing points on the edge of our work. For reconnaissance survey, such a triangulation must be carried out with the greatest expedition; even if all refinements are sacrificed to speed, it is extraordinary how small the errors will be found to be when a more rigid triangulation is made. Any unorthodox method such as carrying through with a resected point or with an astronomical azimuth may be adopted. A bush will often make a good point to observe to, also piles of bushes with a flag on a reed or stick.     

THE TRANSVERSE MERCATOR PROJECTION OF THE SPHEROID

Abstract

In January 1940, in a paper entitled “The Transverse Mercator Projection: A Critical Examination” (E.S.R., v, 35, 285), the late Captain G. T. McCaw obtained expressions for the co-ordinates of a point on the Transverse Mercator projection of the spheroid which appeared to cast suspicion on the results originally derived by Gauss. McCaw considered, in fact, that his expressions gave the true measures of the co-ordinates, and that the Gauss method contained some invalidity. He requested readers to report any flaw that might be discovered in his work, but apparently no such flaw had been detected at the time of his death. It can be shown, however, that the invalidities are in McCaw's methods, and there seems no reason for doubting the results derived by the Gauss method.     

A NOTE ON AZIMUTH DETERMINATION

Abstract

The purpose of this note is twofold; first, to criticize the “azimuth” section of the paper “Some Notes on Astronomy as Applied to Surveying”, by R. W. Pring (E.S.R., July 1952, xi, 85, 309–318),and secondly, out of these criticisms to develop an alternative method of making observations for azimuth. It will be apparent that this method owes much to the ideas put forward by Mr. Pring.