Figure 1.
Principle of CBAM.
The need for multi-attribute decision-making brings more and more complexity, and this type of decision-making extends to an ever wider range of areas of life. A recent model that captures many components of decision-making frameworks is the complex q-rung picture fuzzy set (Cq-RPFS), a generalization of complex fuzzy sets and q-rung picture fuzzy sets. From a different standpoint, linguistic terms are very useful to evaluate qualitative information without specialized knowledge. Inspired by the ease of use of the linguistic evaluations by means of 2-tuple linguistic term sets, and the broad scope of applications of Cq-RPFSs, in this paper we introduce the novel structure called 2-tuple linguistic complex q-rung picture fuzzy sets (2TLCq-RPFSs). We argue that this model prevails to represent the two-dimensional information over the boundary of Cq-RPFSs, thanks to the additional features of 2-tuple linguistic terms. Subsequently, some 2TLCq-RPF aggregation operators are proposed. Fundamental cases include the 2TLCq-RPF weighted averaging/geometric operators. Other sophisticated aggregation operators that we propose are based on the Hamacher operator. In addition, we investigate some essential properties of the new operators. These tools are the building blocks of a multi-attribute decision making strategy for problems posed in the 2TLCq-RPFS setting. Furthermore, a numerical instance that selects an optimal machine is given to guarantee the applicability and effectiveness of the proposed approach. Finally, we conduct a comparison with other existing approaches.
Citation: Muhammad Akram, Ayesha Khan, Uzma Ahmad, José Carlos R. Alcantud, Mohammed M. Ali Al-Shamiri. A new group decision-making framework based on 2-tuple linguistic complex q-rung picture fuzzy sets[J]. Mathematical Biosciences and Engineering, 2022, 19(11): 11281-11323. doi: 10.3934/mbe.2022526
| [1] | Faisal Mehmood Butt, Lal Hussain, Anzar Mahmood, Kashif Javed Lone . Artificial Intelligence based accurately load forecasting system to forecast short and medium-term load demands. Mathematical Biosciences and Engineering, 2021, 18(1): 400-425. doi: 10.3934/mbe.2021022 |
| [2] | Xingyu Tang, Peijie Zheng, Yuewu Liu, Yuhua Yao, Guohua Huang . LangMoDHS: A deep learning language model for predicting DNase I hypersensitive sites in mouse genome. Mathematical Biosciences and Engineering, 2023, 20(1): 1037-1057. doi: 10.3934/mbe.2023048 |
| [3] | Wen Li, Xuekun Yang, Guowu Yuan, Dan Xu . ABCNet: A comprehensive highway visibility prediction model based on attention, Bi-LSTM and CNN. Mathematical Biosciences and Engineering, 2024, 21(3): 4397-4420. doi: 10.3934/mbe.2024194 |
| [4] | Xin Jing, Jungang Luo, Shangyao Zhang, Na Wei . Runoff forecasting model based on variational mode decomposition and artificial neural networks. Mathematical Biosciences and Engineering, 2022, 19(2): 1633-1648. doi: 10.3934/mbe.2022076 |
| [5] | Wenbo Yang, Wei Liu, Qun Gao . Prediction of dissolved oxygen concentration in aquaculture based on attention mechanism and combined neural network. Mathematical Biosciences and Engineering, 2023, 20(1): 998-1017. doi: 10.3934/mbe.2023046 |
| [6] | Guanghua Fu, Qingjuan Wei, Yongsheng Yang . Bearing fault diagnosis with parallel CNN and LSTM. Mathematical Biosciences and Engineering, 2024, 21(2): 2385-2406. doi: 10.3934/mbe.2024105 |
| [7] | Ying Chang, Lan Wang, Yunmin Zhao, Ming Liu, Jing Zhang . Research on two-class and four-class action recognition based on EEG signals. Mathematical Biosciences and Engineering, 2023, 20(6): 10376-10391. doi: 10.3934/mbe.2023455 |
| [8] | Xintong Wu, Yingyi Geng, Xinhong Wang, Jucheng Zhang, Ling Xia . Continuous extraction of coronary artery centerline from cardiac CTA images using a regression-based method. Mathematical Biosciences and Engineering, 2023, 20(3): 4988-5003. doi: 10.3934/mbe.2023231 |
| [9] | Vasileios E. Papageorgiou, Georgios Petmezas, Pantelis Dogoulis, Maxime Cordy, Nicos Maglaveras . Uncertainty CNNs: A path to enhanced medical image classification performance. Mathematical Biosciences and Engineering, 2025, 22(3): 528-553. doi: 10.3934/mbe.2025020 |
| [10] | Keruo Jiang, Zhen Huang, Xinyan Zhou, Chudong Tong, Minjie Zhu, Heshan Wang . Deep belief improved bidirectional LSTM for multivariate time series forecasting. Mathematical Biosciences and Engineering, 2023, 20(9): 16596-16627. doi: 10.3934/mbe.2023739 |
The need for multi-attribute decision-making brings more and more complexity, and this type of decision-making extends to an ever wider range of areas of life. A recent model that captures many components of decision-making frameworks is the complex q-rung picture fuzzy set (Cq-RPFS), a generalization of complex fuzzy sets and q-rung picture fuzzy sets. From a different standpoint, linguistic terms are very useful to evaluate qualitative information without specialized knowledge. Inspired by the ease of use of the linguistic evaluations by means of 2-tuple linguistic term sets, and the broad scope of applications of Cq-RPFSs, in this paper we introduce the novel structure called 2-tuple linguistic complex q-rung picture fuzzy sets (2TLCq-RPFSs). We argue that this model prevails to represent the two-dimensional information over the boundary of Cq-RPFSs, thanks to the additional features of 2-tuple linguistic terms. Subsequently, some 2TLCq-RPF aggregation operators are proposed. Fundamental cases include the 2TLCq-RPF weighted averaging/geometric operators. Other sophisticated aggregation operators that we propose are based on the Hamacher operator. In addition, we investigate some essential properties of the new operators. These tools are the building blocks of a multi-attribute decision making strategy for problems posed in the 2TLCq-RPFS setting. Furthermore, a numerical instance that selects an optimal machine is given to guarantee the applicability and effectiveness of the proposed approach. Finally, we conduct a comparison with other existing approaches.
A typhoon is a kind of tropical cyclone. When the sustained wind speed in the center of a tropical cyclone reaches category 12 to 13 (32.7–41.4 meters per second), it is called a typhoon [1]. Typhoons usually form on the sea surface in tropical areas 3–5 latitudes away from the equator, such as the North and South Pacific, the North Atlantic, and the Indian Ocean. Its movement is mainly affected by large-scale weather systems, and eventually dissipates at sea, becomes an extratropical cyclone, or dissipates after arriving on land. A typhoon is sudden and destructive; it will bring strong winds and rainstorms wherever it goes. The storm surge disaster caused by typhoons is one of the most serious natural disasters in the world. Sun et al. [2] found that typhoons will become stronger, larger, and more destructive in the context of global warming. Kossin [3] mentioned in his study that the translation speeds of tropical-cyclones have decreased globally by 10% and decreased by 30% in Asia during 1949–2016. The lower the velocity of a cyclone, the more rain falls on the location, and the higher the risk of flooding and other disasters. Therefore, slower translation speeds will lead to more serious disasters. In order to reduce the economic losses and casualties caused by typhoons, it is urgent to establish accurate and timely prediction methods.
The motion of a typhoon is highly complex. Unlike many time series problems, this cannot be solved by using the characteristics of periodicity [4]. A meteorological department mainly uses numerical simulation to achieve typhoon track prediction. This requires a lot of typhoon observational data. In order to establish complex atmospheric equations to simulate the motion and variation of the atmosphere, researchers need to fully understand the internal structural changes, key physical processes, and the mutation mechanism of the atmosphere [5]. South Korea, for example, conducts numerical simulations on a Cray XC40 supercomputer with 139, 329 CPUs. Due to its difficulty, numerical simulation requires a lot of calculation time. Even so, numerical simulation cannot obtain the analytical solution, only the numerical solution. Also, expensive hardware means high maintenance costs.
A deep learning method can use a single graphics processing unit (GPU). The prediction results can be generated within seconds by updating the input data of the prediction model. Compared with using hundreds of thousands of CPUs in a simulation, this method is obviously quicker and cheaper. There have been many investigations into the use of deep learning methods for typhoon track prediction. Song et al. [6] used the spatial location information and meteorological factors of typhoons to train a neural network which is constructed based on convolutional neural networks (CNNs) and the gated recurrent neural network (GRU). They state that the proposed model is better than the traditional prediction methods for typhoon track prediction. Song et al. [7] used BiGRU with an attention mechanism to predict typhoon tracks. Gao et al. [8] trained a LSTM with typhoon track data. In addition to describing the relevant physical quantities of the typhoon process, satellite images can also effectively reflect the spatial information.
Satellite images can provide the distribution of clouds, and the changes between satellite images can reflect the evolution of weather systems. With the help of satellite images, more accurate weather forecasts can be made. Lee et al. [9] were probably the first to consider using satellite images of tropical cyclones for typhoon track prediction with the help of neural networks. However, the images did not function as input data for the LSTM. Therefore, they extract information such as Dvorak number, maximum wind speed, or spatial location of the cyclone from satellite images and input these into the network for time series prediction. Kovordanyi et al. [10] used satellite images as neural network input data for the first time. They predicted the shape of cloud clusters and the direction of movement of the typhoon via a neural network. Rüttgers et al. [11] marked the typhoon center on the satellite images and used generative adversarial networks (GANs) [12] to generate the satellite images to predict the shape of the cloud cluster and the typhoon track.
The prediction of typhoon tracks is affected by many factors such as the typhoon model structure, terrain, initial field, and boundary conditions of the model [13]. For example, the movement of a typhoon is affected by the steering flow, which is determined by the environment around the storm. When a typhoon is affected by multiple weather systems at the same time, it is difficult to predict the steering flow [14]. Therefore, a single neural network cannot accurately predict the typhoon track. Therefore, the fusion of multiple neural networks is more suitable for the task of typhoon track prediction rather than a single neural network [7]. Xu et al. [15] modeled the spatiotemporal relationship. They attached the trained generator to the output of the LSTM to generate the cloud images. Hong et al. [16] applied a convolutional LSTM network (ConvLSTM) to typhoon track prediction.
At present, research based on deep learning methods combined with satellite images mostly focuses on typhoon eye tracking [17] or typhoon intensity prediction [18], and less on typhoon track prediction. We believe that it is necessary to extract the temporal and spatial features of time series satellite images for typhoon path prediction. Therefore, we propose a new typhoon track prediction model based on a CNN and LSTM: DeepTyphoon. The specific contributions of this paper are as follows: We mark the satellite images obtained from the Japan Meteorological Agency to create the data set. According to the characteristics of the images in the data set, we establish a prediction network. The features of the clouds and marked squares in the image are extracted by dilated convolution. In addition, a convolutional block attention module (CBAM) is used to enhance the feature extraction ability of the network from two dimensions: Channel and spatial.
The attention module can focus on the local information of the network and suppress invalid information, so as to enhance the expression of effective features. The attention mechanism was initially applied to natural language processing and achieved good results. However, there are few studies on the application of attention mechanisms to time series prediction in meteorology. The convolutional block attention module is a new convolutional attention module proposed by Sanghyun et al. in 2018 [19]. The convolution operation extracts information features by mixing cross channel information and spatial information. Therefore, they innovatively proposed the attention mechanism of integrating channel attention and spatial attention and applied a channel attention module and a spatial attention module for the two respective dimensions. The principle of CBAM is shown in Figure 1.
In this paper, CBAM is applied to the encoder of the neural network to enhance the feature extraction ability of the convolution network, and extracts the features of marked squares and clouds from the two dimensions.
The long short-term memory network is a special recurrent neural network. It was first proposed by Hochreiter and Schmidhuber [20] and later improved by Grave [21], which can address the problems of gradient disappearance and gradient explosion in the process of long sequence training. Like RNNs, LSTMs can also connect previous information to the current task, such as using the previous video frame information to understand the current frame. Compared with RNNs, LSTMs perform better in learning long-term dependencies.
The operation of LSTM is mainly realized by the input gate it, forget gate ft, and output gate ot. First, the LSTM uses the forget gate ft to decide which information in ct-1 to delete, and the activation state of the forget gate is determined by the activation function σ(·):
| ft=σ(Wf⋅[ht−1,xt]+bf) | (1) |
where Wf is the weight matrix of the forget gate; bf is the offset term of the forget gate; ft represents the forget gate at time t; xt is the input of network at time t. Second, the LSTM uses the input gate it to determine the information stored in the new memory unit ct:
| it=σ(Wi⋅[ht−1,xt]+bi) | (2) |
| ~ct=tanh(Wc⋅[ht−1,xt]+bc) | (3) |
where Wi is the weight matrix of the input gate; Wc is the weight moment of the cell state; bi is the offset term of the input gate; bc is the bias term of the cell state; it is the input gate at time t; tanh is the hyperbolic tangent function; ~ct indicates the candidate value ct added to the new memory unit status. Then, the LSTM updates the status of the previous memory unit ct-1 to the status ct:
| ct=ft×ct−1+it×~ct | (4) |
where ct is the input status of the memory unit. Finally, the output gate ot is used to calculate the hidden state ht. The specific formula is as follows:
| ot=σ(Wo⋅[ht−1,xt]+bo) | (5) |
| ht=ot⊙tanh(ct) | (6) |
where Wo is the weight matrix of the output gate; bo is the output gate offset term; σ is the sigmoid function; ⊙ represents the multiplication of matrix elements.
The dilated residual network [22] was proposed to improve the receptive field of convolutions. It maintains the width and height of input features by setting dilation rates. However, the setting of dilation rates leads to the "gridding issue"; when the dilation rate increases, the input sampling will become very sparse, which will be detrimental to the learning of features. Therefore, Wang et al. [23] proposed hybrid dilated convolution (HDC) which alleviates the gridding issue by setting the dilation rates for a series of convolution layers and connecting them. In this paper, the size of the marked square on the satellite images is 3 × 3 pixels. As shown in Figure 2, the three-layer HDC layer with a dilation rate of 1/2/3 can ensure that the input feature of the square is not lost.
The self-built typhoon image data set used in this paper is provided by the National Institute of Information (NII). The satellite image data are taken from Himawari-1~8 series and GOES meteorological satellites. Of these, Himawari-8 can achieve a 10-minute time resolution and 500 m spatial resolution. It has an improved imager observing five bands (one for visible (0.64 μm) and four for infrared (10.4, 12.4, 6.2 and 3.9 μm)). GOES is the geostationary orbit operational satellite series of NOAA in the United States. It also has an improved imager observing five bands (one for visible (0.55–0.75 μm) and four for infrared (10.20–11.20, 11.50–12.50, 6.50–7.00 and 3.80–400 μm)). Satellite images in the recent 40 years are selected as data samples, which are all infrared images (10.4 μm) with a size of 128 × 128 pixels, representing an actual range of 2600 × 2600 km. The middle of the image is the typhoon center, and NII provides the longitude and latitude of the typhoon position corresponding to each image.
Due to the different life cycles of typhoons, the typhoon processes with typhoon cycle less than 60 hours are screened out through data cleaning in order to ensure that the neural network can obtain sufficient information in the training process. We selected satellite images located in the Northwest Pacific and Southwest Pacific to demonstrate the effectiveness of the model. Of these, 513 typhoon processes located in the Northwest Pacific from 1980 to 2020, that is, 513 typhoon sequences, and a total of 13400 satellite images are selected to construct the training set, verification set, and test set, which are 338, 103 and 72 typhoon sequences, respectively, corresponding to 9110, 2680 and 1610 satellite images. Also, 198 typhoon processes located in the Southwest Pacific from 2000 to 2020, a total of 8420 satellite images, are selected to construct the training set, verification set, and test set, which are 108, 60 and 30 typhoon sequences, respectively, corresponding to 5470, 2110 and 840 satellite images.
The original infrared image is a single channel gray image, and the range of digital number of satellite image is 0–255. Typhoons are sequential from start to finish. Moreover, the same typhoon has different forms in different lifecycle stages, and the spiral cloud system shown on the satellite cloud map is also different [24]. Figure 3 shows some images in the lifecycle of typhoon Damrey, which lasted for 11 days and 18 hours.
As described in Section 2.2.1, the typhoon center in the satellite image provided by NII is located in the middle of the image. However, the typhoon center on the satellite images is an absolute position. We need the relative positional information between satellite image sequences to judge the movement of clouds. In order to determine the direction of movement and distance of the typhoon through the satellite image, for each satellite image, we mark the typhoon position at the previous time (for example, six hours ago) on the satellite image at the current time. The marking method is as follows: The spatial coordinate of the typhoon at the current time is (φ, λ). The horizontal offset Dφ and vertical offset Dλ (in km) between the coordinates of the previous time and the coordinates of the current time is calculated according to the Havers formula [25]; the earth radius R is taken from the real coordinate system:
| Dφ=2Rarcsin√sin2(φ′−φ2) | (7) |
| Dλ=2Rarcsin√cosφcosφ′sin2(λ′−λ2) | (8) |
Then, according to Dφ and Dλ, we calculate the horizontal and vertical offset Pφ and Pλ at the pixel level on the image using the reference system. This determines the position of the typhoon at the previous time on the satellite image at the current time.
We use black squares to mark satellite images (gray image) for the experiments. Compared with using RGB image (image data marked with red squares) as the input of neural network, the use of gray images as the input of neural network means a shorter training time which is convenient for the experiments to adjust the parameters of the neural network. In the experimentation, it is found that using a red square marked image can eliminate the interference of the gray background, so as to determine the position of the typhoon more accurately. As shown in Figure 4: Figure 4(a) is the original satellite image, Figure 4(b) is the satellite image marked with a black square, and Figure 4(c) is the satellite image marked with a red square. In Figure 4(b), (c), the center of the image is the typhoon center at the current time, and the marked square represents the typhoon center at the previous time. The size of the marked square is 3 × 3 pixels. There are 13400 single channel gray images in the typhoon image data set with black square marks, and 13, 400 RGB three channel images in the typhoon image data set with red square marks. The range of digital number of image data marked with square (black or red) are both 0–255.The size of the dataset does not change as a result of marking. The data set size of the Northwest Pacific is 9110 images for the training set, 2680 images for the verification set, and 1610 images for the test set. For the Southwest Pacific data set, this is 5470, 2110 and 840, respectively.
In the research of using deep learning to predict typhoon tracks, most researchers use reanalysis data as the input of neural network. For example, the work in [26] converts the reanalysis data into a density map and uses it as the input of a neural network to determine the spatial position of the typhoon according to the predicted density map. However, compared with the thermodynamic diagram transformed according to physical quantity data, satellite images contain more intuitive information. In order to consider the temporal and spatial characteristics of the typhoon process at the same time, we propose the DeepTyphoon prediction model for satellite image training and prediction. DeepTyphoon combines the ability of a CNN to extract advanced semantic features and the utility of LSTM for time series prediction. It takes the marked satellite image as the input, learns the relative position between the square mark and the cloud cluster and the direction of movement of the cloud cluster in the training process, and finally predicts the position of the typhoon according to the square mark on the predicted image.
The prediction process of DeepTyphoon is shown in Figure 5. The whole neural network model includes a convolution neural network (as encoder), an LSTM, and a deconvolution neural network (as decoder). Taking the satellite images marked with black squares as an example, the prediction process is as follows: First, the marked satellite images (denoted by m) are input into the encoder. We use the hybrid dilated convolution layer to extract the features of the satellite image m. At the same time, we add the channel attention module and the spatial attention module to further improve the feature extraction ability of the coding network. After each image of a satellite image sequence is input into the encoder, the encoder outputs the feature vector of m (denoted by t) as the input of LSTM. The output of the LSTM network t' is the predicted value of t. Finally, we deconvolute t' into an image m' in order to determine the typhoon center according to the prediction result t' of LSTM. Due to the black square mark in the input image m, the final output of the model, the predicted image m', will also have a predicted black square mark.
The black square in the predicted image is detected by applying a filter on the predicted image and a convolution with the size of the black square on the filtered image. The pixel with the highest value in the convolved image gives the location of the predicted typhoon center. After identifying the pixel coordinates of the predicted black square, they are converted back to coordinates of latitude (φ) and longitude (λ) through the geographic reference system.
The experimental environment of this paper is mainly based on the PyTorch framework and built on the Ubuntu operating system. The specific configuration is shown in Table 1 below. In order to fully verify the performance, the experiments in this paper are completed in the same experimental environment. The experimental hyperparameters are set as follows: The batch size is 100, the learning rate is 0.0001, and the number of epochs is 100. The Adam algorithm is used to optimize the network parameters, and the weight attenuation coefficient is set to 0.00002.
| Name | Configuration |
| Operating System | Ubuntu 18.04.5 |
| Framework | PyTorch 1.8.0 |
| Memory | 60 Gb |
| GPU | GeForce RTX 3080Ti |
| CUDA | CUDA 10.1 |
DownLoad:
CSV
In order to test the influence of the number of convolution layers on the experimental results, we set different convolution layers for testing. First, we set the dilation rates of the three continuous convolution layers to 1/2/3, respectively, and take the three convolution layers as one HDC layer. We use the satellite image marked with the red square as the input of the model, and calculate the mean average error (MAE) of the predicted results when one, two, and three HDC layers are used.
As shown in Table 2, we achieved the minimum MAE when using two HDC layers. That is to say, when we use six convolution layers and set the dilation rates to 1/2/3/1/2/3 respectively, the prediction result of the neural network is the best.
| Number of HDC layers | MAE (km) |
| 1 | 78.22 |
| 2 | 66.38 |
| 3 | 69.43 |
DownLoad:
CSV
An important hyperparameter of the model is the number of input image sequences. If the number of input image sequences (denoted by n) is too large, the memory usage will be increased, and the training time will be increased. Moreover, the lifecycle of a typhoon is not too long, and too large a value of n will make most typhoon processes unable to be used as training data. However, too small a value of n will lead to insufficient input information and affect the prediction performance. We use the satellite images marked with red squares as the input. The time interval of the satellite image is 6 hours. After training and testing, for the three cases where n is 5, 10 and 20, the MAE after 6 hours is calculated.
As shown in Table 3: In the case of n = 5, the change of the typhoon seems to have not been fully captured. In the case of n = 20, due to the reduction of the number of training data sets, the results are not good compared with n = 10. For example, considering typhoon Mitag with 29 sequences, using n = 10 allows to start training from sequence 11, and provides a number of 19 training sequences. However, when n = 20, this only allows to start training from sequence 21, that is, 9 training sequences.
| n | MAE (km) |
| 5 | 96.05 |
| 10 | 68.49 |
| 20 | 79.43 |
DownLoad:
CSV
During the experiment, it is found that different convolution kernel sizes will affect the shape of the final generated square. In order to accurately predict the square in the predicted images, the convolution kernel size is tested. We select the satellite image marked with a black square as the input, which can save calculation time.
Figure 6 shows the prediction results corresponding to different convolution kernels. Figure 6(a) shows the prediction results when the convolution kernel size is 3. It is obvious from the prediction image that a 3 × 3 pixel black square is generated. Figure 6(b), (c) correspond to the shape of the predicted square in the generated image when the convolution kernel size is 2 and 4, respectively. Although it can be seen from the predicted image that the marks are generated, they are blurred to varying degrees and the central position cannot be accurately determined by convolution. Finally, the size of the convolution kernel is set to 3 so that the exact position of the square can be obtained.
For the typhoon satellite image at the later stage of the typhoon lifecycle, as shown in Figure 7(a), if the black square is used to mark the satellite image, the square in the predicted image almost coincides with the gray background of the image. It is impossible to accurately detect the mark position in the prediction image through convolution. Therefore, we choose the satellite image marked by the red square as the model input so that we can accurately detect the location of the marked pixels. The channel attention module is added after the convolution layer to enhance the weight of the pixel value of the image R channel (red). In the predicted image, only the pixel value on the R channel in the generated image is detected, and the position of the marked square can be obtained accurately. As shown in Figure 7(b), the red square in the figure is located and redrawn after convolution on the R channel.
The parameters of the neural network are finally determined through experiments, as shown in Table 4. Of these, the number of HDC layers is obtained from Table 2, which is six convolution layers with different dilation rates. In order to decode the prediction vector into an image, the number of deconvolution layers is the same as that of convolution layers. The kernel size is determined according to Figure 6. Based on previous experience, we set the number of hidden layers of LSTM to two.
| Name | Parameter |
| HDC layers | 2 |
| Deconvolution layers | 6 |
| LSTM layers | 2 |
| Kernel | 3 × 3 |
DownLoad:
CSV
In order to evaluate the performance of the prediction model, the root mean square error (RMSE) and the mean absolute error are selected to measure the prediction accuracy of the model.
The absolute error (E) describes the distance between the predicted coordinates (φpred, λpred) and the real coordinates (φreal, λreal). It is calculated in kilometers (km) by applying the haversine formula, and the earth radius R is taken from the real coordinate position:
| E=2Rarcsin√sin2(φpred−φreal2)+cosφpredcosφpredsin2(λpred−λreal2) | (9) |
The mean absolute error formula is as follows, in km, where E represents the absolute error of the typhoon:
| MAE=1m∑i=1m|E| | (10) |
The root mean square error formula is as follows, in km, where E represents the absolute error of the typhoon:
| RMSE=√1mm∑i=1(Ei)2 | (11) |
As an example of a predicted image, Figure 8 shows the prediction results of typhoon Gray sequence 11, in which Figure 8(a) is the ground truth and Figure 8(b) is a predicted image.
It is obvious that the predicted image does not have the same definition as the real image; in particular, the clouds are blurred compared with the ground truth. However, we focus on accurately predicting the position of the square in the predicted images, so as to accurately predict the spatial position of the typhoon rather than improve the quality of the generated images.
By calculating the relative position between the square and the center of the predicted image, the spatial position of the predicted typhoon is determined. We compare the position of the predicted typhoon with the position of the real typhoon, and finally we can obtain the prediction accuracy. We select an example in the Northwest Pacific and the Southwest Pacific to show the prediction results of DeepTyphoon. Each typhoon is predicted in chronological order, and the sequence is named by the date and time of the corresponding real image UTC. For example, "1993072718" represents 6 p.m. on July 27, 1993 (UTC).
Table 5 shows the coordinate prediction results and prediction errors of Typhoon Mitag. As can be seen from the table, this typhoon has taken 29 steps from occurrence to extinction; each time step is 6 hours, giving a total of 174 hours. Take ten points as a sequence as the input of the prediction model, the predicted values of Mitag sequence 11 to sequence 29 can be obtained in order. After calculation, the average error of Typhoon Mitag from sequence 11 to sequence 29 is 55.037 km.
| Sequence | Real | Predicted | E (km) |
| 11(2007112300) | 13.7, 127.2 | 13.6, 126.8 | 44.629 |
| 12(2007112306) | 13.7, 127.0 | 14, 127.1 | 35.062 |
| 13(2007112312) | 14.0, 126.7 | 14.5, 126.6 | 56.632 |
| 14(2007112318) | 14.1, 126.5 | 14.3, 126.1 | 48.516 |
| 15(2007112400) | 14.3, 126.2 | 14.8, 126.4 | 59.619 |
| 16(2007112406) | 14.5, 125.9 | 15, 126.2 | 64.278 |
| 17(2007112412) | 14.8, 125.5 | 15, 124.9 | 68.201 |
| 18(2007112418) | 15.1, 125.2 | 15.4, 125.5 | 46.353 |
| 19(2007112500) | 15.5, 124.5 | 16.1, 124.7 | 70.065 |
| 20(2007112506) | 16.4, 123.9 | 16.1, 124.1 | 39.606 |
| 21(2007112512) | 16.9, 122.9 | 16.8, 122.5 | 43.997 |
| 22(2007112518) | 17.3, 121.5 | 17, 122.1 | 71.951 |
| 23(2007112600) | 18.3, 121.0 | 18.7, 121.4 | 61.298 |
| 24(2007112606) | 19.0, 120.7 | 19.4, 120.5 | 49.187 |
| 25(2007112612) | 19.5, 120.7 | 19.2, 121.1 | 53.609 |
| 26(2007112618) | 20.0, 121.1 | 20.6, 121.4 | 73.688 |
| 27(2007112700) | 20.1, 121.5 | 20.4, 122.6 | 56.719 |
| 28(2007112706) | 20.5, 122.8 | 20.6, 123.6 | 23.619 |
| 29(2007112712) | 20.6, 124.0 | 20, 124.4 | 78.685 |
| MAE (km) |  |
 |
