Assessment of UAS Photogrammetry and Planet Imagery for Monitoring Water Levels around Railway Tracks

Figure 1.
Relative location of test sites along the railway track in eastern Ontario for assessing change in surface water levels. Railway layer source: Ontario Ministry of Natural Resources and Forestry—Provincial Mapping Unit.
Figure 1.
Relative location of test sites along the railway track in eastern Ontario for assessing change in surface water levels. Railway layer source: Ontario Ministry of Natural Resources and Forestry—Provincial Mapping Unit.
Figure 2.
Example UAS photographs showing the variability of water near the railway track right-of-way. (A) track crossing a wetland, June 2021; (B) track crossing a creek, October 2021.
Figure 2.
Example UAS photographs showing the variability of water near the railway track right-of-way. (A) track crossing a wetland, June 2021; (B) track crossing a creek, October 2021.

Figure 3.
(A) Photograph of the M300 with P1 camera, circular scanning radar and D-RTK2 base station. (B) Photograph of the M600P with the X5 camera and integrated D-RTK unit.
Figure 3.
(A) Photograph of the M300 with P1 camera, circular scanning radar and D-RTK2 base station. (B) Photograph of the M600P with the X5 camera and integrated D-RTK unit.

Figure 4.
Flowchart of the SfM-MVS workflow for generating the final products from the UAS RGB photographs. Final products are shown in blue. Processes are coloured based on the software used where green = Pix4D Mapper, orange = CloudCompare and red = Quick Terrain Modeler.
Figure 4.
Flowchart of the SfM-MVS workflow for generating the final products from the UAS RGB photographs. Final products are shown in blue. Processes are coloured based on the software used where green = Pix4D Mapper, orange = CloudCompare and red = Quick Terrain Modeler.

Figure 5.
PlanetScope Dove classic (PS) satellite imagery for sites 1A, 1B, and 2 used for water surface classification at the landscape level. Images were clipped to show the sites.
Figure 5.
PlanetScope Dove classic (PS) satellite imagery for sites 1A, 1B, and 2 used for water surface classification at the landscape level. Images were clipped to show the sites.

Figure 6.
Examples of aquatic vegetation that are important to consider within the water class when mapping water surface in the study area. (A) Pondweed (Potamogeton nodosus), (B) water lily (Nymphea odorata), (C) duckweed (Lemna minor), (D) water lettuce (Pistia stratiotes), (E) UAS photograph illustrating the aquatic plants in the context of open water near the track.
Figure 6.
Examples of aquatic vegetation that are important to consider within the water class when mapping water surface in the study area. (A) Pondweed (Potamogeton nodosus), (B) water lily (Nymphea odorata), (C) duckweed (Lemna minor), (D) water lettuce (Pistia stratiotes), (E) UAS photograph illustrating the aquatic plants in the context of open water near the track.

Figure 7.
(A) Example 4-band Planet Dove spectrum from July illustrating the similarity of the wetland class to the tree class; (B) Example multi-temporal spectra from the 20-layer stacked image illustrating wetland, tree, and an agricultural field; (C) 20-band multi-temporal image stack; band combination shown is R: September band 3, G: July band 2, B: June band 2.
Figure 7.
(A) Example 4-band Planet Dove spectrum from July illustrating the similarity of the wetland class to the tree class; (B) Example multi-temporal spectra from the 20-layer stacked image illustrating wetland, tree, and an agricultural field; (C) 20-band multi-temporal image stack; band combination shown is R: September band 3, G: July band 2, B: June band 2.

Figure 8.
Example of a final 3D point cloud (Site 3A—August) showing (A) RGB true color representation and (B) height above ellipsoid. The distribution along the vertical axis of the color bar shows the spread of the height values.
Figure 8.
Example of a final 3D point cloud (Site 3A—August) showing (A) RGB true color representation and (B) height above ellipsoid. The distribution along the vertical axis of the color bar shows the spread of the height values.

Figure 9.
Example orthomosaic (resampled to 3 cm pixel size) for Site 1A from 2 June 2021. Inset shows open water surface next to the culvert.
Figure 9.
Example orthomosaic (resampled to 3 cm pixel size) for Site 1A from 2 June 2021. Inset shows open water surface next to the culvert.

