Next Article in Journal
High Sensitivity Refractometer Based on Reflective Smf-Small Diameter No Core Fiber Structure
Previous Article in Journal
Root System Water Consumption Pattern Identification on Time Series Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Disciplinary Approach to Remote Sensing through Low-Cost UAVs

1
Departamento de Ciencias de la Computación e Inteligencia Artificial, Universidad del País Vasco UPV/EHU, 20018 Donostia-San Sebastián, Spain
2
Centro Universitario de los Valles, Carretera Guadalajara - Ameca Km. 45.5, CP 46600 Ameca, Jalisco, México
3
Centro de Investigación en Matemáticas, Jalisco SN, Col. Valenciana, CP 36240, Guanajuato, México
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1411; https://doi.org/10.3390/s17061411
Submission received: 31 March 2017 / Revised: 6 June 2017 / Accepted: 9 June 2017 / Published: 16 June 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
The use of Unmanned Aerial Vehicles (UAVs) based on remote sensing has generated low cost monitoring, since the data can be acquired quickly and easily. This paper reports the experience related to agave crop analysis with a low cost UAV. The data were processed by traditional photogrammetric flow and data extraction techniques were applied to extract new layers and separate the agave plants from weeds and other elements of the environment. Our proposal combines elements of photogrammetry, computer vision, data mining, geomatics and computer science. This fusion leads to very interesting results in agave control. This paper aims to demonstrate the potential of UAV monitoring in agave crops and the importance of information processing with reliable data flow.

1. Introduction

Remote sensing (RS) through Unmanned Aerial Vehicles (UAVs), is not only a new topic of research in the civil field, but also an alternative to conventional platforms, for the acquisition of data with infinite possibilities. Some examples of UAVs are: vegetation monitoring [1], forest inspection, mapping of territorial coverage [2], disaster response, construction monitoring [3], industrial and residential inspection, three-dimensional photogrammetric models, hydrocarbon pipeline monitoring and coastal surveillance [4]. Regarding the advantages of UAV, we also can highlight the importance to facilitate activities that have a detrimental effect on humans. Currently, we can remotely execute risk tasks, such as flying over contaminated areas, or places with high levels of radiation or in danger of collapse.
The support provided to agriculture through UAVs can be used to create alternatives with greater versatility and low cost. UAV technology in conjunction with other disciplines and fields of research are generating new applications in agriculture, such as crop identification, monitoring and mapping of cultivated areas, pest detection, crop yield estimation and prediction of anomalies. The check schedule in the field is more adjustable for the user. Therefore, monitoring crops through UAV can be a good tool for decision-making, management and planning of public policies in the agriculture. As the satellite sensors, UAV also allows for obtaining reliable data but in a more economical way.
Another important fact that has a direct impact on agriculture monitoring is climate change. This fact generates the need for low-cost and multitemporal monitoring. The increase of CO 2 leads to a remarkable change in the growth and maturation of vegetation [5], which causes problems in the crops. Through a UAV remote sensor, it is possible to detect the location of the weed within an agricultural field, and, consequently, the right herbicide in the correct amount can be applied in a specific area. UAV sensors with high spatial resolution also generate data, which allows for discriminating different types of weeds, crops and textures.
Satellite imagery and aerial photography play an important role in agricultural monitoring; they are a robust research tool for monitoring and evaluating large scale crop inventory; however, there is a limitation: free sources of satellite imagery do not provide images with enough spatial resolution [6] as those given by the UAVs.
Precision farming requires frequent information on crop conditions. UAV is an adequate tool to systematically provide information with a high resolution and low-cost in real time. However, an important fact to be taken into account is that the farmer is often not very familiar with the use of images [6]. Another aspect to consider is that the flight session of the UAVs is limited by the battery charge, which leads to short periods of use and limited areas. In this respect, it is worth mentioning their vulnerability to cyber attacks and being knocked down, and their high risk of mismanagement by the users [7]. Nowadays, there is an increasing use of UAVs in the monitoring area, and the market for ultra-light cameras with spectral and hyperspectral ranges [8,9] has grown. The use of technologies such as special sensors [10,11], synthetic aperture radar and thermal sensors in UAVs is already remarkable. The growth of digital image processing tools is also observed (for example, Skycatch [12], dronedeploy [13], among others). A report about UAV sales estimated about 13.22 billion dollars in 2016 year and speculated an increase of up to 28.27 billion dollars by 2022 [14].
The relationship between computer vision, data mining [15], geomatics, computer science and photogrammetry are generating a very interesting multidisciplinary environment in the extraction of information and its treatment by the means of UAVs [16].
Remote sensing technology allows a feasible solution for crop monitoring. This technology together with computer science leads to monitoring, detecting and classifying different types of crops. In Ref. [17], six different supervised classifiers are implemented with the purpose of monitoring and classifying corn crops through RapidEye images. The studied supervised techniques were: Maximum Likelihood [18,19], Mahalanobis and Minimum distances [20], Spectral Angle Mapper [21,22], the Parallelepiped method [23] and SVM (support vector machine) [24,25,26]. Winner-takes-all (WTA) is the final classifier [27], which is an ensemble based classifier, and, for the research in [17], WTA considers the voting derived from the six mentioned techniques. In [28,29,30], an artificial neural network (ANN) is used to predict vegetation parameters and crop yields using data from satellite images. Vegetation indices are widely used to enhance the vegetation information in satellite images with the purpose of monitoring [31,32]. The mentioned works face the precision problem during the monitoring, detection and classification process of the crops. The problems are not only due to selected classifiers, but also due to the spatial resolution of the satellite sensor, and, even though this parameter is improving constantly because of the new technology, the access to high spatial resolution data also is high-cost. Therefore, achieving a high resolution is a problem for satellite images users in the agriculture. On the contrary, through UAVs, it is possible to achieve high spatial resolution at a low-cost in comparison with satellite sensors.
Agriculture is a very important item for the economy of any country. In Mexico, agriculture plays an important role, so the monitoring of the crops using the new technology is a useful tool nowadays for the farmer’s job.
Agave is one of the important crops of Mexico. Agave is a plant with an American origin, with long leaves in a rosette shape (see Figure 1). Native to the hot and arid regions of Mexico and the Southwestern United States, the Agave Tequilana Weber plant is the one used to elaborate tequila, an original Mexican drink [33].
Ref. [34] reports the study of Agave plants through LANDSAT 7 imagery using supervised and unsupervised classification techniques. They used regression and classification tree methods and achieved 70% of Agave accuracy. The authors in Garnica et al. [34] faced the following problems: plant density in agave crops is not homogeneous, agave occurs in many soil types, humidity resulted in a problem for satellite imagery classification, and the agave is camouflaged with other covers like low tropical forest and grasslands, mainly, [34]. To mitigate the influence of all of these problems in the classification process of agave, a possible solution is the use of UAVs, enabled with high resolution cameras.
In this work, we focus on developing a solution for the monitoring of agave crops, taking advantage of the opportunity to obtain a high spatial resolution, which is provided through low-cost UAVs. For the classification, we use an unsupervised approach: k-means. During the process, we obtain ortho-mosaics, which allow us to separate plants from other elements in the cultivation. The proposed approach let us perform an inspection of agave crops; in this way, detailed monitoring that helps agave farmers in their daily work can be done. It is worth mentioning that the proposed approach could be used as well for other kinds of plants.
The rest of the paper is organized as follows: Materials and Methods are introduced in Section 2; the proposed Method is presented in Section 3; Section 4 presents the evaluation of the proposed approach, and in the final Section 5, conclusions are given and lines of future work are envisaged.

