CELE Surveying (Geomatics) — Area and Volume Computation (Earthworks)Study Notes

The **shoelace formula** (also called the **coordinate method** or **Gauss area formula**) is the most direct approach for computing polygon areas when vertex coordinates are known. For a closed polygon with n vertices listed in order (either clockwise or counterclockwise), the formula is: **A = (1/2) |Σ(xᵢ × yᵢ₊₁ − xᵢ₊₁ × yᵢ)|** (from i = 1 to n, with vertex n+1 returning to vertex 1) Alternatively, this can be written as: **A = (1/2) |Σ xᵢ(yᵢ₊₁ − yᵢ₋₁)|** Both forms are equivalent. The absolute value ensures the area is positive regardless of the direction (clockwise vs. counterclockwise) in which vertices are ordered. **Why it works:** The method sums signed areas of triangles formed between each edge and a reference axis. The cross-product terms capture both magnitude and orientation. **Systematic approach:** 1. List vertices in order (preferably counterclockwise, which is standard in surveying). 2. Set up two columns: one for (xᵢ × yᵢ₊₁) and one for (xᵢ₊₁ × yᵢ). 3. Sum each column separately. 4. Compute the difference: Σ(xᵢyᵢ₊₁) − Σ(xᵢ₊₁yᵢ). 5. Divide by 2 and take the absolute value. **Alternative: Double Meridian Distance (DMD) method** — Uses latitudes (Δy) and departures (Δx) of each course, with the DMD accumulating the east-west offset. This is older but still appears in some exam questions. The result is equivalent to the shoelace method. **Accuracy considerations:** - The shoelace method assumes exact coordinate values (no field errors). - In practice, coordinates come from traverse closure, so closure error should be resolved first. - Rounding errors are minimal if intermediate products are carried to sufficient decimal places.

When the boundary of a parcel is irregular (curved or jagged) and cannot be easily described mathematically, surveying practice is to take perpendicular offsets at regular intervals from a base line. The area can then be estimated using **numerical integration rules**. **Setup:** - Base line (often called the "datum line") is divided into n equal segments of length h. - At each division point (station), a perpendicular offset (ordinate) yᵢ is measured. - We have n+1 ordinates: y₀, y₁, y₂, ..., yₙ. **Trapezoidal Rule** Assumes each strip between consecutive offsets is a trapezoid: **A_trap = h × [(y₀ + yₙ)/2 + y₁ + y₂ + ... + yₙ₋₁]** Or equivalently: **A_trap = h × [y₀/2 + y₁ + y₂ + ... + yₙ₋₁ + yₙ/2]** This formula sums the average of the end ordinates plus all intermediate ordinates, then multiplies by the spacing h. **Characteristics:** - Simple and always applicable (no restrictions on the number of intervals). - Approximates curved boundaries as straight lines between points (slight under-estimation for concave curves). - Less accurate than Simpson's rule but sufficient for many practical applications. - Useful as a quick check or when ordinates are unevenly spaced (though less precise). **Simpson's One-Third Rule** Assumes the boundary follows a parabolic arc between consecutive pairs of intervals: **A_Simpson = (h/3) × [(y₀ + yₙ) + 4(y₁ + y₃ + y₅ + ...) + 2(y₂ + y₄ + y₆ + ...)]** **Critical requirement:** The number of intervals must be **EVEN**, which means you must have an **ODD** number of ordinates (e.g., 3, 5, 7, 9 ordinates). **Pattern:** - First and last ordinates: coefficient 1 - Odd-indexed ordinates (y₁, y₃, y₅, ...): coefficient 4 - Even-indexed ordinates (y₂, y₄, y₆, ...): coefficient 2 **Characteristics:** - More accurate than the trapezoidal rule for curved boundaries. - Requires even number of intervals. - If you have an even number of ordinates (odd number of intervals), split off the last interval as a trapezoid and apply Simpson's to the remaining portion. - Standard method for exam problems involving irregular boundaries. **When Simpson's cannot be applied directly:** If you have 6 ordinates (5 intervals, which is odd), apply Simpson's to the first 5 intervals (4 intervals, even) and the trapezoidal rule to the last interval, then sum the results. **Common exam error:** Forgetting the even-interval requirement and applying Simpson's to an odd number of intervals—this produces an incorrect answer.

