# Ship Classification with High Resolution TerraSAR-X Imagery Based on Analytic Hierarchy Process

Authors: Zhi Zhao; Kefeng Ji; Xiangwei Xing; Wenting Chen; Huanxin Zou
Journal: International Journal of Antennas and Propagation
 (2013)

---

## Abstract

Ship surveillance using space-borne synthetic aperture radar (SAR), taking advantages of high resolution over wide swaths and all-weather working capability, has attracted worldwide attention. Recent activity in this field has concentrated mainly on the study of ship detection, but the classification is largely still open. In this paper, we propose a novel ship classification scheme based on analytic hierarchy process (AHP) in order to achieve better performance. The main idea is to apply AHP on both feature selection and classification decision. On one hand, the AHP based feature selection constructs a selection decision problem based on several feature evaluation measures (e.g., discriminability, stability, and information measure) and provides objective criteria to make comprehensive decisions for their combinations quantitatively. On the other hand, we take the selected feature sets as the input of KNN classifiers and fuse the multiple classification results based on AHP, in which the feature sets’ confidence is taken into account when the AHP based classification decision is made. We analyze the proposed classification scheme and demonstrate its results on a ship dataset that comes from TerraSAR-X SAR images.


---

## Body

## 1. Introduction

With the rapid development of worldwide cooperation in ocean economics and military affairs, ship surveillance has attracted many countries’ attention. Among various sensors, space-borne synthetic aperture radar (SAR) is considered to be an operational choice for its high resolution over wide swaths and all-weather working capability. Ship surveillance based on space-borne SAR has been widely applied to ocean traffic monitoring, environmental protection, antipiracy operation [1], and so forth. There exist many maritime surveillance systems and projects including ocean monitoring workstation (OMW), IMPAST (improving fisheries monitoring through integrating passive and active satellite based technologies), and DECLIMS (detection and classification of marine traffic form space) [2]. Ship detection in space-borne SAR [3, 4] has been investigated with great progress. Although SAR automatic target recognition (ATR) techniques [5–8] have been investigated a lot, ship classification still has a long way to go. With the rapid development of space-borne SAR, the research on ship classification using high resolution SAR imagery attracts much attention.In the recent work on SAR ship classification, how to extract features that describes ship accurately and effectively is emphasized. Knapskog et al. [9] investigated the classification of ships based on 3D model and the contour feature. Margarit and Mallorqui [10] extracted scattering center feature for ship classification. They also applied localized RCS features to ship classification and used ships’ structural information with ENVISAT data [2, 11]. Margarit and Tabasco [12] presented a new ship classification methodology that uses single-pol synthetic aperture radar (SAR) images to categorize targets based on a fuzzy logic decision rule, via radar signature isolation, parametric vector estimation, and decision steps. Premaratne and Safaei [13] extracted superstructure moment invariants feature from ISAR images for ship classification. Pilcher and Khotanzad [14] classified ships by nonlinear joint classifiers based on multiple feature including length, centroid, and peak value features. Touzi et al. [15] described ships based on the permanent symmetric scatters to achieve better classification results. Margarit et al. [16] investigated ship scattering in polarimetric SAR for an improved classification under high clutter conditions and evaluated the potentialities of polarimetric ship scattering for classification methods that provide reasonable performance within cluttered scenes. To utilize the polarimetric interferometric SAR capabilities, they combined Pauli CTD with interferometric-based height retrieval to provide an accurate 3D representation of the ship geometry in terms of a set of persistent polarimetric scatterers [17]. Wu et al. [18] analyzed three types of features from fine quad polarization Radarsat-2 SAR image to characterize the ships and prepare for ship classification. Lakshmi et al. [19] analyzed the polarimetric SAR data for identification and characterization of ships, and the polarimetric signature of interaction of the metallic structure of the ship with water is clearly demarcated for the first time through a particular scattering mechanism representing a dipole scatterer. In the related researches, multifeature is much preferred than single feature and achieves better performance. However, the state-of-the-art feature selection is usually based on the experiences or rough evaluations, and lacks valid comprehensive evaluations.Feature selection is an essential and fundamental problem for classification. Different feature always leads to various classification performances. To choose an optimal feature subset is significant for classification. There are kinds of feature selection methods (i.e., filter and wrapper). State-of-the-art evaluation criteria for filter methods [20–22] include discriminability measure, stability measure, and information measure. Single criteria based feature selection can be easily implemented, but it lacks quantitative and objective evaluation for multifeature combination. Therefore, it is necessary to establish objective and comprehensive evaluation criteria for feature selection.Analytic hierarchy process (AHP) is a structured technique for organizing and analyzing complex decisions. In this paper, we firstly introduce AHP to feature selection. It accounts for the comprehensive performances of multiple evaluation measures and provides valid and objective criteria to make comprehensive decisions for their combinations quantitatively. Moreover, similar to the feature selection procedure, we also introduce AHP to the classification decision to fuse the results of different feature sets. The contribution of the proposed AHP based ship classification scheme lies on two aspects: (1) a feature selection method is derived based on the AHP; (2) a classification decision rule is proposed based on AHP to fuse the classification results of different selected feature sets. The proposed classification scheme has been tested on a ship dataset that comes from high resolution TerraSAR-X SAR images. Preliminary results illustrate that the feature selection method is effective and the classification scheme achieves a ratio of positive classifications around 87.7%.This paper is organized as follows. Section2 provides a brief review of the analytic hierarchy process. Section 3 includes the description of the classification scheme, with special concern on the procedure of feature selection and classification decision based on AHP, respectively. Then, Section 4 validates the method with high resolution TerraSAR-X SAR images, including the detail implementation of AHP for feature selection and classification decision. Finally, the paper is concluded in Section 5.


## 2. Analytic Hierarchy Process