2. Materials and Methods

2.1. Work-Flow

In order to describe our research, the photogrammetric process is first presented in general terms, taking into account the state of the art. Then, a new methodology for classifying agave plants is explained. The methodology is based on photogrammetry and a k-means algorithm. Figure 2 illustrates all steps of our methodology.

2.2. Study Areas and UAV Flight Plan

This study looked at four agave areas managed by the Tequila Regulatory Council (CRT). These areas contain information of agave plants of different sizes, concentration and years of age. The area a, represented in red in Figure 3, has 3.2 ha approximately and its over flight coordinates are 20 56 4.76 N, 104 6 23.70 W. The area b, in purple, on the top right-hand side, has 5.8 ha and coordinates of 20 51 28.24 N, 103 46 48.92 W . The areas in blue and green are indicated with c, on the bottom right of the Figure 3, and they represent around 4.6 ha of land with over flight coordinates 20 45 14.24 N, 103 39 57.03 W. The immediate areas, of the previous areas, are included in this study.
The UAV flight plan was checked in advance via Google Earth (Google Inc-DigitalGlobe 2016, Mountain View, CA, USA), the take-off and landing area were specified. The flight time was about 15 min. We performed flight at different altitudes between 40 mts and 100 mts. The flight of the UAV allowed us to obtain the collection of images and their approximate coordinates. It is known that the conditions of the area, weather and the global positioning system generate errors in the image. For this reason, geodetic control points (GCPs) were distributed in the region of interest using the differential GPS.
The GCPs provide precision and at least three points are required. However, in our case, we used eight different control points for each studied area, with the aim of minimizing the error in georeferencing [35] (see Figure 4).

2.3. Description of the Sensor

For the image acquisition, we used a quad-copter Phantom 4 (DJI, Shenzhen, China), see Figure 5. It is low cost equipment, and it has an obstacle detection system of 0.7 to 15 m. Its operating environment must have good illumination to meet the objective. It has Global Positioning System-Global Navigation Satellite System, (GPS-GLONASS) mode, stabilization of 3 axes with a degree of inclination of 90 to +30 and axes of horizontal movement, vertical and rotation. The operating distance is about 3 km, and the quad-copter always flies within a clear line of sight for safety reasons. In addition to the automatic flight plan, all of the members of the team also have a manual control of the UAV and therefore the level of skills and knowledge about the manual control should be high, due to the setbacks that can arise directly in the acquisition area.
The quad-copter has a sensor sensitive to Red (R), Green (G), Blue (B) light, (RGB sensor), which allows the capture of image size of 4000 × 3000, from a height predefined by the user. The use of stabilizers allows absorbing the vibration and stabilizes the position of the sensor. The stabilizer is mounted on a gimbal platform that allows for obtaining the searched-for nadir in the images. Table 1 summarizes the main characteristics of our UAV.

2.4. Camera Calibration

