Precision Agriculture Precision Agriculture is about

0 downloads 0 Views 5MB Size Report
Precision Agriculture (Satellite Farming) has faced a breakthrough shift since risings of UAVs ... market in Malaysia to be hit by UAVs. Next. Prev ... Small size.
Redmond Ramin Shamshiri, PhD [email protected]

Next

 Precision Agriculture (Satellite Farming) has faced a breakthrough shift since risings of UAVs  Conventional Remote Sensing platforms are being replaced by integrated UAVs  UAVs & Robotics are the future of Precision Ag  Precision Agriculture of Oil Palm is one of the largest market in Malaysia to be hit by UAVs

Precision Agriculture

GPS GIS

Wireless Sensor Networks

Automation Control

Robotics

PA is a farming Management concept based on Sensing, Measuring and Assessment

RGB Infrared NIR Hyperspectral Multispectral Thermal Mapping software, Mobile Apps

Sensing Components

Sensors , Cameras o o o o o o o

Typical Applications in Agriculture UAV

Remote Sensing

VRT

           

As a flexible remote sensing platform Crop/ Tree Scouting Health/growth assessment Inventory management Yield estimation / Monitoring Weed and disease detection Mapping (2D, 3D, GIS, NDVI) Risk/Hazard/Safety management Soil condition assessment VRT Robotics harvesting Academic and Research application

Prev

RGB Hyperspectral Multispectral Thermal NDVI Infrared NIR

Platform

Airborne

Ground-based Handheld

Satellite

Piloted Airplanes

Next

Vehicle mounted

UAV

Fixed wing

Multi-rotor

Prev

Traditional Scouting o o o o o o o o o o

Next

Precision Agriculture is about optimizing returns on inputs while preserving resources

Integrated with Commercial sensors/Cameras/Software       

Prev

Next

Plant densities are an important and limiting factor for growth, Nutritional status, fruiting and hence for a plantation’s yield.

Traditional scouting requires spending hours and hours of visualizing Involves manual operations, ineffective techniques Limitation of human/labor resources (Not on demand) Repeated dull tasks (i.e., Palm census) Not a pleasant environment to work (hot, humid) Hazard/ Safety (Falling from trees, bugs, snakes, etc) Generates Inaccurate/biased statistics Requires expert knowledge/Post processing (i.e., lab analysis of data) Generates limited information (Does not provide comprehensive result) Ignored parameters due to measurement’s difficulties (i.e., tree height,

Oil palm with good plant density

Oil palm with high mortality

 Optimal plant densities depend on different factors, such as cultivars, climate, soil characteristics, land preparation…  Refilling of canopy gaps and correction of non-optimal plant densities are of high priority for a good plantation management

canopy diameter, tasks that involves climbing trees)

Satellite remote sensing o Cost o Low resolution o Difficulties of access (Not on-demand)

Yield reduction due to high density palm areas that causes Etiolation*.

Ground sensing o Time consuming o Limited field of view

Need for Automatic identification of potentially etiolated palms or high density areas Conventional method, solely based on visual observation, inaccurate, particularly when coverage is large and dominant topography is hillocky. * Etiolation is a process in flowering plants grown in partial or complete absence of light.

Prev

Next

Prev

Next

1

Main Objective

Specific Objectives

To develop/adapt a flexible remote sensing platform by integrating UAV, Sensors, and robust machine vision system for Smart Management of Oil Palm Plantations

1. To setup and operate a multi-rotor UAV remote sensing system for Oil Palm plantations 2. To produce 2D and 3D visual maps of the fields under study

Smart Management involves:

(including Blocks, rivers, roads, boundaries)

3. To produce NDVI (Vegetation stress) and GIS maps for health and growth assessment

A: Smart Inventory Management

4. To develop a robust, real-time machine vision system for the following tasks

B: Smart Growth/Health Assessment

(4.1) Inventory management Palm census and density, track and record Crown diameter estimation, Canopy size Palm height measurements Plantable, Unplantable, Overplanted areas Palms distance

“Smart” refers to: Autonomous monitoring, data processing and decision making

(4.2) Yield mapping system Detection/Quantification of fresh fruit bunches from UAV images Yield per Palm and Yield per ha Prev

Next

Prev

Next

Prev

Next

Prev

Next

Prev

Next

Prev

Next

Consideration for a flexible design            