Analytic hierarchy process (AHP) [23–30] is a structured technique for organizing and analyzing complex decisions. Since presented by Saaty in the 1970s, it has been applied in a wide variety of decision situations. The AHP provides a comprehensive and rational framework for structuring a decision problem, for representing and quantifying its elements, for relating those elements to overall goals, and for evaluating alternative solutions. When a hierarchy is built, the decision-makers systematically evaluate various elements by comparing them to one another two at a time, with respect to their impact on an element above them in the hierarchy.AHP has unique advantages when important elements of the decision are difficult to quantify or compare. It is applicable in the condition that there exists multiple-goal decision on evaluating the schemes based on multiple influencing factors. When AHP is applied, four premises need to be satisfied: the elements on every layer are known and the hierarchy is clear; the elements on the same layer are uncorrelated; the elements on the lowest layer could be quantitated; the relations among different elements on different layers need to be definite. The procedure for the AHP can be summarized as follows.(i) Build the hierarchy model. Classify the correlative factors into three or more layers: Top layer (the purpose to solve problem); middle layer (the measures and criteria to use, also called strategy layer or criteria layer), whose number could be more than one; bottom layer (the schemes used to resolve problems).Constitute the pair-wise comparison matrix for the middle layer based on the pair-wise comparison method and 1–9 comparison measure. Priorities among the elements of the hierarchy are achieved by making a series of judgments based on pair-wise comparisons. When comparing the importance of theith factor with the jth factor on the middle layer, it is described by the quantitative relative weighting aij (1~9 and their reciprocals). 1~9 represents different levels for importance; see Table 1. Assume there are M factors to be compared for the middle layer; the pair-wise comparison matrix is AM×M=(aij)M×M.Table 1
Proportional measurements.
Importance evaluation
Coequal importance
Little importance
More importance
Strong importance
Extreme importance
Median value
Numerical value
1
3
5
7
9
2, 4, 6, 8(ii) Calculate the weight vector ofA and check the consistency of the judgments. Calculate the maximum eigenvalue and corresponding eigenvector of A, and then check the consistency. If the check succeeds, the normalized eigenvector is assigned as the weight vector s. If not, the comparison matrix needs to be constituted again. In detail, the steps are illustrated below. Firstly, calculate the indexCIto assess the nonconsistency of pair-wise comparison matrix A:
(1)CI=λmax(A)-MM-1.
And, then, calculate the average random consistency indexRI. It is acquired by building pair-wise comparison matrix A randomly to get the average of maximum eigenvalues with large enough random samples. It is validated to be related to matrix rank numberMonly (see Table 2).Table 2
Standard value of the average random consistency index RI.
Matrix rank number
1
2
3
4
5
6
7
8
9
10
RI
0
0
0.58
0.90
1.12
1.24
1.32
1.41
1.45
1.49Thirdly, calculate the random consistency ratioCR:
(2)CR=CIRI.
WhenCR<0.1, the consistency of matrix A is satisfied. Otherwise, A needs to be reconstructed.(iii) Corresponding to each element of the middle layer, a pair-wise comparison matrix for the bottom layer could be derived. So there areMpair-wise comparison matrix for the bottom layer. Assume that the number of the elements for the bottom layer isN; theN×Npair-wise comparison matrix could be expressed as Bi(i=1,2,…,M). And then, check the consistency of them using (1) and (2). If the check succeeds (CR<0.1), the maximum eigenvalue and corresponding eigenvector are derived for each of them. The weight vectors are calculated by normalizing the maximum eigenvalue vectors. So there arexweight vectors vi(i=1,2,…,M), and they could be concatenated into a weight-vector matrix VM×N.(iv) Come to a final decision based on the results of this process. There is a need needs to calculate the total scores for every scheme of bottom layer and give the orders of them. The final decision is made based on the order ofs·V.


## 3. Ship Classification Based on AHP

In this section, we present the proposed ship classification scheme, with special concern on the AHP based feature selection and classification decision procedures.


### 3.1. Classification Scheme

The proposed scheme for ship classification is illustrated in Figure1. The basic process includes four main parts: feature extraction, feature selection, classifier training and classification decision. The kernel idea of the proposed model is to apply analytic hierarchy process to the feature selection and classification decision. In this paper, the evaluation for feature selection adopts filter measures. Ship classification is mainly based on the single classifier and takes the KNN classifier as an example to illustrate the problem. It is worth noting that the feature extraction and classifier training have been widely investigated [31–33] and this paper will not pay additional attention to these processes.Figure 1
Ship classification model.


### 3.2. Feature Selection Based on AHP

Features derived from ship slices are used as the input of classifier for ship classification. Considering the same classifier, different feature set always leads to different classification result. The feature selection is necessary to improve the performance. The hierarchy model for feature selection is presented in Figure2.Figure 2
Hierarchy model for feature selection.Assume that the evaluation measure vector is(e1,e2,…,eM). The weightiness of them could be expressed as(g1,g2,…,gM). The relative weightinessaijcould be expressed as
(3)aij=φ(gigj),i,j∈{1,2,…,M},
whereφ(·)is the mapping from accuracy rate domain{gi/gj}to the discrete weighting domain illustrated in Table 1. So the pair-wise comparison matrix is:
(4)A=[1a12a13⋯a1Ma211a23⋯a2Ma31a321⋯a3M⋮⋮⋮1⋮aM1aM2aM3⋯1]=ϕ([1g1g2g1g3⋯g1gMg2g11g2g3⋯g2gMg3g1g3g21⋯g3gM⋮⋮⋮1⋮gMg1gMg2gMg3⋯1]).After calculation for the maximum eigenvalue of matrixA, the check of consistency can be performed through comparison betweenCRand 0.1. In addition, the maximum eigenvalue vector can also be derived:
(5)s*=(m1,m2,m3,…,mM).It needs to be normalized to a standardization vectors with the conditions that
(6)mi>0,∑imi=1.As to one of the evaluation measures, the pair-wise comparison matrix for the scheme layer is(7)Bi=[1b12b13⋯b1Nb211b23⋯b2Nb31b321⋯b3N⋮⋮⋮1⋮bN1bN2bN3⋯1],i∈{1,2,…,M}.The weight vectors of{Bi,i=1,2,3,…,M}can be calculated by the normalization of the maximum eigenvalue vectors:
(8)v1=(x1(1),x2(1),x3(1),…,xN(1))v2=(x1(2),x2(2),x3(2),…,xN(2))⋮⋮vM=(x1(M),x2(M),x3(M),…,xN(M)).And then{vi,i=1,2,…,M} constitutes a matrix V:
(9)V=(v1,v2,v3,…,vM)T.The final decision on feature selection can be made through comparisons of total scores of the features to be chosen:(10)Sort(s·V).


