Assessing the Accuracy and Feasibility of Using Close-Range Photogrammetry to Measure Channelized Erosion with a Consumer-Grade Camera

: Water-induced channel is one of the main forms of soil erosion in cultivated fields. Channelized erosion is often measured by the volume of the channels. Traditionally, the measurements were conducted with rulers or measuring tapes. However, these traditional methods are generally time- and labor-consuming and can cause soil surface disturbance. Close-range photogrammetry with a Consumer-Grade Camera (CGC-CRP) provides an alternative way of measuring channel volume and can overcome limitations of traditional methods and provides much higher spatial resolution. However, quantitative information on the accuracy of this technique is rare. In this study, the accuracy of the CGC-CRP method under different settings were examined with an in-house experiment and validated with a field experiment. In the in-house experiment, a wood board surface with Artificial Channels (AC) of different shapes, orientations, and sizes were built. These ACs were surveyed using the CGC-CRP method with a series of settings of shooting angles


Introduction
Water erosion is one of the main forms of soil erosion in cultivated fields, causing soil loss and crop yield reduction on site and environmental concerns off-site [1][2][3][4]. Water erosion is often shown as morphological changes or resultant morphological features. At the beginning of a water erosion event, erosion can be caused by raindrop splashing and surface runoff (in the form of sheet flow). It results in relatively uniform soil loss across the slope, which is termed sheet erosion. As the runoff depth increases, it starts to concentrate and incise into soil surface to form channels. Channels firstly appear as randomly located parallel small rills. As a rill develops, runoff cuts deeper into the rill bottom and widens the channel. Big channels that can be filled by regular tillage operation are called ephemeral gullies. The locations of these channels are often defined by topography so that they usually reappear at roughly the same locations each time even they are filled by tillage every year. When channels become too big to be erased by normal tillage operations, they are termed classic gullies [5][6][7].
Many methods have been developed to assess channelized water erosion by measuring the geometry and volume of channels. Traditionally, accurate assessment of channelized erosion is based on the measurement of multiple cross sections along the channel using a ruler and a measuring tape or a micro-topographic profiler [8,9]. These two methods are simple, low-cost, and widely used [10]. For the first method, a measuring tape is used to measure the horizontal distance across the channel; at discrete locations, a ruler is used to measure depths, based on which the cross section of the channel can be established. For the second method, an ordinary micro-topographic profiler is used. It usually consists of an array of dropping pins arranged along a wooden frame at 2 mm to 5 mm intervals. The wooden frame is placed at the ground level along the cross-sectional transect, and the dropping pins are lowered to the channel bottom, and the profile is recorded using a camera [11]. The traditional method can achieve high accuracy in some cases. For example, a field assessment of erosion damage by Herweg [12] indicated that a carefully implemented survey with a ruler and a measuring tape usually has errors less than 15%. However, the volume calculation is based on measures of discrete cross sections. Errors can be high when cross section changes erratically along the channel, which is very common due to the erratic nature of surface runoff. Moreover, these methods are generally timeand labor-consuming and, therefore, not practical for surveying large areas. As direct contact is required, they may also cause disturbance to the soil surface.
Many studies have used Light Detection and Ranging (Lidar) technology, such as aerial Lidar and terrestrial Lidar, for topography measurement [13][14][15][16][17][18]. Lidar data can be used to generate a Digital Elevation Model (DEM) for the soil surface. Channel volume can be calculated from the DEM so that it overcomes the limitations of discrete cross section measurements. For the aerial Lidar, the DEMs produced usually have a horizontal resolution greater 0.5 m [19,20] and an accuracy of about 0.3 m or better [21]. This is sufficient for detecting general surface form or classic gullies, but not for measuring small erosion channels (i.e., rills), which are usually smaller than 0.3 m wide and 0.3 m deep. The aerial Lidar is also limited by weather conditions and vegetation cover. Besides, aerial Lidar scan is costly. For the terrestrial Lidar, the resolution and accuracy are much higher than aerial Lidar. However, the equipment is very expensive and heavy to carry in the field. Moreover, the setup of terrestrial Lidar is time-consuming as it requires careful calibration and leveling. Photogrammetry has also been used in water erosion studies, especially for the channelized erosion. Photogrammetry is the technique to extract 3D information of objects through photographs [22]. Compared to the traditional method, photogrammetry has the advantage of minimal soil disturbance and low labor and time requirement. Recently, photogrammetry became more popular with the Structure from Motion (SfM) technique [23]. SfM is a technique used to reconstruct three-dimensional structures from two-dimensional image sequences that can be coupled with motion signals. With the SfM, photogrammetry becomes more reliable due to the robustness and efficiency of the algorithm proposed by Lowe [24]. Photogrammetry is classified based on camera location during photography into aerial photogrammetry and close-range photogrammetry. For the aerial photogrammetry, the camera is usually mounted under an aircraft (e.g., Unmanned Aerial Vehicles (UAVs)), whereas, for the close-range photogrammetry, the camera is usually mounted on a tripod or handheld by a photographer. The aerial photogrammetry has long been used to monitor or measure soil erosion, especially ephemeral gully erosion [25][26][27][28][29]. With the advances in technology, numerous studies have applied aerial photogrammetry for assessing channelized erosion at different scales [30][31][32][33][34][35][36][37][38][39]. Recently, close-range photogrammetry has also been utilized to monitor processes of soil erosion in many studies [40][41][42][43][44][45][46][47][48][49]. Many researchers used Consumer-Grade Cameras (CGC) in their studies. To achieve a high level of accuracy, several researchers utilized different types of support tools or metal frames to maintain the stability of the camera for obtaining quality images [43,47]. Those frames added extra weight and were not convenient to carry in the field. Castillo [41] et al. reported a study that required more than one operator for taking the photogrammetric measurements. One operator was responsible for taking the photo, whereas the other operators set up the georeferences. These methods were inconvenient in practice. Therefore, Kaiser et al. [45] pointed out that, from a practical perspective, image capturing with a handheld camera should be the easiest and lightest way to reconstruct gully morphology. However, there is still large uncertainty associated with this method.
Several factors can affect the accuracy of close-range photogrammetric measurement. In general, measurement errors associated with photogrammetry can come from two sources: camera network geometry and channel geometry. The camera network geometry in photogrammetry was associated with the procedure of the photo-taking, including settings such as shooting angle (angle between the camera lens direction and the horizontal surface), image overlapping rate, number of images used, and so on. A study from Dai et al. [50] examined the error sources for photogrammetry in construction engineering applications and summarized that (1) for the shooting angle, the closer to 90° (vertical image), the better quality of the images, thereby the better accuracy, and (2) for the image overlapping rate, the higher rate, the better accuracy. These trends were also reported by Golparvar-Fard et al. [51], who conducted a quantitative comparison regarding cost, accuracy, and ease of use of close-range photogrammetry. With respect to the channel geometry, only a few studies have been conducted. Giménez et al. [33] studied the accuracy of photogrammetric measurements of gullies (gully width ranges from 0.14 m to 17.6 m; gully depth ranges from 0.2 m to 2.7 m) and showed that the photogrammetric measurement was strongly influenced by the gully morphology, such as gully shape (e.g., U-channel versus V-channel). A major cause of the error is the shadow effect. Under natural conditions, the shadow effect occurs when there is shadow from sunlight or when part of the target surface is blocked. The shadow effects are usually associated with narrow and deep gullies. In Marzolff and Poesen's study, the factors affecting the accuracy of gully detection using photogrammetry were similar, including shadow effect; edge effect, resulting from the steep gully wall; and vegetation cover [35]. These previous studies revealed the major factors affecting the accuracy of close-range photogrammetry. However, they are mostly qualitative studies and each focused on one specific aspect of photogrammetry or a sub-group of factors. A systematic assessment on the accuracy of close-range photogrammetry in measuring the volume of water-induced channels is needed.
Overall, there is still a lack of quantitative information on how to use the close-range photogrammetry with a consumer-grade camera (CGC-CRP) to make accurate measurements of water-induced channels. Meanwhile, there is a lack of assessment on time and labor consumption associated with efforts to enhance accuracy with the CGC-CRP method. This information can be used to provide practical guidance for using the CGC-CRP method to measure the volume of waterinduced channels and the size of any geomorphological features at large. Therefore, in this paper, the CGC-CRP method was tested on two types of surfaces, and the objectives are listed below.
1. To estimate the accuracies of the CGC-CRP method under different settings on detecting erosion channels and similar morphological features.

