1958 Jun 433-442
1958 Jun 433-442
INTRODUCTION
electronic digital computer can be
T HE
applied to two major areas of applica-
tion important to the photogrammetric
engineer. The first is concerned with the
reduction of the raw data (the photograph)
to obtain the basic photogrammetric out-
put (the spatial location of points). The
computations associated with analytical
space resection and intersection, aero-
triangulation, transformations, and adjust-
ments would be examples of the first area
of computer application.
The second area of application concerns
the computations associated with the solu-
tion of engineering, scientific, and military
problems which involve the use of photo-
grammetric output data. The determina-
tion of highway earthwork, the determina-
tion of the change in the shape of a glacier,
and the supplying of guidance instructions
for a low level missile would be examples of a
computer utilizing photogrammetric data. C. L. MILLER
and the nature and requirements of his ap- the scan line may be obtained. By setting
plication of the DTM. In any of the sys- the differentiated form equal to zero, the
tems, the density of the points for a given location of the high and low points may be
type of terrain will depend on the accuracy obtained. By integrating the polynomial,
requirements associated with the applica- the area under the curve may be obtained.
tion. Different applications will require The following equations result.
different accuracy levels or degrees of per- dz
missive terrain approximation. In some Slope -
dy
= A + B [ (y - Yl) + (y - yz) .]
engineering problems, a series of DT M's
with progressively higher densities of + C[(y - Yl)(Y - yz) + (y - y,)(y - y,)
points for smaller areas might be advisable + (y - yz) (y - Y3)]
corresponding, for example, to the various Areay"..-y, = (y, - yz) [Zl - AYI - By,yz - CY1YZY3]
stages in the location and design of a high-
way. ()"z _ yzZ)
+ ---- 2
[A - B(y, + yz)
THE MATHEMATICAL TERRAIN MODEL
+ C(YlYZ + y,y, + y,y,) ]
The DT M, like the topographic map, (y,,3 _ yz3)
uses a sample of data to represent the con- + --:3-- [B - C(YI + yz + Y3)]
tinuous surface of the ground. The DT M
uses a sample of the infinite number of + ~~- yz') [C]
surface points, and the topographic map 4
uses a sample of the infinite number of
Interpolation between the scan lines
contour lines of the surface. Just as the
may be accomplished in a similar manner
engineer must interpolate on the topo-
by evaluating the polynomials across the
graphic map, the computer will have to
scan lines and parallel to the scan lines
interpolate with the DTM. And in both
which pass through the desired point from
cases straight line interpolation will often
points previously interpolated on the
be quite satisfactory. polynomials along the scan lines. An alter-
Since with the high-speed electronic
nate approach, and one which shows prom-
computer, much more sophisticated inter-
ise, is to evaluate the equation of a series
polation with the DTM is quite practical,
of surfaces which will fit the DTM points.
it is proposed that the actual model uti-
Considerable work is being performed on
lized by the computer be a mathematical
this subject but it has not yet reached a
model of the surface generated with the
stage far enough advanced for reporting at
data furnished by the DTM. Instead of
this ti me.
connecting each successive pair of points
The mathematical terrain model can be
with a straight line, a third degree poly-
justified if the number of DTM points
nomial will be generated by the computer
necessary to represent an area of interest
(Figure 2). A scan line will then consist of a
can be greatly reduced. A single third de-
number of continuous curves, each good
gree polynomial might, for example,
for a specified range of y val ues. The represent a given profile as accurately as
equation for the section between yz and Y3
would 50 straight line interpolations.
is of the form: (Figure 2) The speed and efficiency of the
z = z, + A (y - Yl) + B (y - Yl) (y - yz) electronic computer permits one to think
+ C(y - Yl)(Y - yz) (y - Y3) in terms of representing the surface of a
project area by thousands or literally tens
where the constants are computed by the of thousands of mathematical equations.
following equations in terms of the yz's of
the four points bracketing the range of DTM DATA INSTRUMENTATION
usefulness of the individual polynomial:
The most important contribution of
(zz - Zl) (Z3 - Zl) - A ()13 - Yl) photogrammetry to the DTM concept is
A = -- B = ::....::---:.:--..,.-"'''--~
(yz - )'1) (Y3 - Yl) (Y3 - yz) the practicality of obtaining the tre·
mendous amount of terrain data neces-
~-~-A~-~-B~-~~-~
C = ..:-c.._2-,--~-:-:-,'--'----:-':--~ sary. It would be quite impractical to ob-
(y, - y,)(y< - yz)(y, - Y3) tain the coordinates of, for example,
By differentiating the equation of the 10,000 points in a one square mile area in
polynomial, the slope at any point along any other way. Using photogrammetry,
DIGITAL TERRAIN MODEL-THEORY AND APPLICATION 437
1040r----,r--r----r--,---r---,---,----,r--r--r--,---,
960
940
'no
900
FIG. 2. This illustrates how terrain which would be reprcscnted by many points, using straight
line interpolation, may be reprcsented by fewer points using a third degree polynomial.
the data procurement problem is reduced lated work on the procurement of digital
to one of obtaining the coordinates of data from maps to operate a numerically
many points from the stereo model in a controlled milling machine for the produc-
fast, efficient, and accurate manner. tion of physical terrain models.
One approach to this problem is to pre- A second manual approach to obtaining
pare the standard topographic map to the DTkf data would be to scale the data
photogrammetric methods. The DTM data directly from the stereoplotter manuscript.
are readily obtainable from the contour Any of the systems of data along scan
map by graphically plotting the scan lines, lines could be selected. \Vi th some plotters
manually scaling and recording the data, it is possible to plot continuous pro/iles
and manually punching the results. Such across the stereo model. This is being ac-
an approach has the advantage that a complished in the M.LT. Photogrammetry
minimum additional investment in instru- Laboratory by means of a Nistri coordina-
mentation is required-the price of a sim- tometer and coordinatograph unit, operat-
ple scale. For limited and special use of the ing on either the Kelsh or Balplex plotter.
DTM system, such an approach might be These plots may then be digitized manu-
economical. However, the totally manual ally or automatically with a scanning and
approach is of course the slowest and most output system.
costly in manpower expenses. It would not The /irst real step in automation of the
be practical for extensive use of the DTM data procurement is obtained by digitizing
concept. Two men, one scaling and one re- one or more of the x, y, and z scanning mo-
cording, will take about 200 points per tions of the stereoplotter. If a rigid grid
hour on the average. Hence, our example system is being used, it would only be nec-
of 10,000 points would take 50 hours at a essary to digitize the z axis, the x and y
cost of 100 manhours plus the cost of positions being furnished by a projected or
punching and verifying the results. plotted grid. If the x increments are con-
Semi-automatic scanning and output stant, and/or plotted, digitizing of the y
instrumentation for operation with a map and z axes would be sufficient. Complete
as the source of data is quite possible and flexibility and automation isonly achieved
is the subject of current research and de- by digi tizing all three axes.
velopment at M.LT. Several different ap- A basic requirement for a scanning sys-
proaches to the problem are being investi- tem is that the scanning motions and the
gated and will be reported on at a later measurements they represent be converted
date. ]. A. Stieber, of the U. S. Naval to a form which can be digitized. This may
Training Device Center, is doing some re- be accomplished by converting the linear
438 PHOTOGRAMMETRIC ENGINEERING
FIG.3. The Nistri Coordinatometer shown on the Balplex Plotter. This unit is used to
obtain the x, Y, and z coordinates of points. rt is also L1sed on the Kelsh Plotter.
motions to equivalent shaft rotations by a static type, meaning that the operator
such means as a lead screw, rack and must stop for each point, and press a
pinion, continuous wire or chain, or odom- button to accomplish the automatic read-
eter type wheel. There are a number of out and recording. Two such systems have
commercial components for converting been designed and are being developed in
shaft rotations to digital output. A system the M.LT. Photogrammetry Laboratory,
of measuring and digitizing based on the one using electronic counters (Figure 4)
use of a diffraction grating might also be and the other, high speed relays. Although
used. A scanning, measuring, and digitiz- there are several commercial readout sys-
ing system of unique design based on send- tems available which are adaptable to
ing sound waves down a wire, has been de- photogrammetry, the M.LT. systems are
veloped by the Physics Department at being built to obtain a number of special
M.LT. for photographic data reduction experimental features for research work.
and shows potential for stereoplotter ap- The second basic type of readout sys-
plication. As a matter of convenience, it is tem which might be used would permit
desirable to have the scanning axes parallel readout "on-the-f1y" meaning that the
to the DTM data coordinate axes in order operator could scan the model with a con-
to directly read DT M coordinates. There- stant or variable but continuous drive. The
fore, either the stereo model or the scan- coordinates of points at equal y, z, or yz
ning unit axes should be capable of being increments would automatically be read
rotated about a vertical axes. out and recorded. Although the initial cost
For experimental work, the M.LT. would be greater than for a static readout
Photogrammetry Laboratory is using the system, the increased speed would eco-
previously mentioned Nistri unit as a nomically justify such a system in many
three-dimensional scanning unit for double cases. The actual data recording unit can
projection plotters (Figure 3). The unit is take many commercial forms such as a
rotated in stereo model space to align the tape punch, card punch, electric type-
axes. wri ter, magnetic tape recorder, or combi-
An odometer type scanning unit for use nation of such units.
with maps or plotters is being built. The The final degree of au tomation would be
Nistri unit may also be used for map scan- achieved by substituting an automatic
ning. scanning system for the human operator.
Two basic types of readout systems may Vlie can also expect to see radical changes
be used. The slowest but simplest would be in the stereoplotter itself as more advan-
DIGITAL TERRAIN MODEL-THEORY AND APPLICATION 439
FIG. 4. The console, in the center of the photograph, gives a visual readout of the coordinates
obtained from the Nistri Coordinatometer. The readout is by means of electronic counting tubes
which also drive the Bendix tape punch, shown in the lower right-hand corner of the photograph.
Three basic electronic computer pro- to select the highway profile or grade line
grams have been written at M.LT. for to meet specified optimization conditions
applying the DTM to highway engineering. and according to controls on gradient,
The Phase I program gives a general curvature, sight distance, design practice,
horizontal alignment solution for the prob- and similar geometric control specifica-
lem of relating any alignment to the DTM tions. AI though the mathematical formu-
data. The input data are the coordinates lation and programming of such problems
and radius of curvature at each point of will be quite complex and costly, they are
intersection-the P.L-and an origin of entirely practical and economically justi-
centerline stationing. The input data can fied.
be referenced to a coordinate system dif-
ferent from that of the DT M data and MULTIPLE VARIABLE EVALUATION
usually will be stated plane coordi na tes. The
The usefulness of the DT M approach
output of the Phase I program is the D.T M
can be greatly extended by simultaneously
coordinates and centerline stationing of
considering additional variables related to
selected points along the alignment, usu-
the problem at hand. This can be accom-
ally at stated intervals such as every 50
plished by attaching classification infor-
feet. The program also serves as a general
mation and quantitative data describing
solution to the profile problem, furnishing
other variables at each point. For example,
a terrain profile defined by the intersection
the total data associated with each terrain
of a series of plane and cylindrical vertical
point in the M.I.T. experimental work and
surfaces ,with the model.
automatically recorded by the output in-
The Phase II program gives a solution
strumentation take the form:
to the vertical alignment problem and re-
lates the profile reference line of the high- pppp xxxxxx YYYYY ZZZZZ CCNN
way surface to the DTM data. The input
to Phase II consists of the station, eleva- comprised of four identification digits, six
tion, and length of parabolic curve associ- digits of the x coordinate, five of y, five of
ates with each P.I. The output is the refer- z, a classification and two digits of quanti-
ence elevations at the same points com- tative data. Six different formatsof the out-
puted in Phase I. put data are presently possible by means of
After the reference line for the proposed a six position switch for different IBM
highway surface has been completely card, Remington Rand card, and Bendix
fixed in three dimensional space and re- tape punch data formats.
lated to the DT M data, the Phase III pro- The classification and quantitative data
gram generates the highway surface, com- associated with a given point can take
putes intersections of the surface and the many forms and includes a number of dif-
model, and determines the enclosed areas ferent types of variables. For example, the
and volumes. classification might be the type of soil at
In the present program, the surface gen- the point, and the quantitative digits the
erated or the cross-section templet can be depth of overburden to bedrock. Such data
composed of a series of 20 planes, all of might be obtained by airphoto analysis
which are variable from job to job, and and geophysical methods. The computer
six of which are variables selected by the program using such DT M data would com-
computer. Although these programs have pute the volumes of each type of material
immediate practical application to a vari- included in the construction requirements.
ety of problems, they represent only the Alternately, the quantitative data might
crude beginning of the programs which will represent the unit right-of-way cost and
ultimately be used with the DTM. the computer would be expected to deter-
The programs described above require mine the relative land acquisition cost for
that the engineer specify to a great extent each alignment evaluated. The solution of
the shape and location of his surface of in- the ul timate forms of such problems will be
terest. In many cases it will be possible based on the digital cost model previously
and desirable for the computer to assist in proposed by the senior author. In solving
selecting the location of the engineering the problems, the computer will evaluate
surface, according to limits and controls the most economical solution considering
set by the engineer. An example would be a all cost and benefit variables. Admittedly,
problem in which the computer is required such problems become quite complex and
DIGITAL TERRAIN MODEL-THEORY AND APPLICATION 441