Enhancing land use classification with fusing dual-polarized TerraSAR-X and multispectral RapidEye data

Abstract. The contribution of dual-polarized synthetic aperture radar (SAR) to optical data for the accuracy of land use classification is investigated. For this purpose, different image fusion algorithms are implemented to achieve spatially improved images while preserving the spectral information. To compare the performance of the fusion techniques, both the microwave X-band dual-polarized TerraSAR-X data and the multispectral (MS) optical image RapidEye data are used. Our test site, Gediz Basin, covers both agricultural fields and artificial structures. Before the classification phase, four data fusion approaches: (1) adjustable SAR-MS fusion, (2) Ehlers fusion, (3) high-pass filtering, and (4) Bayesian data fusion are applied. The quality of the fused images was evaluated with statistical analyses. In this respect, several methods are performed for quality assessments. Then the classification performances of the fused images are also investigated using the support vector machines as a kernel-based method, the random forests as an ensemble learning method, the fundamental k-nearest neighbor, and the maximum likelihood classifier methods comparatively. Experiments provide promising results for the fusion of dual polarimetric SAR data and optical data in land use/cover mapping.


Introduction
A wide variety of remote sensing satellite sensors provide data with diverse spectral and spatial resolution for the observation of many phenomena on the Earth.Land use and cover mapping require both the high spectral and spatial resolution for an accurate analysis and interpretation.The data fusion is a key preprocessing method to integrate multisensor and multiresolution images, which has advantages over the result of each individual data set. 1,2Image fusion is an active topic for researching the performance of image fusion techniques for different sensors using both qualitative and quantitative analyses. 3Previous studies in the literature proved that the fusion of the synthetic aperture radar (SAR) and multispectral (MS) images improves the spatial information while preserving the spectral information. 4,5s part of image processing, several approaches of data fusion methods were proposed, and the contributions of the fusion techniques in image classification accuracies were studied. 3,6For different applications of the SAR and MS data fusion, various satellite images were used and various results have been achieved for each study. 4,7,8A generalized intensity modulation for the fusion of MS LANDSAT and ERS-1 SAR images is addressed. 7The benefit of fusion is demonstrated by the maximum likelihood classifier (MLC).In this study, the classification accuracy of vegetation did not improve when using the SAR data, but the discrimination of urban areas was enhanced.The high-resolution spotlight modes of TerraSAR-X and RapidEye are fused using principal component (PC) substitution, Ehlers fusion (EF), Gram-Schmidt (GS), highpass filtering (HPF), modified intensity hue saturation (M-IHS), and Wavelet algorithms. 8In both visual and statistical analyses, the HPF method gave better results.Another study suggest that PC, color normalization, GS, or the University of New Brunswick methods should only be used for a single date and a single sensor dataset. 4Ehlers method is the only one that preserved the spectral information of the MS data, which is suitable for classification. 4Dual-polarized HH-HV (H-Horizontal, V-Vertical) RADARSAT and PALSAR data are fused with LANDSAT-TM data for the comparison of land cover classification. 2The study indicated that among the discrete wavelet transform (DWT), HPF, principal component analysis (PCA), and normalized multiplication methods, only DWT improved the overall accuracy of the MLC.In a previous study, authors used PALSAR and RADARSAT images to fuse them with an Satellite for Observation of Earth (SPOT) image employing five fusion approaches (IHS, PCA, DWT, HPF, and EF). 9 In the results of the fused images, EF and HPF gave satisfactory results for agricultural areas.As a second study, different combinations of the optical-SAR and SAR-SAR fusion results were compared.The EF demonstrated better visual and statistical results. 9his paper extends the previous study of Ref. 10, which focuses on the fusion of the RapidEye data with VV polarized TerraSAR-X SAR data.For the analysis, three pixelbased fusion methods, namely adjustable SAR-MS fusion (ASMF), EF, and HPF were examined.The ASMF method is able to fuse high-resolution SAR data with low-resolution MS data and vice versa.It is possible to use the scaled values for each of the SAR and MS images. 11,12The EF is a hybrid approach, which uses IHS transform with a Fourier domain filtering.A low-pass filter is used to filter the intensity and an inverse high-pass filter is used on high-resolution images. 4In the HPF method, the high-pass filtered high-resolution image is included in each MS band of the low-resolution image. 2,5At the end, visual and statistical analyses were presented and compared.In the previous study, four metrics were used for the statistical analysis.Among all methods, HPF conserved better spectral information.In addition to the previous study, in this study, the VH image of the TerraSAR-X data is also considered.The contribution of the dual-polarized (VV-VH) TerraSAR-X SAR data to RapidEye over agricultural land types is investigated using different image fusion methods.Furthermore, the statistical analyses were extended by adding quality metrics.In total, seven metrics were applied.The bias of the mean (BM) is the difference of the original image and the fused image relative to the original MS data, 13 the difference in variance (DIV) measure is the difference of variance values relative to the original MS data, 13 entropy is a measure that indicates the additional information in the fused image, 13 relative average spectral error (RASE) provides a value for the average performance of the fusion approach, 14 and correlation coefficient (CC) gives the correlation between the original MS image and the fused image.The universal image quality index (UIQI) measures the combination of various factors, such as luminance distortion, contrast distortion, and loss of correlation. 15Relative Dimensionless Global Error in Synthesis-Erreur Relative Globale Adimensionnelle de Synthese (ERGAS) is a global metric which calculates the spectral distortion in a fused image. 16Additionally, Bayesian data fusion (BDF) is also applied for both the TerraSAR-X VV and VH polarized data.BDF was rarely used in the literature.It was performed successfully for high-resolution IKONOS images and recommended as a promising technique for optical/SAR image fusion. 17It is adapted within the ORFEO Toolbox. 18he BDF method allows the user to adjust the images during the fusion processes for emphasizing spectral information via selecting a small weighting coefficient. 17s a part of this study, it is also intended to investigate the contribution of different polarizations of the SAR data to MS data via various fusion techniques for land use image classification.The land use classification needs robust classification methods, which can help in the accurate mapping of land use or land cover classes.There are many studies in the literature with SAR and MS data used separately or together using support vector machines (SVMs) 19 and MLC 20 classification methods.A decision fusion strategy for joint classification of multiple segmentation levels with the SAR and optical data was also evaluated in the literature. 21Ensemble learning and kernel-based classification methods have been confirmed to improve the land use and land cover classification accuracy in remote sensing areas.An assessment of the effectiveness of random forests (RF) for land-cover classification was investigated. 22A comprehensive analysis for the choice of the kernel function and its parameters for the SVMs were presented for the land cover classification. 23In another study, the multiclass SVMs were used and compared with MLC and artificial neural network classifiers for land cover classification. 24The classification of multitemporal SAR and MS data was achieved by the fusion of SVMs. 25 In this study, the original outputs of each discriminant function were used instead of fusing the final classification outputs.
In our proposed study, the contribution of the dual-polarized SAR to optical data for the accuracy of the land use classification is investigated.The fusion methods' effect on classification accuracies are explored by comparing the SVMs as a kernel-based learning method, RF as an ensemble learning method, k-nearest neighbor (k-NN) as a fundamental machine learning classifier, and MLC as a statistical model.
2 Study Area and Data