### 3.3. Classification Decision Based on AHP

Because the evaluation measures including discriminability, stability, and information are calculated based on certain amount of samples (nonsufficient theoretically), the errors are unavoidable. According to the calculated orders of features, feature selection needs to consider the error tolerance. For example, if the optimal feature set is made up of five features, six or seven features (the former six or seven of the descending order) need to be taken into account, so there areC65 or C75 feature sets {FM} that should be considered.To set the feature sets derived above as the input for KNN classifier, respectively, multiple classification measurement value vectors (not the final labels) could be acquired. And then the final classification decision is made using the AHP based comprehensive processing of these measurement value vectors.Firstly, construct the hierarchy model. When a ship is classified into three kinds (e.g., cargo ship, oil tanker, and container ship) based on the KNN classifier, feature sets{FM} make up the evaluation evidence (criteria layer). The hierarchy model is shown in Figure 3. Unlike the feature selection, the pair-wise comparison matrix for the scheme layer is constructed based on the classification measurement value vectors of a feature set, neither the value of the feature set itself nor the final classification labels (cargo ship, oil tanker, or container ship).Figure 3
Hierarchy model for classification decision based on single classifier.Secondly, the pair-wise comparison matrix can then be constructed based on the comparisons of different feature sets’ contributions. Due to the stochastic errors, the feature sets could be considered with the same importance (the comparison matrix of all ones). Another method to determine the weighting among the feature sets is based on the feature’s priority in the feature set. For example, there are two feature sets{f1,f2,f3,f4,f5}and{f1,f2,f3,f4,f6}; assume that the feature priorities are f1>f2>f3>f4>f5>f6; then the first feature set is judged to be prior to the second. The priorities for multiple feature sets could be transformed into their confidence, so aij is the confidence ratio between feature sets i andj. The comparison matrix AM×Mcould be constituted. After calculation for the maximum eigenvalue of matrix AM×M, the check of consistency can be performed through comparison betweenCRand 0.1. In addition, the maximum eigenvalue vector can also be derived:
(11)s*=(m1,m2,m3⋯mM).It also needs to be transformed into a standardization vectors. There are three schemes (type I: cargo ship, type II: oil tanker, and type III: container ship) for classification under every feature set, so the Mcorresponding pair-wise comparison matrix Bi,i∈{1,2,…,M}, could be produced:
(12)IIIIIIIIIIII[1b21b31b121b32b13b231].The weight coefficients are calculated by(13)bij=φ(gigj),i,j∈{1,2,3},
where gi is the measurement value of the classifier’s output.Take the KNN classifier [34] as an example; when a ship chipxis classified, the discriminant of the ith type of ship is
(14)di(x)=ki,
wherekiis the number ofKnearest neighbor samples which belongs to the ith type of ship. So the measurement value is
(15)gi=di(x).After checking it consistency, if succeeded, the weight vectors matrixV={viT}ofBcould be calculated. The total scores of the three types are calculated by the products of the weight vector of the middle layer pair-wise comparison matrix and the weight vector of the bottom layer pair-wise comparison matrix. The final decision on classification for a single classifier based on multiple feature sets can be made through comparisons of the total scores:
(16)Sort(s·V).


## 3.1. Classification Scheme

The proposed scheme for ship classification is illustrated in Figure1. The basic process includes four main parts: feature extraction, feature selection, classifier training and classification decision. The kernel idea of the proposed model is to apply analytic hierarchy process to the feature selection and classification decision. In this paper, the evaluation for feature selection adopts filter measures. Ship classification is mainly based on the single classifier and takes the KNN classifier as an example to illustrate the problem. It is worth noting that the feature extraction and classifier training have been widely investigated [31–33] and this paper will not pay additional attention to these processes.Figure 1
Ship classification model.


## 3.2. Feature Selection Based on AHP

Features derived from ship slices are used as the input of classifier for ship classification. Considering the same classifier, different feature set always leads to different classification result. The feature selection is necessary to improve the performance. The hierarchy model for feature selection is presented in Figure2.Figure 2
Hierarchy model for feature selection.Assume that the evaluation measure vector is(e1,e2,…,eM). The weightiness of them could be expressed as(g1,g2,…,gM). The relative weightinessaijcould be expressed as
(3)aij=φ(gigj),i,j∈{1,2,…,M},
whereφ(·)is the mapping from accuracy rate domain{gi/gj}to the discrete weighting domain illustrated in Table 1. So the pair-wise comparison matrix is:
(4)A=[1a12a13⋯a1Ma211a23⋯a2Ma31a321⋯a3M⋮⋮⋮1⋮aM1aM2aM3⋯1]=ϕ([1g1g2g1g3⋯g1gMg2g11g2g3⋯g2gMg3g1g3g21⋯g3gM⋮⋮⋮1⋮gMg1gMg2gMg3⋯1]).After calculation for the maximum eigenvalue of matrixA, the check of consistency can be performed through comparison betweenCRand 0.1. In addition, the maximum eigenvalue vector can also be derived:
(5)s*=(m1,m2,m3,…,mM).It needs to be normalized to a standardization vectors with the conditions that
(6)mi>0,∑imi=1.As to one of the evaluation measures, the pair-wise comparison matrix for the scheme layer is(7)Bi=[1b12b13⋯b1Nb211b23⋯b2Nb31b321⋯b3N⋮⋮⋮1⋮bN1bN2bN3⋯1],i∈{1,2,…,M}.The weight vectors of{Bi,i=1,2,3,…,M}can be calculated by the normalization of the maximum eigenvalue vectors:
(8)v1=(x1(1),x2(1),x3(1),…,xN(1))v2=(x1(2),x2(2),x3(2),…,xN(2))⋮⋮vM=(x1(M),x2(M),x3(M),…,xN(M)).And then{vi,i=1,2,…,M} constitutes a matrix V:
(9)V=(v1,v2,v3,…,vM)T.The final decision on feature selection can be made through comparisons of total scores of the features to be chosen:(10)Sort(s·V).


## 3.3. Classification Decision Based on AHP