Small size Light-weight Affordable Autonomous Stable (against wind and other disturbances) Shifting between multiple sensors Payload Flight time Safety (For operator, environment, and platform) Low-altitude flight On-Demand flight (Ease of access operation) Repair and Maintenance costs

2

Prev

Next

Billion US Dollar

Oil Palm contribute to 12 billion USD of Malaysian economy 20 18 16 14 12 10 8 6 4 2 0

Potential for monitoring oil-palm plantations in such a great detail has been never possible

Prev

Next

This research will provide growers/managers with a tool for:  Automatic palm detection, counting, size measurements, etc

Palm Oil Export ($)

 Calculation of planted areas (for replanting or thinning)  Analyzing Palm status based on orthomosaics and digital elevation models  Generating valuable information based on each and every individual palm  Classification of palms based on crown size, height, vegetation indices, etc  Such information can be used for appropriate management decisions (maximize yields)

60

% of the total world

50

Yield Estimation Model development

Palm Oil Export (% of the total)

Correlation between palm height (𝒙𝟏 ) , crown size (𝒙𝟐 ), age (𝒙𝟑 ), vegetation index (𝒙𝟒 ) , …, and yield

40 30

𝒀𝒊𝒆𝒍𝒅 = 𝒇𝒖𝒏𝒄(𝒙𝟏 , 𝒙𝟐 , 𝒙𝟑 , 𝒙𝟒 … )

20 10

 A model that is based on a comprehensive information of each palm location, size, and health, will provide managers with an estimation of yield, and make decisions for sustainable practices methods for production increase without necessary needs for expanding the plantation into natural forests

0

Prev

                      

Prev

Next

Prev

Next

Methodology steps

Research phases

Accurate planted area Creation of Palm trees inventory database for specific plot Total land use Palm distances to specific areas Canopy diameter estimation Tree height measurements Calculating palm density for specific plot Creation of 2D, 3D, GIS, NDVI maps for plantation Monitoring Healthy/Unhealthy palms (Stress assessment) Monitoring exposed soil (VRT application) Quantification of FFB, Estimation of mature fruits Calculating yields for each palm from the acquired images Yield monitoring Creation of yield maps Chlorophyll analysis Drought assessment Biomass indication Leaf area index Growth monitoring Weed detection Inventory management decision support systems Yield Model Development Academic and Research application

1. 2. 3. 4.

Platform setup, integrating UAV, sensors and software Creating high quality 2D and 3D maps of the area Developing custom-built programs/algorithms for smart inventory management Developing custom-built programs/algorithms for smart Health/growth assessment

DATA COLLECTION Data Processing Mapping Modeling

UAV images need to be photogrammetrically processed and translated into accurate 2D orthomosaics and maps, 3D models and surface models, and other GIS datasets

UAV Setup Flight Preparation Test and trials Mission planning

Image acquisition Video streaming Camera/Sensor setups Calibration

Image processing Creating 2D, 3D, GIS, NDVI Maps (Pix4Dmapper, Agisoft)

Results / Reports

Image/Map interpretation GIS analysis Custom software

Data analyzing

Management strategies Decision makings

Correlation analysis

Reports generation

Orthomosaics, 3D Models and Digital Surface Models, 3D Flythrough Videos, Multispectral Image Mosaics, Index Maps, (i.e., NDVI) needs to be processed/interpreted

Next

Prev

Next

3

DJI Phantom, a small hobbyist quadcopter

Purchase, Adapt, Or Build from Scratch? Sensors (Cameras)

UAV Platform

Hyperspectral

Application Software

Programming (Processing)

RGB

Multi-rotor

Fixed wing

Factors to be considered       

NIR Onboard GPS

Auto flight Controller

LiDAR

Multi-rotor UAVs  launch and land vertically  are favored where space is tight Fixed-wing UAV Requires suitable space to launch and land  Can provide longer flight duration  Can carry a heavier payload.

Built from scratch UAV

 Commercial UAVs are expensive,  Not designed for operation inside oil-palm plantations  Might need protection frames

NDVI Thermal

Quadcopter micro UAV 'microdrones md4-200' with the prototype MSMS multispectral sensor

Setting up a UAV Multispectral

Fixed Wing, ebee

The quadrocopter UAV, model md4-1000

