Purpose: updated the core chapter
authorElena Fedotova <no@email>
Sun, 10 Apr 2011 19:37:17 +0000 (19:37 +0000)
committerElena Fedotova <no@email>
Sun, 10 Apr 2011 19:37:17 +0000 (19:37 +0000)
modules/core/doc/operations_on_arrays.rst

index 84655fd..ef4ff7f 100644 (file)
@@ -27,7 +27,8 @@ abs
     
 The output matrix has the same size and the same type as the input one (except for the last case, where ``C`` will be ``depth=CV_8U`` ).
 
-See Also: :ref:`MatrixExpressions`, :func:`absdiff`
+See Also: :ref:`MatrixExpressions`, 
+:func:`absdiff`
 
 .. index:: absdiff
 
@@ -83,11 +84,11 @@ add
 
     :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param sc: A scalar. This is the second input parameter.
+    :param sc: Scalar. This is the second input parameter.
 
-    :param dst: The destination array of the same size and type as  ``src1`` . See  ``Mat::create`` .
+    :param dst: Destination array of the same size and type as  ``src1`` . See  ``Mat::create`` .
     
-    :param mask: An optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
+    :param mask: Optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
 
 The functions ``add`` compute:
 
@@ -116,7 +117,11 @@ The first function in the list above can be replaced with matrix expressions: ::
 In case of multi-channel arrays, each channel is processed independently.
 
 See Also:
-:func:`subtract`,:func:`addWeighted`,:func:`scaleAdd`,:func:`convertScale`,:ref:`MatrixExpressions`
+:func:`subtract`,
+:func:`addWeighted`,
+:func:`scaleAdd`,
+:func:`convertScale`,
+:ref:`MatrixExpressions`
 
 .. index:: addWeighted
 
@@ -136,9 +141,9 @@ addWeighted
     
     :param beta: Weight for the second array elements.
 
-    :param dst: The destination array of the same size and type as  ``src1`` .
+    :param dst: Destination array of the same size and type as  ``src1`` .
     
-    :param gamma: A scalar added to each sum.
+    :param gamma: Scalar added to each sum.
 
 The functions ``addWeighted`` calculate the weighted sum of two arrays as follows:
 
@@ -156,7 +161,11 @@ The first function can be replaced with a matrix expression: ::
 In case of multi-channel arrays, each channel is processed independently.
 
 See Also:
-:func:`add`,:func:`subtract`,:func:`scaleAdd`,:func:`convertScale`,:ref:`MatrixExpressions`
+:func:`add`,
+:func:`subtract`,
+:func:`scaleAdd`,
+:func:`convertScale`,
+:ref:`MatrixExpressions`
 
 .. index:: bitwise_and
 
@@ -174,11 +183,11 @@ bitwise_and
 
     :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param sc: A scalar. This is the second input parameter.
+    :param sc: Scalar. This is the second input parameter.
 
-    :param dst: The destination array of the same size and type as  ``src1`` . See  ``Mat::create`` .    
+    :param dst: Destination array of the same size and type as  ``src1`` . See  ``Mat::create`` .    
     
-    :param mask: An optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
+    :param mask: Optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
 
 The functions ``bitwise_and`` compute the per-element bit-wise logical conjunction:
 
@@ -210,11 +219,11 @@ bitwise_not
 
     Inverts every bit of an array.
 
-    :param src1: The source array.
+    :param src1: Source array.
 
-    :param dst: The destination array. It is reallocated to be of the same size and type as  ``src`` . See  ``Mat::create`` .
+    :param dst: Destination array. It is reallocated to be of the same size and type as  ``src`` . See  ``Mat::create`` .
     
-    :param mask: An optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
+    :param mask: Optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
 
 The functions ``bitwise_not`` compute per-element bit-wise inversion of the source array:
 
@@ -240,11 +249,11 @@ bitwise_or
 
     :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param sc: A scalar. This is the second input parameter.
+    :param sc: Scalar. This is the second input parameter.
 
-    :param dst: The destination array. It is reallocated to be of the same size and type as  ``src1`` . See  ``Mat::create`` .
+    :param dst: Destination array. It is reallocated to be of the same size and type as  ``src1`` . See  ``Mat::create`` .
     
-    :param mask: An optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
+    :param mask: Optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
 
 The functions ``bitwise_or`` compute the per-element bit-wise logical disjunction:
 
@@ -280,11 +289,11 @@ bitwise_xor
 
     :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param sc: A scalar. This is the second input parameter.
+    :param sc: Scalar. This is the second input parameter.
 
-    :param dst: The destination array. It is reallocated to be of the same size and type as  ``src1`` . See  ``Mat::create`` .
+    :param dst: Destination array. It is reallocated to be of the same size and type as  ``src1`` . See  ``Mat::create`` .
     
-    :param mask: An optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
+    :param mask: Optional operation mask, 8-bit single channel array, that specifies elements of the destination array to be changed.
 
 The functions ``bitwise_xor`` compute the per-element bit-wise logical "exclusive or" operation:
 
@@ -352,7 +361,9 @@ calcCovarMatrix
 The functions ``calcCovarMatrix`` calculate the covariance matrix and, optionally, the mean vector of the set of input vectors.
 
 See Also:
-:func:`PCA`,:func:`mulTransposed`,:func:`Mahalanobis`
+:func:`PCA`,
+:func:`mulTransposed`,
+:func:`Mahalanobis`
 
 .. index:: cartToPolar
 
@@ -427,11 +438,11 @@ compare
 
     :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param value: A scalar value to compare each array element with.
+    :param value: Scalar value to compare each array element with.
 
-    :param dst: The destination array of the same size as  ``src1``  and type= ``CV_8UC1`` .
+    :param dst: Destination array of the same size as  ``src1``  and type= ``CV_8UC1`` .
     
-    :param cmpop: The flag specifying the relation between the elements to be checked.
+    :param cmpop: Flag specifying the relation between the elements to be checked.
 
             * **CMP_EQ** :math:`\texttt{src1}(I) = \texttt{src2}(I)`  or  :math:`\texttt{src1}(I) = \texttt{value}`
             * **CMP_GT** :math:`\texttt{src1}(I) > \texttt{src2}(I)`  or  :math:`\texttt{src1}(I) > \texttt{value}`
@@ -453,7 +464,11 @@ The comparison operations can be replaced with the equivalent matrix expressions
 
 
 See Also:
-:func:`checkRange`,:func:`min`,:func:`max`,:func:`threshold`,:ref:`MatrixExpressions`
+:func:`checkRange`,
+:func:`min`,
+:func:`max`,
+:func:`threshold`,
+:ref:`MatrixExpressions`
 
 .. index:: completeSymm
 
@@ -480,7 +495,8 @@ The function ``completeSymm`` copies the lower half of a square matrix to its an
     :math:`\texttt{mtx}_{ij}=\texttt{mtx}_{ji}`     for
     :math:`i < j`     if ``lowerToUpper=true``
     
-See Also: :func:`flip`,:func:`transpose`
+See Also: :func:`flip`,
+:func:`transpose`
 
 .. index:: convertScaleAbs
 
@@ -493,13 +509,13 @@ convertScaleAbs
 
     Scales, computes absolute values, and converts the result to 8-bit.
 
-    :param src: The source array.
+    :param src: Source array.
 
-    :param dst: The destination array.
+    :param dst: Destination array.
 
-    :param alpha: An optional scale factor.
+    :param alpha: Optional scale factor.
 
-    :param beta: An optional delta added to the scaled values.
+    :param beta: Optional delta added to the scaled values.
 
 On each element of the input array, the function ``convertScaleAbs`` performs three operations sequentially: scaling, taking an absolute value, conversion to an unsigned 8-bit type:
 
@@ -518,7 +534,8 @@ In case of multi-channel arrays, the function processes each channel independent
 
 
 See Also:
-:func:`Mat::convertTo`,:func:`abs`
+:func:`Mat::convertTo`,
+:func:`abs`
 
 .. index:: countNonZero
 
@@ -540,7 +557,11 @@ The function ``cvCountNonZero`` returns the number of non-zero elements in ``mtx
     \sum _{I: \; \texttt{mtx} (I) \ne0 } 1
 
 See Also:
-:func:`mean`,:func:`meanStdDev`,:func:`norm`,:func:`minMaxLoc`,:func:`calcCovarMatrix`
+:func:`mean`,
+:func:`meanStdDev`,
+:func:`norm`,
+:func:`minMaxLoc`,
+:func:`calcCovarMatrix`
 
 .. index:: cubeRoot
 
@@ -622,7 +643,12 @@ The last parameter, ``coiMode`` , specifies how to deal with an image with COI s
 :func:`insertImageCOI` , respectively).
 
 See Also:
-:func:`cvGetImage`,:func:`cvGetMat`,:func:`cvGetMatND`,:func:`extractImageCOI`,:func:`insertImageCOI`,:func:`mixChannels` 
+:func:`cvGetImage`,
+:func:`cvGetMat`,
+:func:`cvGetMatND`,
+:func:`extractImageCOI`,
+:func:`insertImageCOI`,
+:func:`mixChannels` 
 
 .. index:: dct
 
@@ -632,21 +658,21 @@ dct
 -------
 .. c:function:: void dct(const Mat& src, Mat& dst, int flags=0)
 
-    Performs a forward or inverse discrete cosine transform of 1D or 2D array
+    Performs a forward or inverse discrete Cosine transform of 1D or 2D array.
 
-    :param src: The source floating-point array
+    :param src: The source floating-point array.
 
-    :param dst: The destination array; will have the same size and same type as  ``src``
+    :param dst: The destination array of the same size and type as  ``src`` .
     
-    :param flags: Transformation flags, a combination of the following values
+    :param flags: Transformation flags, a combination of the following values:
 
-            * **DCT_INVERSE** do an inverse 1D or 2D transform instead of the default forward transform.
+            * **DCT_INVERSE** Perform an inverse 1D or 2D transform instead of the default forward transform.
 
-            * **DCT_ROWS** do a forward or inverse transform of every individual row of the input matrix. This flag allows user to transform multiple vectors simultaneously and can be used to decrease the overhead (which is sometimes several times larger than the processing itself), to do 3D and higher-dimensional transforms and so forth.
+            * **DCT_ROWS** Perform a forward or inverse transform of every individual row of the input matrix. This flag enables you to transform multiple vectors simultaneously and can be used to decrease the overhead (which is sometimes several times larger than the processing itself) to perform 3D and higher-dimensional transforms and so forth.
 
-The function ``dct`` performs a forward or inverse discrete cosine transform (DCT) of a 1D or 2D floating-point array:
+The function ``dct`` performs a forward or inverse discrete Cosine transform (DCT) of a 1D or 2D floating-point array:
 
-Forward Cosine transform of 1D vector of
+* Forward Cosine transform of a 1D vector of
 :math:`N` elements:
 
 .. math::
@@ -663,24 +689,24 @@ and
 :math:`\alpha_0=1`,:math:`\alpha_j=2` for
 :math:`j > 0` .
 
-Inverse Cosine transform of 1D vector of N elements:
+* Inverse Cosine transform of a 1D vector of N elements:
 
 .. math::
 
     X =  \left (C^{(N)} \right )^{-1}  \cdot Y =  \left (C^{(N)} \right )^T  \cdot Y
 
 (since
-:math:`C^{(N)}` is orthogonal matrix,
+:math:`C^{(N)}` is an orthogonal matrix,
 :math:`C^{(N)} \cdot \left(C^{(N)}\right)^T = I` )
 
-Forward Cosine transform of 2D
+Forward Cosine transform of 2D
 :math:`M \times N` matrix:
 
 .. math::
 
     Y = C^{(N)}  \cdot X  \cdot \left (C^{(N)} \right )^T
 
-Inverse Cosine transform of 2D vector of
+* Inverse Cosine transform of a 2D vector of
 :math:`M \times N` elements:
 
 .. math::
@@ -690,29 +716,32 @@ Inverse Cosine transform of 2D vector of
 The function chooses the mode of operation by looking at the flags and size of the input array:
 
 *
-    if ``(flags & DCT_INVERSE) == 0``     , the function does forward 1D or 2D transform, otherwise it is inverse 1D or 2D transform.
+    If ``(flags & DCT_INVERSE) == 0``     , the function does a forward 1D or 2D transform. Otherwise, it is an inverse 1D or 2D transform.
 
 *
-    if ``(flags & DCT_ROWS) :math:`\ne` 0``     , the function performs 1D transform of each row.
+    If ``(flags & DCT_ROWS) :math:`\ne` 0``     , the function performs a 1D transform of each row.
 
 *
-    otherwise, if the array is a single column or a single row, the function performs 1D transform
+    If the array is a single column or a single row, the function performs a 1D transform.
 
 *
-    otherwise it performs 2D transform.
+    If none of the above is true, the function performs a 2D transform.
 
-**Important note**
-: currently dct supports even-size arrays (2, 4, 6 ...). For data analysis and approximation you can pad the array when necessary.
+**Note**
+: 
+Currently ``dct`` supports even-size arrays (2, 4, 6 ...). For data analysis and approximation, you can pad the array when necessary.
 
-Also, the function's performance depends very much, and not monotonically, on the array size, see
-:func:`getOptimalDFTSize` . In the current implementation DCT of a vector of size ``N`` is computed via DFT of a vector of size ``N/2`` , thus the optimal DCT size
+Also, the function performance depends very much, and not monotonically, on the array size (see
+:func:`getOptimalDFTSize` ). In the current implementation DCT of a vector of size ``N`` is computed via DFT of a vector of size ``N/2`` . Thus, the optimal DCT size
 :math:`\texttt{N}^*\geq\texttt{N}` can be computed as: ::
 
     size_t getOptimalDCTSize(size_t N) { return 2*getOptimalDFTSize((N+1)/2); }
 
 
 See Also:
-:func:`dft`,:func:`getOptimalDFTSize`,:func:`idct`
+:func:`dft`,
+:func:`getOptimalDFTSize`,
+:func:`idct`
 
 .. index:: dft
 
@@ -723,24 +752,24 @@ dft
 
 .. c:function:: void dft(const Mat& src, Mat& dst, int flags=0, int nonzeroRows=0)
 
-    Performs a forward or inverse Discrete Fourier transform of 1D or 2D floating-point array.
+    Performs a forward or inverse Discrete Fourier transform of 1D or 2D floating-point array.
 
-    :param src: The source array, real or complex
+    :param src: Source array that could be real or complex.
 
-    :param dst: The destination array, which size and type depends on the  ``flags``
+    :param dst: Destination array whose size and type depends on the  ``flags`` .
     
-    :param flags: Transformation flags, a combination of the following values
+    :param flags: Transformation flags representing a combination of the following values:
 
-            * **DFT_INVERSE** do an inverse 1D or 2D transform instead of the default forward transform.
+            * **DFT_INVERSE** Perform an inverse 1D or 2D transform instead of the default forward transform.
 
-            * **DFT_SCALE** scale the result: divide it by the number of array elements. Normally, it is combined with  ``DFT_INVERSE``             .
-            * **DFT_ROWS** do a forward or inverse transform of every individual row of the input matrix. This flag allows the user to transform multiple vectors simultaneously and can be used to decrease the overhead (which is sometimes several times larger than the processing itself), to do 3D and higher-dimensional transforms and so forth.
+            * **DFT_SCALE** Scale the result: divide it by the number of array elements. Normally, it is combined with  ``DFT_INVERSE`` .             .
+            * **DFT_ROWS** Perform a forward or inverse transform of every individual row of the input matrix. This flag enables you to transform multiple vectors simultaneously and can be used to decrease the overhead (which is sometimes several times larger than the processing itself) to perform 3D and higher-dimensional transforms and so forth.
 
-            * **DFT_COMPLEX_OUTPUT** then the function performs forward transformation of 1D or 2D real array, the result, though being a complex array, has complex-conjugate symmetry ( *CCS* ), see the description below. Such an array can be packed into real array of the same size as input, which is the fastest option and which is what the function does by default. However, you may wish to get the full complex array (for simpler spectrum analysis etc.). Pass the flag to tell the function to produce full-size complex output array.
+            * **DFT_COMPLEX_OUTPUT** Perform a forward transformation of 1D or 2D real array. The result, though being a complex array, has complex-conjugate symmetry ( *CCS* ).?? See the description below for details. Such an array can be packed into a real array of the same size as input, which is the fastest option and which is what the function does by default. However, you may wish to get a full complex array (for simpler spectrum analysis, and so on). Pass the flag to enable the function to produce a full-size complex output array.
 
-            * **DFT_REAL_OUTPUT** then the function performs inverse transformation of 1D or 2D complex array, the result is normally a complex array of the same size. However, if the source array has conjugate-complex symmetry (for example, it is a result of forward transformation with  ``DFT_COMPLEX_OUTPUT``  flag), then the output is real array. While the function itself does not check whether the input is symmetrical or not, you can pass the flag and then the function will assume the symmetry and produce the real output array. Note that when the input is packed real array and inverse transformation is executed, the function treats the input as packed complex-conjugate symmetrical array, so the output will also be real array
+            * **DFT_REAL_OUTPUT** Perform an inverse transformation of 1D or 2D complex array. The result is normally a complex array of the same size. However, if the source array has conjugate-complex symmetry (for example, it is a result of forward transformation with  ``DFT_COMPLEX_OUTPUT``  flag), the output is a real array. While the function itself does not check whether the input is symmetrical or not, you can pass the flag and then the function will assume the symmetry and produce the real output array. Note that when the input is packed into a real array and inverse transformation is executed, the function treats the input as a packed complex-conjugate symmetrical array. So, the output will also be a real array.
 
-    :param nonzeroRows: When the parameter  :math:`\ne 0` , the function assumes that only the first  ``nonzeroRows``  rows of the input array ( ``DFT_INVERSE``  is not set) or only the first  ``nonzeroRows``  of the output array ( ``DFT_INVERSE``  is set) contain non-zeros, thus the function can handle the rest of the rows more efficiently and thus save some time. This technique is very useful for computing array cross-correlation or convolution using DFT
+    :param nonzeroRows: When the parameter  :math:`\ne 0` , the function assumes that only the first  ``nonzeroRows``  rows of the input array ( ``DFT_INVERSE``  is not set) or only the first  ``nonzeroRows``  of the output array ( ``DFT_INVERSE``  is set) contain non-zeros. Thus, the function can handle the rest of the rows more efficiently and save some time. This technique is very useful for computing array cross-correlation or convolution using DFT.
 
 Forward Fourier transform of 1D vector of N elements:
 
@@ -771,37 +800,37 @@ Inverse Fourier transform of 2D vector of
 
     \begin{array}{l} X'=  \left (F^{(M)} \right )^*  \cdot Y  \cdot \left (F^{(N)} \right )^* \\ X =  \frac{1}{M \cdot N} \cdot X' \end{array}
 
-In the case of real (single-channel) data, the packed format called
+In case of real (single-channel) data, the packed format is called
 *CCS*