Because the evaluation measures including discriminability, stability, and information are calculated based on certain amount of samples (nonsufficient theoretically), the errors are unavoidable. According to the calculated orders of features, feature selection needs to consider the error tolerance. For example, if the optimal feature set is made up of five features, six or seven features (the former six or seven of the descending order) need to be taken into account, so there areC65 or C75 feature sets {FM} that should be considered.To set the feature sets derived above as the input for KNN classifier, respectively, multiple classification measurement value vectors (not the final labels) could be acquired. And then the final classification decision is made using the AHP based comprehensive processing of these measurement value vectors.Firstly, construct the hierarchy model. When a ship is classified into three kinds (e.g., cargo ship, oil tanker, and container ship) based on the KNN classifier, feature sets{FM} make up the evaluation evidence (criteria layer). The hierarchy model is shown in Figure 3. Unlike the feature selection, the pair-wise comparison matrix for the scheme layer is constructed based on the classification measurement value vectors of a feature set, neither the value of the feature set itself nor the final classification labels (cargo ship, oil tanker, or container ship).Figure 3
Hierarchy model for classification decision based on single classifier.Secondly, the pair-wise comparison matrix can then be constructed based on the comparisons of different feature sets’ contributions. Due to the stochastic errors, the feature sets could be considered with the same importance (the comparison matrix of all ones). Another method to determine the weighting among the feature sets is based on the feature’s priority in the feature set. For example, there are two feature sets{f1,f2,f3,f4,f5}and{f1,f2,f3,f4,f6}; assume that the feature priorities are f1>f2>f3>f4>f5>f6; then the first feature set is judged to be prior to the second. The priorities for multiple feature sets could be transformed into their confidence, so aij is the confidence ratio between feature sets i andj. The comparison matrix AM×Mcould be constituted. After calculation for the maximum eigenvalue of matrix AM×M, the check of consistency can be performed through comparison betweenCRand 0.1. In addition, the maximum eigenvalue vector can also be derived:
(11)s*=(m1,m2,m3⋯mM).It also needs to be transformed into a standardization vectors. There are three schemes (type I: cargo ship, type II: oil tanker, and type III: container ship) for classification under every feature set, so the Mcorresponding pair-wise comparison matrix Bi,i∈{1,2,…,M}, could be produced:
(12)IIIIIIIIIIII[1b21b31b121b32b13b231].The weight coefficients are calculated by(13)bij=φ(gigj),i,j∈{1,2,3},
where gi is the measurement value of the classifier’s output.Take the KNN classifier [34] as an example; when a ship chipxis classified, the discriminant of the ith type of ship is
(14)di(x)=ki,
wherekiis the number ofKnearest neighbor samples which belongs to the ith type of ship. So the measurement value is
(15)gi=di(x).After checking it consistency, if succeeded, the weight vectors matrixV={viT}ofBcould be calculated. The total scores of the three types are calculated by the products of the weight vector of the middle layer pair-wise comparison matrix and the weight vector of the bottom layer pair-wise comparison matrix. The final decision on classification for a single classifier based on multiple feature sets can be made through comparisons of the total scores:
(16)Sort(s·V).


## 4. Experiments and Analysis

### 4.1. Data Set

To validate the proposed methods, the experiments are carried out using TerraSAR-X SAR imagery with VV polarization, azimuth resolution 3.3 meters, and ground range resolution 1.9 meters [35]. Ship chips extracted from the imagery make up the experiment database.In SAR images, the appearance of ship varies with many factors, such as atmospheric phenomena, image resolution, ship speed, ship orientation, ship size, and the material of the ship. In this paper, the high resolution TerraSAR-X images are captured from the near coastal areas; the sea state is considered as stable. For simplification, the atmospheric effects are neglected. For the cases of nonmetallic or small ship, the ship scattering is not clear and difficult to be identified even by experts. It is a challenging work for the state-of-the-art methods to recognize these kinds of ships. Additionally, the ships with severe smearing and occlusion that result from ship cruising and various orientations are discarded when the data set is constituted. Therefore, we focus on the metallic ships with relative large dimensions and clear appearances on the SAR images.The number of TerraSAR-X ship chips collected is 286. Ship classification was firstly done by the human operator. There are 56 oil tanker ships, 160 cargo ships, and 70 container ships, respectively. Figure4 shows two examples of each kind. We choose 60% of each kind of ships as the training set and the left 40% as the test set, whereas the training set includes 96 cargo ship samples, 42 container ship samples, and 34 oil tanker ship samples, and the test set is composed of 64 cargo ship samples, 28 container ship samples, and 22 oil tanker ship samples.Figure 4
Samples of ship chips.


### 4.2. Feature Selection