Low-cost: affordable by oil palm growers Modular sensor system (shifting different sensors) Higher flight time, more ha coverage Payload Mobile applications Data sharing Wireless networks

MōVI M10 - Digital 3-Axis GyroStabilized , Price: $4,995

The camera system on board the Oktokopter

In-flight stability, flight-time duration and payloads are major paramount concerns

Prev

Next

Visible (RGB)

Prev

Next

Normalized Difference Vegetation Index (NDVI): measurement of the amount of live vegetation in an area

Infrared\ Near Infrared

NDVI - TWIN LUMIX LX7 low cost, rugged, high resolution imaging solution for NDVI, agricultural and archeological data analysis. provides true R and IR information from different sensors, providing the clean photogrammetric information required for vegetation anaysis processing. Combined with a 5cm resolution at 400ft this sensor provides a host of benefits for the agronomist and archeaologist . True Resolution:10.1Mp, Aperture Range:F1.4 - F2.3 Max Shutter Speed: 1 second, ISO:80 - 6400 ISO Image Quality: 82, Resolution (GSD): 5cm @400ft Format:RAW, JPG, Image Stabilisation:roll and anti-shake

SUN

Sony A6000 Visible Sensor

Lumix LX7 Sensor Details

Lumix LX7 Infrared Sensor

Focus Points: 179 True Resolution: 24.3Mp Pixel Size: 15.1 µm² ISO: 1347 ISO Image Quality: 82 Dynamic Range: 13.1 EV Colour Depth:24.1 bits Resolution (GSD):1.5cm – 4.5cm

True Resolution:10.1Mp Aperture Range:F1.4 - F2.3 Max Shutter Speed: 1 second ISO:80 - 6400 ISO Image Quality: 82 Resolution (GSD): 5cm @400ft Format:RAW, JPG Image Stabilisation:roll & anti-shake

True Resolution:10.1Mp Aperture Range:F1.4 - F2.3 Max Shutter Speed: 1 second ISO:80 - 6400 ISO Image Quality: 82 Resolution (GSD): 5cm @400ft Format:RAW, JPG Image Stabilisation:roll and anti-shake

X%

Canon S110-NIR, 12 MP, adapted to be controlled by drones autopilot Acquires image data in the NIR band

Y%

Resolution: 12 MP Ground resolution at 100 m: 3.5 cm/px Sensor size: 7.44 x 5.58 mm Pixel pitch: 1.86 um Image format: JPEG and/or RAW

𝑋% − 𝑌% ≤ +1 𝑋% + 𝑌%

NDVI0.66 Very healthy plants

Healthy plants have a strong near infrared reflectivity, called the "Red Edge".

Prev

−1 ≤ 𝑁𝐷𝑉𝐼 =

RGB

NIR

NDVI Canon PowerShot SX260 12.1 MegaPixel Digital Camera Models: XNiteCanonSX260: UV+Visible+IR XNiteCanonSX260: IR- Only XNiteCanonSX260NDVI: 3-Band Vegetation Stress Remote Sensing Camera

Next

Prev

Next

Object-Oriented Software Framework for Hyperspectral Imaging Thermal Sensors

Multispectral Sensor    

High Resolution Uncooled Thermal Camera: Tau 640



To identify Oil Palm stress factors, soil types, fertilizers, or insecticides To identify differentiate plant species or recognize other plant (weeds…) To identify soil or chemical conditions that are, in each case, able to be identified by their unique spectral signature. To graphically illustrate vegetation indices such as NDVI that are defined by relationships of specific narrow-band wavelengths. With each exposure, 4 or 6 separate bands of visible or near-infrared radiation move through each camera's lens and filter to form a separate monochromatic image on the camera's sensor.

Hyperspectral Imaging

RGB Color Image Pixel 3 dimensional data at R, G, B waveband

High Resolution thermal imaging can assist disease detection and water stress in Oil Palms, or for scouting at nights, fire hazard alarm

Optris PI450 Temperature range: -20°C to 900°C Spectral range: 7,5 bis 13 µm Optical resolution: 382 x 288px Frame rate 80 Hz capture single images at a rate one per 2 seconds. Each pixel from each image has an exact temperature associated with it.

Prev

Next

Prev

Next

4

Hyperspectral image contains enormous information. Pixel dozens or hundreds dimensional data