2.
To assess the applicability of the CGC-CRP method under scenarios with various spatial, temporal, and resource constraints.

Study Surfaces
Two surface types were tested: artificial surface and soil surface. The artificial surface was used to exam accuracies of different settings, whereas the soil surface was used to validate the recommended settings, which has a good performance on detecting artificial surface.
The experiment was carried out in the Fredericton Research and Development Centre, New Brunswick, Canada. The artificial surface (Figure 1) was established in a storage building of the center. It was approximately 2. 4 m by 2.0 m in size. The base of this artificial surface was made from two large plywood boards, on which two-by-four (approximately 3.6 cm x 8.6 cm cross-sectional area) lumber and plywood board cuts were placed to make eight Artificial Channels (AC), mimicking the channel types observed in the field ( Table 1). The ACs were divided into two major groups: the Uchannels and V-channels. For the U-channels, width and depth were kept constant along the channels. In order to test the impact of size and shape, three types of U-channels were installed: regular (10.0 cm wide and 3.9 cm deep), deep (10.0 cm wide and 7.6 cm deep), and wide (20.1 cm wide and 3.9 cm deep) U-channel, respectively. For the V-channels, they were shaped to tape off at one end and, therefore, the width (from 0 cm to ~24 cm) and depth (from 0 cm to ~15 cm) varied along the length of the channels. The design of the V-channel provided the possibility to examine the effects of channel size. In addition, to test the influences of channel orientation on the accuracy, channels were installed in three orientations, referred to as vertical, angled, and horizontal channel.  The soil surface was located in the experimental farm. The field was located on an east-north facing slope with ~10% slope gradient. The soil at the site is classified as Fredericton Orthic Humo-Ferric Podzol, and the surface soil is mostly sandy loam in texture [52]. The field was under terracing practice without any crop residue from the previous year. Therefore, the bare soil surface was susceptible to water erosion. In May 2019, a ~3.5 m long rill ( Figure 2) was created by snowmelt water on the upper slope area. Unlike the U-channels on the artificial surface, the size varied along the rill. The depth of the rill ranges from 0.05 m to 0.23 m, while the width of the rill ranges from 0.30 m to 0.70 m. The rill became wider and shallower at the downslope position. The soil surface, which is a rill, provided an opportunity to validate the performance of the CGC-CRP method for detecting the geometry of rills.