In this experiment, 21 features [31–33] were extracted from each ship from the database (N=21). They include lengthf1, widthf2, length to width ratiof3, areaf4, perimeterf5, shape complexityf6, centroidf7, moment of inertiaf8, massf9, the mean intensityf10, variance coefficientf11, weighted-rank fill ratiof12, standard deviationf13, fractal dimensionf14, and Hu momentf15~f21. These features belong to the geometrical structural features and gray intensity statistical features.In addition, electromagnetic scattering features are also important to ship classification, but they are normally difficult to extract and require much higher image resolution. They increase the information quantity using both the magnitude and phase information of the complex SAR image, relate closely to the ships’ scattering characteristics, and reflect the targets’ essential information. In this paper, in order to illustrate the problem concisely, we only choose the 21 features presented above.Three evaluation measures including feature discriminability, feature stability, and BIF (best individual feature) are considered in this experiment. Correlation measure is one of the important measures for feature selection to reduce the redundancy. The paper selects three measures just to illustrate the problem, and actually more measures could be adopted including the correlation measure. The definition of the features and the evaluation measures are beyond the scope of this paper, and they can be referenced in the literatures [20–22].Firstly, the comparison matrix for the three evaluation measures could be set according to the experiences. If it is not easy to decide the weighting, it could also be set as all ones (with the same importance). In this paper, it is set by experience. Discriminability (measure I) is considered as the most important one, and stability (measure II) is more important than BIF (measure III). SoM=3, and A could be set as
(17)IIIIIIIII13III123[1131512331131312353131123].And then, the maximum eigenvalue ofA is λmax(A)=3.0554, CI=(λmax(A)-n)/(n-1)=0.0277. We look up Table 2 to get the value of RI(RI=0.58), soCR=CI/RI=0.0478<0.1, and the consistency of A is satisfied.The standardization form of the maximum eigenvector is(18)s=[0.6054,0.2915,0.1031]T.Based on the training samples, the three evaluation measures are calculated and the results are illustrated in Figure5, respectively.Analysis of feature attribute. (a) Discriminability (the ratio of intraclass distance and interclass distance), (b) stability (normalized variance coefficient), and (c) BIF (mutual information).
(a)
(b)
(c)According to the results above, the comparison matrix for discriminability, stability, and BIF could be derived:Bi|21×21,i∈{discriminability,stability,BIF}. And then, the weighting vectors of them are acquired:
(19)vdiscriminaloility=[0.0747,0.0305,0.0765,0.0623,0.0676,0.0596,0.0305,0.0623,0.0605,0.0489,0.0182,0.0114,0.0279,0.0308,0.0209,0.0525,0.0311,0.0605,0.0418,0.0640,0.0676]T,vstability=[0.0385,0.0200,0.0696,0.0623,0.0197,0.0623,0.1026,0.0360,0.0385,0.0238,0.0360,0.0360,0.0360,0.0342,0.0281,0.0623,0.1026,0.0733,0.0360,0.0623,0.0197]T,vBIF=[0.0760,0.0472,0.0553,0.0461,0.0461,0.0472,0.0392,0.0461,0.0461,0.0461,0.0403,0.0415,0.0449,0.0438,0.0484,0.0484,0.0472,0.0484,0.0472,0.0461,0.0484]T.The total score for each feature is shown in Figure6.Figure 6
Total score for every feature.Thus the descending order of them isf3>f1>f18>f20>f4>f6>f16>f17>f8>f9>f7>f21>f5>f10>f19>f14>f13>f2>f15>f11>f12.It can be concluded that all the geometrical structural features have better performances except featuref2 (width). Geometrical structural features relate closely to the ships’ size, shape, and structure. In the metallic large ship cases (cargo, oil tanker, and container ship), geometrical structural features perform better theoretically. Three typical types of ships differ from each other in terms of size, shape, and structure, and the results indicate that these features have better performance, especially for feature f3 (length to width ratio) which is important for ship classification. Most of the Hu moment features perform better in the experiment. They describe the characteristics of ships’ gray intensity distribution and edge shapes and reflect the ships’ essential attributes. Gray intensity statistical features including weighted-rank fill ratio f12, standard deviation f13, and fractal dimension f14 describe the regularity and structure information of the pixels’ gray intensity distribution. They are preferable and also widely used in the ship classification; however, the experimental results illustrate that they perform worse than the geometrical structural features as a whole for ship classification of these three types.Feature selection could be implemented according to the results above. The features with higher scores or orders should be considered to make up the optimal feature set, but the number is still unknown. To derive it, 20 feature sets are constructed based on the derived order and ship classification based on the test set is carried out to evaluate the accuracy. The first set is{f3,f1}, adding one feature based on the order once a time, so the second one is{f3,f1,f18}and the third is{f3,f1,f18,f20}. The rest may be deduced by analogy. The KNN based classification result is shown in Figure 7.Figure 7
Classification accuracy for the feature sets.Figure7 illustrates that the accuracy curve presents an ascending-stable-descending trend as a whole (red dot line). The accuracy rises firstly (from 61.4% to 85%) and maintains a relatively stable level (around 85%); then it falls below 80% primarily due to the features’ redundancy. From the fifth to the tenth feature set, the accuracy varies slightly and keeps constant approximately, so the fifth feature set {f3,f1,f18,f20,f4,f6}could be considered as a turning point (intersection of the red solid line and the curve), and it is set as the optimal feature set here (the number of features is six). In practice, the first seven features {f3,f1,f18,f20,f4,f6,f16} are selected to constructC76feature sets to enlarge the scope of tolerance error, and then the derived feature sets are
(20)F1={f3,f1,f18,f20,f4,f6},F2={f3,f1,f18,f20,f4,f16},F3={f3,f1,f18,f20,f6,f16},F4={f3,f1,f18,f4,f6,f16},F5={f3,f1,f20,f4,f6,f16},F6={f3,f18,f20,f4,f6,f16},F7={f1,f18,f20,f4,f6,f16}.
According to the feature scores, the feature sets’ priorities are derived:
(21)F1>F2>F3>F4>F5>F6>F7.


### 4.3. Ship Classification Based on Single Classifier

