一种鸢乌贼渔情数据处理方法、系统、装置和存储介质与流程

文档序号:18943360发布日期:2019-10-23 01:20阅读:294来源:国知局
本发明涉及渔情数据分析
技术领域
:,尤其是一种鸢乌贼渔情数据处理方法、系统、装置和存储介质。
背景技术
::鸢乌贼作为南海最重要的中上层渔业资源捕捞对象。据中国水产科学研究院南海水产研究所调查评估显示,南海鸢乌贼资源年可捕量在400万吨以上,而我国年捕捞量仅10万吨,明显处于未开发状态。为了指导鸢乌贼捕捞生产,需要对鸢乌贼渔情进行分析。目前对南海外海鸢乌贼渔情预报系统的研究,只处于资源量的评估及模型研究阶段,这依赖于使用渔船等设备进行实地考察,分析结果的实时性、精度和可靠性均得不到保证。技术实现要素:为了解决上述技术问题,本发明的目在于提供一种鸢乌贼渔情数据处理方法、系统、装置和存储介质。一方面,本发明包括一种鸢乌贼渔情数据处理方法,包括以下步骤:从渔捞日志中获取产量数据;从数据库中获取环境数据;所述环境数据包括海表面温度、海平面高度、海表面叶绿素a浓度和海表面盐度;对所述产量数据和环境数据进行空间尺度统一化处理;使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化;根据经过标准化后的所述产量数据计算单因子栖息地指数;根据所述单因子栖息地指数计算综合栖息地指数。进一步地,所述对所述产量数据和环境数据进行空间尺度统一化处理这一步骤,具体包括:使用公式对所述产量数据进行空间分辨率转换,式中cpuei,l,k,j为在第i年l月在k经度和j纬度的每艘船每天捕捞产量,∑catchi,l,k,j为在第i年l月在k经度和j纬度的总产量,∑ei,l,k,j为在第i年l月在k经度和j纬度的总作业天数;使用公式对所述海表面温度进行空间分辨率转换,式中ssti,l,k,j为在第i年l月在k经度和j纬度的渔区内平均海表面温度,sstx为在第i年l月在k经度和j纬度的渔区内单个海表面温度,n为在第i年l月在k经度和j纬度的渔区内海表面温度测量点个数;使用公式对所述海表面叶绿素a浓度进行空间分辨率转换,式中chlai,l,k,j为在第i年l月在k经度和j纬度的渔区内平均海表面叶绿素a浓度,chlax为在第i年l月在k经度和j纬度的渔区内单个海表面叶绿素a浓度,n为在第i年l月在k经度和j纬度的渔区内海表面叶绿素a浓度测量点个数;使用三次线条插值方法对所述海平面高度和所述海表面盐度进行空间分辨率转换。进一步地,所述使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化这一步骤,所使用的公式为:ln(cpuei,l,k,j+1)=k+a1yeari+a2monthl+a3lonk+a4latj+a5sst+a6chla+a6sss+a7ssh+εi,l,k,j或ln(cpuei,l,k,j+1)=s(yeari)+s(monthl)+s(lonk)+s(latj)+s(sst)+s(chla)+s(sss)+s(ssh)+εi,l,k,j;式中,cpuei,l,k,j为在第i年l月在k经度和j纬度的每艘船每天捕捞产量,εi,l,k,j为与cpuei,l,k,j对应的误差项,yeari为第i年,monthl为第l月,lonk为k经度,latj为j纬度,a1、a2、a3、a4、a5、a6、a7和k分别为自定义的参数,s()为样条平滑函数。进一步地,所述根据经过标准化后的所述产量数据计算单因子栖息地指数这一步骤,所使用的公式为:式中,cpuei为经过标准化后的一个每艘船每天捕捞产量,cpuemax为经过标准化后的所有每艘船每天捕捞产量中的最大值,sii为计算得到的单因子栖息地指数。进一步地,所述根据所述单因子栖息地指数计算综合栖息地指数这一步骤,所使用的公式为:hsi=minsii,hsi=maxsii,或式中,hsi为计算得到的综合栖息地指数,sii为单因子栖息地指数,n为单因子栖息地指数的个数。进一步地,所述一种鸢乌贼渔情数据处理方法还包括以下步骤:将所述产量数据和环境数据按照所属的渔区进行归类;对归属同一渔区的产量数据和环境数据建立对应关系。进一步地,所述一种鸢乌贼渔情数据处理方法还包括以下步骤:分别计算各渔区对应的综合栖息地指数;根据各所述综合栖息地指数对各渔区进行排名;将排名最前的至少一个渔区标记为中心渔区;输出与所述中心渔区相关的信息。另一方面,本发明实施例还包括一种鸢乌贼渔情数据处理系统,包括:产量数据模块,用于从渔捞日志中获取产量数据;环境数据模块,用于从数据库中获取环境数据;所述环境数据包括海表面温度、海平面高度、海表面叶绿素a浓度和海表面盐度;尺度统一模块,用于对所述产量数据和环境数据进行空间尺度统一化处理;标准化模块,用于使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化;单因子栖息地指数计算模块,用于根据经过标准化后的所述产量数据计算单因子栖息地指数;综合栖息地指数计算模块,用于根据所述单因子栖息地指数计算综合栖息地指数。另一方面,本发明实施例还包括一种鸢乌贼渔情数据处理装置,包括存储器和处理器,所述存储器用于存储至少一个程序,所述处理器用于加载所述至少一个程序以执行本发明方法。另一方面,本发明实施例还包括一种存储介质,其中存储有处理器可执行的指令,所述处理器可执行的指令在由处理器执行时用于执行本发明方法。本发明的有益效果是:通过对积累的产量数据和环境数据进行处理分析,所计算出的综合栖息地指数能够及时反映鸢乌贼的捕捞价值,从而向渔船推送新的渔区,缩短渔船寻找渔场的时间、减少耗费的成本,提高鸢乌贼捕捞产量和捕捞效率。附图说明图1为本发明实施例中鸢乌贼渔情数据处理方法的流程图;图2为本发明实施例中鸢乌贼渔情数据处理系统的结构框图。具体实施方式本实施例包括一种鸢乌贼渔情数据处理方法,参照图1,包括以下步骤:s1.从渔捞日志中获取产量数据;s2.从数据库中获取环境数据;所述环境数据包括海表面温度、海平面高度、海表面叶绿素a浓度和海表面盐度;s3.对所述产量数据和环境数据进行空间尺度统一化处理;s4.使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化;s5.根据经过标准化后的所述产量数据计算单因子栖息地指数;s6.根据所述单因子栖息地指数计算综合栖息地指数。步骤s1中,可以从水产技术推广站等地获取其管理的渔船的渔捞日志,所述渔捞日志通常记录了各渔船的船名、船主、主机功率、作业日期、作业时间、作业经纬度位置和渔获量等数据,这些数据的时间分辨率为天,空间分辨率为0.5°×0.5°,即每0.5°×0.5°所划定的空间作为一个渔区,每个渔区都分别对应位于其内部的渔船的船名、船主、主机功率、作业日期、作业时间、作业经纬度位置和渔获量等数据。例如经度110°-110.5°e,纬度10.5°-11°n确定了一个0.5°×0.5°的范围,那么就可以将110°-110.5°e,10.5°-11°n所确定的范围标记为一个渔区,用(110.25°e,10.75°n)表示。本实施例中,所述产量数据是指渔捞日志中的渔获量数据。步骤s2中,从美国nasa的modisaqua卫星传感器数据库(https://oceandata.sci.gsfc.nasa.gov/)获取海表面温度(seasurfacetemperature,sst)和海表面叶绿素a浓度(seasurfacechlorophyllaconcentration,chla)等环境数据,海表面温度和海表面叶绿素a浓度的时间分辨率为月,空间分辨率为0.0417°×0.0417°。从哥白尼海洋环境管理服务数据库(cmems,copernicusmarineenvironmentmanagementservice)(http://marine.copernicus.eu)中的全球海洋物理再分析产品(globaloceanphysicalreanalysisproduct)中获取海平面高度(seasurfaceheight,ssh)和海表面盐度(seasurfacesalinity,sss)等环境数据,海平面高度和海表面盐度的时间分辨率为月,空间分辨率为1°×1°。步骤s1和步骤s2中所获取的产量数据和环境数据具有不同的空间分辨率,这不利于后续的分析过程。步骤s3中,对产量数据和环境数据进行空间尺度统一化处理,使得它们具有相同的空间分辨率。步骤s4中,使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化,从而在步骤s5和s6中进一步计算出单因子栖息地指数和综合栖息地指数。步骤s4-s5使用一元非线性回归来建立各个显著性因子(即具有较大影响的产量数据和环境数据)与单因子栖息地指数之间的关系模型,从而计算出单因子栖息地指数。最后计算出的综合栖息地指数可以反映渔区环境对鸢乌贼生活的适宜程度,即综合栖息地指数越大,渔区环境越适宜鸢乌贼生活繁殖,表明该渔区具有捕捞价值的概率越大,据此指导渔船进行鸢乌贼捕捞工作。进一步作为优选的实施方式,所述步骤s3,即对所述产量数据和环境数据进行空间尺度统一化处理这一步骤,具体包括:s301.使用公式对所述产量数据进行空间分辨率转换,式中cpuei,l,k,j为在第i年l月在k经度和j纬度的每艘船每天捕捞产量,∑catchi,l,k,j为在第i年l月在k经度和j纬度的总产量,∑ei,l,k,j为在第i年l月在k经度和j纬度的总作业天数;s302.使用公式对所述海表面温度进行空间分辨率转换,式中ssti,l,k,j为在第i年l月在k经度和j纬度的渔区内平均海表面温度,sstx为在第i年l月在k经度和j纬度的渔区内单个海表面温度,n为在第i年l月在k经度和j纬度的渔区内海表面温度测量点个数;s303.使用公式对所述海表面叶绿素a浓度进行空间分辨率转换,式中chlai,l,k,j为在第i年l月在k经度和j纬度的渔区内平均海表面叶绿素a浓度,chlax为在第i年l月在k经度和j纬度的渔区内单个海表面叶绿素a浓度,n为在第i年l月在k经度和j纬度的渔区内海表面叶绿素a浓度测量点个数;s304.使用三次线条插值方法对所述海平面高度和所述海表面盐度进行空间分辨率转换。步骤s301中,∑catchi,l,k,j和∑ei,l,k,j均为从捕捞日志中获取到的产量数据。通过执行步骤s301,所得到的cpuei,l,k,j的空间分辨率被统一为0.5°×0.5°。步骤s302和s303中,sstx和chlax均为从数据库中获取到的环境数据,并且它们对应的经纬度划分到相应的渔区。通过执行步骤s302和s303,海表面温度ssti,l,k,j和海表面叶绿素a浓度chlai,l,k,j的空间分辨率被统一为0.5°×0.5°。步骤s304中,所述海平面高度ssh和所述海表面盐度sss均为从数据库中获取到的环境数据,并且它们对应的经纬度划分到相应的渔区。通过执行步骤s302和s303,海平面高度ssh和所述海表面盐度sss的空间分辨率被统一为0.5°×0.5°。使用计算机程序执行步骤s304所使用的三次线条插值方法时,可以运行以下代码:clear;clc;file='文件位置';%数据读入data=xlsread(file);%赋值lat=data(:,1);%其中的1代表第1列lon=data(:,2);%chl=data(:,3);%其中的3代表第3列datanum=length(chl);%数据的长度monnum=3;%月份的数量[latuni,latindex]=unique(lat);%纬度唯一值[lonuni,lonindex]=unique(lon);%经度唯一值lonuni=sort(lonuni);%经度升序latuni=sort(latuni);%纬度降序,'descend'm1=length(latuni);%行数是433n1=length(lonuni);%列数是529chln1=reshape(chl,[datanum/monnum,monnum]);%重塑矩阵chlapril=reshape(chln1(:,3),[n1,m1]);[m2,n2]=size(chlapril);%size(行数,列数)[latn1,lonn1]=meshgrid(latuni,lonuni);%lonn是列数,latn是行数latunit=min(latuni):0.5:max(latuni);latunit=latunit-0.25;lonunit=min(lonuni):0.5:max(lonuni);lonunit=lonunit-0.25;[latn2,lonn2]=meshgrid(latunit,lonunit);%lonn是列数,latn是行数chlinterp=interp2(latn1,lonn1,chlapril,latn2,lonn2,'spline');%三次样条插值[m3,n3]=size(chlinterp);chlwrite=reshape(chlinterp,[m3*n31]);latwrite=reshape(latn2,[m3*n31]);lonwrite=reshape(lonn2,[m3*n31]);variablename={'lat','lon','wind'};t=table(latwrite,lonwrite,chlwrite);t=table2cell(t);writefile=('c:\users\xps13\desktop\');%文件读出xlswrite(writefile,[variablename;t]);进一步作为优选的实施方式,所述步骤s4,即使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化这一步骤,所使用的公式为:ln(cpuei,l,k,j+1)=k+a1yeari+a2monthl+a3lonk+a4latj+a5sst+a6chla+a6sss+a7ssh+εi,l,k,j或ln(cpuei,l,k,j+1)=s(yeari)+s(monthl)+s(lonk)+s(latj)+s(sst)+s(chla)+s(sss)+s(ssh)+εi,l,k,j;式中,cpuei,l,k,j为在第i年l月在k经度和j纬度的每艘船每天捕捞产量,εi,l,k,j为与cpuei,l,k,j对应的误差项,yeari为第i年,monthl为第l月,lonk为k经度,latj为j纬度,a1、a2、a3、a4、a5、a6、a7和k分别为自定义的参数,s()为样条平滑函数。步骤s4中所使用的glm模型假设响应变量的期望值与解释变量呈线性关系,例如:式中g为关联函数,μi=e(yi),xi为第i个响应变量的解释变量,β为模型估计参数,yi为第i个响应变量。假设cpue服从对数正态分布,因此glm模型表示为公式ln(cpuei,l,k,j+1)=k+a1yeari+a2monthl+a3lonk+a4latj+a5sst+a6chla+a6sss+a7ssh+εi,l,k,j,式中,cpuei,l,k,j为在第i年l月在k经度和j纬度的每艘船每天捕捞产量,εi,l,k,j为与cpuei,l,k,j对应的误差项,yeari为第i年,monthl为第l月,lonk为k经度,latj为j纬度,a1、a2、a3、a4、a5、a6、a7和k分别为自定义的参数。glm模型中,将时间(年、月)、空间(经度、纬度)、环境(sst、chla、ssh、salinity)因素作为解释变量,其中变量年、月、经度、纬度为分类离散变量,其他变量为连续变量。cpue加上常数1,再作对数变换后,作为响应变量,以解决cpue为0的情况。步骤s4中所使用的gam模型为glm模型的延伸,可以用来表示响应变量和解释变量的非线性关系,即:式中,g为链接函数,μi=e(yi),xi为第i个响应变量的解释变量,yi为第i个响应变量。α为适合函数的截距;ε表示残差;fi(xi)表示各自变量的任意单变量函数,。因此gam模型表示为公式ln(cpuei,l,k,j+1)=s(yeari)+s(monthl)+s(lonk)+s(latj)+s(sst)+s(chla)+s(sss)+s(ssh)+εi,l,k,j,式中s()为样条平滑函数(splinesmoothing)。将解释变量依次加入gam模型,得到包含不同个数解释变量的gam模型。aic值最小的模型,为最佳模。aic值的计算如下:aic=-2lnl(p1,…pm,σ2)+2m式中,m为模型中参数的个数。本实施例中,利用r软件rodbc和mgcv包进行统计分析。利用逐步回归(stepwisemethod)选择对模型影响显著的变量,确定gam的表达形式。进一步作为优选的实施方式,所述步骤s5,即根据经过标准化后的所述产量数据计算单因子栖息地指数这一步骤,所使用的公式为:式中,cpuei为经过步骤s4标准化后的一个每艘船每天捕捞产量,cpuemax为经过步骤s4标准化后的所有每艘船每天捕捞产量中的最大值,sii为计算得到的单因子栖息地指数。进一步作为优选的实施方式,所述步骤s6,所述根据所述单因子栖息地指数计算综合栖息地指数这一步骤,所使用的公式为:hsi=minsii,hsi=maxsii,或式中,hsi为计算得到的综合栖息地指数,sii为单因子栖息地指数,n为单因子栖息地指数的个数。本实施例中,使用公式进行计算的方法被称为连乘法,使用公式hsi=minsii进行计算的方法被称为最小值法,使用公式hsi=maxsii进行计算的方法被称为最大值法,使用公式进行计算的方法被称为算术平均法,使用公式进行计算的方法被称为几何平均法。通过连乘法、最小值法、最大值法、算术平均法或几何平均法,均可根据多个单因子栖息地指数计算出一个综合栖息地指数,该综合栖息地指数综合地考虑了所有单因子栖息地指数的作用。进一步作为优选的实施方式,所述鸢乌贼渔情数据处理方法还包括以下步骤:s2a.将所述产量数据和环境数据按照所属的渔区进行归类;s2b.对归属同一渔区的产量数据和环境数据建立对应关系。本实施例中,所述步骤s2a和s2b在步骤s2之后、s3之前执行。从渔捞日志中获取得到的产量数据和从数据库中获取得到环境数据的都对应有经纬度信息,可以确定产量数据和环境数据各自对应的渔区,从而对产量数据和环境数据进行归类。将归属同一渔区的产量数据和环境数据建立对应关系,具体地,可以使用excel2010软件中的vlookup执行,使得通过经纬度可以查找到相应的产量数据和环境数据,便于对产量数据和环境数据进行管理。进一步作为优选的实施方式,所述鸢乌贼渔情数据处理方法还包括以下步骤:s7.分别计算各渔区对应的综合栖息地指数;s8.根据各所述综合栖息地指数对各渔区进行排名;s9.将排名最前的至少一个渔区标记为中心渔区;s10.输出与所述中心渔区相关的信息。步骤s7中,针对各渔区分别执行步骤s1-s6,从而计算出各渔区对应的综合栖息地指数。步骤s8中,按照对应的综合栖息地指数从大到小的顺序,对各渔区进行排名。由于综合栖息地指数越大表示相应渔区越适合鸢乌贼生存,该渔区对鸢乌贼捕捞产生的经济价值越大,因此排名越靠前的渔区越可能成为渔船的目标渔区。步骤s9中,可以选择将排名最前的一个、三个或五个等数量的渔区标记为中心渔区,从而提示这些渔区具有较高的经济价值。步骤s10中,通过卫星电话或者4g等通信手段,将中心渔区的经纬度等相关信息发送到渔船,从而指导渔船到这些中心渔区进行捕捞生产。综上所述,本实施例中的鸢乌贼渔情数据处理方法具有以下优点:通过对积累的产量数据和环境数据进行处理分析,所计算出的综合栖息地指数能够及时反映鸢乌贼的捕捞价值,从而向渔船推送新的渔区,缩短渔船寻找渔场的时间、减少耗费的成本,提高鸢乌贼捕捞产量和捕捞效率。本实施例还包括一种鸢乌贼渔情数据处理系统,参照图2,包括:产量数据模块,用于从渔捞日志中获取产量数据;环境数据模块,用于从数据库中获取环境数据;所述环境数据包括海表面温度、海平面高度、海表面叶绿素a浓度和海表面盐度;尺度统一模块,用于对所述产量数据和环境数据进行空间尺度统一化处理;标准化模块,用于使用glm模型或gam模型对所述环境数据进行筛选,根据筛选结果对所述产量数据进行标准化;单因子栖息地指数计算模块,用于根据经过标准化后的所述产量数据计算单因子栖息地指数;综合栖息地指数计算模块,用于根据所述单因子栖息地指数计算综合栖息地指数。所述产量数据模块、环境数据模块、尺度统一模块、标准化模块、单因子栖息地指数计算模块和综合栖息地指数计算模块可以是计算机系统上具有相应功能的硬件模块或软件模块。本实施例还包括一种鸢乌贼渔情数据处理装置,包括存储器和处理器,所述存储器用于存储至少一个程序,所述处理器用于加载所述至少一个程序以执行本发明鸢乌贼渔情数据处理方法。本实施例还包括一种存储介质,其中存储有处理器可执行的指令,所述处理器可执行的指令在由处理器执行时用于执行本发明鸢乌贼渔情数据处理方法。本实施例中的鸢乌贼渔情数据处理系统、装置和存储介质,可以执行本发明的鸢乌贼渔情数据处理方法,可执行方法实施例的任意组合实施步骤,具备该方法相应的功能和有益效果。需要说明的是,如无特殊说明,当某一特征被称为“固定”、“连接”在另一个特征,它可以直接固定、连接在另一个特征上,也可以间接地固定、连接在另一个特征上。此外,本公开中所使用的上、下、左、右等描述仅仅是相对于附图中本公开各组成部分的相互位置关系来说的。在本公开中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式,除非上下文清楚地表示其他含义。此外,除非另有定义,本实施例所使用的所有的技术和科学术语与本
技术领域
:的技术人员通常理解的含义相同。本实施例说明书中所使用的术语只是为了描述具体的实施例,而不是为了限制本发明。本实施例所使用的术语“和/或”包括一个或多个相关的所列项目的任意的组合。应当理解,尽管在本公开可能采用术语第一、第二、第三等来描述各种元件,但这些元件不应限于这些术语。这些术语仅用来将同一类型的元件彼此区分开。例如,在不脱离本公开范围的情况下,第一元件也可以被称为第二元件,类似地,第二元件也可以被称为第一元件。本实施例所提供的任何以及所有实例或示例性语言(“例如”、“如”等)的使用仅意图更好地说明本发明的实施例,并且除非另外要求,否则不会对本发明的范围施加限制。应当认识到,本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现,其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而,若需要,该程序可以以汇编或机器语言实现。在任何情况下,该语言可以是编译或解释的语言。此外,为此目的该程序能够在编程的专用集成电路上运行。此外,可按任何合适的顺序来执行本实施例描述的过程的操作,除非本实施例另外指示或以其他方式明显地与上下文矛盾。本实施例描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机系统的控制下执行,并且可作为共同地在一个或多个处理器上执行的代码(例如,可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。进一步,所述方法可以在可操作地连接至合适的任何类型的计算平台中实现,包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本发明的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现,无论是可移动的还是集成至计算平台,如硬盘、光学读取和/或写入存储介质、ram、rom等,使得其可由可编程计算机读取,当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外,机器可读代码,或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时,本实施例所述的发明包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本发明所述的方法和技术编程时,本发明还包括计算机本身。计算机程序能够应用于输入数据以执行本实施例所述的功能,从而转换输入数据以生成存储至非易失性存储器的输出数据。输出信息还可以应用于一个或多个输出设备如显示器。在本发明优选的实施例中,转换的数据表示物理和有形的对象,包括显示器上产生的物理和有形对象的特定视觉描绘。以上所述,只是本发明的较佳实施例而已,本发明并不局限于上述实施方式,只要其以相同的手段达到本发明的技术效果,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。在本发明的保护范围内其技术方案和/或实施方式可以有各种不同的修改和变化。当前第1页12当前第1页12
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1