Test Site
The study area Menemen Plain is located on the west of Turkey in Izmir Province as shown in Fig. 1.The Aegean Sea lies on the west and Izmir Bay on the south, which shapes the border of the study area.The area covers mostly agricultural fields and is approximately 50 km 2 .The crop species depend on the harvesting period and the characteristics of the soil.In this study area, fields were covered with summer crops such as corn, cotton, watermelon, and meadow when the RapidEye data were acquired.There are also some residential areas and small bodies of water in the region.The topographic relief of the study area is lower than 1%, which reduces the effects of topography in image processing.

Data Set
In this experiment, the dual-polarized (VV and VH) TerraSAR-X SAR data and MS RapidEye data were used.The TerraSAR-X image has 8 × 8 m 2 ground resolution, and it was preprocessed to the Enhanced Ellipsoid Corrected product type (e.g., radiometrically enhanced).It was acquired on August 29, 2010, in an ascending pass direction.Data were taken on the Strip Map mode.Detailed specifications of the data sets are given in Table 1.
The RapidEye image was acquired on August 10, 2010, as an L3A format, which was a radiometrically calibrated and orthorectified data resampled to 5 × 5 m 2 ground resolution (WGS 84 Datum UTM projection, zone 35).The RapidEye provides five optical spectral bands, which range between 400 and 850 nm.The RapidEye differs significantly from the standard high-resolution MS satellite sensors (e.g., IKONOS, QuickBird, SPOT), having an extra spectrum called the red edge (690 to 730 nm).For the classification analysis, fieldwork was carried out and crop types were determined using handheld GPSs in agricultural lands on the same date of the RapidEye MS acquisition.The crop types were defined carefully as representing all types of crops in the test site (i.e., for each crop type, ground-truth data were collected from 5 to 15 different fields).