Advantages • We can analyze the target in detail. • We can combine both spatial and spectral analysis. • We can apply it to the wide area of agricultural sensing.

[-] Intensity [-] Intensity

250

200 200

Disadvantages • Data structure is complicated and special. • Data size is large. • There are few software applications and libraries. • We have to develop the software by ourselves. • We have to understand complicated data structure.

150 150 100 100 50 50 0

0 0

25 100

50

150

50 200

Waveband Spectral Band [-] [-]

Prev

Next

Prev

Next

Mechanism of Hyperspectral Camera (ImSpector)

400 nm from http://www.specim.fi/ functions as line sensor

Camera

Wavelength spectral axis

Camera

Camera axis

1000 nm

Target

Spatial axis by Electric moving pan head

by Driving vehicle

Prev

Next

Prev

Next

Dataflow of Sampling Pixel Spectral Data Sensing line frame

Spectral image A

B B-D

A-C

y (x, y)

x Pixel spectral data 250 Intensity [-]

200

C

D

150 100 50 0 0

50

100

150

200

Spectral Band [-]

Prev

Next

Prev

Next

5

Extraction of Desired Waveband Value ("WavebandProcessor" object) Outdoor Fields Hyperspectral Camera

Target

Digital Video Recorder

Indoor Laboratory

IEEE 1394 Interface

Personal Computer

Sugar beet Waveband No.15(500nm)

Prev

Next

Composition of 3 Waveband Values as False Color

("FalseColorProcessor" object)

Sugar beet Waveband No.35(750nm)

Prev

Next

Distinction between Plant and Soil ("PlantSoilProcessor" object)

1: Plant

0: Soil Sugar beet NIR color (NIR, R, G)

Soybean (crop row) Visible rays(R, G, B)

Sugar beet

Prev

Next

Estimation of Chlorophyll Content (SPAD Value) ("SpadProcessor" object)

24.7

Soybean (crop row)

Prev

Next

Plant Classification ("PlantClassificationProcessor" object)

50 - 55 45 -

1: Sugar beet

40 35 -

26.0

0: Soil

30 25 -

36.9

40.2

20 15 -

2: Sugina

10 5-

24.7

0-

Estimated SPAD values of sugar beet

Prev

Estimated plant ID number (Soil, Sugar beet, Sugina) Next

Prev

Next

6

Light Detection and Ranging (LiDAR)         

Topographic Survey and Mapping Contour Mapping Cross Section / Longitudinal Analysis 3D Mapping Floodplain Mapping Vegetation Mapping Shoreline Analysis Corridor / Route Studies Slope Analysis

3D Point Cloud of LiDAR

Ortho-mosaic of Aerial Photo

Routescene LidarPod

Potential Application Palm height measurements Prev

Next

Price 8700 USD unlimited use in time 3500 USD /year (rent) 350 USD /month (rent)

Convert thousands of UAV images into  Geo-referenced 2D maps  Orthomosaics models  3D surface models  point clouds Orthomosaics

Prev

3D Models

Next

Price: £400.00

Vegetation, Health, growth analysis SampleVideoAgisoft

 Process thousands of aerial images  Suitable for a non-specialist operator  Generate high- resolution Geo-Ref orthophotos  Exceptionally detailed Geo-Ref DEMs*  Fully automated workflow  Easy integration to the Q-Pods system  Create projects using more than one camera and process imagery together (i.e., NIR and RGB)

SampleVideoPIX4D

* DEM: Digital Elevation Model

Prev

Next

Prev

Next

Prev

Next

Creating smart data for decision support systems GIS tasks  Essentials of GIS & Aerial Image Interpretation  Map Generation  Creating Workflows  Management of Spatial Data  Natural Color Images  Multi-spectral Images  Digital Elevation Models  Multi-temporal Images Images source: Adapted from Terracentra

Prev

Images source: Adapted from Terracentra

Next

7

Plantation Infrastructure Inventory, road mapping, inventory and monitoring is very important for efficient plantation management.

Estimated/Precise Palm tree counts in a selected area of interest

Fertility Mapping, detecting and mapping of oil palm fertility or palm vigorous growth level

Identifying unhealthy palms

Plot size: (ha) Palm counts: 2500 Palm density: 100 trees/ha

Prev

Next

Prev

Next

Identification of poor spots examples

Finding average distance of a plot to the river/road

Prev

Next