Data Collection
For the artificial surface, the cross-sectional area was calculated based on the width and depth measurements to characterize the geometry of ACs. The parameter was first collected manually and was then obtained using the CGC-CRP method with different settings. The manual measurement was used as the references (control), and the differences between the manual method and the CGC-CRP method were used to assess the errors associated with the CGC-CRP method. Similarly, for the soil surface, the cross-sectional area was first measured manually and was then assessed using the CGC-CRP method.

Manual Measurement
For the artificial surface experiment, maximum channel width and depth were measured manually using a measuring tape at 10 locations in 10 cm interval along each channel. The accuracy of the manual measurements was at approximately 1 mm level. At each location, the maximum channel width was measured at the top of the channel along a transect perpendicular to the length of the channel, and maximum channel depth was measured from a meter stick resting on the channel bank to the bottom of a U-channel or the wedging tip of a V-channel. The cross-sectional area was calculated from the measured width and depth values. For the U-channels, the differences between the ten locations were minimal, so parameters determined at the 10 locations were treated as repeated measurements, and the average values were used as the reference values characterizing the channel as a whole. For the V-channels, these parameters varied along the channel, so they were treated as individual measurements at each location.
For the soil surface experiment, ten pairs of plastic sticks were inserted on both sides of the channel to indicate the positions of the transect ( Figure 2). The pairs of sticks were installed at an interval ranging from 0.2 m to 0.6 m along the channel. At each transect location, a straight ruler was placed at the top of the channel. The channel depth was measured from a vertical meter stick resting on the channel bottom to the straight ruler on the top. The depth measurement was recorded starting from the left side to the right side of the stick at 0.01 m interval along the transect. Each depth measurement represented the nearby one cm wide area. Therefore, the cross-sectional area (in cm 2 ) for each transect was calculated by summing all the depth measurements together.

Consumer-Grade Camera Based Photogrammetry (CGC-CRP)
For the photogrammetry method, a Nikon D7000 camera (Table 2) and an AF-S Nikkor 18-200 mm f/3.5-5.6G ED VR II lens were used in this study. For the artificial surface, in order to assess the impact of different factors on channel detection, the camera position and shooting angle need to be stable throughout the test, so the camera was attached to an aluminum post at approximately 1.92 m above the ground, and the shooting angle was adjusted with a tripod adapter ( Figure 3). The 1.92 m height was chosen because it is a height that is convenient to the operator and there is an original drilling hole at this height on the post. Before data collection, the lens was autofocused to the board surface firstly. Then, the camera was switched to the manual focusing mode so that the focusing distance will not change during one set of photo shooting. Calibration for lens interior orientation and lens distortion was performed by the PhotoModeler Scanner® software automatically using a built-in algorithm with multiple (more than 8) overlapped photos. The coordinate system and scale were defined using five plastic targets fixed on the board surface. Four of the five targets were placed at the middle point of each edge, and the last one was placed in the center of the board. The distance between the targets ranged from 0.9 m to 2.2 m. Their locations were precisely measured with a measuring tape. These targets can be automatically detected by the software, facilitating fast data processing. All photos were taken along the two long edges of the board surface. For all photos, the photographer was located ~0.2 m to ~0.5 m away from the edges of the board, and the camera was facing toward the surface. Six shooting angle treatments (90°, 80°, 70°, 60°, 45°, and 30°) were applied. With each shooting angle, the first photo was taken at the corner of the board, then the photographer moved along the long edge, and subsequent photos were taken at 10 cm interval. The same procedure was repeated along the other edge for the other set of photos. In total, 26 photos were taken to cover the entire surface. With 10 cm interval, the image overlapping rate was 95%. For testing the effect of overlapping rate, photos at 20 cm, 50 cm, 70 cm, 100 cm, and 120 cm intervals were selected to make datasets with image overlapping rate of 90%, 80%, 70%, 60%, and 50%, respectively. For comparing different shooting angles, data with image overlapping rate of 80% was presented because it is a common overlapping rate used in the close-range photogrammetry [53]. For comparing different image overlapping rates, data with shooting angle of 70° was presented because the 70° shooting angle was found to perform the best in the above-mentioned shooting angle test. Moreover, it is a shooting angle can cover a relatively large area and is easy to operate by a photographer in the field.
For the soil surface, the consumer-grade camera was handheld by the photographer without any support post. Similarly, before data collection, the lens was autofocused to the channel. Then, the camera was switched to the manual focusing mode. The coordinate system was defined with a reference device placed on the ground. The device has three 15 cm long arms, which are perpendicular to each other, to form a coordinate system on its own. At the end of each arm, a target was installed. As the locations of the three targets are fixed and accurately measured, when the device was captured in the photo, it can be used as a reference for length and height measurement. All photos were taken along the two long edges of the channel. The photographer was located ~0.3 m away from the edges, and the camera was facing toward the channel at ~70° shooting angle. The first photo was taken at the top of the channel, then the photographer moved along both sides of the channel, and subsequent photos were taken with ~70% overlap.