-(complex-conjugate-symmetrical) that was borrowed from IPL and used to represent the result of a forward Fourier transform or input for an inverse Fourier transform:
+(complex-conjugate-symmetrical). It was borrowed from IPL and used to represent the result of a forward Fourier transform or input for an inverse Fourier transform:
 
 .. math::
 
     \begin{bmatrix} Re Y_{0,0} & Re Y_{0,1} & Im Y_{0,1} & Re Y_{0,2} & Im Y_{0,2} &  \cdots & Re Y_{0,N/2-1} & Im Y_{0,N/2-1} & Re Y_{0,N/2}  \\ Re Y_{1,0} & Re Y_{1,1} & Im Y_{1,1} & Re Y_{1,2} & Im Y_{1,2} &  \cdots & Re Y_{1,N/2-1} & Im Y_{1,N/2-1} & Re Y_{1,N/2}  \\ Im Y_{1,0} & Re Y_{2,1} & Im Y_{2,1} & Re Y_{2,2} & Im Y_{2,2} &  \cdots & Re Y_{2,N/2-1} & Im Y_{2,N/2-1} & Im Y_{1,N/2}  \\ \hdotsfor{9} \\ Re Y_{M/2-1,0} &  Re Y_{M-3,1}  & Im Y_{M-3,1} &  \hdotsfor{3} & Re Y_{M-3,N/2-1} & Im Y_{M-3,N/2-1}& Re Y_{M/2-1,N/2}  \\ Im Y_{M/2-1,0} &  Re Y_{M-2,1}  & Im Y_{M-2,1} &  \hdotsfor{3} & Re Y_{M-2,N/2-1} & Im Y_{M-2,N/2-1}& Im Y_{M/2-1,N/2}  \\ Re Y_{M/2,0}  &  Re Y_{M-1,1} &  Im Y_{M-1,1} &  \hdotsfor{3} & Re Y_{M-1,N/2-1} & Im Y_{M-1,N/2-1}& Re Y_{M/2,N/2} \end{bmatrix}
 
-in the case of 1D transform of real vector, the output will look as the first row of the above matrix.
+In case of 1D transform of a real vector, the output looks like the first row of the matrix above.
 
-So, the function chooses the operation mode depending on the flags and size of the input array:
+So, the function chooses an operation mode depending on the flags and size of the input array:
 
- * if ``DFT_ROWS`` is set or the input array has single row or single column then the function performs 1D forward or inverse transform (of each row of a matrix when ``DFT_ROWS`` is set, otherwise it will be 2D transform.
+ * If ``DFT_ROWS`` is set or the input array has a single row or single column, the function performs a 1D forward or inverse transform of each row of a matrix when ``DFT_ROWS`` is set. Otherwise, it performs a 2D transform.
 
- * if input array is real and ``DFT_INVERSE`` is not set, the function does forward 1D or 2D transform:
+ * If the input array is real and ``DFT_INVERSE`` is not set, the function performs a forward 1D or 2D transform:
 
-    * when ``DFT_COMPLEX_OUTPUT`` is set then the output will be complex matrix of the same size as input.
+    * When ``DFT_COMPLEX_OUTPUT`` is set, the output is a complex matrix of the same size as input.
 
-    * otherwise the output will be a real matrix of the same size as input. in the case of 2D transform it will use the packed format as shown above; in the case of single 1D transform it will look as the first row of the above matrix; in the case of multiple 1D transforms (when using ``DCT_ROWS``         flag) each row of the output matrix will look like the first row of the above matrix.
+    * When ``DFT_COMPLEX_OUTPUT`` is not set, the output is a real matrix of the same size as input. In case of 2D transform, it uses the packed format as shown above. In case of a single 1D transform, it looks like the first row of the matrix above. In case of multiple 1D transforms (when using the ``DCT_ROWS``         flag), each row of the output matrix looks like the first row of the matrix above.
 
- * otherwise, if the input array is complex and either ``DFT_INVERSE``     or ``DFT_REAL_OUTPUT``     are not set then the output will be a complex array of the same size as input and the function will perform the forward or inverse 1D or 2D transform of the whole input array or each row of the input array independently, depending on the flags ``DFT_INVERSE`` and ``DFT_ROWS``.
+ * If the input array is complex and either ``DFT_INVERSE``     or ``DFT_REAL_OUTPUT``     are not set, the output is a complex array of the same size as input. The function performs a forward or inverse 1D or 2D transform of the whole input array or each row of the input array independently, depending on the flags ``DFT_INVERSE`` and ``DFT_ROWS``.
 
- * otherwise, i.e. when ``DFT_INVERSE`` is set, the input array is real, or it is complex but ``DFT_REAL_OUTPUT``     is set, the output will be a real array of the same size as input, and the function will perform 1D or 2D inverse transformation of the whole input array or each individual row, depending on the flags ``DFT_INVERSE`` and ``DFT_ROWS``.
+ * When ``DFT_INVERSE`` is set and the input array is real, or it is complex but ``DFT_REAL_OUTPUT``     is set, the output is a real array of the same size as input. The function performs a 1D or 2D inverse transformation of the whole input array or each individual row, depending on the flags ``DFT_INVERSE`` and ``DFT_ROWS``.
 
-The scaling is done after the transformation if ``DFT_SCALE`` is set.
+If ``DFT_SCALE`` is set, the scaling is done after the transformation.
 
 Unlike
-:func:`dct` , the function supports arrays of arbitrary size, but only those arrays are processed efficiently, which sizes can be factorized in a product of small prime numbers (2, 3 and 5 in the current implementation). Such an efficient DFT size can be computed using
+:func:`dct` , the function supports arrays of arbitrary size. But only those arrays are processed efficiently, whose sizes can be factorized in a product of small prime numbers (2, 3, and 5 in the current implementation). Such an efficient DFT size can be computed using the
 :func:`getOptimalDFTSize` method.
 
-Here is the sample on how to compute DFT-based convolution of two 2D real arrays: ::
+Here is a sample illustrating how to compute a DFT-based convolution of two 2D real arrays: ::
 
     void convolveDFT(const Mat& A, const Mat& B, Mat& C)
     {
@@ -833,7 +862,7 @@ Here is the sample on how to compute DFT-based convolution of two 2D real arrays
 
         // transform the product back from the frequency domain.
         // Even though all the result rows will be non-zero,
-        // we need only the first C.rows of them, and thus we
+        // you need only the first C.rows of them, and thus you
         // pass nonzeroRows == C.rows
         dft(tempA, tempA, DFT_INVERSE + DFT_SCALE, C.rows);
 
@@ -844,20 +873,28 @@ Here is the sample on how to compute DFT-based convolution of two 2D real arrays
     }
 
 
-What can be optimized in the above sample?
+To optimize this sample, consider the following approaches:
 
 *
-    since we passed :math:`\texttt{nonzeroRows} \ne 0`     to the forward transform calls and since we copied ``A``     / ``B``     to the top-left corners of ``tempA``     / ``tempB``     , respectively, it's not necessary to clear the whole ``tempA``     and ``tempB``     ; it is only necessary to clear the ``tempA.cols - A.cols``     ( ``tempB.cols - B.cols``     ) rightmost columns of the matrices.
+    Since :math:`\texttt{nonzeroRows} \ne 0`     is passed to the forward transform calls and since  ``A``     / ``B``     is copied to the top-left corners of ``tempA``     / ``tempB``     , respectively, it is not necessary to clear the whole ``tempA``     and ``tempB``     . It is only necessary to clear the ``tempA.cols - A.cols``     ( ``tempB.cols - B.cols``     ) rightmost columns of the matrices.
 
-* this DFT-based convolution does not have to be applied to the whole big arrays, especially if ``B``     is significantly smaller than ``A``     or vice versa. Instead, we can compute convolution by parts. For that we need to split the destination array ``C``     into multiple tiles and for each tile estimate, which parts of ``A``     and ``B``     are required to compute convolution in this tile. If the tiles in ``C``     are too small, the speed will decrease a lot, because of repeated work - in the ultimate case, when each tile in ``C``     is a single pixel, the algorithm becomes equivalent to the naive convolution algorithm. If the tiles are too big, the temporary arrays ``tempA``     and ``tempB``     become too big and there is also slowdown because of bad cache locality. So there is optimal tile size somewhere in the middle.
+* This DFT-based convolution does not have to be applied to the whole big arrays, especially if ``B``     is significantly smaller than ``A``     or vice versa. Instead, you can compute convolution by parts. To do this, you need to split the destination array ``C``     into multiple tiles. For each tile, estimate which parts of ``A``     and ``B``     are required to compute convolution in this tile. If the tiles in ``C``     are too small, the speed will decrease a lot because of repeated work. In the ultimate case, when each tile in ``C``     is a single pixel, the algorithm becomes equivalent to the naive convolution algorithm. If the tiles are too big, the temporary arrays ``tempA``     and ``tempB``     become too big and there is also a slowdown because of bad cache locality. So, there is an optimal tile size somewhere in the middle.
 
 *
-    if the convolution is done by parts, since different tiles in ``C``     can be computed in parallel, the loop can be threaded.
+    If different tiles in ``C``     can be computed in parallel and, thus, the convolution is done by parts, the loop can be threaded.
 
-All of the above improvements have been implemented in :func:`matchTemplate` and :func:`filter2D` , therefore, by using them, you can get even better performance than with the above theoretically optimal implementation (though, those two functions actually compute cross-correlation, not convolution, so you will need to "flip" the kernel or the image around the center using :func:`flip` ).
+All of the above improvements have been implemented in :func:`matchTemplate` and :func:`filter2D` . Therefore, by using them, you can get the performance even better than with the above theoretically optimal implementation. Though, those two functions actually compute cross-correlation, not convolution, so you need to "flip" the kernel or the image around the center using :func:`flip` .
 
 See Also:
-:func:`dct`,:func:`getOptimalDFTSize`,:func:`mulSpectrums`,:func:`filter2D`,:func:`matchTemplate`,:func:`flip`,:func:`cartToPolar`,:func:`magnitude`,:func:`phase`
+:func:`dct`,
+:func:`getOptimalDFTSize`,
+:func:`mulSpectrums`,
+:func:`filter2D`,
+:func:`matchTemplate`,
+:func:`flip`,
+:func:`cartToPolar`,
+:func:`magnitude`,
+:func:`phase`
 
 .. index:: divide
 
@@ -875,13 +912,13 @@ divide
 
     Performs per-element division of two arrays or a scalar by an array.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array; should have the same size and same type as  ``src1``
+    :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param scale: Scale factor
+    :param scale: Scalar factor.
 
-    :param dst: The destination array; will have the same size and same type as  ``src2``
+    :param dst: The destination array of the same size and type as  ``src2`` .
     
 The functions ``divide`` divide one array by another:
 
@@ -895,10 +932,13 @@ or a scalar by array, when there is no ``src1`` :
 
     \texttt{dst(I) = saturate(scale/src2(I))}
 
-The result will have the same type as ``src1`` . When ``src2(I)=0``,``dst(I)=0`` too.
+The result has the same type as ``src1`` . When ``src2(I)=0``,``dst(I)=0`` too.
 
 See Also:
-:func:`multiply`,:func:`add`,:func:`subtract`,:ref:`MatrixExpressions`
+:func:`multiply`,
+:func:`add`,
+:func:`subtract`,
+:ref:`MatrixExpressions`
 
 .. index:: determinant
 
@@ -909,20 +949,24 @@ determinant
 
 .. c:function:: double determinant(const Mat& mtx)
 
-    Returns determinant of a square floating-point matrix.
+    Returns the determinant of a square floating-point matrix.
 
-    :param mtx: The input matrix; must have  ``CV_32FC1``  or  ``CV_64FC1``  type and square size
+    :param mtx: The input matrix that must have  ``CV_32FC1``  or  ``CV_64FC1``  type and square size.
 
-The function ``determinant`` computes and returns determinant of the specified matrix. For small matrices ( ``mtx.cols=mtx.rows<=3`` )
-the direct method is used; for larger matrices the function uses LU factorization.
+The function ``determinant`` computes and returns the determinant of the specified matrix. For small matrices ( ``mtx.cols=mtx.rows<=3`` ),
+the direct method is used. For larger matrices, the function uses LU factorization.
 
-For symmetric positive-determined matrices, it is also possible to compute
+For symmetric positively-determined matrices, it is also possible to compute
 :func:`SVD` :
 :math:`\texttt{mtx}=U \cdot W \cdot V^T` and then calculate the determinant as a product of the diagonal elements of
 :math:`W` .
 
 See Also:
-:func:`SVD`,:func:`trace`,:func:`invert`,:func:`solve`,:ref:`MatrixExpressions`
+:func:`SVD`,
+:func:`trace`,
+:func:`invert`,
+:func:`solve`,
+:ref:`MatrixExpressions`
 
 .. index:: eigen
 
@@ -937,31 +981,30 @@ eigen
 
     Computes eigenvalues and eigenvectors of a symmetric matrix.
 
-    :param src: The input matrix; must have  ``CV_32FC1``  or  ``CV_64FC1``  type, square size and be symmetric:  :math:`\texttt{src}^T=\texttt{src}`
+    :param src: The input matrix that must have  ``CV_32FC1``  or  ``CV_64FC1``  type, square size and be symmetric:  :math:`\texttt{src}^T=\texttt{src}`
     
-    :param eigenvalues: The output vector of eigenvalues of the same type as  ``src`` ; The eigenvalues are stored in the descending order.
+    :param eigenvalues: The output vector of eigenvalues of the same type as  ``src`` . The eigenvalues are stored in the descending order.
 
-    :param eigenvectors: The output matrix of eigenvectors; It will have the same size and the same type as  ``src`` ; The eigenvectors are stored as subsequent matrix rows, in the same order as the corresponding eigenvalues
+    :param eigenvectors: The output matrix of eigenvectors. It has the same size and type as  ``src`` . The eigenvectors are stored as subsequent matrix rows, in the same order as the corresponding eigenvalues.
 
-    :param lowindex: Optional index of largest eigenvalue/-vector to calculate. (See below.)
+    :param lowindex: Optional index of largest eigenvalue/-vector to calculate. See below for more details.
 
-    :param highindex: Optional index of smallest eigenvalue/-vector to calculate. (See below.)
+    :param highindex: Optional index of smallest eigenvalue/-vector to calculate. See below for more details.
 
-The functions ``eigen`` compute just eigenvalues, or eigenvalues and eigenvectors of symmetric matrix ``src`` : ::
+The functions ``eigen`` compute just eigenvalues, or eigenvalues and eigenvectors of the symmetric matrix ``src`` : ::
 
     src*eigenvectors(i,:)' = eigenvalues(i)*eigenvectors(i,:)' (in MATLAB notation)
 
 
-If either low- or highindex is supplied the other is required, too.
-Indexing is 0-based. Example: To calculate the largest eigenvector/-value set
-lowindex = highindex = 0.
-For legacy reasons this function always returns a square matrix the same size
-as the source matrix with eigenvectors and a vector the length of the source
-matrix with eigenvalues. The selected eigenvectors/-values are always in the
-first highindex - lowindex + 1 rows.
+If either ``low-`` or ``highindex`` is supplied, the other one is required, too.
+Indexing is 0-based. For example, to calculate the largest eigenvector/-value set,
+``lowindex = highindex = 0`` .
+For legacy reasons, this function always returns a square matrix of the same size as the source matrix with eigenvectors and a vector of the length of the source matrix with eigenvalues. The selected eigenvectors/-values are always in the first ``highindex`` - ``lowindex`` + 1 rows.
 
 See Also:
-:func:`SVD`,:func:`completeSymm`,:func:`PCA`
+:func:`SVD`,
+:func:`completeSymm`,
+:func:`PCA`
 
 .. index:: exp
 
@@ -976,9 +1019,9 @@ exp
 
     Calculates the exponent of every array element.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param dst: The destination array; will have the same size and same type as  ``src``
+    :param dst: Destination array of the same size and same type as  ``src`` .
 
 The function ``exp`` calculates the exponent of every element of the input array:
 
@@ -992,7 +1035,13 @@ The maximum relative error is about
 :math:`\pm \infty` ) are not handled.
 
 See Also:
-:func:`log`,:func:`cartToPolar`,:func:`polarToCart`,:func:`phase`,:func:`pow`,:func:`sqrt`,:func:`magnitude`
+:func:`log`,
+:func:`cartToPolar`,
+:func:`polarToCart`,
+:func:`phase`,
+:func:`pow`,
+:func:`sqrt`,
+:func:`magnitude`
 
 .. index:: extractImageCOI
 
@@ -1003,20 +1052,27 @@ extractImageCOI
 
 .. c:function:: void extractImageCOI(const CvArr* src, Mat& dst, int coi=-1)
 
-    Extract the selected image channel
+    Extracts the selected image channel.
 
-    :param src: The source array. It should be a pointer to  ``CvMat``  or  ``IplImage``
+    :param src: Source array. It should be a pointer to  ``CvMat``  or  ``IplImage`` .
     
-    :param dst: The destination array; will have single-channel, and the same size and the same depth as  ``src``
+    :param dst: Destination array with a single channel and the same size and depth as  ``src`` .
     
-    :param coi: If the parameter is  ``>=0`` , it specifies the channel to extract; If it is  ``<0`` , ``src``  must be a pointer to  ``IplImage``  with valid COI set - then the selected COI is extracted.
+    :param coi: If the parameter is  ``>=0`` , it specifies the channel to extract. If it is  ``<0`` and ``src``  is a pointer to  ``IplImage``  with a  valid COI set, the selected COI is extracted.
 
-The function ``extractImageCOI`` is used to extract image COI from an old-style array and put the result to the new-style C++ matrix. As usual, the destination matrix is reallocated using ``Mat::create`` if needed.
+The function ``extractImageCOI`` is used to extract an image COI from an old-style array and put the result to the new-style C++ matrix. As usual, the destination matrix is reallocated using ``Mat::create`` if needed.
 
 To extract a channel from a new-style matrix, use
 :func:`mixChannels` or
-:func:`split` See Also:
-:func:`mixChannels`,:func:`split`,:func:`merge`,:func:`cvarrToMat`,:func:`cvSetImageCOI`,:func:`cvGetImageCOI`
+:func:`split` .
+
+See Also:
+:func:`mixChannels`,
+:func:`split`,
+:func:`merge`,
+:func:`cvarrToMat`,
+:func:`cvSetImageCOI`,
+:func:`cvGetImageCOI`
 
 .. index:: fastAtan2
 
@@ -1027,14 +1083,13 @@ fastAtan2
 
 .. c:function:: float fastAtan2(float y, float x)
 
-    Calculates the angle of a 2D vector in degrees
+    Calculates the angle of a 2D vector in degrees.
 
-    :param x: x-coordinate of the vector
+    :param x: x-coordinate of the vector.
 
-    :param y: y-coordinate of the vector
+    :param y: y-coordinate of the vector.
 
-The function ``fastAtan2`` calculates the full-range angle of an input 2D vector. The angle is
-measured in degrees and varies from
+The function ``fastAtan2`` calculates the full-range angle of an input 2D vector. The angle is measured in degrees and varies from
 :math:`0^\circ` to
 :math:`360^\circ` . The accuracy is about
 :math:`0.3^\circ` .
@@ -1045,13 +1100,13 @@ flip
 --------
 .. c:function:: void flip(const Mat& src, Mat& dst, int flipCode)
 
-    Flips a 2D array around vertical, horizontal or both axes.
+    Flips a 2D array around vertical, horizontal, or both axes.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param dst: The destination array; will have the same size and same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
     
-    :param flipCode: Specifies how to flip the array: 0 means flipping around the x-axis, positive (e.g., 1) means flipping around y-axis, and negative (e.g., -1) means flipping around both axes. See Also the discussion below for the formulas.
+    :param flipCode: Flag to specify how to flip the array. 0 means flipping around the x-axis. Positive value (for example, 1) means flipping around y-axis. Negative value (for example, -1) means flipping around both axes. See the discussion below for the formulas.
 
 The function ``flip`` flips the array in one of three different ways (row and column indices are 0-based):
 