Figure 10.
Examples of the point density of the final 3D point clouds for (A) Site 3A—August and (B) Site 2—August.
Figure 10.
Examples of the point density of the final 3D point clouds for (A) Site 3A—August and (B) Site 2—August.

Figure 11.
Comparison of difference in position of points (distributions) representing the railway track and ties, the positions of which are considered to have negligible change throughout the sampling period June–October. The values in cm represent the absolute different between the June–September point clouds in comparison to October for each site. The two July datasets (Sites 1A and 1B) were acquired with the M600P and the X5 camera, all other datasets were acquired with the M300 and the P1 camera. The distributions along the vertical axes of the color bars illustrate the spread of the values.
Figure 11.
Comparison of difference in position of points (distributions) representing the railway track and ties, the positions of which are considered to have negligible change throughout the sampling period June–October. The values in cm represent the absolute different between the June–September point clouds in comparison to October for each site. The two July datasets (Sites 1A and 1B) were acquired with the M600P and the X5 camera, all other datasets were acquired with the M300 and the P1 camera. The distributions along the vertical axes of the color bars illustrate the spread of the values.

Figure 12.
Surface water (including aquatic vegetation) in percent cover across sites and periods based on a GEOBIA classification of the UAS orthomosaics (3 cm pixel size).
Figure 12.
Surface water (including aquatic vegetation) in percent cover across sites and periods based on a GEOBIA classification of the UAS orthomosaics (3 cm pixel size).

Figure 13.
Vertical distance between the rails crossing a culvert or bridge and the water level for each site across sampling dates. The error bars represent the uncertainty calculated as the variability in the rail and water surface heights.
Figure 13.
Vertical distance between the rails crossing a culvert or bridge and the water level for each site across sampling dates. The error bars represent the uncertainty calculated as the variability in the rail and water surface heights.

Figure 14.
(A) Changes in exposed surface water area and the delineation of the wetlands from multi-temporal satellite image stack; (B) Comparison of wetland baseline classification (wetland symbology) to the Ontario Ministry of Natural Resources and Forestry wetland database (red outline). The scale is the same for panels A and B.
Figure 14.
(A) Changes in exposed surface water area and the delineation of the wetlands from multi-temporal satellite image stack; (B) Comparison of wetland baseline classification (wetland symbology) to the Ontario Ministry of Natural Resources and Forestry wetland database (red outline). The scale is the same for panels A and B.

Figure 15.
Comparison of PS imagery from (A) July and (B) October of the railway track crossing the wetland in Site 2. In the July imagery the aquatic vegetation can be seen as bright green, obscuring the water surface below.
Figure 15.
Comparison of PS imagery from (A) July and (B) October of the railway track crossing the wetland in Site 2. In the July imagery the aquatic vegetation can be seen as bright green, obscuring the water surface below.