Data Extraction and Analysis
For the CGC-CRP method, the raw data were photo images. These photo images were imported into PhotoModeler Scanner®, which uses an algorithm based on the scale-invariant feature transform (SIFT) developed by Lowe [24] to extract point cloud with 3D coordinates information in 1 mm resolution. After the images were processed in PhotoModeler Scanner®, more than 5,000,000 3D coordinates were extracted for each raw dataset.
The 3D coordinate data derived from PhotoModeler Scanner® were imported into Surfer®, and the triangulation with linear interpolation method was used to generate grid format DEMs. The triangulation with linear interpolation method was selected because it is an exact interpolator, so it does not extrapolate beyond the elevation value of the data range, representing the real measured elevation values better [54]. All DEMs generated were in 1 mm resolution (cell size = 1 mm). After the DEMs were generated, ten cross-sectional reference transects, matching those measured manually, were extracted from each AC or water-induced channel. The area of the cross section was obtained for each transect using a tool provided in Surfer®. The accuracy of the method was assessed by comparing the value derived from the DEMs to the value derived from manual measurements.
For the U-channels of the artificial surface, the precision of the method was assessed using Standard Deviation (SD), Mean Error (SE), and Root Mean Square Error (RMSE), defined as follows, where Xi is the parameter value for the i th transect, indicates the average parameter value of all N transects, and Ref denotes the reference value from manual measurements. It should be noted that the SD value was not applicable for the V-channels or the soil surface channel examined in this study because the area of cross section varied along the channel, so that all measurements are individual measurements and there were no repeats. Therefore, no mean or SD can be calculated.

Applicability Measures
The applicability of manual method and the CGC-CRP method on detecting channels were tested by comparing equipment, time consumption and detection area (Figure 4). The equipment needed and time consumption of each test were recorded during the experiments.

Visual Observation of the DEMs
In all generated DEMs, all objects were clearly identifiable and in the correct locations ( Figures  5 and 6). For camera shooting angle settings, there was a general trend of increasing noise with the decrease of shooting angle, whereas for image overlapping rate settings, there was a general trend of increasing noise with the decrease of image overlapping rate. At the lowest shooting angle (30°) and lowest image overlapping rates (60% and 50%), some objects were significantly deformed.

U-Channels
For U-channel objects, the channel size maintained the same for each individual object so that the cross-sectional area measurements can be considered repeated measurements. For the vertically oriented U-channel (VU), manually measured cross-sectional area averaged at 39.0 cm 2 with a CV of 1.4%. The CGC-CRP-derived area with 70° shooting angle and 80% overlapping rate (Table 3) ranged from 25.5 cm 2 to 36.2 cm 2 , and averaged at 31.0 cm 2 with a SD of 3.4 cm 2 , a ME of −8.0 cm 2 and a RMSE of 8.6 cm 2 . Table 3. Average, standard deviation (SD), mean error (ME), and root mean square error (RMSE) of the cross-sectional area (cm 2 ) of selected transects derived from the CGC-CRP method for U-channels on the artificial surface. With respect to the different shooting angle setting, 70° shooting angle (Table 3) appears to have the lowest SD (except for comparing to 80° shooting angle), ME, and RMSE values. From the 70° shooting angle, increasing and decreasing of shooting angle both led to increasing SD (except for the 80° shooting angle), ME, and RMSE values. When the shooting angle is low, the angle of view to the ground surface is narrow, which increases the coverage area on the ground captured by the image but reduces its resolution. With fewer points representing the same area, the accuracy is expected to be lower. However, the situation is different for the channel walls. Because the wall is perpendicular to the ground surface, when the shooting angle is high, the angle of view towards the channel wall is narrow, resulting in a shorter length for the height of the wall on the image and, therefore, fewer points on the wall that can be captured by the photogrammetry. A 70° shooting angle appears to be the best shooting angle setting because it balances the errors for the ground surface and those for the channel walls. From 70°, reducing the shooting angle appears to reduce the accuracy of the photogrammetry more than increasing the shooting angle. When the shooting angle dropped from 45° to 30°, there were substantial drops in SD, ME, and RMSE values.
With respect to the image overlapping rate setting, SD, ME, and RMSE all increase with decreasing image overlapping rate, which is expected. In particular, when the overlapping rate dropped from 60% to 50%, there were substantial drops in SD, ME, and RMSE values (Table 3). It should be noted that ME values for all settings were negative, indicating that there was always an underestimation of the cross-sectional area. This could be explained by that it is hard to capture the corner of the channel, which was artificially made at a 90° angle. In fact, under all settings, the channels on the DEMs had more or less a trapezoidal shape, with the corners cut out.
For the Angled U-channel (AU) and the Horizontal U-channel (HU), the trends of accuracy measurements (SD, ME, and RMSE) with camera shooting angle and image overlapping rate were almost the same as those for the VU. Moreover, for both AU and HU, the ME values for all settings were negative, the same as for VU. However, substantial drop in accuracy for the shooting angle occurred from 60° to 45° for both AU and HU, and also from 70° to 80° for HU whereas the substantial drop in accuracy for the overlapping rate occurred from 60% to 45% for AU (the same as for VU) and 70% to 60% for HU (Table 3). This suggested that the camera shooting angle needs to be between 60° and 80° and the image overlapping rate needs to be ≥70% to maintain an acceptable level of accuracy. For AU, with only a few exceptions, the absolute values of SD, ME, and RMSE were greater than those for VU and HU, indicating a low accuracy with AU than with VU and HU (Table 3). This could be due to fact that the coordinate axes of the DEM were set up along the vertical and horizontal directions, and, as such, measurements along these two directions on the DEM were more direct. Measurements along diagonal directions require interpolation and, therefore, often have higher errors.
For the Deep U-channel (DU) and Wide U-channel (WU), the trend of SD measurements with camera shooting angle and image overlapping rate were similar to those for the VU. However, the trend of ME and RMSE measurements was different. For both DU and WU, these two measurements increased as the shooting angle dropped from 90° to 70°. A substantial decrease in measurements (ME and RMSE) for the shooting angle was occurred at 60° and 45° for DU and WU, respectively. Below those angles, the measurements (ME and RMSE) of both DU and WU experienced an increase again. Besides, unlike VU, these two accuracy measurements fluctuate with a decreasing image overlapping rate (Table 3). These different trends observed from the DU and WU are hard to explain. Further investigation might be required to understand the details. As the cross-sectional area of Regular U-channel (VU) was only half of that of Deep U-channel (DU) and Wide U-channel (WU), so the ratio between the absolute ME and true value was a reasonable accuracy parameter for comparing the different channel sizes. A higher ratio indicates a lower accuracy. The ratio of DU was the highest for the shooting angle that occurred from 90° to 70°, whereas that of WU was the smallest. Meanwhile, the ratio of VU was at a similar level as that of DU. This could be explained by that when the shooting angle is relatively high, the angle of view towards the ground surfaces of both DU and VU was narrow; therefore, fewer points on the ground surface could be captured by the photogrammetry, leading to low accuracy of DU and VU. In contrast, the ground surface of WU was wide, leading to the relatively high accuracy of WU. When shooting angle was from 60° to 30°, the accuracy of DU became higher than those of VU, probably because with a relatively low shooting angle, the angle of view towards the channel wall was wide; therefore, more channel wall area of DU could be captured than that of VU by the photogrammetry. Therefore, DU had a higher accuracy level than VU. With respect to the image overlapping rate, the ratio followed the order of DU > VU > WU for most cases. This could be rooted from that the default shooting angle (70°) was relatively high. Therefore, the photogrammetry was easier to capture the ground surface than the channel wall. As a result, the accuracy of WU was better than that of DU.