We used a chessboard pattern approach for camera calibration and we obtained 16 calibration images in different orientations. For this purpose, we used the Camera Calibration Toolbox of Matlab (v. 2012, MathWorks, Inc., Natick, MA, USA) [37,38]. The results of camera calibration process are shown in Table 2.
The calibration parameters allows us to extract the information of the image. The data generated by the calibration process provide a mapping from the image to the real-world dimensions [40].
The parameter that indicates the size of the pixel is called Ground Sample Distance (GSD), and it can be calculated through the Ground Sampling Distance Calculator tool by Pix4D in.
The computation of the size is done according to the following equation:
G S D = S w · H · 100 F r · i m W ,
where GSD is the Ground Sampling Distance (centimeters/pixel) and represents the distance between two consecutive pixel centers, Sw denotes the sensor width of the camera (millimeters), H is the flight height (meters), Fr is the real focal length of the camera (millimeters) and imW is the image width (pixels). In our case, S w = 6 . 25 mm, the average of the flight height was H = 60 m, F r = 3 . 6 mm and the i m W = 4000 pixels, and, therefore, the distance between the centers of two pixels is 2.6 cm. For H [ 40 , 80 ] , the GSD [ 1 . 74 , 3 . 47 ] . The variation in altitude in the previous range did not affect the quality of the agave detection. Therefore, we suggest to use H = 60 m, in order to avoid obstacles during the flight and, in some sense, increase the time flight using the same battery.

2.5. Photogrammetric Flow

In order to obtain a good result in the image processing, a set of processing steps must be carried out [35]. Currently, in the market, there are a variety of photogrammetric software packages that can perform processes on the UAV images. These packages usually use an algorithm called structure from motion that is a set of techniques of photogrammetry and computer vision [41]. In our case, we use the software called Inpho UAS Master 6.0 (Trimble Inc, Sunnyvale, CA, USA) [42] and the application ExifTool that allows to read the metadata of a variety of photographic formats [43].
The starting point for a typical photogrammetric flow are the images set acquired during the flight. In general, all images are georeferenced [44].
In the integration process of photogrammetric flow, the most important phases are:
  • The interior orientation: it refers to the internal geometry of the camera and defines the coordinates of the principal point and focal length.
  • The exterior orientation: [45] It refers to coordinates system projection and attitude (roll, pitch and yaw), which allow for specifying, for each single image, the real position in space. These parameters may be included to Exchangeable Image File Format (EXIF-metadata) [43].
  • The aerial triangulation: it delivers 3D positions of points, measured on images, in a ground control coordinate system. This process consists in generating the correct overlap of each image [46], which, in our case, was in the horizontal of 70% and in the vertical of 30%.
We use the Root Mean Square Error (RMSE), to measure the quality of the aerial triangulation. This indicator is based on the residuals of the image coordinates and the ground coordinates. Taking into account conventional aerial photography, an RMSE of up to 1 pixel is desirable; however, according to Laliberte et al. [47], and due to larger distortion of the imagery obtained with low-cost cameras, an acceptable RMSE error is considerable of 1.5 to 2 pixels from the aerial triangulation for UAV imagery (see Table 3). Once aerial triangulation process is finished, a digital terrain model (DTM) can be generated by a dense image matching. The ortho-mosaic can be generated from UAV-based images with known camera parameters and the obtained DTM (see Figure 6). The accuracy values for DTM were: 0.08 m for the area a, 0.11 m for the area b and 0.07 m for the area c. The described procedure is automated by Inpho UAS Master, (Trimble Inc, Sunnyvale, CA,USA) in order to improve the quality of the image.
As a result of this process, we obtain a georeferenced ortho-mosaic image in GeoTIFF file format.

3. Image Processing

