一种基于变分模态分解和多尺度下采样的光伏功率出力预测方法与流程

1.本发明属于电网运行和调度领域，涉及了一种基于变分模态分解和多尺度下采样的光伏功率出力预测方法。


背景技术：

2.天气、地理位置等外界因素决定了光电具有随机性、间歇性。当光伏发电接入电网后，这种不稳定的特性会给整个电网带来安全隐患、影响电能质量。预测光伏发电功率，提前做好合理的安排与准备，以应对各种随机事件，降低光电并网后带来的不良影响，可以保障电网的安全性、稳定性、和经济性。
3.目前，可以从物理层面和统计层面进行光伏发电功率预测。站在物理层面主要考虑太阳的辐照时间、辐照强度、气温等气象数据以及地理信息等影响因素，然后根据输出功率曲线计算出光伏电站的输出功率，其过程比较复杂。站在统计层面主要考虑历史数据，基于模型的输入、输出，从数学角度找到内在的统计规律，从而进行预测。其不需要考虑太阳辐照的具体过程和光伏组件的特性原理，简化了预测过程。
4.利用灰色理论预测光伏发电功率，虽然不需要大量样本、工作量小，但是灰色理论对历史数据依赖性强，而且没有考虑影响因素之间的联系。基于马尔科夫链的光伏并网发电量预测，能够较好的适用光伏发电的随机性和波动性，但是要求历史数据“无后效性”且状态的划分没有完善的理论。采用线性回归模型预测光伏发电功率，能够充分考虑环境等因素的影响而且预测方法简单，在外在环境相对稳定的情况下预测精度比较理想，但是在影响因素波动性较大的时间点预测误差较大。也有采用组合权重相似日、径向基函数神经网络(rbf)、灰色神经网络组合模型等统计法对光伏发电功率进行直接预测，但其过程比较繁琐，而且在预测精度上并没有突出的提高。
5.以往的研究大多是面向单尺度的光功率预测，本专利从多尺度的角度提升预测精度。事实上，光伏功率在不同的时间尺度上受外界不同时间尺度的因素影响，故多尺度行为对于光伏功率预测至关重要。为了自动学习光伏功率数据中的多尺度信息，本专利提出了一种利用多尺度下采样和模态信息实现光功率预测的方法。该方法分别利用基于变分模态分解和多尺度下采样的两种尺度信息，通过变分模态分解获取同一个特征中的不同时间尺度信息，通过下采样获取不同特征的中的相同尺度信息。同时变分模态特征进行模态卷积网络融合，并使用光伏真实值对模态卷积网络进行预训练，对6个尺度信号通过6个尺度卷积网络进行融合，并使用光伏真实值对尺度卷积网络进行预训练，最后通过全连接网络将模态卷积网络和尺度卷积网络进行训练，并采用真实值对整个网络进行训练，从而多尺度信息可以显著提高光伏预测的性能，且具有较强的泛化性能，具有较好的实用性。


技术实现要素：