Prev

Next

Tree Counting methods

Calculating total unplanted area from Geo-referenced maps Training rules

Commercial Software

Other techniques

S5

S4

S3 S1

S2 S6

Total Unplanted Area =

PLOT B

𝒏 𝒊=𝟏 𝑺𝒊

Images source: Adapted from Terracentra

Prev

Next

Prev

Next

8

Study Area=10.3305 Sqr.Inch Estimated count= ?

Training Area = 0.9 Sqr.Inch Training Count = 4

1 2

7

13

8

14

9 3

15 10

4

16

11

21

24 25 26 27

22

5 6

19 20

17

12

18

23 0.5

30

36 37

31

42 43

38

32

44 39

33 28 29

Banana plantation in Indonesia

open-source software QGIS The Leading Open Source Desktop GIS

40 34 41 35

45

The software allows plant counting, density calculations and the generation of mortality maps by visual inspection of the image products.

46

Pineapple plantation in the Philippines

0.4

Manual count ≅ 47

Study Area

Factor =

______ Training Area

Estimated count = Factor × Training count

10.33 = ---------- = 11.47 0.9

= 45.9 Accuracy: 97.6%

Prev

Next

Prev

Next

ERDAS IMAGINE's Spatial Modeler (NIR thresholding) The overall performance (detection rate) between 0.916 to 0.998.

RGB, NIR

NIR

EDF Image

LPF Image

Threshold Image

Prev

Next

Photogrammetric point clouds Technique

Prev

Next

Photogrammetric point clouds Technique

photogrammetric point clouds

Digital surface model

Local maximum display tree position The mapping accuracy amounts for 86.1% for the entire study area and 98.2% for dense growing palm stands.

Prev

Next

Prev

Next

9

Correlating between palm heights, ages, and yield

Prev

Next

Crown (Canopy) Volume Crown diameter

Correlating between image and mass of FFB

Prev

Next

Geo-referenced

FFB detection/quantification

Prev

Next

 Autopilot and mission control  Visual servo control mechanism for FFB detection  Breakthrough innovative ideas (i.e., Night-mission flights)  Making smarter UAV platforms that can learn while flying  Techniques for Improving accuracy and resolution

Y=Weight (Kg)

 Development of customized sensors with built-in algorithms for specific task  Improvement of low altitude flight mission  Improvement of flight control over actuator limits and noise (i.e., Controller Design for Stabilizing of an Autonomous fixed wing Crop Surveillance Osprey Drone with Actuator Limits and Sensor Noise)

X=% of Indexed pixels, No. of FFB, etc

Prev

Next

Prev

Next

10

     

 UAV can contribute to mechanization of Oil Palm Agriculture

Reducing labor force on field Developing GIS dabase, 2D, 3D, NDVI, and thermal maps Reducing labor hazards Reducing management time Palm tree tagging Monitoring fungal disease with different sensors

 FFB quantification is the first step toward building cost-effective robotic harvesting system for existing palm trees

Potential to be extended to other fields in Malaysia, i.e., rice and rubber Embedded Child Script ROS nodes

Further contribution

Plugins

      

Smart pesticide control Enhanced pollination Constant track and record of growth condition Drastically help growers in decision making Early warnings for disease A ground for autonomous robotic harvesting Academic application

Control Mechanisms Serial port Remote API

Preliminary design Simulation

Re-design

Prev

Next

Nutrient contents, N, P, K, Mg, B

Prev

Field experiments

Next

Link to draft

RGB NDVI

Evaluation Improvement

Prev

Predicted Yield Disease detection Management decisions

NIR

Building the prototype

Rubber plantation

Link to draft

Palm heights, density, Crown Diameter, Volume, FFB quantification, etc

Next

Prev

Next

Picture of an 8 year Oil Palm that produced 535kgs of FFB in one year

PrecisionHawk will presented its drone platform for early disease detection on February 4, 2016 at Dubai Internet City. Prev

Next

Prev

Next

11

Prev

Next

Prev

Next

Prev

Next

Prev

Next

Prev

Next

Prev

Next

5 LITERS $5,299

10 LITERS $8,399

15 LITERS $10,399

12

(UPM) (UPM)

(Univ of Florida) (Univ of Florida) (Univ of Florida) (Wageningen UR) For their insightful suggestions and ideas

Prev

Next

13