The experiment is carried out based on the test set of ship chips. When a ship chip is input for classification, the feature sets{Fn,n=1~7}derived above are set as the inputs for KNN classifier (the neighbor parameterK=3in this experiment). In this experiment, M=7, N=3. According to the feature sets’ priorities derived in the feature selection, the pair-wise comparison matrix can be set as
(22)A=[F1F2F3F4F5F6F7F11234567F212123456F3131212345F41413121234F515141312123F6161514131212F71716151413121].
The maximum eigenvalue ofAisλmax(A)=7.2786, CI=(λmax(A)-n)/(n-1)=0.0464. From Table 2, the value ofRIcould be acquired (RI=1.32), and thenCR=CI/RI=0.0352<0.1, so the consistency ofAis satisfied. The standardization form of the maximum eigenvector is
(23)s=[0.3086,0.2369,0.1745,0.1221,0.0803,0.0490,0.0286]T.Then, ship classification using KNN classifier is based on the seven feature sets{Fi,i=1,2,3,4,5,6,7}. As to each feature set, the comparison matrix for the scheme layer is constructed according to Section 3.3. And then, the comparison matrix for other feature sets could be constituted similarly. The corresponding weighting vectors can be calculated. The total scores for the three kinds (cargo ship, container ship, or oil tanker) can be derived and the final decision (the type of input ship chip) is made.When the classification decision of one ship chip is made, others are processed one by one based on the same method. Once the classification decisions of all the ship chips are made, the classification accuracy (for each type and the total accuracy) is calculated by the ratio of correct-classified number to the total number. Classification accuracies of the comparisons between the single-feature-set based classifications and the presented AHP-based classification are shown in Table3. The corresponding diagram of confusion matrix is presented in Figure 8. KNN classification based on the derived feature sets {Fn,n=1~7}could achieve satisfied performance. Inparticular for the AHP based method, the classification accuracy rates of the cargo and container ship come to 89.1% and 89.3%  respectively, and the comprehensive performance is superior to other methods with the total accuracy of 87.7%.Table 3
Ship classification accuracy based on KNN classifier.
Methods
Ship type
Total accuracy
Cargo ship
Oil tanker
Container ship
F
1
87.5%
77.3%
85.7%
85.1%
F
2
87.5%
72.7%
89.3%
85.1%
F
3
85.9%
77.3%
78.6%
82.5%
F
4
82.8%
68.2%
75%
78.1%
F
5
84.4%
72.7%
78.6%
80.7%
F
6
82.8%
68.2%
78.6%
78.9%
F
7
82.8%
63.6%
75%
77.2%
AHP
89.1%
81.8%
89.3%
87.7%Diagram of confusion matrix: (a)F1, (b) F2, (c) F3, (d) F4, (e) F5, (f) F6, (g) F7, and (h) AHP.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)In addition, some cases in which the vessel was wrongly classified (take the oil tanker case as an example) are presented in Table4 and some analysis of the possible causes is also presented. There are three chips whose oil pipe structure could be identified by human operator but not as clear as the ideal ones due to many interferential brighter scattering points around the oil pipe. However, they are all wrongly classified into cargo ships. In these three cases, the third chip is disturbed by strong scattering points badly and none of the chosen feature sets could be classified correctly, so the final decision based on AHP considers 0% possibility for the oil tanker. From the cases for chips No. 1 and No. 2, we can see that onlyF3andF5classify the first one into oil tanker, and F2, F3, F5, andF6classify the second one correctly. Referring to the detailed features for each set shown in Section 4.2, the possible reasons for these two cases are analyzed below.Table 4
Wrong classification cases for oil tanker.
Oil tanker cases
F
1
×
×
×
F
2
×
✓
×
F
3
✓
✓
×
F
4
×
×
×
F
5
✓
✓
×
F
6
×
✓
×
F
7
×
×
×
AHP
×
×
×
✓: identified; ×: unidentified.F
1classifies the first and second one wrongly; however,F3classifies both correctly andF2classifies the second one correctly. The features affecting their performances focus on the features f4 (area),f6 (shape complexity), and f16 (Hu moment). Considering their detailed performances,f4andf6are more sensitive to the classification of these two chips compared to other features. Actually, the interferential scattering points really affect ships’ shape complexity and the accurate extraction of targets’ area. Meanwhile, we can see that featuref16might be more preferable tof4orf6for the oil tanker classification although Section 4.2 shows that the latter two are a little more superior in the overall performance. Concluded from the performance of feature setF5(compared toF4), possibly Hu momentf20is superior to Hu momentf18in oil tanker classification.Actually ship features are due to the ship scattering, and the relations among them could be considered as being nonlinear, not totally independent. Essentially, the feature set including several features contributes to the classification as a whole. It is quite difficult to clarify their independent contribution clearly.When it refers to the second ship chip classification, four of seven feature sets are classified correctly, however,F2, F3, F5, and F6 have their own superiorities (in this experiment they are assigned 1/2, 1/3, 1/5, and 1/6). When the final decision based on AHP is made, the total score of classification into oil tanker is a little less than the score of classification into cargo. Actually, this is only a decision problem, and it just illustrates the possibility of falling within a certain type of ships. When the superiority factors of the feature sets vary proportionally, the results might vary but the conclusions are roughly consistent.


## 4.1. Data Set

To validate the proposed methods, the experiments are carried out using TerraSAR-X SAR imagery with VV polarization, azimuth resolution 3.3 meters, and ground range resolution 1.9 meters [35]. Ship chips extracted from the imagery make up the experiment database.In SAR images, the appearance of ship varies with many factors, such as atmospheric phenomena, image resolution, ship speed, ship orientation, ship size, and the material of the ship. In this paper, the high resolution TerraSAR-X images are captured from the near coastal areas; the sea state is considered as stable. For simplification, the atmospheric effects are neglected. For the cases of nonmetallic or small ship, the ship scattering is not clear and difficult to be identified even by experts. It is a challenging work for the state-of-the-art methods to recognize these kinds of ships. Additionally, the ships with severe smearing and occlusion that result from ship cruising and various orientations are discarded when the data set is constituted. Therefore, we focus on the metallic ships with relative large dimensions and clear appearances on the SAR images.The number of TerraSAR-X ship chips collected is 286. Ship classification was firstly done by the human operator. There are 56 oil tanker ships, 160 cargo ships, and 70 container ships, respectively. Figure4 shows two examples of each kind. We choose 60% of each kind of ships as the training set and the left 40% as the test set, whereas the training set includes 96 cargo ship samples, 42 container ship samples, and 34 oil tanker ship samples, and the test set is composed of 64 cargo ship samples, 28 container ship samples, and 22 oil tanker ship samples.Figure 4
Samples of ship chips.


## 4.2. Feature Selection