Methodology
In this section, the methodological approach is represented.First of all, preprocessing steps were applied to SAR images and then different fusion methods were applied with RapidEye data.Subsequently, four fusion methods were utilized and the quality assessment was conducted with various quantitative analyses.Lastly, four image classification methods were used to evaluate the contribution of SAR images to the optical images.The workflow is shown in generalized form with Fig. 2.

Preprocessing
Before the application of image fusion methods, image preprocessing steps are necessary.First, the SAR images were filtered using a gamma map filter with 3 × 3 kernel window to reduce speckles.Then, both the VV and VH polarized TerraSAR-X images were registered to the RapidEye image with a less than AE1 pixel root mean square error and resampled to its original pixel size as 8 × 8 m 2 .Since only one optical data are used in the study and the data are not affected severely by atmospheric conditions, an atmospheric correction has not been applied before the fusion process.

Image Fusion Methods
In this study, four different pixel level image fusion approaches have been utilized.They are ASMF, EF, HPF, and BDF.The ASMF method accepts weights for the SAR and MS images separately.In this study, two types of ASMF-based fused images are obtained with different weights.The first image, ASMF-I, is obtained by giving 100% weights to both the SAR and MS images, whereas the second image, ASMF-II, is obtained by giving 50% and 100% weights to the SAR and MS images, respectively.In the Ehlers fused image, the spatial information is improved and spectral characteristics of the MS image are preserved.The HPF method fuses both spatial and spectral information using a band addition approach.In BDF method, the weighting parameter (w) changes between 0 and 1 for the PAN and MS images.In this study, the TerraSAR-X images are first resampled to the same spatial resolution with the MS data.Afterward, two weighting values (0.5 and 0.1) are selected in the fusion processes and fused images, BDF-I and BDF-II, are produced, respectively.

Quality Analysis of Image Fusion Results
The fusion quality assessment is conducted via statistical analyses.A quantitative analysis is applied to evaluate the spectral quality of the fused image by the BM, DIV, entropy, CC, UIQI, RASE, and ERGAS quality metrics.It is expected to be close to zero for BM and DIV metrics.Conversely, while the CC value is close to one, it signifies a better correlation result.Small values of entropy difference between the original MS and the fused image indicate better spectral quality.The higher UIQI values specify better spectral quality in the image.Small values of the RASE and ERGAS indicate a better spectral quality. 16