6.有鉴于此，本发明的第一方面的目的是提出一种基于变分模态分解和多尺度下采
样的光伏功率出力预测方法，从而提高光伏预测的精度。
7.本发明的目的是通过以下技术方案实现的：
8.本发明的一种基于变分模态分解和多尺度下采样的光伏功率出力预测方法，包括以下步骤：
9.步骤s1:进行数据准备
10.对某个历史时刻t0，获取前m个时刻的样本集，m为正整数；其中光伏功率序列为f
power
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]，相应的，所有序列对应的label均为f
power
(t0 1)，样本序列一共有n个；
[0011]
步骤s2：进行模态信息融合
[0012]
包括以下子步骤：
[0013]
步骤s21：进行特征模态分解；
[0014]
步骤s22：进行模态向量拼接；
[0015]
步骤s23:构建模态卷积网络netv；
[0016]
步骤s24：进行模态卷积网络预训练；即以均方误差为损失函数，采用随机梯度下降法训练netv，直至收敛；
[0017]
步骤s25：提取netv的全连接层输入结果，记为
[0018]
步骤s3：进行尺度信息融合，包括：
[0019]
步骤s31:构建不同尺度采样池序列；
[0020]
步骤s32：确定下采样尺度数量；
[0021]
步骤s33:进行尺度卷积网络设计；根据样本序列的数量，设计的卷积神经网络为步骤s33:进行尺度卷积网络设计；根据样本序列的数量，设计的卷积神经网络为
[0022]
步骤s34:进行尺度卷积网络预训练；
[0023]
步骤s35：尺度卷积网络数据拼接；即将步骤s34训练得到的多个网络的最后一层输出结果拼接成一个向量形成
[0024]
步骤s4：进行模态-尺度信息融合，包括：
[0025]
步骤s41:将下采样预测结果和变分模态预测结果拼接成一个新的向量x
joint
，即
[0026]
步骤s42：构建模态-尺度信息融合，该网络为一个全连接网络netf；
[0027]
步骤s43:进行网络训练，采用随机梯度下降法训练netf；
[0028]
步骤s5：包括如下子步骤：
[0029]
步骤s51：获取实时的样本；
[0030]
步骤s52：获取前m个时刻的样本的数值；
[0031]
步骤s53：对前m个时刻的样本的数值进行变分模态分解并拼接成列向量；
[0032]
步骤s54：将样本以及变分模态分解并拼接成列向量都输入网络即可得到下一个时刻的光伏出力预测值。
[0033]
进一步的，所述步骤s1中，所述样本包括但不限于温度、云量、雨量、风速、光照、光伏功率的样本，其中温度序列设定为f
temp
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]，其他样本序列以此
类推。
[0034]
进一步，所述步骤s21中，以温度序列f
temp
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]为对象，进行变分模态分解，包括以下子步骤：
[0035]
步骤s211：设置参数：
[0036]
最大模态数量k＝4；设置阈值ε＝0.01
[0037]
步骤s212:初始化：
[0038]
初始化和λ0，这里k≤k；迭代次数n＝0；(上述参数分别代表何种变量，需要说明)
[0039]
步骤s213：计算模态分量频域值：
[0040]
对进行快速傅里叶变换得到对应的频域分量
[0041]
对λ0进行快速傅里叶变换得到对应的频域分量
[0042]
对f
temp
(t)进行快速傅里叶变换得到对应的频域分量
[0043]
步骤s214：进行第n次迭代，对第k个模态，更新采用如下形式更新傅里叶变换第n次计算结果
[0044][0045]
步骤s215：进行第n次迭代，对第k个模态，更新采用如下形式更新ωk的傅里叶变换第n次计算结果
[0046][0047]
上述形式为积分形式，可以采用离散求和形式逼近。
[0048]
步骤s216：进行第n次迭代，对第k个模态，更新采用如下形式更新λ的傅里叶变换第n次计算结果
[0049][0050]
步骤s217:若满足停止迭代，否则n＝n 1继续迭代；
[0051]
步骤s218：反变换：
[0052]
对第k个模态，将频域信号反变换至模态时域分量uk(t)，并令
[0053]
通过上述步骤如此得到对应的变分模态分解结果：
[0054][0055]
同理可得其他序列变分模态分解结果。
[0056]
进一步，所述s22中，是将各样本序列的变分模态分解结果拼接至行向量，然后将各个行向量拼接成矩阵，记为u。
[0057]
进一步，步骤s23包括以下子步骤：
[0058]
步骤s231:将u进行pool1和pad操作；
[0059]
步骤s232：对u采用cnn进行融合，记该模态卷积网络为netv，netv两层卷积层和两层全连接层组成；
[0060]
步骤s233：将第二池化层结果拉直成一维。
[0061]
进一步，所述步骤s31中，是对历史某个时刻t0，获取前m个时刻的样本集，其中光伏功率序列为f
power
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；将所有序列拼接成维数为m
×
n的总特征矩阵f0(t0),其中n为序列的数量；
[0062]
令δt＝1为下采样单位步长，计s为下采样尺度，s为总的下采样尺度数量，对f0(t0)进行尺度s下的下采样操作：
[0063]fs
(t0)＝{f0(t sδt),t sδt≤t0}。
[0064]
进一步，所述步骤s32中，采用x个时间尺度足以保证下采样序列具有充分的信息，故取s＝x，如此便可得到不同尺度下的采样池序列fc(t0)：
[0065]
fd(t0)＝{fs(t0)},
[0066]
进一步，在步骤s34中，仍以f
power
(t0 1)为前述网络的label，并以均方误差为损失函数，采用随机梯度下降法训练直至收敛；
[0067]
在步骤35中，尺度卷积网络数据拼接是将共n个网络的最后一层输出结果拼接成一个向量形成
[0068]
进一步，所述步骤s43中，冻结netv，共n 1个网络的参数，以f
power
(t0 1)为上述netf网络的label，并以均方误差为损失函数，采用随机梯度下降法训练netf；
[0069]
冻结netf的参数，以f
power
(t0 1)为netf的label，采用随机梯度下降法训练的label，采用随机梯度下降法训练网络；
[0070]
重复直至收敛。
[0071]
本发明的第二方面的目的是提供一种计算机装置，包括存储器、处理器及储存在存储器上并能够在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如前所述的方法。
[0072]
本发明的有益效果是：方法分别利用基于变分模态分解和多尺度下采样的两种尺度信息，通过变分模态分解获取同一个特征中的不同时间尺度信息，通过下采样获取不同
特征的中的相同尺度信息。同时变分模态特征进行模态卷积网络融合，并使用光伏真实值对模态卷积网络进行预训练，对多个尺度信号通过对应数目的尺度卷积网络进行融合，并使用光伏真实值对尺度卷积网络进行预训练，最后通过全连接网络将模态卷积网络和尺度卷积网络进行训练，并采用真实值对整个网络进行训练，从而多尺度信息可以显著提高光伏预测的性能，且具有较强的泛化性能，具有良好的实用性。
[0073]
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和前述的权利要求书来实现和获得。
附图说明
[0074]
为了使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明作进一步的详细描述，其中：
[0075]
附图1为下采样环节示意图；
[0076]
附图2为尺度卷积网络示意图；
[0077]
附图3为单个特征行向量示意图；
[0078]
附图4为温度模态分解和拼接示意图；
[0079]
附图5为模态卷积网络示意图；
[0080]
附图6为模态卷积网络和尺度卷积网络的融合示意图。
具体实施方式
[0081]
以下将参照附图，对本发明的优选实施例进行详细的描述。应当理解，优选实施例仅为了说明本发明，而不是为了限制本发明的保护范围。
[0082]
本发明旨在提出一种基于变分模态分解和多尺度下采样的光伏功率出力预测方法，从而提高光伏预测的精度。
[0083]
为了实现上述技术目的，本发明的步骤如下：
[0084]
为叙述方便，此处先对一些符号和变量做统一说明：本实施例中，一共采用了温度、云量、雨量、风速、光照、光伏功率在内的6个样本，也即是n＝6；其中，
[0085]
记当前时刻为t；
[0086]
计温度序列为f
temp
(t)；
[0087]
云量序列为f
cloud
(t)；
[0088]
雨量序列为f
rain
(t)；
[0089]
风速序列为f
wind
(t)；
[0090]
光照序列为f
light
(t)；
[0091]
光伏功率序列为f
power
(t)
[0092]
本专利的目的是通过上述序列预测下一个时刻的光伏出力f
power
(t 1)，采用的技术手段如下：
[0093]
1、预测影响因素选择
[0094]
光伏发电主要依赖于太阳光照，而不同的地理位置、不同的季节太阳的辐照强度
和辐射时间都有很大的不同。所以在预测时应该根据季节特性把整体划分为不同的预测阶段，在每个阶段中，不仅太阳光照条件类似，而且其他的气象条件也都类似。同时，预测阶段的划分和回归模型的理论特性也相吻合，能够在很大程度上提高预测精度。故本专利除了太阳辐照强度外，还有许多影响光伏发电的因素。对观测数据的分析，还选取其中主要的几个因素：风速、温度、云量、降雨量。
[0095]
2、卷积神经网络
[0096]
卷积神经网络通过增加网络的层数来模拟人类视觉系统的信息处理过程。从最底层到最顶层逐渐提取图像信息的有效抽象化特征表示，即越高层数上的图像特征具有越强的表征能力。卷积神经网络可以归结为输入层、多个隐含层和输出层，其通过共享网络参数来减少网络间的权重，以此来降低训练样本数量和提高训练速度，为了降低特征维数和避免出现过度拟合(over-fitting)，cnns需要进行池化操作。池化操作可以提高特征描述的鉴别力，常用的池化操作是平均池化(mean-pooling)、最大池化(max-pooling)和随机池化(stochastic pooling)操作。
[0097]
3、变分模态分解
[0098]
变分模态分解的目标是将输入信号分解为离散数量的子模态，其中每个模态在频谱域中具有有限的带宽。即每个模态主要紧凑地围绕中心频率k，其具体方法是：
[0099]
(1)设对第k个变分模态分量为uk(t)，利用希尔伯特变换计算相关的解析信号，以获得单边频谱
[0100][0101]
其中f(ω)是第k个变分模态分量为uk(t)平移至ωk对应的傅里叶变换，δ(t)是狄拉克分布函数，t是时间，j是复数单位，下标k是第k个变分模态分量，*是卷积算子。本专利中最大模态数k设置为4，以下为叙述方便，记uk(t)为uk。
[0102]
(2)对于每个模态，通过与调谐到相应估计中心频率的指数混合，将模态的频谱移至“基带”：
[0103][0104]
其中，f(ω-ωk)是第k个变分模态分量为uk(t)平移至ωk对应的傅里叶变换，ωk为第k个变分模态分量对应的中心频率。
[0105]
通过解调信号的平滑度估计带宽。由此产生的约束变分问题如下：
[0106][0107][0108]
式中{uk}＝{u1,...,uk}表示k个模态分量，{ωk}＝{ω1,...,ωk}表示k个模态分量的中心频率，表示微分算子，表示二范数平方和，f表示原始信号。
[0109]
本专利采用admm框架来处理上述优化问题。遵循admm的标准格式，本专利通过lagrange乘子法，将有约束变分优化变为无约束优化，其原理如下：
其中λ(t)是拉格朗日乘子，以下简记λ(t)为λ。
[0110]
推导可得第k个模态第n 1次迭代结果为：
[0111]
将上述等式变换到频域：
[0112]
其中是的傅里叶变换，sgn(ω ωk)是δ(t)的傅里叶变换，是uk(t)的傅里叶变换，是f(t)的傅里叶变换，是λ(t)的傅里叶变换。
[0113]
进一步的推导可得上述优化问题的解具有如下解析形式：
[0114][0115]
同理，对于的最小值问题，在频域通过优化问题解得中心频率为：
[0116][0117]
下采样：对于一个样值序列间隔几个样值取样一次,这样得到新序列就是原序列的下采样。本专利中下需要不同采样率主要是有以下两方面的需求：
[0118]
(1)降低计算量；
[0119]
(2)捕获不同特征在尺度上的信息。
[0120]
综上所述，本专利的光伏预测方法流程如下：
[0121]
步骤s1:数据准备环节
[0122]
对历史某个时刻t0，获取前m个时刻的样本集(本专利取m＝15)：
[0123]
温度序列为f
temp
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0124]
云量序列为f
cloud
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0125]
雨量序列为f
rain
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0126]
风速序列为f
wind
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0127]
光照序列为f
light
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0128]
光伏功率序列为f
power
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]
[0129]
相应的，上述6序列对应的label均为f
power
(t0 1)
[0130]
步骤s2：模态信息融合
[0131]
步骤s21：进行特征模态分解
[0132]
对温度序列f
temp
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]，进行变分模态分解，其计算步骤为：
[0133]
步骤s211：设置参数：本实施例中，设置最大模态数量k＝4；设置阈值ε＝0.01
[0134]
步骤s212:初始化：
[0135]
初始化和λ0，本实施例中，k≤k；迭代次数n＝0；为第k个模态迭代初始值，为第k个模态初始迭代频率，λ0为拉格朗日乘子的初始值
[0136]
步骤s213：计算模态分量频域值：
[0137]
对进行快速傅里叶变换得到对应的频域分量
[0138]
对λ0进行快速傅里叶变换得到对应的频域分量
[0139]
对f
temp
(t)进行快速傅里叶变换得到对应的频域分量
[0140]
步骤s214：进行第n次迭代，对第k个模态，更新采用如下形式更新傅里叶变换第n次计算结果
[0141][0142]
步骤s215：进行第n次迭代，对第k个模态，更新采用如下形式更新ωk的傅里叶变换第n次计算结果
[0143][0144]
上述形式为积分形式，可以采用离散求和形式逼近。
[0145]
步骤s216：进行第n次迭代，对第k个模态，更新采用如下形式更新λ的傅里叶变换第n次计算结果
[0146][0147]
步骤s217:若满足停止迭代，否则n＝n 1继续迭代；
[0148]
步骤s218：反变换：
[0149]
对第k个模态，将频域信号反变换至模态时域分量uk(t)，并令
[0150]
通过上述步骤如此得到对应的变分模态分解结果：
[0151][0152]
同理可可得云量序列变分模态分解结果结果；
[0153][0154]
同理可得雨量序列变分模态分解结果：
[0155][0156]
同理可得风速列进行变分模态分解结果：
[0157][0158]
同理可得光照序列变分模态分解结果：
[0159][0160]
同理可得功率序列变分模态分解结果：
[0161][0162]
步骤ss2：模态向量拼接
[0163]
将拼接至行向量，简记
[0164]
将拼接至行向量，简记
[0165]
将拼接至行向量，简记
[0166]
将拼接至行向量，简记
[0167]
将拼接至行向量，简记
[0168]
将拼接至行向量，简记
[0169]
将拼接成矩阵，记为u；
[0170]
步骤s23：构建模态卷积网络netv:
[0171]
包括子步骤为：
[0172]
s231:将u进行pool1和pad操作，本实施例中，使其大小为32
×
32；需要说明的是，此处的大小是可以根据实际情况进行自行设定的。
[0173]
s232:对u采用cnn(卷积神经网络)进行融合，记该模态卷积网络为netv，netv两层卷积层和两层全连接层组成，具体的结构如下：(需要说明的是，此处的卷积层和全连接层以及具体结构是可以根据实际情况进行自行设定的)
[0174]
第一卷积层conv1输入数据：32
×
32
×
1。
[0175]
第一卷积层conv1滤波器：维度5
×
5，步长1，填充0，滤波器6个。
[0176]
第一池化层pool1输入数据：28
×
28
×
6。
[0177]
第一池化层pool1滤波器：维度2
×
2，步长2，填充0。
[0178]
第二卷积层conv2输入数据：14
×
14
×
6。
[0179]
第二卷积层conv2滤波器：维度5
×
5，步长1，填充0，滤波器16个。
[0180]
第二池化层pool2输入数据：10
×
10
×
16。
[0181]
第二池化层pool2滤波器：维度2
×
2，步长2，填充0。
[0182]
s233:将第二池化层结果拉直成一维；其中，
[0183]
全连接层fc3输入数据：400
×
1，输出数据120
×1[0184]
全连接层fc4输入数据：120
×
1，输出数据84
×1[0185]
softmax回归输入数据：84
×
1，输出数据10
×1[0186]
步骤s24：模态卷积网络预训练
[0187]
以f
power
(t0 1)为上述网络的label，并以均方误差为损失函数，采用随机梯度下降法训练netv，直至收敛；
[0188]
步骤s25：提取netv全连接层fc4输入结果，记为
[0189]
步骤s3：尺度信息融合
[0190]
步骤s31：构建不同尺度采样池序列：
[0191]
对历史某个时刻t0，获取前m个时刻的样本集(本专利取m＝15)：
[0192]
计温度序列为f
temp
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0193]
云量序列为f
cloud
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0194]
雨量序列为f
rain
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0195]
风速序列为f
wind
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0196]
光照序列为f
ligh
t(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0197]
光伏功率序列为f
power
(t),t∈[t-m 1
,
…
,t-2
,t-1
,t0]；
[0198]
将上述6个序列拼接成维数为m
×
6的总特征矩阵f0(t0)
[0199]
令δt＝1为下采样单位步长，计s为下采样尺度，s为总的下采样尺度数量，对f0(t0)进行尺度s下的下采样操作：
[0200]fs
(t0)＝{f0(t sδt),t sδt≤t0}
[0201]
步骤s32：确定下采样尺度数量
[0202]
光伏预测问题中，一般采样时间为15分钟，故6个时间尺度足以保证下采样序列具有充分的信息，故本专利取s＝6，如此便可得到不同尺度下的采样池序列fc(t0)：
[0203]
fd(t0)＝{fs(t0)},
[0204]
步骤s33：尺度卷积网络设计
[0205]
不同的尺度应该对应不同的网络结构，但不同的网络训练计算量大，收敛性差，故本专利设计了三种不同的卷积神经网络，具体的：
[0206]
对尺度1和2，本专利设计的网络和结构如下：
[0207]
类型滤波器数量滤波器尺寸strip步长卷积层147*72最大池化层142*21卷积层257*71最大池化层252*21全连接层2561*1无
[0208]
需要注意的是，和只是结构相同，但参数并不相同，下同。
[0209]
对尺度3和4，本专利设计的网络和结构如下：
[0210]
类型滤波器数量滤波器尺寸strip步长卷积层1165*51最大池化层142*23卷积层283*31最大池化层243*31全连接层641*1无
[0211]
对于尺度5和6，本专利设计的网络和结构如下：
[0212]
类型滤波器数量滤波器尺寸strip步长卷积层1532*322最大池化层182*21卷积层2616*161最大池化层268*81全连接层321*1无
[0213]
步骤s34：尺度卷积网络预训练：
[0214]
仍以f
power
(t0 1)为上述网络6个的label，并以均方误差为损失函数，采用随机梯度下降法训练直至收敛；
[0215]
步骤s35：尺度卷积网络数据拼接：
[0216]
将6个网络的最后一层输出结果拼接成一个向量形成
[0217]
步骤s4：模态-尺度信息融合
[0218]
步骤s41：将下采样预测结果和变分模态预测结果拼接成一个新的向量x
joint
，即
[0219]
步骤s42：构建模态-尺度信息融合，该网络为一个全连接网络netf，其网络结构为：
[0220]
类型隐藏节点数是否有bn操作全连接层132有全连接层215有全连接层38无全连接层41无
[0221]
步骤s43：网络训练
[0222]
冻结netv，7个网络的参数，以f
power
(t0 1)为上述netf网络的label，并以均方误差为损失函数，采用随机梯度下降法训练netf；
[0223]
冻结netf的参数，以f
power
(t0 1)为netf的label，采用随机梯度下降法训练netv，
7个网络；
[0224]
重复直至收敛。
[0225]
五：在线预测
[0226]
步骤1：获取实时温度、云量、雨量、风速、光照以及光功率；
[0227]
步骤2：获取前m个时刻的温度、云量、雨量、风速、光照以及光功率；本实施例中，m取15，即取前15个时刻的温度、云量、雨量、风速、光照以及光功率的数值。
[0228]
步骤3：对前m个时刻的温度、云量、雨量、风速、光照以及光功率的进行变分模态分解并拼接成列向量；
[0229]
步骤4：将实时温度、云量、雨量、风速、光照以及光功率以及变分模态分解并拼接成列向量都输入网络即可得到下一个时刻的光伏出力预测值。
[0230]
应当认识到，本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现，其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而，若需要，该程序可以以汇编或机器语言实现。在任何情况下，该语言可以是编译或解释的语言。此外，为此目的该程序能够在编程的专用集成电路上运行。
[0231]
此外，可按任何合适的顺序来执行本文描述的过程的操作，除非本文另外指示或以其他方式明显地与上下文矛盾。本文描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机系统的控制下执行，并且可作为共同地在一个或多个处理器上执行的代码(例如，可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。
[0232]
进一步，所述方法可以在可操作地连接至合适的任何类型的计算平台中实现，包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本发明的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现，无论是可移动的还是集成至计算平台，如硬盘、光学读取和/或写入存储介质、ram、rom等，使得其可由可编程计算机读取，当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外，机器可读代码，或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时，本文所述的发明包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本发明所述的的方法和技术编程时，本发明还包括计算机本身。
[0233]
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。