Next Article in Journal
Factors Driving Changes in Vegetation in Mt. Qomolangma (Everest): Implications for the Management of Protected Areas
Previous Article in Journal
Knowledge-Aided Ground Moving Target Relocation for Airborne Dual-Channel Wide-Area Radar by Exploiting the Antenna Pattern Information
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seamless Mosaicking of UAV-Based Push-Broom Hyperspectral Images for Environment Monitoring

1
School of Geoscience and Surveying Engineering, China University of Mining and Technology-Beijing, Ding 11 Xueyuan Road, Haidian District, Beijing 100083, China
2
Department of Geography and Planning, University of Toronto, Toronto, ON M5S 3G3, Canada
3
School of Geographical Science, Fujian Normal University, No. 8 Shangsan Road, Cangshan District, Fuzhou 350007, China
4
Key Laboratory of Computational Optics Imaging Technology, Aerospace Information Research Institute, Chinese Academy of Sciences, No.9 Dengzhuang South Road, Haidian District, Beijing 100094, China
5
School of Opto-Electronics, University of Chinese Academy of Sciences, No.19(A) Yuquan Road, Shijingshan District, Beijing 100049, China
6
Nanjing Institute of Software Technology, Institute of software Chinese Academy of Sciences (ISCAS), Building 4, Artificial Intelligence Industrial Park, No. 266, Chuangyan Road, Kirin Science and Technology Park, Nanjing 210000, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(22), 4720; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224720
Submission received: 15 October 2021 / Revised: 16 November 2021 / Accepted: 19 November 2021 / Published: 22 November 2021

Abstract

:
This paper proposes a systematic image mosaicking methodology to produce hyperspectral image for environment monitoring using an emerging UAV-based push-broom hyperspectral imager. The suitability of alternative methods in each step is assessed by experiments of an urban scape, a river course and a forest study area. First, the hyperspectral image strips were acquired by sequentially stitching the UAV images acquired by push-broom scanning along each flight line. Next, direct geo-referencing was applied to each image strip to get initial geo-rectified result. Then, with ground control points, the curved surface spline function was used to transform the initial geo-rectified image strips to improve their geometrical accuracy. To further remove the displacement between pairs of image strips, an improved phase correlation (IPC) and a SIFT and RANSAC-based method (SR) were used in image registration. Finally, the weighted average and the best stitching image fusion method were used to remove the spectral differences between image strips and get the seamless mosaic. Experiment results showed that as the GCPs‘ number increases, the mosaicked image‘s geometrical accuracy increases. In image registration, there exists obvious edge information that can be accurately extracted from the urban scape and river course area; comparative results can be achieved by the IPC method with less time cost. However, for the ground objects with complex texture like forest, the edges extracted from the image is prone to be inaccurate and result in the failure of the IPC method, and only the SR method can get a good result. In image fusion, the best stitching fusion method can get seamless results for all three study areas. Whereas, the weighted average fusion method was only useful in eliminating the stitching line for the river course and forest areas but failed for the urban scape area due to the spectral heterogeneity of different ground objects. For different environment monitoring applications, the proposed methodology provides a practical solution to seamlessly mosaic UAV-based push-broom hyperspectral images with high geometrical accuracy and spectral fidelity.

1. Introduction

Hyperspectral images have hundreds of narrow and nearly continuous spectral bands that represent the reflectance signals of the observed objects [1]. In early years the imaging spectrometers were often used in multispectral mode due to the technical limitations of the instruments. Fine spatial and spectral resolution creates problems with the signal-to-noise ratio of the sensor used. The number of spectral bands had to be reduced in order to get smaller pixel size as reading the data from the sensor takes some time and using many spectral bands would have resulted in pixels elongated towards the flight direction because plane speed cannot be reduced below a certain critical limit [2]. Recent developments in the miniaturization of electro-optical sensors and unmanned aerial vehicles (UAVs) have established a new era of hyperspectral imaging [3]. Due to the high flexibility of UAV flight planning, hyperspectral images taken by UAVs can be widely used in environment monitoring applications [4,5,6,7,8,9,10,11], such as precision agriculture, species classification [11,12,13], surface parameter retrieval [14,15], etc.
The wider use of UAVs has also pushed also instrument developers. There are two main types of lightweight hyperspectral sensors for UAVs: push-broom scanners [16,17,18] and snapshot imagers [18,19,20,21]. Push-broom hyperspectral scanners are based on a linear array that captures scenes while operating in a push-broom mode and the acquired images have a high spectral sensibility [16]. Snapshot imaging sensors can obtain bundles of frames at the same time to form an image with two spatial dimensions and tens to hundreds of bands [18]. The push-broom sensor system outperforms snapshot image capturing approaches, as the latter systems require a compromise between spatial coverage, spatial resolution and spectral resolution [20]. Spectral resolution of push-broom hyperspectral scanners has improved in recent years. Most of the sensors have 5 nm or better spectral resolution and can be flown with 1 m or smaller pixel sizes. A UAV carrying a push-broom hyperspectral imager that is light weight and has a small size gradually becomes a payload for acquiring multi-temporal and high spatial resolution hyperspectral images.
In hyperspectral data acquisition, the UAV is usually operated at low altitudes, and the implementations always result in some geometric and radiometric distortions within the acquired data such as pixel displacement errors and spectral noises caused by smile effects. Thus there exists a lot of challenges in the post-processing of hyperspectral images for spectrally complex environments [3]. It is often necessary to obtain multiple UAV images to cover a study area with different overlaps between images. Mosaicking the UAV images into geometrical and radiometric high quality hyperspectral images is one of the most essential procedures to meet the application needs of remote sensing environment monitoring [22,23,24]. However, as the push-broom sensor system is an emerging technology that is of high cost, the research work on it is not popularized in science. Only a few research works have been published to mosaic the push-broom hyperspectral images mounted on a UAV [16,17,18].
The ZK-VNIR-FPG480 hyperspectral imaging system is an emerging push-broom sensor with independent property rights in China that offers a good combination of spatial and spectral resolution and has been proved effective in environment monitoring applications [25]. The image has a total of 270 bands, the spectral range is 400–1000 nm, the spectral resolution is 2.8 nm, and the spatial resolution is 0.9m at 1 km flight height. While no comprehensive tools are available for producing mosaics of UAV hyperspectral images acquired by this push-broom scanner, how to process and mosaic the images to get a seamless hyperspectral image is studied in this research.
Generally, the geometric and radiometric correction of the UAV hyperspectral images are two key steps in image mosaicking. In geometric correction, images obtained with different positions and attitudes should be mapped to an identical reference plane. In radiometric correction, the inconsistency in brightness between geometrically corrected images should be minimized, and then the radiometric corrected multiple images should be fused into a mosaicked image. The geometric quality of the mosaicked image is dependent on the geometric correction accuracy, while the radiometric quality of the image is impacted by the sensor calibration, radiometric correction and image fusion procedures.

1.1. Geometric Correction