Earthwork volumes arise when computing cut-and-fill quantities for road construction, dam design, or site grading. The volume between two cross-sections (perpendicular slices of the ground) separated by distance L is calculated using the **end-area method**. **End-Area (or Average-End-Area) Formula:** **V_end = (L/2) × (A₁ + A₂)** Where: - A₁ = area of the first cross-section (m²) - A₂ = area of the second cross-section (m²) - L = distance between the two sections (m) - V_end = volume (m³) **How it works:** The formula assumes the volume between two sections is approximated by the average of the two end areas multiplied by the distance. Geometrically, this is equivalent to a prism whose base is the average of the two cross-sectional areas. **Characteristics:** - Simple and quick to calculate. - Works for both "cut" (excavation) and "fill" (placement) sections. - Slightly **over-estimates** the volume for non-prismatic shapes (e.g., when the sides slope). - Commonly used in practice because it is fast and reasonably accurate for typical earthwork situations. - Suitable when cross-section spacing is small relative to the changes in area. **Sign convention (important for cut/fill accounting):** - **Cut sections:** Positive area (excavation) - **Fill sections:** Negative area (material placement) - Mixed sections (part cut, part fill) must be split and handled separately. **Practical note:** In real projects, cross-sections are typically taken at 20–50 m intervals along a centerline. Each pair of adjacent sections yields a volume, which are then summed to get the total earthwork quantity.

The **prismoidal method** provides a more accurate volume estimate than the end-area method, especially for sections that are not simple prisms. **Prismoidal Formula:** **V_prism = (L/6) × (A₁ + 4A_m + A₂)** Where: - A₁ = area of the first cross-section (m²) - A_m = area of the **middle (mid-section)** located at L/2 (m²) - A₂ = area of the second cross-section (m²) - L = distance between sections A₁ and A₂ (m) - V_prism = more accurate volume (m³) **Key distinction:** A_m is NOT the average of A₁ and A₂; it is the actual cross-sectional area measured at the midpoint distance L/2 between the two end sections. **Why it's more accurate:** The prismoidal formula is derived from Simpson's one-third rule applied to volume. It accounts for the variation in cross-sectional area along the centerline and is exact for prismatic or parabolically varying sections. For earthwork with sloping sides (common in road and canal design), this formula typically reduces error by 50–90% compared to the end-area method. **Prismoidal Correction:** Often, engineers have only the end-area result and wish to improve it without computing the mid-section area. The **prismoidal correction** is: **Correction = V_end − V_prism = (L/12) × (A₁ − 2A_m + A₂)** This can be rewritten as: **V_prism = V_end − (L/12) × (A₁ − 2A_m + A₂)** If A_m is known, this allows quick refinement of an end-area calculation. **When to use each method:** - **End-area:** Quick estimates, small spacing, gentle slopes, or when only two section areas are available. - **Prismoidal:** Accurate calculations, large spacing, significant slope variations, or design-stage volume estimation. Always use for exam problems that provide the mid-section area. **Accuracy comparison:** For a typical fill with linear slopes from A₁ to A₂, if the mid-section area A_m differs from (A₁ + A₂)/2, the prismoidal method is significantly more accurate. Exam problems typically expect both methods to be calculated and compared.

The **mass diagram** (also called a **haul diagram** or **cumulative volume curve**) is a graphical tool that plots cumulative cut or fill volumes along the project centerline. It is essential for: 1. **Planning material haul** — determining where cut material can be used as fill 2. **Identifying borrow and waste locations** — where additional material must be imported (borrow) or excess material disposed of (waste) 3. **Optimizing haul distance** — minimizing cost by matching cut and fill volumes over the shortest distance 4. **Cross-haul analysis** — detecting inefficient material movement (e.g., hauling cut material past a fill section) **Construction:** - **Horizontal axis:** Station (distance along the project) - **Vertical axis:** Cumulative volume (m³), starting at zero - **Ordinate values:** - Increase when volume is cut (positive) - Decrease when volume is fill (negative) - **Slope of the curve:** Represents the local cut or fill rate **Key features:** - **Peaks and valleys:** Indicate transitions from cut to fill - **Rising section:** Net cut (cut exceeds fill) - **Falling section:** Net fill (fill exceeds cut) - **Horizontal tangent:** Cut and fill are in balance locally - **Intersection with horizontal reference line:** Cumulative cut = cumulative fill **Material balance:** At the completion of a project, if total cut equals total fill, the mass diagram returns to zero (or a specified final elevation). If not: - **Positive final ordinate:** Excess cut (waste material) - **Negative final ordinate:** Deficit (borrow material required) **Haul calculation:** The area between the mass diagram curve and a horizontal reference line (called the "balance line") represents the **station-metres** or **volume-metres** of haul. Dividing by the hauling distance gives the total haul quantity in cubic metres. **Exam perspective:** While full mass diagram construction is time-intensive, exam questions typically ask to: - Identify regions of cut vs. fill from ordinate slopes - Calculate cumulative volumes at specific stations - Determine borrow/waste quantities - Interpret the diagram for haul planning