Image Classification Methods
SVM is very popular and powerful kernel-based learning method.SVMs are introduced as a kernel-based classification algorithm in machine learning society. 26Kernel-based learning aims to separate data in a high-dimensional feature space by mapping the data points with a kernel function.The SVMs are defined for the binary separation problem of samples with n-dimensional feature vector x i and binary class label y i , which can be expressed mathematically as ðx 1 ; y 1 Þ; ðx 2 ; y 2 Þ; : : : ; The SVM creates a decision surface between the samples of the different classes by finding the optimal hyperplane that is closest to the deciding training samples (support vectors).In this way, an optimal classification can be achieved for linearly separable classes.In the linearly inseparable cases, kernel versions of the SVM are defined.The main purpose of the kernel approach in the SVM is to transform the data to  a higher-dimensional space ðϕ∶R n → R h ; h > nÞ, where binary classification can again be achieved linearly. 27The SVM utilizes a kernel function that corresponds to the inner product in the higher-dimensional space.The deciding support vectors can be found with the optimization problem that maximizes the Eq. ( 1) subject to Eq. ( 2): where N denotes the number of training samples, C is the penalty parameter, Kðx u ; x v Þ is the kernel function, and α is the Lagrange multiplier coefficients.
C controls the number of misclassified training samples to allow for a better margin maximization.The SVM does not require the explicit definition of the transformation function ðϕÞ, but is based on the definition of the inner product result in the high-dimensional space as Each nonzero α u value corresponds to a support vector.Given all support vectors ðN SV Þ, the nonlinear classification result can be given in Eq. ( 3) for an arbitrary sample x f ¼ sgn where N SV denotes the number of support vectors and sgn stands for the sign function.
Widely used kernel functions for the SVM can be given as linear kernel Kðx u ; sÞ d with a degree parameter ðdÞ, a scaling factor ðγÞ, and radial basis function kernel with a scaling factor ðγÞ.The multiclass problems can be broken down into several one-against-one binary problems.For an m class problem, a total of mðm − 1Þ∕2 one-against-one SVMs are calculated.The majority vote from the one-against-one classification for each sample decides for the final result. 28F is a supervised ensemble learning technique, which has received highlighted interest in machine learning and pattern recognition society. 29The RF is a powerful "classification and regression tree" (CART) type classifier and it is popular in the remote sensing society with applications on MS and hyperspectral images, as well as kernel-based classification methods.In essence, the RF algorithm builds an ensemble of tree-based classifiers and makes use of bagging or bootstraps aggregating to form an ensemble of CART-like classifiers.In the ensemble of classifiers, each classifier produces a single vote for the appointment of the most frequent class label to the input vector x and fhðx; θ i Þ; i ¼ 1; : : : ; g, where the θ i is independent identically distributed random vectors. 30The random term in the RF refers to the way each tree is trained.Thus, each tree is chosen from a random subset of the features in the training data.In this way, the classifier becomes more robust against the minor variations in input data.The RF also tries to minimize the correlation between the tree classifiers in the ensemble in order for the forest to represent the independent identical distributions in each of the classifiers.Because of the nonparametric structure, high classification accuracy, and feature importance capability, the RF is a promising method in the remote sensing area. 22-NN is a well-known, nonparametric method in CART tasks.The k-NN is an instance-(or sample) based cornerstone classifier in machine learning.It basically predicts labels of the test samples by using labels of the k-nearest training samples in the feature space.Consequently, the majority voting of closest neighbors defines the label which is assigned to the test sample. 31he MLC is a well known and widely used parametric supervised classification method in machine learning and pattern recognition.MLC assumes a multidimensional normal distribution for each class, and computes the probability of a test pixel based on this distribution model for the classification task.In the MLC, the likelihood of a new sample x belonging to a class ω c can be calculated with the following discriminant equation: where x is an n-dimensional data, pðω c Þ is the probability that class ω c occurs in the image, Σ c is the covariance matrix for class ω c , and μ c is the mean vector for class ω c .The key feature of the MLC is the inclusion of the covariance in the normal class distributions. 32,33Results and Discussions