Table 1.
Acquisition mode, system, and base station used for data collection at railway sites.
Table 1.
Acquisition mode, system, and base station used for data collection at railway sites.
Date | Sampling Period | Sites | Acquisition Mode | System | Base Station | Altitude (m) |
---|---|---|---|---|---|---|
June 2 | 1 | 1A, 1B, 2 | Smart oblique | M300 + P1 | RS2 & Smartnet NTRIP | 100 |
July 30 | 2 | 1A, 1B | Double Grid | M600P + X5 | DRTK1 * | 50 |
August 6 | 2 | 2, 3A | Smart oblique | M300 + P1 | DRTK2, RS2 & Smartnet NTRIP | 80 |
September 9 | 3 | 1A, 1B | Smart oblique | M300 + P1 | DRTK2, RS2 & Smartnet NTRIP | 80 |
September 10 | 3 | 2, 3A | Smart oblique | M300 + P1 | DRTK2, RS2 & Smartnet NTRIP | 80 |
October 19 | 4 | 1A, 1B, 2, 3A | Smart oblique | M300 + P1 | DRTK2, RS2 & Smartnet NTRIP | 80 |
Table 2.
Confusion matrix for the September Site 1B classification. The water class includes aquatic vegetation as described above. The vegetation class does not include the aquatic species. The infrastructure class includes the manmade objects (e.g., rails, railway ties, bridges).
Table 2.
Confusion matrix for the September Site 1B classification. The water class includes aquatic vegetation as described above. The vegetation class does not include the aquatic species. The infrastructure class includes the manmade objects (e.g., rails, railway ties, bridges).
Infrastructure Reference |
Vegetation Reference |
Water Reference |
User’s Accuracy (%) |
|
---|---|---|---|---|
Infrastructure classification |
22 | 0 | 0 | 100 |
Vegetation classification |
0 | 84 | 0 | 100 |
Water classification |
0 | 5 | 33 | 86.8 |
Producer’s accuracy (%) |
100 | 94.3 | 100 | OA (%) = 96.5 |
Table 3.
Confusion matrix for the wetland baseline from the multi-temporal PlanetScope Dove image classification.
Table 3.
Confusion matrix for the wetland baseline from the multi-temporal PlanetScope Dove image classification.
Wetland Reference |
Other Reference |
User’s Accuracy (%) |
|
---|---|---|---|
Wetland classification |
52 | 1 | 98.1 |
Other classification |
23 | 81 | 77.8 |
Producer’s accuracy (%) |
69.3 | 98.7 | OA (%) = 84.7 |
Table 4.
UAS data acquisition and SfM-MVS product characteristics.
Table 4.
UAS data acquisition and SfM-MVS product characteristics.
No. Photos |
Area Orthomosaic (ha) |
Total No. Points 3D Cloud (MM) |
Avg. Point Density (pts/m2) | Point Cloud File Size (GB) | |
---|---|---|---|---|---|
Site 1A | |||||
2 June 2021 | 659 | 3.11 | 53.8 | 1664 | 1.37 |
30 July 2021 | 605 | 3.11 | 35.9 | 1140 | 1.19 |
9 September 2021 | 1017 | 3.11 | 120.0 | 3709 | 3.98 |
19 October 2021 | 1019 | 3.11 | 115.6 | 3574 | 4.74 |
Site 1B | |||||
2 June 2021 | 360 | 1.44 | 26.9 | 1786 | 0.89 |
30 July 2021 | 405 | 1.44 | 16.7 | 1106 | 0.55 |
9 September 2021 | 485 | 1.44 | 37.7 | 2501 | 1.25 |
19 October 2021 | 483 | 1.44 | 37.6 | 2495 | 1.25 |
Site 2 | |||||
2 June 2021 | 908 | 4.88 | 95.5 | 1900 | 3.92 |
6 August 2021 | 1262 | 4.88 | 119.7 | 2381 | 4.91 |
10 September 2021 | 1576 | 4.88 | 136.9 | 2722 | 5.61 |
19 October 2021 | 1499 | 4.88 | 135.1 | 2687 | 4.75 |
Site 3A | |||||
6 August 2021 | 529 | 1.79 | 33.2 | 1781 | 1.10 |
10 September 2021 | 590 | 1.79 | 45.1 | 2418 | 1.50 |
19 October 2021 | 589 | 1.79 | 50.8 | 2724 | 2.08 |
Total | 11,986 | 43.09 | 39.09 |
Table 5.
Total area (ha) of surface water extent mapped from the UAS orthomosaics.
Table 5.
Total area (ha) of surface water extent mapped from the UAS orthomosaics.
Site 1A | Site 1B | Site 2 | Site 3A | |
---|---|---|---|---|
June | 0.15 | 0.25 | 0.79 | – |
July/August | 0.09 | 0.17 | 0.58 | 0.08 |
September | 0.02 | 0.20 | 0.56 | 0.09 |
October | 0.07 | 0.17 | 0.75 | 0.11 |
Disasters Expo USA, is proud to be supported by Inergency for their next upcoming edition on March 6th & 7th 2024!
The leading event mitigating the world’s most costly disasters is returning to the Miami Beach
Inergency is celebrating 4 years as partner with the Disasters Expo USA 😍 😍