Area Measurement: Planimeters & Green's Theorem

Introduction

    Determining the areas of irregular objects is a common and important problem. How large is the lake? How large is the drainage basin? What is the surface area of the leaf? What is the area under the peak obtained from a gas chromatograph? How much sealer do I need for the deck? How much sod do I need for my yard?
    Perhaps the most basic examples come from surveying: How do you determine the size of a parcel of land? Many fields are bordered by creeks, rivers, and lakes, thus  irregular boundaries occur frequently. The best practical information on measuring area is not found in math books, but in surveying books.
    If you had a map of a lake and wanted to determine its surface area, what would you do? You could xerox the map onto a piece of graph paper and count squares. You could xerox the map and cut out the lake and weigh the paper! What else? If you had to do this type of thing very often, you would probably invest in an electronic or mechanical planimeter. What are these devices, and how do they work?    

Electronic planimetry

    Consider again measuring the area of a lake, given a good map. One could be a little more sophisticated by dividing the lake into geometric objects of known area, other than just squares. For example, one might divide the lake into a combination of squares, rectangles, trapezoids, triangles, etc., and add up the area of all these. Generalizing a little bit more, the boundary can be approximated by an arbitrary number of straight line segments. These form a polygon.
    The area of an arbitrary polygon can be found exactly by dividing it into trapezoids.

[Graphics:Images/area_gr_1.gif]

The area of polygon ABCDEF is the area between ABCD and the x axis minus the area between DEFA and the x axis. Note that the area between ABCD and the x axis is the sum of the three trapezoids ABB'A', BCC'B', and CDD'C'. Similarly, the area between DEFA and the x axis is the sum of three trapezoids. The area of a trapezoid is the average height times the base, thus the area of ABB'A' is (AA'+BB')/2 times B'A'. If we use subscript one for the x and y coordinates of point A, and subscript 2 for the coordinates of point B, then the area of this trapezoid in terms of its coordinates is

[Graphics:Images/area_gr_2.gif]

which is positive. If we start and end at the same point (say A), and go around the polygon in a counter-clockwise direction, the signed area of the polygons underneath the closed curve will automatically be negative. For example, the polygon FAA'F' will have signed area determined by

[Graphics:Images/area_gr_3.gif]

in this example. Hence if we add up the trapezoids all the way around and make the ending point the same as our starting point we will get the area of the enclosed polygon. You might think the area is

[Graphics:Images/area_gr_4.gif]
[Graphics:Images/area_gr_5.gif]

but this isn't quite right because we need to set

[Graphics:Images/area_gr_6.gif]

then we get the correct result:

[Graphics:Images/area_gr_7.gif]
[Graphics:Images/area_gr_8.gif]

This simplifies to

[Graphics:Images/area_gr_9.gif]
[Graphics:Images/area_gr_10.gif]

Electronic planimeters calculate area using methods based on this formula. Edward Donley in Finding Areas with the Gauss-Green Formula gives an excellent example using Mathematica.

In summary, the area of a polygon with n vertices can be found from

[Graphics:Images/area_gr_11.gif]

but we have to be careful to "close the path" by replacing[Graphics:Images/area_gr_12.gif] . Also note that we need a simple closed curve which is a curve that doesn't intersect itself. Otherwise parts of the curve might not be traversed in a counter-clockwise direction which would give the wrong signed area. As an exercise these formulas should be compared to those for numerical integration using the Trapezoidal Rule. How and why do they differ?
    A planimeter is an instrument for determining the area of a plane figure. Pictures and descriptions of electronic planimeters are available from manufacturers such as Lasico . See Appendix A for a Mathematica program for the coordinate area of a polygon, and Appendix B for a description of how to use Mathematica itself as a digitizer (your own electronic planimeter!).

Green's Theorem

    Green's Theorem is in essence a form of the Fundamental Theorem of Calculus. It relates the double integral over a closed region to a line integral over its boundary:    

[Graphics:Images/area_gr_13.gif]

Applications include converting line integrals to double integrals or vice versa, and calculating areas. Green first published the theorem in 1828, but it did not become well known until 1846 when it was republished by Lord Kelvin.
    In the last section we derived the coordinate formula for the area of a polygon using only geometry:

[Graphics:Images/area_gr_14.gif]

Is it possible to derive this formula from Green's Theorem? The objective is to pick Q and P such that

[Graphics:Images/area_gr_15.gif]

then the integral on the right is the area A which may be obtained from the line integral on the left. One possible choice for P and Q is

[Graphics:Images/area_gr_16.gif]

First just consider the line integral from [Graphics:Images/area_gr_17.gif] to [Graphics:Images/area_gr_18.gif] . This can be represented parametrically by

[Graphics:Images/area_gr_19.gif]

as t goes from 0 to 1. Then

[Graphics:Images/area_gr_20.gif]
[Graphics:Images/area_gr_21.gif]

from which it follows that

[Graphics:Images/area_gr_22.gif]
[Graphics:Images/area_gr_23.gif]

Our closed path is the sum of all the straight line segments making up the polygon or

[Graphics:Images/area_gr_24.gif]

and we have derived the coordinate formula for the area of a polygon using Green's Theorem.

Mechanical planimeters

     These ingenious mechanical devices, invented by the Swiss mathematician Jacob Amsler in 1854, are a source of fascination to many. Although there isn't a great deal of readily available information on mechanical planimeters in books or journals, there is on the web. For general information on mechanical planimeters see Larry's Planimeter Planet by Larry Leinweber and planimeters by Robert Foote.
    An excellent geometric explanation of these devices has been given by Tanya Leise .
    These devices can also be explained using Green's Theorem, although it is a bit more difficult than it was for the electronic planimeter (see Tanya Leise's planimeter site ). A simpler explanation using Green's Theorem, but taking a slightly different approach is in the forthcoming calculus book by Brian Blank and Steven G. Krantz. Another reference containing both approaches (using Green's Theorem) is Calculus of One & Several Variable by Robert T. Seeley, Foresman & Co., 1973.

Conclusion

    One of the fascinating things about mathematics is seeing and finding relationships between things that at first glance may appear unrelated. We have considered electronic planimeters in detail where we found they could be easily understood from the geometric formula for the area of a polygon. The interesting part was that they could also be explained using Green's Theorem from multi-variable calculus. A similar situation exists for mechanical planimeters - they can be explained geometrically or using Green's theorem. Please pursue the links and references to see for yourself!
    Note that Amsler invented the polar planimeter in 1854 and Green's Theorem became widely known in 1846. Did Amsler know about Green's Theorem? If so, what role did it play in his invention of the planimeter?    

Appendix A: Mathematica program for coordinate area of a polygon

Appendix B: Mathematica as a Digitizer

Bruce Atwood's Math page .


Converted by Mathematica      May 19, 2000