In our approach, we process the information corresponding to the regions located between 380 nm and 780 nm of the electromagnetic spectrum, i.e., the red (R), green (G) and blue (B) bands. The RGB ortho-mosaic is transformed into the International Commission on Illumination (Commission Internationale de l’éclairage), CIE L*a*b* color space. The CIE L*a*b* was developed by the International Commission on Illumination (CIE - Commission International de lÉlairage). CIE color spaces have the capacity to represent perceived color differences across Euclidean distance and are considered as an approximation of the human visual system [48]. For that reason, CIE color spaces are perceptually uniform. In order to convert from the RGB color space to the CIE space L*a*b*, it is first necessary to obtain the so-called artificial primaries, denoted as X, Y, Z [48]. The CIE XYZ space is the result of direct measurements on the human eye made in the late 1920s by W. David Wright [49] and John Guild [50] and serves as the basis for other color representations. The values of XYZ are calculated by means of linear transformation of the RGB given by the Expression (2):
X Y Z = 0 . 4125 0 . 3576 0 . 1804 0 . 2125 0 . 7154 0 . 0721 0 . 0193 0 . 1192 0 . 9502 R G B .
In (2), the values of R, G and B are in the interval [0,1]. The elements of the transformation matrix depend on the type of selected reference white [48,51], and these values are tabulated in [48]. We considered D65 reference white [48], which is usually used for standard RGB monitors (sRGB) [52]. The values in the space L*a*b* are calculated from the XYZ, by a non-linear transformation, see Equations (3)–(5):
L * = 116 Y Y 0 3 16 ,
a * = 500 X X 0 3 Y Y 0 3 ,
b * = 200 Y Y 0 3 Z Z 0 3 .
In Equations (3)–(5) X 0 , Y 0 and Z 0 are the values corresponding to the RGB vector [1,1,1], i.e., the white color in RGB color space. For details of the implementation, see the information described in http://www.brucelindbloom.com/.
In Figure 7, the examples of the results are shown corresponding to the color space CIE L*a*b*.
After the color space transformation, the k-means algorithm is applied on the CIE L*a*b* ortho-mosaic. This approach is a non-supervised learning algorithm, which allows for generating different class groups. k-means uses the distance criterion as a measure of similarity, and it is widely used in scientific classification schemes and in the field of pattern recognition [53]. The criterion distance justifies the use of the CIE L*a*b*.
According to the research in [34], an unsupervised approach is a feasible strategy for agave monitoring. Supervised algorithms require good samples and enough samples for the training step, and, in the case of agave study, it is very difficult to have training samples without other land covers. This being the reason, in our proposal, we use an unsupervised algorithm in order to separate the plantations of agave in relation to other land covers.
The k-means algorithm [54] allows us to create two segmented layers: agave plant and weeds.
Some authors addressed the computational limitations of k-means [55]. In order to improve the performance of the k-means, we use a parallel approach [56].
In Figure 8, an example of the results of classification through k-means is depicted. The image represented in Figure 8 a,b corresponds to regions located in study areas b and c, indicated in Figure 3.
After the classification step through k-means, we create a copy of the geographic data of our ortho-mosaic. The geographic data are extracted from the GeoTIFF file [57,58]. The created copy is annexed to the file created by means of k-means. We used the Matlab (v.2012) implementation of k-means and GeoTIFF procedures. The programs that we elaborated in Matlab allows users to fix all necessary parameters. We carried out several experiments in order to find the best number of classes, k, for the k-means algorithm. According to our results, k = 3 was the best value of k, because, with this value, the agave plants and weeds were best discriminated. The third group detected regions not relevant for our application. In our study, we use the computer workstation with a high performance processors: Intel® Xeon® (Intel Corporation, Santa Clara, CA, USA) E3-1280 v5 3.7 GHz, up to 4 GHz with Intel Turbo Boost Technology, 8 MB cache, 4 cores, with Ram memory 32 GB DDR4 and with a Serial ATA, hard drive (SATA technology, Beaverton, OR, USA) with 2 TB storage [55].

4. Evaluation of Methodology

The accuracy of the processing in the described methodology depends mainly on three aspects: on the resolution of the UAV sensor, on the photogrammetric process and on the georeference. In order to evaluate our proposal, 25 samples were taken at different sites around the area of interest: 10 of them represent weeds and 15 agave plants. Each sample was obtained with relative accuracy planimetric [59] through the georeferenced ortho-mosaic, comparing this image with segmented images of agave plants and weed on a geographic information system QGIS (Quantum Geographic Information System v2.162, Project of the Open Source Geospatial Foundation, Beaverton, OR, USA). With this procedure, we gather the information about the position of plants or weed areas into conformance with the Universal Transverse Mercator (UTM) map projection [60].
Figure 9 and Figure 10 depict how we can validate the results obtained from the segmentation of agave plants and weed.
First, we create random polygons of agave and weed. Then, we apply the identify tool in QGIS on every single polygon, and we verify the information output and the attributes for both agave and weed (see tables in Figure 9b and Figure 10b).
After checking the spatial information, the segmentation results of the selected area are verified in the field by the Tequila Regulatory Council (CRT) in Mexico.
It is worth mentioning that this process has been applied to all of the acquired images. In total, we processed four ortho-mosaics, and, for all of them, we obtained a valuable result. Table 4 contains the numerical evaluation of the segmentation of the agave plants by k-means. As it can be seen, obtained results are all over 99 . 999 % in accuracy when compared to the human made process, which has been considered by the Agave Regulation Agency as a very good result. Study areas in column 1 correspond to the areas described in Section 2.2.
Figure 11 and Figure 12 illustrate an example of the segmentation results. In both figures, (a) represents the studied land part, and (b), (c) and (d) represent the overlap between the original and segmented images.

5. Conclusions

In this work, we proposed a methodology for agave crop monitoring. The methodology combines remote sensing through low-cost UAV, photogrammetry, computer vision, data mining, geomatics and computer science. This study has demonstrated the potential development of low-cost unmanned aerial vehicles in the area of agave monitoring. We achieved excellent detection results, which is demonstrated by the obtained precision value of 99%. The monitoring of the vegetation through UAV will allow, in the near future, the generation of very important data for the study of plants such as agave. The results of this study is the base for the geospatial database, which we are building to analyze the behavior of the agave plants. At the moment, we work together with the Tequila Regulatory Council in Mexico. To the best of our knowledge, this is the first application that integrates remote sensing based on low cost UAV, image processing and pattern recognition techniques for georeferenced images for agave crop monitoring.
As future work, an extension of the presented approach is envisaged, in order to apply it to wider areas of agave and help farmers in other places different to those used in the experimental phase. The presented approach could be applied as well to supervise other types of plants; an improvement of the model is needed to this end, in order to adapt to the characteristics of the plant of interest.

Acknowledgments

We wish to acknowledge the Consejo Nacional de Ciencia y Tecnologia (CONACYT) for its financial support to the PhD studies of Gabriela Calvario. We are grateful to Cubo Geoespacial S.A .de C.V. and special to Ing. Jordan Martinez for the stimulus to this work, more information about this Company is available at: http://www.cubogeoespacial.com/. In addition, we are grateful to the support of the Tequila Regulatory Council (CRT), which has allowed us to monitor several crops. This paper has been supported by the Spanish Ministerio de Economía y Competitividad, contract TIN2015-64395-R (MINECO/FEDER, UE), as well as by the Basque Government, contract IT900-16. This work was also supported in part by CONACYT (Mexico), Grant 258033.