In this experiment, 21 features [31–33] were extracted from each ship from the database (N=21). They include lengthf1, widthf2, length to width ratiof3, areaf4, perimeterf5, shape complexityf6, centroidf7, moment of inertiaf8, massf9, the mean intensityf10, variance coefficientf11, weighted-rank fill ratiof12, standard deviationf13, fractal dimensionf14, and Hu momentf15~f21. These features belong to the geometrical structural features and gray intensity statistical features.In addition, electromagnetic scattering features are also important to ship classification, but they are normally difficult to extract and require much higher image resolution. They increase the information quantity using both the magnitude and phase information of the complex SAR image, relate closely to the ships’ scattering characteristics, and reflect the targets’ essential information. In this paper, in order to illustrate the problem concisely, we only choose the 21 features presented above.Three evaluation measures including feature discriminability, feature stability, and BIF (best individual feature) are considered in this experiment. Correlation measure is one of the important measures for feature selection to reduce the redundancy. The paper selects three measures just to illustrate the problem, and actually more measures could be adopted including the correlation measure. The definition of the features and the evaluation measures are beyond the scope of this paper, and they can be referenced in the literatures [20–22].Firstly, the comparison matrix for the three evaluation measures could be set according to the experiences. If it is not easy to decide the weighting, it could also be set as all ones (with the same importance). In this paper, it is set by experience. Discriminability (measure I) is considered as the most important one, and stability (measure II) is more important than BIF (measure III). SoM=3, and A could be set as
(17)IIIIIIIII13III123[1131512331131312353131123].And then, the maximum eigenvalue ofA is λmax(A)=3.0554, CI=(λmax(A)-n)/(n-1)=0.0277. We look up Table 2 to get the value of RI(RI=0.58), soCR=CI/RI=0.0478<0.1, and the consistency of A is satisfied.The standardization form of the maximum eigenvector is(18)s=[0.6054,0.2915,0.1031]T.Based on the training samples, the three evaluation measures are calculated and the results are illustrated in Figure5, respectively.Analysis of feature attribute. (a) Discriminability (the ratio of intraclass distance and interclass distance), (b) stability (normalized variance coefficient), and (c) BIF (mutual information).
(a)
(b)
(c)According to the results above, the comparison matrix for discriminability, stability, and BIF could be derived:Bi|21×21,i∈{discriminability,stability,BIF}. And then, the weighting vectors of them are acquired:
(19)vdiscriminaloility=[0.0747,0.0305,0.0765,0.0623,0.0676,0.0596,0.0305,0.0623,0.0605,0.0489,0.0182,0.0114,0.0279,0.0308,0.0209,0.0525,0.0311,0.0605,0.0418,0.0640,0.0676]T,vstability=[0.0385,0.0200,0.0696,0.0623,0.0197,0.0623,0.1026,0.0360,0.0385,0.0238,0.0360,0.0360,0.0360,0.0342,0.0281,0.0623,0.1026,0.0733,0.0360,0.0623,0.0197]T,vBIF=[0.0760,0.0472,0.0553,0.0461,0.0461,0.0472,0.0392,0.0461,0.0461,0.0461,0.0403,0.0415,0.0449,0.0438,0.0484,0.0484,0.0472,0.0484,0.0472,0.0461,0.0484]T.The total score for each feature is shown in Figure6.Figure 6
Total score for every feature.Thus the descending order of them isf3>f1>f18>f20>f4>f6>f16>f17>f8>f9>f7>f21>f5>f10>f19>f14>f13>f2>f15>f11>f12.It can be concluded that all the geometrical structural features have better performances except featuref2 (width). Geometrical structural features relate closely to the ships’ size, shape, and structure. In the metallic large ship cases (cargo, oil tanker, and container ship), geometrical structural features perform better theoretically. Three typical types of ships differ from each other in terms of size, shape, and structure, and the results indicate that these features have better performance, especially for feature f3 (length to width ratio) which is important for ship classification. Most of the Hu moment features perform better in the experiment. They describe the characteristics of ships’ gray intensity distribution and edge shapes and reflect the ships’ essential attributes. Gray intensity statistical features including weighted-rank fill ratio f12, standard deviation f13, and fractal dimension f14 describe the regularity and structure information of the pixels’ gray intensity distribution. They are preferable and also widely used in the ship classification; however, the experimental results illustrate that they perform worse than the geometrical structural features as a whole for ship classification of these three types.Feature selection could be implemented according to the results above. The features with higher scores or orders should be considered to make up the optimal feature set, but the number is still unknown. To derive it, 20 feature sets are constructed based on the derived order and ship classification based on the test set is carried out to evaluate the accuracy. The first set is{f3,f1}, adding one feature based on the order once a time, so the second one is{f3,f1,f18}and the third is{f3,f1,f18,f20}. The rest may be deduced by analogy. The KNN based classification result is shown in Figure 7.Figure 7
Classification accuracy for the feature sets.Figure7 illustrates that the accuracy curve presents an ascending-stable-descending trend as a whole (red dot line). The accuracy rises firstly (from 61.4% to 85%) and maintains a relatively stable level (around 85%); then it falls below 80% primarily due to the features’ redundancy. From the fifth to the tenth feature set, the accuracy varies slightly and keeps constant approximately, so the fifth feature set {f3,f1,f18,f20,f4,f6}could be considered as a turning point (intersection of the red solid line and the curve), and it is set as the optimal feature set here (the number of features is six). In practice, the first seven features {f3,f1,f18,f20,f4,f6,f16} are selected to constructC76feature sets to enlarge the scope of tolerance error, and then the derived feature sets are
(20)F1={f3,f1,f18,f20,f4,f6},F2={f3,f1,f18,f20,f4,f16},F3={f3,f1,f18,f20,f6,f16},F4={f3,f1,f18,f4,f6,f16},F5={f3,f1,f20,f4,f6,f16},F6={f3,f18,f20,f4,f6,f16},F7={f1,f18,f20,f4,f6,f16}.
According to the feature scores, the feature sets’ priorities are derived:
(21)F1>F2>F3>F4>F5>F6>F7.


## 4.3. Ship Classification Based on Single Classifier