Spectral Quality of Fused Images
For each fused image, seven quality metrics were calculated and an average over five bands were scored, as in Table 2.In the previous work, 10 five bands of each fused image were analyzed separately.The last row of Table 2 shows the ideal image values of each metric.It must be noted that each fusion technique gives different results for each quantitative metric.Due to these different results, each metric was ranked for the fusion algorithms.Then the averages of seven metrics were taken, and Table 3 shows the fusion algorithms ranked by their quality scores.The results of the previous study 10 indicated that comparing five bands of fused images with the original image was not an efficient way of understanding the comparison, due to the problem of having so many bands.We conclude that the TSX VH fused BDF-II image is the best fused image, in that it preserves the spectral information of the RapidEye data better than the other results.In the BDF, the selection of weights plays an important role for the quality of fusion.Assigning a higher weight value decreases the preservation of spectral characteristics for an optical image.
In the previous study, 10 using only the VV image showed that the HPF result of the fusion image gave higher statistical results as compared to the EF and ASMF methods.When comparing the previous results with the recent study, the results show that using the BDF method increased the correlation between the original optical data and the SAR data.Moreover, adding the VH polarized image increased the correlation and introduced the best result among all the fused images using the BDF fusion method.In reference to the polarization of the SAR images, the TSX VV fused images yield slightly better results than the VH fused images (Table 2).The use of more quality metrics and a ranking system demonstrated a better understanding and provided an efficient way of comparison which improved previous results.The best three image fusion results according to quality metrics overall ranking are given in Fig. 3, zoomed into a particular patch for easier inspection.

Classification Results
The class information table of the Menemen test site for seven classes, namely corn type-1, corn type-2, cotton, water, bare soil, artificial structures, and orchards, is shown in Fig.   category on the entire area as shown in Fig. 4. As indicated in Fig. 4, a higher the spatial coverage means a greater number of training samples.The classification results are obtained for the following classifiers representing major classification paradigms in machine learning: the k-NN, RF, SVM, and MLC.In the SVM, linear, polynomial, and radial basis kernel functions (SVM-lin, SVM-poly, and SVM-rbf, respectively) are used.All of the classification experiments are realized in the MATLAB environment.The number of neighbors is chosen as k ¼ 1 and k ¼ 9 for the k-NN classifier.The maximum number of trees is set to 500 in the RF, which is composed of several decision trees.The best parameter optimization of the SVM classifiers is accomplished using the grid search method.The penalty parameter of the SVM C is evaluated between ½1 − 100, with a step size increment of 2 for all kernel function types.In SVM-rbf, the γ parameter of the kernel function is tested between ½0.01 − 10, with a step size increment of 0.1.The polynomial degree d is selected as 3, and the γ parameter of the polynomial kernel is evaluated between ½0.01 − 10, with a step size increment of 0.1.The best parameter settings are obtained by one-against-one multiclass classifier modeling.In the MLC, the classification task has been realized according to the discriminant function given in Eq. ( 4).
Initially, the classification results for the original SAR and MS images are obtained as the baseline results to compare the fused images.In Table 4, the classification accuracies for two polarization bands for the SAR image (TSX VH and TSX VV), five band RapidEye MS images, and a VH-VV polarization layer-stacked SAR image (TSX VV-VH) are given separately.According to the results from Table 4, the VV polarization of the SAR data produces better classification accuracies as compared to the VH polarization.The RapidEye MS data have a classification accuracy up to 94.89% with the SVM-rbf kernel function classifier.
As the second step of the experiment, the classification accuracies of the single polarization SAR and five band MS fused images are given with fusion methods designated for the rows, and classifiers with the columns in Table 5.The overall best classification accuracy is obtained for the ASMF-II_TSX VV with an SVM-rbf as 95.32%.The ASMF-II_TSX VV fused image has the highest classification accuracies of all the proposed classifiers, except the SVM-poly, against other fusion methods.The second most successful case is the ASMF-I_TSX VV fusion, which has the best classification accuracies in the SVM-lin and SVM-poly.The EF_TSX VH fusion method also gives the best result for RF together with the ASMF-II_TSX VV.It must be noted that there is at least one fusion method that surpasses the original RapidEye data in classification accuracy for all the proposed classifiers except RF.It can be concluded that with the SAR and MS fusion, the classification accuracies can be improved as compared to the original MS image.The highest classification accuracies are given in bold in Table 5 for each classification method.According to Table 5, the ASMF-II_TSX VV has the best classification accuracy (95.32%) as compared to Table 4, which includes all the original SAR and MS data.
As the third step of the experiment, the fused images are used in a layer-stacked structure with two SAR polarizations cascaded together as in Table 6.In this way, all the SAR and MS information can be used collectively for the scene.The highest classification accuracies are given in bold in Table 6 for each classification method.It can be seen from Table 6 that the layer-stacked   fusion results reveal higher accuracies as compared to their single polarization counterparts as in Table 5.The only exception to our conclusion is the best SVM-rbf results for Tables 5 and 6.In these tables, equal results (95.32%) are obtained for the ASMF-II_TSX VV in Table 5 and for the HPF_TSX VV-VH layer-stacked fusion in Table 6.The highest classification accuracy for the test site is obtained by the ASMF-II_TSX VV-VH and BDF-II_TSX VV-VH layer-stacked fusion images with 95.74% in Table 6.

