AvicMat 笛卡尔坐标系转换为极坐标系的算法实现

该算法实现了将二维平面上的笛卡尔坐标系转换为极坐标系的功能。输入参数'x'和'y'分别表示二维平面上的x坐标和y坐标，'magnitude'和'angle'分别为输出参数，用于保存转换后的极坐标的大小和角度信息。'angleInDegrees'参数用于指定是否将角度以度数的形式表示。

算法流程如下：

1. 获取输入矩阵'x'和'y'的通道数'cn'，以及四个指针'ptrs[4]'，分别指向'x'、'y'、'magnitude'和'angle'矩阵的数据。
2. 定义两个缓冲区'buf[2]'，用于存储计算过程中的临时数据。
3. 计算输入矩阵的大小'size'和总元素数'total'。
4. 定义循环变量'i'和'j'，循环变量'i'用于控制并行计算的次数，循环变量'j'用于遍历每个输入矩阵的元素。
5. 在并行区域中，对于每个输入矩阵的元素，将'x'和'y'指针转换为float类型的指针，并获取'mag'和'angle'指针。
6. 调用'avic_magnitude_32f'函数计算输入矩阵的大小，并将结果保存到'mag'指针指向的内存中。
7. 调用'avic_fastAtan2_32f'函数计算输入矩阵的角度，并将结果保存到'angle'指针指向的内存中。
8. 更新指针'ptrs[0]'、'ptrs[1]'、'ptrs[2]'和'ptrs[3]'，将它们分别指向下一个元素的位置。
9. 循环结束后，返回1表示算法执行成功。

int avic_cartToPolar(const AvicMat* x, const AvicMat* y, AvicMat* magnitude, AvicMat* angle, bool angleInDegrees)
{

    int cn = AVICCV_MAT_CN(x->type);
    unsigned char* ptrs[4] = {x->data.ptr, y->data.ptr, magnitude->data.ptr, angle->data.ptr };

    float* buf[2] = { 0, 0 };
    size_t size = x->cols * x->rows;
    int j, k, total = (int)(size * cn), blockSize = AVICCV_MIN(total, ((AVIC_BLOCK_SIZE + cn - 1) / cn) * cn);
    size_t esz1 = AVICCV_ELEM_SIZE1(x->type);
    size_t i;
#ifdef WITH_OMP
    #pragma omp parallel for private(j) schedule(guided ,8)
#endif
    for (i = 0; i < 1; i++)
    {
        for (j = 0; j < total; j += blockSize)
        {
            int len = AVICCV_MIN(total - j, blockSize);

            const float* x = (const float*)ptrs[0], * y = (const float*)ptrs[1];
            float* mag = (float*)ptrs[2], * angle = (float*)ptrs[3];
            avic_magnitude_32f(x, y, mag, len);
            avic_fastAtan2_32f(y, x, angle, len, angleInDegrees);
            

            ptrs[0] += len * esz1;
            ptrs[1] += len * esz1;
            ptrs[2] += len * esz1;
            ptrs[3] += len * esz1;
        }
    }

	return 1;
}