Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm Using Multispectral Imagery Based UAV
Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm Using Multispectral Imagery Based UAV
Abstract: Oil palm is a commodity which contributes to the largest foreign exchange. In Indonesia, the area of oil palm
plantations has a large area compared to other commodities. Proper and efficient fertilization is needed to
reduce production costs. This study aims to estimate the nutrient content of calcium, magnesium, and sulfur
in oil palm plants using UAV-based multispectral cameras. The method used is divided into three stages, data
preparation, pre-processing, and data analysis. At the data preparation stage, the things done are leaf sampling,
sample coordinate points, and multispectral image capture. In the pre-data processing stage the things done
are stitching, georeferencing, and digitizing. The last stage is data analysis like multiple linear regression
analysis to get a model of the relationship of multispectral images and actual nutrition of lab test results. The
results of this study obtained a model for predicting of calcium content Ca = 0.994+0.00723*GREEN-
0.00863*RED EDGE, for magnesium content is Mg = 0.693+0.00531*RED-0.00541*RED EDGE, and for
sulfur content is S = -0.222+0.00338*RED EDGE. The results of overall accuracy using a confusion matrix
of 66.7% in the calcium model, 63.3% in the magnesium model, and 36.6% in the sulfur model.
127
Suyuthi, M., Seminar, K. and Sudradjat, .
Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm using Multispectral Imagery based UAV.
DOI: 10.5220/0009978700002833
In Proceedings of the 2nd SEAFAST International Seminar (2nd SIS 2019) - Facing Future Challenges: Sustainable Food Safety, Quality and Nutrition, pages 127-134
ISBN: 978-989-758-466-4
Copyright
c 2022 by SCITEPRESS – Science and Technology Publications, Lda. All rights reserved
2nd SIS 2019 - SEAFAST International Seminar
128
Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm using Multispectral Imagery based UAV
2.3.1 Data Preparation Stage to be used as a reference point when uniting images
into a complete map. (Adillah 2018). A tool to draw
In the data preparation stage, several things are GCP coordinates can be by using a handheld GPS.
carried out, namely literature study, method
determination, discussion, and data collection. After 2.3.2 Pre-processing Stage
the research planning is done, the data collection
stage is carried out. The data collection phase is done The objectives of the pre-processing stage such as
by three things, namely sampling, drone imagery, and normalization and noise reduction are to produce
GCP. Leaf sampling is carried out in a spread where clean and ready-to-use data (Nanda et al. 2019; Nanda
the sample plants are proportionally determined. et al. 2018a; Nanda et al. 2018b). At this stage, the
Guidelines for sampling leaves on oil palm plants stitching and georeferencing process is performed
based on Winarna et al. (2005) where the leaves used using pix4Dmapper software. The stitching process is
as the main sample must meet a number of provisions, basically a combination or combination of two or
namely they are not the mains of inserts, grow more different images to create or form one image
normally, do not lie adjacent to roads or ditches/rivers, called a panorama (Kale and Singh 2015). Before
do not coexist with insertion trees and are not attacked stitching, georeferencing needs to be done, namely
by pests or diseases. the process of giving geographic references to raster
Leaf samples taken were leaves from the 17th or images that do not yet have a coordinate system
midrib. According to Chapman and Gray (1949) in reference. The coordinate reference used is WGS 84 /
Pahan (2006) said that the leaves of the 17th midrib UTM zone 48S with EPSG: 32748.
are the most sensitive leaves because they show the Then do bordering of the garden and canopy of the
greatest difference in nutrient levels. In addition, sample plants done using QGIS software. The results
nutrient status on the 17th leaf has a better correlation of the bordering of the sample plant canopy are used
to crop production when compared to other younger to extract the reflectance value in the form of a digital
leaves. The leaves of the 17th midrib are taken by six number at each pixel. The tool used to extract these
leaflets (three strands on the left and three strands on values is Zonal Statistics on the Raster menu.
the right at the meeting point of the two sides of the
midrib). Leaves that have been taken are stored in 2.3.3 Data Analysis Stage
envelopes that have been labeled according to the
location of the sample. The selected sample plants Estimation using multiple linear regression analysis
were given raffia to indicate the tree was a sample with the stepwise method is performed by Minitab
plant. Then the sample plants are marked on GPS which aims to determine the factors that effect and
which will be used to correct the geometry between make a estimation model of calcium, magnesium, and
the map and the image results. Samples that have been sulfur nutrition. At this stage two types of variables
obtained were analyzed for nutrients of calcium, are needed, namely the independent variable and the
magnesium, and sulfur in the Testing Laboratory of dependent variable. The independent variable
the Department of Agronomy and Horticulture. contains the reflectance value data which is the
Multispectral image capture using a Phantom 4 average value of the drone image in the digitizing
drone with a Parrot Sequoia camera. Drone flight attribute of the sample plant. While the dependent
planning automatically uses the Pix4D Mapper variable contains the actual calcium, magnesium, and
application which is adjusted to the taking land area sulfur nutrition data from the 17th midrib nutrient
and hours of time that can be taken by the drone. analysis results in the sample plants.
According to Kasih (2012) the optimum shooting is The model obtained is then evaluated to see the
done in the morning because the effect of reflected strength of the estimator model. The method used to
light from the sun is still weak. Besides the wind evaluate the model is the Mean Absolute Percentage
speed in the morning still tends to be low, thereby Error (MAPE). MAPE shows how big the difference
reducing the risk of UAV shake while shooting which between the actual results and the predicted results.
can cause poor quality captured images. Table 2 shows the criteria used to decide the
GCP determination aims to reduce errors or predictive power of an estimator model.
changes in position when integrating the results of the
image into a map that can make the data change. GCP
can be in the form of objects, buildings, or forms of
certain locations that can be clearly seen on the image
so that the coordinates of the object can be measured
129
2nd SIS 2019 - SEAFAST International Seminar
Table 2: Criteria for estimating the strength of the model. plantations and sample canopy are done manually
MAPE (%) Power of prediction using the add feature tool after creating a new
<10 Very Good polygon layer in QGIS software. Figure 1 and Figure
10 – 20 Good 2 show the results of bordering the sample plant
20-50 Moderate canopy digitized to extract the average value of the
>50 Bad digital number (DN) at each pixel. The extraction
Source: Wang et al. (2012) process is carried out on all layers of stitches using
Zonal Statistics in QGIS. Extraction results can be
After that, the nutrient estimation layer is made seen in appendix 1.
which contains the estimated nutrient value for each
pixel. The estimated nutritional value is the
application of the model obtained using the raster
calculator tool in the QGIS application. Then the
estimate value is classified based on the criteria in
Table 3.
130
Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm using Multispectral Imagery based UAV
comparison of actual and predicted nutrition. optimum conditions, light blue at above optimum
Evaluate the model using Mean Absolute Percentage conditions, and dark blue at excess nutrient conditions.
Error (MAPE). Evaluation of the model is done to see Figure 3 shows the results of the calcium nutrient
how strong the estimation of the model is made. The estimation layer, Figure 4 for magnesium nutrition,
results of the model evaluation obtained the value of and Figure 6 for sulfur nutrition.
MAPE in the Ca model, Mg model, and the S model The result of nutrient estimation layer is
respectively 22.72%, 11.62%, and 46.5%. Based on performed by accuracy test using a confusion matrix.
the strength estimation criteria of the model according In the confusion matrix the predicted nutrition of the
to Wang et al (2012) in table 2, the evaluation results model results is compared with the actual nutrition of
show that the Ca and S models are categorized as the lab test results. From the results of the accuracy-
medium and the Mg model is categorized as good. test using a confusion matrix, Overall Accuracy
The model obtained is then interpreted in the form of obtained in the nutrients of calcium, magnesium, and
an estimated nutrient layer that has been processed sulfur was 66.6%, 63.3%, and 36.7%. According to
and classified according to table 3. Classification of Jaya (2015) a good accuracy value is an accuracy
nutritional assumptions is presented in the form of a value that has reached a score of> 85%, so the
color that matches the nutrient content. The red color estimation of calcium and magnesium nutrients has
indicates that the oil palm tree on the grid lacks moderate accuracy and for sulfur has poor accuracy.
nutrients, yellow under optimal conditions, green at
Table 4: Comparison of actual nutrition and prediction of the 17th midrib leaf laboratory results.
Magnesium (%DM) Kalsium (%DM) Sulfur (%DM)
Sample
Actual Prediction Actual Prediction Actual Prediction
1A 0.36 0.296 0.51 0.643 0.15 0.231
1B 0.34 0.296 0.69 0.594 0.17 0.282
1C 0.19 0.304 0.53 0.616 0.17 0.266
1D 0.32 0.316 0.82 0.668 0.13 0.193
1E 0.32 0.284 0.69 0.622 0.43 0.264
2A 0.32 0.383 0.39 0.621 0.18 0.193
2B 0.42 0.356 0.59 0.642 0.14 0.218
2C 0.32 0.378 0.69 0.649 0.13 0.184
2D 0.41 0.390 0.91 0.658 0.40 0.167
2E 0.39 0.353 0.79 0.651 0.17 0.203
2F 0.35 0.351 0.68 0.638 0.17 0.233
2G 0.34 0.379 0.67 0.683 0.07 0.156
2H 0.44 0.386 0.42 0.668 0.09 0.181
3A 0.40 0.357 0.48 0.632 0.46 0.220
3B 0.38 0.370 0.68 0.603 0.40 0.265
3C 0.36 0.424 0.82 0.671 0.16 0.125
3D 0.30 0.375 0.79 0.615 0.14 0.236
3E 0.40 0.385 0.61 0.633 0.22 0.192
3F 0.41 0.389 0.66 0.661 0.34 0.157
3G 0.47 0.355 0.62 0.647 0.23 0.202
3H 0.45 0.416 0.55 0.629 0.20 0.182
3I 0.41 0.461 0.60 0.582 0.23 0.218
3J 0.49 0.458 0.51 0.583 0.28 0.221
3K 0.44 0.410 0.70 0.570 0.19 0.286
3L 0.33 0.375 0.67 0.591 0.36 0.263
4A 0.48 0.362 0.56 0.608 0.26 0.247
4B 0.50 0.423 0.72 0.634 0.10 0.147
4C 0.37 0.387 0.87 0.625 0.17 0.169
4D 0.43 0.413 0.45 0.619 0.15 0.165
4E 0.40 0.382 0.52 0.577 0.15 0.243
131
2nd SIS 2019 - SEAFAST International Seminar
132
Estimation of Calcium, Magnesium and Sulfur Content in Oil Palm using Multispectral Imagery based UAV
133
2nd SIS 2019 - SEAFAST International Seminar
Utilization to Support the Production of Remote Sudradjat. (2016). Kelapa Sawit: Peningkatan
Sensing Information). Inderaja, 5(7), pp. 24-31. Produktivitas(Palm Oil: Increased Productivity).
Goh, K.J., Hardter, R. (2003). General Oil Palm Nutrition. Bogor, Indonesia: IPB Press, p. 68.
In: Fairhust T.H., Hardter, ed., Oil Palm: Management Von Uexkull, H.R. and Fairhurst, T.H. (1991). Fertilizing
for Large and Sustainable Yield. Switzerland: PPI., pp. for High Yield and Quality. Basel, Switzerland:
191-230. International Potash Institute.
Gerendas, J. and Heng, A. (2010). Oil Palm Fertilization – Wang, F., Hsiao, Y. and Chang, K. (2012). Combining
Sharing Some Perspectives. Plantation Industry: Diffusion and Grey Models Based on Evolutionary
Competitive Strategies In Achieving A Sustainable Optimization Algorithms to Forecast Motherboard
Future. Selangor, Malaysia, pp. 1-6. Shipments. Journal of Mathematical Problems in
Jaya, I.N.S. (2015). Analisis Citra Digital: Perspektif Engineering, pp 1-10.
Penginderaan Jauh untuk Pengelolaan Sumberdaya Winarna, Sutarta, E.S. and Sugiono. (2005). Pedoman
Alam (Digital Image Analysis: A Remote Sensing Pengambilan Contoh Daun dan Tanah pada Tanaman
Perspective for Natural Resource Management). Bogor, Kelapa Sawit (Guidelines for Sampling Leaves and
Indonesia: IPB Press. Soils in Palm Oil Plants). Medan, Indonesia: Pusat
Kale, P. and Singh, K.R. (2015). A Technical Analysis of Penelitian Kelapa Sawit.
Image Stitching Algorithm. International Journal of
Computer Science and Information Technologies, 6(1),
pp. 284-288.
Kasih, L.S.B. (2012). Penentuan Iluminansi dan Ketinggian APPENDIX
Terbang Pesawat Optimum untuk Pemetaan Tingkat
Warna Daun Padi (Determination of Optimum Aircraft Digital number values of digitization of the sample plant
Illuminance and Altitude for Mapping the Color of Rice canopy.
Leaves). Bachelor. Bogor, Indonesia: IPB University.
Sample Green NIR Red Red Edge
Nanda, M. A., Seminar, K. B., Nandika, D., & Maddu, A.
(2019). Development of Termite Detection System 1A 97.696 150.290 55.234 126.499
Based on Acoustic and Temperature Signals. 1B 99.226 162.430 60.269 128.360
Measurement, 147, 106902. 1C 107.011 158.916 62.020 128.823
Nanda, M. A., Seminar, K. B., Nandika, D., & Maddu, A. 1D 100.055 149.172 54.509 120.953
(2018a). A Comparison Study of Kernel Functions in 1E 101.391 162.627 55.288 125.934
the Support Vector Machine and Its Application for 2A 92.608 131.776 71.301 127.583
Termite Detection. Information, 9(1), 5. 2B 110.319 123.704 75.585 139.960
Nanda, M. A., Seminar, K. B., Nandika, D., & Maddu, A. 2C 103.668 124.700 73.264 132.114
(2018b). Discriminant Analysis as a Tool for Detecting 2D 104.833 123.975 73.732 129.730
the Acoustic Signals of Termites Coptotermes 2E 107.398 130.823 70.143 133.121
curvignathus (Isoptera: Rhinotermitidae). International 2F 118.497 135.772 75.499 137.183
Journal of Technology, 9(5): pp. 840-851. 2G 112.261 130.371 69.891 125.834
Pahan, I. (2006). Panduan Lengkap Kelapa Sawit 2H 122.389 129.942 77.669 131.776
(Complete Palm Oil Guide). Jakarta, Indonesia: 3A 107.613 133.104 73.033 134.595
Penebar Swadaya. 3B 125.211 138.418 85.341 141.739
Parrot, S.A. (2017). Parrot Sequoia Capture the Iinvisible 3C 102.690 124.847 73.835 121.059
Monitor Your Crops. [Online]. Avaiable at: 3D 114.509 134.506 80.459 137.144
https://www.parrot.com/us/business-solutions/parrot-
3E 103.420 133.144 73.862 128.594
sequoia [Accessed 1 January 2018].
3F 101.884 131.221 69.520 123.551
Pusat Data dan Sistem Informasi Pertanian. (2014). Outlook
3G 105.717 133.299 69.451 131.241
Komoditi Kelapa Sawit (Palm Oil Commodity
Outlook). Jakarta, Indonesia: Pusat Data dan Sistem 3H 108.704 131.464 80.579 128.028
Informasi Pertanian Sekretariat Jenderal-Kementerian 3I 121.342 135.507 97.446 132.733
Pertanian. 3J 122.091 132.585 98.171 135.082
Shofiyati, R. (2011). Teknologi Pesawat Tanpa Awak untuk 3K 134.299 137.072 99.472 146.368
Pemetaan dan Pemantauan Tanaman dan Lahan 3L 114.651 134.315 84.555 141.575
Pertanian (Unmanned Aerial Vehicles Technology for 4A 111.116 137.428 77.785 136.967
Mapping and Monitoring Crops and Agricultural Land). 4B 90.891 131.026 72.818 119.278
Informatika Pertanian, 20(2):,pp. 58–64. 4C 79.805 129.379 66.129 122.606
Starks, P.J., Zhao, D., Phillips, W.A. and Coleman, S.W. 4D 85.741 126.997 73.204 123.770
(2006). Development of Canopy Reflectance 4E 96.247 141.433 77.157 131.028
Algorithms for Realtime Prediction of Bermuda Grass
Pasture Biomass and Nutritive Values. J. Crop Sc,.
46(2), pp. 927-934.
134