Comparison of Results
The statistical significance of the fusion methods' accuracies compared to the RapidEye classification accuracy is evaluated by McNemar's test for each classification method.A contingency table is formed based on the agreements and disagreements of the RapidEye classification results and the evaluated fusion method's classification results.This approach enables us to distinguish the fusion methods' statistical significance over RapidEye classification accuracies.McNemar's test score is evaluated as in Eq. ( 5) according to the contingency matrix given in Table 7 McNemar 0 s χ The results of McNemar's test should be compared to the χ 2 table value for 1 degree of freedom. 34McNemar's test scores for the classification accuracy of each fusion method as compared to the RapidEye classification accuracy per classification method is given in Table 8.The statistically significant results that yield better accuracies as compared to the RapidEye results, with a 90% confidence, are given in italic and bold, whereas the better accuracies compared to the RapidEye results, with a 95% confidence, are given in bold only in Table 8.
The classification maps enable better visual assessment of the classification accuracies for all the fusion methods.Hence, the classification maps are generated based on the predictions of the whole scene with the trained classifier models.The corresponding classification maps for each fusion method that produces the highest classification accuracy are generated.In Fig. 5, the classification maps for the TSX VH, TSX VV, and RapidEye are given, respectively, according to the obtained training models in Table 4. Since the SAR images are heavily contaminated with speckle noise, the whole scene classification results produced a relatively noisy map.On the other hand, the classification map of the whole scene with the RapidEye data shows mostly homogeneous regions, whereas some boundary regions are lost.
The classification maps for the single polarization fused images are shown in Fig. 6.The first column of Fig. 6 provides the results for two types of the ASMF fusion for the MS and different SAR polarizations.The second column of Fig. 6 shows the results for both types of the BDF fusion for the MS and different SAR polarizations.The third column of Fig. 6 presents the results for the EF and HPF fusion methods for the MS and different SAR polarizations.If the results are examined, it can be seen that the MS image fusion with the VH polarized SAR data forms more salt and pepper noise-like labeling in all the fusion methods.The visual inspection of Fig. 6 also shows that the BDF type fusion causes deficiencies in the boundary regions.According to Fig. 6, it can also be concluded that the EF spreads the fused features spatially, so the resulting labels are highly affected by the neighboring labels.
In Fig. 7, the classification maps for the layer-stacked fused images for all the SAR polarizations are presented.The first column in Fig. 7 provides the results for the layer-stacked ASMF fusion images.The second column shows the results for the layer-stacked BDF fusion images, and finally, the third column gives the classification maps for the EF and HPF layer-stacked images.
In Turkey, agricultural statistical data such as yield and acreage are collected by the local technical staff of the Ministry of Agriculture and based on the declarations of farmers.The government gives subsidies to farmers based on declared crop types, yield, and acreage.These data sometimes could be considered as unreliable due to false declarations from farmers.Therefore, environmental decision makers and local authorities will benefit from the detailed information regarding the crop pattern that can be used for strategic planning and sustainable management of agricultural resources.
RapidEye data itself are quite successful; however, there is still the contribution of TerraSAR-X to land use/cover classification.This is an indication that when the RapidEye data are not available (e.g., cloudiness, temporal resolution, etc.) the fusion of TerraSAR-X with any optical imagery might be still considered for land use/cover classification with a higher accuracy.

