visual and thermal lowvisual and thermal low--altitude ...€¦ · saeed yahyanejad...

36
Visual and thermal low Visual and thermal low-altitude aerial imagery for altitude aerial imagery for real real-time monitoring time monitoring real real time monitoring time monitoring by: Saeed Yahyanejad

Upload: others

Post on 15-Sep-2020

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Visual and thermal lowVisual and thermal low--altitude aerial imagery for altitude aerial imagery for realreal--time monitoringtime monitoringrealreal time monitoringtime monitoring

by: Saeed Yahyanejad

Page 2: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Outline

1 Motivation1. Motivation 2. Challenges3 Problem definition3. Problem definition4. Orthorectified mosaicking of UAV images

(A quick overview)(A quick overview)5. Multispectral mosaicking

(Main focus of this presentation)(Main focus of this presentation)6. Conclusion 7 Future work7. Future work

2Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 3: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Motivation

Advances in computer vision and aerial imagery.

Small-scale UAVs - easy to deploy. y p y

Provide overview images similar to Google Maps and Bing g g p gMaps, but with higher temporal and spatial resolution.

Applications of aerial images: environmental monitoring, surveillance, border control, detection, and disaster

3Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

management.

Page 4: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Goal:  Generate overview image

4Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 5: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Challenges of aerial image mosaicking

Non-planar scene and low altitude P i f i i i i i i Point of view in image acquisition Overlapping ratio

5Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Panorama: extension to FOV Mosaic: extension to POV

Page 6: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Problem definition:  (How to mosaick?)

6Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Geosense. http://www.geosense.com.my/sample_uav_image/Index.htm

Page 7: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Problem definition(A quick overview)(A quick overview)

1 H t th t d t (IMU & GPS) f i ki ?1. How to use the metadata (IMU & GPS) for mosaicking? Contribution: A hybrid approach: find adjacent images + appro imate positions and orientationapproximate positions and orientation

2 H t d th tifi d i ith t i th2. How to produce orthorectified mosaics without using the metadata?

Contribution: Study the sources of error + minimize theseContribution: Study the sources of error + minimize these accumulated errors (low overlap no global optimization)

7Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 8: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Problem definition(main focus of this presentation)(main focus of this presentation)

3. How to perform a robust interspectral registration3. How to perform a robust interspectral registration between images taken from heterogeneous sensors (e.g., thermal and visual cameras)?

Contributions:

• General lens distortion correction of thermal camerasGeneral lens distortion correction of thermal cameras• Introduction of feature descriptors for robustly identifying correspondences between images of different spectrumsp g p• Registration of image mosaics• Registration based on depth mapsRegistration based on depth maps

8Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 9: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Orthorectified mosaicking of UAV images

Mosaic On is constructed incrementallyn y

where O0 is a zero matrix &Ii’ is a transformation of original image Ii by:Ii is a transformation of original image Ii by:

Projectivetransformation:

9Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 10: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Orthorectified mosaicking of UAV images

Incremental mosaicking:

10Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 11: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

1) Hybrid method

Reduce the search space:

1.Place the images based on camera posebased on camera pose

2.Use image-processing to fine tune

11Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 12: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

2) Loop‐independent mosaicking

Major sources of error in pairwise mosaicking:

1. Unleveled features

There is a loop if there is a simple cycle in the

2. Lens distortion

a simple cycle in the adjacency graph.3. Projection and

transformation model

12Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 13: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Unleveled features ‐ Disparity implies the depth

13Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 14: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Loop‐independent mosaicking result

No error control With considering sources of error

14Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 15: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Multispectral aerial imaging

multiple heterogeneous sensors - needs alignment

Visual cameraThermal camera

15Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 16: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Multispectral aerial imaging

General goal: Register visual and thermal imagesChallenges:

• Different features because of daylight, IR radiation, fog, y g gsmoke, water and vegetation

Partial goal: Correct thermal lens distortionChallenges:

• Conventional patterns (e.g. chessboard) is not easily applicable for thermal

16Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 17: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Correct thermal lens distortion

ti lconventional chessboard

pattern

Easy toEasy to reimplement

&Accurate

IR radiation,

Accurate

conveniently generated by a heat lamp

Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt 17

Page 18: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Correct thermal lens distortion

N l i f hNormal image of the pattern:

Result of our assembly:

Now it is possible to feed it in existing calibration SW

Adaptive thresholding:Erosion + Convex hull:Dilation + Convex hull:

Centers of blobs:

Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt 18

Page 19: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration ‐ State of the artInterspectral registration  State of the artInterspectral registration (e.g., finding correspondences between thermal and visual images)between thermal and visual images)

R. Istenic, et al. Thermal and visual image registration in Hough parameter space.14th International Workshop on Systems, Signals and Image Processing (IWSSIP), pages 106-109, 2007.

G. Schaefer et al. User Centered Design for Medical VisualizationUser Centered Design for Medical Visualization, chapter Automated overlay of infrared and visual medical images, pages 174-183. IGI Global, 2008.

19Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 20: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

State of the art

E. Corias, J. Santamaria, and C. Miravet. A Segment-based Registration Techniquefor Visual-IR Images. O ti l E i i (1) 1 29 2000Optical Engineering, (1):1-29, 2000.

J W J J W Ch i d D L Ch J. W. Joo, J. W. Choi, and D. L. Cho. Robust registration in two heterogeneous sequenceimages on moving objects. In Proceedings of Sixth International Conference ofIn Proceedings of Sixth International Conference of Information Fusion, pages 277-282, 2003.

20Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 21: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration pipeline

Interspectral RegistrationSingle pair

Multiple pairsMultiple pairs

Without overlap With overlap

Split to individual single pairs

Extract Robust Features (e.g., SIFT, SURF, or RFAE)

21Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Register the images

Page 22: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration – single pair

Scale-, rotation-, and illumination-invariant Manipulate existing descriptors (e g SIFT or SURF) Manipulate existing descriptors (e.g., SIFT or SURF) Robust Features Along the Edges (RFAE)

where B is the binary operator and S is the Sobel operator

8 SURF79 RFAEfeatures

Failed Registration

features

SuccessfulRegistrationRegistrationRegistration

22Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 23: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Single pair registration – Satellite sample

RFAE Successful RegistrationSURF Failed Registration

23Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Bands 1 and 4 of the Landsat satellite image of Iowa state (image source: NASA/USGS).

Page 24: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Single pair registration – Human sample

SURF Failed RegistrationRFAE Successful Registrationg g

Visual and thermal images of humans

24Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 25: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration evaluation – single pairSuccessful registration ratio

SIFT SURF U-SURF RFAE SURF+ RFAE

Satellite (low deviation) 24/24 24/24 24/24 24/24 24/24 Satellite (high deviation) 7/10 6/10 9/10 10/10 10/10

Human 0/16 0/16 2/16 12/16 12/16 Surveillance 1/14 2/14 4/14 12/14 12/14

UAV 13/20 14/20 2/20 13/20 17/20

Average quality (Q values are calculated based on coverage and sparsity of features)SIFT SURF U-SURF RFAE SURF+ RFAE

Satellite (low deviation) 72 % 66 % 71 % 53 % 68 %

Satellite (high deviation) 45 % 41 % 49 % 46 % 52 %

Human 0 % 0 % 1 % 8 % 8 %

25Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Surveillance 1 % 3 % 4 % 18 % 19 % UAV 9 % 10 % 1 % 6 % 11 %

Page 26: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration pipeline

Interspectral RegistrationSingle pair

Multiple pairsMultiple pairs

Without overlap With overlap

Split to individual 

Alternative 1 Alternative 2

single pairsMosaic the images from each spectrum together.  Output is p g p

a single pair of mosaics.

Extract Robust Features (e.g., SIFT, SURF, or RFAE)

26Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Register the images

Page 27: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration – multiple pairs

Register the mosaics if they have enough overlaps

27Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 28: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration pipeline

Interspectral RegistrationSingle pair

Multiple pairsMultiple pairs

Without overlap With overlap

Split to individual 

Alternative 1 Alternative 2

single pairsMosaic the images from each spectrum together.  Output is 

Construct a depth map from each spectrum.  Output is a p g p

a single pair of mosaics.p p

single pair of depth maps.

Extract Robust Features (e.g., SIFT, SURF, or RFAE)

28Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Register the images

Page 29: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration – multiple pairs

Exploit the 3D structure – rough depth map

+ =

29Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Depth map registration

Page 30: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration – multiple pairs

Exploit the 3D structure – full reconstruction

3D model from 25 visual images

3D model from 25 thermal images

30Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 31: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Interspectral registration – multiple pairs

Project 3D to 2D – equivalent depth map

Perpendicular base point on plane Distance of the point to the plane

SURFover

d thdepth map

31Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 32: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Multiple pair registration ‐ depth map sample

Visual Thermal

Successful registration using 3D structure- different seasons and sensorsThermalWinter

Summer WinterThermalSummer

32Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Surf failed - same sensor but different seasons

Successful registration using 3D structure different seasons and sensorsWinter Summer

Page 33: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Alignment sample – Multispectral image

33Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 34: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

ConclusionWe construct high-resolution mosaic images with high geometric accuracy from a set of images taken from small-scale UAVs.

Present a hybrid approach which combines metadata and the image data (increased accuracy and 70% computational time reduction)g ( y % p )

Quantify the influence of different parameters - studying the sources f i i i i kiof error in pairwise mosaicking

(increased accuracy with the same complexity)

Robust interspectral image registration• A general method (RFAE) which exploits the existing scale-

invariant feature extraction methods such as SIFT and SURFinvariant feature extraction methods such as SIFT and SURF• Register the image mosaics (complexity n times more, for n images )

• Use depth maps of a target scene for the feature extraction

34Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

p p g(For a rough depth map 4n times more, for n images)

Page 35: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Future work

Advanced UAVs: handle higher resolution or videos

Hybrid approach: improve optimization (heuristic)

Loop-independent mosaicking: generalize in presence of loop

Interspectral registration: more complex features (e g classified areas or object can be used SGM)(e.g., classified areas or object can be used, SGM)

System level future works: Underlying architectures for deployment ofSystem level future works: Underlying architectures for deployment of multi-UAV, autonomous planning and deployment, communication structure, and optimal coverage methods

35Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt

Page 36: Visual and thermal lowVisual and thermal low--altitude ...€¦ · Saeed Yahyanejad ‐Alpen‐Adria‐UniversitätKlagenfurt 8. Orthorectified mosaicking of UAV images Mosaic O n

Questions?Q

fThank you for your attention!

36Saeed Yahyanejad ‐ Alpen‐Adria‐Universität Klagenfurt