面向灰度的直方圖

Question

對於一個項目，我正在編寫一些代碼來計算一些圖像的灰度值，但我堅持使用atan2函數時我的方向只在0〜90度之間。 我猜測，出現這種問題是由於的OpenCV的filter2D功能，但我不知道這是什麼原因或者說我在做別的事情不對：

Vector<Vector<Mat_<float>>> HoG(Mat image) { 
Mat img_x; 
Mat img_y; 
IplImage img = image; 

Mat kern_x = (Mat_<char>(1, 3) << -1, 0, 1); 
Mat kern_y = (Mat_<char>(3, 1) << -1, 0, 1); 

filter2D(image, img_x, image.depth(), kern_x); 
filter2D(image, img_y, image.depth(), kern_y); 

Vector<Vector<Mat_<float>>> histograms; 

for(int y = 0; y < image.rows - size; y += size) { 
    Vector<Mat_<float>> temp_hist; 
    for(int x = 0; x < image.cols - size; x += size) { 
     float total_mag = 0; 
     Mat hist = Mat::zeros(1, 8, CV_32FC1); 
     for(int i = y; i < y + size; ++i) { 
      for(int j = x; j < x + size; ++j) { 
       float grad_x = (float)img_x.at<uchar>(i, j); 
       float grad_y = (float)img_y.at<uchar>(i, j); 
       double ori = myatan2(grad_x, grad_y); 
       float mag = sqrt(pow(grad_x, 2) + pow(grad_y, 2)); 
       int bin = round(ori/45); 
       hist.at<float>(0, (bin - 1 < 0 ? 7 : bin - 1)) += - (float)(ori - ((round(ori/45) - 1) * 45.0 + 22.5))/45.0f; 
       hist.at<float>(0, bin) += -(float)(ori - ((round(ori/45) - 1) * 45.0 + 22.5))/45.0f; 
       total_mag += mag; 
      } 
     } 
     // Normalize the histogram 
     for(int i = 0; i < 8; ++i) { 
      hist.at<float>(0, i) = hist.at<float>(0, i)/total_mag; 
     } 
     temp_hist.push_back(hist); 
    } 
    histograms.push_back(temp_hist); 
} 

return histograms; 
} 

如果您有任何其他提示在我的代碼或其他東西中加快速度，當然也是受歡迎的。

Answer 1

我注意到這一點：

float grad_x = (float)img_x.at<uchar>(i, j); 
float grad_y = (float)img_y.at<uchar>(i, j); 

你似乎使用uchar。這不應該是char？