55.037 |
DownLoad:
CSV
As shown in Figure 9, in order to display the forecast results more intuitively, they are marked in the map according to the real track (yellow) and predicted track (blue) of Typhoon Mitag.
Table 6 shows the coordinate prediction results and prediction errors of Typhoon Marcia. Typhoon Marcia has taken 22 steps from occurrence to extinction; each time step is 6 hours, giving a total of 132 hours. Take ten points as a sequence as the input of the prediction model, the predicted values of Marcia sequence 11 to sequence 22 can be obtained in order.
| Sequence | Real | Predicted | E (km) |
| 11(2015021909) | –20.5, 150.6 | –20.4, 150.3 | 33.175 |
| 12(2015021912) | –20.8, 150.5 | –20.6, 150.9 | 47.177 |
| 13(2015021915) | –21.1, 150.5 | –21.5, 150.8 | 54.261 |
| 14(2015021918) | –21.5, 150.5 | –21.5, 151.1 | 70.413 |
| 15(2015021921) | –22.1, 150.5 | –21.9, 150.9 | 46.853 |
| 16(2015022000) | –22.7, 150.5 | –22.5, 150.6 | 24.494 |
| 17(2015022003) | –23.1, 150.5 | –22.7, 150.2 | 54.061 |
| 18(2015022006) | –23.7, 150.6 | –23.3, 150.9 | 53.983 |
| 19(2015022009) | –24.2, 150.9 | –24.7, 150.2 | 90.065 |
| 20(2015022012) | –24.5, 151.3 | –24.9, 150.8 | 67.302 |
| 21(2015022015) | –24.9, 151.6 | –25.1, 151 | 64.426 |
| 22(2015022018) | –25.3, 151.8 | –25.4, 151.4 | 41.705 |
| MAE (km) |  |
 |
