1.一种基于三维梯度幅度的立体图像质量客观评价方法，其特征在于包括以下步骤：

①令Sorg表示原始的无失真的立体图像，令Sdis表示待评价的失真的立体图像，将Sorg的左视点图像记为{Lorg(x,y)}，将Sorg的右视点图像记为{Rorg(x,y)}，将Sdis的左视点图像记为{Ldis(x,y)}，将Sdis的右视点图像记为{Rdis(x,y)}，其中，(x,y)表示左视点图像和右视点图像中的像素点的坐标位置，1≤x≤W，1≤y≤H，W表示左视点图像和右视点图像的宽度，H表示左视点图像和右视点图像的高度，Lorg(x,y)表示{Lorg(x,y)}中坐标位置为(x,y)的像素点的像素值，Rorg(x,y)表示{Rorg(x,y)}中坐标位置为(x,y)的像素点的像素值，Ldis(x,y)表示{Ldis(x,y)}中坐标位置为(x,y)的像素点的像素值，Rdis(x,y)表示{Rdis(x,y)}中坐标位置为(x,y)的像素点的像素值；

②根据{Lorg(x,y)}中的每个像素点和{Rorg(x,y)}中对应坐标位置的像素点在多个视差值下的视差空间值，获得Sorg的视差空间图，记为{DSIorg(x,y,d)}，并根据{Ldis(x,y)}中的每个像素点和{Rdis(x,y)}中对应坐标位置的像素点在多个视差值下的视差空间值，获得Sdis的视差空间图，记为{DSIdis(x,y,d)}，其中，DSIorg(x,y,d)表示{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的视差空间值，DSIdis(x,y,d)表示{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的视差空间值，0≤d≤dmax，dmax表示最大视差值；

③计算{DSIorg(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度和视点方向梯度，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的水平方向梯度记为gxorg(x,y,d)，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的垂直方向梯度记为gyorg(x,y,d)，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的视点方向梯度记为gdorg(x,y,d)；

同样，计算{DSIdis(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度和视点方向梯度，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的水平方向梯度记为gxdis(x,y,d)，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的垂直方向梯度记为gydis(x,y,d)，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的视点方向梯度记为gddis(x,y,d)；

④根据{DSIorg(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度

和视点方向梯度，计算{DSIorg(x,y,d)}中的每个像素点的三维梯度幅度，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的三维梯度幅度记为morg(x,y,d)，同样，根据{DSIdis(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度和视点方向梯度，计算{DSIdis(x,y,d)}中的每个像素点的三维梯度幅度，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的三维梯度幅度记为mdis(x,y,d)，⑤根据{DSIorg(x,y,d)}和{DSIdis(x,y,d)}中的每个像素点的三维梯度幅度，计算{DSIdis(x,y,d)}中的每个像素点的客观评价度量值，将{DSIdis(x,y,d)}中 坐 标 位 置 为(x,y,d) 的 像 素 点 的 客 观 评 价 度 量 值 记 为QDSI(x,y,d)，其中，C为控制参数；

⑥根据{DSIdis(x,y,d)}中的每个像素点的客观评价度量值，计算Sdis的图像质量客观评价预测值，记为Q， 其中，Ω表示{DSIdis(x,y,d)}中的所有像素点的坐标位置的集合，N表示{DSIdis(x,y,d)}中包含的像素点的总个数。

2.根据权利要求1所述的一种基于三维梯度幅度的立体图像质量客观评价方法，其特征在于所述的步骤②中Sorg的视差空间图的获取过程为：

②-a1、将{Lorg(x,y)}中当前待处理的像素点定义为当前第一像素点，将{Rorg(x,y)}中当前待处理的像素点定义为当前第二像素点；

②-a2、假设当前第一像素点为{Lorg(x,y)}中坐标位置为(x1,y1)的像素点，假设当前第二像素点为{Rorg(x,y)}中坐标位置为(x1,y1)的像素点，取视差值d0=0，然后计算当前第一像素点和当前第二像素点在该视差值d0下的视差空间值，记为DSIorg(x1,y1,d0)，DSIorg(x1,y1,d0)=|Lorg(x1,y1)-Rorg(x1-d0,y1)|，其中，1≤x1≤W，1≤y1≤H，0≤d0≤dmax，dmax表示最大视差值，Lorg(x1,y1)表示{Lorg(x,y)}中坐标位置为(x1,y1)的像素点的像素值，Rorg(x1-d0,y1)表示{Rorg(x,y)}中坐标位置为(x1-d0,y1)的像素点的像素值，“||”为取绝对值符号；

②-a3、选取dmax个与d0各不相同的视差值，分别记为 然后分别计

算当前第一像素点和当前第二像素点在该dmax个各不相同的视差值下的视差空间值，对应的分别记为

DSIorg(x1,y1,d1)=|Lorg(x1,y1)-Rorg(x1-d1,y1)|，DSIorg(x1,y1,d2)=|Lorg(x1,y1)-Rorg(x1-d2,y1)|，DSIorg(x1,y1,di)=|Lorg(x1,y1)-Rorg(x1-di,y1)|，其中，1≤i≤dmax，di=d0+i，DSIorg(x1,y1,d1)表示当前第一像素点和当前第二像素点在视差值d1下的视差空间值，DSIorg(x1,y1,d2)表示当前第一像素点和当前第二像素点在视差值d2下的视差空间值，DSIorg(x1,y1,di)表示当前第一像素点和当前第二像素点在视差值di下的视差空间值，表示当前第一像素点和当前第二像素点在视差值 下的视差空间值，Rorg(x1-d1,y1)表示{Rorg(x,y)}中坐标位置为(x1-d1,y1)的像素点的像素值，Rorg(x1-d2,y1)表示{Rorg(x,y)}中坐标位置为(x1-d2,y1)的像素点的像素值，Rorg(x1-di,y1)表示{Rorg(x,y)}中坐标位置为(x1-di,y1)的像素点的像素值， 表示{Rorg(x,y)}中坐标位置为 的像素点的像素值；

②-a4、将{Lorg(x,y)}中下一个待处理的像素点作为当前第一像素点，将{Rorg(x,y)}中下一个待处理的像素点作为当前第二像素点，然后返回步骤②-a2继续执行，直至{Lorg(x,y)}和{Rorg(x,y)}中的所有像素点处理完毕，获得Sorg的视差空间图，记为{DSIorg(x,y,d)}，其中，DSIorg(x,y,d)表示{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的视差空间值，DSIorg(x,y,d)的值为{Lorg(x,y)}中坐标位置为(x,y)的像素点和{Rorg(x,y)}中坐标位置为(x,y)的像素点在视差值d下的视差空间值，所述的步骤②中Sdis的视差空间图的获取过程为：

②-b1、将{Ldis(x,y)}中当前待处理的像素点定义为当前第一像素点，将{Rdis(x,y)}中当前待处理的像素点定义为当前第二像素点；

②-b2、假设当前第一像素点为{Ldis(x,y)}中坐标位置为(x1,y1)的像素点，假设当前第二像素点为{Rdis(x,y)}中坐标位置为(x1,y1)的像素点，取视差值d0=0，然后计算当前第一像素点和当前第二像素点在该视差值d0下的视差空间值，记为DSIdis(x1,y1,d0)，DSIdis(x1,y1,d0)=|Ldis(x1,y1)-Rdis(x1-d0,y1)|，其中，1≤x1≤W，1≤y1≤H，0≤d0≤dmax，dmax表示最大视差值，Ldis(x1,y1)表示{Ldis(x,y)}中坐标位置为(x1,y1)的像素点的像素值，Rdis(x1-d0,y1)表示{Rdis(x,y)}中坐标位置为(x1-d0,y1)的像素点的像素值，“||”为取绝对值符号；

②-b3、选取dmax个与d0各不相同的视差值，分别记为 然后分别计

算当前第一像素点和当前第二像素点在该dmax个各不相同的视差值下的视差空间值，对应的分别记为

DSIdis(x1,y1,d1)=|Ldis(x1,y1)-Rdis(x1-d1,y1)|，DSIdis(x1,y1,d2)=|Ldis(x1,y1)-Rdis(x1-d2,y1)|，DSIdis(x1,y1,di)=|Ldis(x1,y1)-Rdis(x1-di,y1)|，其中，1≤i≤dmax，di=d0+i，DSIdis(x1,y1,d1)表示当前第一像素点和当前第二像素点在视差值d1下的视差空间值，DSIdis(x1,y1,d2)表示当前第一像素点和当前第二像素点在视差值d2下的视差空间值，DSIdis(x1,y1,di)表示当前第一像素点和当前第二像素点在视差值di下的视差空间值，表示当前第一像素点和当前第二像素点在视差值 下的视差空间值，Rdis(x1-d1,y1)表示{Rdis(x,y)}中坐标位置为(x1-d1,y1)的像素点的像素值，Rdis(x1-d2,y1)表示{Rdis(x,y)}中坐标位置为(x1-d2,y1)的像素点的像素值，Rdis(x1-di,y1)表示{Rdis(x,y)}中坐标位置为(x1-di,y1)的像素点的像素值， 表示{Rdis(x,y)}中坐标位置为 的像素点的像素值；

②-b4、将{Ldis(x,y)}中下一个待处理的像素点作为当前第一像素点，将{Rdis(x,y)}中下一个待处理的像素点作为当前第二像素点，然后返回步骤②-b2继续执行，直至{Ldis(x,y)}和{Rdis(x,y)}中的所有像素点处理完毕，获得Sdis的视差空间图，记为{DSIdis(x,y,d)}，其中，DSIdis(x,y,d)表示{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的视差空间值，DSIdis(x,y,d)的值为{Ldis(x,y)}中坐标位置为(x,y)的像素点和{Rdis(x,y)}中坐标位置为(x,y)的像素点在视差值d下的视差空间值，

3.根据权利要求1或2所述的一种基于三维梯度幅度的立体图像质量客观评价方法，其特征在于所述的步骤③中{DSIorg(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度和视点方向梯度的获取过程为：③-a1、采用水平梯度算子对{DSIorg(x,y,d)}进行卷积，得到{DSIorg(x,y,d)}中的每个像素点的水平方向梯度，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的水平方向梯度记为gxorg(x,y,d)，其中，DSIorg(u,v,j)表示{DSIorg(x,y,d)}中坐标位置为(u,v,j)的像素点的视差空间值；

③-a2、采用垂直梯度算子对{DSIorg(x,y,d)}进行卷积，得到{DSIorg(x,y,d)}中的每个像素点的垂直方向梯度，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的垂直方向梯度记为gyorg(x,y,d)，③-a3、采用视点梯度算子对{DSIorg(x,y,d)}进行卷积，得到{DSIorg(x,y,d)}中的每个像素点的视点方向梯度，将{DSIorg(x,y,d)}中坐标位置为(x,y,d)的像素点的视点方向梯度记为gdorg(x,y,d)， 其中，sign()为阶跃函数，

上述步骤③-a1至步骤③-a3中，如果u<1，则DSIorg(u,v,j)的值由DSIorg(1,v,j)的值替代，如果u>W，则DSIorg(u,v,j)的值由DSIorg(W,v,j)的值替代，如果v<1，则DSIorg(u,v,j)的值由DSIorg(u,1,j)的值替代，如果v>H，则DSIorg(u,v,j)的值由DSIorg(u,H,j)的值替代，如果j<0，则DSIorg(u,v,j)的值由DSIorg(u,v,0)的值替代，如果j>dmax，则DSIorg(u,v,j)的值由DSIorg(u,v,dmax)的值替代，DSIorg(1,v,j)表示{DSIorg(x,y,d)}中坐标位置为(1,v,j)的像素点的视差空间值，DSIorg(W,v,j)表示{DSIorg(x,y,d)}中坐标位置为(W,v,j)的像素点的视差空间值，DSIorg(u,1,j)表示{DSIorg(x,y,d)}中坐标位置为(u,1,j)的像素点的视差空间值，DSIorg(u,H,j)表示{DSIorg(x,y,d)}中坐标位置为(u,H,j)的像素点的视差空间值，DSIorg(u,v,0)表示{DSIorg(x,y,d)}中坐标位置为(u,v,0)的像素点的视差空间值，DSIorg(u,v,dmax)表示{DSIorg(x,y,d)}中坐标位置为(u,v,dmax)的像素点的视差空间值；

所述的步骤③中{DSIdis(x,y,d)}中的每个像素点的水平方向梯度、垂直方向梯度和视点方向梯度的获取过程为：

③-b1、采用水平梯度算子对{DSIdis(x,y,d)}进行卷积，得到{DSIdis(x,y,d)}中的每个像素点的水平方向梯度，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的水平方向梯度记为gxdis(x,y,d)，其中，DSIdis(u,v,j)表示{DSIdis(x,y,d)}中坐标位置为(u,v,j)的像素点的视差空间值；

③-b2、采用垂直梯度算子对{DSIdis(x,y,d)}进行卷积，得到{DSIdis(x,y,d)}中的每个像素点的垂直方向梯度，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的垂直方向梯度记为gydis(x,y,d)，③-b3、采用视点梯度算子对{DSIdis(x,y,d)}进行卷积，得到{DSIdis(x,y,d)}中的每个像素点的视点方向梯度，将{DSIdis(x,y,d)}中坐标位置为(x,y,d)的像素点的视差方向梯度记为gddis(x,y,d)， 其中，sign()为阶跃函数，

上述步骤③-b1至步骤③-b3中，如果u<1，则DSIdis(u,v,j)的值由DSIdis(1,v,j)的值替代，如果u>W，则DSIdis(u,v,j)的值由DSIdis(W,v,j)的值替代，如果v<1，则DSIdis(u,v,j)的值由DSIdis(u,1,j)的值替代，如果v>H，则DSIdis(u,v,j)的值由DSIdis(u,H,j)的值替代，如果j<0，则DSIdis(u,v,j)的值由DSIdis(u,v,0)的值替代，如果j>dmax，则DSIdis(u,v,j)的值由DSIdis(u,v,dmax)的值替代，DSIdis(1,v,j)表示{DSIdis(x,y,d)}中坐标位置为(1,v,j)的像素点的视差空间值，DSIdis(W,v,j)表示{DSIdis(x,y,d)}中坐标位置为(W,v,j)的像素点的视差空间值，DSIdis(u,1,j)表示{DSIdis(x,y,d)}中坐标位置为(u,1,j)的像素点的视差空间值，DSIdis(u,H,j)表示{DSIdis(x,y,d)}中坐标位置为(u,H,j)的像素点的视差空间值，DSIdis(u,v,0)表示{DSIdis(x,y,d)}中坐标位置为(u,v,0)的像素点的视差空间值，DSIdis(u,v,dmax)表示{DSIdis(x,y,d)}中坐标位置为(u,v,dmax)的像素点的视差空间值。