Geometric Quality Analysis of AVHRR Orthoimages
1
Department of Geomatics Engineering, Hacettepe University, Ankara 06800, Turkey
2
Department of Geography, University of Bern, CH-3012 Bern, Switzerland
3
ETH Zurich, Institute of Geodesy and Photogrammetry, CH-8093 Zurich, Switzerland
*
Author to whom correspondence should be addressed.
†
The work has been done while she was at ETH Zurich.
Remote Sens. 2015, 7(3), 3293-3319; https://doi.org/10.3390/rs70303293
Submission received: 15 September 2014
/
Revised: 6 March 2015
/
Accepted: 9 March 2015
/
Published: 23 March 2015
Abstract
:The geometric accuracy of 2008 AVHRR orthoimages from Metop-A, NOAA-17 and NOAA-18 satellites over Switzerland have been investigated here. The methods employed in the study are fully automated, with an accuracy of 0.1–0.2 pixels, however, blunders do occur and this requests a careful blunder detection approach. The investigations include analysis of relative, absolute and band-to-band registration (BBR) accuracy. Regarding relative accuracy, thousands of points are matched between Metop-A, NOAA-17 and NOAA-18 images of the same day. The accuracy is quite high with mean shifts between 0.2 and 0.4 pixels. Systematic stripes have been observed when NOAA-18 images are involved in matching. In spite of many efforts to find the source of this error, no explanation could be found. In addition, large shifts up to 2.9 pixels on some days between September and December 2008 were observed. Regarding absolute accuracy, digitized lakes as reference polygons have been used and a subpixel lake matching method has been applied. The mean shifts generally fulfilled EUMETSAT and GCOS specifications, although some partial results exceed them, especially for Metop-A. Regarding BBR accuracy, six multispectral bands have been compared, also with point matching. The EUMETSAT specification is 0.1 km, however, this specification refers to original images, not orthoimages. Taking also into account the matching errors of 0.1 km, the EUMETSAT specifications are in principle fulfilled in all cases except matching of Metop-A and NOAA-17 Band-2 images with Bands 4 and 5. The overall work showed that although, in general, accuracies are high and fulfill specifications, errors exceeding the specifications can occur and vary depending on the satellite used, time and location. Such errors influence subsequent geometric or thematic processing; thus, an automated and permanent quality control of such images should be executed.
1. Introduction
1.1. Motivation and Background
A research project has been carried out at ETH Zurich in cooperation with the Swiss GCOS (Global Climate Observing System) office at MeteoSwiss, which coordinates the GCOS activities in Switzerland [1,2]. The GCOS is a long-term program aiming at providing systematic observations required for monitoring and improving the understanding of the climate system, detecting climate changes, and assessing their impacts and sponsored by WMO (World Meteorological Organization), the Intergovernmental Oceanographic Commission (IOC) of the UNESCO, the UNEP (UN Environmental Programme), and the ICSU (International Council of Science) [3]. GCOS has defined the Essential Climate Variables (ECVs) and the target requirements for their geometric accuracy [4]. The ECVs for which satellite observations make a significant contribution are listed by WMO [4].
Possible geometric errors influence the subsequent extraction of both geometric and thematic information from satellite images, and in this project especially the location of ECVs. Even small geometric errors can have a significant influence on retrieval of thematic information (e.g., see [5]). The project aims at contributing to the GCOS objectives by analyzing the geometric accuracy of three commonly used satellite sensors especially for meteorological applications: SEVIRI (Spinning Enhanced Visible and Infrared Imager on board of the European meteorological satellites of the Meteosat Second Generation (MSG) Series), MODIS (Moderate Resolution Imaging Spectroradiometer on board the NASA satellites Terra and Aqua), and AVHRR (Advanced Very High Resolution Radiometer on board the NOAA and Metop satellite series). The sensor data and the product level (in case of AVHRR, orthoimages) were decided by the project initiator, MeteoSwiss, as they are routinely used by them for further processing. Due to time and reference data acquisition constraints, the project area is limited to Switzerland and surroundings.
The main aim of this paper is to present the investigations of the relative, absolute and band-to-band registration (BBR) accuracies of the AVHRR orthoimages and to provide a comparison to the product specifications. The relative accuracy refers to the geometric accuracy and stability between two images (usually acquired with small time difference, and often from the same sensor) and also consistency of the accuracy within an image. In previous MSG tests, daily, consecutive images of SEVIRI with 15 min time difference have been compared. Since such timewise dense data do not exist for AVHRR, daily AVHRR orthoimages (usually three) from three satellites have been compared. The absolute accuracy refers to the coordinate accuracy of image points or features with respect to their ground coordinates, i.e., absolute georeferencing (geolocation). The BBR accuracy refers to geometric registration errors between different, simultaneously acquired spectral bands of a sensor. In all three accuracy investigations, corresponding features were matched using two shifts. One reason to use only shifts is the fact that the data used here are orthoimages; thus, scale, rotation and other higher-order deformations between the images are negligible. Furthermore, since the number of features matched is 10–20 for absolute accuracy, and several hundreds for relative and BBR accuracy, higher order transformations can be estimated from the local feature shifts. The results have been compared with the accuracy specifications given by EUMETSAT and GCOS (latter for map products). The specifications of EUMETSAT [6] for absolute and BBR accuracy are 1 km and 0.1 km, respectively. Thereby, it is not clarified whether they refer to radial distance or tolerances for each of the involved coordinates. The BBR specifications refer to the original spectral bands, not orthoimages as used in this case, where additional errors due to the used DSM (Digital Surface Model) and sensor geometric calibration and orientation may influence the product accuracy (however, the orthoimages used here have the same DSM, geometric calibration and orientation errors, thus these tests are equivalent to testing original images). The GCOS specifications for map products (including orthoimages) are 1/3 pixel.
The methodology was initially developed and tested with the MSG-SEVIRI Level 1.5 images acquired from Meteosat-8 and Meteosat-9. The preliminary results of MSG-SEVIRI investigations in terms of relative and absolute accuracy were published in [7]. More comprehensive analysis on relative accuracy was given in [8] and the BBR accuracy evaluation results were presented in [9]. The tests with the MSG-SEVIRI images have shown that the accuracy of the point matching method (with Kanade-Lukas-Tomasi (KLT) tracking) obtained in the tests is about 0.1 pixels and for lake matching about 0.2 pixels. The methods have been adapted for the AVHRR orthoimages. The analyses were done using data acquired in 2008 over Switzerland and surroundings.
1.2. AVHRR Sensor and Data Characteristics
The Advanced Very High Resolution Radiometer/3 (AVHRR/3) is a multipurpose imaging instrument used for global monitoring of cloud cover, sea surface temperature, ice, snow and vegetation cover characteristics and is/was aboard the NOAA satellite series and Metop-A [6]. The AVHRR operates at six different bands simultaneously, but only five bands are transmitted to ground. It provides cloud and surface information with about 1.1 km resolution near nadir [6] and about 4 km at maximum off-nadir angle. The band information of AVHRR is provided in Table 1.
AVHRR orthoimages with 1 km spatial resolution have been investigated in this project. The orthorectification procedure has been developed by Khlopenkov et al. [10] and accuracy better than 1/3 GSD (Ground Sampling Distance) has been achieved in their tests with historical AVHRR images over Canada. This current study employed the AVHRR imagery processed with an augmented system (software package SAPS (Science Systems and Applications, Inc. AVHRR Processing System)) which was developed for the region of Europe including Switzerland based on principles described in [10]. The data have been acquired from three different polar-orbiting satellites, NOAA-17, NOAA-18 and Metop-A. The raw AVHRR images have been georeferenced at the University of Bern. Processing system employed 250 m GSD MODIS monthly composites and feature matching between MODIS and AVHRR images to acquire GCPs over the region of Europe. The images were then orthorectified using a 500 m grid spacing SRTM (Shuttle Radar Topography Mission) DSM, with voids filled by other DSM sources. The orthoimages are produced with 1 km GSD in Lambert Conical Projection (E0: −1700000.0, N0: 8600000.0, false origin at: 0°/10°, standard parallels at: 40°/55°) on WGS84 ellipsoid. The same georeferencing information is used for all six bands of an image acquisition, while the DSM used is the same for all data used here. Visual checks were performed on the accuracy of the orthoimages by the University of Bern and provided very stable and accurate impressions except for the very off-nadir acquisitions (personal communications). Cloud masks were also produced by the University of Bern and the methodology was under testing at the time of data delivery.
Table 1.
AVHRR/3 band characteristics [11].
| Band Number | Resolution at Nadir (km) | Wavelength (µm) | Typical Use |
|---|---|---|---|
| 1 | 1.09 | 0.58–0.68 | Daytime cloud and surface mapping |
| 2 | 1.09 | 0.725–1.00 | Land-water boundaries |
| 3A | 1.09 | 1.58–1.64 | Snow and ice detection |
| 3B | 1.09 | 3.55–3.93 | Night cloud mapping, sea surface temperature |
| 4 | 1.09 | 10.30–11.30 | Night cloud mapping, sea surface temperature |
| 5 | 1.09 | 11.50–12.50 | Sea surface temperature |
The data from 24 days—two days from each month in 2008—has been received from the University of Bern together with the cloud masks. The days are selected according to minimum cloud coverage. Since the main aim is to investigate the geometric accuracy over Switzerland, only parts of the orthoimages have been used for the tests (will mainly be referred as images from here on). Mainly, three close to nadir acquisitions (one per satellite) over Switzerland are available per day. However, more images have been provided for some days where the project area is covered by off-nadir images (close to the image strip borders). Although the preliminary investigations were performed including off-nadir images (the part of the AVHRR image corresponding to a large viewing zenith angle of about >45°), the very off-nadir images have been excluded from the final analysis due to much larger GSD and image blurring.
2. Methodology
The methods have been previously developed for the relative, absolute, and BBR accuracy investigations of MSG-SEVIRI images and implemented as standalone software using Python programming language. These methods have been described in detail in [7,8,9] and have been adapted for the AVHRR data in this study. The major modifications applied to the processes are described in the following subsections. The tests with the MSG-SEVIRI images have shown that the accuracy of the point matching method (with KLT tracking) obtained in the tests and excluding blunders is about 0.1 pixels and for lake matching about 0.2 pixels. The image preprocessing, cloud conversion, the KLT and the lake matching, and the statistical analysis and outlier detection methods run automatically for all images without user interaction.
The KLT tracking method used in relative and BBR accuracy extracts interest points and matches them (estimating two shifts) using least squares matching, by matching two image patches [12]. The KLT Tracking function of OpenCV) [13] has been employed for image matching and a statistical outlier detection has been used iteratively. A point is excluded if Δx or Δy > 3σ, with Δx or Δy, the difference between x-, y-shift of each point and the mean x-, y-shift for all image points. For BBR accuracy, a stricter threshold of 2σ was used. For relative and BBR accuracy, an additional outlier detection method using cross-correlation was employed, e.g. for relative accuracy, if after matching the cross-correlation value is <0.8, the match point is rejected. For BBR, a stricter threshold of 0.9 was used.
The results for each image pair are evaluated via statistical analysis of point/lake shifts (i.e., minimum, maximum, mean, standard deviation, median (MED) and median absolute deviation from the median (MAD)) and additional visual checks of the shift plots. MAD is a robust estimator of standard deviation, when multiplied by a constant, which for normal distribution is 1.4826. The MAD values given in tables are without this multiplication and this should be taken into account for comparison of MAD with standard deviation. The number of points used in relative accuracy analysis were a few thousand, while in BBR accuracy several hundreds, except for bands 3B, 4 and 5 where they were about 50–100. However, the point distribution, excluding cloud and water areas, was not always homogeneous. Areas with poor texture had too little points, while in mountains more matching errors occurred and these points were mostly rejected by the blunder detection methods. For absolute accuracy, in general, 10–20 lakes were used.
Preprocessing adapted to the type of accuracy analysis was always performed. Especially to reduce large radiometric differences and enhance contrast, a set of preprocessing techniques are applied. The following steps are applied to all images and cloud mask data:
- Load 16-bit images
- Linear stretching (the grey values usually cover <2% of the histogram)
- Conversion to 8-bit (due to software implementation issues, easier visualization and less data, without significant loss of information)
- Extract cloud mask data (different masks for KLT tracking method and lake matching)
- Crop both the image and the cloud mask data
These processing steps have also been implemented in Python with the help of OpenCV and Python imaging libraries [13,14]. Additional techniques have been applied to all images for the BBR accuracy investigations as explained in Section 2.3.
2.1. Relative Accuracy
The KLT tracking method has been used for matching of large numbers of image points between two images. A patch size of 15 × 15 pixels is used for matching. The relative accuracy has been investigated using the AVHRR Band-2 images of the same day. Figure 1 shows one full scene, its part used for the relative accuracy evaluation, and the corresponding cloud mask. The cloud mask used in the matching process with KLT tracking method includes: (a) clouds; (b) water; (c) cloud shadows over water, land, and snow; (d) pixels adjacent to cloud over water, land, and snow; and (e) no data (e.g., image borders with no information). The cloud masks of both images have also been combined for the matching of the pairs. Figure 2 shows parts from Band-2 images from two consecutive acquisitions from Metop-A and NOAA-17 on the same day. The points extracted from Metop-A images and successfully matched in NOAA-17 images are shown in both images. The average shift between the two images are 0.2 and 0.1 pixels in x and y, respectively, with standard deviations of 0.1 pixels in both directions.
Figure 1.
AVHRR Band-2 orthoimage acquired from NOAA-18 on 7 May 2008 at 12:08 (left). The subarea used for relative accuracy tests is marked with red rectangle on the full image and shown in upper right image. The cloud mask used for relative accuracy tests is given in the lower right.
Figure 1.
AVHRR Band-2 orthoimage acquired from NOAA-18 on 7 May 2008 at 12:08 (left). The subarea used for relative accuracy tests is marked with red rectangle on the full image and shown in upper right image. The cloud mask used for relative accuracy tests is given in the lower right.
Figure 2.
AVHRR Band-2 orthoimages acquired on 10 July 2008 from Metop-A (left) at 09:36 and NOAA-17 (right) at 09:45. The matched points between the orthoimages are shown in both.
Figure 2.
AVHRR Band-2 orthoimages acquired on 10 July 2008 from Metop-A (left) at 09:36 and NOAA-17 (right) at 09:45. The matched points between the orthoimages are shown in both.
2.2. Absolute Accuracy
The absolute accuracy of the AVHRR imagery has been evaluated based on lakes by estimating the 2D shifts of the lake polygons in image space (i.e., orthoimage coordinate system). Band-2 was used, since lakes are dark in this band and demonstrate a good contrast compared to their surroundings. The evaluation is done based on lakes, since the detection of shorelines in the images is relatively easy and reliable [7]. In addition, with a resolution of 1 km, there are hardly any other stable image features in our test region that are well defined and measurable apart from lake and sea shorelines (and to a much lesser extent forest borders). An alternative would be to use orthorectified images of higher accuracy and resolution and match them with AVHRR orthoimages; however, this matching would be difficult due to large differences (spectral, multitemporal etc.) and was not feasible in this case. Each lake is matched in the images based on the minimization of normalized sum of pixel intensities inside the lake polygons [7]. A total of 20 lakes is used for the evaluations (Figure 3). The lakes were digitized from Landsat-5 orthoimages with 25 m resolution and their locations in the AVHRR images have been computed via coordinate transformations.
Regarding the cloud mask data, only no data and cloud pixels are used as mask in lake matching. Lakes having a cloud cover percentage over 60% are excluded from the process since the matching is usually not successful. The statistical parameters of the shifts are derived from all successfully matched lakes (equally weighted) within an image.
2.3. Band-to-Band Registration Accuracy
Band-to-band registration (BBR) accuracies of the six AVHRR bands are also evaluated using the KLT tracking method [12]. Band-2 images are used as reference for the other five bands (four bands per acquisition, since only 3A or 3B are transmitted to the ground). Band-2 was selected since it was used in the absolute accuracy and thus, through the BBR results, the absolute accuracy of all remaining bands could be determined. Due to large radiometric differences, the following preprocessing techniques are applied to the images prior to the matching:
- Image inverse: Radiometric inverse has been applied to the Bands 3A, 3B, 4 and 5 since the appearance of the mountain areas are similar to those in Band-2 only after the inverse is taken. The functions in Python imaging library is used for this purpose [14].
- Wallis filtering: applied to the image inverses to increase the local contrast and the radiometric similarity of the images. The method has been described in [15] and the software implementation has been done by the Photogrammetry and Remote Sensing Group at ETH Zurich. The same filtering parameters are used for all images (i.e., block size (20 × 20), target mean (127), target standard deviation (50), brightness enforcing constant (1), contrast enforcing constant (0.95)).
- Sobel filtering: applied to the Wallis-filtered images to detect edges in the images. The filter has been implemented here using Python programming language.
- Thresholding: reduces weak edges, often due to noise, from Sobel-filtered images. The threshold is selected as 50% of the cumulative histogram of edge gradients for all images.
Figure 3.
AVHRR Band-2 image part acquired from Metop-A on 7 May 2008 at 10:00 and the lakes used for absolute accuracy investigations.
Figure 3.
AVHRR Band-2 image part acquired from Metop-A on 7 May 2008 at 10:00 and the lakes used for absolute accuracy investigations.
After preliminary tests, a patch size of 9 × 9 pixels is selected as it led to the best matching results. The cloud masks, which have been used in relative accuracy evaluations and described above, have been expanded by 4 pixels (half the patch size) in order to prevent selection of candidate points close to cloud areas and thus reduce the matching errors caused by overlap of the matching patch with clouds. This problem was realized after the relative accuracy investigations. A stricter blunder detection, compared to relative accuracy, was also applied since correct matching of a smaller number of points (due to quasi-identical acquisition times and same orientation) would be sufficient to estimate band-misregistration errors, while the large spectral differences of the matched images caused more matching blunders, needing stricter detection. Figure 4 shows Band-2 and Band-4 images of one acquisition (from Metop-A on 23 January 2008) with all intermediate preprocessing steps applied prior to the matching (i.e., a: 8-bit image; b: inverse; c: Wallis filtered; d: Sobel filtered with thresholding). In the last one (Figure 4d), a smaller area is zoomed in, in order to provide better insight into the process. Figure 5 shows the remaining band images of the same acquisition together with the cloud mask.
Figure 6 shows parts from Band-2 and Band-5 images from Meteop-2 on 10 July 2008 at 09:36. The matching is performed on preprocessed (i.e., Wallis and Sobel filtered with thresholding) images. The candidate points were extracted from Band-2 images and matched in Band-5 images. The plots given in the Figure show the successfully matched points both on preprocessed images used for matching and also on original (linear stretched) images. The average shifts between the two images are 0.3 in both directions, with standard deviations of 0.3 and 0.1 pixels in x and y, respectively.
Figure 4.
AVHRR Band-4 and Band-2 images acquired from Metop-A on 23 January 2008. (a) 8-bit image; (b) inverse ; (c) Wallis filtered; (d) Sobel filtered with thresholding.
Figure 4.
AVHRR Band-4 and Band-2 images acquired from Metop-A on 23 January 2008. (a) 8-bit image; (b) inverse ; (c) Wallis filtered; (d) Sobel filtered with thresholding.
Figure 5.
AVHRR Band-1, 5, and 3A images acquired from Metop-A on 23 January 2008 and the cloud mask used in the matching.
Figure 5.
AVHRR Band-1, 5, and 3A images acquired from Metop-A on 23 January 2008 and the cloud mask used in the matching.
Figure 6.
Parts of AVHRR Band-2 (top) and Band-5 (bottom) images acquired on 10 July 2008 at 09:36 from Metop-A. The matched points between the preprocessed images (right) are shown in both and presented on linear stretched images (left) as well.
Figure 6.
Parts of AVHRR Band-2 (top) and Band-5 (bottom) images acquired on 10 July 2008 at 09:36 from Metop-A. The matched points between the preprocessed images (right) are shown in both and presented on linear stretched images (left) as well.
3. Results and Discussion
The relative and absolute accuracy have been investigated using the AVHRR Band-2 images from 24 different days in 2008. For the relative accuracy evaluations, only the images acquired on the same day have been matched. Up to five acquisitions from all three sensors per day are available. However, due to the inferior GSD of the large off-nadir images, only three images, which cover Switzerland at nadir or near-nadir position, have been used. On the other hand, a total of four images for some days has been evaluated when two images taken from the same sensor have similar off-nadir angles.
Since systematic errors have been observed during the relative accuracy tests (see Section 3.1.2), further analysis has been performed using additional data: i.e., all spectral bands for some days, images processed with two orthoimage generation software versions (SAPS v. 2.0 and 3.0), data of one day in 2009 including images from NOAA-19, and matching original Level 1B reflective images. The data used for these tests are described in the corresponding sections below.
3.1. Relative Accuracy Results
3.1.1. Results (Overall)
A total of 75 images from 24 days was used for matching. The images acquired on the same day from different sensors are matched by pairing the images with shorter time interval in order to reduce matching errors caused by illumination and cloud differences. For example, on the 23 January, three images acquired from Metop-A at 09:33, from NOAA-17 at 09:59, and from NOAA-18 at 11:56 are matched by pairing Metop-A-NOAA-17 and NOAA-17-NOAA-18 images. The matching results of a total of 51 image pairs are given in Table 2 separately for different satellite combinations. P0 and P1 denote the number of candidate and matched points. The Min, Max, Mean, Median (MED), and the standard deviation (σ) values in the second column (Parameter) are derived from all image pairs obtained from the respective satellites given in the first column. The Shift (x, y), MED (x, y), σ(x, y) and MAD(x, y) values are computed from all successfully matched points in an image pair. The Shift, MED, σ and MAD values are presented in the table as absolute shifts (not signed). In this table, only absolute values have been presented to highlight the magnitude of the errors, while with signed shifts the signed mean could be small although the individual shifts large. The Max shifts obtained from the matching of Metop-A and NOAA-17 images are 1.8 (x) and 2.9 (y) pixels. Between the Metop-A-NOAA-18 images, the Max shifts are 1.1 (x) and 0.9 (y) pixels. Between NOAA-17-NOAA-18, the Max shifts are 1.6 (x) and 0.6 (y) pixels. Two NOAA-18 pairs are also evaluated and the Max shifts are 1.1 (x) and 0.1 (y) pixels.
The differences between P0 and P1 are caused by the matching error elimination strategy based on the statistical evaluation (see Section 2). The σ values obtained in the matching are higher than those previously obtained from MSG-SEVIRI images. This can be explained by the larger illumination and atmospheric differences between the images (longer image acquisition intervals), matching images of different satellite sensors and different viewing angles, and existence of local systematic errors in the images. As explained in the next section, systematic errors in stripe layout are observed in all NOAA-18 images.
Figure 7 and Figure 8 (split in two due to space constraints) show the mean x-, y-shifts for all pairs. It can be seen in Figure 8 that the larger shifts occurred after August 2008 on some days only, involving pairs of all three satellites. The green dotted lines in Figure 8 denote the range of vertical axis of Figure 7 for comparison.
Table 2.
Relative accuracy evaluation; statistical values of shifts (x) and (y) from 51 image pairs acquired on 24 different days in 2008. All results except the number of points P0 and P1 are given in pixels and are derived from absolute values.
| Satellites | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Metop-A & NOAA-17 (25 pairs) | Min | 1490 | 1173 | 0.0 | 0.0 | 0.1 | 0.1 | 0.0 | 0.0 | 0.1 | 0.1 |
| Max | 5000 | 4455 | 1.8 | 2.9 | 1.0 | 1.2 | 1.7 | 2.9 | 0.6 | 0.9 | |
| Mean | 4623 | 3701 | 0.4 | 0.4 | 0.4 | 0.4 | 0.3 | 0.4 | 0.3 | 0.3 | |
| MED | 5000 | 3960 | 0.2 | 0.2 | 0.4 | 0.3 | 0.2 | 0.2 | 0.2 | 0.2 | |
| σ | 844 | 733 | 0.5 | 0.7 | 0.2 | 0.3 | 0.4 | 0.7 | 0.1 | 0.2 | |
| Metop-A & NOAA-18 (9 pairs) | Min | 2848 | 1979 | 0.0 | 0.0 | 0.3 | 0.5 | 0.1 | 0.1 | 0.2 | 0.4 |
| Max | 5000 | 4427 | 1.1 | 0.9 | 0.8 | 1.0 | 1.1 | 0.7 | 0.5 | 0.7 | |
| Mean | 4698 | 3632 | 0.2 | 0.3 | 0.5 | 0.6 | 0.2 | 0.2 | 0.3 | 0.5 | |
| MED | 5000 | 3814 | 0.1 | 0.2 | 0.4 | 0.5 | 0.1 | 0.2 | 0.3 | 0.4 | |
| σ | 718 | 834 | 0.3 | 0.2 | 0.2 | 0.2 | 0.3 | 0.2 | 0.1 | 0.1 | |
| NOAA-17 & NOAA-18 (15 pairs) | Min | 1535 | 372 | 0.0 | 0.0 | 0.3 | 0.5 | 0.0 | 0.0 | 0.2 | 0.3 |
| Max | 5000 | 4651 | 1.6 | 0.6 | 1.0 | 0.7 | 1.6 | 0.5 | 0.6 | 0.5 | |
| Mean | 4437 | 2938 | 0.3 | 0.2 | 0.5 | 0.5 | 0.3 | 0.2 | 0.3 | 0.4 | |
| MED | 5000 | 3093 | 0.1 | 0.2 | 0.5 | 0.5 | 0.2 | 0.2 | 0.3 | 0.4 | |
| σ | 1005 | 1097 | 0.4 | 0.1 | 0.2 | 0.1 | 0.4 | 0.1 | 0.1 | 0.0 | |
| NOAA-18 & NOAA-18 (2 pairs) | Min | 3989 | 2558 | 0.7 | 0.0 | 0.8 | 0.7 | 0.6 | 0.0 | 0.4 | 0.5 |
| Max | 5000 | 3179 | 1.1 | 0.1 | 1.4 | 0.9 | 1.1 | 0.1 | 0.7 | 0.6 | |
| Mean | 4495 | 2869 | 0.9 | 0.1 | 1.1 | 0.8 | 0.8 | 0.1 | 0.6 | 0.5 | |
| MED | 4495 | 2869 | 0.9 | 0.1 | 1.1 | 0.8 | 0.8 | 0.1 | 0.6 | 0.5 | |
| σ | 715 | 439 | 0.3 | 0.1 | 0.4 | 0.1 | 0.3 | 0.1 | 0.2 | 0.0 |
Figure 7.
The mean x-, y-shifts obtained from 32 image pairs obtained from Metop-A, NOAA-17 and NOAA-18 between January and August 2008. Horizontal axis: Time, Vertical axis: pixels.
Figure 7.
The mean x-, y-shifts obtained from 32 image pairs obtained from Metop-A, NOAA-17 and NOAA-18 between January and August 2008. Horizontal axis: Time, Vertical axis: pixels.
Figure 8.
The mean x-, y-shifts obtained from 19 image pairs obtained from Metop-A, NOAA-17 and NOAA-18 between September and December 2008. Horizontal axis: Time, Vertical axis: pixels.
Figure 8.
The mean x-, y-shifts obtained from 19 image pairs obtained from Metop-A, NOAA-17 and NOAA-18 between September and December 2008. Horizontal axis: Time, Vertical axis: pixels.
3.1.2. Systematic Errors
During the relative accuracy tests, systematic stripes have been observed in some of the image pairs. After visual inspections on the shift plots of all 51 pairs, it has been seen that this effect exists in all pairs where NOAA-18 images are involved (i.e., Metop-A/NOAA-18, NOAA-17/NOAA-18, and NOAA-18/NOAA-18 pairs). These systematic errors have been investigated further using additional data as mentioned before Section 3.1.1.
Figure 9, Figure 10 and Figure 11 show the shift plots of three image pairs acquired on 26 December 2008. In the figures of this Section, the vectors are not shown with background of the original images, since the vectors would be hardly visible, but on images that were especially preprocessed to suppress the background, keeping water surfaces dark for better visualization (thus, what appears in the background as being clouds, is not really so). In Figure 10, the image pair, which contains NOAA-18 image, shows systematic errors as stripes across the flight direction of NOAA-18. This sort of systematic errors can be seen also in the matching of two NOAA-18 images from the same day (Figure 11), but due to phase differences the stripes may appear differently. Stripes have also been observed along the whole length of an AVHRR image (in N. Europe and N. Africa).
The striping errors occur for NOAA-18 and for NOAA-19 (for the latter only one day in 2009 has been checked). Both satellites take images while ascending, however no reason could been seen why ascending image acquisition should have such systematic errors. The stripes are perpendicular to the satellite track and alternating with a width of about 30 pixels (corresponding to 5 s of data acquisition), i.e., one stripe has shifts of about 1 pixel in flight direction, the next stripe has no shifts, etc. The stripes appear in all spectral bands. The matched images were from the same day and using the same DSM as input to the orthorectification. Although the GCPs were different for the various AVHRR images, their number and distribution were quite similar for the same day, with differences mainly due to varying clouds and less illumination. The stripes appear over the whole length of the original AVHRR images (from Scandinavia to N. Africa). This investigation has been done since in the AVHRR orientation with SAPS, the image is divided into three almost equal segments, and orientation errors may differ for each segment. They also have a very high frequency, which cannot be explained by errors in sensor orientation or DSM. The error is assumed to be existing in the input data for orthorectification (Level 1B), due perhaps to errors of geometric calibration, scanning or time synchronization between the system’s position and attitude determination units. An additional indication that the systematic errors exist in the Level 1B data is the inclination of the stripes (see Figure 10, Figure 11, Figure 12 and Figure 13). The stripe direction (long dimension) is not aligned with North/South used in the orthoimage generation but is perpendicular to the flight direction of ascending NOAA-18/19 satellites. Thus, the Level 1B images have also been checked.
Figure 9.
Image space shifts obtained from matching of NOAA-17 and Metop-A images acquired on 26 December 2008. The axes denote the image coordinates.
Figure 9.
Image space shifts obtained from matching of NOAA-17 and Metop-A images acquired on 26 December 2008. The axes denote the image coordinates.
Figure 10.
Image space shifts obtained from matching of Metop-A and NOAA-18 images acquired on 26 December 2008. The axes denote the image coordinates.
Figure 10.
Image space shifts obtained from matching of Metop-A and NOAA-18 images acquired on 26 December 2008. The axes denote the image coordinates.
Figure 11.
Image space shifts obtained from matching of two NOAA-18 images acquired on 26 December 2008. The axes denote the image coordinates.
Figure 11.
Image space shifts obtained from matching of two NOAA-18 images acquired on 26 December 2008. The axes denote the image coordinates.
Additional data acquired on 6 August 2009 has been provided by the University of Bern, where one acquisition from NOAA-19 satellite is also included. Similar striping errors in across-track direction are also observed in the NOAA-19 image. The shift plots from the matching of NOAA-17 and NOAA-18 images with the NOAA-19 image are provided in Figure 12a,b.
Figure 12.
Image space shifts obtained from matching of NOAA-18-NOAA-19 (a) and NOAA-17-NOAA-19 (b) images acquired on 6 August 2009. The axes denote the image coordinates.
Figure 12.
Image space shifts obtained from matching of NOAA-18-NOAA-19 (a) and NOAA-17-NOAA-19 (b) images acquired on 6 August 2009. The axes denote the image coordinates.
The analysis of NOAA-18 images has been extended by using the bands for the some days. As can be seen in the Figure 13a,b, the stripes are clearly visible in the matching of two different bands (Band-1 and Band-5) acquired on 8 February 2008.
Figure 13.
Image space shifts obtained from the matching of the Band-1 (a) and Band-5 (b) images acquired from NOAA-17 and NOAA-18 on February 8, 2008. The axes denote the image coordinates.
Figure 13.
Image space shifts obtained from the matching of the Band-1 (a) and Band-5 (b) images acquired from NOAA-17 and NOAA-18 on February 8, 2008. The axes denote the image coordinates.
All data were generated by SAPS v. 2.0. We also processed some data generated by SAPS v3.0 and the results were practically identical. To reduce the risk that the systematic stripe problem was due to a software problem, the University of Bern analysed all outputs of the SAPS processing software and checked the quality of the georeferencing and matching between MODIS and AVHRR. No problem could be detected, even using different versions of the SAPS software. The author of SAPS was also informed about this systematic error but could not provide any explanation (Khlopenkov, personal communications). Furthermore, A.P. Cracknell (personal communications) who has worked extensively with AVHRR was contacted for a possible explanation of the systematic stripes, but he could not provide one.
Figure 14.
Matching of Band-2 L1B images (input to orthorectification). Top: left NOAA-17 after being co-registered with NOAA-18, right NOAA-18. Bottom: shifts after matching NOAA-17 with NOAA-18. The large vectors correspond to matching errors due to multitemporal differences (e.g., clouds) and significant other differences in grey values due to different illumination conditions, viewing angle, etc. The systematic stripes are clearly visible. They have a width of about 30 pixels, and alternating shifts of about 1 and 0 pixels. The axes in the plot (bottom) denote the image coordinates.
Figure 14.
Matching of Band-2 L1B images (input to orthorectification). Top: left NOAA-17 after being co-registered with NOAA-18, right NOAA-18. Bottom: shifts after matching NOAA-17 with NOAA-18. The large vectors correspond to matching errors due to multitemporal differences (e.g., clouds) and significant other differences in grey values due to different illumination conditions, viewing angle, etc. The systematic stripes are clearly visible. They have a width of about 30 pixels, and alternating shifts of about 1 and 0 pixels. The axes in the plot (bottom) denote the image coordinates.
Original Level-1B AVHRR reflective images (radiometrically calibrated) have also been used to find out if the systematic errors exist also in the original images. After the finalization of this project, we used accurate coregistration methods to perform the matching between L1B images within another project [16]. These methods computed a co-registration with 3rd degree polynomial (20 parameters) between NOAA-17 and NOAA-18 images by using about 350 hundred points found and matched with SIFT, and manually edited to exclude blunders. As can be seen in Figure 14, the systematic stripes are visible in the results. The matching between the transformed NOAA-17 and the original NOAA-18 image are dense on a regular grid, instead of using good texture points as with KLT; thus, the systematic pattern is more visible. The matching algorithm used between the co-registered images was a Fourier-Transform based Cross-Correlation (FTCC) as explained in [16]. Thus, we can show that this error exists in the L1B data and is not due to the orthorectification software. This error maybe due to errors of calibration, scanning or time synchronization between the system’s main three units. We will do further future investigations and contact experts of the NOAA-18 sensors to try to find the cause.
3.2. Absolute Accuracy Results
A total of 72 images acquired on 23 different days (one day excluded due to total cloud coverage) has been processed using lake matching. Figure 15 shows the plots of two lake polygons at their initial and the matched positions in the images. The results are summarized in Table 3, Table 4 and Table 5 for Metop-A, NOAA-17 and NOAA-18. The shifts in Table 3, Table 4 and Table 5 show that up to 4 pixels, 1.6 pixels and 1 pixel shifts are present in Metop-A, NOAA-17 and NOAA-18 images, respectively. The shifts are shown with signs, where +Δy indicates a shift to the North and +Δx indicates a shift to the East.
Figure 15.
Lake matching results from Lake Leman in a NOAA-17 image (left) and Lake Garda in a Metop-A image (right). The green polygons represent the initial positions of the lake points and the red polygons represent the matched positions.
Figure 15.
Lake matching results from Lake Leman in a NOAA-17 image (left) and Lake Garda in a Metop-A image (right). The green polygons represent the initial positions of the lake points and the red polygons represent the matched positions.
The standard deviations (σx, σy) given for each date show the variability in the shifts of all lakes in an image. The larger standard deviations, which occurred in the shifts of all lakes within an image, have been observed in the more off-nadir images. The standard deviations increase even more (up to 0.7 pixels) at the off-nadir images, which are excluded from the analysis in the following tables (images closer to the nadir are preferred for each day). However, in each table, there are days with two acquisitions, where both images were acquired from a relatively off-nadir position.
When the mean and σ values (given in the lower part of the Tables) of Shift (x) and Shift (y) parameters are compared for all three satellites, it can be said that NOAA-17 and NOAA-18 images have better absolute accuracy than Metop-A. The mean (x,y) shifts obtained from all 24 images of NOAA-17 are (0.2, −0.2) pixels with shift standard deviations of (0.3, 0.2) pixels. For NOAA-18, the mean (x, y) shifts are (0.1, −0.3) pixels with standard deviations of (0.2, 0.2) pixels. For Metop-A images, although the mean (x, y) shifts are small (−0.1, −0.4 pixels), the standard deviations of shifts are high (0.4, 1.1 pixels), which reflects existence of large shifts, especially in y, on three days. The results of these days were visually checked, but no matching errors could be found.
Table 3.
Absolute accuracy evaluation; statistical summary from the lake shifts (equally weighted) in AVHRR images acquired from Metop-A. All values are in pixels except the number of lakes and are derived from signed values.
| Date | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) | No. of Lakes |
|---|---|---|---|---|---|---|---|---|---|
| 23.01.2008 | −0.4 | −0.3 | 0.2 | 0.3 | −0.5 | −0.3 | 0.2 | 0.2 | 20 |
| 24.01.2008 | −0.1 | −0.8 | 0.3 | 0.3 | −0.2 | −0.8 | 0.2 | 0.2 | 20 |
| 08.02.2008 | −0.2 | −0.6 | 0.2 | 0.2 | −0.2 | −0.6 | 0.2 | 0.2 | 20 |
| 09.02.2008 | −0.4 | −1.5 | 0.3 | 0.3 | −0.4 | −1.6 | 0.3 | 0.2 | 20 |
| 09.02.2008 | 0.3 | −0.1 | 0.3 | 0.4 | 0.4 | −0.1 | 0.3 | 0.3 | 20 |
| 06.03.2008 | −0.1 | −0.2 | 0.2 | 0.3 | −0.1 | −0.2 | 0.2 | 0.2 | 20 |
| 29.03.2008 | 0.2 | −0.4 | 0.2 | 0.3 | 0.2 | −0.5 | 0.1 | 0.3 | 18 |
| 01.04.2008 | 0.0 | −0.4 | 0.1 | 0.1 | 0.0 | −0.4 | 0.1 | 0.1 | 7 |
| 07.05.2008 | 0.1 | 0.0 | 0.2 | 0.2 | 0.0 | 0.0 | 0.1 | 0.1 | 20 |
| 08.05.2008 | 0.0 | −0.3 | 0.1 | 0.1 | 0.0 | −0.3 | 0.1 | 0.1 | 19 |
| 21.06.2008 | 0.0 | −0.2 | 0.2 | 0.2 | 0.0 | −0.2 | 0.1 | 0.0 | 20 |
| 25.06.2008 | −0.1 | −0.1 | 0.1 | 0.1 | −0.1 | −0.1 | 0.1 | 0.1 | 19 |
| 10.07.2008 | 0.0 | −0.1 | 0.2 | 0.2 | 0.0 | −0.1 | 0.1 | 0.1 | 19 |
| 15.07.2008 | −0.2 | 0.4 | 0.1 | 0.3 | −0.2 | 0.4 | 0.1 | 0.2 | 20 |
| 18.08.2008 | 0.0 | −0.2 | 0.2 | 0.2 | 0.0 | −0.3 | 0.1 | 0.1 | 19 |
| 30.08.2008 | −0.1 | −0.1 | 0.4 | 0.3 | −0.1 | −0.1 | 0.2 | 0.2 | 16 |
| 09.09.2008 | −1.7 | 2.0 | 0.4 | 0.5 | −1.6 | 1.8 | 0.3 | 0.3 | 20 |
| 28.09.2008 | 0.1 | −0.1 | 0.4 | 0.4 | 0.2 | −0.1 | 0.2 | 0.3 | 14 |
| 05.10.2008 | 0.0 | −0.4 | 0.2 | 0.2 | 0.0 | −0.4 | 0.1 | 0.2 | 20 |
| 11.10.2008 | −0.3 | 0.0 | 0.2 | 0.3 | −0.3 | 0.0 | 0.1 | 0.2 | 11 |
| 15.11.2008 | 0.4 | −4.0 | 0.4 | 0.4 | 0.5 | −4.2 | 0.3 | 0.2 | 11 |
| 27.11.2008 | −0.4 | 0.1 | 0.2 | 0.4 | −0.4 | 0.1 | 0.2 | 0.3 | 14 |
| 08.12.2008 | 0.3 | −2.3 | 0.3 | 0.4 | 0.3 | −2.4 | 0.2 | 0.4 | 13 |
| 26.12.2008 | 0.0 | −0.4 | 0.3 | 0.3 | 0.0 | −0.5 | 0.2 | 0.3 | 18 |
| Min | −1.7 | −4.0 | 0.1 | 0.1 | −1.6 | −4.2 | 0.1 | 0.0 | 7 |
| Max | 0.4 | 2.0 | 0.4 | 0.5 | 0.5 | 1.8 | 0.3 | 0.4 | 20 |
| Mean | −0.1 | −0.4 | 0.2 | 0.3 | −0.1 | −0.5 | 0.2 | 0.2 | 17.4 |
| MED | 0.0 | −0.2 | 0.2 | 0.3 | 0.0 | −0.3 | 0.2 | 0.2 | 19.0 |
| σ | 0.4 | 1.1 | 0.1 | 0.1 | 0.4 | 1.1 | 0.1 | 0.1 | 3.7 |
Table 4.
Absolute accuracy evaluation; statistical summary from the lake shifts in AVHRR images acquired from NOAA-17. All values are in pixels except the number of lakes and are derived from signed values.
| Date | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) | No. of Lakes |
|---|---|---|---|---|---|---|---|---|---|
| 23.01.2008 | 0.0 | −0.2 | 0.2 | 0.3 | 0.0 | −0.2 | 0.2 | 0.3 | 20 |
| 24.01.2008 | 0.2 | −0.4 | 0.2 | 0.2 | 0.2 | −0.5 | 0.1 | 0.3 | 20 |
| 08.02.2008 | −0.3 | 0.0 | 0.2 | 0.3 | −0.3 | −0.1 | 0.1 | 0.3 | 20 |
| 09.02.2008 | 0.0 | −0.3 | 0.2 | 0.2 | 0.0 | −0.3 | 0.1 | 0.2 | 20 |
| 06.03.2008 | 0.1 | −0.2 | 0.2 | 0.2 | 0.1 | −0.2 | 0.1 | 0.2 | 19 |
| 29.03.2008 | 0.4 | −0.1 | 0.2 | 0.2 | 0.3 | −0.2 | 0.1 | 0.1 | 18 |
| 01.04.2008 | 0.0 | −0.2 | 0.1 | 0.2 | 0.0 | −0.3 | 0.1 | 0.1 | 9 |
| 07.05.2008 | 0.2 | −0.1 | 0.2 | 0.3 | 0.2 | −0.2 | 0.1 | 0.1 | 20 |
| 08.05.2008 | −0.1 | −0.2 | 0.4 | 0.4 | −0.3 | −0.3 | 0.3 | 0.2 | 19 |
| 08.05.2008 | 0.1 | −0.2 | 0.4 | 0.4 | 0.1 | −0.1 | 0.2 | 0.3 | 19 |
| 21.06.2008 | 0.2 | 0.1 | 0.3 | 0.2 | 0.2 | 0.0 | 0.2 | 0.1 | 17 |
| 25.06.2008 | 0.2 | 0.0 | 0.2 | 0.1 | 0.3 | 0.0 | 0.1 | 0.0 | 19 |
| 10.07.2008 | 0.1 | 0.0 | 0.2 | 0.2 | 0.1 | 0.0 | 0.1 | 0.2 | 20 |
| 15.07.2008 | 0.2 | 0.0 | 0.2 | 0.2 | 0.3 | 0.1 | 0.1 | 0.1 | 20 |
| 18.08.2008 | 0.0 | −0.1 | 0.2 | 0.1 | 0.1 | 0.0 | 0.1 | 0.1 | 20 |
| 30.08.2008 | 0.1 | −0.1 | 0.2 | 0.1 | 0.0 | 0.0 | 0.1 | 0.1 | 17 |
| 09.09.2008 | 0.0 | −0.2 | 0.1 | 0.2 | 0.0 | −0.2 | 0.1 | 0.1 | 20 |
| 28.09.2008 | 0.3 | −0.2 | 0.5 | 0.3 | 0.4 | −0.1 | 0.2 | 0.2 | 15 |
| 05.10.2008 | 0.3 | −0.2 | 0.2 | 0.2 | 0.2 | −0.2 | 0.1 | 0.2 | 19 |
| 11.10.2008 | 1.6 | 0.3 | 0.4 | 0.3 | 1.6 | 0.3 | 0.3 | 0.2 | 12 |
| 15.11.2008 | 0.2 | −0.6 | 0.3 | 0.4 | 0.2 | −0.6 | 0.2 | 0.3 | 13 |
| 27.11.2008 | −0.2 | −0.6 | 0.3 | 0.3 | −0.2 | −0.6 | 0.2 | 0.4 | 12 |
| 08.12.2008 | 0.0 | −0.3 | 0.3 | 0.3 | −0.1 | −0.3 | 0.1 | 0.3 | 13 |
| 26.12.2008 | 0.1 | −0.4 | 0.3 | 0.3 | 0.1 | −0.4 | 0.1 | 0.3 | 18 |
| Min | −0.3 | −0.6 | 0.1 | 0.1 | −0.3 | −0.6 | 0.1 | 0.0 | 9 |
| Max | 1.6 | 0.3 | 0.5 | 0.4 | 1.6 | 0.3 | 0.3 | 0.4 | 20 |
| Mean | 0.2 | −0.2 | 0.3 | 0.2 | 0.1 | −0.2 | 0.1 | 0.2 | 17.5 |
| MED | 0.1 | −0.2 | 0.2 | 0.2 | 0.1 | −0.2 | 0.1 | 0.2 | 19.0 |
| σ | 0.3 | 0.2 | 0.1 | 0.1 | 0.4 | 0.2 | 0.1 | 0.1 | 3.3 |
Table 5.
Absolute accuracy evaluation; statistical summary from the lake shifts in AVHRR images acquired from NOAA-18. All values are in pixels except the number of lakes and are derived from signed values.
| Date | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) | No. of Lakes |
|---|---|---|---|---|---|---|---|---|---|
| 23.01.2008 | 0.1 | −0.6 | 0.2 | 0.5 | 0.2 | −0.6 | 0.2 | 0.5 | 17 |
| 24.01.2008 | 0.2 | −0.5 | 0.3 | 0.5 | 0.3 | −0.8 | 0.2 | 0.1 | 16 |
| 08.02.2008 | −0.2 | −0.5 | 0.2 | 0.3 | −0.2 | −0.5 | 0.2 | 0.2 | 20 |
| 09.02.2008 | 0.0 | −0.6 | 0.2 | 0.3 | 0.1 | −0.6 | 0.1 | 0.3 | 20 |
| 06.03.2008 | 0.1 | −0.2 | 0.2 | 0.4 | 0.1 | −0.1 | 0.2 | 0.2 | 18 |
| 29.03.2008 | 0.1 | −0.4 | 0.2 | 0.3 | 0.1 | −0.4 | 0.1 | 0.2 | 20 |
| 01.04.2008 | −0.1 | 0.1 | 0.4 | 0.3 | 0.0 | 0.1 | 0.2 | 0.3 | 9 |
| 07.05.2008 | 0.0 | −0.1 | 0.2 | 0.3 | 0.0 | −0.1 | 0.2 | 0.2 | 17 |
| 08.05.2008 | 0.1 | −0.2 | 0.3 | 0.2 | 0.2 | −0.2 | 0.1 | 0.1 | 19 |
| 21.06.2008 | 0.2 | −0.1 | 0.2 | 0.3 | 0.2 | −0.1 | 0.1 | 0.2 | 18 |
| 25.06.2008 | 0.2 | −0.2 | 0.2 | 0.3 | 0.3 | −0.3 | 0.2 | 0.3 | 12 |
| 10.07.2008 | 0.2 | −0.1 | 0.2 | 0.3 | 0.3 | 0.0 | 0.1 | 0.2 | 20 |
| 15.07.2008 | 0.0 | −0.4 | 0.3 | 0.5 | −0.1 | −0.7 | 0.2 | 0.2 | 20 |
| 18.08.2008 | 0.1 | −0.1 | 0.2 | 0.4 | 0.1 | −0.1 | 0.1 | 0.4 | 19 |
| 30.08.2008 | 0.2 | −0.2 | 0.2 | 0.4 | 0.2 | −0.2 | 0.1 | 0.3 | 19 |
| 09.09.2008 | 0.2 | −0.2 | 0.2 | 0.3 | 0.2 | −0.1 | 0.1 | 0.2 | 17 |
| 28.09.2008 | 0.1 | −0.2 | 0.1 | 0.3 | 0.1 | −0.2 | 0.1 | 0.2 | 17 |
| 05.10.2008 | 0.3 | −0.2 | 0.2 | 0.3 | 0.3 | −0.3 | 0.1 | 0.2 | 8 |
| 11.10.2008 | 0.3 | −0.6 | 0.3 | 0.4 | 0.4 | −0.6 | 0.2 | 0.1 | 12 |
| 15.11.2008 | 0.0 | −1.0 | 0.3 | 0.4 | 0.0 | −1.1 | 0.2 | 0.3 | 14 |
| 27.11.2008 | 0.0 | −0.3 | 0.4 | 0.3 | 0.0 | −0.3 | 0.2 | 0.3 | 17 |
| 08.12.2008 | 0.1 | −0.4 | 0.2 | 0.5 | 0.0 | −0.5 | 0.1 | 0.3 | 13 |
| 08.12.2008 | 0.7 | −0.7 | 0.3 | 0.3 | 0.6 | −0.7 | 0.2 | 0.3 | 16 |
| 26.12.2008 | 0.0 | −0.5 | 0.4 | 0.5 | −0.1 | −0.7 | 0.3 | 0.1 | 20 |
| Min | −0.2 | −1.0 | 0.1 | 0.2 | −0.2 | −1.1 | 0.1 | 0.1 | 8 |
| Max | 0.7 | 0.1 | 0.4 | 0.5 | 0.6 | 0.1 | 0.3 | 0.5 | 20 |
| Mean | 0.1 | −0.3 | 0.2 | 0.4 | 0.1 | −0.4 | 0.1 | 0.2 | 16.6 |
| MED | 0.1 | −0.3 | 0.2 | 0.3 | 0.1 | −0.3 | 0.2 | 0.2 | 17.0 |
| σ | 0.2 | 0.2 | 0.1 | 0.1 | 0.2 | 0.3 | 0.1 | 0.1 | 3.5 |
3.3. BBR Accuracy Results
The images of 63 acquisitions from 24 different days have been processed for the BBR accuracy analysis. The plots of the shift vectors for two different bands of the NOAA-18 acquisition on 27 November 2008 are given in Figure 16. The matching results are provided in Table 6, Table 7, Table 8, Table 9 and Table 10, one per band, and given separately per satellite. P0 and P1 indicate the numbers of the candidate and the matched points. These values are in general much lower than for the relative accuracy investigations, due to matching difficulties and mainly stricter blunder detection criteria (thresholds). Given that the same georeferencing parameters and DSM are used for all bands, the test results provided here should represent the BBR accuracy of the original bands plus the accuracy of matching, the latter having its influence even if matching the original images. The Shift, MED, σ and MAD values are computed using the absolute shifts (not signed).
As can be seen in Table 6, the mean (x,y) shifts are (0.07, 0.03) pixels for Metop-A, (0.04, 0.03) pixels for NOAA-17 and (0.08, 0.09) pixels for NOAA-18 between Bands 1 and 2. The maximum (x,y) shifts are (0.12, 0.08) pixels for Metop-A, (0.11, 0.12) pixels for NOAA-17 and (0.15, 0.21) pixels for NOAA-18. For Band-3A (Table 7), the mean (x, y) shifts are (0.17, 0.05) pixels for Metop-A and (0.15, 0.04) pixels for NOAA-17. The maximum (x, y) shifts are (0.28, 0.14) pixels for Metop-A and (0.35, 0.10) pixels for NOAA-17. For Band-3B (Table 8), only the data from NOAA-18 is available in the given dataset. It can be seen that the mean x-shift is comparable to those of Band-3A (0.16 pixels). The mean y-shift is slightly higher (0.13 pixels) than that of Band-3A. The maximum (x, y) shifts are (0.42, 0.24) pixels. The larger shifts are observed in the more off-nadir acquisitions.
Table 9 and Table 10 show the results obtained from the matching of Band-2 images with Bands 4 and 5, respectively. It can be seen in the mean, minimum, and maximum values of each sensor that the results are very similar for these two bands. The images are visually very similar as well. Regarding Band-4, the mean (x, y) shifts are (0.30, 0.07) pixels for Metop-A, (0.26, 0.09) pixels for NOAA-17 and (0.17, 0.11) pixels for NOAA-18. The maximum (x, y) shifts are (0.41, 0.20) pixels for Metop-A, (0.40, 0.21) pixels for NOAA-17 and (0.31, 0.27) pixels for NOAA-18. Regarding Band-5 (Table 10), the mean (x,y) shifts are (0.29, 0.09) pixels for Metop-A, (0.27, 0.09) pixels for NOAA-17 and (0.13, 0.13) pixels for NOAA-18. The maximum (x, y) shifts are (0.46, 0.24) pixels for Metop-A, (0.39, 0.26) pixels for NOAA-17 and (0.27, 0.29) pixels for NOAA-18. Similar to Band-3B, the larger shifts occur in the more off-nadir acquisitions. Results for Band-3A, even more for Band-3B and especially Bands 4 and 5 are worse than for Band-1, due to more spectral difference to the reference Band-2 and thus more difficult matching.
Figure 16.
Image space shifts obtained from the matching of the Band-1 (left) and Band-4 (right) images with Band-2 image acquired from NOAA-18 on 27 November 2008. The axes denote the image coordinates.
Figure 16.
Image space shifts obtained from the matching of the Band-1 (left) and Band-4 (right) images with Band-2 image acquired from NOAA-18 on 27 November 2008. The axes denote the image coordinates.
Table 6.
BBR accuracy evaluation results between AVHRR Bands 1 and 2; statistical summary from the point matching results of all images for the three satellites. The given statistics are in pixels and are derived from the absolute values.
| Satellite | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Metop-A (21 image pairs) | Min | 5763 | 38 | 0.03 | 0.00 | 0.04 | 0.04 | 0.03 | 0.00 | 0.03 | 0.03 |
| Max | 10,000 | 1419 | 0.12 | 0.08 | 0.22 | 0.23 | 0.12 | 0.10 | 0.17 | 0.18 | |
| Mean | 8503 | 321 | 0.07 | 0.03 | 0.12 | 0.11 | 0.07 | 0.04 | 0.09 | 0.08 | |
| MED | 9223 | 209 | 0.07 | 0.03 | 0.11 | 0.10 | 0.07 | 0.03 | 0.09 | 0.07 | |
| σ | 1650 | 363 | 0.02 | 0.02 | 0.06 | 0.05 | 0.03 | 0.03 | 0.05 | 0.04 | |
| NOAA-17 (22 image pairs) | Min | 4267 | 25 | 0.00 | 0.00 | 0.04 | 0.04 | 0.00 | 0.00 | 0.03 | 0.03 |
| Max | 10,000 | 1323 | 0.11 | 0.12 | 0.26 | 0.24 | 0.12 | 0.11 | 0.20 | 0.19 | |
| Mean | 8080 | 383 | 0.04 | 0.03 | 0.12 | 0.11 | 0.04 | 0.04 | 0.10 | 0.09 | |
| MED | 8881 | 184 | 0.03 | 0.03 | 0.12 | 0.10 | 0.03 | 0.03 | 0.08 | 0.08 | |
| σ | 2076 | 435 | 0.03 | 0.03 | 0.07 | 0.06 | 0.04 | 0.03 | 0.06 | 0.05 | |
| NOAA-18 (20 image pairs) | Min | 3594 | 21 | 0.00 | 0.00 | 0.06 | 0.05 | 0.00 | 0.00 | 0.00 | 0.03 |
| Max | 10,000 | 1436 | 0.15 | 0.21 | 0.27 | 0.22 | 0.15 | 0.21 | 0.18 | 0.21 | |
| Mean | 7421 | 364 | 0.08 | 0.09 | 0.14 | 0.13 | 0.08 | 0.09 | 0.08 | 0.10 | |
| MED | 7416 | 152 | 0.08 | 0.06 | 0.13 | 0.13 | 0.08 | 0.06 | 0.09 | 0.07 | |
| σ | 2326 | 457 | 0.05 | 0.06 | 0.07 | 0.05 | 0.05 | 0.06 | 0.05 | 0.06 |
Table 7.
BBR accuracy evaluation results between AVHRR Bands 2 and 3A; statistical summary from the point matching results of all images for the two satellites. The given statistics are in pixels and are derived from the absolute values.
| Satellite | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Metop-A (21 image pairs) | Min | 5763 | 45 | 0.11 | 0.00 | 0.06 | 0.04 | 0.11 | 0.00 | 0.05 | 0.03 |
| Max | 10,000 | 1692 | 0.28 | 0.14 | 0.17 | 0.21 | 0.28 | 0.14 | 0.13 | 0.14 | |
| Mean | 8503 | 781 | 0.17 | 0.05 | 0.11 | 0.13 | 0.17 | 0.05 | 0.09 | 0.10 | |
| MED | 9223 | 727 | 0.17 | 0.05 | 0.11 | 0.13 | 0.17 | 0.05 | 0.08 | 0.10 | |
| σ | 1650 | 510 | 0.04 | 0.04 | 0.02 | 0.03 | 0.04 | 0.04 | 0.02 | 0.02 | |
| NOAA-17 (22 image pairs) | Min | 4267 | 34 | 0.01 | 0.01 | 0.05 | 0.10 | 0.03 | 0.01 | 0.04 | 0.07 |
| Max | 10,000 | 1854 | 0.35 | 0.10 | 0.20 | 0.17 | 0.35 | 0.10 | 0.17 | 0.14 | |
| Mean | 8080 | 891 | 0.16 | 0.04 | 0.12 | 0.13 | 0.16 | 0.05 | 0.10 | 0.10 | |
| MED | 8881 | 980 | 0.15 | 0.04 | 0.12 | 0.13 | 0.15 | 0.05 | 0.09 | 0.10 | |
| σ | 2076 | 562 | 0.08 | 0.03 | 0.03 | 0.02 | 0.07 | 0.03 | 0.03 | 0.02 |
Table 8.
BBR accuracy evaluation results between AVHRR Bands 2 and 3B; statistical summary from the point matching results of all NOAA-18 satellites. The given statistics are in pixels and are derived from the absolute values.
| Satellite | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| NOAA-18 (15 image pairs) | Min | 4301 | 12 | 0.00 | 0.01 | 0.11 | 0.08 | 0.01 | 0.01 | 0.09 | 0.06 |
| Max | 10000 | 153 | 0.42 | 0.24 | 0.42 | 0.19 | 0.42 | 0.25 | 0.35 | 0.17 | |
| Mean | 7945 | 59 | 0.16 | 0.13 | 0.24 | 0.15 | 0.18 | 0.14 | 0.19 | 0.11 | |
| MED | 8712 | 47 | 0.16 | 0.14 | 0.25 | 0.15 | 0.18 | 0.16 | 0.18 | 0.12 | |
| σ | 2131 | 45 | 0.13 | 0.07 | 0.09 | 0.04 | 0.12 | 0.07 | 0.08 | 0.03 |
Table 9.
BBR accuracy evaluation results between AVHRR Bands 2 and 4; statistical summary from the point matching results of all images for the three satellites. The given statistics are in pixels and are derived from the absolute values.
| Satellite | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Metop-A (20 image pairs) | Min | 5763 | 11 | 0.18 | 0.00 | 0.07 | 0.10 | 0.18 | 0.00 | 0.04 | 0.08 |
| Max | 10,000 | 183 | 0.41 | 0.20 | 0.30 | 0.24 | 0.40 | 0.21 | 0.29 | 0.19 | |
| Mean | 8499 | 60 | 0.30 | 0.07 | 0.15 | 0.15 | 0.30 | 0.08 | 0.12 | 0.12 | |
| MED | 9330 | 50 | 0.31 | 0.05 | 0.14 | 0.16 | 0.31 | 0.08 | 0.11 | 0.12 | |
| σ | 1692 | 42 | 0.05 | 0.06 | 0.06 | 0.04 | 0.05 | 0.06 | 0.06 | 0.03 | |
| NOAA-17 (21 image pairs) | Min | 4267 | 10 | 0.02 | 0.02 | 0.08 | 0.09 | 0.10 | 0.00 | 0.06 | 0.06 |
| Max | 10,000 | 252 | 0.40 | 0.21 | 0.45 | 0.34 | 0.42 | 0.21 | 0.27 | 0.35 | |
| Mean | 8097 | 80 | 0.26 | 0.09 | 0.16 | 0.17 | 0.28 | 0.10 | 0.12 | 0.13 | |
| MED | 8948 | 68 | 0.26 | 0.09 | 0.14 | 0.16 | 0.28 | 0.11 | 0.10 | 0.11 | |
| σ | 2126 | 66 | 0.08 | 0.06 | 0.09 | 0.06 | 0.06 | 0.06 | 0.05 | 0.07 | |
| NOAA-18 (22 image pairs) | Min | 3339 | 10 | 0.02 | 0.00 | 0.12 | 0.06 | 0.01 | 0.00 | 0.05 | 0.05 |
| Max | 10,000 | 235 | 0.31 | 0.27 | 0.41 | 0.33 | 0.34 | 0.26 | 0.29 | 0.27 | |
| Mean | 7052 | 89 | 0.17 | 0.11 | 0.20 | 0.16 | 0.18 | 0.12 | 0.15 | 0.13 | |
| MED | 6532 | 88 | 0.17 | 0.09 | 0.19 | 0.16 | 0.17 | 0.11 | 0.15 | 0.13 | |
| σ | 2515 | 70 | 0.07 | 0.07 | 0.07 | 0.06 | 0.08 | 0.07 | 0.05 | 0.05 |
Table 10.
BBR accuracy evaluation results between AVHRR Bands 2 and 5; statistical summary from the point matching results of all images for the three satellites. The given statistics are in pixels and are derived from the absolute values.
| Satellite | Parameter | P0 | P1 | Shift (x) | Shift (y) | σ (x) | σ (y) | MED (x) | MED (y) | MAD (x) | MAD (y) |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Metop-A (20 image pairs) | Min | 5763 | 13 | 0.13 | 0.01 | 0.04 | 0.12 | 0.10 | 0.01 | 0.04 | 0.08 |
| Max | 10,000 | 144 | 0.46 | 0.24 | 0.32 | 0.22 | 0.44 | 0.25 | 0.30 | 0.16 | |
| Mean | 8499 | 54 | 0.29 | 0.09 | 0.15 | 0.16 | 0.30 | 0.10 | 0.12 | 0.12 | |
| MED | 9330 | 42 | 0.30 | 0.08 | 0.15 | 0.16 | 0.31 | 0.09 | 0.11 | 0.13 | |
| σ | 1692 | 35 | 0.08 | 0.07 | 0.06 | 0.03 | 0.08 | 0.07 | 0.05 | 0.02 | |
| NOAA-17 (19 image pairs) | Min | 4267 | 13 | 0.12 | 0.00 | 0.07 | 0.06 | 0.12 | 0.01 | 0.05 | 0.05 |
| Max | 10,000 | 217 | 0.39 | 0.26 | 0.25 | 0.26 | 0.40 | 0.25 | 0.16 | 0.20 | |
| Mean | 8400 | 86 | 0.27 | 0.09 | 0.15 | 0.15 | 0.26 | 0.11 | 0.11 | 0.12 | |
| MED | 9211 | 76 | 0.26 | 0.09 | 0.14 | 0.16 | 0.26 | 0.10 | 0.10 | 0.13 | |
| σ | 1992 | 58 | 0.06 | 0.07 | 0.05 | 0.05 | 0.06 | 0.07 | 0.03 | 0.04 | |
| NOAA-18 (20 image pairs) | Min | 3370 | 10 | 0.03 | 0.01 | 0.12 | 0.05 | 0.04 | 0.00 | 0.09 | 0.04 |
| Max | 10,000 | 200 | 0.27 | 0.29 | 0.30 | 0.21 | 0.29 | 0.31 | 0.25 | 0.17 | |
| Mean | 7410 | 94 | 0.13 | 0.13 | 0.20 | 0.15 | 0.13 | 0.14 | 0.15 | 0.12 | |
| MED | 7416 | 96 | 0.13 | 0.12 | 0.20 | 0.15 | 0.13 | 0.13 | 0.15 | 0.12 | |
| σ | 2346 | 69 | 0.06 | 0.07 | 0.05 | 0.03 | 0.07 | 0.08 | 0.04 | 0.04 |
4. Conclusions
The relative, absolute and band-to-band registration (BBR) accuracies of the AVHRR sensor aboard Metop-A, NOAA-17 and NOAA-18 have been investigated in this study. A set of methods, previously developed for similar investigations of MSG-SEVIRI images, was adapted for AVHRR orthoimages. The proposed methods provide fully automated processing, with feature matching at subpixel accuracy. Extensive tests have been performed and problem cases have been visually checked.
Regarding the relative accuracy evaluation of AVHRR, images from Metop-A, NOAA-17 and NOAA-18 acquired on the same day are matched. The results show a good overall accuracy of the system. The mean values of the shifts (x, y) range from 0.1 to 0.4 pixels between all acquisitions from the three different satellites. However, blunders and systematic errors exist in the datasets. Systematic stripes have been observed in the image pairs, where one NOAA-18 (or NOAA-19) image is involved. Mean shifts, which are high and up to 2.9 pixels, have been observed for some days between September and December 2008. The mean of the standard deviations obtained in the matching is higher in comparison to those of SEVIRI HRV images. This can be explained by the greater illumination and atmospheric differences between the images (longer image acquisition intervals), matching images of different satellite sensors and different viewing angles, and existence of local systematic errors.
The lake matching results show in principle both EUMETSAT and GCOS specifications have been fulfilled and the mean values of the 2D shift errors in the images range from 0.1 to 0.4 pixels (Table 11). However, there are blunders in the images, especially in the images from Metop-A, as can be seen from the σ(x, y) values. The shifts are up to 4 pixels for Metop-A, 1.6 pixels in NOAA-17, and 1 pixel in NOAA-18 images in the worst cases. These values are worse than the GCOS accuracy specification. However, all NOAA-18 images stay within the EUMETSAT specification (1 km), whereas one NOAA-17 image exceeds it (1.6 km). For Metop-A images, although the mean (x, y) of shifts are small (−0.1, −0.4 pixels), the standard deviations are high (0.4, 1.1 pixels), which reflects the existence of large shifts, due to unknown reasons, on some days.
The BBR accuracy investigations of AVHRR show that co-registration of Bands 1, 2, 3A and 3B is good (mean shifts usually <0.1 pixels), although in some cases, shifts of up to 0.4 pixels have been observed. The larger shifts are observed in the more off-nadir acquisitions. Regarding Bands 4 and 5, the matching results are very similar for all images, with worse mean shifts of up to about 0.3 pixels, mainly due to more difficult matching.
Some results are summarized in Table 11 and compared to the specifications.
Table 11.
Geometric accuracy specifications and achieved performances for AVHRR. The achieved accuracy values are presented as absolute values.
| Accuracy Type | Specification EUMETSAT/GCOS | Achieved | Achieved | Achieved |
|---|---|---|---|---|
| Metop-A (km) | NOAA-17 (km) | NOAA-18 (km) | ||
| Absolute accuracy | < 1.0 km / < 1/3 pixel | Max(x): 1.7 | Max(x): 1.6 | Max(x): 0.7 |
| Max(y): 4.0 | Max(y): 0.6 | Max(y): 1.0 | ||
| Mean(x): 0.1 | Mean(x): 0.2 | Mean(x): 0.1 | ||
| Mean(y): 0.4 | Mean(y): 0.2 | Mean(y): 0.3 | ||
| σ(x): 0.4 | σ(x): 0.3 | σ(x): 0.2 | ||
| σ(y): 1.1 | σ(y): 0.2 | σ(y): 0.2 | ||
| BBR accuracy | 0.1mrad (~0.1 km) / - | 2&1: Mean(x): 0.07 Mean(y): 0.03 2&3A: Mean(x): 0.17 Mean(y): 0.05 2&4: Mean(x): 0.30 Mean(y): 0.07 2&5: Mean(x): 0.29 Mean(y): 0.09 | 2&1: Mean(x): 0.04 Mean(y): 0.03 2&3A: Mean(x): 0.16 Mean(y): 0.05 2&4: Mean(x): 0.36 Mean(y): 0.09 2&5: Mean(x): 0.27 Mean(y): 0.09 | 2&1: Mean(x): 0.08 Mean(y): 0.09 2&3B: Mean(x): 0.16 Mean(y): 0.13 2&4: Mean(x): 0.17 Mean(y): 0.11 2&5: Mean(x): 0.13 Mean(y): 0.13 |
Overall, the used methods have proven their value. However, the feature distribution in relative and BBR accuracy is not optimal and may require modifications. The blunder detection, as well as fine-tuning of some parameters, could also be improved. The overall work showed that although, in general, accuracies are high and fulfill specifications, errors exceeding the specifications can occur and vary depending on the satellite used, time and location. Sometimes, the reasons for errors, as for the systematic stripes, are not clear. Thus, to ensure the high quality of the image products and information derived thereof, fully automated, comprehensive and permanent quality control procedures should be applied to detect major problems, especially deviations from the specifications.
Acknowledgments
This study is funded by the Swiss GCOS Office at the Federal Office of Meteorology and Climatology MeteoSwiss. The authors thank to Gabriela Seiz and Fabio Fontana from the Swiss GCOS Office and Fabia Hüsler and Stefan Wunderle from the University of Bern for providing data and their valuable inputs. Special thanks to Konstantin Khlopenkov and Alexander Trishchenko for their feedback and Charis Lanaras from ETH Zurich for providing his support during the investigations of the systematic errors.
Author Contributions
Sultan Kocaman Aksakal, Emmanuel Baltsavias and Konrad Schindler provided the overall conception of this research, designed the methodologies and experiments, performed the analysis and wrote the majority of the manuscript. Christoph Neuhaus provided expertise and support regarding the low-level processing as representative of the University of Bern and investigated potential systematic error sources in the original data and software.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Seiz, G.; Foppa, N. National Climate Observing System (GCOS Switzerland). Available online: http://www.gcos.ch (accessed on 15 September 2014).
- Seiz, G.; Foppa, N.; Meier, M.; Paul, F. The role of satellite data within GCOS Switzerland. Remote Sens. 2011, 3, 767–780. [Google Scholar] [CrossRef] [Green Version]
- Global Climate Observing System (GCOS). Available online: http://www.wmo.int/pages/prog/gcos/ (accessed on 14 September 2014).
- World Meteorological Organization (WMO). Systematic Observation Requirements for Satellite-based Products for Climate. Available online: http://www.wmo.int/pages/prog/gcos/Publications/gcos-107.pdf (accessed on 15 September 2014).
- Arnold, G.T; Hubanks, P.A.; Platnick, S.; King, M.D.; Bennartz, R. Impact of Aqua misregistration on MYD06 cloud retrieval properties. In Proceeding of MODIS Science Team Meeting, Washington, DC, USA, 26–28 January 2010.
- EUMETSAT. AVHRR Level 1b Product Guide. Available online: http://oiswww.eumetsat.org/WEBOPS/eps-pg/AVHRR/AVHRR-PG-0TOC.htm (accessed on 21 January 2011).
- Kocaman Aksakal, S. Geometric accuracy investigations of SEVIRI HRV Level 1.5 imagery. Remote Sens. 2013, 5, 2475–2491. [Google Scholar]
- Kocaman Aksakal, S.; Baltsavias, E.; Schindler, K. Analysis of the geometric accuracy of MSG-SEVIRI imagery with focus on estimation of climate variables. In Proceeding of 34th Asian Conference on Remote Sensing, Bali, Indonesia, 20–24 October 2013.
- Kocaman Aksakal, S.; Baltsavias, E.; Schindler, K. Geometric Accuracy Assessment of MSG-SEVIRI Level 1.5 Imagery. Available online: http://www.igp.ethz.ch/photogrammetry/publications/pdf_folder/KocamanEUMETSAT-AMS_13.pdf (accessed on 15 September 2014).
- Khlopenkov, K.; Trishchenko, A.P.; Luo, Y. Achieving subpixel georeferencing accuracy in the Canadian AVHRR processing system. IEEE Trans. Geosci. Remote Sens. 2010, 48, 2150–2161. [Google Scholar] [CrossRef]
- NOAA. Advanced Very High Resolution Radiometer—AVHRR. Available online: http://noaasis.noaa.gov/NOAASIS/ml/avhrr.html (accessed on 15 September 2014).
- Lucas, B.; Kanade, T. An iterative image registration technique with an application to stereo vision. In Proceedings of 7th International Joint Conference on Artificial Intelligence (IJCAI), Vancouver, BC, Canada, 24–28 August 1981; pp. 674–679.
- OpenCV. Available online: http://opencv.org/ (accessed on 30 November 2014).
- Python Imaging Library. Available online: https://pypi.python.org/pypi/PIL (accessed on 30 November 2014).
- Wallis, R. An approach to the space variant restoration and enhancement of images. In Proceedings of Symp. on Current Mathematical Problems in Image Science, Naval Postgraduate School, Monterey, CA, USA, 10–12 November 1976.
- Lanaras, C.; Baltsavias, E.; Schindler, K. A comparison and combination of methods for co-registration of multi-modal images. In Proceedings of 35th Asian Conference on Remote Sensing, Nay Pyi Taw, Myanmar, 27–31 October 2014.
© 2015 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
MDPI and ACS Style
Aksakal, S.K.; Neuhaus, C.; Baltsavias, E.; Schindler, K. Geometric Quality Analysis of AVHRR Orthoimages. Remote Sens. 2015, 7, 3293-3319. https://doi.org/10.3390/rs70303293
AMA Style
Aksakal SK, Neuhaus C, Baltsavias E, Schindler K. Geometric Quality Analysis of AVHRR Orthoimages. Remote Sensing. 2015; 7(3):3293-3319. https://doi.org/10.3390/rs70303293
Chicago/Turabian StyleAksakal, Sultan Kocaman, Christoph Neuhaus, Emmanuel Baltsavias, and Konrad Schindler. 2015. "Geometric Quality Analysis of AVHRR Orthoimages" Remote Sensing 7, no. 3: 3293-3319. https://doi.org/10.3390/rs70303293