Exam success in area and volume computation depends not only on knowing the formulas but also on recognizing problem types and avoiding systematic errors. This section summarizes the most frequent mistakes and strategies to avoid them. **Pitfall 1: Simpson's Rule Interval Count** **Error:** Applying Simpson's one-third rule to an odd number of intervals. **Example:** You have 6 ordinates (5 intervals). You try to apply Simpson's directly and get a wrong answer. **Solution:** Check that you have an **even** number of intervals (i.e., odd number of ordinates: 3, 5, 7, 9, ...). If not, split the last interval using trapezoid and apply Simpson's to the remainder. **Exam strategy:** Always count the ordinates first. If the count is even, you have an odd number of intervals—apply trapezoidal or split. If the count is odd, you have an even number of intervals—Simpson's is safe. **Pitfall 2: Confusing A_m in the Prismoidal Formula** **Error:** Using A_m = (A₁ + A₂)/2 (the average area) instead of the actual mid-section area. **Example:** Problem says "the section midway between stations has area 28 m²." You incorrectly average the end areas instead of using 28 m². **Solution:** A_m is **always** the **physically measured** area at the geometric midpoint distance L/2. It is NOT calculated as an average unless the problem explicitly says so. **Exam strategy:** Highlight or underline the phrase "the mid-section area is..." in the problem. That value goes directly into the 4A_m term. **Pitfall 3: Sign Errors in Shoelace Formula** **Error:** Forgetting to take the absolute value, or incorrectly calculating the cross-product terms. **Example:** You compute Σ(x_i y_{i+1}) − Σ(x_{i+1} y_i) = −50, then report area = −25 m² (negative). **Solution:** **Always take the absolute value.** A = (1/2)|−50| = 25 m². Areas are never negative. **Exam strategy:** Set up two clear columns: one for (x_i × y_{i+1}), one for (x_{i+1} × y_i). Sum each column separately. Subtract them. Take the absolute value. Divide by 2. **Pitfall 4: Mixing Units or Forgetting to Scale** **Error:** Ordinates given in metres, spacing given in feet (or vice versa). Final area computed but not converted to the required units. **Example:** h = 10 feet, offsets in metres, area computed but not converted to m². **Solution:** **Standardize all inputs to SI units first.** Convert all distances to metres, all areas to m². The formula then automatically produces the correct unit (m² for area, m³ for volume). **Exam strategy:** At the start of every problem, write a "Unit Check" line: "All inputs in SI (metres): spacing h = ___ m, ordinates in m, etc." **Pitfall 5: Applying End-Area Method When Prismoidal Is Required** **Error:** Problem provides mid-section area A_m, but you use only end areas and report the end-area volume without the correction. **Solution:** If the problem gives or allows you to compute A_m, **always** report both end-area and prismoidal volumes. Compare them. The prismoidal is more accurate. **Exam strategy:** Read the problem statement carefully. If it mentions "the middle section," "the mid-area," or "at station halfway," prismoidal is expected. **Pitfall 6: Incorrect Vertex Ordering in Shoelace** **Error:** Listing vertices out of sequence (not in a closed loop around the polygon). **Example:** You list vertices as A, C, B, D instead of A, B, C, D. The formula produces a nonsensical result. **Solution:** **Trace the polygon boundary in one direction (clockwise or counterclockwise) and list vertices in that order.** Close the loop by returning to the starting vertex. **Exam strategy:** Before applying the formula, sketch the polygon (even roughly) and number the vertices in the order you will use them. Verify that you form a closed path. **Pitfall 7: Forgetting the L/2 or L/6 Factor** **Error:** Computing V = A₁ + A₂ (omitting the L/2 factor in end-area) or V = A₁ + 4A_m + A₂ (omitting L/6 in prismoidal). **Solution:** **Include the spacing factor every time.** Write out the formula explicitly before substituting numbers. **Exam strategy:** Use a template: V_end = (L/2)(A₁ + A₂) = (___/2)(___+___) = ___ V_prism = (L/6)(A₁ + 4A_m + A₂) = (___/6)(___+4(___)+___) = ___ **Pitfall 8: Misidentifying Cut vs. Fill in Mass Diagrams** **Error:** Not recognizing that the mass diagram should rise (positive slope) for cut and fall (negative slope) for fill. **Solution:** Remember: **Cut → cumulative volume increases (curve rises)**; **Fill → cumulative volume decreases (curve falls)**. **Exam strategy:** Write "Cut = +, Fill = −" at the top of your mass diagram work. Verify each cumulative value makes sense relative to the previous one. **Strategic Problem-Solving Checklist:** 1. **Read carefully:** Identify the method required (shoelace, Simpson's, end-area, prismoidal, mass diagram). 2. **Check inputs:** Verify all units are consistent (SI), and count ordinates or sections. 3. **Set up explicitly:** Write the formula in words and symbols before substituting numbers. 4. **Compute step-by-step:** Show intermediate results. Use a calculator, but write down partial sums. 5. **Verify reasonableness:** Does the answer make physical sense? Is the area positive? Is the volume positive and non-zero? 6. **State the result:** Include units (m², m³, station-metres, etc.) and a brief interpretation if applicable.