Author Contributions

Gabriela Calvario conceived the idea of the paper; Basilio Sierra, Teresa E. Alarcón and Oscar Dalmau helped in the experimental design and writing process; Gabriela Calvario performed the experiments; Carmen Hernández analyzed the data and the results of the experiments. All the authors have taken part in the elaboration of this document.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aasen, H.; Burkart, A.; Bolten, A.; Bareth, G. Generating 3D hyperspectral information with lightweight UAV snapshot cameras for vegetation monitoring: From camera calibration to quality assurance. ISPRS J. Photogramm. Remote Sens. 2015, 108, 245–259. [Google Scholar] [CrossRef]
  2. Horcher, A.; Visser, R.J. Unmanned aerial vehicles: Applications for natural resource management and monitoring. In Proceedings of the Council on Forest Engineering Conference, Hot Springs, AR, USA, 27–30 April 2004. [Google Scholar]
  3. Tsach, S.; Penn, D.; Levy, A. Advanced technologies and approaches for next generation UAVs. In Proceedings of the International Council of the Aeronautical Sciences, Congress ICAS, Toronto, ON, Canada, 8–13 September 2002. [Google Scholar]
  4. Gonçalves, J.; Henriques, R. UAV photogrammetry for topographic monitoring of coastal areas. ISPRS J. Photogramm. Remote Sens. 2015, 104, 101–111. [Google Scholar] [CrossRef]
  5. Kimball, B.; Idso, S. Increasing atmospheric CO2: Effects on crop yield, water use and climate. Agric. Water Manag. 1983, 7, 55–72. [Google Scholar] [CrossRef]
  6. Seelan, S.K.; Laguette, S.; Casady, G.M.; Seielstad, G.A. Remote sensing applications for precision agriculture: A learning community approach. Remote Sens. Environ. 2003, 88, 157–169. [Google Scholar] [CrossRef]
  7. Weibel, R.; Hansman, R.J. Safety considerations for operation of different classes of uavs in the nas. In Proceedings of the AIAA 3rd “Unmanned Unlimited” Technical Conference, Workshop and Exhibit, Chicago, IL, USA, 20–22 September 2004; p. 6421. [Google Scholar]
  8. Tetracam INC. Available online: http://www.tetracam.com/Products1.htm (accessed on 22 March 2017).
  9. Resonon Systems. Available online: http://www.resonon.com (accessed on 25 March 2017).
  10. Zhou, G.; Yang, J.; Li, X.; Yang, X. Advances of flash LiDAR development onboard UAV. ISPRS-Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, 39, 193–198. [Google Scholar]
  11. Geodetics Incorporated. Available online: http://geodetics.com/wp-content/uploads/2016/10/Geo-MMS.pdf (accessed on 22 March 2017).
  12. Skycatch. Available online: https://www.skycatch.com/ (accessed on 22 March 2017).
  13. Dronedeploy. Available online: https://www.dronedeploy.com/ (accessed on 22 March 2017).
  14. marketsandmarkets.com. Unmanned Aerial Vehicle (UAV) Market, by Application, Class (Mini, Micro, Nano, Tactical, MALE, HALE, UCAV), SubSystem (GCS, Data Link, Software), Energy Source, Material Type, Payload and Region—Global Forecast to 2022. Available online: http://www.marketsandmarkets.com/Market-Reports/unmanned-aerial-vehicles-uav-market-662.html (accessed on 22 December 2016).
  15. Andersson, P.J. Hazard: A Framework Towards Connecting Artificial Intelligence And Robotics. Available online: https://www.cs.auckland.ac.nz/courses/compsci777s2c/ijcai05.pdf#page=5 (accessed on 15 June 2017).
  16. Colomina, I.; de la Tecnologia, P.M. Towards A New Paradigm for High-Resolution Low-Cost Photogrammetryand Remote Sensing. In Proceedings of the International Society for Photogrammetry and Remote Sensing, (ISPRS) XXI Congress, Beijing, China, 3–11 July 2008; pp. 1201–1206. [Google Scholar]
  17. Ustuner, M.; Esetlili, M.T.; Sanli, F.B.; Abdikan, S.; Kurucu, Y. Comparison of crop classification methods for the sustanaible agriculture management. J. Environ. Prot. Ecol. 2016, 17, 648–655. [Google Scholar]
  18. Gould, W.; Pitblado, J.; Sribney, W. Maximum Likelihood Estimation With Stata; Stata Press: Station, TX, USA, 2006. [Google Scholar]
  19. Mather, P.; Tso, B. Classification Methods for Remotely Sensed Data, Second Edition; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  20. Richards, J.A. Remote Sensing Digital Image Analysis; Springer: Berlin, Germany, 2013. [Google Scholar]
  21. Kruse, F.A.; Lefkoff, A.B.; Boardman, J.W.; Heidebrecht, K.B.; Shapiro, A.T.; Barloon, P.J.; Goetz, A.F.H. The Spectral Image Processing System (SIPS)—Interactive Visualization and Analysis of Imaging Spectrometer Data. Remote Sens. Environ. 1993, 44, 145–163. [Google Scholar] [CrossRef]
  22. Kruse, F.A.; Lefkoff, A.B.; Boardman, J.W.; Heidebrecht, K.B.; Shapiro, A.T.; Barloon, P.J.; Goetz, A.F.H. Comparison of Spectral Angle Mapper and Artificial Neural Network Classifiers Combined with Landsat TM Imagery Analysis for Obtaining Burnt Area Mapping. Sensors 2010, 10, 1967–1985. [Google Scholar]
  23. Mather, P.M.; Koch, M. Computer Processing of Remotely-sensed Images: An Introduction, 4th ed.; Wiley-Blackwell: Chichester, UK, 2011. [Google Scholar]
  24. Cortes, C.; Vapnik, V. Support-Vector Networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar]
  25. Huang, C.; Davis, L.S.; Townshend, J.R.G. An Assessment of Support Vector Machines for Land Cover Classification. Int. J. Remote Sens. 2002, 23, 725–749. [Google Scholar]
  26. Melgani, F.; Bruzzone, L. Classification of Hyperspectral Remote Sensing Images with Support Vector Machines. Geosci. Remote Sens. IEEE Trans. 2004, 42, 1778–1790. [Google Scholar] [CrossRef]
  27. Zhang, C.; Xie, Z. Data fusion and classifier ensemble techniques for vegetation mapping in the coastal Everglades. Geocarto Int. 2014, 29, 228–243. [Google Scholar] [CrossRef]
  28. Pinter, P.J., Jr.; Hatfield, J.L.; Schepers, J.S.; Barnes, E.M.; Moran, M.S.; Daughtry, C.S.; Upchurch, D.R. Remote sensing for crop management. Photogramm. Eng. Remote Sens. 2003, 69, 647–664. [Google Scholar] [CrossRef]
  29. Yu, K.; Li, F.; Gnyp, M.L.; Miao, Y.; Bareth, G.; Chen, X. Remotely detecting canopy nitrogen concentration and uptake of paddy rice in the Northeast China Plain. ISPRS J. Photogramm. Remote Sens. 2013, 78, 102–115. [Google Scholar] [CrossRef]
  30. Farifteh, J.; Van der Meer, F.; Atzberger, C.; Carranza, E. Quantitative analysis of salt-affected soil reflectance spectra: A comparison of two adaptive methods (PLSR and ANN). Remote Sens. Environ. 2007, 110, 59–78. [Google Scholar] [CrossRef]
  31. Motohka, T.; Nasahara, K.N.; Oguma, H.; Tsuchida, S. Applicability of green-red vegetation index for remote sensing of vegetation phenology. Remote Sens. 2010, 2, 2369–2387. [Google Scholar] [CrossRef]
  32. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef]
  33. Bautista-Justo, M.; García-Oropeza, L.; Barboza-Corona, J.; Parra-Negrete, L. El Agave Tequilana Weber y la Producción de Tequila; Red Acta Universitaria: Guanajuato, Mexico, 2000. [Google Scholar]
  34. Garnica, J.F.; Reich, R.; Zuñiga, E.T.; Bravo, C.A. Using Remote Sensing to Support different Approaches to identify Agave (Agave tequilana Weber) CROPS. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, XXXVII, 941–944. [Google Scholar]
  35. Vallet, J.; Panissod, F.; Strecha, C.; Tracol, M. Photogrammetric performance of an ultra light weight swinglet UAV. In Proceedings of the AV-g (Unmanned Aerial Vehicle in Geomatics) Conference, Zurich, Switzerland, 14–16 September 2011. [Google Scholar]
  36. DJI. Available online: https://www.dji.com/phantom-4/info (accessed on 23 February 2017).
  37. Bouguet, J.Y. Camera Calibration Toolbox for Matlab. Available online: http://www.vision.caltech.edu/bouguetj/calib_doc/ (accessed on 25 April 2017).
  38. MathWorks. Single Camera Calibration App. Available online: https://es.mathworks.com/help/vision/ug/single-camera-calibrator-app.html (accessed on 25 April 2017).
  39. Gašparović, M.; Gajski, D. Two-step camera calibration method developed for micro UAV’s. In Proceedings of the XXIII ISPRS Congress, Prague, Czech Republic, 12–19 July 2016. [Google Scholar]
  40. Bath, W.; Paxman, J. UAV localisation & control through computer vision. In Proceedings of the Australasian Conference on Robotics and Automation, Sydney, Australia, 7 December 2005. [Google Scholar]
  41. Snavely, N.; Seitz, S.M.; Szeliski, R. Modeling the world from internet photo collections. Int. J. Comput. Vis. 2008, 80, 189–210. [Google Scholar] [CrossRef]
  42. Sani, S.; Carlos, J. Vehículos aéReos no Tripulados-UAV Para la Elaboración de Cartografía Escalas Grandes Referidas al Marco de Referencia Sirgas-Ecuador. Ph.D. Thesis, Universidad de las Fuerzas Armadas ESPE, Sangolquí, Ecuador, 2015. [Google Scholar]
  43. Cohen, K. Digital still camera forensics. Small Scale Digit. Device Forensics J. 2007, 1, 1–8. [Google Scholar]
  44. Laliberte, A.S.; Winters, C.; Rango, A. UAS remote sensing missions for rangeland applications. Geocarto Int. 2011, 26, 141–156. [Google Scholar] [CrossRef]
  45. Gašparović, M.; Jurjević, L. Gimbal Influence on the Stability of Exterior Orientation Parameters of UAV Acquired Images. Sensors 2017, 17, 401. [Google Scholar] [CrossRef] [PubMed]
  46. García, J.L.L. Aerotriangulación: Cálculo y Compensación de un Bloque Fotogramétrico; Universidad Politécnica de Valencia: Valencia, Spain, 1999. [Google Scholar]
  47. Laliberte, A.S.; Herrick, J.E.; Rango, A.; Winters, C. Acquisition, orthorectification, and object-based classification of unmanned aerial vehicle (UAV) imagery for rangeland monitoring. Photogramm. Eng. Remote Sens. 2010, 76, 661–672. [Google Scholar] [CrossRef]
  48. Wyszecki, G.; Stiles, W. Color Science: Concepts and Methods, Quantitative Data and Formulae, 2nd ed.; Wiley: New York, NY, USA, 1982. [Google Scholar]
  49. Wright, W.D. A re-determination of the trichromatic coeficients of spectral colours. Trans. Opt. Soc. 1929, 30, 141. [Google Scholar] [CrossRef]
  50. Guild, J. The colorimetric propierties of the spectrum. Philos. Trans. R. Soc. London 1932, 230, 149–187. [Google Scholar] [CrossRef]
  51. Hutson, G.H. Teoría de la Television en Color, 2nd ed.; S.A. MARCOMBO: Barcelona, Spain, 1984. [Google Scholar]
  52. Alarcón, T.E.; Marroquín, J.L. Linguistic color image segmentation using a hierarchical Bayesian approach. Color Res. Appl. 2009, 34, 299–309. [Google Scholar] [CrossRef]
  53. Jain, A.K. Data clustering: 50 years beyond K-means. Pattern Recognit. Lett. 2010, 31, 651–666. [Google Scholar] [CrossRef]
  54. Martins, O.; Braz Junior, G.; Corrêa Silva, A.; Cardoso de Paiva, A.; Gattass, M. Detection of masses in digital mammograms using K-means and support vector machine. ELCVIA: Electron. Lett. Comput. Vis. Image Anal. 2009, 8, 39–50. [Google Scholar]
  55. Grama, A. Introduction to Parallel Computing; Pearson Education: Harlow, UK, 2003. [Google Scholar]
  56. Zhao, W.; Ma, H.; He, Q. Parallel k-means clustering based on mapreduce. In IEEE International Conference on Cloud Computing; Springer: Berlin/Heidelberg, Germany, 2009; pp. 674–679. [Google Scholar]
  57. Ritter, N.; Ruth, M.; Grissom, B.B.; Galang, G.; Haller, J.; Stephenson, G.; Covington, S.; Nagy, T.; Moyers, J.; Stickley, J.; et al. GeoTIFF Format Specification GeoTIFF Revision 1.0; SPOT Image Corp.: Reston, VA, USA, 28 December 2000. [Google Scholar]
  58. Mahammad, S.S.; Ramakrishnan, R. GeoTIFF-A standard image file format for GIS applications. Map India 2003, 28–31. [Google Scholar]
  59. Bonilla Romero, J.H. Método Para Generar Modelos Digitales de Terreno con Base en Datos de Escáner láser Terrestre. Ph.D. Thesis, Universidad Nacional de Colombia-Sede Bogotá, Bogotá, Colombia, 2016. [Google Scholar]
  60. Langley, R.B. The UTM grid system. GPS World 1998, 9, 46–50. [Google Scholar]