V-Channels
For V-channel objects, the channel size changed along the channel so that the cross-sectional area measurements varied. For the vertically oriented V-channel (VV), the manual measurements of the cross-sectional area (Table 4) ranged from 0.1 cm 2 to 113.9 cm 2 , whereas the CGC-CRP-derived value with 70° shooting angle and 80% overlapping rate ranged from 0.0 cm 2 to 114.2 cm 2 . Unlike the Uchannel, there was no systematic under-or overestimation with the V-channel (Figure 7). Moreover, the absolute values of error for V-channel were approximately 15% less than those of U-channel. This could be rooted from that the sidewalls of V-channel were angled, not straight down and perpendicular to the ground. Therefore, the errors due to cutting out corners were expected to be much lower. With respect to the different shooting angle settings, the accuracy measurements of VV showed a similar trend as those for the U-channel. The 70° shooting angle appears to have the lowest absolute ME and RMSE values for VV. For different image overlapping settings, the absolute values of ME and RMSE for VV experienced a slight decrease with overlapping rate dropping from 90% to 70%, but increased slightly again when the overlapping rate was less than 70%. However, the largest discrepancy of the RMSE was only 1.7 cm 2 ( Table 4). With respect to the different channel orientations, for the Angled V-channel (AV) and the Horizontal V-channel (HV), the trends of ME and RMSE with camera shooting angle were similar to those for the VV in most cases. However, the RMSE values for both AV and HV were mostly double or triple of those for VV (Table 4). For the overlapping rate, unlike VV, a substantial drop in accuracy was observed from 60% to 50% for AV, and for HV, the V-channels were not well-formed in the DEMs (Table 4; Figure 6). This result suggested that, with a low image overlapping rate, the orientation affected the accuracy significantly. Overall, the results also suggest that the camera shooting angle needs to be ≥70° and the image overlapping rate needs to be ≥60% to maintain a high level of accuracy.
With respect to channel size, for the VV with 70° shooting angle and 80% image overlapping rate, the ME value was close to 0 when the cross-sectional area was less than 10 cm 2 . When the crosssectional area was greater than 10 cm 2 , the ME value ranged from −3.0 cm 2 to 1.4 cm 2 (Figure 7). This trend indicated that the accuracy was maintained at a high level under this setting, and the channel size did not have a significant impact on the accuracy level. However, when the shooting angle increased or decreased from 70°, the ME values became higher when the cross-sectional area was greater than 10 cm 2 . Moreover, in most cases, the ME values increased with increasing channel size. In addition, the channel orientation had an impact on the ME values. The ME values followed the order of HV > AV > VV for all settings (Figure 7). This could be explained by the shade effect. In the experiment, the source of the light was perpendicular but angled towards the direction of HV. As a result, there were shades on HV. There were also shades on AV, but the sizes were smaller. In addition, the width of the shade varied with channel size. This also explains why the error decreased with a decreasing channel size.
Based on the results of both U-channel and V-channel derived from the artificial surface test, we suggest using a 70° shooting angle and ≥70% image overlapping rate because these settings could obtain acceptable levels of accuracy for both U-and V-channel in different orientations. Therefore, this setting was applied for the soil surface validation test below.

Soil Surface
For the soil surface, manually measured cross-sectional areas for the ten transects ranged from 376.2 cm 2 to 853.1 cm 2 , whereas the values derived from the CGC-CRP methods ranged from 385.9 cm 2 to 877.9 cm 2 (Table 5, Figure 8). There was no systematic under-or overestimation. The ME ranged from −7.8 cm 2 to 17.7 cm 2 , which can be translated to −1.6% to 3.1% of the measurement. The overall RMSE for the ten transects was 8.6 cm 2 , which could be translated to ~1.6% of the averaged cross-sectional area. This agrees with Castillo et al. [42], who reported that close-range photogrammetry and total station methods could produce cross-sectional area values with an error less than 4%. Table 5. Error and RMSE of cross-section area (cm 2 ) of selected transects derived from the CGC-CRP method for the water-induced channel on the soil surface. In comparison to the board surface test, the accuracy of the soil surface test was much higher. This could be explained by the lack of pattern contrast on the artificial surface. For the artificial surface, the channels were made from wood board and lumbers which were white in color. There was no clear pattern observed on these materials. As a result, the photogrammetric software was unable to pick out matching points in adjacent images correctly or obtain enough matching points, which led to higher errors. With the real soil surface, the pattern of soil particles and aggregates was clearly identifiable so that it was very easy for the software to find matching points. Therefore, the accuracy of channels on the real soil surface was much higher than that for the artificial surface. In other words, the accuracy measures derived from the artificial surface can be considered as the worst scenario values when using the CGC-CRP for real soil surface morphology measurement.

Applicability Analysis
A reasonably good quality consumer-grade camera is the major equipment used with the CGC-CRP method. The price varies with brand and region but is not very expensive. The one used in this study was approximately $800. The dimension of the camera was 132 × 105 × 77 mm, and the weight was ~690 g, which makes it easy to carry in a field tour. Moreover, with modern technology, camera adjustment is largely automatic, so no sophisticated photo shooting skills are required for taking a clear picture. The camera used in this study can be operated in various weather conditions. With respect to time consumption, for all settings, it took no more than 15 minutes to finish raw data collections (Table 6). However, it took much longer (~2 hours) for data processing. The detection area of the CGC-CRP method varied with the settings (Table 6). However, there was a trade-off between the detection area and accuracy. The detection area increased but accuracy decreased with the decrease of shooting angle. On the other hand, with a lower overlapping rate, the detection area increased but the accuracy decreased substantially.