Conclusion
In this study, the influence of the fusion techniques for the dual-polarized SAR and MS images on the classification performance is presented.Several fusion techniques are utilized in the experiments using the microwave X-band dual-polarized TerraSAR-X data and the MS optical image RapidEye data of the Menemen (Izmir) Plain data.The classification performances are investigated using the SVM, RF, k-NN, and MLC methods in a comparative manner.
The study shows some insight into the improvement of the individual sensor results of land use/cover monitoring using the interoperability of multisensor data.Within this context, even though not all of the fusion methods improved the classification result using the single polarimetric SAR data, the study provides promising results of the fusion of dual polarimetric SAR data and the optical data for the mapping of land use/cover types.The results also confirm that using both dual polarimetric SAR data and MS fused images yielded higher classification results as compared to the single polarimetric SAR data fused and the original images.
Quality metrics may give different results which can cause misunderstanding of results.A ranking method for different quality metrics is more efficient than using them individually.The ranking method could be used with any other data set and region, which makes it a universal application for the assessment of fusion quality.
A single fusion technique is inadequate to improve the image quality for land use mapping.Different image fusion techniques should be compared.For the contribution of SAR characteristics to MS images, we recommend the use of image fusion approaches which are developed to combine radar and the MS dataset.In this study, not only a single polarization of an SAR fused image, but also a stack of dual polarization of data were applied to investigate the individual (VV and VH fused images separately) and dual polarization data (VV and VH fused images together) performances.Using the stack of fused dual-polarized SAR data gave better classification accuracies than both the original MS RapidEye data and the single polarized SAR fused data.
It is suggested that quality metrics should not be the only way for the interpretation of fused images.Image classification should be applied on all fused images although statistical analysis of an image that has a worse result can give the highest classification accuracy.Although fused images' classification accuracies are slightly lower than the accuracy of RapidEye (except the one given in Table 5), it is concluded that VV polarized images have a higher accuracy than VH polarized images.
Using this dataset on an agriculturally dominant area, it is concluded that some image fusion methods performed better than others and improved the results.The selection of the fusion approach and the classification method can play an important role.Although we suggest this appropriate methodology for other study areas, generalizations have not been performed yet for other applications.As further research, we will investigate merging full polarimetric SAR data with optical data to explore the contribution of polarimetry to other remote sensing applications.We planned to test these fusion algorithms for other topographically heterogeneous areas to determine the contributions of fusion techniques.

4 .
The class of artificial structures represents residential areas and other constructions such as roads and airports.Corn type-1 and 2 represent the crop being in one of the two distinct biological periods that is the result of different planting times.The orchards class represents different kinds of planted trees.The ground-truth data, which are required for training and reference/test samples, have been collected by systematic sampling based on ground control points by the fieldwork.The number of training samples of each class category is based on the ratio of the coverage of the class

Table 1
Specifications of data set.

Table 2
Statistical results of quality metrics.

Table 3
Rank values of fused images.

Class Names and Map Colour Number of Fields Ground-truth Collected Number of Pixels Used as Samples
Class information table for seven classes in the Menemen data.

Table 8
McNemar's test scores for difference of agreement compared to RapidEye classification, for fused and layer-stacked fused images, χ 2 critical values at 1 degree of freedom for 90% confidence: 2.71, for 95% confidence: 3.84.