In geometric correction, the geo-rectification and registration of UAV hyperspectral images are two key steps and have been active research fields recently [17,26,27,28].
In geo-rectification, each image acquired from an imaging system is tied to spatial locations relative to a desired reference frame with a geometric control. This control can be established using GCPs and/or the implementation of a Position and Orientation System (POS) [29]. Images can be rectified by global transformations derived from the relationships between GCPs and the corresponding image points [30] or by a direct Geo-reference scheme [31]. Push-broom scanners have an inherently unstable imaging geometry due to wind-related 3-dimensional motions of the imaging platform during data acquisition, and therefore their post-processing ideally requires accurate navigation data which are typically measured by a high-end direct georeferencing system based on global navigation satellite system (GNSS)/Inertial navigation system (INS) [32,33]. There are a lot of GNSS and INS systems developed by UAV companies all over the world, such as the Sky2 airborne GNSS receiver (Hi-target company in China, https://en.hi-target.com.cn/precision-uav-kits, accessed on 1 November 2021), mosaic-X5 GNSS module (septentrio company in European, https://www.septentrio.com/en/applications/uavandrobotics, accessed on 1 November 2021) etc. Choosing the right module that is optimal for UAV use is the key to successful products. The Sky2 differential GPS module adopts a high dynamic receiver board, which can achieve cm level differential positioning results even at a high speed of 100 km/h. Mosaic-X5 proved to be the ideal PPK-capable module to deliver the accuracy for survey-grade results. Given that the accurate POS information provides the possibility to get the direct Geo-referenced images, the location errors are unavoidable. Using GCPs can increase the geometric accuracy, however, the displacement between pairs of images always exists.
To get a seamless mosaicked image, the accurate tie points in the overlapped area between pairs of images should be extracted and the geometric differences of the tie points should be minimized during image registration. An algorithm was proposed to automatically combine multiple overlapping images for a scene by a 2-dimensional fast fourier transform to estimate the displacement between image pairs by accounting for image rotation and scale variation. In this research work, the selection of a suitable hyperspectral band is critical for best co-registration to deal with the data volume and processing complexity of hundreds of bands; due to the limitation of the computer’s computational capabilities, the Pearson correlation coefficient (PCC) was computed for each band in the hypercube pair by the iterating step first, then the band with the highest PCC value was retained to the mosaicking process [34]. A modified Speeded-Up Robust Features matching method was used to find the tie points between the partially rectified hyperspectral images and the RGB orthophoto to mosaic the UAV images for agricultural monitoring [28]. The above two methods were validated on frame images for which the overlap rate is high. Whereas, in push-broom hyperspectral image acquisition, sufficient overlaps between adjacent images may not be ensured, especially when data collection efficiency is a priority for a large study area with the constraint of UAV battery duration. Moreover, in some natural environments where the target objects are mainly water and forest, there exist even less texture or distinct features that can be detected as the tie points due to the homogeneity of the spectral signatures. Although existing relative transformation models and resampling methods have been analyzed for narrow overlapped image mosaicking, the lack of tie points and weak stereo geometry makes it difficult to get accurately geo-registered images [23].

1.2. Radiometric Correction

Due to the influences of the light variation during a flight, caused by viewing geometry or wind, as well as instrument noise, spectral signatures consecutively acquired over the same area can vary significantly, and hence a stitching line would appear in mosaics of geometrically corrected images. Stitching lines are therefore an issue that should be addressed in order to obtain seamlessly mosaicked hyperspectral images. It is therefore highly useful to develop methods to eliminate the spectral differences and improve the radiometric accuracy of mosaicked images [27].
One solution is to ensure the image data collected at different strips are radiometrically consistent and the obtained surface reflectance estimates are reliable. In recent research, the authors performed experiments to assess the accuracy of the Parrot Sequoia camera and the sunshine sensor to get an indication of whether the quality of the data collected is sufficient to create accurate reflectance maps. The developed workflow is evaluated against data from a handheld spectroradiometer, giving the highest correlation (R2 = 0.99) for the normalized difference vegetation index (NDVI). For the individual wavelength bands, R2 was 0.80–0.97 for the red-edge, near-infrared, and red bands. They stated that surface reflectance is not directly measured by the UAV-based cameras. When an image is captured, the image sensor in the camera records the radiant (light) received by each pixel as a digital number (DN). To convert the DN to surface reflectance, a radiometric correction must be performed by (1) applying sensor-related corrections to obtain the radiance received by the camera from the DN, and (2) converting the radiance received by the camera to surface reflectance [35].
Although the above rigorous radiometric correction method is preferred, in practice, another solution for eliminating the spectral distortions, namely relative radiometric correction, is used more frequently. For example, a radiometric block adjustment method was proposed for the UAV image datasets by processing of close-range 2D frame format images with a comprehensive weighting scheme for both solar elevations and reflectance panel observations. However, this method was sensitive to a priori values, in particular to the illumination level, so it is important to devise more accurate measurement methods for the incident irradiance [36]. In other research, the hyperspectral imagery was obtained using a snapshot hyperspectral sensor mounted on a UAV, the radiometric response linearity and radiometric response variations were quantitatively evaluated and a method to correct the radiometric response variation effect was proposed. Comparison with spectra measured with an ASD Field Spec Pro spectrometer (Analytical Spectral Devices, Boulder, CO, USA) was done, and the results show that the radiometric precision was 5% for bands between 500 and 945 nm, and the reflectance curve in the infrared spectral region did not decrease [19].
The focus of this study is on the processing of the push-broom hyperspectral UAV images to acquire a seamless image that has a high geometric accuracy and spectral fidelity. A practical image mosaicking chain of methods that are applicable to different types of land surfaces was proposed. Several methods for image rectification, image registration and image fusion were compared and assessed by using a new set of UAV data acquired in urban scape, river course and forest areas in experiments. Results show that the proposed image mosaicking methodology can not only produce a hyperspectral image with a high image geometrical accuracy but also ensure the spectral fidelity before and after image mosaicking. The following paper is organized as follows: Section 2 introduces the data acquisition and preprocessing principles of the UAV-mounted hyperspectral imagery used in this study and proposes the image mosaicking methodology with comprehensive processing steps including geometric rectification, image registration and fusion. Section 3 presents experimental results and discusses the best strategy for each study area. Section 4 draws conclusions and proposes future work.

2. Materials and Methods

2.1. Data Acquisition

The ZK-VNIR-FPG480 hyperspectral imaging system (ZKYD Data Technology Co., Ltd., Beijing, China) (Figure 1a) includes an six-rotor UAV platform (DJI M600 Pro), an image stabilization platform, a hyperspectral imager, a high-speed data acquisition controller, and a POS. The hyperspectral imager and the Sky2 airborne GNSS receiver were pointed vertically downward on the drone DJI M600 Pro (Figure 1b). The imager is at the center of rotation of the gimbal, which is used to keep the sensor horizontal when the UAV is moving. Table 1 shows a summary of the UAV hyperspectral imaging system properties. The sensor covers the spectral range of 400–1000 nm with a spectral resolution of 2.8 nm and a spatial resolution of 0.9 m at 1 km flying height. The max framerate, which is the maximum lines of data that are acquired in one second during flight, is 100 frames per second.
The UAV-based hyperspectral images are acquired at low altitudes in sunny weather. While the influence of the atmosphere on the images can be almost entirely ignored, a whiteboard with known reflectance data was put as the calibration target to facilitate the conversion of UAV image DN to surface reflectance [26]. To obtain as large an area of image data as possible during one flight, the flight lines were set parallel to acquire the hyperspectral image strips. In each flight line, linear push-broom detectors were used to scan the ground surface along the flight direction with 480 spatial pixels across the flight line, and the spectral signals were recorded and decomposed into 270 bands to construct the arrays of each hyperspectral image. At the same time of imaging, the Sky2 and GNSS satellite receiving antenna records the GNSS base’s real-time kinematics (RTK) signals that are used to derive the UAV external orientation elements. The UAV turned around after reaching the terminal and flew along the next line until the entire study area was scanned.
In different areas, there are different limitations for UAV flights in China; the legislation for the UAV flights in China is published in 2015 by the Civil Aviation Administration of China [37]. As the flight height becomes higher, the photograph scale becomes smaller, the coverage of the area that can be captured in one scanning strip is larger and the number of strips of image that is needed to cover the study area become fewer. In the experiment, data collection of the three typical study areas are implemented. To reduce the time that is needed for data collection as well as consider the latter convenience of image mosaicking, different flight parameters and side overlap rates are set according to the characteristics of the study area before the flight implementation (Table 2).
As shown in Figure 2, the locations of the study areas are labelled using the red polygon. The urban scape area (960 m∗300 m) is a place near Erhai lake in Yunnan Province which is located in the southwest of China. When collecting the hyperspectral images, the UAV flight height is set at 200 m. From Figure 3a, it shows that the urban scape area is composed of buildings, roads, trees, farmland and a lake; it is easy to extract multiple high-precision tie points for image mosaicking, so the side overlap rate is set as 30%. The river course (1100 m∗68 m) is a part of Zhang wei xin river in DeZhou distinct of Shandong Province, which is located in the east of China, the flight height is set as 120m. The area is dominated by water (Figure 3b) where the spectral reflectance is homogeneous so it is difficult to find enough high-precision tie points. Moreover, the width of the river course is about 76% of the width of image strip; therefore, the side overlap rate was set at 50% to make sure there exist some distinct feature points that can be used in image registration. The forest area (570 m∗135 m) is an area that is located in Zhan jiang mangrove forest protection area in Guang dong Province, which is located in the south of China. As shown in Figure 3c, it is representative of mangrove, riparian plants and water. To get a high spatial resolution image for mangrove forest monitoring, the flight height is set as 90 m, so the coverage of each image strip is relatively small. In this situation, the side overlap rate was set at 20% to improve the data collection efficiency.

2.2. Data Preprocessing

To achieve the image strips that were used in the latter image mosaicking, three preprocessing steps were implemented to the original data set. First, the images that have a large distortion at the beginning and end of each flight line were removed from the original data set. Then, the images that were retained were stitched along the flight line to get the image strips. Finally, the sensor-related calibration and reflectance calculation were implemented for each hyperspectral image strip to convert the digital number to reflectance data. In sensor-related calibration, there were three key steps: spectral calibration, relative radiometric calibration and absolute radiometric calibration [38]. Spectral calibration ensures an accurate instrumental spectral response to surface reflectance, while relative radiometric calibration ensures a uniform output of the sensor across space and time. Absolute radiometric calibration determines the conversion from DN to radiance and to reflectance units. First in spectral calibration, the central wavelength and the full width at half maximum of each spectral band are determined. In radiometric calibration, the corresponding relationship between the DN of each spectral band and the true radiance is established and recorded in the calibration files [19]. Next, the radiance of the pixels in a hyperspectral image can be calculated using the corresponding calibration files. Finally, the spectral reflectance of each pixel Rijλ can be calculated using the following equation:
R i j λ = S i j λ D i j λ T λ D i j λ × R r λ
where i, j, λ represent the pixel row number, column number and wavelength, Sijλ is the radiance of the pixel (i, j) in the hyperspectral image, Dijλ is the spectral density of dark current, and Tλ is the radiance of the reference of a calibration whiteboard in the hyperspectral image which has the reflectance R.
The reflectance data of the urban scape, river course and forest area are shown in Figure 3a,b. The data size of Urban scape, river course and forest area is 2.6, 6.81 and 3.06 Gigabytes respectively. After the preprocessing of UAV hyperspectral images, the distortion that occurs in the push-broom sensor’s scanning procedure along the flight line can be mostly reduced because a gimbal was used in the whole data capture process. However, there are still scale, rotation and translation distortions between the narrow overlapped image strips as shown in Figure 3.

2.3. Image Mosaicking Methodology

To get seamless mosaicked image bands, an image mosaicking chain was proposed as shown in Figure 4. Here, the hyperspectral images are the reflectance data calculated from the image strips’ DN data, as shown in Figure 3. The geometric rectification (Section 2.3.1) is first implemented. Three kinds of geometric rectification methods based on the GCPs, POS and a combined use of POS and GCPs can be chosen. While geometric errors are unavoidable, the image registration (Section 2.3.2) method was then implemented to adjust the coordinates of the tie points on different image strips to ensure that the neighboring image strips can be stitched together by overlaying with no dislocation. Two image registration methods of the SIFT and RANSAC methods (Section 2.3.2.1) and the improved phase correlation method (Section 2.3.2.2) can be chosen. Finally, the weighted averaging and best stitching image fusion methods can be chosen to mosaic the neighboring image strips into a large seamless image (Section 2.3.3).

2.3.1. Geometrical Rectification

Geometric rectification of each strip can reduce the differences and provide the geometric rectified data for image registration. In the POS-based method, the exterior parameters are derived from the POS data, including (Xs, Ys), which is the location of the light beam center, the pitch, yaw, roller which are the three angles that describe the flight attitude. The focal length ƒ is derived from the camera calibration report. The ground is assumed as a flat plain because the area coverage is small, and then the ZA coordinate of each ground point A is set as ZS-h; h represents the flight height. The coordinates of each ground point A (XA, YA) is calculated by direct geo-referencing based on the collinear condition equation [39] as shown in Equation (2). In the equation, the a1, a2, a3, b1, b2, b3, c1, c2, c3 are calculated based on the three angles of pitch, yaw, roller to construct the rotation matrix. By a bilinear interpolation method, each image strip can be rectified and resampled to be the coarsely rectified image with geographical coordinates. In this process, the self-developed hyperspetral image processing software, which is programmed by Visual C++, was used.
{ X A X S = ( Z A Z S ) a 1 x + a 2 y a 3 f c 1 x + c 2 y c 3 f Y A Y S = ( Z A Z S ) b 1 x + b 2 y b 3 f c 1 x + c 2 y c 3 f
Due to the limited positioning and attitude accuracy of the POS system on board, there will be an obvious misalignment for the neighboring strips after direct geo-referencing. To improve the accuracy of geometrical rectification, GCP-based image rectification using the curved surface spline function (CSSF) [40] is applied as shown in Equation (3). Each image strip is assigned with new geographical coordinates after the GCP-based image rectification by resampling it with a bilinear interpolation method. As a non-linear transformation function, CSSF only needs three GCPs to simulate a high-order differentiable smooth surface. After the two image strips were geometrically rectified, the smooth surface passes through the same GCPs, so the difference between the pair of rectified image strips becomes smaller. Here, a reference image in the same area was used to extract the GCPs, and ArcGIS 10.1 software was used to implement the CSSF-based image rectification method.
W ( x , y ) = a 0 + a 1 x + a 2 y + i = 1 n F i r i 2 ln ( r i 2 + ε )
where, a 0 , a 1 , a 2 , F i ( i = 1 , 2 , k n ) are parameters to be estimated, r i 2 = ( x x i ) 2 + ( y y i ) 2 , and ε is the parameter that is used to adjust the surface curvature. Normally ε = 1 ~ 10 2 when the surfaces are mildly distorted, and ε = 10 5 ~ 10 6 when the surfaces are significantly distorted.

2.3.2. Image Registration of Neighboring Image Strips

Because noise in the hyperspectral image of one band would result in image matching errors, it is preferred to choose an optimal band that has the least noise and the best image quality to increase the image matching accuracy. Moreover, to make sure that the selected band is suitable for complex scene feature detection, the local variance method based on edge block elimination [41] is used to calculate the image signal-to-noise ratio (SNR) first. Then, the band with the highest SNR was used in the latter image registration.

2.3.2.1. SIFT and RANSAC-Based Method

The traditional feature point-based technical process of image registration is as follows: first the feature points on two images are extracted; then the feature point pairs are matched through feature similarity measurements to obtain tie points; next, the image space coordinate transformation function is established using the tie points; and finally, image registration is performed by coordinate transformation [42].
The image feature point extraction is fundamental for the success of the subsequent image registration. There are many traditional image matching methods, such as scale-invariant feature transform (SIFT) [43] and Speeded-Up Robust Features (SURF) algorithms [42]. The SIFT algorithms run fast and have a high stability in extracting feature point pairs from images with multiple changes such as translation, rotation, and zooming scale. In using a SIFT algorithm, there exist many image matching strategies such as the brute-force matching [44], the K-nearest neighbor matching [45] and the best-bin-first matching methods [46]. For simplicity, the nearest neighbor ratio method [47] is used in image matching first. However, there are some inevitably wrong matched points. So the random sample consensus (RANSAC) algorithm, which is effective and fast in removing some outliers and mismatched points are used for purifying the tie points [48]. After extracting the tie points, the spatial similarity transformation [39] with 7 parameters (3 rotation parameters, 3 translation parameters, 1 scale parameter) was applied to the right image strip with the left image strip as the reference image.

2.3.2.2. The Improved Phase Correlation Method

With consideration that the phase correlation method [49] based on the two-dimensional Fourier transform is scene-independent, and can accurately align two-dimensional images with small translation, an improved phase correlation method was proposed and compared with the SIFT and RANSAC methods. The advantages of this method include its simplicity, its fast running speed, its insensitivity to brightness changes, and its strong anti-interference ability.
First, to deal with the image strips with rotation and scaling, the phase correlation method was improved by transforming the images from the Cartesian coordinate system to the polar coordinate system first, then the rotation relationship between the images is represented as θ = θ + θ 0 , where θ represents the scaling relationship represented as r = r 0 r and transferred to an additive relationship by logarithm operation [50].
The main steps of the phase correlation method accounting for the rotation and scaling differences are:
(1) Transfer the image from the Cartesian coordinate system to the polar coordinate system.
In the Cartesian coordinate system, suppose the relative rotation angle is θ0 and scaling factor is k, then the spatial relationship between the neighboring image strip can be represented in Equation (4):
f 2 ( x , y ) = f 1 [ k ( x cos θ 0 + y sin θ 0 ) , k ( x sin θ 0 + y cos θ 0 ) ]
(2) Then the Fourier transformer of the images F1 and F2 satisfies Equation (5):
F 2 ( u , v ) = 1 k 2 e i 2 π ( u x 0 + v y 0 ) F 1 [ 1 k ( u cos θ 0 + v sin θ 0 ) , 1 k ( u sin θ 0 + v cos θ 0 ) ]
(3) Next, the amplitude spectra of the two images have the following relationship:
M 2 ( u , v ) = 1 k 2 M 1 [ 1 k ( u cos θ 0 + v sin θ 0 ) , 1 k ( u sin θ 0 + v cos θ 0 ) ]
Set { u = ρ cos θ v = ρ sin θ , Equation (6) can be simplified as:
M 2 ( ρ cos θ , ρ sin θ ) = 1 k 2 M 1 [ ρ k cos ( θ θ 0 ) , ρ k sin ( θ θ 0 ) ]
(4) Based on Equation (7), the parameter ρ k can be transformed by the logarithmic transformation lg ρ k = lg ρ lg k , then Equation (7) can be transformed to Equation (8):
M 2 ( θ , lg ρ ) = 1 k 2 M 1 [ ( θ θ 0 ) , lg ρ lg k ]
(5) Using Equation (8), the phase correlation method can be used to calculate the rotation angle θ 0 and scaling parameter k .
Secondly, considering that in fast Fourier transform, if the image size N is the power of 2 (N = 2p, p is an integer), its calculation efficiency is highest, sub-images (Figure 5), which are cut from the original image with the image size of 2-power, can be used to calculate the transformation parameters first. Then, the translation offset for the original images can be calculated according to the location of sub-images. There are two overlap patterns including the left–right overlap as shown in Figure 5a and the up-down overlap as shown in Figure 5b. Using the 2-power sub-image, the image size becomes smaller, the overlapped rate is increased and the efficiency can be promoted.
Thirdly, when the image has many distortions or noises, the image registration result may be inaccurate [51]. Since image edges are the most basic and relatively stable feature, edges can be used to promote the efficiency and accuracy of the phase correlation method.
The steps of the improved phase correlation (ICP) method includes:
  • The 2-power sub-images of the reference image F 1 ( x , y ) and the registered image F 2 ( x , y ) are set as f 1 ( x , y ) and f 2 ( x , y ) respectively. The Canny edge detector is used to extract the edge of f 1 ( x , y ) and f 2 ( x , y ) , which are represented as f 1 ( x , y ) and f 2 ( x , y ) .
  • The Fourier transformation is implemented to f 1 ( x , y ) and f 2 ( x , y ) respectively, and then the amplitude spectra are obtained as M 1 ( u , v ) and M 2 ( u , v ) . The amplitude spectra is converted from the Cartesian coordinate system to the polar coordinate system, and M 1 ( θ , lg ρ ) and M 2 ( θ , lg ρ ) were got by logarithmic operation.
  • The traditional phase correlation method is applied to M 1 ( θ , lg ρ ) and M 2 ( θ , lg ρ ) to get the rotation angle θ 0 and the scale parameter k . Then, the θ 0 and k are used to do the inverse transformation to the 2-power sub-image of the registered image f 2 ( x , y ) to get a transition image f 3 ( x , y ) with only translation offsets.
  • Canny edge detector is used to get the edge image of f 3 ( x , y ) as f 3 ( x , y ) . Then, the phase correction method is used to calculate the relative translation offset ( x 0 , y 0 ) between the registered image F 2 ( x , y ) and f 3 ( x , y ) .
  • Based on the calculated θ 0 , k and ( x 0 , y 0 ) , the spatial similarity transformation is used and incorporated in a bilinear interpolation resampling procedure is used in image registration to produce the geometrically rectified hyperspectral image strips for the latter seamless mosaicking.

2.3.3. Image Fusion

Hyperspectral images are wete used for water monitoring, forestry analysis, mineral exploration, etc. For these applications, spectra fidelity of mosaicked images must be ensured. To eliminate the stitching line [52], the weighted average [53] and the best stitch fusion methods [54] were used. To analyze the spectral fidelity of the mosaicked image, the spectral reflectance curves of the different types of ground objects’ typical pixel in the overlapped area on the original, the geometrically corrected and the mosaicked images were compared. The spectral angle cosine (SAC), spectral correlation (SC), spectral information divergence (SID) and Euclidean distance (ED) [55] were used to measure the similarity of spectral curves in the overlapped area before and after image fusion. The greater the SAC value and the SC value, the greater the spectral similarity is, and the smaller the SID value and the ED value, the greater the spectral similarity is.

3. Experimental Results

As shown in Figure 4, different methods can be chosen in geometrical rectification, image registration and image fusion. In the following, the experiment results of different methods are presented. First, because few distinct GCPs can be extracted in the reference image in the river course area, the river course area is chosen as the representative study area to show the results of the different geometric rectification methods with or without GCPs. Secondly, the results of the two image registration methods of all the three study areas are presented. Thirdly, because the urban scape area consists of many types of ground object with different spectral feature characteristics, the two image fusion strategies are compared using the urban scape area as the representative study area. Moreover, the weighted average fusion results of the river course and forest areas are presented.

3.1. Geometric Rectification Results of River Course Area

Because a large portion of the river course area is dominated by water, the ground control points were derived from google earth (GE) images rather than from field surveys in the research. Eleven GCPs shown in Figure 6a were used for geometrical rectification.
Ten other GCPs in the GE image were used to assess the geometrical accuracy. As shown in Table 3, geometric rectification with POS data and GCPs is a compromise between efficiency and accuracy, and improves the accuracy compared to geometric rectification using only POS data. The medium error of the rectified image was 1.76 m, while the spatial resolution of the GE image was 2.5 m. The geometrical rectification error was less than 1 pixel, in accordance with the reference image.

3.2. Image Registration Results

As shown in Figure 7 and Figure 8, the ICP image registration method can get comparative results in the urban scape and river course areas; however, it failed in the forest area.
To compare the image registration methods, the testing points derived from the GE image for image rectification accuracy assessment were used to evaluate the accuracy for each study area after image registration. The time cost and accuracy are first analyzed for the river course and urban scape areas. The computer setting is as follows: Intel(R) Xeon (R) CPU E5-2630 v4 @ 2.20 GHz 2.20 GHz, the memory storage is 32G.
As shown in Table 4 and Table 5, for the urban scape and river course areas, the SIFT and RANSAC methods cost more time to establish the transformation model with a similar accuracy. The IPC method first extracts the edge image from the initial rectified image‘s 2-power sub-image, and then used Fourier transformation to convert the edge image from the spatial space to the frequency space; the transform parameters are calculated by the traditional phase correlation method, and finally the transform parameters of the original images are calculated. The time cost is much less than the traditional SIFT- and RANSAC-based methods.
For the forest data, the SIFT and RANSAC-based methods get good image registration results, with the medium error of 1.9669m, less than one pixel of the GE image, as shown in Table 6. However, the ICP method failed, with no forest result in Figure 8. For the spectral complex environment such as forest, the high spatial and spectral resolution image was characterized by the spectral heterogeneity and texture, and led to the difficulty of extracting accurate edges from the images, so the IPC method cannot get the accurate spatial similarity transformation parameters that are needed to co-register the neighboring image strips.

3.3. Image Fusion Results

As shown in Figure 7 and Figure 8, the spectral differences between the left and right geo-registered images are obvious. To get seamless image mosaicking results, the weighted average and the best stitching image fusion methods were used and compared in experiments.

3.3.1. Urban Scape Image Fusion Results

The weighted average fusion result and the best stitching fusion result of the urban scape areas are shown in Figure 9 and Figure 10. As shown in Figure 9, for images with complex features, the weighted average fusion is not suitable, and the red house has obvious ghosting problems. Since the best stitching line fusion takes into account the color and structure of the image, the stitching line avoids the ground objects with large differences in color and geometric structure, such as passing the houses, so the ghosting problem can be avoided and a better fusion effect is achieved for the urban scape area in Figure 10.
To further evaluate the best stitching fusion method for the urban scape area, the five typical ground objects including buildings, water, soil, farmland and roads in the overlapped area between the left and right image strips are compared with their spectral curves between the original, geo-registered and the mosaicked hyperspectral images (Figure 11). The spectral curve was extracted from one typical pixel of each type of ground object in the overlapped area between the neighboring image strips. From Figure 11, it can be seen that the difference between the spectra of the left and right original image strips is obvious especially for soil, farmland and road ground objects. These are mainly due to the following reasons.
First, the illumination difference of the round-trip flight strips; secondly, only one white board is set in the experimental area, and there may exist errors in the radiation calibration of images. Moreover, as the original images are re-sampled after the geometric rectification and image registration, the change in the original spectrum makes the left and right rectified images have larger spectral reflectance differences, especially for soil (600–700 nm), with a maximum difference of about 0.125, farmland (750–850 nm) with a maximum difference about 0.008, and road (550–850 nm) with a maximum difference about 0.03. After image fusion, the spectral curves of different ground objects (red curves) on the mosaic image is closer to that on the left registered (Rec_left) image (green curves) so the stitching line will disappear after image fusion, as shown in Figure 10. To find out whether the spectral reflectance before and after image mosaicking changes a lot, the spectral similarity assessment indices were calculated between the spectral curve of the mosaicked image and the original left and right images and the Rec_left image and Rec_right image as shown in Table 7. As can be seen from Table 7, the typical pixel of each type of ground object extracted from the mosaicked image has similar spectral reflectance to the Original left and Rec_ left images, which indicates that the spectral similarity of the pixels of the overlapped area before and after the mosaicking is high. Image fusion can reduce the spectral difference of the overlapped area’s ground object.

3.3.2. The Weighted Average Fusion Results of River Course and Forest Area

To further evaluate the image fusion methods, the river course and forest images before and after the weighted average fusion are shown in Figure 12 and Figure 13. As shown in the red rectangle on Figure 12 and Figure 13a, the spectral differences that are obvious in Figure 7c and Figure 8b disappear. Figure 13b shows the mosaicked image of multiple image strips of the mangrove forest. The mosaicked images have a good spectral smoothing quality in the overlapped area. It can be seen that the weighted average fusion method can get good results for the spectral homogeneous environments, such as river course and forest areas. With the consideration that the weighted average fusion method is more convenient than the best stitching fusion method, it is suggested to use the weighted average fusion method for the water and vegetation-dominated areas when the data volume is huge. However, when the computation ability is essential, the best stitching fusion method can perform better than the weighted average fusion method.

4. Discussion

UAV-based image products are generally produced by stitching together hundreds of partially overlapping image strips or frames captured in certain geometrical patterns [56]. However, when “matching” any two images, the transformation and reprojection using various algorithms routinely introduce localized distortions. While the positional accuracy of an individual image may be on the order of a few centimeters, the accuracy of the completed mosaic may increase to the decimeter range as the positional errors accumulate through the merging process. Accurate overlapping of swaths requires a good number of matching GCPs to avoid further distortion in the final mosaic [30]. To improve the geometrical accuracy, some researchers used UAV orthomaps as the reference data in image mosaicking. The identification of tie points in high resolution imagery, especially those exhibiting different geometric and spectral characteristics, is mainly established by applying a feature-based approach such as the modified SURF algorithm [28]. This method needs an accurate reference image in advance and much computational support for improving the automation and efficiency of its processing procedures [17]. Some other researchers used multispectral UAV imagery to obtain Ortho-mosaicked NDVI images for precision agriculture [57]. They used a common RGB camera and a modified camera sensitive to near-infrared radiation to gather images. The OpenDroneMap (ODM) was used to merge the two sets of images into ortho-mosaics. This work provided an efficient way to acquire large-area data for agriculture research using UAV. This method used the common aerial photogrammetry technique that required high image overlapping rates.
In comparison, our work is focused on the push-broom hyperspectral image mosaicking between image strips. The proposed chain of image mosaicking mainly includes three steps of geometrical rectification, image registration and image fusion.

4.1. Geometric Rectification

For the study areas that are dominated by water or forest there doesn’t exist so many distinct feature points that can be easily collected as GCPs. The traditional image matching approaches cannot extract sufficient tie points. To improve the efficiency and accuracy of geometrical correction, the different strategies of image rectification and image registration were compared in experiments. First, as shown in the image rectification experiments of the river course, the medium error of the POS-based method is 4.9 m, POS alone achieved a low quality of UAV georeferencing, especially in push-broom spectrometers [18]. Navigation grade GPS sensors do not facilitate hyperspectral push-broom scanning on a UAV [58]. The use of ground control points (GCPs) is absolutely necessary for improving georeferencing. Recent advances in direct georeferencing and imaging technologies—like the GNSS and INS—enable precise mapping using a minimal number of GCPs [30]. The medium error of the CSSF method with GCPs is 1.76 m, while the medium error of the POS-based method with less GCPs is 2.37 m. Geometric rectification based on the combination of POS data and GCPs is a compromise between efficiency and accuracy, and it not only reduces the demand for a large number of GCPs, but also improves the accuracy compared to geometric rectification using only POS data. In the experiments, the CSSF method is used to calculate the coordinates of rectified images with GCPs. By setting appropriate parameters, the modeled curved surface passes through each known GCP on the rectified images, and it is useful in eliminating the local deformation of images around GCPs and will benefit subsequent image registration. In all, the geometric correction accuracy of the mosaicked image depends mainly on the accuracy and the number of GCPs. In the study, GCPs were manually extracted from the GE images, so the GCP location errors may impact the geometric rectification of hyperspectral images, and the medium errors of geo-rectification results are almost 2 m for different study areas. In future, with more precise GCPs and as the amount of GCPs becomes larger, the geometric accuracy of the mosaic will increase.

4.2. Image Registration

To select an optimal band from the hyperspectral image before image registration, the local variance method based on edge block elimination [41] was used to calculate SNR for data quality characterization assessment of each band with high efficiency and precision [59]. In the image registration, the band with the highest SNR was selected in calculating the transform parameters between neighboring images. The SIFT and RANSAC-based method and the improved phase correlation method were compared from the aspects of time cost and geometric accuracy.
In image registration, the traditional phase correlation method is improved by using a 2-power sub-image and the edge information in the calculation of transformation parameters. The image registration experiment results of three types of ground areas show that the improved phase correlation method can achieve comparable results to the SIFT and RANSAC-based image registration method for study areas with distinct features such as urban scape and river course areas with less time cost. However, for the forest area, the improved phase correlation method failed and the SIFT and RANSAC-based method is suitable. To improve the effectiveness of the improved phase correlation method in areas with complex textures such as forest, further work needs to be done to find ways to calculate transformation parameters correctly.

4.3. Image Fusion

Another contribution of our work is that the spectral fidelity of each kind of ground object on the image before and after image mosaicking is ensured. The two image fusion methods were compared by both qualitative and quantitative assessments. In study areas such as the forest and the river course, the spectra of ground objects in the overlapped areas of images are almost homogeneous, the weighted average fusion method is suggested due to its simplicity for calculation. However, the urban scape area consists of various land cover types, and its spectra are variable in the overlapped area. The best stitching fusion method is suggested because it can avoid the ghosting problem that occurs when the weighted average fusion method was used.

5. Conclusions

This paper proposed a simple and effective chain of methods to produce geometrically corrected and mosaicked images with high spectral fidelity from UAV push-broom hyperspectral images. The proposed workflow generates a mosaicked output progressively by three major steps, including image rectification, registration and fusion. In the image mosaicking, the design of each step should be considered according to the application purpose. In our research work, the spectral homogenous area such as the river course and forest area always needs GCPs to get higher geometric accuracy. However, in practice, it is hard to get GCPs in such places by a field survey. Using Google earth images or another source of images as reference data, the proposed image mosaicking chain integrates the POS information, GCPs and the tie points together to refine the geometric accuracy of the mosaic. Researchers that are interested in using the UAV-based push-broom images and multi-source images together in environment monitoring can implement the proposed methodology to get a co-registered image data set. The image mosaic is of a high spectral fidelity. While it is essential to acquire data that can preserve scientific rigor, the radiometric calibration of the mosaic is beyond the scope of our research. There still remains some further works that can be done in the future:
(1) In the geometric rectification, the influence of the number and distribution of GCPs on the correction accuracy needs to be further studied. Since GCPs also have inherent errors, how to develop a method that can take into account the uncertainty of GCPs remains an issue that needs further investigation. In POS-based geometric rectification, this paper assumes that the scene area is a flat plain, which is not true for some study areas. Methods are yet to be developed to apply POS-based geometric rectification to non-flat terrains.
(2) In selecting the optimal band for image registration, SNR is used. However, the band with a higher SNR may have less information. This will also affect the accuracy of subsequent image registration and the subsequent processing steps. Therefore, other band selection methods may be explored to find an optimal band with both a high SNR and a high information content.
(3) Differences in reflectance between tie points in the original image strips lead to stitching lines. We mainly studied image fusion methods to eliminate the spectral difference in the overlapped area between the image strips, but did not calibrate the spectral reflectance of hyperspectral images against ground-based spectral measurements.
(4) The number of bands of the hyperspectral images used in this study is large, and with the continuous improvement in image resolution, the amount of data will become even larger in the future. It is necessary to study fast mosaicking methods or use multiple parallel computing systems to improve the efficiency of data processing.

Author Contributions

Conceptualization, L.Y.; methodology, G.Z.; validation, X.X.; data resources, X.M. and W.G.; data curation X.X.; writing—original draft preparation, L.Y.; writing—review and editing, J.M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the State Scholarship Fund of China (No. 201906435005), Chinese Academy of Sciences Strategic Leading Science and Technology Project (Class A) (XDA13020506), National Natural Science Foundation of China (61405204), China’s National Key R&D Programme (2018YFB0504903, 2016YFB0501402), the High Resolution Remote Sensing, Surveying and Mapping Application Demonstration System (Phase II) (42-Y30B04-9001-19/21), the Key Program of Marine Economy Development Special Foundation of the Department of Natural Resources of Guangdong Province (GDNRC [2020]012) and in part by the Key Deployment of Chinese Academy of Sciences (Y9F0600Z2F) and by the Fundamental Research Funds for the Central Universities(2021YQDC06) of China university of mining and technology-Beijing.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goetz, A.F.H. Three decades of hyperspectral remote sensing of the Earth: A personal view. Remote Sens. Environ. 2009, 113, S5–S16. [Google Scholar] [CrossRef]
  2. Kutser, T.; Hedley, J.; Giardino, C.; Roelfsema, C.; Brando, V.E. Remote sensing of shallow waters—A 50 year retrospective and future directions. Remote Sens. Environ. 2020, 240, 111619. [Google Scholar] [CrossRef]
  3. Banerjee, B.P.; Raval, S.; Cullen, P.J. UAV-hyperspectralimaging of spectrally complex environments. Int. J. Remote Sens. 2020, 41, 4136–4159. [Google Scholar] [CrossRef]
  4. Avta, R.; Watanabe, T. Unmanned Aerial Vehicle: Applications in Agriculture and Environment; Springer Nature: Basingstoke, UK, 2020. [Google Scholar]
  5. Jang, G.; Kim, J.; Yu, J.-K.; Kim, H.-J.; Kim, Y.; Kim, D.-W.; Kim, K.-H.; Lee, C.W.; Chung, Y.S. Review: Cost-Effective Unmanned Aerial Vehicle (UAV) Platform for Field Plant Breeding Application. Remote Sens. 2020, 12, 998. [Google Scholar] [CrossRef] [Green Version]
  6. Zhang, S.M.; Zhao, G.X.; Lang, K.; Su, B.W.; Chen, X.N.; Xi, X.; Zhang, H.B. Integrated Satellite, Unmanned Aerial Vehicle (UAV) and Ground Inversion of the SPAD of Winter Wheat in the Reviving Stage. Sensors 2019, 19, 1485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Honkavaara, E.; Rosnell, T.; Oliveira, R.; Tommaselli, A. Band registration of tuneable frame format hyperspectral UAV imagers in complex scenes. ISPRS J. Photogramm. Remote Sens. 2017, 134, 96–109. [Google Scholar] [CrossRef]
  8. Wu, W.B.; Zhang, Z.B.; Zheng, L.J.; Han, C.Y.; Wang, X.M.; Xu, J.; Wang, X.R. Research Progress on the Early Monitoring of Pine Wilt Disease Using Hyperspectral Techniques. Sensors 2020, 20, 3729. [Google Scholar] [CrossRef]
  9. Zarco-Tejada, P.J.; González-Dugo, V.; Berni, J.A.J. Fluorescence, temperature and narrow-band indices acquired from a UAV platform for water stress detection using a micro-hyperspectral imager and a thermal camera. Remote Sens. Environ. 2012, 117, 322–337. [Google Scholar] [CrossRef]
  10. Stuart, M.B.; Mcgonigle, A.J.S.; Willmott, J.R. Hyperspectral Imaging in Environmental Monitoring: A Review of Recent Developments and Technological Advances in Compact Field Deployable Systems. Sensors 2019, 19, 3071. [Google Scholar] [CrossRef] [Green Version]
  11. Cao, J.; Leng, W.; Liu, K.; Liu, L.; He, Z.; Zhu, Y. Object-Based Mangrove Species Classification Using Unmanned Aerial Vehicle Hyperspectral Images and Digital Surface Models. Remote Sens. 2018, 10, 89. [Google Scholar] [CrossRef] [Green Version]
  12. Lu, B.; He, Y. Optimal spatial resolution of Unmanned Aerial Vehicle (UAV)-acquired imagery for species classification in a heterogeneous grassland ecosystem. GIScience Remote Sens. 2018, 55, 205–220. [Google Scholar] [CrossRef]
  13. Lu, B.; He, Y. Species classification using Unmanned Aerial Vehicle (UAV)-acquired high spatial resolution imagery in a heterogeneous grassland. ISPRS J. Photogramm. Remote Sens. 2017, 128, 73–85. [Google Scholar] [CrossRef]
  14. Lu, B.; He, Y.; Liu, H.H.T. Mapping vegetation biophysical and biochemical properties using unmanned aerial vehicles-acquired imagery. Int. J. Remote Sens. 2017, 39, 5265–5287. [Google Scholar] [CrossRef]
  15. Büttner, A.; Röser, H.-P. Hyperspectral remote sensing with the UAS ‘‘Stuttgarter Adler”—system setup, calibration and first results. Photogramm. Fernerkund. Geoinf. 2014, 4, 265–274. [Google Scholar] [CrossRef] [PubMed]
  16. Barreto, M.A.P.; Johansen, K.; Angel, Y.; McCabe, M.F. Radiometric Assessment of a UAV-Based Push-Broom Hyperspectral Camera. Sensors 2019, 19, 4699. [Google Scholar] [CrossRef] [Green Version]
  17. Angel, Y.; Turner, D.; Parkes, S.; Malbeteau, Y.; Lucieer, A.; McCabe, M.F. Automated Georectification and Mosaicking of UAV-Based Hyperspectral Imagery from Push-Broom Sensors. Remote Sens. 2020, 12, 34. [Google Scholar] [CrossRef] [Green Version]
  18. Habib, A.; Xiong, W.; He, F.; Yang, H.L.; Crawford, M. Improving Orthorectification of UAV-Based Push-Broom Scanner Imagery Using Derived Orthophotos From Frame Cameras. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 10, 262–276. [Google Scholar] [CrossRef]
  19. Yang, G.; Li, C.; Wang, Y.; Yuan, H.; Feng, H.; Xu, B.; Yang, X. The DOM Generation and Precise Radiometric Calibration of a UAV-Mounted Miniature Snapshot Hyperspectral Imager. Remote Sens. 2017, 9, 642. [Google Scholar] [CrossRef] [Green Version]
  20. 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]
  21. Hagen, N.; Kudenov, M.W. Review of snapshot spectral imaging technologies. Opt. Eng. 2013, 52, 090901. [Google Scholar] [CrossRef] [Green Version]
  22. Suomalainen, J.; Anders, N.; Iqbal, S.; Roerink, G.; Franke, J.; Wenting, P.; Hünniger, D.; Bartholomeus, H.; Becker, R.; Kooistra, L. A Lightweight Hyperspectral Mapping System and Photogrammetric Processing Chain for Unmanned Aerial Vehicles. Remote Sens. 2014, 6, 11013–11030. [Google Scholar] [CrossRef] [Green Version]
  23. Kim, J.-I.; Kim, T.; Shin, D.; Kim, S. Fast and robust geometric correction for mosaicking UAV images with narrow overlaps. Int. J. Remote Sens. 2017, 38, 2557–2576. [Google Scholar] [CrossRef]
  24. Faraji, M.R.; Qi, X.; Jensen, A. Computer vision–based orthorectification and georeferencing of aerial image sets. J. Appl. Remote Sens. 2016, 10, 036027. [Google Scholar] [CrossRef]
  25. Shen, X.; Cao, L.; Coops, N.C.; Fan, H.; Wu, X.; Liu, H.; Wang, G.; Cao, F. Quantifying vertical profiles of biochemical traits for forest plantation species using advanced remote sensing approaches. Remote Sens. Environ. 2020, 250, 112041. [Google Scholar] [CrossRef]
  26. Aasen, H.; Honkavaara, E.; Lucieer, A.; Zarco-Tejada, P.J. Quantitative Remote Sensing at Ultra-High Resolution with UAV Spectroscopy: A Review of Sensor Technology, Measurement Procedures, and Data Correction Workflows. Remote Sens. 2018, 10, 1091. [Google Scholar] [CrossRef] [Green Version]
  27. Koshy, G.; Vishnukumar, S. A hybrid approach to generate visually seamless aerial mosaicks from unmanned aerial vehicles. J. Intell. Fuzzy Syst. 2019, 36, 2075–2083. [Google Scholar]
  28. Habib, A.; Han, Y.; Xiong, W.; He, F.; Zhang, Z.; Crawford, M. Automated Ortho-Rectification of UAV-Based Hyperspectral Data over an Agricultural Field Using Frame RGB Imagery. Remote Sens. 2016, 8, 796. [Google Scholar] [CrossRef] [Green Version]
  29. Kirsch, M.; Lorenz, S.; Zimmermann, R.; Tusa, L.; Möckel, R.; Hödl, P.; Booysen, R.; Khodadadzadeh, M.; Gloaguen, R. Integration of Terrestrial and Drone-Borne Hyperspectral and Photogrammetric Sensing Methods for Exploration Mapping and Mining Monitoring. Remote Sens. 2018, 10, 1366. [Google Scholar] [CrossRef] [Green Version]
  30. Adão, T.; Hruška, J.; Pádua, L.; Bessa, J.; Peres, E.; Morais, R.; Sousa, J.J. Hyperspectral imaging: A review on UAV-based sensors, data processing and applications for agriculture and forestry. Remote Sens. 2017, 9, 1110. [Google Scholar] [CrossRef] [Green Version]
  31. Jaud, M.; Le Dantec, N.; Ammann, J.; Grandjean, P.; Constantin, D.; Akhtman, Y.; Barbieux, K.; Allemand, P.; Delacourt, C.; Merminod, B. Direct Georeferencing of a Pushbroom, Lightweight Hyperspectral System for Mini-UAV Applications. Remote Sens. 2018, 10, 204. [Google Scholar] [CrossRef] [Green Version]
  32. Li, X.Y. Principle, Method and Practice of IMU/DGPS Based Photogrammetry. Ph.D. Thesis, Information Engineering University, Zhengzhou, China, 2005. [Google Scholar]
  33. Yuan, X.X.; Zhang, X.P. Theoretical accuracy of direct georeferencing with position and orientation system in aerial photogrammetry. Int. Arch. Photogramm. Remote Sens. Spat. Sci. 2008, XXXVII Pt B1, 617–622. [Google Scholar]
  34. Moroni, M.; Dacquino, C.; Cenedese, A. Mosaicing of Hyperspectral Images: The Application of a Spectrograph Imaging Device. Sensors 2012, 12, 10228–10247. [Google Scholar] [CrossRef] [PubMed]
  35. Olsson, P.-O.; Vivekar, A.; Adler, K.; Garcia Millan, V.E.; Koc, A.; Alamrani, M.; Eklundh, L. Radiometric Correction of Multispectral UAS Images: Evaluating the Accuracy of the Parrot Sequoia Camera and Sunshine Sensor. Remote Sens. 2021, 13, 577. [Google Scholar] [CrossRef]
  36. Honkavaara, E.; Khoramshahi, E. Radiometric Correction of Close-Range Spectral Image Blocks Captured Using an Unmanned Aerial Vehicle with a Radiometric Block Adjustment. Remote Sens. 2018, 10, 256. [Google Scholar] [CrossRef] [Green Version]
  37. Flight Rules of Light-Weighted UAV, Civil Aviation Adminstration of China. 2015. Available online: http://www.caac.gov.cn/XXGK/XXGK/GFXWJ/201601/P020170527591647559640.pdf (accessed on 13 November 2021).
  38. Hruska, R.; Mitchell, J.; Anderson, M.; Glenn, N.F. Radiometric and Geometric Analysis of Hyperspectral Imagery Acquired from an Unmanned Aerial Vehicle. Remote Sens. 2012, 4, 2736–2752. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, S.G. The Principal and Applications on Photogrammetry, 1st ed.; Wuhan University Press: Wuhan, China, 2009. [Google Scholar]
  40. Harder, R.L.; Desmarais, R.N. Interpolation using surface splines. J. Aircr. 1972, 9, 189–191. [Google Scholar] [CrossRef]
  41. Zhang, L.; Huang, W.M.; Zhang, Y.; Xu, Y.M.; Zhou, C. Estimation of Signal-Noise-Ratio for HJ-1/CCD Data. Geospatial Inf. 2013, 3, 73–75. [Google Scholar]
  42. Tan, X.; Sun, C.; Sirault, X.; Furbank, R.; Pham, T.D. Feature matching in stereo images encouraging uniform spatial distribution. Pattern Recognit. 2015, 48, 2530–2542. [Google Scholar] [CrossRef]
  43. Abdel-Hakim, A.E.; Aly, A.F. Csift: A sift descriptor with color invariant characteristics. In Proceedings of the 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR′06), New York, NY, USA, 17–22 June 2006. [Google Scholar]
  44. Jakubovic, A.; Velagic, J. Image Feature Matching and Object Detection Using Brute-Force Matchers. In Proceedings of the 2018 International Symposium ELMAR, Zadar, Croatia, 16–19 September 2018. [Google Scholar]
  45. Li, T.T.; Jiang, B.; Tu, Z.Z.; Luo, B.; Tang, J. Image Matching Using Mutual k-Nearest Neighbor Graph. In Intelligent Computation in Big Data Era; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  46. Liu, H.; Deng, M.; Xiao, C.; Deng, M. An improved best bin first algorithm for fast image registration. In Proceedings of the 2011 International Conference on Electronic & Mechanical Engineering and Information Technology, Harbin, China, 12–14 August 2011; Volume 1, pp. 355–358. [Google Scholar]
  47. Bellavia, F.; Colombo, C. Is There Anything New to Say About SIFT Matching? Int. J. Comput. Vis. 2020, 128, 1847–1866. [Google Scholar] [CrossRef] [Green Version]
  48. Li, J.; Wang, H.; Zhang, L.; Wang, Z.; Wang, M. The Research of Random Sample Consensus Matching Algorithm in PCA-SIFT Stereo Matching Method. In Proceedings of the 2019 Chinese Control And Decision Conference (CCDC), Nanchang, China, 3–5 June 2019; pp. 3338–3341. [Google Scholar]
  49. Zhou, M.L.; Bai, Z.W.; Yan, X.J. Design of the phase correction based image stitching system. Foreign Electron. Meas. Technol. 2015, 5, 31–33. (In Chinese) [Google Scholar]
  50. Reddy, B.; Chatterji, B. An FFT-based technique for translation, rotation, and scale-invariant image registration. IEEE Trans. Image Process. 1996, 5, 1266–1271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Isoardi, R.A.; Osorio, A.R.; Mato, G. Medical Image Registration with Fourier basis Functions; INTECH Open Access Publisher: Metro Manila, Philippines, 2001. [Google Scholar]
  52. Szeliski, R. Image Alignment and Stitching: A Tutorial. Found. Trends® Comput. Graph. Vis. 2007, 2, 1–104. [Google Scholar] [CrossRef]
  53. Luo, Y.T.; Wang, Y.; Zhang, H.M. Image-stitching Algorithm by Combining the Optimal Seam and an Improved Gradual Fusion Method. Infrared Technol. 2018, 40, 382–387. [Google Scholar]
  54. Gu, Y.; Zhou, G.; Ren, G. Image stitching by combining optimal seam and multi-resolution fusion. J. Image Graphics. 2017, 22, 842–851. [Google Scholar]
  55. Zhang, J.Z.; Zhu, W.Q.; Zheng, Z.T. Research on the similarity of hyperspectral image. Sci. Surv. Mapp. 2013, 38, 33–36. [Google Scholar]
  56. Jakob, S.; Zimmermann, R.; Gloaguen, R. The Need for Accurate Geometric and Radiometric Corrections of Drone-Borne Hyperspectral Data for Mineral Exploration: MEPHySTo—A Toolbox for Pre-Processing Drone-Borne Hyperspectral Data. Remote Sens. 2017, 9, 88. [Google Scholar] [CrossRef] [Green Version]
  57. de Oca, A.M.; Arreola, L.; Flores, A.; Sanchez, J.; Flores, G. Low-cost multispectral imaging system for crop monitoring. In Proceedings of the 2018 International Conference on Unmanned Aircraft Systems (ICUAS), Dallas, TX, USA, 12–15 June 2018; pp. 443–451. [Google Scholar]
  58. Lucieer, A.; Malenovský, Z.; Veness, T.; Wallace, L. HyperUAS-Imaging Spectroscopy from a Multirotor Unmanned Aircraft System: HyperUAS-Imaging Spectroscopy from a Multirotor Unmanned. J. Field Robot 2014, 31, 571–590. [Google Scholar] [CrossRef] [Green Version]
  59. Richter, R.; Schläpfer, D.; Muller, A. Operational Atmospheric Correction for Imaging Spectrometers Accounting for the Smile Effect. IEEE Trans. Geosci. Remote Sens. 2011, 49, 1772–1780. [Google Scholar] [CrossRef]
Figure 1. UAV hyperspectral system used in the study: (a) M600 Pro UAV platform with the hyperspectral imaging system onboard; (b) The ZK-VNIR-FPG480 hyperspectral imaging spectrometry system.
Figure 1. UAV hyperspectral system used in the study: (a) M600 Pro UAV platform with the hyperspectral imaging system onboard; (b) The ZK-VNIR-FPG480 hyperspectral imaging spectrometry system.
Remotesensing 13 04720 g001
Figure 2. Location of the study areas in China.
Figure 2. Location of the study areas in China.
Remotesensing 13 04720 g002
Figure 3. True color images of two neighboring image strips (bands of 103, 63 and 27): (a) Urban scape; (b) River course; (c) Forest.
Figure 3. True color images of two neighboring image strips (bands of 103, 63 and 27): (a) Urban scape; (b) River course; (c) Forest.
Remotesensing 13 04720 g003
Figure 4. The image mosaicking chain.
Figure 4. The image mosaicking chain.
Remotesensing 13 04720 g004
Figure 5. Two-power sub-images (the blue and red line represent the boundaries of 2-power sub-images): (a) the left-right overlap 2-power sub-images; (b) the up-down overlap 2-power sub-images.
Figure 5. Two-power sub-images (the blue and red line represent the boundaries of 2-power sub-images): (a) the left-right overlap 2-power sub-images; (b) the up-down overlap 2-power sub-images.
Remotesensing 13 04720 g005
Figure 6. River course GCPs and geometric rectified results: (a) GCPs on GE image; (b) CSSF method with 11 GCPs; (c) POS-based method without GCP; (d) POS-based geometric rectification method with 4 GCPs that labeled as 2, 6, 8 and 11 on Figure 5a.
Figure 6. River course GCPs and geometric rectified results: (a) GCPs on GE image; (b) CSSF method with 11 GCPs; (c) POS-based method without GCP; (d) POS-based geometric rectification method with 4 GCPs that labeled as 2, 6, 8 and 11 on Figure 5a.
Remotesensing 13 04720 g006
Figure 7. SIFT and RANSAC-based image registration results: (a) Urban scape; (b) River course; (c) Forest.
Figure 7. SIFT and RANSAC-based image registration results: (a) Urban scape; (b) River course; (c) Forest.
Remotesensing 13 04720 g007
Figure 8. Image registration results using the improved phase correlation method: (a) Urban scape; (b) River course.
Figure 8. Image registration results using the improved phase correlation method: (a) Urban scape; (b) River course.
Remotesensing 13 04720 g008
Figure 9. The weighted average fusion result of a urban scape area.
Figure 9. The weighted average fusion result of a urban scape area.
Remotesensing 13 04720 g009
Figure 10. The best stitching image fusion result of an urban scape area.
Figure 10. The best stitching image fusion result of an urban scape area.
Remotesensing 13 04720 g010
Figure 11. Comparison of typical spectral curves before and after image mosaicking for the urban scape area: the spectral curves of objects of (a) building; (b) water; (c) soil; (d) farmland; (e) road are represented in five colors, the sources of the spectral curves are shown in (f) using the texts in the corresponding colors, where the original_left means the original left image, original_right means the original right image; mosaic means the mosaicked image; Rec_left means the geo-registered left image, Rec_right means the geo-registered right image.
Figure 11. Comparison of typical spectral curves before and after image mosaicking for the urban scape area: the spectral curves of objects of (a) building; (b) water; (c) soil; (d) farmland; (e) road are represented in five colors, the sources of the spectral curves are shown in (f) using the texts in the corresponding colors, where the original_left means the original left image, original_right means the original right image; mosaic means the mosaicked image; Rec_left means the geo-registered left image, Rec_right means the geo-registered right image.
Remotesensing 13 04720 g011
Figure 12. The weighted average fusion result of river course area.
Figure 12. The weighted average fusion result of river course area.
Remotesensing 13 04720 g012
Figure 13. Forest experiment results: (a) the fused image by the weighted average method; and (b) the fused image of multiple forest image strips.
Figure 13. Forest experiment results: (a) the fused image by the weighted average method; and (b) the fused image of multiple forest image strips.
Remotesensing 13 04720 g013
Table 1. Technical indicators of ZK-VNIR-FPG480 hyperspectral imager.
Table 1. Technical indicators of ZK-VNIR-FPG480 hyperspectral imager.
ParameterIndicators
Spectral range400 nm~1000 nm
Number of spectral channels270
Spectral resolution2.8 nm
Line array width 480
Spatial resolution9 cm at 100 m
Field of View26° at 35 mm (lens focal length)
A/D conversion12 bits
Max framerate100 fps
Table 2. Flight parameters of different study areas.
Table 2. Flight parameters of different study areas.
Study AreaFlight HeightSide Overlap Rate
Urban scape200 m30%
River course120 m50%
Forest90 m20%
Table 3. Geometric rectification accuracy (units in meters) of the river course area.
Table 3. Geometric rectification accuracy (units in meters) of the river course area.
Mean Residual Medium ErrorMean Absolute DeviationMedian Absolute DeviationStandard DeviationMaximum Residual
CSSF method with 11 GCPs x 0.08620.69540.60680.86210.69001.0909
y −0.11981.62071.24742.23331.61632.8232
x o y 1.7636 1.7574
POS-based method without GCPs x −4.13544.35931.15551.17441.37926.8187
y −0.09612.24692.10941.90822.24493.8767
x o y 4.9043 2.6347
POS-based method with 4 GCPs x −1.25531.74131.0260.50311.20683.372
y 0.21621.62081.36271.03671.60633.1469
x o y 2.3789 2.0091
Table 4. The image registration accuracy—urban scape (units in meters).
Table 4. The image registration accuracy—urban scape (units in meters).
Mean Residual Medium ErrorMean Absolute DeviationMedian Absolute DeviationStandard DeviationMaximum ResidualTime Cost
The improved phase correlation method x −0.670251.16641.59920.68991.10193.54980.329s
y 0.247541.08131.3071.55170.86183.4393
x o y 1.5905 1.3989
The SIFT and RANSAC-based method x −1.67311.21040.956942.84781.13263.06295.395s
y −0.821381.15520.971381.07080.81372.7707
x o y 1.6731 1.3945
Table 5. The image registration accuracy—river course (units in meters).
Table 5. The image registration accuracy—river course (units in meters).
Mean Residual Medium ErrorMean Absolute DeviationMedian Absolute DeviationStandard DeviationMaximum ResidualTime Cost
The improved phase correlation method x −2.79751.48561.75810.956151.60693.77091.169s
y −0.954351.24700.731390.412350.950642.8895
1.9395 1.8670
The SIFT and RANSAC-based method x −1.15301.21911.75201.01801.59614.22474.082s
y −0.974211.26130.636850.195850.801092.5976
x o y 1.7541 1.7858
Table 6. The SIFT and RANSAC-based image registration accuracy—forest (units in meters).
Table 6. The SIFT and RANSAC-based image registration accuracy—forest (units in meters).
Mean Residual Medium ErrorMean Absolute DeviationMedian Absolute DeviationStandard DeviationMaximum Residual
x 0.08640.90650.79771.40190.90241.3841
y 0.21661.74551.37502.62051.73213.0199
x o y 1.9669 1.9530
Table 7. Evaluation of similarity of spectral curves of different ground objects in the overlapped area before and after image mosaicking (urban scape).
Table 7. Evaluation of similarity of spectral curves of different ground objects in the overlapped area before and after image mosaicking (urban scape).
Ground ObjectsThe Compared ImageSACSCSIDED
buildingRec_left image0.99990.99980.000020.0377
Rec_right image0.99970.99840.000220.1114
Original left image0.99990.99970.000040.0407
Original right image0.99980.99900.000230.1313
waterRec_left image0.99890.99640.000950.03670
Rec_right image0.99340.98630.004790.05112
Original left image0.99230.98940.005480.05837
Original right image0.98570.97760.010660.07115
soilRec_left image0.99510.99210.003960.3190
Rec_right image0.99680.99690.006700.3277
Original left image0.99710.99300.007450.4289
Original right image0.99920.99680.001430.3315
farmlandRec_left image0.99910.99840.002090.1664
Rec_right image0.99970.99950.000710.1927
Original left image0.99860.99730.002310.2403
Original right image0.99950.99870.000870.2628
roadRec_left image0.99990.99920.000010.0457
Rec_right image0.99860.91110.001050.1207
Original left image0.99860.95740.001220.2362
Original right image0.99760.90910.002170.1511
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yi, L.; Chen, J.M.; Zhang, G.; Xu, X.; Ming, X.; Guo, W. Seamless Mosaicking of UAV-Based Push-Broom Hyperspectral Images for Environment Monitoring. Remote Sens. 2021, 13, 4720. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224720

AMA Style

Yi L, Chen JM, Zhang G, Xu X, Ming X, Guo W. Seamless Mosaicking of UAV-Based Push-Broom Hyperspectral Images for Environment Monitoring. Remote Sensing. 2021; 13(22):4720. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224720

Chicago/Turabian Style

Yi, Lina, Jing M. Chen, Guifeng Zhang, Xiao Xu, Xing Ming, and Wenji Guo. 2021. "Seamless Mosaicking of UAV-Based Push-Broom Hyperspectral Images for Environment Monitoring" Remote Sensing 13, no. 22: 4720. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224720

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