The experiment is carried out based on the test set of ship chips. When a ship chip is input for classification, the feature sets{Fn,n=1~7}derived above are set as the inputs for KNN classifier (the neighbor parameterK=3in this experiment). In this experiment, M=7, N=3. According to the feature sets’ priorities derived in the feature selection, the pair-wise comparison matrix can be set as
(22)A=[F1F2F3F4F5F6F7F11234567F212123456F3131212345F41413121234F515141312123F6161514131212F71716151413121].
The maximum eigenvalue ofAisλmax(A)=7.2786, CI=(λmax(A)-n)/(n-1)=0.0464. From Table 2, the value ofRIcould be acquired (RI=1.32), and thenCR=CI/RI=0.0352<0.1, so the consistency ofAis satisfied. The standardization form of the maximum eigenvector is
(23)s=[0.3086,0.2369,0.1745,0.1221,0.0803,0.0490,0.0286]T.Then, ship classification using KNN classifier is based on the seven feature sets{Fi,i=1,2,3,4,5,6,7}. As to each feature set, the comparison matrix for the scheme layer is constructed according to Section 3.3. And then, the comparison matrix for other feature sets could be constituted similarly. The corresponding weighting vectors can be calculated. The total scores for the three kinds (cargo ship, container ship, or oil tanker) can be derived and the final decision (the type of input ship chip) is made.When the classification decision of one ship chip is made, others are processed one by one based on the same method. Once the classification decisions of all the ship chips are made, the classification accuracy (for each type and the total accuracy) is calculated by the ratio of correct-classified number to the total number. Classification accuracies of the comparisons between the single-feature-set based classifications and the presented AHP-based classification are shown in Table3. The corresponding diagram of confusion matrix is presented in Figure 8. KNN classification based on the derived feature sets {Fn,n=1~7}could achieve satisfied performance. Inparticular for the AHP based method, the classification accuracy rates of the cargo and container ship come to 89.1% and 89.3%  respectively, and the comprehensive performance is superior to other methods with the total accuracy of 87.7%.Table 3
Ship classification accuracy based on KNN classifier.
Methods
Ship type
Total accuracy
Cargo ship
Oil tanker
Container ship
F
1
87.5%
77.3%
85.7%
85.1%
F
2
87.5%
72.7%
89.3%
85.1%
F
3
85.9%
77.3%
78.6%
82.5%
F
4
82.8%
68.2%
75%
78.1%
F
5
84.4%
72.7%
78.6%
80.7%
F
6
82.8%
68.2%
78.6%
78.9%
F
7
82.8%
63.6%
75%
77.2%
AHP
89.1%
81.8%
89.3%
87.7%Diagram of confusion matrix: (a)F1, (b) F2, (c) F3, (d) F4, (e) F5, (f) F6, (g) F7, and (h) AHP.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)In addition, some cases in which the vessel was wrongly classified (take the oil tanker case as an example) are presented in Table4 and some analysis of the possible causes is also presented. There are three chips whose oil pipe structure could be identified by human operator but not as clear as the ideal ones due to many interferential brighter scattering points around the oil pipe. However, they are all wrongly classified into cargo ships. In these three cases, the third chip is disturbed by strong scattering points badly and none of the chosen feature sets could be classified correctly, so the final decision based on AHP considers 0% possibility for the oil tanker. From the cases for chips No. 1 and No. 2, we can see that onlyF3andF5classify the first one into oil tanker, and F2, F3, F5, andF6classify the second one correctly. Referring to the detailed features for each set shown in Section 4.2, the possible reasons for these two cases are analyzed below.Table 4
Wrong classification cases for oil tanker.
Oil tanker cases
F
1
×
×
×
F
2
×
✓
×
F
3
✓
✓
×
F
4
×
×
×
F
5
✓
✓
×
F
6
×
✓
×
F
7
×
×
×
AHP
×
×
×
✓: identified; ×: unidentified.F
1classifies the first and second one wrongly; however,F3classifies both correctly andF2classifies the second one correctly. The features affecting their performances focus on the features f4 (area),f6 (shape complexity), and f16 (Hu moment). Considering their detailed performances,f4andf6are more sensitive to the classification of these two chips compared to other features. Actually, the interferential scattering points really affect ships’ shape complexity and the accurate extraction of targets’ area. Meanwhile, we can see that featuref16might be more preferable tof4orf6for the oil tanker classification although Section 4.2 shows that the latter two are a little more superior in the overall performance. Concluded from the performance of feature setF5(compared toF4), possibly Hu momentf20is superior to Hu momentf18in oil tanker classification.Actually ship features are due to the ship scattering, and the relations among them could be considered as being nonlinear, not totally independent. Essentially, the feature set including several features contributes to the classification as a whole. It is quite difficult to clarify their independent contribution clearly.When it refers to the second ship chip classification, four of seven feature sets are classified correctly, however,F2, F3, F5, and F6 have their own superiorities (in this experiment they are assigned 1/2, 1/3, 1/5, and 1/6). When the final decision based on AHP is made, the total score of classification into oil tanker is a little less than the score of classification into cargo. Actually, this is only a decision problem, and it just illustrates the possibility of falling within a certain type of ships. When the superiority factors of the feature sets vary proportionally, the results might vary but the conclusions are roughly consistent.


## 5. Discussion and Conclusion

This paper investigates the application of AHP in the feature selection and classification decision. Ship classification based on AHP could make more accurate decisions. It considers various contributions which is beneficial to the final decision and optimizes the process of decision. The experiment results illustrated that the AHP-based classification achieves higher accuracy.Similar to the KNN classifier, ship classification based on Bayes classifier, BP-NN classifier, SVM classifier, and the combination of multiple classifiers could also be implemented using AHP. The comparison matrix for theN classifiers is aN×Nmatrix. Because there exit no absolute assessments for the classifiers’ performance, the weighting coefficients need to be determined according to the different conditions. For example, when there are less samples comparatively, SVM has better performance, but the classification accuracy of BP-NN is difficult to be improved, so SVM has higher priority than BP-NN. When it is not easy to decide, the coefficients could be set as all ones (with coequal importance). The following steps are similar to ship classification based on single classifier. The only difference is the evaluation criteria. AHP based feature selection and classification are easy to be affected by the subjective factors, so the classification results may vary according to the different evaluations on the importance of criteria. How to establish objective evaluation criteria to constitute reasonable pair-wise comparison matrix needs to be emphasized in the future.Essentially, ship classification using space-borne SAR imagery is a quite difficult problem to deal with. This paper primarily researches on the decision problems for feature selection and classification. Ship classification is based on scattering characteristics affected by SAR monitoring, environmental conditions, and ship status, such as sea state, atmospheric phenomena, ship speed, and ship orientation. Emphasizing the illustration the AHP decision method concisely, this paper simplifies these experimental conditions properly. These problems are important for accurate classification and need to be considered in the future research. Ground truth could also be retrieved via transponder polls, that is, automatic identification system (AIS) [36], and provides references for the ship classification. Multisource data fusion (including SAR, optical satellite, infrared, coastal radar, and AIS, etc.) would provide more valid references and benefit ship classification [37]. These effective validation approaches are necessary to aid ship classification based on SAR imagery.



---

*Source: 698370-2013-12-10.xml*