Guidance for Using a Consumer-Grade Camera to Measure Water-Induced Channels
Results obtained from this study can be used to guide field measurement of water-induced channels using a consumer-grade camera. Lower shooting angle leading to lower accuracy in photogrammetry has been reported in many previous studies [45][46][47]. However, our results indicate that the optimum shooting angle ranges from 60° to 80°, not 90°, and minor changes in shooting angle will not affect the accuracy significantly. This provides much flexibility for shooting angle setting. In fact, it is difficult to achieve a 90° shooting angle or any other specific shooting angle. The range of 60° to 80° shooting angle allows the user to hold a camera in a relatively natural position. This range can be achieved by adjusting the height of the camera and the distance to the center of the area of interest. For example, if the photographer extends her or his arms and holds the camera 0.5 m away from the body, and 1.75 m above the ground, to achieve a shooting angle between 60° and 80°, the photographer needs to stand 0.8-1.5 m away from the center of the area of interest while pointing the camera lens at it. Apparently, this range of distance will change when the position of the camera changes. However, a photographer can practice to hold the camera at a relatively fixed position based on her or his height and arm length. The parameters associated with this position (0.5 m and 1.75 m in the previous example) can be measured before fieldwork and used to back calculate the range of distance between the photographer and the area of interest. The distance on the ground can be measured on site with a measuring tape. For image overlapping rate, it is also well known that the accuracy of photogrammetry increases with overlapping rate. Our study shows that enhancement in accuracy with increasing overlapping rate is limited when the overlapping rate is greater than 70%. This sets a lower limit on the shifting distance from one photo to its adjacent photo. The shifting distance is determined by camera parameters and its position. Assuming we are using a Nikon 7000 with the focal length set at 18 mm and the camera is held at 1.75 m above the ground with a 75° shooting angle, in order to obtain a 70% or greater overlapping rate, photos need to be taken at an interval less than 50 cm. In other words, along the channel to be measured, the maximum distance we can move from one photo to its adjacent one is 50 cm.
For camera setting, the focal length of the lens must be fixed in each survey. This is to ensure the proper calibration of camera distortion, which was done by the software automatically. For the software to quantify the distortion, it requires all paired photos to be taken at the same focal length. Fixed focal length sets the acceptable range of distance between the camera and the object (depending on the camera's depth of field), beyond which the picture will be blurry. This means that while taking multiple photos to make a DEM, the distance between the camera and the survey area needs to be kept the same as much as possible. However, this does mean that the focal length has to be kept the same all the time. Separate survey areas can use different focal lengths. Channel orientation also plays an essential role in controlling the measurement accuracy. For channel survey, taking photos from the upslope or downslope position while facing the channel vertically is the best. However, photos need to be taken from both sides of the channel to capture the walls better. Taking photos at a large angle or horizontally is not recommended as it may increase the measurement error. The CGC-CRP method can be used in various different weathers. However, bright sunshine is one of the weather conditions that needs to be avoided because it can create strong shadow effect. If the survey is conducted in a sunny day, shading the entire survey area might be necessary. Light rain and wind will not strongly affect the CGC-CRP survey. However, heavy rainfall and strong wind are to be avoided because most cameras are not waterproof and it is hard to hold the camera still in strong wind.

Comparison to Other Methods
For monitoring channelized erosion, the CGC-CRP method has several advantages over the traditional manual method (with rulers and measuring tapes). First, the CGC-CRP method requires much less time and labor in the field. In the soil surface survey conducted in the study, raw data collection in the field with the CGC-CRP method took less than five minutes as the operator only needed to take a series of photos, whereas for the traditional method, the fieldwork took two people one and a half hours. In this case, it was 36 times quicker with the CGC-CRP method. Saving fieldwork time is important because fieldwork is usually much more costly, time-sensitive, and demanding than office work, such as data analysis on a computer. Second, the CGC-CRP method can generate a DEM representing the whole detecting surface, whereas the traditional method only provides measurements for a limited number of cross-sectional areas. With a DEM, one can typically obtain more accurate channel volume measurements than with just cross-sectional area measurements because it needs interpolation to calculate channel volume. Moreover, with a DEM, one can extract many other valuable information (such as the sinuosity of the channel and the position of the thalweg in a channel).
The CGC-CRP method also has a few advantages over the aerial photogrammetry method on channel survey. First, the price of the major equipment used in the CGC-CRP method, a consumergrade camera, is lower than the major equipment used in the aerial photogrammetry method (e.g., an UAV). Second, aerial photogrammetry is strongly affected by field conditions. For example, it is not recommended to be used in forested area for channel survey because an UAV is difficult to fly in the woods and if it is flying above the tree tops, the canopy and tree trunks will block the view [55][56][57][58]. Third, the use of the aerial photogrammetry method is restricted by weather conditions. It is usually very difficult to operate UAVs in windy or rainy conditions. In contrast, the effects of weather on the CGC-CRP method are minor, although it also has its limits as discussed in the previous session. Last, for the CGC-CRP method, the camera location is close to the detection surface (usually <2 m), whereas for the aerial photogrammetry, the camera is much further away (usually >30 m). Consequently, the resolution of the DEM derived from the CGC-CRP is typically much higher than those derived from the aerial photogrammetry. Higher resolution often leads to higher accuracy. However, the higher the location of the camera, the larger the detection area. Therefore, there is a trade-off between the detection area and accuracy. Compared to the CGC-CRP method, aerial photogrammetry can usually measure a much larger area in one run but its accuracy is often lower.

Applications and Future Studies
Our study showed that the CGC-CRP method is an easy-to-use method for topographic survey of a small area. The accuracy of CGC-CRP can reach the centimeter level. For different forms of erosion, sheet erosion is normally at the millimeter level (erosion rate of 10 Mg ha -1 yr -1 is equivalent to 0.8 mm yr -1 with a soil bulk density of 1250 kg m -3 ), so the CGC-CRP method may not be accurate enough and the traditional CRP method with more strict setup and error control is probably needed. However, centimeter accuracy is sufficient for rill and gully erosion studies, as rills and gullies are normally wider and deeper than 1 cm. Many field surveyors carry a consumer-grade camera for their field trips so the equipment for the CGC-CRP method is readily available. This makes the CGC-CRP method well suited for any type of field survey. It can be a survey designed for monitoring erosion channels, a soil survey or a reconnaissance type of survey for any geoscience or environment-related topics.
The robustness of the CGC-CRP method also makes it particularly suited for some special conditions (e.g., limiting weather conditions or vegetation conditions making it difficult for other methods). For example, in Atlantic Canada, a critical time for gully erosion in agricultural fields is the snowmelt period in early spring. There is a short time window for measuring those gullies. However, in this time period, fields are wet and hard to walk on so that it is hard to carry heavy equipment into the field. This makes it very difficult for the use of terrestrial Lidar. Moreover, rainy and windy days are common in this period, which limits the applicability of the aerial photogrammetry and other airborne remote sensing methods. Further, as mentioned in the previous session, the CGC-CRP method is uniquely suited for channel survey in forested areas.
The discussion for the use of the CGC-CRP method is focused on measuring channelized erosion in this study. However, the method can be used to measure the shape and size of any objects. There are many other uses (e.g., archaeology and biology) for this method. The objects studied in other disciplines are diverse in shape (i.e., not necessary like a channel), therefore the results obtained from this study may not be directly applicable. Future studies are recommended to examine the effects of different settings on other types of shapes. Also, most cell phones nowadays are equipped with a very good camera. However, the cell phone camera typically has less control and is less precise than a typical consumer-grade camera. Future studies are recommended to test the limit of using a cell phone camera for a CRP survey.

Conclusions
Our study demonstrates that the accuracy of the measurement with the CGC-CRP method was strongly affected by camera shooting angle; image overlapping rate; as well as the shape, size, and orientation of the channels. As expected, the accuracy of the CGC-CRP method increased with image overlapping rate. To maintain an acceptable accuracy level, the image overlapping rate was suggested to be ≥ 70%. On the other hand, contrary to the popular belief on higher accuracy with higher camera shooting angle, we found that the optimum shooting angle is in the range of 60° to 80° for measuring channels and the accuracy will drop either with higher or lower shooting angles. With respect to different channel shapes, the accuracy of V-channel was ~15% higher than that of Uchannel. With respect to channel orientation, for the U-channel, the impact of the orientation on the accuracy was not significant when the shooting angle was relatively high, whereas for the V-channel, the vertically oriented channel had higher accuracy than channels oriented horizontally or in an angle. Last, channel size did not strongly affect accuracy when the channel was vertically orientated when the shooting angle and image overlapping rate was set in the optimum range. However, when the shooting angle or image overlapping rate was low, or when the channel was angled or horizontally orientated, the accuracy was lower with larger channel size. The test on a soil surface in the field proved that the CGC-CRP method could provide centimeter-level accuracy for morphological measurements on channels and other field objects. Therefore, in soil erosion studies, the CGC-CRP method is probably accurate enough for measuring rill and gully erosion (which normally result in above centimeter-level elevation change or difference). The CGC-CRP method has a few advantages over the traditional method and the aerial photogrammetry method, which makes it uniquely suitable under various conditions for field survey of erosion channels. Funding: Funds for this study were provided through Agriculture and Agri-Food Canada (AAFC) A-base projects "Sustainability measures to monitor and analyze the environmental impact of Canadian agriculture" and "Landscape Integrated Soil and Water Conservation (LISWC) on sloping fields under potato production in