Figure 1. Agave plants and crops.
Figure 1. Agave plants and crops.
Sensors 17 01411 g001
Figure 2. Work-flow system.
Figure 2. Work-flow system.
Sensors 17 01411 g002
Figure 3. Selected study areas. Regions in red and blue depict areas a and b, respectively, regions in blue and green correspond to area c.
Figure 3. Selected study areas. Regions in red and blue depict areas a and b, respectively, regions in blue and green correspond to area c.
Sensors 17 01411 g003
Figure 4. Ground Control Points, (GCPs), generation.
Figure 4. Ground Control Points, (GCPs), generation.
Sensors 17 01411 g004
Figure 5. Low cost equipment and field work in selected study areas. (a) Agave crop area; (b) Unmanned Aerial Vehicle, (UAV); (c) Remote Controller.
Figure 5. Low cost equipment and field work in selected study areas. (a) Agave crop area; (b) Unmanned Aerial Vehicle, (UAV); (c) Remote Controller.
Sensors 17 01411 g005
Figure 6. Ortho-mosaics processed with Inpho UAS Master corresponding to: (a) area a; (b) area b. Images in panels (c,d) correspond to polygons in the study area c.
Figure 6. Ortho-mosaics processed with Inpho UAS Master corresponding to: (a) area a; (b) area b. Images in panels (c,d) correspond to polygons in the study area c.
Sensors 17 01411 g006
Figure 7. Ortho-mosaics represented in CIE L*a*b color space: (a) Ortho-mosaic corresponding to area a; (b) Ortho-mosaic corresponding to area b. Images in panels (c,d) correspond to polygons in the study area c.
Figure 7. Ortho-mosaics represented in CIE L*a*b color space: (a) Ortho-mosaic corresponding to area a; (b) Ortho-mosaic corresponding to area b. Images in panels (c,d) correspond to polygons in the study area c.
Sensors 17 01411 g007
Figure 8. Example of classification results by k-means. (a) Image sample of the study area b; (b) Image sample of the study area c; (c) Detected agave plants in the sample represented in panel (a); (d) Detected agave plants in the image represented in panel (b); (e) Detected weeds in the sample represented in panel (a); (f) Detected weeds in the image represented in panel (b).
Figure 8. Example of classification results by k-means. (a) Image sample of the study area b; (b) Image sample of the study area c; (c) Detected agave plants in the sample represented in panel (a); (d) Detected agave plants in the image represented in panel (b); (e) Detected weeds in the sample represented in panel (a); (f) Detected weeds in the image represented in panel (b).
Sensors 17 01411 g008
Figure 9. (a) example of sample weed, showed in the QGIS system (polygon in red), (b) table with the object attributes described by QGIS, (c) segmentation of the weed sample by k-means (polygon in red), (d) full study area b.
Figure 9. (a) example of sample weed, showed in the QGIS system (polygon in red), (b) table with the object attributes described by QGIS, (c) segmentation of the weed sample by k-means (polygon in red), (d) full study area b.
Sensors 17 01411 g009
Figure 10. (a) position of agave plants in the QGIS system; (b) table with the object attributes described by QGIS; (c) sample of detected plants; and (d) full study area b described in Section 2.2, plants’ segmentation, shown in Google earth.
Figure 10. (a) position of agave plants in the QGIS system; (b) table with the object attributes described by QGIS; (c) sample of detected plants; and (d) full study area b described in Section 2.2, plants’ segmentation, shown in Google earth.
Sensors 17 01411 g010
Figure 11. Detection of agave plants and weed corresponding to the study area c Described in Section 2.2: (a) Original ortho-mosaic; (b) Agave detection bottom–top; (c) Agave detection left–right; and (d) Detected weed.
Figure 11. Detection of agave plants and weed corresponding to the study area c Described in Section 2.2: (a) Original ortho-mosaic; (b) Agave detection bottom–top; (c) Agave detection left–right; and (d) Detected weed.
Sensors 17 01411 g011
Figure 12. Detection of agave plants and Weed corresponding to the study area b Described in Section 2.2: (a) Original ortho-mosaic; (b) Agave detection top–bottom; (c) Agave detection left–right; and (d) Detected weed.
Figure 12. Detection of agave plants and Weed corresponding to the study area b Described in Section 2.2: (a) Original ortho-mosaic; (b) Agave detection top–bottom; (c) Agave detection left–right; and (d) Detected weed.
Sensors 17 01411 g012
Table 1. General characteristics of Unmanned Aerial Vehicle (UAV) [36].
Table 1. General characteristics of Unmanned Aerial Vehicle (UAV) [36].
ParameterValue
Sensor RGB6.25 mm × 4.68 mm
Weight25 grams
Sensor12.4 Megapixels
LensFOV 94
Focal length20 mm (35 mm format equivalent) f/2.8 focus at
Pixel size1.5625 μ m
Measurement of image4000 × 3000
Image TypeJPEG, DNG (RAW)
Temperature 0 to 40
Table 2. Camera calibration parameter values [39].
Table 2. Camera calibration parameter values [39].
ParameterValues
Focal length(2.2495 × 10 3 , 2.2498 × 10 3 )
Principal point coordinates(2.0159 × 10 3 , 1.5088 × 10 3 )
Skew−7.2265
Lens distortion
Tangential Distortion coefficients(0.0011, 5.6749 × 10 4 )
Radial distortion coefficients(−0.0160, −0.0336)
Num. Patterns16
Table 3. Root Mean Square Error (RMSE) aerial triangulation.
Table 3. Root Mean Square Error (RMSE) aerial triangulation.
Study AreasImages Low-Cost CameraNumber of Flight LinesImage Resolution (cm)RMSE (Pixels/cm)
a146102.601.4
b14081.631.7
c266182.101.6
Table 4. Plant detection accuracy.
Table 4. Plant detection accuracy.
Study AreasPrecisionOverall Accuracy
a0.999950.99994
b0.999980.99998
c0.999610.99961
d0.999910.99998

Share and Cite

MDPI and ACS Style

Calvario, G.; Sierra, B.; Alarcón, T.E.; Hernandez, C.; Dalmau, O. A Multi-Disciplinary Approach to Remote Sensing through Low-Cost UAVs. Sensors 2017, 17, 1411. https://doi.org/10.3390/s17061411

AMA Style

Calvario G, Sierra B, Alarcón TE, Hernandez C, Dalmau O. A Multi-Disciplinary Approach to Remote Sensing through Low-Cost UAVs. Sensors. 2017; 17(6):1411. https://doi.org/10.3390/s17061411

Chicago/Turabian Style

Calvario, Gabriela, Basilio Sierra, Teresa E. Alarcón, Carmen Hernandez, and Oscar Dalmau. 2017. "A Multi-Disciplinary Approach to Remote Sensing through Low-Cost UAVs" Sensors 17, no. 6: 1411. https://doi.org/10.3390/s17061411

APA Style

Calvario, G., Sierra, B., Alarcón, T. E., Hernandez, C., & Dalmau, O. (2017). A Multi-Disciplinary Approach to Remote Sensing through Low-Cost UAVs. Sensors, 17(6), 1411. https://doi.org/10.3390/s17061411

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop