• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    VGG-CovidNet:Bi-Branched Dilated Convolutional Neural Network for Chest X-Ray-Based COVID-19 Predictions

    2021-12-11 13:32:58MuhammedBinsawadMarwanAlbaharandAbdullahBinSawad
    Computers Materials&Continua 2021年8期

    Muhammed Binsawad,Marwan Albahar and Abdullah Bin Sawad

    1Information Systems,King Abdulaziz University,Jeddah,Saudi Arabia

    2Department of Science,Umm Al Qura University,Mecca,Saudi Arabia

    Abstract: The coronavirus disease 2019 (COVID-19) pandemic has had a devastating impact on the health and welfare of the global population.A key measure to combat COVID-19 has been the effective screening of infected patients.A vital screening process is the chest radiograph.Initial studies have shown irregularities in the chest radiographs of COVID-19 patients.The use of the chest X-ray(CXR),a leading diagnostic technique,has been encouraged and driven by several ongoing projects to combat this disease because of its historical effectiveness in providing clinical insights on lung diseases.This study introduces a dilated bi-branched convoluted neural network (CNN)architecture,VGG-COVIDNet,to detect COVID-19 cases from CXR images.The front end of the VGG-COVIDNet consists of the first 10 layers of VGG-16, where the convolutional layers in these layers are reduced to two to minimize latency during the training phase.The last two branches of the proposed architecture consist of dilated convolutional layers to reduce the model’s computational complexity while retaining the feature maps’spatial information.The simulation results show that the proposed architecture is superior to all the state-of-the-art architecture in accuracy and sensitivity.The proposed architecture’s accuracy and sensitivity are 96.5% and 96%,respectively,for each infection type.

    Keywords:Coronavirus disease 2019; prognosis; X-ray images; deep learning; artificial intelligence

    1 Introduction

    Coronavirus disease 2019 (COVID-19) has been spreading exponentially; thus, it has become a pandemic.Its transmission process has not been fully comprehended.The virus often causes few or no symptoms; however, in 2%-8% of infected individuals, rapidly advancing and, frequently,deadly pneumonia can occur [1-3].The precise prevalence, mortality, and transmission approaches are still not completely understood partly because of the distinct issues posed by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection.They include peak infectiousness at or immediately preceding symptom onset and a poorly understood multi-organ pathophysiology with major features and fatality within the lungs [4].Because of the rapid transmission rate of infection, healthcare systems worldwide have been affected by the unavailability of key protective gear and qualified providers [5].In addition, practitioners have varied access to point-of-care testing techniques, such as reverse transcription-polymerase chain reaction (RT-PCR) [6].Despite the increase in the availability of rapid RT-PCR testing, challenges persist.They include high rates of false negatives, delays in processing, changes in testing methodologies, and sensitivity rates that are often found to be low [7,8].

    Radiography is a technique that has been applied to COVID-19 screening.Chest radiography(e.g., computed tomography [CT] or [CXR]) is performed and read by radiologists to visually identify infection-related features.Initial studies have shown irregularities in patients’chest radiographs.This is typically the case in those infected with COVID-19 [9-11].Studies have indicated that radiographic assessments could serve as the main technique for COVID-19 screening in epidemic regions [12].For example, Guan et al.[11] found that radiographic irregularities, including bilateral abnormalities, interstitial irregularities, and ground-glass opacity in CT and CXR images, were present in COVID-19 positive cases.Huang et al.[10] found bilateral radiographic inconsistencies in the CXR images of most COVID-19 positive cases.Most new studies have focused on CT imaging because of the higher image detail.However, the use of CXR imaging for COVID-19 screening during the pandemic has several benefits, especially for resource-limited and highly affected regions.

    Rapid training:Through CXR images, patients with suspected COVID-19 infection can be rapidly trained.This can be performed along with viral testing (the results of which take time)to bring relief to the large number of patients, particularly those in the highly affected regions(e.g., Spain and Italy) where hospitals have exceeded capacity.CXRs can also be performed alone when viral testing is not possible because of low supplies.In addition, CXR imaging is highly effective for triaging in geographical regions (e.g., in New York City) where patients are given instructions to stay at home until they start experiencing more advanced symptoms.The reason is that abnormalities are frequently observed in patients presenting with suspected COVID-19 at clinics [13].

    Availability and accessibility:Many clinics and imaging centers have easily available and accessible CXR imaging because most healthcare systems consider it standard equipment.Particularly in developing countries, CT imaging is less available than CXR imaging because of the high acquisition and maintenance costs.

    Portability:Imaging can be performed in isolation rooms because of the availability of portable CXR systems.This considerably decreases the risk of COVID-19 transmission in rooms with fixed imaging systems or during patient transport to these rooms [13].

    Therefore, it is possible to perform radiographic assessments more rapidly.In addition, these systems have become more widely available, as is clear from the availability of chest radiographic systems and portable equipment in modern healthcare systems.Hence, computer-assisted diagnostic systems that can improve the accuracy of radiologists’interpretations of images to determine the status of COVID-19 cases are highly desirable.Several models have been proposed and implemented.Most have limitations, such as computational overhead, a limited number of training or testing sets, the requirement for a separate training session, and the need to train a large number of parameters.Thus, the development of accurate diagnostic methods is important.

    Consequently, the present study has proposed VGG-COVIDNet, a dilated bi-branched convoluted neural network (CNN) system for COVID-19 detection on the basis of CXR images.The motivation is the growing need for solutions to combat the pandemic.A review of the literature indicates that COVIDx is the most extensive open-access benchmark dataset on the number of publicly available COVID-19 positive cases.Furthermore, the predictive effectiveness of VGG-COVIDNet was investigated by leveraging the transfer-learning capability of the first 10 layers of VGG-16 and dilated convolution to capture the feature maps’spatial information to reduce the complexity of the model.

    With the rapid increase in the number of suspected cases, artificial intelligence methods may play a role in the identification and characterization of COVID-19 through imagery.CXR gives this process a clear and expedited window.Deep learning from large CXR multinationals could provide automated reproducible biomarkers to classify and to quantify COVID-19.This study aimed to develop and to evaluate a CNN-based model to classify three types of infections.The focus was the high-accuracy detection of normal, COVID-19, and non-COVID-19 cases.The rationale behind the choice of these categories was the opportunity for clinicians to prioritize patients for COVID-19 testing and to determine the appropriate care on the basis of the cause of infection because of the differences in the treatment of COVID19 and non-COVID-19 infections.The main contributions of this work are as follows:

    1.A dilated convolution and branching strategy model are proposed.The dilated convolution aims to grab spatial information from the feature maps by avoiding pooling layers, which increase complexity.

    2.Given the nature of the problem, the first 10 layers of the VGG-16 model [14] are used as the front end of the architecture because of the model’s flexibility and transfer-learning capability.These first 10 layers extract and use the high-level feature and provide the generated outputs to the two branches to further extract and to incorporate low-level features throughout the training process.

    This paper is organized as follows.Section 2 discusses the methods employed in the detection of COVID-19 and other diseases.Section 3 explains the use of the dataset for infection detection.In Section 4, the benefits of the proposed architecture and dilated convolution are discussed.Section 5 explains the efficiency of the architecture through evaluation matrices and provides a comparison with the state-of-the-art methods.Finally, conclusions are drawn, and future directions are discussed in Section 6.

    2 Related Work

    Studies have examined the potential of CXR and CT scans in the detection of COVID-19 [15-17].Most have sought to generate automatic COVID-19 classification systems that typically employ CNNs [18].Xu et al.[19] initially used a pre-trained three-dimensional CNN to extract potentially infected regions via CT scans.The samples were then added to another CNN and sorted with 86.7% accuracy into three categories:non-infection, influenza viral pneumonia, and COVID-19.Wang et al.[20] first extracted candidates by employing a thresholdbased strategy.Next, for every case, a random selection of two or three regions was used to create the dataset, which was used to fine-tune a pre-trained CNN.Finally, features were removed from the CNN and added to a group of classifiers for predicting COVID-19 with 88% accuracy.The authors in [21] used CT scans to identify positive COVID-19 cases in which all the slices were separately fed into the model and a max-pooling function was used to combine all the outputs,thereby achieving 90% sensitivity.In [22], a pretraining CNN model was initially trained on the ImageNet dataset [23].Next, it was refined through the use of a CXR image dataset to classify the samples as normal, bacterial, COVID-19 viral infection, and non-COVID-19 viral infection with 83.5% accuracy.Sethy et al.[24] conducted similar research in which they trained different CNN models on CXR images, then an SVM classifier used to detect positive COVID-19 cases with 95.38% accuracy.In [25], capsule networks, COVID-CAPS, were proposed.This model, which includes capsule and convolution layers, attained an accuracy of 95.7%, a sensitivity of 90%, a specificity of 95.8%, and an area under the curve (AUC) of 0.97%.

    CXRs have been widely used to identify COVID-19 positive cases [22,24,26,27].Chen et al.[28] found bilateral pneumonia in patient CXRs in a sample of 99 cases.Interestingly, the link between CXRs and CT images was noted by Yoon et al.[26] in a study of nine COVID-19 positive cases.Other studies have focused on neural network-tailored deep-learning (DL) models,such as ResNet50 [24] and COVID-Net [22].ResNet50 was tested on 25 positive COVID-19 cases,and COVID-Net was used on just 31 positive cases.A classic deep-learning model was employed by Zhang et al.[27] on 100 COVID-19 samples to detect positive cases.This method exhibited a high accuracy of 90%.To detect COVID-19 positive cases from CXRs, Ozturk et al.[29]developed a deep neural network (DNN)-based method that achieved 98.08% accuracy in 125 positive cases.Moreover, the model also achieved 87.02% accuracy by employing the multi-class scenario.Narin et al.[16] used three types of DNN architecture (Inception-ResNetV2, ResNet50,and InceptionV3) to identify COVID-19 positive cases from a total of 50 samples.They achieved accuracy levels of 97%, 98%, and 87%, respectively.Using a sample of 155 COVID-19 positive cases, Mangal et al.[30] proposed a DL-based system, CovidAID, to determine positive cases from CXRs.The study achieved 90.5% accuracy and 100% sensitivity.

    To detect COVID-19 cases from CT images with weak labels, Wang et al.[31] proposed a DLbased method, which was trained by using 499 volumes and tested with an additional 131 volumes.The authors achieved specificity and sensitivity values of 0.911 and 0.907, respectively.Farooq et al.[17] proposed a DL-based approach to distinguish COVID-19 cases on the basis of CXRs.Using eight COVID-19 positive cases, the study achieved 96.23% accuracy.A DL-based method was used by Hall et al.[32] to determine COVID-19 cases from CXRs.From 135 COVID-19 positive cases, the authors achieved 89.2% accuracy, exhibiting a true positivity rate of 0.8039 and an AUC of 0.95.An ensemble-based method that was used to test a group of 33 CXRs achieved 91.24% accuracy and an AUC of 0.7879.The authors in [33] employed CXRs and CT scans to test a proposed CNN-tailored DNN, with 96.28% accuracy (false ?negativerate=0.0208,AUC=0.9808).While the authors in [34] used traditional machine learning and DL to facilitate COVID-19 detection.The authors claimed that this was the first approach to combine both types of learning to detect COVID-19 infection by learning CXR image representations of highly discriminatory features.They stated that the DL models had achieved the optimum accuracy of 98.8% in the ResNet50 model.In [35], the authors proposed a parallel architecture namely COVID-CheXNet to provide a high degree of confidence in discriminating between healthy and COVID-19 infected individuals.The study used two pre-trained DL models based on the incorporation of a ResNet34 and a high-resolution network model.Given the 99.99% accuracy rate,the authors claimed that the proposed system had correctly diagnosed the COVID-19 patients.

    3 COVIDx Dataset

    The present study used the COVIDx dataset proposed by [22].It consisted of 13,975 CXR images across 13,870 patient cases.The dataset was generated through the integration and modification of five publicly available data archives:

    1.COVID-19 Image Data Collection [36]

    2.COVID-19 Chest X-Ray Dataset Initiative [37]

    3.ActualMed COVID-19 Chest X-Ray Dataset Initiative [38]

    4.RSNA Pneumonia Detection Challenge Dataset [39]

    5.COVID-19 Radiography Database [40].

    Some images from the COVIDx dataset are shown in Fig.1.All the combined images from the above archives were divided into the following three categories.

    Figure 1:Sample images from the COVIDx dataset.(1) COVID-19 images, (2) non-COVID-19 images, (3) normal images

    1.Normal:This comprised 7,983 normal images, i.e., the absence of chest disease.

    2.Non-COVID 19:This contained 5,475 images with signs of non-COVID-19-related chest disease.

    3.COVID-19:This included 517 images of COVID-19 patients.

    The data distribution by infection type is presented in Fig.2.It illustrates the paucity of COVID-19 patient data in the publicly available repositories.Fig.2 also indicates the need for additional COVID-19 data to improve the results and to reduce the false-negative and falsepositive rates in the testing phase.Tab.1 presents the distribution of the images in the training and testing sets.

    Figure 2:Distribution of COVIDx dataset in training and testing datasets.The values are represented on the y-axis after taking log with base 10

    Table 1:Image distribution in training and testing sets

    4 Methods

    This study presents a hybrid bi-branched CNN model based on dilated convolution to classify CXR images in three infection types:COVID-19, non-COVID-19, and normal.The first 10 layers of the VGG-16 architecture were used as the front end of the proposed model, VGG-COVIDNet.The back end was based on dilated kernels, which are explained in the succeeding sections.The main idea of the proposed solution was the deployment of a deeper CNN to produce high-level features without losing detailed information in the images.In addition, it has large receptive fields and classifies the CXRs in the respective categories without the negative exploitation of network complexity.This section describes the proposed architecture and presents the training method.

    4.1 Pre-Processing

    Initially, the input images were converted to the JPEG format and resized to 480 × 480 ×3 dimension.Similarly, the normalization of images was performed in the 0 to 1 range based on the formulation given in Eq.(1).Next, the power law transformationwas performed in accordance with Eq.(2) withγ= 0.6 to make the dark areas of the images more prominent to reveal additional information about the corresponding disease.

    Here,IoutandIinare the output and input images, respectively, wherec=1 andγ=0.6.

    4.2 Architecture of the Proposed VGG-CovidNet Model

    The proposed model consists of VGG-16 [14] and some CNN layers with expanded convolutional layers to classify the CXR images in the respective classes.VGG-16 was chosen as the front-end architecture on the basis of [41-43] because of its strong transfer learning ability and the flexibility to easily concatenate the back end for the desired task.In the absence of a back end,the VGG-16 performs supportive tasks without boosting performance, and this results in very low accuracy.In the present study, the fully connected VGG-16 layers (classification part) are removed,and the proposed VGG-CovidNet is made by appending the convolutional layers in VGG-16.The front-end network’s output size is one-eighth of the original input size, and upon stacking multiple layers, the information in the images is lost.In turn, it affects classification performance.Inspired by the work of [41-43], the study incorporated dilated convolutional layers as the back end to extract deeper patency information, which is crucial for classification tasks (Fig.3).

    Figure 3:Proposed VGG-CovidNet hybrid architecture.The parameters of the convolutional layers are denoted as Conv-(kernel size), and the second line represents (number of filters) ?(dilation rate)

    4.3 Dilated Convolution

    In the proposed model, one of the crucial components is the dilated convolution inspired by [44].The mathematical formulation of the dilated convolution can be given in Eq.(3).

    O(m,n)is the resultant output of the dilated convolution, whereI(.)and kernelK(i,j)are the inputs to convolution with the length and width ofmandn, respectively.The stridesdetermines the dilation size, where fors=1, the convolution becomes the normal dilation.An alternative to pooling layers, dilated convolution for segmentation tasks has been found to yield significant improvements in accuracy [39-41].Pooling is used to maintain invariance and to avoid overfitting in the learning model.However, it reduces the spatial resolution of the input, and this in turn causes the spatial information of the feature map to be lost.This information loss can be ameliorated through the deconvolutional layers but at the cost of increased complexity and latency, which may not be suitable in every case.Dilated convolution is the solution to avoid the disadvantages of using the pooling and deconvolutional layer.In the dilated convolutional layers, sparse kernels that replace the pooling and deconvolutional layers are used by enlarging the receptive field to avoid information loss and to control complexity.In addition, they reduce the number of operations by performing down-sampling followed by up-sampling.Thus, it facilitates the flexible aggregation of the multiscale contextual information with the same resolution.The dilated convolution demonstration is illustrated in Fig.4, which shows an enlargement of a smallp×pkernel top+(p1)(s1).

    Figure 4:3 × 3 dilated convolution kernel with stride s=1, 2, and 3

    4.4 Network Configuration

    This subsection describes the proposed VGG-CovidNet architecture, which has the same front end, which consists of the first 10 layers of the VGG-16 model.The back end is two branches with different dilation rates.Each branch extracts pieces of spatial information missed by the other processes during dilation.In the end, both outputs are combined to enhance classification accuracy.Moreover, the pooling layers in VGG-16 are reduced to two instead of three layers to overcome the detrimental effects on output accuracy.

    4.5 Implementation Details

    To train the VGG-CovidNet, the procedure of [22] was followed.However, instead of being pretrained on the ImageNet [45] dataset, the VGG-CovidNet was trained afresh with the Adam optimizer The hyperparameter values were adjusted, as indicated below:learning rate =2e?4,number of epochs=500, and batch size=64.Data augmentation was performed through the use of the following augmentation types:intensity shift, horizontal flip, translation, zoom, and rotation.The proposed VGG-CovidNet architecture was implemented with the highly recommended PyTorch [46] framework.The training process for the model was performed on a system with 64 GB memory, a quad-core processor, and a GeForce GTX 1080 Ti graphics card with a large 11 GB frame buffer.The duration of the process was approximately 161 h, and the testing phase lasted 3.5 h.

    4.6 Hyper-Parameter Adjustment

    The hyper-parameters of a neural network are variables that represent the network structure and training methods on the conditions provided.To avoid the fluctuation of the weight values and direct them toward the minimum solution, the learning rate was set to 2e?4.The dropout rate at the last layers was set to 20% because the high value would lead to underfitting, and the low value would have a minimal effect.Based on the network results during training, the number of layers after VGG-16 was set to 6.The activation function is also an important factor in network architecture.Therefore, in the hidden layers, the rectified linear unit activation function was used to avoid negative values during the process, and in the last layer, the softmax entropy was used for multi-class predictions.

    5 Experimental Results

    To evaluate the effectiveness of the proposed architecture, the results were analyzed to better observe the classification performance and to evaluate it against various types of DNN architecture.

    5.1 Network Configuration

    In this section, the effects of the VGG-CovidNet architecture are analyzed in terms of the performance balance in the classification process on the resulting network architecture.For comparison, an in-depth analysis and evaluation of the performance of the following DNN architecture were performed:

    VGG-19[14]:A DNN architecture with low architectural diversity that does not use residual design principles.

    ResNet-50[47]:A DNN architecture that uses lightweight patterns and moderate architectural diversity with minimal long-range connectivity (it does not leverage lightweight PEPX patterns).

    COVID-Net [22]:A deep convolutional multistage neural network design.In the first stage,the design principle was adopted to enable reliable high-performance trainable neural network structures.In the second stage, a design pattern for projection expansion and projection was used to provide a better representation while retaining computational efficiency.

    5.2 Analysis of Results

    To analyze the quantitative performance of the proposed model, the test accuracy, positive predictive value (PPV), and sensitivity were recorded for each infection type on the COVIDx dataset.In addition to these performance measures, the architectural complexity regarding the parameters and computation complexity in terms of the multiplier-accumulator (MAC) numbers were also measured (Tab.1).The values of the number of parameters and MAC of VGG-16 of all 16 layers [14] were 138 M and 15.4 G, respectively, because of the large number and channels of the convolutional layers.However, in the present study, they were reduced to 13.24 M and 14.3 G, respectively, by decreasing the number of convolutions in each layer to two instead of three and by using just 10 layers of VGG-16.One of the main reasons for these lower values was dilated convolution, which reduces the number of parameters and MACs.The accuracy achieved through the proposed model is 96.5%, which is superior to that of all other types of architecture(Tab.2).However, the architectural and computational complexity of the proposed model are slightly higher than that of [22].Unlike the model in [22], this model did not require continuous retraining whenever new data were collected.Tab.2 indicates that the complexity of the proposed model was the result of the higher transfer learning based on VGG-16.Therefore, it will incur higher computational costs than other models during retraining.However, if the first 10 layers of the model are kept frozen, then the computational complexity in retraining the model will become reasonable.

    To gain deeper insights into the proposed model, the sensitivity and PPV were computed(Tabs.3 and 4, respectively).Figs.5 and 6 illustrate the comparison of the PPV and sensitivity of the proposed model with those of other models.These measures, sensitivity and PPV, were derived from the confusion matrix (Fig.7) and compared with the confusion matrix for COVID-Net in Fig.8, which resulted from the testing phase.The proposed model’s sensitivity is 96%, which is the highest value achieved so far for COVID-19 cases.This high sensitivity suggests that the rate of missed COVID-19 positive cases is lower than that in the other state-of-the-art methods (Tab.3).The PPV value is comparable to that of [22], which shows fewer incorrectly classified COVID-19 cases.As shown in Fig.7, two of the four patients were classified as non-COVID-19 patients, and two were classified as normal.In sum, the results indicate that the proposed model outperformed the state-of-the-art method in terms of accuracy and sensitivity on the COVIDx dataset.

    Table 2:Performance comparison of state-of-the-art architecture

    Table 3:Sensitivity for each infection type

    Table 4:Positive predictive value for each infection type

    Figure 5:Comparison of positive predictive values of VGG-CovidNet and those of other models

    Figure 6:Comparison of sensitivity of VGG-CovidNet and that of other models

    Figure 7:Confusion matrix for the proposed model on the COVIDx dataset

    Figure 8:Confusion matrix for the COVID-Net on the COVIDx dataset [22]

    Table 5:Comparison of the proposed model and other models for COVID-19 detection

    The comparison of VGG-COVIDNet and the latest models to detect COVID-19 is shown in Tab.5.The proposed model is superior to the other models because it identifies the features that distinguish COVID-19 from other classes.The same dataset (COVIDx) and ImageNet were used in most of the models.The DeepCOVIDExplainer method employed two ensemble strategies, and ECOVNet used an EfficientNet base model with ImageNet pre-trained weights.Similar to the present work, PDCOVIDNet is a dilated convolution-based COVID-19 detection network.The accuracy of PDCOVIDNet has been found to be comparable to that of the proposed method.However, in PDCOVIDNet, training and testing were performed on a limited dataset.The other ensemble models have limited applicability because of the computational cost of training multiple DL models for ensemble prediction.COVID-Net used many parameters that produced a high computational overhead; therefore, it is not very practicable.

    5.3 Limitations

    The proposed study has a few limitations, such as model complexity and the unavailability of a large set.The minimal number of images in the dataset restricted system performance in terms of the high false-positive rate.The false positive rate can be reduced with the availability of additional COVID-19 images.The proposed model has a higher degree of complexity than the network proposed by [22], as can be observed from the number of parameters shown in Tab.2.There is a need to reduce the complexity without affecting or degrading performance; however,there is a trade-off.Regardless, model retraining will be required with the availability of more images to further increase performance.

    6 Conclusion

    This paper has introduced VGG-CovidNet, a dilated convolution-based bi-branched architecture for classifying infectious images into corresponding classes.The front end of VGG-CovidNet is integrated with the first 10 layers of the VGG-16 model.In these layers, the convolution is reduced to two instead of three to reduce computational complexity.The back end of VGGCovidNet consists of six parallel convolutional layers with different strides of the convolutional kernel.Using the dilated convolution facilitated the capture of the feature maps’spatial information and the reduction of the model’s complexity.A benchmark dataset (COVIDx) was used to train and to test the proposed model.Further, VGG-CovidNet was quantitatively analyzed on the basis of accuracy, sensitivity, and PPV matrices.Comparisons with the state-of-the-art models indicate that VGG-CovidNet undoubtedly yields promising results with 96.5% accuracy, 96.0%sensitivity, and 98% PPV.

    Funding Statement:This work was funded by Institutional Fund projects (Grant No.IFPHI-255-611-2020).The authors therefore gratefully acknowledge the technical and financial support from the Ministry of Education and King Abdulaziz University, DSR, Jeddah, Saudi Arabia.

    Conflicts of Interest:The authors declare that they have no competing interests.

    欧美国产精品一级二级三级| av在线app专区| 成人漫画全彩无遮挡| 免费观看无遮挡的男女| 免费观看无遮挡的男女| 99久国产av精品国产电影| 欧美成人精品欧美一级黄| 老熟女久久久| 国产成人精品一,二区| 巨乳人妻的诱惑在线观看| 精品视频人人做人人爽| 满18在线观看网站| 一级a做视频免费观看| 中文字幕精品免费在线观看视频 | 波野结衣二区三区在线| 亚洲国产av新网站| 一级毛片我不卡| 国产精品久久久久久久久免| 黄色一级大片看看| 亚洲欧洲国产日韩| 午夜福利视频精品| 九九爱精品视频在线观看| 国产一区二区三区综合在线观看 | 韩国精品一区二区三区 | 一本大道久久a久久精品| 免费黄网站久久成人精品| 国产一区二区三区综合在线观看 | 老女人水多毛片| 欧美激情极品国产一区二区三区 | 丝袜喷水一区| 国产福利在线免费观看视频| 日韩大片免费观看网站| 中文字幕精品免费在线观看视频 | 国产精品一区二区在线不卡| 久久精品久久久久久久性| 午夜免费男女啪啪视频观看| 男的添女的下面高潮视频| 亚洲综合色网址| 99久久精品国产国产毛片| 伊人久久国产一区二区| 精品国产一区二区三区久久久樱花| 国产欧美日韩综合在线一区二区| 男男h啪啪无遮挡| 亚洲av综合色区一区| a级片在线免费高清观看视频| 22中文网久久字幕| 搡老乐熟女国产| 少妇 在线观看| 久久精品aⅴ一区二区三区四区 | 亚洲国产av新网站| 国产在视频线精品| 亚洲国产毛片av蜜桃av| 成人毛片60女人毛片免费| 日韩中文字幕视频在线看片| 黄色视频在线播放观看不卡| 18禁国产床啪视频网站| av在线观看视频网站免费| 国产视频首页在线观看| 超碰97精品在线观看| 欧美激情 高清一区二区三区| 午夜日本视频在线| 日韩成人av中文字幕在线观看| 韩国高清视频一区二区三区| 深夜精品福利| 91国产中文字幕| 精品第一国产精品| 成年av动漫网址| 黄色一级大片看看| 亚洲精品乱久久久久久| 熟女人妻精品中文字幕| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 国产高清国产精品国产三级| 晚上一个人看的免费电影| 精品少妇黑人巨大在线播放| 五月伊人婷婷丁香| 最后的刺客免费高清国语| 伊人久久国产一区二区| 天天操日日干夜夜撸| 一本—道久久a久久精品蜜桃钙片| 一本久久精品| 考比视频在线观看| 久久精品久久久久久噜噜老黄| 中文字幕av电影在线播放| 三上悠亚av全集在线观看| 性高湖久久久久久久久免费观看| 国产麻豆69| 99九九在线精品视频| av黄色大香蕉| 日韩一本色道免费dvd| 欧美少妇被猛烈插入视频| 亚洲五月色婷婷综合| 国产精品人妻久久久久久| 只有这里有精品99| 国产综合精华液| 蜜桃国产av成人99| 熟女电影av网| 久久久精品94久久精品| 欧美 日韩 精品 国产| 岛国毛片在线播放| xxx大片免费视频| 女性被躁到高潮视频| 成年女人在线观看亚洲视频| 亚洲av中文av极速乱| 精品第一国产精品| 新久久久久国产一级毛片| 中文字幕另类日韩欧美亚洲嫩草| 久久久久人妻精品一区果冻| 久久人人97超碰香蕉20202| 久久国产精品大桥未久av| 少妇的丰满在线观看| 久久久久精品性色| 午夜免费男女啪啪视频观看| 免费播放大片免费观看视频在线观看| 亚洲精品国产av成人精品| 97人妻天天添夜夜摸| 久久婷婷青草| 欧美bdsm另类| 国语对白做爰xxxⅹ性视频网站| av黄色大香蕉| 亚洲av中文av极速乱| 黄色怎么调成土黄色| 中文字幕制服av| 90打野战视频偷拍视频| 五月玫瑰六月丁香| 超碰97精品在线观看| 精品酒店卫生间| www.色视频.com| 国产精品久久久久成人av| 青春草视频在线免费观看| 久久免费观看电影| 亚洲情色 制服丝袜| 国产精品蜜桃在线观看| 国产精品久久久久久久电影| 国产精品欧美亚洲77777| 免费看av在线观看网站| 熟女电影av网| 国产精品女同一区二区软件| 国产一区有黄有色的免费视频| 在线看a的网站| 精品一区在线观看国产| 国产深夜福利视频在线观看| 一级片免费观看大全| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 大香蕉久久成人网| 精品久久久精品久久久| 国产不卡av网站在线观看| 一区在线观看完整版| 精品卡一卡二卡四卡免费| 国产xxxxx性猛交| 久久久亚洲精品成人影院| 亚洲国产精品一区三区| 中文字幕亚洲精品专区| 国产女主播在线喷水免费视频网站| 亚洲av电影在线观看一区二区三区| 国产成人欧美| 国产免费一级a男人的天堂| 欧美精品av麻豆av| 国产在线一区二区三区精| 永久免费av网站大全| 亚洲欧洲精品一区二区精品久久久 | 精品人妻在线不人妻| 91aial.com中文字幕在线观看| 91精品伊人久久大香线蕉| 欧美精品高潮呻吟av久久| 男女啪啪激烈高潮av片| videossex国产| 涩涩av久久男人的天堂| 久久99蜜桃精品久久| 国产成人a∨麻豆精品| 乱码一卡2卡4卡精品| 亚洲精品乱码久久久久久按摩| 久久青草综合色| 精品亚洲成a人片在线观看| 91午夜精品亚洲一区二区三区| 久久99热这里只频精品6学生| 亚洲国产精品国产精品| 天天操日日干夜夜撸| av又黄又爽大尺度在线免费看| 侵犯人妻中文字幕一二三四区| 精品人妻熟女毛片av久久网站| 久久久久久伊人网av| 久久久久精品性色| 欧美日韩综合久久久久久| av一本久久久久| 精品99又大又爽又粗少妇毛片| 一级毛片黄色毛片免费观看视频| 人人澡人人妻人| 国产乱人偷精品视频| 国产精品 国内视频| 久久精品久久精品一区二区三区| 国产淫语在线视频| 99香蕉大伊视频| 国产亚洲精品久久久com| 国产成人午夜福利电影在线观看| 久久这里只有精品19| 一二三四在线观看免费中文在 | 男女高潮啪啪啪动态图| 精品人妻一区二区三区麻豆| 男人操女人黄网站| 大香蕉97超碰在线| 岛国毛片在线播放| 国产精品三级大全| 乱码一卡2卡4卡精品| 免费高清在线观看视频在线观看| 咕卡用的链子| 亚洲精品美女久久久久99蜜臀 | 18禁国产床啪视频网站| 亚洲综合色网址| 欧美 日韩 精品 国产| 国产福利在线免费观看视频| 午夜91福利影院| 少妇被粗大的猛进出69影院 | 亚洲成人一二三区av| 熟女电影av网| 午夜福利视频精品| 亚洲,欧美,日韩| 中国三级夫妇交换| 国产精品蜜桃在线观看| 久久99精品国语久久久| 国产精品一区www在线观看| 久久久国产精品麻豆| 一级a做视频免费观看| 熟妇人妻不卡中文字幕| 久久99精品国语久久久| 夜夜骑夜夜射夜夜干| 亚洲成av片中文字幕在线观看 | 色视频在线一区二区三区| 在线观看人妻少妇| 少妇被粗大猛烈的视频| 日韩一区二区三区影片| 高清黄色对白视频在线免费看| 日本av手机在线免费观看| 一级爰片在线观看| 少妇人妻 视频| 国产成人aa在线观看| 国产男女超爽视频在线观看| 久久97久久精品| 欧美变态另类bdsm刘玥| 777米奇影视久久| 国产成人欧美| 国产黄频视频在线观看| 国产老妇伦熟女老妇高清| 免费在线观看完整版高清| 最近最新中文字幕免费大全7| a级毛色黄片| 国产精品久久久久久久电影| 大片免费播放器 马上看| 日韩成人伦理影院| 有码 亚洲区| 亚洲图色成人| 一区在线观看完整版| 国产 精品1| 免费av不卡在线播放| 欧美精品av麻豆av| 秋霞伦理黄片| 国产欧美另类精品又又久久亚洲欧美| 黄网站色视频无遮挡免费观看| 国产在线免费精品| 久久99蜜桃精品久久| 捣出白浆h1v1| 边亲边吃奶的免费视频| 成人漫画全彩无遮挡| 丝袜美足系列| 建设人人有责人人尽责人人享有的| 色5月婷婷丁香| www.色视频.com| 下体分泌物呈黄色| 欧美成人午夜精品| 黄色配什么色好看| 一本—道久久a久久精品蜜桃钙片| 国产精品一区www在线观看| 熟女av电影| 国产一区二区在线观看日韩| 黄色毛片三级朝国网站| 亚洲图色成人| 人成视频在线观看免费观看| 视频中文字幕在线观看| 人妻系列 视频| 午夜免费男女啪啪视频观看| 亚洲av电影在线进入| 日本wwww免费看| 咕卡用的链子| 午夜福利视频在线观看免费| 91午夜精品亚洲一区二区三区| 黄色怎么调成土黄色| av黄色大香蕉| 国产黄频视频在线观看| √禁漫天堂资源中文www| 亚洲成人一二三区av| 国产av码专区亚洲av| 在线观看国产h片| 男人舔女人的私密视频| 久久国产精品大桥未久av| 午夜激情久久久久久久| 最近的中文字幕免费完整| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产高清三级在线| 国产免费福利视频在线观看| 久久久精品94久久精品| 久久人人97超碰香蕉20202| 男的添女的下面高潮视频| 黄色怎么调成土黄色| av黄色大香蕉| 亚洲国产精品一区二区三区在线| 97在线视频观看| 丰满饥渴人妻一区二区三| 黄色配什么色好看| 日韩三级伦理在线观看| 黄色视频在线播放观看不卡| 国产av精品麻豆| 国产精品欧美亚洲77777| 亚洲精品456在线播放app| 久久女婷五月综合色啪小说| 18禁裸乳无遮挡动漫免费视频| 国产 精品1| 欧美激情国产日韩精品一区| av在线观看视频网站免费| 好男人视频免费观看在线| 亚洲三级黄色毛片| 黄片无遮挡物在线观看| 天堂8中文在线网| 日韩视频在线欧美| 老司机影院毛片| 中文精品一卡2卡3卡4更新| 欧美 亚洲 国产 日韩一| a级片在线免费高清观看视频| 99香蕉大伊视频| 久久久久精品性色| 欧美+日韩+精品| 日日撸夜夜添| 亚洲综合精品二区| 久久精品aⅴ一区二区三区四区 | 最近中文字幕2019免费版| 国产视频首页在线观看| 久久久a久久爽久久v久久| 亚洲av中文av极速乱| 亚洲高清免费不卡视频| 久热这里只有精品99| 亚洲天堂av无毛| videossex国产| 18禁国产床啪视频网站| 少妇的丰满在线观看| 久久99热6这里只有精品| 国产在线免费精品| 在线观看国产h片| 精品99又大又爽又粗少妇毛片| 狠狠婷婷综合久久久久久88av| 午夜久久久在线观看| 十八禁高潮呻吟视频| freevideosex欧美| 天堂8中文在线网| 亚洲伊人色综图| 美女脱内裤让男人舔精品视频| 欧美国产精品一级二级三级| 麻豆乱淫一区二区| 亚洲成人av在线免费| 国产男女内射视频| 在线 av 中文字幕| 日本欧美视频一区| 一区二区三区精品91| 日本爱情动作片www.在线观看| 免费黄网站久久成人精品| 国产成人精品久久久久久| av卡一久久| 国产色婷婷99| 国产黄频视频在线观看| 春色校园在线视频观看| 久久久国产精品麻豆| 中国三级夫妇交换| 国产精品蜜桃在线观看| 久久99一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 免费黄网站久久成人精品| 久久99一区二区三区| 欧美精品av麻豆av| 久久精品国产亚洲av天美| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 国产在线一区二区三区精| 只有这里有精品99| 五月天丁香电影| av有码第一页| 亚洲内射少妇av| 黑人巨大精品欧美一区二区蜜桃 | 一二三四中文在线观看免费高清| 极品少妇高潮喷水抽搐| 亚洲av国产av综合av卡| 亚洲欧美清纯卡通| 亚洲经典国产精华液单| 久久ye,这里只有精品| 精品亚洲成a人片在线观看| 91国产中文字幕| 欧美人与性动交α欧美软件 | 天美传媒精品一区二区| 精品亚洲成a人片在线观看| 日本av免费视频播放| 欧美日韩av久久| 黑丝袜美女国产一区| 色视频在线一区二区三区| 中文字幕制服av| 日韩一本色道免费dvd| 精品卡一卡二卡四卡免费| 亚洲性久久影院| 亚洲综合色惰| 美女国产高潮福利片在线看| 午夜免费观看性视频| 日本vs欧美在线观看视频| 秋霞在线观看毛片| 亚洲欧美一区二区三区黑人 | 久久精品国产综合久久久 | 丰满乱子伦码专区| 日韩视频在线欧美| 久久精品国产自在天天线| 国产 精品1| 少妇高潮的动态图| 国产视频首页在线观看| 在线观看一区二区三区激情| 中文字幕人妻熟女乱码| 国产无遮挡羞羞视频在线观看| 国产福利在线免费观看视频| 久久韩国三级中文字幕| 日韩成人伦理影院| 日韩大片免费观看网站| 黄色 视频免费看| 在线观看免费日韩欧美大片| 久久午夜福利片| 免费在线观看完整版高清| 国产精品欧美亚洲77777| 成人综合一区亚洲| 人人妻人人澡人人看| 免费看av在线观看网站| 国产女主播在线喷水免费视频网站| 久久国产亚洲av麻豆专区| 宅男免费午夜| 亚洲精品色激情综合| 美女视频免费永久观看网站| 亚洲一码二码三码区别大吗| 美女福利国产在线| 18+在线观看网站| 亚洲av福利一区| 中文字幕制服av| 只有这里有精品99| 欧美xxxx性猛交bbbb| 香蕉丝袜av| 国产精品一区二区在线不卡| 欧美变态另类bdsm刘玥| 精品久久国产蜜桃| 99久久人妻综合| 日本色播在线视频| 久久久久国产精品人妻一区二区| 精品国产露脸久久av麻豆| 9191精品国产免费久久| 深夜精品福利| 精品久久国产蜜桃| 国产一区有黄有色的免费视频| 亚洲美女视频黄频| 男人舔女人的私密视频| 视频区图区小说| 女的被弄到高潮叫床怎么办| 亚洲精品,欧美精品| 国产无遮挡羞羞视频在线观看| 汤姆久久久久久久影院中文字幕| 九色亚洲精品在线播放| 亚洲一码二码三码区别大吗| 日本-黄色视频高清免费观看| 亚洲图色成人| 蜜臀久久99精品久久宅男| 黄色视频在线播放观看不卡| 你懂的网址亚洲精品在线观看| 亚洲精品av麻豆狂野| 久久鲁丝午夜福利片| 国产成人av激情在线播放| 国产不卡av网站在线观看| 国国产精品蜜臀av免费| av线在线观看网站| 亚洲国产精品国产精品| 国产精品成人在线| 久久人人爽人人片av| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 波野结衣二区三区在线| 美女福利国产在线| 成人免费观看视频高清| 美女福利国产在线| 国产免费现黄频在线看| 日本av手机在线免费观看| 人妻 亚洲 视频| 亚洲国产av新网站| 国产在线一区二区三区精| 精品久久久久久电影网| 18禁动态无遮挡网站| 69精品国产乱码久久久| 在线观看www视频免费| 中文欧美无线码| 国产精品一二三区在线看| 亚洲精品乱码久久久久久按摩| 97精品久久久久久久久久精品| 久久99一区二区三区| 日韩免费高清中文字幕av| 欧美日韩综合久久久久久| 亚洲成人手机| 成年人免费黄色播放视频| 卡戴珊不雅视频在线播放| 免费看光身美女| 色视频在线一区二区三区| 免费在线观看完整版高清| 久久综合国产亚洲精品| 十八禁网站网址无遮挡| av在线播放精品| 国产在线免费精品| 黄色怎么调成土黄色| 久久鲁丝午夜福利片| 亚洲成人手机| 黑丝袜美女国产一区| 亚洲欧美日韩卡通动漫| 90打野战视频偷拍视频| 中文字幕精品免费在线观看视频 | av一本久久久久| 日韩视频在线欧美| av免费在线看不卡| 国产亚洲av片在线观看秒播厂| 日韩中文字幕视频在线看片| av在线老鸭窝| 日韩中文字幕视频在线看片| 久久国产精品男人的天堂亚洲 | 在线观看免费视频网站a站| 99久久中文字幕三级久久日本| 国产爽快片一区二区三区| 亚洲欧洲国产日韩| 晚上一个人看的免费电影| 国产一区二区三区av在线| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 三级国产精品片| 波多野结衣一区麻豆| 天天躁夜夜躁狠狠躁躁| 国产一区二区在线观看av| av线在线观看网站| 国产一区二区在线观看日韩| 亚洲图色成人| 一本—道久久a久久精品蜜桃钙片| 亚洲国产毛片av蜜桃av| 亚洲成国产人片在线观看| 亚洲国产看品久久| 欧美另类一区| 18+在线观看网站| 综合色丁香网| 国产精品一区二区在线观看99| 亚洲成人手机| 色网站视频免费| 国产亚洲精品久久久com| 一级片免费观看大全| 黄色视频在线播放观看不卡| 日韩人妻精品一区2区三区| 黄网站色视频无遮挡免费观看| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 日产精品乱码卡一卡2卡三| 丰满饥渴人妻一区二区三| 国产精品秋霞免费鲁丝片| 巨乳人妻的诱惑在线观看| 国产精品三级大全| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 精品国产露脸久久av麻豆| 亚洲av免费高清在线观看| 亚洲av中文av极速乱| 香蕉精品网在线| 中文字幕人妻熟女乱码| 黑人巨大精品欧美一区二区蜜桃 | 日韩电影二区| 久久久国产一区二区| 国产av一区二区精品久久| 一级a做视频免费观看| kizo精华| 人人妻人人爽人人添夜夜欢视频| 91在线精品国自产拍蜜月| 99国产综合亚洲精品| av电影中文网址| 国产成人精品久久久久久| 亚洲国产精品国产精品| 国产精品成人在线| 黄色 视频免费看| 一区二区av电影网| 精品99又大又爽又粗少妇毛片| 国产精品国产三级国产av玫瑰| 国产淫语在线视频| 国产高清国产精品国产三级| 高清视频免费观看一区二区| 精品卡一卡二卡四卡免费| 中文天堂在线官网| 欧美日韩亚洲高清精品| 久久久久视频综合| 在线观看免费视频网站a站| 两个人看的免费小视频| 国产在线视频一区二区| 亚洲高清免费不卡视频| 亚洲av日韩在线播放| 欧美成人午夜精品| 一区二区av电影网| 久久国产精品男人的天堂亚洲 | 久久国产精品男人的天堂亚洲 | 欧美老熟妇乱子伦牲交| 欧美日韩国产mv在线观看视频| 成人毛片a级毛片在线播放| 一边摸一边做爽爽视频免费| 少妇被粗大猛烈的视频| 涩涩av久久男人的天堂| 国产精品欧美亚洲77777| 国产一区二区三区综合在线观看 | 日韩,欧美,国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 亚洲一码二码三码区别大吗|