@@ -1061,26 +1116,28 @@ The function ``flip`` flips the array in one of three different ways (row and co
     { \texttt{src} _{i, \texttt{src.cols} -j-1}}{if  \texttt{flipCode} > 0}
     { \texttt{src} _{ \texttt{src.rows} -i-1, \texttt{src.cols} -j-1}}{if  \texttt{flipCode} < 0}
 
-The example scenarios of function use are:
+The example scenarios of using the function are the following:
 
 *
-    vertical flipping of the image (
-    :math:`\texttt{flipCode} = 0`     ) to switch between top-left and bottom-left image origin, which is a typical operation in video processing in Windows.
+    Vertical flipping of the image (
+    :math:`\texttt{flipCode} = 0`     ) to switch between top-left and bottom-left image origin. This is a typical operation in video processing on Microsoft Windows* OS.
 
 *
-    horizontal flipping of the image with subsequent horizontal shift and absolute difference calculation to check for a vertical-axis symmetry (
-    :math:`\texttt{flipCode} > 0`     )
+    Horizontal flipping of the image with the subsequent horizontal shift and absolute difference calculation to check for a vertical-axis symmetry (
+    :math:`\texttt{flipCode} > 0`     ).
 
 *
-    simultaneous horizontal and vertical flipping of the image with subsequent shift and absolute difference calculation to check for a central symmetry (
-    :math:`\texttt{flipCode} < 0`     )
+    Simultaneous horizontal and vertical flipping of the image with the subsequent shift and absolute difference calculation to check for a central symmetry (
+    :math:`\texttt{flipCode} < 0`     ).
 
 *
-    reversing the order of 1d point arrays (
+    Reversing the order of 1D point arrays (
     :math:`\texttt{flipCode} > 0`     or
-    :math:`\texttt{flipCode} = 0`     )
+    :math:`\texttt{flipCode} = 0`     ).
 
-See Also: :func:`transpose`,:func:`repeat`,:func:`completeSymm`
+See Also: :func:`transpose`,
+:func:`repeat`,
+:func:`completeSymm`
 
 .. index:: gemm
 
@@ -1093,17 +1150,17 @@ gemm
 
     Performs generalized matrix multiplication.
 
-    :param src1: The first multiplied input matrix; should have  ``CV_32FC1`` , ``CV_64FC1`` , ``CV_32FC2``  or  ``CV_64FC2``  type
+    :param src1: The first multiplied input matrix that should have  ``CV_32FC1`` , ``CV_64FC1`` , ``CV_32FC2`` , or  ``CV_64FC2``  type.
 
-    :param src2: The second multiplied input matrix; should have the same type as  ``src1``
+    :param src2: The second multiplied input matrix of the same type as  ``src1`` .
     
-    :param alpha: The weight of the matrix product
+    :param alpha: Weight of the matrix product.
 
-    :param src3: The third optional delta matrix added to the matrix product; should have the same type as  ``src1``  and  ``src2``
+    :param src3: The third optional delta matrix added to the matrix product. It should have the same type as  ``src1``  and  ``src2`` .
     
-    :param beta: The weight of  ``src3``
+    :param beta: Weight of  ``src3`` .
     
-    :param dst: The destination matrix; It will have the proper size and the same type as input matrices
+    :param dst: The destination matrix. It has the proper size and the same type as input matrices.
 
     :param flags: Operation flags:
 
@@ -1117,13 +1174,15 @@ The function performs generalized matrix multiplication and similar to the corre
 
     \texttt{dst} =  \texttt{alpha} \cdot \texttt{src1} ^T  \cdot \texttt{src2} +  \texttt{beta} \cdot \texttt{src3} ^T
 
-The function can be replaced with a matrix expression, e.g. the above call can be replaced with: ::
+The function can be replaced with a matrix expression. For example, the above call can be replaced with: ::
 
     dst = alpha*src1.t()*src2 + beta*src3.t();
 
 
 See Also:
-:func:`mulTransposed`,:func:`transform`,:ref:`MatrixExpressions`
+:func:`mulTransposed`,
+:func:`transform`,
+:ref:`MatrixExpressions`
 
 .. index:: getConvertElem
 
@@ -1140,26 +1199,28 @@ getConvertElem
 
 .. c:function:: typedef void (*ConvertScaleData)(const void* from, void* to, int cn, double alpha, double beta)
 
-    Returns conversion function for a single pixel
+    Returns a conversion function for a single pixel.
 
-    :param fromType: The source pixel type
+    :param fromType: The source pixel type.
 
-    :param toType: The destination pixel type
+    :param toType: The destination pixel type.
 
-    :param from: Callback parameter: pointer to the input pixel
+    :param from: Callback parameter: pointer to the input pixel.
 
     :param to: Callback parameter: pointer to the output pixel
 
-    :param cn: Callback parameter: the number of channels; can be arbitrary, 1, 100, 100000, ...
+    :param cn: Callback parameter: the number of channels. It can be arbitrary, 1, 100, 100000, ...
 
-    :param alpha: ConvertScaleData callback optional parameter: the scale factor
+    :param alpha: ``ConvertScaleData`` callback optional parameter: the scale factor.
 
-    :param beta: ConvertScaleData callback optional parameter: the delta or offset
+    :param beta: ``ConvertScaleData`` callback optional parameter: the delta or offset.
 
 The functions ``getConvertElem`` and ``getConvertScaleElem`` return pointers to the functions for converting individual pixels from one type to another. While the main function purpose is to convert single pixels (actually, for converting sparse matrices from one type to another), you can use them to convert the whole row of a dense matrix or the whole matrix at once, by setting ``cn = matrix.cols*matrix.rows*matrix.channels()`` if the matrix data is continuous.
 
 See Also:
-:func:`Mat::convertTo`,:func:`MatND::convertTo`,:func:`SparseMat::convertTo`
+:func:`Mat::convertTo`,
+:func:`MatND::convertTo`,
+:func:`SparseMat::convertTo`
 
 .. index:: getOptimalDFTSize
 
@@ -1170,14 +1231,14 @@ getOptimalDFTSize
 
 .. c:function:: int getOptimalDFTSize(int vecsize)
 
-    Returns optimal DFT size for a given vector size.
+    Returns the optimal DFT size for a given vector size.
 
-    :param vecsize: Vector size
+    :param vecsize: Vector size.
 
-DFT performance is not a monotonic function of a vector size, therefore, when you compute convolution of two arrays or do a spectral analysis of array, it usually makes sense to pad the input data with zeros to get a bit larger array that can be transformed much faster than the original one.
-Arrays, which size is a power-of-two (2, 4, 8, 16, 32, ...) are the fastest to process, though, the arrays, which size is a product of 2's, 3's and 5's (e.g. 300 = 5*5*3*2*2), are also processed quite efficiently.
+DFT performance is not a monotonic function of a vector size. Therefore, when you compute convolution of two arrays or perform the spectral analysis of an array, it usually makes sense to pad the input data with zeros to get a bit larger array that can be transformed much faster than the original one.
+Arrays whose size is a power-of-two (2, 4, 8, 16, 32, ...) are the fastest to process. Though, the arrays whose size is a product of 2's, 3's, and 5's (for example, 300 = 5*5*3*2*2) are also processed quite efficiently.
 
-The function ``getOptimalDFTSize`` returns the minimum number ``N`` that is greater than or equal to ``vecsize`` , such that the DFT
+The function ``getOptimalDFTSize`` returns the minimum number ``N`` that is greater than or equal to ``vecsize``  so that the DFT
 of a vector of size ``N`` can be computed efficiently. In the current implementation
 :math:`N=2^p \times 3^q \times 5^r` , for some
 :math:`p`,:math:`q`,:math:`r` .
@@ -1187,7 +1248,11 @@ The function returns a negative number if ``vecsize`` is too large (very close t
 While the function cannot be used directly to estimate the optimal vector size for DCT transform (since the current DCT implementation supports only even-size vectors), it can be easily computed as ``getOptimalDFTSize((vecsize+1)/2)*2`` .
 
 See Also:
-:func:`dft`,:func:`dct`,:func:`idft`,:func:`idct`,:func:`mulSpectrums`
+:func:`dft`,
+:func:`dct`,
+:func:`idft`,
+:func:`idct`,
+:func:`mulSpectrums`
 
 .. index:: idct
 
@@ -1198,17 +1263,20 @@ idct
 
 .. c:function:: void idct(const Mat& src, Mat& dst, int flags=0)
 
-    Computes inverse Discrete Cosine Transform of a 1D or 2D array
+    Computes the inverse Discrete Cosine Transform of a 1D or 2D array.
 
-    :param src: The source floating-point single-channel array
+    :param src: The source floating-point single-channel array.
 
-    :param dst: The destination array. Will have the same size and same type as  ``src``
+    :param dst: The destination array of the same size and type as  ``src`` .
     
     :param flags: The operation flags.
     
 ``idct(src, dst, flags)`` is equivalent to ``dct(src, dst, flags | DCT_INVERSE)``.
 
-See Also: :func:`dct`,:func:`dft`,:func:`idft`,:func:`getOptimalDFTSize`
+See Also: :func:`dct`,
+:func:`dft`,
+:func:`idft`,
+:func:`getOptimalDFTSize`
 
 .. index:: idft
 
@@ -1219,23 +1287,29 @@ idft
 
 .. c:function:: void idft(const Mat& src, Mat& dst, int flags=0, int outputRows=0)
 
-    Computes inverse Discrete Fourier Transform of a 1D or 2D array
+    Computes the inverse Discrete Fourier Transform of a 1D or 2D array.
 
-    :param src: The source floating-point real or complex array
+    :param src: The source floating-point real or complex array.
 
-    :param dst: The destination array, which size and type depends on the  ``flags``
+    :param dst: The destination array whose size and type depend on the  ``flags`` .
     
-    :param flags: The operation flags. See  :func:`dft`
+    :param flags: The operation flags. See  :func:`dft` .
     
-    :param nonzeroRows: The number of  ``dst``  rows to compute. The rest of the rows will have undefined content. See the convolution sample in  :func:`dft`  description
+    :param nonzeroRows: The number of  ``dst``  rows to compute. The rest of the rows have undefined content. See the convolution sample in  :func:`dft`  description.
     
 ``idft(src, dst, flags)`` is equivalent to ``dct(src, dst, flags | DFT_INVERSE)`` .
 
 See :func:`dft` for details.
-Note, that none of ``dft`` and ``idft`` scale the result by default.
+**Note**:
+
+None of ``dft`` and ``idft`` scale the result by default.
 Thus, you should pass ``DFT_SCALE`` to one of ``dft`` or ``idft`` explicitly to make these transforms mutually inverse.
 
-See Also: :func:`dft`,:func:`dct`,:func:`idct`,:func:`mulSpectrums`,:func:`getOptimalDFTSize`
+See Also: :func:`dft`,
+:func:`dct`,
+:func:`idct`,
+:func:`mulSpectrums`,
+:func:`getOptimalDFTSize`
 
 .. index:: inRange
 
@@ -1254,27 +1328,31 @@ inRange
 
     Checks if array elements lie between the elements of two other arrays.
 
-    :param src: The first source array
+    :param src: The first source array.
 
-    :param lowerb: The inclusive lower boundary array of the same size and type as  ``src``
+    :param lowerb: Inclusive lower boundary array of the same size and type as  ``src`` .
     
-    :param upperb: The exclusive upper boundary array of the same size and type as  ``src``
+    :param upperb: Exclusive upper boundary array of the same size and type as  ``src`` .
     
-    :param dst: The destination array, will have the same size as  ``src``  and  ``CV_8U``  type
+    :param dst: Destination array of the same size as  ``src``  and  ``CV_8U``  type.
 
-The functions ``inRange`` do the range check for every element of the input array:
+The functions ``inRange`` check the range as follows:
+
+* for every element of the input array:
 
 .. math::
 
     \texttt{dst} (I)= \texttt{lowerb} (I)_0  \leq \texttt{src} (I)_0 <  \texttt{upperb} (I)_0
 
-for single-channel arrays,
+* for single-channel arrays:
 
 .. math::
 
     \texttt{dst} (I)= \texttt{lowerb} (I)_0  \leq \texttt{src} (I)_0 <  \texttt{upperb} (I)_0  \land \texttt{lowerb} (I)_1  \leq \texttt{src} (I)_1 <  \texttt{upperb} (I)_1
 
-for two-channel arrays and so forth. ``dst`` (I) is set to 255 (all ``1`` -bits) if ``src`` (I) is within the specified range and 0 otherwise.
+* for two-channel arrays and so forth. 
+
+``dst`` (I) is set to 255 (all ``1`` -bits) if ``src`` (I) is within the specified range and 0 otherwise.
 
 .. index:: invert
 
@@ -1285,32 +1363,33 @@ invert
 
 .. c:function:: double invert(const Mat& src, Mat& dst, int method=DECOMP_LU)
 
-    Finds the inverse or pseudo-inverse of a matrix
+    Finds the inverse or pseudo-inverse of a matrix.
 
-    :param src: The source floating-point  :math:`M \times N`  matrix
+    :param src: The source floating-point  :math:`M \times N`  matrix.
 
-    :param dst: The destination matrix; will have  :math:`N \times M`  size and the same type as  ``src``
+    :param dst: The destination matrix of  :math:`N \times M`  size and the same type as  ``src`` .
     
     :param flags: The inversion method :
 
-            * **DECOMP_LU** Gaussian elimination with optimal pivot element chosen
+            * **DECOMP_LU** Gaussian elimination with the optimal pivot element chosen.
 
-            * **DECOMP_SVD** Singular value decomposition (SVD) method
+            * **DECOMP_SVD** Singular value decomposition (SVD) method.
 
-            * **DECOMP_CHOLESKY** Cholesky decomposion. The matrix must be symmetrical and positively defined
+            * **DECOMP_CHOLESKY** Cholesky decomposion. The matrix must be symmetrical and positively defined.
 
-The function ``invert`` inverts matrix ``src`` and stores the result in ``dst`` .
-When the matrix ``src`` is singular or non-square, the function computes the pseudo-inverse matrix, i.e. the matrix ``dst`` , such that
+The function ``invert`` inverts the matrix ``src`` and stores the result in ``dst`` .
+When the matrix ``src`` is singular or non-square, the function computes the pseudo-inverse matrix (the ``dst`` matrix) so that
 :math:`\|\texttt{src} \cdot \texttt{dst} - I\|` is minimal.
 
-In the case of ``DECOMP_LU`` method, the function returns the ``src`` determinant ( ``src`` must be square). If it is 0, the matrix is not inverted and ``dst`` is filled with zeros.
+In case of the ``DECOMP_LU`` method, the function returns the ``src`` determinant ( ``src`` must be square). If it is 0, the matrix is not inverted and ``dst`` is filled with zeros.
 
-In the case of ``DECOMP_SVD`` method, the function returns the inversed condition number of ``src`` (the ratio of the smallest singular value to the largest singular value) and 0 if ``src`` is singular. The SVD method calculates a pseudo-inverse matrix if ``src`` is singular.
+In case of the ``DECOMP_SVD`` method, the function returns the inverse condition number of ``src`` (the ratio of the smallest singular value to the largest singular value) and 0 if ``src`` is singular. The SVD method calculates a pseudo-inverse matrix if ``src`` is singular.
 
-Similarly to ``DECOMP_LU`` , the method ``DECOMP_CHOLESKY`` works only with non-singular square matrices. In this case the function stores the inverted matrix in ``dst`` and returns non-zero, otherwise it returns 0.
+Similarly to ``DECOMP_LU`` , the method ``DECOMP_CHOLESKY`` works only with non-singular square matrices. In this case, the function stores the inverted matrix in ``dst`` and returns non-zero. Otherwise, it returns 0.
 
 See Also:
-:func:`solve`,:func:`SVD`
+:func:`solve`,
+:func:`SVD`
 
 .. index:: log
 
@@ -1325,9 +1404,9 @@ log
 
     Calculates the natural logarithm of every array element.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param dst: The destination array; will have the same size and same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
     
 The function ``log`` calculates the natural logarithm of the absolute value of every element of the input array:
 
@@ -1335,14 +1414,20 @@ The function ``log`` calculates the natural logarithm of the absolute value of e
 
     \texttt{dst} (I) =  \fork{\log |\texttt{src}(I)|}{if $\texttt{src}(I) \ne 0$ }{\texttt{C}}{otherwise}
 
-Where ``C`` is a large negative number (about -700 in the current implementation).
+where ``C`` is a large negative number (about -700 in the current implementation).
 The maximum relative error is about
 :math:`7 \times 10^{-6}` for single-precision input and less than
 :math:`10^{-10}` for double-precision input. Special values (NaN,
 :math:`\pm \infty` ) are not handled.
 
 See Also:
-:func:`exp`,:func:`cartToPolar`,:func:`polarToCart`,:func:`phase`,:func:`pow`,:func:`sqrt`,:func:`magnitude`
+:func:`exp`,
+:func:`cartToPolar`,
+:func:`polarToCart`,
+:func:`phase`,
+:func:`pow`,
+:func:`sqrt`,
+:func:`magnitude`
 
 .. index:: LUT
 
@@ -1355,11 +1440,11 @@ LUT
 
     Performs a look-up table transform of an array.
 
-    :param src: Source array of 8-bit elements
+    :param src: Source array of 8-bit elements.
 
-    :param lut: Look-up table of 256 elements. In the case of multi-channel source array, the table should either have a single channel (in this case the same table is used for all channels) or the same number of channels as in the source array
+    :param lut: Look-up table of 256 elements. In case of multi-channel source array, the table should either have a single channel (in this case the same table is used for all channels) or the same number of channels as in the source array.
 
-    :param dst: Destination array; will have the same size and the same number of channels as  ``src`` , and the same depth as  ``lut``
+    :param dst: Destination array of the same size and the same number of channels as  ``src`` , and the same depth as  ``lut`` .
     
 The function ``LUT`` fills the destination array with values from the look-up table. Indices of the entries are taken from the source array. That is, the function processes each element of ``src`` as follows:
 
@@ -1385,22 +1470,25 @@ magnitude
 
 .. c:function:: void magnitude(const Mat& x, const Mat& y, Mat& magnitude)
 
-    Calculates magnitude of 2D vectors.
+    Calculates the magnitude of 2D vectors.
 
-    :param x: The floating-point array of x-coordinates of the vectors
+    :param x: Floating-point array of x-coordinates of the vectors.
 
-    :param y: The floating-point array of y-coordinates of the vectors; must have the same size as  ``x``
+    :param y: Floating-point array of y-coordinates of the vectors. It must have the same size as  ``x`` .
     
-    :param dst: The destination array; will have the same size and same type as  ``x``
+    :param dst: Destination array of the same size and type as  ``x`` .
     
-The function ``magnitude`` calculates magnitude of 2D vectors formed from the corresponding elements of ``x`` and ``y`` arrays:
+The function ``magnitude`` calculates the magnitude of 2D vectors formed from the corresponding elements of ``x`` and ``y`` arrays:
 
 .. math::
 
     \texttt{dst} (I) =  \sqrt{\texttt{x}(I)^2 + \texttt{y}(I)^2}
 
 See Also:
-:func:`cartToPolar`,:func:`polarToCart`,:func:`phase`,:func:`sqrt`
+:func:`cartToPolar`,
+:func:`polarToCart`,
+:func:`phase`,
+:func:`sqrt`
 
 .. index:: Mahalanobis
 
@@ -1413,11 +1501,11 @@ Mahalanobis
 
     Calculates the Mahalanobis distance between two vectors.
 
-    :param vec1: The first 1D source vector
+    :param vec1: The first 1D source vector.
 
-    :param vec2: The second 1D source vector
+    :param vec2: The second 1D source vector.
 
-    :param icovar: The inverse covariance matrix
+    :param icovar: Inverse covariance matrix.
 
 The function ``cvMahalonobis`` calculates and returns the weighted distance between two vectors:
 
@@ -1427,7 +1515,7 @@ The function ``cvMahalonobis`` calculates and returns the weighted distance betw
 
 The covariance matrix may be calculated using the
 :func:`calcCovarMatrix` function and then inverted using the
-:func:`invert` function (preferably using DECOMP_SVD method, as the most accurate).
+:func:`invert` function (preferably using the ``DECOMP_SVD`` method, as the most accurate).
 
 .. index:: max
 
@@ -1450,17 +1538,17 @@ max
 
 .. c:function:: void max(const MatND& src1, double value, MatND& dst)
 
-    Calculates per-element maximum of two arrays or array and a scalar
+    Calculates per-element maximum of two arrays or array and a scalar.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array of the same size and type as  ``src1``
+    :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param value: The real scalar value
+    :param value: Real scalar value.
 
-    :param dst: The destination array; will have the same size and type as  ``src1``
+    :param dst: Destination array of the same size and type as  ``src1`` .
     
-The functions ``max`` compute per-element maximum of two arrays:
+The functions ``max`` compute the per-element maximum of two arrays:
 
 .. math::
 
@@ -1475,10 +1563,14 @@ or array and a scalar:
 In the second variant, when the source array is multi-channel, each channel is compared with ``value`` independently.
 
 The first 3 variants of the function listed above are actually a part of
-:ref:`MatrixExpressions` , they return the expression object that can be further transformed, or assigned to a matrix, or passed to a function etc.
+:ref:`MatrixExpressions` . They return an expression object that can be further either transformed/ assigned to a matrix, or passed to a function, and so on.
 
 See Also:
-:func:`min`,:func:`compare`,:func:`inRange`,:func:`minMaxLoc`,:ref:`MatrixExpressions`
+:func:`min`,
+:func:`compare`,
+:func:`inRange`,
+:func:`minMaxLoc`,
+:ref:`MatrixExpressions`
 
 .. index:: mean
 
@@ -1495,11 +1587,11 @@ mean
 
 .. c:function:: Scalar mean(const MatND& mtx, const MatND& mask)
 
-    Calculates average (mean) of array elements
+    Calculates an average (mean) of array elements.
 
-    :param mtx: The source array; it should have 1 to 4 channels (so that the result can be stored in  :func:`Scalar` )
+    :param mtx: Source array that should have from 1 to 4 channels so that the result can be stored in  :func:`Scalar` .
 
-    :param mask: The optional operation mask
+    :param mask: Optional operation mask.
 
 The functions ``mean`` compute mean value ``M`` of array elements, independently for each channel, and return it:
 
@@ -1510,7 +1602,10 @@ The functions ``mean`` compute mean value ``M`` of array elements, independently
 When all the mask elements are 0's, the functions return ``Scalar::all(0)`` .
 
 See Also:
-:func:`countNonZero`,:func:`meanStdDev`,:func:`norm`,:func:`minMaxLoc`
+:func:`countNonZero`,
+:func:`meanStdDev`,
+:func:`norm`,
+:func:`minMaxLoc`
 
 .. index:: meanStdDev
 
@@ -1523,30 +1618,35 @@ meanStdDev
 
 .. c:function:: void meanStdDev(const MatND& mtx, Scalar& mean, Scalar& stddev, const MatND& mask=MatND())
 
-    Calculates mean and standard deviation of array elements
+    Calculates mean and standard deviation of array elements.
 
-    :param mtx: The source array; it should have 1 to 4 channels (so that the results can be stored in  :func:`Scalar` 's)
+    :param mtx: Source array that should have from 1 to 4 channels so that the results can be stored in  :func:`Scalar` 's.
 
-    :param mean: The output parameter: computed mean value
+    :param mean: Output parameter: computed mean value.
 
-    :param stddev: The output parameter: computed standard deviation
+    :param stddev: Output parameter: computed standard deviation.
 
-    :param mask: The optional operation mask
+    :param mask: Optional operation mask.
 
-The functions ``meanStdDev`` compute the mean and the standard deviation ``M`` of array elements, independently for each channel, and return it via the output parameters:
+The functions ``meanStdDev`` compute the mean and the standard deviation ``M`` of array elements independently for each channel and return it via the output parameters:
 
 .. math::
 
     \begin{array}{l} N =  \sum _{I, \texttt{mask} (I)  \ne 0} 1 \\ \texttt{mean} _c =  \frac{\sum_{ I: \; \texttt{mask}(I) \ne 0} \texttt{src} (I)_c}{N} \\ \texttt{stddev} _c =  \sqrt{\sum_{ I: \; \texttt{mask}(I) \ne 0} \left ( \texttt{src} (I)_c -  \texttt{mean} _c \right )^2} \end{array}
 
 When all the mask elements are 0's, the functions return ``mean=stddev=Scalar::all(0)`` .
-Note that the computed standard deviation is only the diagonal of the complete normalized covariance matrix. If the full matrix is needed, you can reshape the multi-channel array
+**Note**:
+The computed standard deviation is only the diagonal of the complete normalized covariance matrix. If the full matrix is needed, you can reshape the multi-channel array
 :math:`M \times N` to the single-channel array
 :math:`M*N \times \texttt{mtx.channels}()` (only possible when the matrix is continuous) and then pass the matrix to
 :func:`calcCovarMatrix` .
 
 See Also:
-:func:`countNonZero`,:func:`mean`,:func:`norm`,:func:`minMaxLoc`,:func:`calcCovarMatrix`
+:func:`countNonZero`,
+:func:`mean`,
+:func:`norm`,
+:func:`minMaxLoc`,
+:func:`calcCovarMatrix`
 
 .. index:: merge
 
@@ -1565,11 +1665,11 @@ merge
 
     Composes a multi-channel array from several single-channel arrays.
 
-    :param mv: The source array or vector of the single-channel matrices to be merged. All the matrices in  ``mv``  must have the same size and the same type
+    :param mv: Source array or vector of the single-channel matrices to be merged. All the matrices in  ``mv``  must have the same size and the same type.
 
-    :param count: The number of source matrices when  ``mv``  is a plain C array; must be greater than zero
+    :param count: Number of source matrices when  ``mv``  is a plain C array. It must be greater than zero.
 
-    :param dst: The destination array; will have the same size and the same depth as  ``mv[0]`` , the number of channels will match the number of source matrices
+    :param dst: Destination array of the same size and the same depth as  ``mv[0]`` . The number of channels matches the number of source matrices.
 
 The functions ``merge`` merge several single-channel arrays (or rather interleave their elements) to make a single multi-channel array.
 
@@ -1578,9 +1678,13 @@ The functions ``merge`` merge several single-channel arrays (or rather interleav
     \texttt{dst} (I)_c =  \texttt{mv} [c](I)
 
 The function
-:func:`split` does the reverse operation and if you need to merge several multi-channel images or shuffle channels in some other advanced way, use
-:func:`mixChannels` See Also:
-:func:`mixChannels`,:func:`split`,:func:`reshape`
+:func:`split` does the reverse operatio. If you need to merge several multi-channel images or shuffle channels in some other advanced way, use
+:func:`mixChannels` .
+
+See Also:
+:func:`mixChannels`,
+:func:`split`,
+:func:`reshape`
 
 .. index:: min
 
@@ -1603,17 +1707,17 @@ min
 
 .. c:function:: void min(const MatND& src1, double value, MatND& dst)
 
-    Calculates per-element minimum of two arrays or array and a scalar
+    Calculates per-element minimum of two arrays or array and a scalar.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array of the same size and type as  ``src1``
+    :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param value: The real scalar value
+    :param value: Real scalar value.
 
-    :param dst: The destination array; will have the same size and type as  ``src1``
+    :param dst: Destination array of the same size and type as  ``src1`` .
     
-The functions ``min`` compute per-element minimum of two arrays:
+The functions ``min`` compute the per-element minimum of two arrays:
 
 .. math::
 
@@ -1627,11 +1731,15 @@ or array and a scalar:
 
 In the second variant, when the source array is multi-channel, each channel is compared with ``value`` independently.
 
-The first 3 variants of the function listed above are actually a part of
-:ref:`MatrixExpressions` , they return the expression object that can be further transformed, or assigned to a matrix, or passed to a function etc.
+The first three variants of the function listed above are actually a part of
+:ref:`MatrixExpressions` . They return the expression object that can be further either transformed/assigned to a matrix, or passed to a function, and so on.
 
 See Also:
-:func:`max`,:func:`compare`,:func:`inRange`,:func:`minMaxLoc`,:ref:`MatrixExpressions`
+:func:`max`,
+:func:`compare`,
+:func:`inRange`,
+:func:`minMaxLoc`,
+:ref:`MatrixExpressions`
 
 .. index:: minMaxLoc
 
@@ -1646,38 +1754,44 @@ minMaxLoc
 
 .. c:function:: void minMaxLoc(const SparseMat& src, double* minVal, double* maxVal, int* minIdx=0, int* maxIdx=0)
 
-    Finds global minimum and maximum in a whole array or sub-array
+    Finds the global minimum and maximum in a whole array or sub-array.
 
-    :param src: The source single-channel array
+    :param src: Source single-channel array.
 
-    :param minVal: Pointer to returned minimum value;  ``NULL``  if not required
+    :param minVal: Pointer to the returned minimum value.  ``NULL`` is used if not required.
 
-    :param maxVal: Pointer to returned maximum value;  ``NULL``  if not required
+    :param maxVal: Pointer to the returned maximum value.  ``NULL`` is used if not required.
 
-    :param minLoc: Pointer to returned minimum location (in 2D case);  ``NULL``  if not required
+    :param minLoc: Pointer to the returned minimum location (in 2D case).  ``NULL`` is used if not required.
 
-    :param maxLoc: Pointer to returned maximum location (in 2D case);  ``NULL``  if not required
+    :param maxLoc: Pointer to the returned maximum location (in 2D case).  ``NULL`` is used if not required.
 
-    :param minIdx: Pointer to returned minimum location (in nD case); ``NULL``  if not required, otherwise must point to an array of  ``src.dims``  elements and the coordinates of minimum element in each dimensions will be stored sequentially there.
+    :param minIdx: Pointer to the returned minimum location (in nD case). ``NULL`` is used if not required. Otherwise, it must point to an array of  ``src.dims``  elements. The coordinates of minimum element in each dimensions will be stored sequentially there.
 
-    :param maxIdx: Pointer to returned maximum location (in nD case);  ``NULL``  if not required
+    :param maxIdx: Pointer to the returned maximum location (in nD case).  ``NULL`` is used if not required.
 
-    :param mask: The optional mask used to select a sub-array
+    :param mask: Optional mask used to select a sub-array.
 
-The functions ``ninMaxLoc`` find minimum and maximum element values
-and their positions. The extremums are searched across the whole array, or,
+The functions ``ninMaxLoc`` find minimum and maximum element values and their positions. The extremums are searched across the whole array or,
 if ``mask`` is not an empty array, in the specified array region.
 
 The functions do not work with multi-channel arrays. If you need to find minimum or maximum elements across all the channels, use
-:func:`reshape` first to reinterpret the array as single-channel. Or you may extract the particular channel using
-:func:`extractImageCOI` or
-:func:`mixChannels` or
+:func:`reshape` first to reinterpret the array as single-channel. Or you may extract the particular channel using either
+:func:`extractImageCOI` or
+:func:`mixChannels` or
 :func:`split` .
 
-in the case of a sparse matrix the minimum is found among non-zero elements only.
+In case of a sparse matrix, the minimum is found among non-zero elements only.
 
 See Also:
-:func:`max`,:func:`min`,:func:`compare`,:func:`inRange`,:func:`extractImageCOI`,:func:`mixChannels`,:func:`split`,:func:`reshape` .
+:func:`max`,
+:func:`min`,
+:func:`compare`,
+:func:`inRange`,
+:func:`extractImageCOI`,
+:func:`mixChannels`,
+:func:`split`,
+:func:`reshape` 
 
 .. index:: mixChannels
 
@@ -1694,18 +1808,17 @@ mixChannels
 
 .. c:function:: void mixChannels(const vector<MatND>& srcv, vector<MatND>& dstv, const int* fromTo, int npairs)
 
-    Copies specified channels from input arrays to the specified channels of output arrays
+    Copies specified channels from input arrays to the specified channels of output arrays.
 
-    :param srcv: The input array or vector of matrices.
-        All the matrices must have the same size and the same depth
+    :param srcv: Input array or vector of matrices. All the matrices must have the same size and the same depth.
 
-    :param nsrc: The number of elements in  ``srcv``
+    :param nsrc: Number of elements in  ``srcv`` .
     
-    :param dstv: The output array or vector of matrices. All the matrices  *must be allocated* , their size and depth must be the same as in  ``srcv[0]``
+    :param dstv: Output array or vector of matrices. All the matrices  *must be allocated* . Their size and depth must be the same as in  ``srcv[0]`` .
         
-    :param ndst: The number of elements in  ``dstv``
+    :param ndst: Number of elements in  ``dstv`` .
     
-    :param fromTo: The array of index pairs, specifying which channels are copied and where. ``fromTo[k*2]``  is the 0-based index of the input channel in  ``srcv``  and ``fromTo[k*2+1]``  is the index of the output channel in  ``dstv`` . Here the continuous channel numbering is used, that is, the first input image channels are indexed from  ``0``  to  ``srcv[0].channels()-1`` , the second input image channels are indexed from  ``srcv[0].channels()``  to ``srcv[0].channels() + srcv[1].channels()-1``  etc., and the same scheme is used for the output image channels. As a special case, when  ``fromTo[k*2]``  is negative, the corresponding output channel is filled with zero. ``npairs``
+    :param fromTo: Array of index pairs specifying which channels are copied and where. ``fromTo[k*2]``  is a 0-based index of the input channel in  ``srcv`` . ``fromTo[k*2+1]``  is an index of the output channel in  ``dstv`` . The continuous channel numbering is used: the first input image channels are indexed from  ``0``  to  ``srcv[0].channels()-1`` , the second input image channels are indexed from  ``srcv[0].channels()``  to ``srcv[0].channels() + srcv[1].channels()-1``  and so on. The same scheme is used for the output image channels. As a special case, when  ``fromTo[k*2]``  is negative, the corresponding output channel is filled with zero ``npairs`` .
     
 The functions ``mixChannels`` provide an advanced mechanism for shuffling image channels.
     
@@ -1713,14 +1826,13 @@ The functions ``mixChannels`` provide an advanced mechanism for shuffling image
 :func:`merge` and some forms of
 :func:`cvtColor` are partial cases of ``mixChannels`` .
 
-As an example, this code splits a 4-channel RGBA image into a 3-channel
-BGR (i.e. with R and B channels swapped) and separate alpha channel image: ::
+As an example, this code splits a 4-channel RGBA image into a 3-channel BGR (with R and B channels swapped) and separate alpha channel image: ::
 
     Mat rgba( 100, 100, CV_8UC4, Scalar(1,2,3,4) );
     Mat bgr( rgba.rows, rgba.cols, CV_8UC3 );
     Mat alpha( rgba.rows, rgba.cols, CV_8UC1 );
 
-    // forming array of matrices is quite efficient operations,
+    // forming an array of matrices is a quite efficient operation,
     // because the matrix data is not copied, only the headers
     Mat out[] = { bgr, alpha };
     // rgba[0] -> bgr[2], rgba[1] -> bgr[1],
@@ -1729,11 +1841,15 @@ BGR (i.e. with R and B channels swapped) and separate alpha channel image: ::
     mixChannels( &rgba, 1, out, 2, from_to, 4 );
 
 
-Note that, unlike many other new-style C++ functions in OpenCV (see the introduction section and
-:func:`Mat::create` ), ``mixChannels`` requires the destination arrays be pre-allocated before calling the function.
+**Note**
+
+Unlike many other new-style C++ functions in OpenCV (see the introduction section and
+:func:`Mat::create` ), ``mixChannels`` requires the destination arrays to be pre-allocated before calling the function.
 
 See Also:
-:func:`split`,:func:`merge`,:func:`cvtColor`
+:func:`split`,
+:func:`merge`,
+:func:`cvtColor`
 
 .. index:: mulSpectrums
 
@@ -1744,23 +1860,23 @@ mulSpectrums
 
 .. c:function:: void mulSpectrums(const Mat& src1, const Mat& src2, Mat& dst, int flags, bool conj=false)
 
-    Performs per-element multiplication of two Fourier spectrums.
+    Performs the per-element multiplication of two Fourier spectrums.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array; must have the same size and the same type as  ``src1``
+    :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param dst: The destination array; will have the same size and the same type as  ``src1``
+    :param dst: Destination array. It has the same size and type as  ``src1`` .
     
-    :param flags: The same flags as passed to  :func:`dft` ; only the flag  ``DFT_ROWS``  is checked for
+    :param flags: The same flags as passed to  :func:`dft` . Only the flag  ``DFT_ROWS``  is checked for.??
 
-    :param conj: The optional flag that conjugate the second source array before the multiplication (true) or not (false)
+    :param conj: Optional flag that conjugates the second source array before the multiplication (true) or not (false).
 
-The function ``mulSpectrums`` performs per-element multiplication of the two CCS-packed or complex matrices that are results of a real or complex Fourier transform.
+The function ``mulSpectrums`` performs the per-element multiplication of the two CCS-packed or complex matrices that are results of a real or complex Fourier transform.
 
 The function, together with
 :func:`dft` and
-:func:`idft` , may be used to calculate convolution (pass ``conj=false`` ) or correlation (pass ``conj=false`` ) of two arrays rapidly. When the arrays are complex, they are simply multiplied (per-element) with optional conjugation of the second array elements. When the arrays are real, they assumed to be CCS-packed (see
+:func:`idft` , may be used to calculate convolution (pass ``conj=false`` ) or correlation (pass ``conj=false`` ) of two arrays rapidly. When the arrays are complex, they are simply multiplied (per element) with an optional conjugation of the second-array elements. When the arrays are real, they are assumed to be CCS-packed (see
 :func:`dft` for details).
 
 .. index:: multiply
@@ -1774,15 +1890,15 @@ multiply
 
 .. c:function:: void multiply(const MatND& src1, const MatND& src2, MatND& dst, double scale=1)
 
-    Calculates the per-element scaled product of two arrays
+    Calculates the per-element scaled product of two arrays.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array of the same size and the same type as  ``src1``
+    :param src2: The second source array of the same size and the same type as  ``src1`` .
     
-    :param dst: The destination array; will have the same size and the same type as  ``src1``
+    :param dst: Destination array of the same size and type as  ``src1`` .
     
-    :param scale: The optional scale factor
+    :param scale: Optional scale factor.
 
 The function ``multiply`` calculates the per-element product of two arrays:
 
@@ -1791,14 +1907,23 @@ The function ``multiply`` calculates the per-element product of two arrays:
     \texttt{dst} (I)= \texttt{saturate} ( \texttt{scale} \cdot \texttt{src1} (I)  \cdot \texttt{src2} (I))
 
 There is also
-:ref:`MatrixExpressions` -friendly variant of the first function, see
+:ref:`MatrixExpressions` -friendly variant of the first function. See
 :func:`Mat::mul` .
 
-If you are looking for a matrix product, not per-element product, see
+For a not-per-element matrix product, see
 :func:`gemm` .
 
 See Also:
-:func:`add`,:func:`substract`,:func:`divide`,:ref:`MatrixExpressions`,:func:`scaleAdd`,:func:`addWeighted`,:func:`accumulate`,:func:`accumulateProduct`,:func:`accumulateSquare`,:func:`Mat::convertTo`
+:func:`add`,
+:func:`substract`,
+:func:`divide`,
+:ref:`MatrixExpressions`,
+:func:`scaleAdd`,
+:func:`addWeighted`,
+:func:`accumulate`,
+:func:`accumulateProduct`,
+:func:`accumulateSquare`,
+:func:`Mat::convertTo`
 
 .. index:: mulTransposed
 
@@ -1811,17 +1936,17 @@ mulTransposed
 
     Calculates the product of a matrix and its transposition.
 
-    :param src: The source matrix
+    :param src: Source matrix.
 
-    :param dst: The destination square matrix
+    :param dst: Destination square matrix.
 
-    :param aTa: Specifies the multiplication ordering; see the description below
+    :param aTa: Flag specifying the multiplication ordering. See the description below.
 
-    :param delta: The optional delta matrix, subtracted from  ``src``  before the multiplication. When the matrix is empty ( ``delta=Mat()`` ), it's assumed to be zero, i.e. nothing is subtracted, otherwise if it has the same size as  ``src`` , then it's simply subtracted, otherwise it is "repeated" (see  :func:`repeat` ) to cover the full  ``src``  and then subtracted. Type of the delta matrix, when it's not empty, must be the same as the type of created destination matrix, see the  ``rtype``  description
+    :param delta: Optional delta matrix subtracted from  ``src``  before the multiplication. When the matrix is empty ( ``delta=Mat()`` ), it is assumed to be zero, that is, nothing is subtracted. If it has the same size as  ``src`` , it is simply subtracted. Otherwise, it is "repeated" (see  :func:`repeat` ) to cover the full  ``src``  and then subtracted. Type of the delta matrix, when it is not empty, must be the same as the type of created destination matrix. See the  ``rtype``  description.
 
-    :param scale: The optional scale factor for the matrix product
+    :param scale: Optional scale factor for the matrix product.
 
-    :param rtype: When it's negative, the destination matrix will have the same type as  ``src`` . Otherwise, it will have  ``type=CV_MAT_DEPTH(rtype)`` , which should be either  ``CV_32F``  or  ``CV_64F``
+    :param rtype: When it is negative, the destination matrix has the same type as  ``src`` . Otherwise, it has  ``type=CV_MAT_DEPTH(rtype)`` that should be either  ``CV_32F``  or  ``CV_64F`` .
     
 The function ``mulTransposed`` calculates the product of ``src`` and its transposition:
 
@@ -1835,12 +1960,15 @@ if ``aTa=true`` , and
 
     \texttt{dst} = \texttt{scale} ( \texttt{src} - \texttt{delta} ) ( \texttt{src} - \texttt{delta} )^T
 
-otherwise. The function is used to compute covariance matrix and with zero delta can be used as a faster substitute for general matrix product
+otherwise. The function is used to compute the covariance matrix. With zero delta, it can be used as a faster substitute for general matrix product
 :math:`A*B` when
 :math:`B=A^T` .
 
 See Also:
-:func:`calcCovarMatrix`,:func:`gemm`,:func:`repeat`,:func:`reduce`
+:func:`calcCovarMatrix`,
+:func:`gemm`,
+:func:`repeat`,
+:func:`reduce`
 
 .. index:: norm
 
@@ -1863,17 +1991,17 @@ norm
 
 .. c:function:: double norm( const SparseMat& src, int normType )
 
-    Calculates absolute array norm, absolute difference norm, or relative difference norm.
+    Calculates an absolute array norm, an absolute difference norm, or a relative difference norm.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array of the same size and the same type as  ``src1``
+    :param src2: The second source array of the same size and the same type as  ``src1`` .
     
-    :param normType: Type of the norm; see the discussion below
+    :param normType: Type of the norm. See the details below.
 
-    :param mask: The optional operation mask
+    :param mask: Optional operation mask.
 
-The functions ``norm`` calculate the absolute norm of ``src1`` (when there is no ``src2`` ):
+The functions ``norm`` calculate an absolute norm of ``src1`` (when there is no ``src2`` ):
 
 .. math::
 
@@ -1899,9 +2027,9 @@ or
 
 The functions ``norm`` return the calculated norm.
 
-When there is ``mask`` parameter, and it is not empty (then it should have type ``CV_8U`` and the same size as ``src1`` ), the norm is computed only over the specified by the mask region.
+When the ``mask`` parameter is used and it is not empty (has the type ``CV_8U`` and the same size as ``src1`` ), the norm is computed only over the region specified by the mask.
 
-A multiple-channel source arrays are treated as a single-channel, that is, the results for all channels are combined.
+A multi-channel source arrays are treated as a single-channel, that is, the results for all channels are combined.
 
 .. index:: normalize
 
@@ -1916,21 +2044,21 @@ normalize
 
 .. c:function:: void normalize( const SparseMat& src, SparseMat& dst, double alpha, int normType )
 
-    Normalizes array's norm or the range
+    Normalizes an array norm or a range.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param dst: The destination array; will have the same size as  ``src``
+    :param dst: Destination array of the same size as  ``src`` .
     
-    :param alpha: The norm value to normalize to or the lower range boundary in the case of range normalization
+    :param alpha: Norm value to normalize to or the lower range boundary in case of the range normalization.
 
-    :param beta: The upper range boundary in the case of range normalization; not used for norm normalization
+    :param beta: Upper range boundary in case ofthe range normalization. It is not used for the norm normalization.
 
-    :param normType: The normalization type, see the discussion
+    :param normType: Normalization type. See the discussion.
 
-    :param rtype: When the parameter is negative, the destination array will have the same type as  ``src`` , otherwise it will have the same number of channels as  ``src``  and the depth ``=CV_MAT_DEPTH(rtype)``
+    :param rtype: When the parameter is negative, the destination array has the same type as  ``src`` . Otherwise, it has the same number of channels as  ``src``  and the depth ``=CV_MAT_DEPTH(rtype)`` .
     
-    :param mask: The optional operation mask
+    :param mask: Optional operation mask.
 
 The functions ``normalize`` scale and shift the source array elements, so that
 
@@ -1947,15 +2075,18 @@ The functions ``normalize`` scale and shift the source array elements, so that
 
 when ``normType=NORM_MINMAX`` (for dense arrays only).
 
-The optional mask specifies the sub-array to be normalize, that is, the norm or min-n-max are computed over the sub-array and then this sub-array is modified to be normalized. If you want to only use the mask to compute the norm or min-max, but modify the whole array, you can use
+The optional mask specifies a sub-array to be normalized. This means that the norm or min-n-max are computed over the sub-array, and then this sub-array is modified to be normalized. If you want to only use the mask to compute the norm or min-max but modify the whole array, you can use
 :func:`norm` and
 :func:`Mat::convertScale` /
 :func:`MatND::convertScale` /cross{SparseMat::convertScale} separately.
 
-in the case of sparse matrices, only the non-zero values are analyzed and transformed. Because of this, the range transformation for sparse matrices is not allowed, since it can shift the zero level.
+In case of sparse matrices, only the non-zero values are analyzed and transformed. Because of this, the range transformation for sparse matrices is not allowed, since it can shift the zero level.
 
 See Also:
-:func:`norm`,:func:`Mat::convertScale`,:func:`MatND::convertScale`,:func:`SparseMat::convertScale`
+:func:`norm`,
+:func:`Mat::convertScale`,
+:func:`MatND::convertScale`,
+:func:`SparseMat::convertScale`
 
 .. index:: PCA
 
@@ -1993,21 +2124,21 @@ Class for Principal Component Analysis ::
     };
 
 
-The class ``PCA`` is used to compute the special basis for a set of vectors. The basis will consist of eigenvectors of the covariance matrix computed from the input set of vectors. And also the class ``PCA`` can transform vectors to/from the new coordinate space, defined by the basis. Usually, in this new coordinate system each vector from the original set (and any linear combination of such vectors) can be quite accurately approximated by taking just the first few its components, corresponding to the eigenvectors of the largest eigenvalues of the covariance matrix. Geometrically it means that we compute projection of the vector to a subspace formed by a few eigenvectors corresponding to the dominant eigenvalues of the covariation matrix. And usually such a projection is very close to the original vector. That is, we can represent the original vector from a high-dimensional space with a much shorter vector consisting of the projected vector's coordinates in the subspace. Such a transformation is also known as Karhunen-Loeve Transform, or KLT. See
+The class ``PCA`` is used to compute a special basis for a set of vectors. The basis will consist of eigenvectors of the covariance matrix computed from the input set of vectors. The class ``PCA`` can also transform vectors to/from the new coordinate space defined by the basis. Usually, in this new coordinate system, each vector from the original set (and any linear combination of such vectors) can be quite accurately approximated by taking just the first few its components, corresponding to the eigenvectors of the largest eigenvalues of the covariance matrix. Geometrically it means that you compute a projection of the vector to a subspace formed by a few eigenvectors corresponding to the dominant eigenvalues of the covariation matrix. And usually such a projection is very close to the original vector. That is, you can represent the original vector from a high-dimensional space with a much shorter vector consisting of the projected vector's coordinates in the subspace. Such a transformation is also known as Karhunen-Loeve Transform, or KLT. See
 http://en.wikipedia.org/wiki/Principal\_component\_analysis
-The following sample is the function that takes two matrices. The first one stores the set of vectors (a row per vector) that is used to compute PCA, the second one stores another "test" set of vectors (a row per vector) that are first compressed with PCA, then reconstructed back and then the reconstruction error norm is computed and printed for each vector. ::
+The sample below is the function that takes two matrices. The first one stores the set of vectors (a row per vector) that is used to compute PCA. The second one stores another "test" set of vectors (a row per vector) that are first compressed with PCA, then reconstructed back, and then the reconstruction error norm is computed and printed for each vector. ::
 
     PCA compressPCA(const Mat& pcaset, int maxComponents,
                     const Mat& testset, Mat& compressed)
     {
         PCA pca(pcaset, // pass the data
-                Mat(), // we do not have a pre-computed mean vector,
+                Mat(), // there is no pre-computed mean vector,
                        // so let the PCA engine to compute it
                 CV_PCA_DATA_AS_ROW, // indicate that the vectors
                                     // are stored as matrix rows
                                     // (use CV_PCA_DATA_AS_COL if the vectors are
                                     // the matrix columns)
-                maxComponents // specify, how many principal components to retain
+                maxComponents // specify how many principal components to retain
                 );
         // if there is no test data, just return the computed basis, ready-to-use
         if( !testset.data )
@@ -2033,7 +2164,11 @@ The following sample is the function that takes two matrices. The first one stor
 
 
 See Also:
-:func:`calcCovarMatrix`,:func:`mulTransposed`,:func:`SVD`,:func:`dft`,:func:`dct`
+:func:`calcCovarMatrix`,
+:func:`mulTransposed`,
+:func:`SVD`,
+:func:`dft`,
+:func:`dct`
 
 .. index:: PCA::PCA
 
@@ -2047,17 +2182,17 @@ PCA::PCA
 
     PCA constructors
 
-    :param data: the input samples, stored as the matrix rows or as the matrix columns
+    :param data: Input samples stored as matrix rows or matrix columns.
 
-    :param mean: the optional mean value. If the matrix is empty ( ``Mat()`` ), the mean is computed from the data.
+    :param mean: Optional mean value. If the matrix is empty ( ``Mat()`` ), the mean is computed from the data.
 
-    :param flags: operation flags. Currently the parameter is only used to specify the data layout.
+    :param flags: Operation flags. Currently the parameter is only used to specify the data layout.
 
-        * **CV_PCA_DATA_AS_ROWS** Indicates that the input samples are stored as matrix rows.
+        * **CV_PCA_DATA_AS_ROWS** Indicate that the input samples are stored as matrix rows.
 
-        * **CV_PCA_DATA_AS_COLS** Indicates that the input samples are stored as matrix columns.
+        * **CV_PCA_DATA_AS_COLS** Indicate that the input samples are stored as matrix columns.
 
-    :param maxComponents: The maximum number of components that PCA should retain. By default, all the components are retained.
+    :param maxComponents: Maximum number of components that PCA should retain. By default, all the components are retained.
 
 The default constructor initializes empty PCA structure. The second constructor initializes the structure and calls
 :func:`PCA::operator ()` .
@@ -2073,19 +2208,19 @@ PCA::operator ()
 
     Performs Principal Component Analysis of the supplied dataset.
 
-    :param data: the input samples, stored as the matrix rows or as the matrix columns
+    :param data: Input samples stored as the matrix rows or as the matrix columns.
 
-    :param mean: the optional mean value. If the matrix is empty ( ``Mat()`` ), the mean is computed from the data.
+    :param mean: Optional mean value. If the matrix is empty ( ``Mat()`` ), the mean is computed from the data.
 
-    :param flags: operation flags. Currently the parameter is only used to specify the data layout.
+    :param flags: Operation flags. Currently the parameter is only used to specify the data layout.
 
-        * **CV_PCA_DATA_AS_ROWS** Indicates that the input samples are stored as matrix rows.
+        * **CV_PCA_DATA_AS_ROWS** Indicate that the input samples are stored as matrix rows.
 
-        * **CV_PCA_DATA_AS_COLS** Indicates that the input samples are stored as matrix columns.
+        * **CV_PCA_DATA_AS_COLS** Indicate that the input samples are stored as matrix columns.
 
-    :param maxComponents: The maximum number of components that PCA should retain. By default, all the components are retained.
+    :param maxComponents: Maximum number of components that PCA should retain. By default, all the components are retained.
 
-The operator performs PCA of the supplied dataset. It is safe to reuse the same PCA structure for multiple dataset. That is, if the  structure has been previously used with another dataset, the existing internal data is reclaimed and the new ``eigenvalues``,``eigenvectors`` and ``mean`` are allocated and computed.
+The operator performs PCA of the supplied dataset. It is safe to reuse the same PCA structure for multiple datasets. That is, if the  structure has been previously used with another dataset, the existing internal data is reclaimed and the new ``eigenvalues``,``eigenvectors`` , and ``mean`` are allocated and computed.
 
 The computed eigenvalues are sorted from the largest to the smallest and the corresponding eigenvectors are stored as ``PCA::eigenvectors`` rows.
 
@@ -2100,13 +2235,13 @@ PCA::project
 
 .. c:function:: void PCA::project(const Mat& vec, Mat& result) const
 
-    Project vector(s) to the principal component subspace
+    Projects vector(s) to the principal component subspace.
 
-    :param vec: the input vector(s). They have to have the same dimensionality and the same layout as the input data used at PCA phase. That is, if  ``CV_PCA_DATA_AS_ROWS``  had been specified, then  ``vec.cols==data.cols``  (that's vectors' dimensionality) and  ``vec.rows``  is the number of vectors to project; and similarly for the  ``CV_PCA_DATA_AS_COLS``  case.
+    :param vec: Input vector(s). They must have the same dimensionality and the same layout as the input data used at PCA phase. That is, if  ``CV_PCA_DATA_AS_ROWS``  are specified, then  ``vec.cols==data.cols``  (vectors' dimensionality) and  ``vec.rows``  is the number of vectors to project. The same is true for the  ``CV_PCA_DATA_AS_COLS``  case.
 
-    :param result: the output vectors. Let's now consider  ``CV_PCA_DATA_AS_COLS``  case. In this case the output matrix will have as many columns as the number of input vectors, i.e.  ``result.cols==vec.cols``  and the number of rows will match the number of principal components (e.g.  ``maxComponents``  parameter passed to the constructor).
+    :param result: Output vectors. In case of  ``CV_PCA_DATA_AS_COLS``  , the output matrix has as many columns as the number of input vectors. This means that  ``result.cols==vec.cols``  and the number of rows match the number of principal components (for example,  ``maxComponents``  parameter passed to the constructor).
 
-The methods project one or more vectors to the principal component subspace, where each vector projection is represented by coefficients in the principal component basis. The first form of the method returns the matrix that the second form writes to the result. So the first form can be used as a part of expression, while the second form can be more efficient in a processing loop.
+The methods project one or more vectors to the principal component subspace, where each vector projection is represented by coefficients in the principal component basis. The first form of the method returns the matrix that the second form writes to the result. So the first form can be used as a part of expression while the second form can be more efficient in a processing loop.
 
 .. index:: PCA::backProject
 
@@ -2119,14 +2254,14 @@ PCA::backProject
 
 .. c:function:: void PCA::backProject(const Mat& vec, Mat& result) const
 
-    Reconstruct vectors from their PC projections.
+    Reconstructs vectors from their PC projections.
 
     :param vec: Coordinates of the vectors in the principal component subspace. The layout and size are the same as of  ``PCA::project``  output vectors.
 
-    :param result: The reconstructed vectors. The layout and size are the same as of  ``PCA::project``  input vectors.
+    :param result: Reconstructed vectors. The layout and size are the same as of  ``PCA::project``  input vectors.
 
 The methods are inverse operations to
-:func:`PCA::project` . They take PC coordinates of projected vectors and reconstruct the original vectors. Of course, unless all the principal components have been retained, the reconstructed vectors will be different from the originals, but typically the difference will be small is if the number of components is large enough (but still much smaller than the original vector dimensionality) - that's why PCA is used after all.
+:func:`PCA::project` . They take PC coordinates of projected vectors and reconstruct the original vectors. Of course, unless all the principal components have been retained, the reconstructed vectors are different from the originals. But typically, the difference is small if the number of components is large enough (but still much smaller than the original vector dimensionality). As a result, PCA is used.
 
 .. index:: perspectiveTransform
 
@@ -2136,17 +2271,15 @@ perspectiveTransform
 --------------------
 .. c:function:: void perspectiveTransform(const Mat& src, Mat& dst, const Mat& mtx )
 
-    Performs perspective matrix transformation of vectors.
+    Performs the perspective matrix transformation of vectors.
 
-    :param src: The source two-channel or three-channel floating-point array;
-                    each element is 2D/3D vector to be transformed
+    :param src: Source two-channel or three-channel floating-point array. Each element is a 2D/3D vector to be transformed.
 
-    :param dst: The destination array; it will have the same size and same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
     
-    :param mtx: :math:`3\times 3`  or  :math:`4 \times 4`  transformation matrix
+    :param mtx: :math:`3\times 3`  or  :math:`4 \times 4`  transformation matrix.
 
-The function ``perspectiveTransform`` transforms every element of ``src``,by treating it as 2D or 3D vector, in the following way (here 3D vector transformation is shown; in the case of 2D vector transformation the
-:math:`z` component is omitted):
+The function ``perspectiveTransform`` transforms every element of ``src``by treating it as a 2D or 3D vector, in the following way:
 
 .. math::
 
@@ -2164,13 +2297,21 @@ and
 
     w =  \fork{w'}{if $w' \ne 0$}{\infty}{otherwise}
 
-Note that the function transforms a sparse set of 2D or 3D vectors. If you want to transform an image using perspective transformation, use
-:func:`warpPerspective` . If you have an inverse task, i.e. want to compute the most probable perspective transformation out of several pairs of corresponding points, you can use
+Here a 3D vector transformation is shown. In case of a 2D vector transformation, the
+:math:`z` component is omitted.
+
+**Note**:
+
+The function transforms a sparse set of 2D or 3D vectors. If you want to transform an image using perspective transformation, use
+:func:`warpPerspective` . If you have an inverse task, that is, you want to compute the most probable perspective transformation out of several pairs of corresponding points, you can use
 :func:`getPerspectiveTransform` or
 :func:`findHomography` .
 
 See Also:
-:func:`transform`,:func:`warpPerspective`,:func:`getPerspectiveTransform`,:func:`findHomography`
+:func:`transform`,
+:func:`warpPerspective`,
+:func:`getPerspectiveTransform`,
+:func:`findHomography`
 
 .. index:: phase
 
@@ -2181,15 +2322,15 @@ phase
 
 .. c:function:: void phase(const Mat& x, const Mat& y, Mat& angle, bool angleInDegrees=false)
 
-    Calculates the rotation angle of 2d vectors
+    Calculates the rotation angle of 2d vectors.
 
-    :param x: The source floating-point array of x-coordinates of 2D vectors
+    :param x: Source floating-point array of x-coordinates of 2D vectors.
 
-    :param y: The source array of y-coordinates of 2D vectors; must have the same size and the same type as  ``x``     
+    :param y: Source array of y-coordinates of 2D vectors. It must have the same size and the same type as  ``x``  .   
     
-    :param angle: The destination array of vector angles; it will have the same size and same type as  ``x``
+    :param angle: Destination array of vector angles. It has the same size and same type as  ``x`` .
     
-    :param angleInDegrees: When it is true, the function will compute angle in degrees, otherwise they will be measured in radians
+    :param angleInDegrees: When it is true, the function computes the angle in degrees. Otherwise, they are measured in radians.
 
 The function ``phase`` computes the rotation angle of each 2D vector that is formed from the corresponding elements of ``x`` and ``y`` :
 
@@ -2214,17 +2355,17 @@ polarToCart
 
     Computes x and y coordinates of 2D vectors from their magnitude and angle.
 
-    :param magnitude: The source floating-point array of magnitudes of 2D vectors. It can be an empty matrix ( ``=Mat()`` ) - in this case the function assumes that all the magnitudes are =1. If it's not empty, it must have the same size and same type as  ``angle``
+    :param magnitude: Source floating-point array of magnitudes of 2D vectors. It can be an empty matrix ( ``=Mat()`` ). In this case, the function assumes that all the magnitudes are =1. If it is not empty, it must have the same size and type as  ``angle`` .
     
-    :param angle: The source floating-point array of angles of the 2D vectors
+    :param angle: Source floating-point array of angles of the 2D vectors.
 
-    :param x: The destination array of x-coordinates of 2D vectors; will have the same size and the same type as  ``angle``     
+    :param x: Destination array of x-coordinates of 2D vectors. It has the same size and type as  ``angle`` .
     
-    :param y: The destination array of y-coordinates of 2D vectors; will have the same size and the same type as  ``angle``     
+    :param y: Destination array of y-coordinates of 2D vectors. It has the same size and type as  ``angle`` .
     
-    :param angleInDegrees: When it is true, the input angles are measured in degrees, otherwise they are measured in radians
+    :param angleInDegrees: When it is true, the input angles are measured in degrees. Otherwise. they are measured in radians.
 
-The function ``polarToCart`` computes the cartesian coordinates of each 2D vector represented by the corresponding elements of ``magnitude`` and ``angle`` :
+The function ``polarToCart`` computes the Cartesian coordinates of each 2D vector represented by the corresponding elements of ``magnitude`` and ``angle`` :
 
 .. math::
 
@@ -2234,7 +2375,13 @@ The relative accuracy of the estimated coordinates is
 :math:`\sim\,10^{-6}` .
 
 See Also:
-:func:`cartToPolar`,:func:`magnitude`,:func:`phase`,:func:`exp`,:func:`log`,:func:`pow`,:func:`sqrt`
+:func:`cartToPolar`,
+:func:`magnitude`,
+:func:`phase`,
+:func:`exp`,
+:func:`log`,
+:func:`pow`,
+:func:`sqrt`
 
 .. index:: pow
 
@@ -2249,11 +2396,11 @@ pow
 
     Raises every array element to a power.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param p: The exponent of power
+    :param p: Exponent of power.
 
-    :param dst: The destination array; will have the same size and the same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
 
 The function ``pow`` raises every element of the input array to ``p`` :
 
@@ -2261,17 +2408,21 @@ The function ``pow`` raises every element of the input array to ``p`` :
 
     \texttt{dst} (I) =  \fork{\texttt{src}(I)^p}{if \texttt{p} is integer}{|\texttt{src}(I)|^p}{otherwise}
 
-That is, for a non-integer power exponent the absolute values of input array elements are used. However, it is possible to get true values for negative values using some extra operations, as the following example, computing the 5th root of array ``src`` , shows: ::
+That is, for a non-integer power exponent, the absolute values of input array elements are used. However, it is possible to get true values for negative values using some extra operations. In the example below, computing the 5th root of array ``src``  shows: ::
 
     Mat mask = src < 0;
     pow(src, 1./5, dst);
     subtract(Scalar::all(0), dst, dst, mask);
 
 
-For some values of ``p`` , such as integer values, 0.5, and -0.5, specialized faster algorithms are used.
+For some values of ``p`` , such as integer values, 0.5 and -0.5, specialized faster algorithms are used.
 
 See Also:
-:func:`sqrt`,:func:`exp`,:func:`log`,:func:`cartToPolar`,:func:`polarToCart`
+:func:`sqrt`,
+:func:`exp`,
+:func:`log`,
+:func:`cartToPolar`,
+:func:`polarToCart`
 
 .. index:: RNG
 
@@ -2280,7 +2431,7 @@ See Also:
 RNG
 ---
 
-Random number generator class. ::
+Random number generator class ::
 
     class CV_EXPORTS RNG
     {
@@ -2291,7 +2442,7 @@ Random number generator class. ::
         RNG();
         RNG(uint64 state);
 
-        // returns 32-bit unsigned random number
+        // returns 32-bit unsigned random number
         unsigned next();
 
         // return random numbers of the specified type
@@ -2311,10 +2462,10 @@ Random number generator class. ::
         float uniform(float a, float b);
         double uniform(double a, double b);
 
-        // returns Gaussian random number with zero mean.
+        // returns the Gaussian random number with zero mean.
             double gaussian(double sigma);
 
-        // fills array with random numbers sampled from the specified distribution
+        // fills thec array with random numbers sampled from the specified distribution
         void fill( Mat& mat, int distType, const Scalar& a, const Scalar& b );
         void fill( MatND& mat, int distType, const Scalar& a, const Scalar& b );
 
@@ -2323,9 +2474,9 @@ Random number generator class. ::
     };
 
 
-The class ``RNG`` implements random number generator. It encapsulates the RNG state (currently, a 64-bit integer) and  has methods to return scalar random values and to fill arrays with random values. Currently it supports uniform and Gaussian (normal) distributions. The generator uses Multiply-With-Carry algorithm, introduced by G. Marsaglia (
+The class ``RNG`` implements the random number generator. It encapsulates the RNG state (currently, a 64-bit integer) and  has methods to return scalar random values and to fill arrays with random values. Currently it supports uniform and Gaussian (normal) distributions. The generator uses Multiply-With-Carry algorithm, introduced by G. Marsaglia (
 http://en.wikipedia.org/wiki/Multiply-with-carry
-). Gaussian-distribution random numbers are generated using Ziggurat algorithm (
+). Gaussian-distribution random numbers are generated using the Ziggurat algorithm (
 http://en.wikipedia.org/wiki/Ziggurat_algorithm
 ), introduced by G. Marsaglia and W. W. Tsang.
 
@@ -2341,9 +2492,9 @@ RNG::RNG
 
     RNG constructors
 
-    :param state: the 64-bit value used to initialize the RNG
+    :param state: 64-bit value used to initialize the RNG.
 
-These are the RNG constructors. The first form sets the state to some pre-defined value, equal to ``2**32-1`` in the current implementation. The second form sets the state to the specified value. If the user passed ``state=0`` , the constructor uses the above default value instead, to avoid the singular random number sequence, consisting of all zeros.
+These are the RNG constructors. The first form sets the state to some pre-defined value, equal to ``2**32-1`` in the current implementation. The second form sets the state to the specified value. If you passed ``state=0`` , the constructor uses the above default value instead to avoid the singular random number sequence, consisting of all zeros.
 
 .. index:: RNG::next
 
@@ -2353,9 +2504,9 @@ RNG::next
 -------------
 .. c:function:: unsigned RNG::next()
 
-    Returns the next random number
+    Returns the next random number.
 
-The method updates the state using MWC algorithm and returns the next 32-bit random number.
+The method updates the state using the MWC algorithm and returns the next 32-bit random number.
 
 .. index:: RNG::operator T
 
@@ -2378,9 +2529,9 @@ RNG::operator T
 
 .. cpp:function:: RNG::operator double()
 
-    Returns the next random number of the specified type
+    Returns the next random number of the specified type.
 
-Each of the methods updates the state using MWC algorithm and returns the next random number of the specified type. In the case of integer types the returned number is from the whole available value range for the specified type. In the case of floating-point types the returned value is from ``[0,1)`` range.
+Each of the methods updates the state using the MWC algorithm and returns the next random number of the specified type. In case of integer types, the returned number is from the available value range for the specified type. In case of floating-point types, the returned value is from ``[0,1)`` range.
 
 .. index:: RNG::operator ()
 
@@ -2392,12 +2543,12 @@ RNG::operator ()
 
 .. c:function:: unsigned RNG::operator ()(unsigned N)
 
-    Returns the next random number
+    Returns the next random number.
 
-    :param N: The upper non-inclusive boundary of the returned random number
+    :param N: Upper non-inclusive boundary of the returned random number.
 
-The methods transforms the state using MWC algorithm and returns the next random number. The first form is equivalent to
-:func:`RNG::next` , the second form returns the random number modulo ``N`` , i.e. the result will be in the range ``[0, N)`` .
+The methods transform the state using the MWC algorithm and return the next random number. The first form is equivalent to
+:func:`RNG::next` . The second form returns the random number modulo ``N`` , which means that the result is in the range ``[0, N)`` .
 
 .. index:: RNG::uniform
 
@@ -2411,34 +2562,34 @@ RNG::uniform
 
 .. c:function:: double RNG::uniform(double a, double b)
 
-    Returns the next random number sampled from the uniform distribution
+    Returns the next random number sampled from the uniform distribution.
 
-    :param a: The lower inclusive boundary of the returned random numbers
+    :param a: Lower inclusive boundary of the returned random numbers.
 
-    :param b: The upper non-inclusive boundary of the returned random numbers
+    :param b: Upper non-inclusive boundary of the returned random numbers.
 
-The methods transforms the state using MWC algorithm and returns the next uniformly-distributed random number of the specified type, deduced from the input parameter type, from the range ``[a, b)`` . There is one nuance, illustrated by the following sample: ::
+The methods transform the state using the MWC algorithm and return the next uniformly-distributed random number of the specified type, deduced from the input parameter type, from the range ``[a, b)`` . There is a nuance illustrated by the following sample: ::
 
     RNG rng;
 
-    // will always produce 0
+    // always produces 0
     double a = rng.uniform(0, 1);
 
-    // will produce double from [0, 1)
+    // produces double from [0, 1)
     double a1 = rng.uniform((double)0, (double)1);
 
-    // will produce float from [0, 1)
+    // produces float from [0, 1)
     double b = rng.uniform(0.f, 1.f);
 
-    // will produce double from [0, 1)
+    // produces double from [0, 1)
     double c = rng.uniform(0., 1.);
 
-    // will likely cause compiler error because of ambiguity:
+    // may cause compiler error because of ambiguity:
     //  RNG::uniform(0, (int)0.999999)? or RNG::uniform((double)0, 0.99999)?
     double d = rng.uniform(0, 0.999999);
 
 
-That is, the compiler does not take into account type of the variable that you assign the result of ``RNG::uniform`` to, the only thing that matters to it is the type of ``a`` and ``b`` parameters. So if you want a floating-point random number, but the range boundaries are integer numbers, either put dots in the end, if they are constants, or use explicit type cast operators, as in ``a1`` initialization above.
+The compiler does not take into account the type of the variable to which you assign the result of ``RNG::uniform`` . The only thing that matters to the compiler is the type of ``a`` and ``b`` parameters. So, if you want a floating-point random number, but the range boundaries are integer numbers, either put dots in the end, if they are constants, or use explicit type cast operators, as in the ``a1`` initialization above.
 
 .. index:: RNG::gaussian
 
@@ -2448,11 +2599,11 @@ RNG::gaussian
 -----------------
 .. c:function:: double RNG::gaussian(double sigma)
 
-    Returns the next random number sampled from the Gaussian distribution
+    Returns the next random number sampled from the Gaussian distribution.
 
-    :param sigma: The standard deviation of the distribution
+    :param sigma: Standard deviation of the distribution.
 
-The methods transforms the state using MWC algorithm and returns the next random number from the Gaussian distribution ``N(0,sigma)`` . That is, the mean value of the returned random numbers will be zero and the standard deviation will be the specified ``sigma`` .
+The method transforms the state using the MWC algorithm and returns the next random number from the Gaussian distribution ``N(0,sigma)`` . That is, the mean value of the returned random numbers is zero and the standard deviation is the specified ``sigma`` .
 
 .. index:: RNG::fill
 
@@ -2464,18 +2615,18 @@ RNG::fill
 
 .. c:function:: void RNG::fill( MatND& mat, int distType, const Scalar& a, const Scalar& b )
 
-    Fill arrays with random numbers
+    Fills arrays with random numbers.
 
     :param mat: 2D or N-dimensional matrix. Currently matrices with more than 4 channels are not supported by the methods. Use  :func:`reshape`  as a possible workaround.
 
-    :param distType: The distribution type, ``RNG::UNIFORM``  or  ``RNG::NORMAL``
+    :param distType: Distribution type, ``RNG::UNIFORM``  or  ``RNG::NORMAL`` .
     
-    :param a: The first distribution parameter. In the case of uniform distribution this is inclusive lower boundary. In the case of normal distribution this is mean value.
+    :param a: The first distribution parameter. In case of the uniform distribution, this is an inclusive lower boundary. In case of the normal distribution, this is a mean value.
 
-    :param b: The second distribution parameter. In the case of uniform distribution this is non-inclusive upper boundary. In the case of normal distribution this is standard deviation.
+    :param b: The second distribution parameter. In case of the uniform distribution, this is a non-inclusive upper boundary. In case of the normal distribution, this is a standard deviation.
 
-Each of the methods fills the matrix with the random values from the specified distribution. As the new numbers are generated, the RNG state is updated accordingly. In the case of multiple-channel images every channel is filled independently, i.e. RNG can not generate samples from multi-dimensional Gaussian distribution with non-diagonal covariation matrix directly. To do that, first, generate matrix from the distribution
-:math:`N(0, I_n)` , i.e. Gaussian distribution with zero mean and identity covariation matrix, and then transform it using
+Each of the methods fills the matrix with the random values from the specified distribution. As the new numbers are generated, the RNG state is updated accordingly. In case of multiple-channel images, every channel is filled independently, which means that RNG cannot generate samples from the multi-dimensional Gaussian distribution with non-diagonal covariation matrix directly. To do that, first, generate matrix from the distribution
+:math:`N(0, I_n)` (Gaussian distribution with zero mean and identity covariation matrix) and then transform it using
 :func:`transform` and the specific covariation matrix.
 
 .. index:: randu
@@ -2489,15 +2640,15 @@ randu
 
 .. c:function:: void randu(Mat& mtx, const Scalar& low, const Scalar& high)
 
-    Generates a single uniformly-distributed random number or array of random numbers
+    Generates a single uniformly-distributed random number or an array of random numbers.
 
-    :param mtx: The output array of random numbers. The array must be pre-allocated and have 1 to 4 channels
+    :param mtx: Output array of random numbers. The array must be pre-allocated and have 1 to 4 channels.
 
-    :param low: The inclusive lower boundary of the generated random numbers
+    :param low: Inclusive lower boundary of the generated random numbers.
 
-    :param high: The exclusive upper boundary of the generated random numbers
+    :param high: Exclusive upper boundary of the generated random numbers.
 
-The template functions ``randu`` generate and return the next uniformly-distributed random value of the specified type. ``randu<int>()`` is equivalent to ``(int)theRNG();`` etc. See
+The template functions ``randu`` generate and return the next uniformly-distributed random value of the specified type. ``randu<int>()`` is an equivalent to ``(int)theRNG();`` , and so on. See
 :func:`RNG` description.
 
 The second non-template variant of the function fills the matrix ``mtx`` with uniformly-distributed random numbers from the specified range:
@@ -2507,7 +2658,9 @@ The second non-template variant of the function fills the matrix ``mtx`` with un
     \texttt{low} _c  \leq \texttt{mtx} (I)_c <  \texttt{high} _c
 
 See Also:
-:func:`RNG`,:func:`randn`,:func:`theRNG` .
+:func:`RNG`,
+:func:`randn`,
+:func:`theRNG` 
 
 .. index:: randn
 
@@ -2518,19 +2671,20 @@ randn
 
 .. c:function:: void randn(Mat& mtx, const Scalar& mean, const Scalar& stddev)
 
-    Fills array with normally distributed random numbers
+    Fills the array with normally distributed random numbers.
 
-    :param mtx: The output array of random numbers. The array must be pre-allocated and have 1 to 4 channels
+    :param mtx: Output array of random numbers. The array must be pre-allocated and have 1 to 4 channels.
 
-    :param mean: The mean value (expectation) of the generated random numbers
+    :param mean: Mean value (expectation) of the generated random numbers.
 
-    :param stddev: The standard deviation of the generated random numbers
+    :param stddev: Standard deviation of the generated random numbers.
 
 The function ``randn`` fills the matrix ``mtx`` with normally distributed random numbers with the specified mean and standard deviation.
-is applied to the generated numbers (i.e. the values are clipped)
+?? is applied to the generated numbers (that is, the values are clipped)
 
 See Also:
-:func:`RNG`,:func:`randu`
+:func:`RNG`,
+:func:`randu`
 
 .. index:: randShuffle
 
@@ -2541,16 +2695,19 @@ randShuffle
 
 .. c:function:: void randShuffle(Mat& mtx, double iterFactor=1., RNG* rng=0)
 
-    Shuffles the array elements randomly
+    Shuffles the array elements randomly.
 
-    :param mtx: The input/output numerical 1D array
+    :param mtx: Input/output numerical 1D array.
 
-    :param iterFactor: The scale factor that determines the number of random swap operations. See the discussion
+    :param iterFactor: Scale factor that determines the number of random swap operations. See the details below.
 
-    :param rng: The optional random number generator used for shuffling. If it is zero, :func:`theRNG` () is used instead
+    :param rng: Optional random number generator used for shuffling. If it is zero, :func:`theRNG` () is used instead.
 
-The function ``randShuffle`` shuffles the specified 1D array by randomly choosing pairs of elements and swapping them. The number of such swap operations will be ``mtx.rows*mtx.cols*iterFactor`` See Also:
-:func:`RNG`,:func:`sort`
+The function ``randShuffle`` shuffles the specified 1D array by randomly choosing pairs of elements and swapping them. The number of such swap operations will be ``mtx.rows*mtx.cols*iterFactor`` .
+
+See Also:
+:func:`RNG`,
+:func:`sort`
 
 .. index:: reduce
 
@@ -2561,27 +2718,27 @@ reduce
 
 .. c:function:: void reduce(const Mat& mtx, Mat& vec, int dim, int reduceOp, int dtype=-1)
 
-    Reduces a matrix to a vector
+    Reduces a matrix to a vector.
 
-    :param mtx: The source 2D matrix
+    :param mtx: Source 2D matrix.
 
-    :param vec: The destination vector. Its size and type is defined by  ``dim``  and  ``dtype``  parameters
+    :param vec: Destination vector. Its size and type is defined by  ``dim``  and  ``dtype``  parameters.
 
-    :param dim: The dimension index along which the matrix is reduced. 0 means that the matrix is reduced to a single row and 1 means that the matrix is reduced to a single column
+    :param dim: Dimension index along which the matrix is reduced. 0 means that the matrix is reduced to a single row. 1 means that the matrix is reduced to a single column.
 
-    :param reduceOp: The reduction operation, one of:
+    :param reduceOp: Reduction operation that could be one of the following:
 
-            * **CV_REDUCE_SUM** The output is the sum of all of the matrix's rows/columns.
+            * **CV_REDUCE_SUM** The output is the sum of all rows/columns of the matrix.
 
-            * **CV_REDUCE_AVG** The output is the mean vector of all of the matrix's rows/columns.
+            * **CV_REDUCE_AVG** The output is the mean vector of all rows/columns of the matrix.
 
-            * **CV_REDUCE_MAX** The output is the maximum (column/row-wise) of all of the matrix's rows/columns.
+            * **CV_REDUCE_MAX** The output is the maximum (column/row-wise) of all rows/columns of the matrix.
 
-            * **CV_REDUCE_MIN** The output is the minimum (column/row-wise) of all of the matrix's rows/columns.
+            * **CV_REDUCE_MIN** The output is the minimum (column/row-wise) of all rows/columns of the matrix.
 
-    :param dtype: When it is negative, the destination vector will have the same type as the source matrix, otherwise, its type will be  ``CV_MAKE_TYPE(CV_MAT_DEPTH(dtype), mtx.channels())``
+    :param dtype: When it is negative, the destination vector will have the same type as the source matrix. Otherwise, its type will be  ``CV_MAKE_TYPE(CV_MAT_DEPTH(dtype), mtx.channels())`` .
     
-The function ``reduce`` reduces matrix to a vector by treating the matrix rows/columns as a set of 1D vectors and performing the specified operation on the vectors until a single row/column is obtained. For example, the function can be used to compute horizontal and vertical projections of an raster image. In the case of ``CV_REDUCE_SUM`` and ``CV_REDUCE_AVG`` the output may have a larger element bit-depth to preserve accuracy. And multi-channel arrays are also supported in these two reduction modes.
+The function ``reduce`` reduces the matrix to a vector by treating the matrix rows/columns as a set of 1D vectors and performing the specified operation on the vectors until a single row/column is obtained. For example, the function can be used to compute horizontal and vertical projections of a raster image. In case of ``CV_REDUCE_SUM`` and ``CV_REDUCE_AVG`` , the output may have a larger element bit-depth to preserve accuracy. And multi-channel arrays are also supported in these two reduction modes.
 
 See Also: :func:`repeat`
 
@@ -2596,15 +2753,15 @@ repeat
 
 .. c:function:: Mat repeat(const Mat& src, int ny, int nx)
 
-    Fill the destination array with repeated copies of the source array.
+    Fills the destination array with repeated copies of the source array.
 
-    :param src: The source array to replicate
+    :param src: Source array to replicate.
 
-    :param dst: The destination array; will have the same type as  ``src``
+    :param dst: Destination array of the same type as  ``src`` .
     
-    :param ny: How many times the  ``src``  is repeated along the vertical axis
+    :param ny: Flag to specify how many times the  ``src``  is repeated along the vertical axis.
 
-    :param nx: How many times the  ``src``  is repeated along the horizontal axis
+    :param nx: Flag to specify how many times the  ``src``  is repeated along the horizontal axis.
 
 The functions
 :func:`repeat` duplicate the source array one or more times along each of the two axes:
@@ -2614,8 +2771,11 @@ The functions
     \texttt{dst} _{ij}= \texttt{src} _{i \mod \texttt{src.rows} , \; j \mod \texttt{src.cols} }
 
 The second variant of the function is more convenient to use with
-:ref:`MatrixExpressions` See Also:
-:func:`reduce`,:ref:`MatrixExpressions`
+:ref:`MatrixExpressions` . 
+
+See Also:
+:func:`reduce`,
+:ref:`MatrixExpressions`
 
 .. index:: saturate_cast
 
@@ -2640,24 +2800,28 @@ saturate_cast
 
 .. c:function:: template<typename _Tp> inline _Tp saturate_cast(double v)
 
-    Template function for accurate conversion from one primitive type to another
+    Template function for accurate conversion from one primitive type to another.
 
-    :param v: The function parameter
+    :param v: Function parameter.
 
-The functions ``saturate_cast`` resembles the standard C++ cast operations, such as ``static_cast<T>()`` etc. They perform an efficient and accurate conversion from one primitive type to another, see the introduction. "saturate" in the name means that when the input value ``v`` is out of range of the target type, the result will not be formed just by taking low bits of the input, but instead the value will be clipped. For example: ::
+The functions ``saturate_cast`` resemble the standard C++ cast operations, such as ``static_cast<T>()`` and others. They perform an efficient and accurate conversion from one primitive type to another (see the introduction chapter). ``saturate`` in the name means that when the input value ``v`` is out of the range of the target type, the result is not formed just by taking low bits of the input, but instead the value is clipped. For example: ::
 
     uchar a = saturate_cast<uchar>(-100); // a = 0 (UCHAR_MIN)
     short b = saturate_cast<short>(33333.33333); // b = 32767 (SHRT_MAX)
 
 
-Such clipping is done when the target type is ``unsigned char, signed char, unsigned short or signed short`` - for 32-bit integers no clipping is done.
+Such clipping is done when the target type is ``unsigned char`` , ``signed char`` , ``unsigned short`` or ``signed short`` . For 32-bit integers, no clipping is done.
 
-When the parameter is floating-point value and the target type is an integer (8-, 16- or 32-bit), the floating-point value is first rounded to the nearest integer and then clipped if needed (when the target type is 8- or 16-bit).
+When the parameter is floating-point value and the target type is an integer (8-, 16- or 32-bit), the floating-point value is first rounded to the nearest integer and then clipped if needed (when the target type is 8- or 16-bit).
 
-This operation is used in most simple or complex image processing functions in OpenCV.
+This operation is used in the simplest or most complex image processing functions in OpenCV.
 
 See Also:
-:func:`add`,:func:`subtract`,:func:`multiply`,:func:`divide`,:func:`Mat::convertTo`
+:func:`add`,
+:func:`subtract`,
+:func:`multiply`,
+:func:`divide`,
+:func:`Mat::convertTo`
 
 .. index:: scaleAdd
 
@@ -2672,13 +2836,13 @@ scaleAdd
 
     Calculates the sum of a scaled array and another array.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param scale: Scale factor for the first array
+    :param scale: Scale factor for the first array.
 
-    :param src2: The second source array; must have the same size and the same type as  ``src1``
+    :param src2: The second source array of the same size and type as  ``src1`` .
     
-    :param dst: The destination array; will have the same size and the same type as  ``src1``
+    :param dst: Destination array of the same size and type as  ``src1`` .
     
 The function ``scaleAdd`` is one of the classical primitive linear algebra operations, known as ``DAXPY`` or ``SAXPY`` in `BLAS <http://en.wikipedia.org/wiki/Basic_Linear_Algebra_Subprograms>`_. It calculates the sum of a scaled array and another array:
 
@@ -2694,7 +2858,12 @@ The function can also be emulated with a matrix expression, for example: ::
 
 
 See Also:
-:func:`add`,:func:`addWeighted`,:func:`subtract`,:func:`Mat::dot`,:func:`Mat::convertTo`,:ref:`MatrixExpressions`
+:func:`add`,
+:func:`addWeighted`,
+:func:`subtract`,
+:func:`Mat::dot`,
+:func:`Mat::convertTo`,
+:ref:`MatrixExpressions`
 
 .. index:: setIdentity
 
@@ -2705,11 +2874,11 @@ setIdentity
 
 .. c:function:: void setIdentity(Mat& dst, const Scalar& value=Scalar(1))
 
-    Initializes a scaled identity matrix
+    Initializes a scaled identity matrix.
 
-    :param dst: The matrix to initialize (not necessarily square)
+    :param dst: Matrix to initialize (not necessarily square).
 
-    :param value: The value to assign to the diagonal elements
+    :param value: Value to assign to diagonal elements.
 
 The function
 :func:`setIdentity` initializes a scaled identity matrix:
@@ -2725,7 +2894,11 @@ The function can also be emulated using the matrix initializers and the matrix e
 
 
 See Also:
-:func:`Mat::zeros`,:func:`Mat::ones`,:ref:`MatrixExpressions`,:func:`Mat::setTo`,:func:`Mat::operator=`
+:func:`Mat::zeros`,
+:func:`Mat::ones`,
+:ref:`MatrixExpressions`,
+:func:`Mat::setTo`,
+:func:`Mat::operator=`
 
 .. index:: solve
 
@@ -2738,25 +2911,25 @@ solve
 
     Solves one or more linear systems or least-squares problems.
 
-    :param src1: The input matrix on the left-hand side of the system
+    :param src1: Input matrix on the left-hand side of the system.
 
-    :param src2: The input matrix on the right-hand side of the system
+    :param src2: Input matrix on the right-hand side of the system.
 
-    :param dst: The output solution
+    :param dst: Output solution.
 
-    :param flags: The solution (matrix inversion) method
+    :param flags: Solution (matrix inversion) method.
 
-            * **DECOMP_LU** Gaussian elimination with optimal pivot element chosen
+            * **DECOMP_LU** Gaussian elimination with optimal pivot element chosen.
 
-            * **DECOMP_CHOLESKY** Cholesky  :math:`LL^T`  factorization; the matrix  ``src1``  must be symmetrical and positively defined
+            * **DECOMP_CHOLESKY** Cholesky  :math:`LL^T`  factorization. The matrix  ``src1``  must be symmetrical and positively defined.
 
-            * **DECOMP_EIG** Eigenvalue decomposition; the matrix  ``src1``  must be symmetrical
+            * **DECOMP_EIG** Eigenvalue decomposition. The matrix  ``src1``  must be symmetrical.
 
-            * **DECOMP_SVD** Singular value decomposition (SVD) method; the system can be over-defined and/or the matrix  ``src1``  can be singular
+            * **DECOMP_SVD** Singular value decomposition (SVD) method. The system can be over-defined and/or the matrix  ``src1``  can be singular.
 
-            * **DECOMP_QR** QR factorization; the system can be over-defined and/or the matrix  ``src1``  can be singular
+            * **DECOMP_QR** QR factorization. The system can be over-defined and/or the matrix  ``src1``  can be singular.
 
-            * **DECOMP_NORMAL** While all the previous flags are mutually exclusive, this flag can be used together with any of the previous. It means that the normal equations  :math:`\texttt{src1}^T\cdot\texttt{src1}\cdot\texttt{dst}=\texttt{src1}^T\texttt{src2}`  are solved instead of the original system  :math:`\texttt{src1}\cdot\texttt{dst}=\texttt{src2}`
+            * **DECOMP_NORMAL** While all the previous flags are mutually exclusive, this flag can be used together with any of the previous. It means that the normal equations  :math:`\texttt{src1}^T\cdot\texttt{src1}\cdot\texttt{dst}=\texttt{src1}^T\texttt{src2}`  are solved instead of the original system  :math:`\texttt{src1}\cdot\texttt{dst}=\texttt{src2}` .
             
 The function ``solve`` solves a linear system or least-squares problem (the latter is possible with SVD or QR methods, or by specifying the flag ``DECOMP_NORMAL`` ):
 
@@ -2765,14 +2938,18 @@ The function ``solve`` solves a linear system or least-squares problem (the latt
     \texttt{dst} =  \arg \min _X \| \texttt{src1} \cdot \texttt{X} -  \texttt{src2} \|
 
 If ``DECOMP_LU`` or ``DECOMP_CHOLESKY`` method is used, the function returns 1 if ``src1`` (or
-:math:`\texttt{src1}^T\texttt{src1}` ) is non-singular and 0 otherwise; in the latter case ``dst`` is not valid. Other methods find some pseudo-solution in the case of singular left-hand side part.
+:math:`\texttt{src1}^T\texttt{src1}` ) is non-singular. Otherwise, it returns 0. In the latter case, ``dst`` is not valid. Other methods find a pseudo-solution in case of a singular left-hand side part.
+
+**Note** 
 
-Note that if you want to find unity-norm solution of an under-defined singular system
+If you want to find a unity-norm solution of an under-defined singular system
 :math:`\texttt{src1}\cdot\texttt{dst}=0` , the function ``solve`` will not do the work. Use
 :func:`SVD::solveZ` instead.
 
 See Also:
-:func:`invert`,:func:`SVD`,:func:`eigen`
+:func:`invert`,
+:func:`SVD`,
+:func:`eigen`
 
 .. index:: solveCubic
 
@@ -2784,25 +2961,25 @@ solveCubic
 
     Finds the real roots of a cubic equation.
 
-    :param coeffs: The equation coefficients, an array of 3 or 4 elements
+    :param coeffs: Equation coefficients, an array of 3 or 4 elements.
 
-    :param roots: The destination array of real roots which will have 1 or 3 elements
+    :param roots: Destination array of real roots that has 1 or 3 elements.
 
 The function ``solveCubic`` finds the real roots of a cubic equation:
 
-(if coeffs is a 4-element vector)
+(if ``coeffs`` is a 4-element vector)
 
 .. math::
 
     \texttt{coeffs} [0] x^3 +  \texttt{coeffs} [1] x^2 +  \texttt{coeffs} [2] x +  \texttt{coeffs} [3] = 0
 
-or (if coeffs is 3-element vector):
+or (if ``coeffs`` is a 3-element vector):
 
 .. math::
 
     x^3 +  \texttt{coeffs} [0] x^2 +  \texttt{coeffs} [1] x +  \texttt{coeffs} [2] = 0
 
-The roots are stored to ``roots`` array.
+The roots are stored in the ``roots`` array.
 
 .. index:: solvePoly
 
@@ -2813,13 +2990,13 @@ solvePoly
 
 .. c:function:: void solvePoly(const Mat& coeffs, Mat& roots, int maxIters=20, int fig=100)
 
-    Finds the real or complex roots of a polynomial equation
+    Finds the real or complex roots of a polynomial equation.
 
-    :param coeffs: The array of polynomial coefficients
+    :param coeffs: Array of polynomial coefficients.
 
-    :param roots: The destination (complex) array of roots
+    :param roots: Destination (complex) array of roots.
 
-    :param maxIters: The maximum number of iterations the algorithm does
+    :param maxIters: Maximum number of iterations the algorithm does.
 
     :param fig:
 
@@ -2838,26 +3015,27 @@ sort
 
 .. c:function:: void sort(const Mat& src, Mat& dst, int flags)
 
-    Sorts each row or each column of a matrix
+    Sorts each row or each column of a matrix.
 
-    :param src: The source single-channel array
+    :param src: Source single-channel array.
 
-    :param dst: The destination array of the same size and the same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
     
-    :param flags: The operation flags, a combination of the following values:
+    :param flags: Operation flags, a combination of the following values:
 
-            * **CV_SORT_EVERY_ROW** Each matrix row is sorted independently
+            * **CV_SORT_EVERY_ROW** Each matrix row is sorted independently.
 
-            * **CV_SORT_EVERY_COLUMN** Each matrix column is sorted independently. This flag and the previous one are mutually exclusive
+            * **CV_SORT_EVERY_COLUMN** Each matrix column is sorted independently. This flag and the previous one are mutually exclusive.
 
-            * **CV_SORT_ASCENDING** Each matrix row is sorted in the ascending order
+            * **CV_SORT_ASCENDING** Each matrix row is sorted in the ascending order.
 
-            * **CV_SORT_DESCENDING** Each matrix row is sorted in the descending order. This flag and the previous one are also mutually exclusive
+            * **CV_SORT_DESCENDING** Each matrix row is sorted in the descending order. This flag and the previous one are also mutually exclusive.
 
 The function ``sort`` sorts each matrix row or each matrix column in ascending or descending order. If you want to sort matrix rows or columns lexicographically, you can use STL ``std::sort`` generic function with the proper comparison predicate.
 
 See Also:
-:func:`sortIdx`,:func:`randShuffle`
+:func:`sortIdx`,
+:func:`randShuffle`
 
 .. index:: sortIdx
 
@@ -2868,23 +3046,23 @@ sortIdx
 
 .. c:function:: void sortIdx(const Mat& src, Mat& dst, int flags)
 
-    Sorts each row or each column of a matrix
+    Sorts each row or each column of a matrix.
 
-    :param src: The source single-channel array
+    :param src: Source single-channel array.
 
-    :param dst: The destination integer array of the same size as  ``src``
+    :param dst: Destination integer array of the same size as  ``src`` .
     
-    :param flags: The operation flags, a combination of the following values:
+    :param flags: Operation flags that could be a combination of the following values:
 
-            * **CV_SORT_EVERY_ROW** Each matrix row is sorted independently
+            * **CV_SORT_EVERY_ROW** Each matrix row is sorted independently.
 
-            * **CV_SORT_EVERY_COLUMN** Each matrix column is sorted independently. This flag and the previous one are mutually exclusive
+            * **CV_SORT_EVERY_COLUMN** Each matrix column is sorted independently. This flag and the previous one are mutually exclusive.
 
-            * **CV_SORT_ASCENDING** Each matrix row is sorted in the ascending order
+            * **CV_SORT_ASCENDING** Each matrix row is sorted in the ascending order.
 
-            * **CV_SORT_DESCENDING** Each matrix row is sorted in the descending order. This flag and the previous one are also mutually exclusive
+            * **CV_SORT_DESCENDING** Each matrix row is sorted in the descending order. This flag and the previous one are also mutually exclusive.
 
-The function ``sortIdx`` sorts each matrix row or each matrix column in ascending or descending order. Instead of reordering the elements themselves, it stores the indices of sorted elements in the destination array. For example: ::
+The function ``sortIdx`` sorts each matrix row or each matrix column in the ascending or descending order. Instead of reordering the elements themselves, it stores the indices of sorted elements in the destination array. For example: ::
 
     Mat A = Mat::eye(3,3,CV_32F), B;
     sortIdx(A, B, CV_SORT_EVERY_ROW + CV_SORT_ASCENDING);
@@ -2894,7 +3072,8 @@ The function ``sortIdx`` sorts each matrix row or each matrix column in ascendin
 
 
 See Also:
-:func:`sort`,:func:`randShuffle`
+:func:`sort`,
+:func:`randShuffle`
 
 .. index:: split
 
@@ -2911,21 +3090,25 @@ split
 
 .. c:function:: void split(const MatND& mtx, vector<MatND>& mv)
 
-    Divides multi-channel array into several single-channel arrays
+    Divides a multi-channel array into several single-channel arrays.
 
-    :param mtx: The source multi-channel array
+    :param mtx: Source multi-channel array.
 
-    :param mv: The destination array or vector of arrays; The number of arrays must match  ``mtx.channels()`` . The arrays themselves will be reallocated if needed
+    :param mv: Destination array or vector of arrays. The number of arrays must match  ``mtx.channels()`` . The arrays themselves are reallocated, if needed.
 
-The functions ``split`` split multi-channel array into separate single-channel arrays:
+The functions ``split`` split multi-channel array into separate single-channel arrays:
 
 .. math::
 
     \texttt{mv} [c](I) =  \texttt{mtx} (I)_c
 
-If you need to extract a single-channel or do some other sophisticated channel permutation, use
-:func:`mixChannels` See Also:
-:func:`merge`,:func:`mixChannels`,:func:`cvtColor`
+If you need to extract a single channel or do some other sophisticated channel permutation, use
+:func:`mixChannels` .
+
+See Also:
+:func:`merge`,
+:func:`mixChannels`,
+:func:`cvtColor`
 
 .. index:: sqrt
 
@@ -2938,16 +3121,17 @@ sqrt
 
 .. c:function:: void sqrt(const MatND& src, MatND& dst)
 
-    Calculates square root of array elements
+    Calculates a quare root of array elements.
 
-    :param src: The source floating-point array
+    :param src: Source floating-point array.
 
-    :param dst: The destination array; will have the same size and the same type as  ``src``
+    :param dst: Destination array of the same size and type as  ``src`` .
     
-The functions ``sqrt`` calculate square root of each source array element. in the case of multi-channel arrays each channel is processed independently. The accuracy is approximately the same as of the built-in ``std::sqrt`` .
+The functions ``sqrt`` calculate a square root of each source array element. In case of multi-channel arrays, each channel is processed independently. The accuracy is approximately the same as of the built-in ``std::sqrt`` .
 
 See Also:
-:func:`pow`,:func:`magnitude`
+:func:`pow`,
+:func:`magnitude`
 
 .. index:: subtract
 
@@ -2972,42 +3156,42 @@ subtract
 
 .. c:function:: void subtract(const Scalar& sc, const MatND& src2, MatND& dst, const MatND& mask=MatND())
 
-    Calculates per-element difference between two arrays or array and a scalar
+    Calculates the per-element difference between two arrays or array and a scalar.
 
-    :param src1: The first source array
+    :param src1: The first source array.
 
-    :param src2: The second source array. It must have the same size and same type as  ``src1``
+    :param src2: The second source array. It must have the same size and same type as  ``src1`` .
     
-    :param sc: Scalar; the first or the second input parameter
+    :param sc: Scalar that could be the first or the second input parameter.
 
-    :param dst: The destination array; it will have the same size and same type as  ``src1`` ; see  ``Mat::create``     
+    :param dst: Destination array of the same size and type as  ``src1`` . See  ``Mat::create``  .   
     
-    :param mask: The optional operation mask, 8-bit single channel array; specifies elements of the destination array to be changed
+    :param mask: Optional operation mask. This is an 8-bit single channel array that specifies elements of the destination array to be changed.
 
-The functions ``subtract`` compute
+The functions ``subtract`` compute the following values:
 
 *
-    the difference between two arrays
+    Difference between two arrays
 
     .. math::
 
         \texttt{dst} (I) =  \texttt{saturate} ( \texttt{src1} (I) -  \texttt{src2} (I)) \quad \texttt{if mask} (I) \ne0
 
 *
-    the difference between array and a scalar:
+    Difference between an array and a scalar:
 
     .. math::
 
         \texttt{dst} (I) =  \texttt{saturate} ( \texttt{src1} (I) -  \texttt{sc} ) \quad \texttt{if mask} (I) \ne0
 
 *
-    the difference between scalar and an array:
+    Difference between a scalar and an array:
 
     .. math::
 
         \texttt{dst} (I) =  \texttt{saturate} ( \texttt{sc} -  \texttt{src2} (I)) \quad \texttt{if mask} (I) \ne0
 
-where ``I`` is multi-dimensional index of array elements.
+where ``I`` is a multi-dimensional index of the array elements.
 
 The first function in the above list can be replaced with matrix expressions: ::
 
@@ -3016,7 +3200,11 @@ The first function in the above list can be replaced with matrix expressions: ::
 
 
 See Also:
-:func:`add`,:func:`addWeighted`,:func:`scaleAdd`,:func:`convertScale`,:ref:`MatrixExpressions`,.
+:func:`add`,
+:func:`addWeighted`,
+:func:`scaleAdd`,
+:func:`convertScale`,
+:ref:`MatrixExpressions`
 
 .. index:: SVD
 
@@ -3040,8 +3228,8 @@ Class for computing Singular Value Decomposition ::
         // decomposes A into u, w and vt.
         SVD& operator ()( const Mat& A, int flags=0 );
 
-        // finds such vector x, norm(x)=1, so that A*x = 0,
-        // where A is singular matrix
+        // finds vector x, norm(x)=1, so that A*x = 0,
+        // where A is singular matrix
         static void solveZ( const Mat& A, Mat& x );
         // does back-subsitution:
         // x = vt.t()*inv(w)*u.t()*rhs ~ inv(A)*rhs
@@ -3053,11 +3241,14 @@ Class for computing Singular Value Decomposition ::
     };
 
 
-The class ``SVD`` is used to compute Singular Value Decomposition of a floating-point matrix and then use it to solve least-square problems, under-determined linear systems, invert matrices, compute condition numbers etc.
-For a bit faster operation you can pass ``flags=SVD::MODIFY_A|...`` to modify the decomposed matrix when it is not necessarily to preserve it. If you want to compute condition number of a matrix or absolute value of its determinant - you do not need ``u`` and ``vt`` , so you can pass ``flags=SVD::NO_UV|...`` . Another flag ``FULL_UV`` indicates that full-size ``u`` and ``vt`` must be computed, which is not necessary most of the time.
+The class ``SVD`` is used to compute Singular Value Decomposition of a floating-point matrix. The Singular Value Decomposition is used to solve least-square problems, under-determined linear systems, invert matrices, compute condition numbers, and so on.
+For a bit faster operation, you can pass ``flags=SVD::MODIFY_A|...`` to modify the decomposed matrix when it is not necessary to preserve it. If you want to compute a condition number of a matrix or an absolute value of its determinant, you do not need ``u`` and ``vt`` . You can pass ``flags=SVD::NO_UV|...`` . Another flag ``FULL_UV`` indicates that full-size ``u`` and ``vt`` must be computed, which is not necessary most of the time.
 
 See Also:
-:func:`invert`,:func:`solve`,:func:`eigen`,:func:`determinant`
+:func:`invert`,
+:func:`solve`,
+:func:`eigen`,
+:func:`determinant`
 
 .. index:: SVD::SVD
 
@@ -3072,17 +3263,17 @@ SVD::SVD
 
     SVD constructors
 
-    :param A: The decomposed matrix
+    :param A: Decomposed matrix.
 
-    :param flags: Operation flags
+    :param flags: Operation flags.
 
-        * **SVD::MODIFY_A** The algorithm can modify the decomposed matrix. It can save some space and speed-up processing a bit
+        * **SVD::MODIFY_A** Use the algorithm to modify the decomposed matrix. It can save some space and speed-up processing a bit.
 
-        * **SVD::NO_UV** Indicates that only the vector of singular values  ``w``  is to be computed, while  ``u``  and  ``vt``  will be set to empty matrices
+        * **SVD::NO_UV** Indicate that only a vector of singular values  ``w``  is to be computed, while  ``u``  and  ``vt``  will be set to empty matrices.
 
         * **SVD::FULL_UV** When the matrix is not square, by default the algorithm produces  ``u``  and  ``vt``  matrices of sufficiently large size for the further  ``A``  reconstruction. If, however, ``FULL_UV``  flag is specified, ``u``  and  ``vt``  will be full-size square orthogonal matrices.
 
-The first constructor initializes empty ``SVD`` structure. The second constructor initializes empty ``SVD`` structure and then calls
+The first constructor initializes an empty ``SVD`` structure. The second constructor initializes an empty ``SVD`` structure and then calls
 :func:`SVD::operator ()` .
 
 .. index:: SVD::operator ()
@@ -3094,19 +3285,19 @@ SVD::operator ()
 
 .. c:function:: SVD& SVD::operator ()( const Mat& A, int flags=0 )
 
-    Performs SVD of a matrix
+    Performs SVD of a matrix.
 
-    :param A: The decomposed matrix
+    :param A: Decomposed matrix.
 
-    :param flags: Operation flags
+    :param flags: Operation flags.
 
-        * **SVD::MODIFY_A** The algorithm can modify the decomposed matrix. It can save some space and speed-up processing a bit
+        * **SVD::MODIFY_A** Use the algorithm to modify the decomposed matrix. It can save some space and speed-up processing a bit.
 
-        * **SVD::NO_UV** Only singular values are needed. The algorithm will not compute  ``u``  and  ``vt``  matrices
+        * **SVD::NO_UV** Use only singular values. The algorithm does not compute  ``u``  and  ``vt``  matrices.
 
-        * **SVD::FULL_UV** When the matrix is not square, by default the algorithm produces  ``u``  and  ``vt``  matrices of sufficiently large size for the further  ``A``  reconstruction. If, however, ``FULL_UV``  flag is specified, ``u``  and  ``vt``  will be full-size square orthogonal matrices.
+        * **SVD::FULL_UV** When the matrix is not square, by default the algorithm produces  ``u``  and  ``vt``  matrices of sufficiently large size for the further  ``A``  reconstruction. If, however, the ``FULL_UV``  flag is specified, ``u``  and  ``vt``  are full-size square orthogonal matrices.
 
-The operator performs singular value decomposition of the supplied matrix. The ``u``,``vt`` and the vector of singular values ``w`` are stored in the structure. The same ``SVD`` structure can be reused many times with different matrices. Each time, if needed, the previous ``u``,``vt`` and ``w`` are reclaimed and the new matrices are created, which is all handled by
+The operator performs the singular value decomposition of the supplied matrix. The ``u``,``vt`` , and the vector of singular values ``w`` are stored in the structure. The same ``SVD`` structure can be reused many times with different matrices. Each time, if needed, the previous ``u``,``vt`` , and ``w`` are reclaimed and the new matrices are created, which is all handled by
 :func:`Mat::create` .
 
 .. index:: SVD::solveZ
@@ -3118,16 +3309,16 @@ SVD::solveZ
 
 .. c:function:: static void SVD::solveZ( const Mat& A, Mat& x )
 
-    Solves under-determined singular linear system
+    Solves an under-determined singular linear system.
 
-    :param A: The left-hand-side matrix.
+    :param A: Left-hand-side matrix.
 
-    :param x: The found solution
+    :param x: Found solution.
 
-The method finds unit-length solution
+The method finds unit-length solution
 **x**
 of the under-determined system
-:math:`A x = 0` . Theory says that such system has infinite number of solutions, so the algorithm finds the unit-length solution as the right singular vector corresponding to the smallest singular value (which should be 0). In practice, because of round errors and limited floating-point accuracy, the input matrix can appear to be close-to-singular rather than just singular. So, strictly speaking, the algorithm solves the following problem:
+:math:`A x = 0` . Theory says that such a system has an infinite number of solutions, so the algorithm finds a unit-length solution as the right singular vector corresponding to the smallest singular value (which should be 0). In practice, because of round errors and limited floating-point accuracy, the input matrix can appear to be close-to-singular rather than just singular. So, strictly speaking, the algorithm solves the following problem:
 
 .. math::
 
@@ -3142,20 +3333,24 @@ SVD::backSubst
 
 .. c:function:: void SVD::backSubst( const Mat& rhs, Mat& x ) const
 
-    Performs singular value back substitution
+    Performs a singular value back substitution.
 
-    :param rhs: The right-hand side of a linear system  :math:`\texttt{A} \texttt{x} = \texttt{rhs}`  being solved, where  ``A``  is the matrix passed to  :func:`SVD::SVD`  or  :func:`SVD::operator ()`
+    :param rhs: Right-hand side of a linear system  :math:`\texttt{A} \texttt{x} = \texttt{rhs}`  to be solved, where  ``A``  is the matrix passed to  :func:`SVD::SVD`  or  :func:`SVD::operator ()` .
     
-    :param x: The found solution of the system
+    :param x: Found solution of the system.
 
-The method computes back substitution for the specified right-hand side:
+The method computes back substitution for the specified right-hand side:
 
 .. math::
 
     \texttt{x} =  \texttt{vt} ^T  \cdot diag( \texttt{w} )^{-1}  \cdot \texttt{u} ^T  \cdot \texttt{rhs} \sim \texttt{A} ^{-1}  \cdot \texttt{rhs}
 
-Using this technique you can either get a very accurate solution of convenient linear system, or the best (in the least-squares terms) pseudo-solution of an overdetermined linear system. Note that explicit SVD with the further back substitution only makes sense if you need to solve many linear systems with the same left-hand side (e.g. ``A`` ). If all you need is to solve a single system (possibly with multiple ``rhs`` immediately available), simply call
-:func:`solve` add pass ``DECOMP_SVD`` there - it will do absolutely the same thing.
+Using this technique you can either get a very accurate solution of the convenient linear system, or the best (in the least-squares terms) pseudo-solution of an overdetermined linear system. 
+
+**Note**
+
+Explicit SVD with the further back substitution only makes sense if you need to solve many linear systems with the same left-hand side (for example, ``A`` ). If all you need is to solve a single system (possibly with multiple ``rhs`` immediately available), simply call
+:func:`solve` add pass ``DECOMP_SVD`` there. It does absolutely the same thing.
 
 .. index:: sum
 
@@ -3168,14 +3363,19 @@ sum
 
 .. c:function:: Scalar sum(const MatND& mtx)
 
-    Calculates sum of array elements
+    Calculates the sum of array elements.
 
-    :param mtx: The source array; must have 1 to 4 channels
+    :param mtx: Source array that must have from 1 to 4 channels.
 
 The functions ``sum`` calculate and return the sum of array elements, independently for each channel.
 
 See Also:
-:func:`countNonZero`,:func:`mean`,:func:`meanStdDev`,:func:`norm`,:func:`minMaxLoc`,:func:`reduce`
+:func:`countNonZero`,
+:func:`mean`,
+:func:`meanStdDev`,
+:func:`norm`,
+:func:`minMaxLoc`,
+:func:`reduce`
 
 .. index:: theRNG
 
@@ -3184,14 +3384,16 @@ theRNG
 
 .. c:function:: RNG& theRNG()
 
-    Returns the default random number generator
+    Returns the default random number generator.
 
-The function ``theRNG`` returns the default random number generator. For each thread there is separate random number generator, so you can use the function safely in multi-thread environments. If you just need to get a single random number using this generator or initialize an array, you can use
+The function ``theRNG`` returns the default random number generator. For each thread, there is a separate random number generator, so you can use the function safely in multi-thread environments. If you just need to get a single random number using this generator or initialize an array, you can use
 :func:`randu` or
-:func:`randn` instead. But if you are going to generate many random numbers inside a loop, it will be much faster to use this function to retrieve the generator and then use ``RNG::operator _Tp()`` .
+:func:`randn` instead. But if you are going to generate many random numbers inside a loop, it is much faster to use this function to retrieve the generator and then use ``RNG::operator _Tp()`` .
 
 See Also:
-:func:`RNG`,:func:`randu`,:func:`randn`
+:func:`RNG`,
+:func:`randu`,
+:func:`randn`
 
 .. index:: trace
 
@@ -3202,9 +3404,9 @@ trace
 
 .. c:function:: Scalar trace(const Mat& mtx)
 
-    Returns the trace of a matrix
+    Returns the trace of a matrix.
 
-    :param mtx: The source matrix
+    :param mtx: Source matrix.
 
 The function ``trace`` returns the sum of the diagonal elements of the matrix ``mtx`` .
 
@@ -3221,15 +3423,15 @@ transform
 
 .. c:function:: void transform(const Mat& src, Mat& dst, const Mat& mtx )
 
-    Performs matrix transformation of every array element.
+    Performs the matrix transformation of every array element.
 
-    :param src: The source array; must have as many channels (1 to 4) as  ``mtx.cols``  or  ``mtx.cols-1``
+    :param src: Source array that must have as many channels (1 to 4) as  ``mtx.cols``  or  ``mtx.cols-1``.
     
-    :param dst: The destination array; will have the same size and depth as  ``src``  and as many channels as  ``mtx.rows``     
+    :param dst: Destination array of the same size and depth as  ``src`` . It has as many channels as  ``mtx.rows``  .   
     
-    :param mtx: The transformation matrix
+    :param mtx: Transformation matrix.
 
-The function ``transform`` performs matrix transformation of every element of array ``src`` and stores the results in ``dst`` :
+The function ``transform`` performs the matrix transformation of every element of the array ``src`` and stores the results in ``dst`` :
 
 .. math::
 
@@ -3243,20 +3445,24 @@ The function ``transform`` performs matrix transformation of every element of ar
 
 (when ``mtx.cols=src.channels()+1`` )
 
-That is, every element of an ``N`` -channel array ``src`` is
-considered as ``N`` -element vector, which is transformed using
-a
+Every element of the ``N`` -channel array ``src`` is
+considered as an ``N`` -element vector that is transformed using
+the
 :math:`\texttt{M} \times \texttt{N}` or
-:math:`\texttt{M} \times \texttt{N+1}` matrix ``mtx`` into
-an element of ``M`` -channel array ``dst`` .
+:math:`\texttt{M} \times \texttt{N+1}` matrix ``mtx`` into??
+an element of the ``M`` -channel array ``dst`` .
 
 The function may be used for geometrical transformation of
 :math:`N` -dimensional
 points, arbitrary linear color space transformation (such as various kinds of RGB
-:math:`\rightarrow` YUV transforms), shuffling the image channels and so forth.
+:math:`\rightarrow` YUV transforms), shuffling the image channels, and so forth.
 
 See Also:
-:func:`perspectiveTransform`,:func:`getAffineTransform`,:func:`estimateRigidTransform`,:func:`warpAffine`,:func:`warpPerspective`
+:func:`perspectiveTransform`,
+:func:`getAffineTransform`,
+:func:`estimateRigidTransform`,
+:func:`warpAffine`,
+:func:`warpPerspective`
 
 .. index:: transpose
 
@@ -3267,11 +3473,11 @@ transpose
 
 .. c:function:: void transpose(const Mat& src, Mat& dst)
 
-    Transposes a matrix
+    Transposes a matrix.
 
-    :param src: The source array
+    :param src: Source array.
 
-    :param dst: The destination array of the same type as  ``src``
+    :param dst: Destination array of the same type as  ``src`` .
     
 The function :func:`transpose` transposes the matrix ``src`` :
 
@@ -3279,5 +3485,6 @@ The function :func:`transpose` transposes the matrix ``src`` :
 
     \texttt{dst} (i,j) =  \texttt{src} (j,i)
 
-Note that no complex conjugation is done in the case of a complex
-matrix, it should be done separately if needed.
+**Note**:
+
+No complex conjugation is done in case of a complex matrix. It it should be done separately if needed.