我尝试使用GpuMat
数据编写自定义内核,以查找图像像素的反余弦值。当GPU上传数据CV_8UC1
时,我可以上传,下载和更改值,但字符不能用于计算反余弦值。但是,当我尝试将GPU转换为CV_32FC1
类型(浮点数)时,在下载部分期间出现非法内存访问错误。这是我的代码:
//.cu code
#include <cuda_runtime.h>
#include <stdlib.h>
#include <iostream>
#include <stdio.h>
__global__ void funcKernel(const float* srcptr, float* dstptr, size_t srcstep, const size_t dststep, int cols, int rows){
int rowInd = blockIdx.y*blockDim.y+threadIdx.y;
int colInd = blockIdx.x*blockDim.x+threadIdx.x;
if(rowInd >= rows || colInd >= cols)
return;
const float* rowsrcptr=srcptr+rowInd*srcstep;
float* rowdstPtr= dstptr+rowInd*dststep;
float val = rowsrcptr[colInd];
if((int) val % 90 == 0)
rowdstPtr[colInd] = -1 ;
else{
float acos_val = acos(val);
rowdstPtr[colInd] = acos_val;
}
}
int divUp(int a, int b){
return (a+b-1)/b;
}
extern "C"
{
void func(const float* srcptr, float* dstptr, size_t srcstep, const size_t dststep, int cols, int rows){
dim3 blDim(32,8);
dim3 grDim(divUp(cols, blDim.x), divUp(rows,blDim.y));
std::cout << "calling kernel from func\n";
funcKernel<<<grDim,blDim>>>(srcptr,dstptr,srcstep,dststep,cols,rows);
std::cout << "done with kernel call\n";
cudaDeviceSynchronize();
}
//.cpp code
void callKernel(const GpuMat &src, GpuMat &dst){
float* p = (float*)src.data;
float* p2 =(float*) dst.data;
func(p,p2,src.step,dst.step,src.cols,src.rows);
}
int main(){
Mat input = imread("cat.jpg",0);
Mat float_input;
input.convertTo(float_input,CV_32FC1);
GpuMat d_frame,d_output;
Size size = float_input.size();
d_frame.upload(float_input);
d_output.create(size,CV_32FC1);
callKernel(d_frame,d_output);
Mat output(d_output);
return 0;
}
当我运行程序时,我的编译器告诉我:
OpenCV错误:Gpu API调用(遇到非法内存访问) 在副本,文件 /home/mobile/opencv-2.4.9/modules/dynamicuda/include/opencv2/dynamicuda/dynamicuda.hpp, 第882行终止在抛出一个实例后调用 &#39; CV ::异常&#39;什么(): /home/mobile/opencv-2.4.9/modules/dynamicuda/include/opencv2/dynamicuda/dynamicuda.hpp:882: 错误:(-217)在函数中遇到非法内存访问 复制
答案 0 :(得分:12)
您可以使用cv::cuda::PtrStp<>
或cv::cuda::PtrStpSz<>
编写自己的内核(因此您不必使用GpuMat的step-Parameter,它可以简化您的代码:D):
内核:
__global__ void myKernel(const cv::cuda::PtrStepSzf input,
cv::cuda::PtrStepSzf output)
{
int x = blockIdx.x * blockDim.x + threadIdx.x;
int y = blockIdx.y * blockDim.y + threadIdx.y;
if (x <= input.cols - 1 && y <= input.rows - 1 && y >= 0 && x >= 0)
{
output(y, x) = input(y, x);
}
}
注意:
cv::cuda::PtrStep<>
:没有尺码信息
cv::cuda::PtrStepSz<>
:尺寸信息
cv::cuda::PtrStepSzb
:对于unsigned char Mats(CV_8U)
cv::cuda::PtrStepSzf
:浮球垫(CV_32F)
cv::cuda::PtrStep<cv::Point2f>
:其他类型的示例
内核调用:
void callKernel(cv::InputArray _input,
cv::OutputArray _output,
cv::cuda::Stream _stream)
{
const cv::cuda::GpuMat input = _input.getGpuMat();
_output.create(input.size(), input.type());
cv::cuda::GpuMat output = _output.getGpuMat();
dim3 cthreads(16, 16);
dim3 cblocks(
static_cast<int>(std::ceil(input1.size().width /
static_cast<double>(cthreads.x))),
static_cast<int>(std::ceil(input1.size().height /
static_cast<double>(cthreads.y))));
cudaStream_t stream = cv::cuda::StreamAccessor::getStream(_stream);
myKernel<<<cblocks, cthreads, 0, stream>>>(input, output);
cudaSafeCall(cudaGetLastError());
}
您可以使用cv::cuda::GpuMat
调用此功能:
callKernel(d_frame, d_output, cv::cuda::Stream());
答案 1 :(得分:6)
您正在将图像step
视为float
偏移。它是从一行到下一行的字节偏移量。
尝试这样的事情:
const float* rowsrcptr= (const float *)(((char *)srcptr)+rowInd*srcstep);
float* rowdstPtr= (float *) (((char *)dstptr)+rowInd*dststep);
step - 每个矩阵行占用的字节数。
将proper cuda error checking添加到您的代码中也是一个好主意(例如func
)。您可以使用cuda-memcheck
运行代码,以查看生成无效读/写的实际内核故障。