gzr2017 / imageprocessing100wen Goto Github PK

问题八十一，第24行，多写了一个字母'j':
"在Hessian角点检测中，$\det{H}$将极大点视为j角点。"
上下文没有提到这个j，感觉应该是多余的。

jpeg图片压缩四个步骤依次为，颜色模式转换及采样、DCT变换、量化、编码。

报错信息如下，原因是y代表的y坐标超过了原图(128*128)本身的大小

解决方案：由于是运算过程中产生可比如127.5最后向上(round)取整得到128，将向上取整改为向下取整(floor)可以解决

• 向上取整：
  

• 向下取整：
  

// get and assign median
std::sort(vs, vs + (kernel_size * kernel_size));
out.atcv::Vec3b(y, x)[c] = (uchar)vs[int(floor(count / 2)) + 1];
既然取的是中值，索引从0开始，不应该+1吧

q4答案链接失效日版有，还有后面几个线性滤波的代码kernel的大小都是动态申请的（cpp），一般现在的vc应该都过不了编译把

你好，当我点击“答案”下方的链接时，
答案
Python >> answers/answer_11.py
C++ >> answers_cpp/answer_11.cpp
跳转出来的页面是“Page not found”，所以一直都是点开上面的文件夹查看答案文件的。

问题集中出现在 仿射变换（28-31）中
img[1:H+1, 1:W+1] = img 这一赋值语句，感觉是超出容限的问题，在pycharm中运行源码报错，请教下该如何解决？

你好，我写了关于这个题目系列的csdn博客文章，其中引用了你的翻译内容，希望得到许可。
还有，方便留个联系方式（email or qq）相互交流吗？

算法步骤的步骤二的第4点，代码中写的是x2，x4，x8。
说明中写的是x2，x4，x6。

答案缺失

In the description, the morphological operation is dilating before eroding. But in your code, you erode before you dilate. And we have the same problem with open operation.
File: ImageProcessing100Wen/Question_41_50/answers_py/answer_49.py, answer_50.py

line 46:val /= sqrt(height * width);
line 73: g = std::abs(val) / sqrt(height * width);
我昨天在做的时候就觉得很奇怪，怎么会能用错误的式子给出正确解答，看到answer_32.cpp就明白了...这愣是把逆变换的中的g = std::abs(val) /(height * width);拆分成了在正逆变换中分别除以sqrt(height * width)的形式...😂

  // get sum
  for (int y = 0; y < height; y++){
    for (int x = 0; x < width; x++){
      for (int c = 0; c < 3; c++){
        val = (float)img.at<cv::Vec3b>(y, x)[c];
        sum += val;
        squared_sum += (val * val);
      }
    }
  }

  // get standard deviation
  m = sum / (height * width * channel);
  s = sqrt(squared_sum / (height * width * channel) - m * m);

这一部分的均值计算是没有错的，但是，在计算标准差时，怎么能直接把所有rgb值的平方加起来算均值减去m^2？
https://pdfs.semanticscholar.org/fee0/d0c91442d465f343af18b5fce8aeff594d02.pdf
参考文章中STANDARD DEVIATION计算，应该是先计算出均值之后算每个值平方减去均值的平方再开根号。

像素间距离公式有问题

• a容易与下面的权重函数中的a混淆
• 如果我没理解错的话, ax/ay是指x/y方向上的放大倍数, 公式像a*x/a*y, 尤其是在渲染后的pdf中

1. 使用仿射变换，输出（1）那样的$x$轴倾斜$30$度的图像（ $d_x=30$），这种变换被称为X-sharing。
2. 使用仿射变换，输出（2）那样的y轴倾斜$30$度的图像（ $d_y=30$），这种变换被称为Y-sharing。
3. 使用仿射变换，输出（3）那样的$x$轴、$y$轴都倾斜$30$度的图像($d_x = 30$，$d_y = 30$)。

以及底下的两个公式，t_x，t_y都应该替换为d_x和d_y

平均池化的输入矩阵大小和输出矩阵大小应该不是一样大，为什么代码中要将均值化后的值赋给每一个像素？

Thank you for your great work.

_{Sent with GitHawk}

打扰了
我想和yoyoyo-yo/Gasyori100knock保持同步的同时忽略掉README文件。但是如果直接merge好像不能忽略掉这些文件，也会被同步成日语。

大家有什么好方法能保持代码同步的同时忽略掉翻译好的REAME文件吗？

Affine Transformations.
很棒的教程！

首先谢谢楼主的翻译!

首先定义p_low为值大于小阈值的像素, p_high为>大阈值点的像素.

楼主实现的canny是: p_low属于edge, 仅当他周围有一个p_high.

但实际上的算法应该是: p_low属于edge, 当他通过一系列p_low连接到了某个p_high. (所以会比楼主的算法产生更多的边缘点). 然后这应该是用dfs实现的吧.

I cant`t run this anser even at colaboratory.
error below:img[1:H+1, 1:W+1] = _img
NameError: name '_img' is not defined

两个min()里面应该减2不然会越界，以及某大神pr的公式错误

我注意到在Adobe Photoshop较新版本中，有类似的功能，名字译为“色调分离”，不知是否具有参考性？

def otsu_binarization(img, th=128):
    max_sigma = 0
    max_t = 0
    
    # determine threshold
    for _t in range(1, 255):
        v0 = out[np.where(out < _t)]
        m0 = np.mean(v0) if len(v0) > 0 else 0.
        w0 = len(v0) / (H * W)
        v1 = out[np.where(out >= _t)]
        m1 = np.mean(v1) if len(v1) > 0 else 0.
        w1 = len(v1) / (H * W)
        sigma = w0 * w1 * ((m0 - m1) ** 2)
        if sigma > max_sigma:
            max_sigma = sigma
            max_t = _t

    # Binarization
    print("threshold =", max_t)
    th = max_t
    out[out < th] = 0
    out[out >= th] = 255

    return out

作者v0 = out[np.where(out < _t)]中的out可能是指函数传入参的img吧，不是太理解

在文件https://github.com/gzr2017/ImageProcessing100Wen/blob/master/Question_01_10/answers_cpp/answer_9.cpp

for (int y = 0; y < height; y++){
    for (int x = 0; x < width; x++){
      for (int c = 0; c < channel; c++){

      v = 0;

      for (int dy = -pad; dy < pad + 1; dy++){
        for (int dx = -pad; dx < pad + 1; dx++){
          if (((x + dx) >= 0) && ((y + dy) >= 0)){
            v += (double)img.at<cv::Vec3b>(y + dy, x + dx)[c] * kernel[dy + pad][dx + pad];
          }
        }
      }
      out.at<cv::Vec3b>(y, x)[c] = v;
      }
    }
  }

这个

if (((x + dx) >= 0) && ((y + dy) >= 0))

是不是有问题。

获取edge时使用的fx、fy为uint8格式，是否会导致fx、fy计算edge时出现数据溢出反而导致edge的数值偏小？

404 not found

膨胀和腐蚀应该是对out进行的而不是mor把，否则N>1时结果一样

line 41:
out.atcv::Vec3b(y, x)[_c] = (uchar)((b - a) / (d - c) * (val - c) + a);

这里先计算(b - a) / (d - c)会损失数据，导致结果错误，若b-a=255, d-c=178, 斜率会得出1，结果错误。

应该改成(b-a)*(val-c)/(d-c)+a