53.992 |
DownLoad:
CSV
As shown in Figure 10, we also display the forecast results for Typhoon Marcia in the map according to Table 6. The real track is represented by yellow points, and the predicted path is represented by blue points.
We tested two test sets located in the Northwest Pacific and Southwest Pacific. Table 7 shows the prediction results of DeepTyphoon for the two test sets. The predicted time range is 6 hours. The average absolute error of the prediction results of DeepTyphoon for the two test sets is less than 70 km. In addition, we recorded the calculation time of the DeepTyphoon prediction model. After 14 hours of training, the DeepTyphoon prediction model can process 1000 satellite images, that is, 100 typhoon points, in one minute.
| Data set | Test set size | RMSE (km) | MAE (km) | Computational time (s) |
| Northwest Pacific | 1610 | 72.72 | 65.39 | 87.42 |
| Southwest Pacific | 840 | 69.58 | 62.41 | 49.80 |
DownLoad:
CSV
Table 8 shows the prediction results for RNN, GAN, AE–GRU, TraijGRU, and DeepTyphoon. The performance indicators include RMSE and MAE.
| Model | Predicted Time Range | RMSE (km) | MAE (km) |
| RNN[27] | 12 h | 164.22 | |
| AE-GRU[28] | 12 h | 138.67 | |
| TraijGRU[29] | 6 h | 66.6 | |
| DeepTyphoon | 6 h | 73.96 | 64.17 |
| 12 h | 134.70 | 129.82 |
DownLoad:
CSV
As shown in Table 8, the evaluation indexes for DeepTyphoon are the best after 6 and 12 hours. Through comparison, it is shown that DeepTyphoon can more accurately predict the location of typhoons using satellite images.
In order to further verify the validity of the proposed model, we calculated the MAEs in longitude and latitude and compared it with SOTA models. As shown in Table 9, our model has the best prediction result at 6 hours and 12 hours. This verifies the validity of the proposed model.
| Model | Latitude | Longitude | ||
| 6 h | 1 h | 6 h | 12 h | |
| DeepFR[30] | 0.9100 | 1.1765 | 0.8910 | 1.2430 |
| Trj-DMFMG[31] | 0.6602 | 0.9275 | 0.8747 | 1.1329 |
| DeepTyphoon | 0.4876 | 0.8812 | 0.6645 | 1.0874 |
DownLoad:
CSV
The experimental results show that marking the typhoon center at the previous time on the satellite images and using CNN and LSTM for feature extraction can effectively learn cloud motions in the satellite image. Also, we can use the deconvolution network to detect the typhoon center at a future time so as to accurately predict the typhoon path in the future.
However, we expect to get more information from the prediction image obtained by the deconvolution network. Therefore, the focus of future work should be to further optimize the neural network and better extract cloud features, so as to obtain more realistic satellite images through prediction. In this way, we can not only detect the typhoon center in the prediction image, but also determine the intensity of the typhoon. In addition, the prediction accuracy of the typhoon path should be further improved. Then, how to ensure the accuracy of typhoon track prediction and the authenticity of the prediction image at the same time is the difficulty in the future.
In this paper, we use a deep learning method combined with satellite images to predict typhoon tracks, and propose the DeepTyphoon prediction model. The typhoon center is marked on the satellite image. Combined with the characteristics of CNN and LSTM, we extract the features of the satellite image by HDC and CBAM, and use LSTM to predict the features. After the predicted feature vector is transformed into an image, the position of the typhoon is detected on the predicted image. More than 20000 time series satellite images in two data sets released by the Japan Meteorological Agency in recent years are used for training and prediction. Experimental results on the Northwest Pacific and Southwest Pacific data sets show that our prediction method is better than existing deep learning methods in typhoon track prediction.
This work was supported by "Shanghai science and technology innovation plan" (Project number: 20dz1203800) and "the Capacity Development for Local College Project" (Project number: 19050502100).
The authors declare there is no conflict of interests.
| [1] |
L. A. Zadeh, Fuzzy sets, Inf. Control, 8 (1965), 338–353. https://doi.org/10.1016/S0019-9958(65)90241-X doi: 10.1016/S0019-9958(65)90241-X
|
| [2] |
S. M. Chen, W. H. Hsiao, W. T. Jong, Bidirectional approximate reasoning based on interval-valued fuzzy sets, Fuzzy Sets Syst., 91 (1997), 339–353. https://doi.org/10.1016/S0165-0114(97)86594-3 doi: 10.1016/S0165-0114(97)86594-3
|
| [3] |
S. M. Chen, W. H. Hsiao, Bidirectional approximate reasoning for rule-based systems using interval-valued fuzzy sets, Fuzzy Sets Syst., 113 (2000), 185–203. https://doi.org/10.1016/S0165-0114(98)00351-0 doi: 10.1016/S0165-0114(98)00351-0
|
| [4] |
S. M. Chen, W. T. Jong, Fuzzy query translation for relational database systems, IEEE Trans. Syst. Man Cybern. Syst. Part B, 27 (1997), 714–721. https://doi.org/10.1109/3477.604117 doi: 10.1109/3477.604117
|
| [5] |
S. M. Chen, S. J. Niou, Fuzzy multiple-attributes group decision-making based on fuzzy preference relations, Expert Syst. Appl., 38 (2011), 3865–3872. https://doi.org/10.1016/j.eswa.2010.09.047 doi: 10.1016/j.eswa.2010.09.047
|
| [6] |
M. I. Ali, J. Zhan, M. J. Khan, T. Mahmood, H. Faizan, Another view on knowledge measures in atanassov intuitionistic fuzzy sets, Soft Comput., 26 (2022), 6507–6517. https://doi.org/10.1007/s00500-022-07127-3 doi: 10.1007/s00500-022-07127-3
|
| [7] |
K. T. Atanassov, Intuitionistic fuzzy sets, Fuzzy Sets Syst., 20 (1986), 87–96. https://doi.org/10.1007/978-3-7908-1870-3-1 doi: 10.1007/978-3-7908-1870-3-1
|
| [8] |
W. Wang, J. Zhan, J. Mi, A three-way decision approach with probabilistic dominance relations under intuitionistic fuzzy information, Inf. Sci., 582 (2022), 114–145. https://doi.org/10.1016/j.ins.2021.09.018 doi: 10.1016/j.ins.2021.09.018
|
| [9] | R. R. Yager, Pythagorean fuzzy subsets, in 2013 Joint IFSA World Congress and NAFIPS Annual Meeting (IFSA/NAFIPS), (2013), 57–61. https://doi.org/10.1109/IFSA-NAFIPS.2013.6608375 |
| [10] |
R. R. Yager, Pythagorean membership grades in multi-criteria decision making, IEEE Trans. Fuzzy Syst., 22 (2013), 958–965. https://doi.org/10.1109/TFUZZ.2013.2278989 doi: 10.1109/TFUZZ.2013.2278989
|
| [11] |
T. Senapati, R. R. Yager, Fermatean fuzzy sets, J. Ambient. Intell. Humaniz. Comput., 11 (2020), 663–674. https://doi.org/10.1007/s12652-019-01377-0 doi: 10.1007/s12652-019-01377-0
|
| [12] |
V. Torra, Hesitant fuzzy sets, Int. J. Intell. Syst., 25 (2010), 529–539. https://doi.org/10.1002/int.20418 doi: 10.1002/int.20418
|
| [13] |
W. Wang, X. Ma, Z. Xu, W. Pedrycz, J. Zhan, A three-way decision method with prospect theory to multi-attribute decision-making and its applications under hesitant fuzzy environments, Appl. Soft Comput., 126 (2022), 109283. https://doi.org/10.1016/j.asoc.2022.109283 doi: 10.1016/j.asoc.2022.109283
|
| [14] |
R. R. Yager, Generalized orthopair fuzzy sets, IEEE Trans. Fuzzy Syst., 26 (2016), 1222–1230. https://doi.org/10.1109/TFUZZ.2016.2604005 doi: 10.1109/TFUZZ.2016.2604005
|
| [15] |
B. C. Cuong, V. Kreinovich, Picture fuzzy sets, J. Comput. Sci. Cybern., 30 (2014), 409–420. https://doi.org/10.1109/WICT.2013.7113099 doi: 10.1109/WICT.2013.7113099
|
| [16] |
L. Li, R. T. Zhang, J. Wang, X. P. Shang, K. Y. Bai, A novel approach to multi-attribute group decision-making with q-rung picture linguistic information, Symmetry, 10 (2018), 172. https://doi.org/10.3390/sym10050172 doi: 10.3390/sym10050172
|
| [17] |
M. Akram, S. Alsulami, F. Karaaslan, A. Khan, q-Rung orthopair fuzzy graphs under Hamacher operators, J. Intell. Fuzzy Syst., 40 (2021), 1367–1390. https://doi.org/10.3233/JIFS-201700 doi: 10.3233/JIFS-201700
|
| [18] |
M. J. Khan, J. C. R. Alcantud, P. Kumam, W. Kumam, A. N. Al-Kenani, An axiomatically supported divergence measures for q-rung orthopair fuzzy sets, Int. J. Intell. Syst., 36 (2021), 6133–6155. https://doi.org/10.1002/int.22545 doi: 10.1002/int.22545
|
| [19] |
L. A. Zadeh, The concept of a linguistic variable and its application to approximate reasoning Part Ⅰ, Inf. Sci., 8 (1975), 199–249. https://doi.org/10.1016/0020-0255(75)90046-8 doi: 10.1016/0020-0255(75)90046-8
|
| [20] |
H. Zhang, Linguistic intuitionistic fuzzy sets and application in MAGDM, J. Appl. Math., 2014 (2014), 1–11. https://doi.org/10.1155/2014/432092 doi: 10.1155/2014/432092
|
| [21] |
H. Garg, Linguistic Pythagorean fuzzy sets and its applications in multiattribute decision-making process, Int. J. Intell. Syst., 33 (2018), 1234–1263. https://doi.org/10.1002/int.21979 doi: 10.1002/int.21979
|
| [22] |
M. Lin, J. Wei, Z. Xu, R. Chen, Multiattribute group decision-making based on linguistic Pythagorean fuzzy interaction partitioned Bonferroni mean aggregation operators, Complexity, 2018 (2018), 1–24. https://doi.org/10.1155/2018/9531064 doi: 10.1155/2018/9531064
|
| [23] |
M. Lin, X. Li, L. Chen, Linguistic q-rung orthopair fuzzy sets and their interactional partitioned Heronian mean aggregation operators, Int. J. Intell. Syst., 35 (2020), 217–249. https://doi.org/10.1155/2018/9531064 doi: 10.1155/2018/9531064
|
| [24] |
M. Akram, S. Naz, S. A. Edalatpanah, R. Mehreen, Group decision-making framework under linguistic q-rung orthopair fuzzy Einstein models, Soft. Comput., 25 (2021), 10309–10334. https://doi.org/10.1007/s00500-021-05771-9 doi: 10.1007/s00500-021-05771-9
|
| [25] |
F. Herrera, L. Martínez, A 2-tuple fuzzy linguistic representation model for computing with words, IEEE Trans. Fuzzy Syst., 8 (2000), 746–752. https://doi.org/10.1109/91.890332 doi: 10.1109/91.890332
|
| [26] |
F. Herrera, L. Martínez, An approach for combining linguistic and numerical information based on the 2-tuple fuzzy linguistic representation model in decision-making, Int. J. Uncertain. Fuzz. Knowl. Based Syst., 8 (2000), 539–562. https://doi.org/10.1142/S0218488500000381 doi: 10.1142/S0218488500000381
|
| [27] |
X. Liu, H. S. Kim, F. Feng, J. C. R. Alcantud, Centroid transformations of intuitionistic fuzzy values based on aggregation operators, Mathematics, 6 (2018), 215. https://doi.org/10.3390/math6110215 doi: 10.3390/math6110215
|
| [28] |
Y. Liu, Y. Qin, F. Liu, Y. Rong, GIBWM-MABAC approach for MAGDM under multi-granularity intuitionistic 2-tuple linguistic information model, J. Ambient Intell. Hum. Comput., 2021 (2021), 1–17. https://doi.org/10.1007/s12652-021-03476-3 doi: 10.1007/s12652-021-03476-3
|
| [29] |
A. Luqman, M. Akram, J. C. R. Alcantud, Digraph and matrix approach for risk evaluations under Pythagorean fuzzy information, Expert Syst. Appl., 170 (2021), 114518. https://doi.org/10.1016/j.eswa.2020.114518 doi: 10.1016/j.eswa.2020.114518
|
| [30] |
Y. Qin, Y. Liu, S. Abdullah, G. Wei, Group decision support methodology based upon the multigranular generalized orthopair 2-tuple linguistic information model, Int. J. Intell. Syst., 36 (2021), 3367–3400. https://doi.org/10.1002/int.22419 doi: 10.1002/int.22419
|
| [31] |
Y. Xu, H. Wang, Approaches based on 2-tuple linguistic power aggregation operators for multiple attribute group decision making under linguistic environment, Appl. Soft Comput., 11 (2011), 3988–3997. https://doi.org/10.1016/j.asoc.2011.02.027 doi: 10.1016/j.asoc.2011.02.027
|
| [32] |
M. A. Dulebenets, J. Pasha, M. Kavoosi, O. F. Abioye, E. E. Ozguven, R. Moses, et al., Multiobjective optimization model for emergency evacuation planning in geographical locations with vulnerable population groups, J. Manag. Eng., 36 (2020), 04019043. https://doi.org/10.1061/(ASCE)ME.1943-5479.0000730 doi: 10.1061/(ASCE)ME.1943-5479.0000730
|
| [33] |
J. Pasha, M. A. Dulebenets, A. M. Fathollahi-Fard, G. Tian, Y. Y. Lau, P. Singh, B. Liang, An integrated optimization method for tactical-level planning in liner shipping with heterogeneous ship fleet and environmental considerations, Adv. Eng. Inform., 48 (2021), 101299. https://doi.org/10.1016/j.aei.2021.101299 doi: 10.1016/j.aei.2021.101299
|
| [34] |
J. Pasha, A. L. Nwodu, A. M. Fathollahi-Fard, G. Tian, Z. Li, H. Wang, et al., Exact and metaheuristic algorithms for the vehicle routing problem with a factory-in-a-box in multi-objective settings, Adv. Eng. Inform., 52 (2022), 101623. https://doi.org/10.1016/j.aei.2022.101623 doi: 10.1016/j.aei.2022.101623
|
| [35] |
H. Zhao, C. Zhang, An online-learning-based evolutionary many-objective algorithm, Inf. Sci., 509 (2020), 1–21. https://doi.org/10.1016/j.ins.2019.08.069 doi: 10.1016/j.ins.2019.08.069
|
| [36] |
M. A. Dulebenets, An adaptive polyploid memetic algorithm for scheduling trucks at a cross-docking terminal, Inf. Sci., 565 (2021), 390–421. https://doi.org/10.1016/j.ins.2021.02.039 doi: 10.1016/j.ins.2021.02.039
|
| [37] |
M. Rabbani, N. Oladzad-Abbasabady, N. Akbarian-Saravi, Ambulance routing in disaster response considering variable patient condition: NSGA-Ⅱ and MOPSO algorithms, J. Ind. Manage. Optim., 18 (2022), 1035. https://doi.org/10.3934/jimo.2021007 doi: 10.3934/jimo.2021007
|
| [38] |
T. Y. Chen, Novel generalized distance measure of Pythagorean fuzzy sets and a compromise approach for multiple criteria decision analysis under uncertainty, IEEE Access, 7 (2019), 58168–58185. https://doi.org/10.1109/ACCESS.2019.2914703 doi: 10.1109/ACCESS.2019.2914703
|
| [39] |
H. Garg, D. Rani, Robust averaging-geometric aggregation operators for complex intuitionistic fuzzy sets and their applications to MCDM process, Arab. J. Sci. Eng. 45 (2020), 2017–2033. https://doi.org/10.1007/s13369-019-03925-4 doi: 10.1007/s13369-019-03925-4
|
| [40] |
C. Li, S. Zeng, T. Pan, L. Zheng, A method based on induced aggregation operators and distance measures to multiple attribute decision making under 2-tuple linguistic environment, J. Comput. Syst. Sci., 80 (2014), 1339–1349. https://doi.org/10.1016/j.jcss.2014.03.004 doi: 10.1016/j.jcss.2014.03.004
|
| [41] |
J. H. Park, J. M. Park, Y. C. Kwun, 2-tuple linguistic harmonic operators and their applications in group decision making, Knowl. Based Syst., 44 (2013), 10–19. https://doi.org/10.1016/j.knosys.2013.01.006 doi: 10.1016/j.knosys.2013.01.006
|
| [42] |
G. W. Wei, Some generalized aggregating operators with linguistic information and their application to multiple attribute group decision making, Comput. Ind. Eng., 61 (2011), 32–38. https://doi.org/10.1016/j.cie.2011.02.007 doi: 10.1016/j.cie.2011.02.007
|
| [43] |
F. Zhou, T. Y. Chen, An integrated multicriteria group decision-making approach for green supplier selection under Pythagorean fuzzy scenarios, IEEE Access, 8 (2020), 165216–165231. https://doi.org/10.1109/ACCESS.2020.3022377 doi: 10.1109/ACCESS.2020.3022377
|
| [44] |
P. Liu, S. Naz, M. Akram, M. Muzammal, Group decision-making analysis based on linguistic q-rung orthopair fuzzy generalized point weighted aggregation operators, Int. J. Mach. Learn. Cybern., 13 (2021), 883–906. https://doi.org/10.1007/s13042-021-01425-2 doi: 10.1007/s13042-021-01425-2
|
| [45] | H. Hamacher, Über logische verknünpfungen unscharfer Aussagen und deren zugehörige Bewertungsfunktionen, in Progress in Cybernetics and Systems Research (Eds. Trappl, Klir, Riccardi), Hemisphere, Washington DC, 3 (1978), 276–288. |
| [46] |
S. Faizi, W. Salabun, S. Nawaz, A. U. Rehman, J. W
|
| [47] |
S. Abdullah, O. Barukab, M. Qiyas, M. Arif, S. A. Khan, Analysis of decision support system based on 2-tuple spherical fuzzy linguistic aggregation information, Appl. Sci., 10 (2020), 276. https://doi.org/10.3390/app10010276 doi: 10.3390/app10010276
|
| [48] |
M. Akram, G. Shahzadi, J. C. R. Alcantud, Multi-attribute decision-making with q-rung picture fuzzy information, Granul. Comput., 7 (2022), 197–215. https://doi.org/10.1007/s41066-021-00260-8 doi: 10.1007/s41066-021-00260-8
|
| [49] |
X. P. Jiang, G. W. Wei, Some Bonferroni mean operators with 2-tuple linguistic information and their application to multiple attribute decision making, J. Intell. Fuzzy Syst., 27 (2014), 2153–2162. https://doi.org/10.3233/IFS-141180 doi: 10.3233/IFS-141180
|
| [50] |
N. Waseem, M. Akram, J. C. R. Alcantud, Multi-attribute decision-making based on m-polar fuzzy Hamacher aggregation operators, Symmetry, 11 (2019), 1498. https://doi.org/10.3390/sym11121498 doi: 10.3390/sym11121498
|
| [51] |
Z. Zhang, F. Wei, S. Zhou, Approaches to comprehensive evaluation with 2-tuple linguistic information, J. Intell. Fuzzy Syst., 28 (2015), 469–475. https://doi.org/10.3233/IFS-141323 doi: 10.3233/IFS-141323
|
| [52] |
H. Jin, S. Ashraf, S. Abdullah, M. Qiyas, M. Bano, S. Zeng, Linguistic spherical fuzzy aggregation operators and their applications in multi-attribute decision making problems, Mathematics, 7 (2019), 413. https://doi.org/10.3390/math7050413 doi: 10.3390/math7050413
|
| [53] |
Y. Ju, A. Wang, J. Ma, H. Gao, E. D. Santibanez Gonzalez, Some q-rung orthopair fuzzy 2-tuple linguistic Muirhead mean aggregation operators and their applications to multiple-attribute group decision making, Int. J. Intell. Syst., 35 (2020), 184–213. https://doi.org/10.1002/int.22205 doi: 10.1002/int.22205
|
| [54] |
X. Deng, J. Wang, G. Wei, Some 2-tuple linguistic Pythagorean Heronian mean operators and their application to multiple attribute decision-making, J. Exp. Theor. Artif. Intell., 31 (2019), 555–574. https://doi.org/10.1080/0952813X.2019.1579258 doi: 10.1080/0952813X.2019.1579258
|
| [55] |
X. Deng, J. Wang, G. Wei, M. Lu, Models for multiple attribute decision making with some 2-tuple linguistic Pythagorean fuzzy Hamy mean operators, Mathematics, 6 (2018), 236. https://doi.org/10.3390/math6110236 doi: 10.3390/math6110236
|
| [56] |
G. W. Wei, 2-tuple intuitionistic fuzzy linguistic aggregation operators in multiple attribute decision making, Iran. J. Fuzzy Syst., 16 (2019), 159–174. https://doi.org/10.22111/IJFS.2019.4789 doi: 10.22111/IJFS.2019.4789
|
| [57] |
M. Lu, G. Wei, F. E. Alsaadi, T. Hayat, A. Alsaedi, Bipolar 2-tuple linguistic aggregation operators in multiple attribute decision making, J. Intell. Fuzzy Syst., 33 (2017), 1197–1207. https://doi.org/10.3233/JIFS-16946 doi: 10.3233/JIFS-16946
|
| [58] |
Y. Zhang, G. Wei, Y. Guo, C. Wei, TODIM method based on cumulative prospect theory for multiple attribute group decision-making under 2-tuple linguistic Pythagorean fuzzy environment, Int. J. Intell. Syst., 36 (2021), 2548–2571. https://doi.org/10.1002/int.22393 doi: 10.1002/int.22393
|
| [59] |
M. Akram, R. Bibi, M. A. Al-Shamiri, A decision-making framework based on 2-tuple linguistic Fermatean fuzzy Hamy mean operators, Math. Probl. Eng., 2022 (2022), 1501880. https://doi.org/10.1155/2022/1501880 doi: 10.1155/2022/1501880
|
| [60] |
M. Akram, A. Khan, U. Ahmad, Extended MULTIMOORA method based on 2-tuple linguistic Pythagorean fuzzy sets for multi-attribute group decision-making, Granul. Comput., 2022 (2022), 1–22. https://doi.org/10.1007/s41066-022-00330-5. doi: 10.1007/s41066-022-00330-5
|
| [61] |
D. Ramot, M. Friedman, G. Langholz, G. A. Kandel, Complex fuzzy logic, IEEE Trans. Fuzzy Syst., 11 (2003), 450–461. https://doi.org/10.1109/TFUZZ.2003.814832 doi: 10.1109/TFUZZ.2003.814832
|
| [62] |
D. Ramot, R. Milo, M. Fiedman, A. Kandel, Complex fuzzy sets, IEEE Trans. Fuzzy Syst., 10 (2002), 171–186. https://doi.org/10.1109/91.995119 doi: 10.1109/91.995119
|
| [63] | A. M. D. J. S. Alkouri, A. R. Salleh, Complex intuitionistic fuzzy sets, in AIP Conference Proceedings, 1482 (2012), 464–470. https://doi.org/10.1063/1.4757515 |
| [64] |
A. U. M. Alkouri, A. R. Salleh, Complex Atanassov's intuitionistic fuzzy relation, Abstr. Appl. Anal., 2013 (2013), 1–18. https://doi.org/10.1155/2013/287382 doi: 10.1155/2013/287382
|
| [65] |
Y. Rong, Y. Liu, Z. Pei, Complex q-rung orthopair fuzzy 2-tuple linguistic Maclaurin symmetric mean operators and its application to emergency program selection, Int. J. Intell. Syst., 35 (2020), 1749–1790. https://doi.org/10.1002/int.22271 doi: 10.1002/int.22271
|
| [66] |
L. Bi, S. Dai, B. Hu, Complex fuzzy geometric aggregation operators, Symmetry, 10 (2018), 251, https://doi.org/10.3390/sym10070251. doi: 10.3390/sym10070251
|
| [67] |
L. Bi, S. Dai, B. Hu, S. Li, Complex fuzzy arithmetic aggregation operators, J. Intell. Fuzzy Syst., 36 (2019), 2765–2771. https://doi.org/10.3233/JIFS-18568 doi: 10.3233/JIFS-18568
|
| [68] |
P. Liu, T. Mahmood, Z. Ali, Complex q-rung orthopair fuzzy aggregation operators and their applications in multiattribute group decision making, Information, 11 (2020), 5. https://doi.org/10.3390/info11010005 doi: 10.3390/info11010005
|
| [69] |
A. Luqman, M. Akram, A. N. Al-Kenani, J. C. R. Alcantud, A study on hypergraph representations of complex fuzzy information, Symmetry, 11 (2019), 1381. https://doi.org/10.3390/sym11111381 doi: 10.3390/sym11111381
|
| [70] |
S. Naz, M. Akram, M. M. A. Al-Shamiri, M. M. Khalaf, G. Yousaf, A new MAGDM method with 2-tuple linguistic bipolar fuzzy Heronian mean operators, Math. Biosci. Eng., 19 (2022), 3843–3878. https://doi.org/10.3934/mbe.2022177 doi: 10.3934/mbe.2022177
|
| [71] |
S. Naz, M. Akram, M. M. A. Al-Shamiri, M. R. Saeed, Evaluation of network security service provider using 2-tuple linguistic complexq-rung orthopair fuzzy COPRAS method, Complexity, 2022 (2022), 1–27. https://doi.org/10.1155/2022/4523287 doi: 10.1155/2022/4523287
|
| [72] |
P. Liu, Z. Ali, T. Mahmood, Generalized complex q-rung orthopair fuzzy Einstein averaging aggregation operators and their application in multi-attribute decision making, Complex Intell. Syst., 7 (2021), 511–538. https://doi.org/10.1007/s40747-020-00197-6 doi: 10.1007/s40747-020-00197-6
|
| [73] |
M. Akram, A. Bashir, S. A. Edalatpanah, A hybrid decision-making analysis under complex q-rung picture fuzzy Einstein averaging operators, Comput. Appl. Math., 40 (2021), 1–35. https://doi.org/10.1007/s40314-021-01651-y doi: 10.1007/s40314-021-01651-y
|
| [74] |
M. Akram, X. Peng, A. Sattar, Multi-criteria decision-making model using complex Pythagorean fuzzy Yager aggregation operators, Arab. J. Sci. Eng., 46 (2021), 1691–1717. https://doi.org/10.1007/s13369-020-04864-1 doi: 10.1007/s13369-020-04864-1
|
| [75] |
H. Garg, D. Rani, Some generalized complex intuitionistic fuzzy aggregation operators and their application to multicriteria decision-making process, Arab. J. Sci. Eng., 44 (2019), 2679–2698. https://doi.org/10.1007/s13369-018-3413-x doi: 10.1007/s13369-018-3413-x
|
| [76] |
P. Liu, Z. Ali, T. Mahmood, Novel complex T-spherical fuzzy 2-tuple linguistic Muirhead mean aggregation operators and their application to multi-attribute decision-making, Int. J. Comput. Intell. Syst., 14 (2021), 295–331. https://doi.org/10.2991/ijcis.d.201207.003 doi: 10.2991/ijcis.d.201207.003
|
| [77] |
M. Akram, S. Naz, F. Feng, A. Shafiq, Assessment of hydropower plants in Pakistan: Muirhead mean-based 2-tuple linguistic T-spherical fuzzy model combining SWARA with COPRAS, Arabian J. Sci. Eng., 2022 (2022), 1–30. https://doi.org/10.1007/s13369-022-07081-0 doi: 10.1007/s13369-022-07081-0
|
| [78] |
M. Akram, N. Ramzan, F. Feng, Extending COPRAS method with linguistic Fermatean fuzzy sets and Hamy mean operators, J. Math., 2022 (2022), 8239263. https://doi.org/10.1155/2022/8239263 doi: 10.1155/2022/8239263
|
| [79] |
T. Mahmood, Z. Ali, A novel approach of complex q-rung orthopair fuzzy Hamacher aggregation operators and their application for cleaner production assessment in gold mines, J. Ambient Intell. Humaniz. Comput., 12 (2021), 8933–8959. https://doi.org/10.1007/s12652-020-02697-2 doi: 10.1007/s12652-020-02697-2
|
| [80] |
D. Rani, H. Garg, Complex intuitionistic fuzzy power aggregation operators and their applications in multicriteria decision-making, Expert Syst., 35 (2018), e12325. https://doi.org/10.1111/exsy.12325 doi: 10.1111/exsy.12325
|
| 1. | Nam Bui Duc, Nguyen Minh Hai, Luu Vu Cam Hoan, Le Dinh Long, On inverse source term for heat equation with memory term, 2024, 57, 2391-4661, 10.1515/dema-2023-0138 | |
| 2. | Aslam Khan, Abdul Ghafoor, Emel Khan, Manzoor Hussain, Haar Wavelet based method of lines for the numerical solutions of inverse problems with space dependent source and noise intensity parameter on the boundary, 2025, 100, 0031-8949, 045251, 10.1088/1402-4896/adc20d |
Figures(7) / Tables(12)
Muhammad Akram, Ayesha Khan, Uzma Ahmad, José Carlos R. Alcantud, Mohammed M. Ali Al-Shamiri. A new group decision-making framework based on 2-tuple linguistic complex q-rung picture fuzzy sets[J]. Mathematical Biosciences and Engineering, 2022, 19(11): 11281-11323. doi: 10.3934/mbe.2022526
| Name | Configuration |
| Operating System | Ubuntu 18.04.5 |
| Framework | PyTorch 1.8.0 |
| Memory | 60 Gb |
| GPU | GeForce RTX 3080Ti |
| CUDA | CUDA 10.1 |
DownLoad:
CSV
| Number of HDC layers | MAE (km) |
| 1 | 78.22 |
| 2 | 66.38 |
| 3 | 69.43 |
DownLoad:
CSV
| n | MAE (km) |
| 5 | 96.05 |
| 10 | 68.49 |
| 20 | 79.43 |
DownLoad:
CSV
| Name | Parameter |
| HDC layers | 2 |
| Deconvolution layers | 6 |
| LSTM layers | 2 |
| Kernel | 3 × 3 |
DownLoad:
CSV
| Sequence | Real | Predicted | E (km) |
| 11(2007112300) | 13.7, 127.2 | 13.6, 126.8 | 44.629 |
| 12(2007112306) | 13.7, 127.0 | 14, 127.1 | 35.062 |
| 13(2007112312) | 14.0, 126.7 | 14.5, 126.6 | 56.632 |
| 14(2007112318) | 14.1, 126.5 | 14.3, 126.1 | 48.516 |
| 15(2007112400) | 14.3, 126.2 | 14.8, 126.4 | 59.619 |
| 16(2007112406) | 14.5, 125.9 | 15, 126.2 | 64.278 |
| 17(2007112412) | 14.8, 125.5 | 15, 124.9 | 68.201 |
| 18(2007112418) | 15.1, 125.2 | 15.4, 125.5 | 46.353 |
| 19(2007112500) | 15.5, 124.5 | 16.1, 124.7 | 70.065 |
| 20(2007112506) | 16.4, 123.9 | 16.1, 124.1 | 39.606 |
| 21(2007112512) | 16.9, 122.9 | 16.8, 122.5 | 43.997 |
| 22(2007112518) | 17.3, 121.5 | 17, 122.1 | 71.951 |
| 23(2007112600) | 18.3, 121.0 | 18.7, 121.4 | 61.298 |
| 24(2007112606) | 19.0, 120.7 | 19.4, 120.5 | 49.187 |
| 25(2007112612) | 19.5, 120.7 | 19.2, 121.1 | 53.609 |
| 26(2007112618) | 20.0, 121.1 | 20.6, 121.4 | 73.688 |
| 27(2007112700) | 20.1, 121.5 | 20.4, 122.6 | 56.719 |
| 28(2007112706) | 20.5, 122.8 | 20.6, 123.6 | 23.619 |
| 29(2007112712) | 20.6, 124.0 | 20, 124.4 | 78.685 |
| MAE (km) | |
|
55.037 |
DownLoad:
CSV
| Sequence | Real | Predicted | E (km) |
| 11(2015021909) | –20.5, 150.6 | –20.4, 150.3 | 33.175 |
| 12(2015021912) | –20.8, 150.5 | –20.6, 150.9 | 47.177 |
| 13(2015021915) | –21.1, 150.5 | –21.5, 150.8 | 54.261 |
| 14(2015021918) | –21.5, 150.5 | –21.5, 151.1 | 70.413 |
| 15(2015021921) | –22.1, 150.5 | –21.9, 150.9 | 46.853 |
| 16(2015022000) | –22.7, 150.5 | –22.5, 150.6 | 24.494 |
| 17(2015022003) | –23.1, 150.5 | –22.7, 150.2 | 54.061 |
| 18(2015022006) | –23.7, 150.6 | –23.3, 150.9 | 53.983 |
| 19(2015022009) | –24.2, 150.9 | –24.7, 150.2 | 90.065 |
| 20(2015022012) | –24.5, 151.3 | –24.9, 150.8 | 67.302 |
| 21(2015022015) | –24.9, 151.6 | –25.1, 151 | 64.426 |
| 22(2015022018) | –25.3, 151.8 | –25.4, 151.4 | 41.705 |
| MAE (km) | |
|
53.992 |
DownLoad:
CSV
| Data set | Test set size | RMSE (km) | MAE (km) | Computational time (s) |
| Northwest Pacific | 1610 | 72.72 | 65.39 | 87.42 |
| Southwest Pacific | 840 | 69.58 | 62.41 | 49.80 |
DownLoad:
CSV
| Name | Configuration |
| Operating System | Ubuntu 18.04.5 |
| Framework | PyTorch 1.8.0 |
| Memory | 60 Gb |
| GPU | GeForce RTX 3080Ti |
| CUDA | CUDA 10.1 |
| Number of HDC layers | MAE (km) |
| 1 | 78.22 |
| 2 | 66.38 |
| 3 | 69.43 |
| n | MAE (km) |
| 5 | 96.05 |
| 10 | 68.49 |
| 20 | 79.43 |
| Name | Parameter |
| HDC layers | 2 |
| Deconvolution layers | 6 |
| LSTM layers | 2 |
| Kernel | 3 × 3 |
| Sequence | Real | Predicted | E (km) |
| 11(2007112300) | 13.7, 127.2 | 13.6, 126.8 | 44.629 |
| 12(2007112306) | 13.7, 127.0 | 14, 127.1 | 35.062 |
| 13(2007112312) | 14.0, 126.7 | 14.5, 126.6 | 56.632 |
| 14(2007112318) | 14.1, 126.5 | 14.3, 126.1 | 48.516 |
| 15(2007112400) | 14.3, 126.2 | 14.8, 126.4 | 59.619 |
| 16(2007112406) | 14.5, 125.9 | 15, 126.2 | 64.278 |
| 17(2007112412) | 14.8, 125.5 | 15, 124.9 | 68.201 |
| 18(2007112418) | 15.1, 125.2 | 15.4, 125.5 | 46.353 |
| 19(2007112500) | 15.5, 124.5 | 16.1, 124.7 | 70.065 |
| 20(2007112506) | 16.4, 123.9 | 16.1, 124.1 | 39.606 |
| 21(2007112512) | 16.9, 122.9 | 16.8, 122.5 | 43.997 |
| 22(2007112518) | 17.3, 121.5 | 17, 122.1 | 71.951 |
| 23(2007112600) | 18.3, 121.0 | 18.7, 121.4 | 61.298 |
| 24(2007112606) | 19.0, 120.7 | 19.4, 120.5 | 49.187 |
| 25(2007112612) | 19.5, 120.7 | 19.2, 121.1 | 53.609 |
| 26(2007112618) | 20.0, 121.1 | 20.6, 121.4 | 73.688 |
| 27(2007112700) | 20.1, 121.5 | 20.4, 122.6 | 56.719 |
| 28(2007112706) | 20.5, 122.8 | 20.6, 123.6 | 23.619 |
| 29(2007112712) | 20.6, 124.0 | 20, 124.4 | 78.685 |
| MAE (km) | |
|
55.037 |
| Sequence | Real | Predicted | E (km) |
| 11(2015021909) | –20.5, 150.6 | –20.4, 150.3 | 33.175 |
| 12(2015021912) | –20.8, 150.5 | –20.6, 150.9 | 47.177 |
| 13(2015021915) | –21.1, 150.5 | –21.5, 150.8 | 54.261 |
| 14(2015021918) | –21.5, 150.5 | –21.5, 151.1 | 70.413 |
| 15(2015021921) | –22.1, 150.5 | –21.9, 150.9 | 46.853 |
| 16(2015022000) | –22.7, 150.5 | –22.5, 150.6 | 24.494 |
| 17(2015022003) | –23.1, 150.5 | –22.7, 150.2 | 54.061 |
| 18(2015022006) | –23.7, 150.6 | –23.3, 150.9 | 53.983 |
| 19(2015022009) | –24.2, 150.9 | –24.7, 150.2 | 90.065 |
| 20(2015022012) | –24.5, 151.3 | –24.9, 150.8 | 67.302 |
| 21(2015022015) | –24.9, 151.6 | –25.1, 151 | 64.426 |
| 22(2015022018) | –25.3, 151.8 | –25.4, 151.4 | 41.705 |
| MAE (km) | |
|
53.992 |
| Data set | Test set size | RMSE (km) | MAE (km) | Computational time (s) |
| Northwest Pacific | 1610 | 72.72 | 65.39 | 87.42 |
| Southwest Pacific | 840 | 69.58 | 62.41 | 49.80 |
| Model | Predicted Time Range | RMSE (km) | MAE (km) |
| RNN[27] | 12 h | 164.22 | |
| AE-GRU[28] | 12 h | 138.67 | |
| TraijGRU[29] | 6 h | 66.6 | |
| DeepTyphoon | 6 h | 73.96 | 64.17 |
| 12 h | 134.70 | 129.82 |
| Model | Latitude | Longitude | ||
| 6 h | 1 h | 6 h | 12 h | |
| DeepFR[30] | 0.9100 | 1.1765 | 0.8910 | 1.2430 |
| Trj-DMFMG[31] | 0.6602 | 0.9275 | 0.8747 | 1.1329 |
| DeepTyphoon | 0.4876 | 0.8812 | 0.6645 | 1.0874 |
