【图像处理】SFR算法详解3

本篇为SFR算法详解系列的第三篇，前面两篇为： 
《SFR算法详解1》：什么是MTF,MTF50,MTF50P 
《SFR算法详解2》：算法过程描述

本篇主要根据第二篇中的算法流程，描述前面的几个流程及其相关原理。这里参考标准中的附录C算法原理及mitre SFR的代码实现，并结合自己需要对其进行优化。

1、获取图像ROI

这里从拍摄的图像（一般为RGB）中截取ROI区域，并将其转换为Gray。

有个问题就是：用于SFR算法的图像大小多大比较合适？ROI的区域大小信息并没有严格定义，但会有影响，具体各个方面的要求可参考如下： 
对于定义ROI，有以下几个要求：（针对垂直的斜边） 
1） 对于垂直的斜边，ROI的宽度要小于高度。 
2）在ROI中，只允许出现一个斜边，并且大致居中。 
3）在底部或者顶部，比例小的那一部分不能低于5个像素 
4）在中心位置，边界到中间斜边的距离需要在20至60个像素之间 
5）高度需要在80到300之间。 
6）ROI中避免污点，保证图像边界基本连续。

而以下的图片已经很能说明问题如何解决。 

2、线性化图像

主要是对获取到的斜边数据进行伽马变换，获取到线性的图像数据。

为什么要获取线性的图像数据？ 
在SFR的标准算法中，建议的算法是用OECF将图像数据转化为线性数据。因此需要讨论下我们获取到的数据是否线性数据。线性跟非线性应该怎样去理解。

通常，sensor获取到的数据是线性的，sensor上每个pixel都是根据物体所呈现出来的颜色强度转化为对应的数值。但是，由于人眼对颜色的感知不是线性的，因此，通常在相机内部处理中，都有对转化后的图像进行伽马变换，变换后的图像就不是线性的，使得图像数据能够适应人眼。而通常取值为2.2的伽马变换。

为了获取到线性的图像数据，需要抵消这个伽马变换，在算法上，我们可以采取1/2.2的伽马变换来抵消。

在imatest的SFR算法中，也是采用的伽马变换来获得线性数据。

• The cropped image is linearized, i.e., the pixel levels are adjusted to remove the gamma encoding applied by the camera. (Gamma is adjustable with a default of 0.5). 
参考：http://www.imatest.com/docs/sharpness/#calc

在这里，我们根据前面的GRAY数据，处理获得线性化数据ϕ(p,r)，其中p表示像素的个数，r表示图像的行数。

3、找到图像中的矩心centroid

这里是为了获取图像每一行的矩心centroid，因为我们是要计算出边缘的回归线，经过矩心的回归线是最佳的。

对于图像数据，在边缘的地方，我们可以认为每一行都是边缘扩散函数的估计值，这些边缘扩散函数的估计值能够通过差分计算获得对应的线扩散函数。这里先计算每一行的矩心： 

在代码中对应的函数为：locate_centroids()，size_y为ROI的行数，shifts[]为每一行上矩心的位置，temp[]为每一行矩心距离图像中心的Y的数值。方向表示数值符号为正的方向。

•                     <li class="tool-item tool-active is-like "><a href="javascript:;"><svg class="icon" aria-hidden="true">
                          <use xlink:href="#csdnc-thumbsup"></use>
                      </svg><span class="name">点赞</span>
                      <span class="count">2</span>
                      </a></li>
                      <li class="tool-item tool-active is-collection "><a href="javascript:;" data-report-click="{&quot;mod&quot;:&quot;popu_824&quot;}"><svg class="icon" aria-hidden="true">
                          <use xlink:href="#icon-csdnc-Collection-G"></use>
                      </svg><span class="name">收藏</span></a></li>
                      <li class="tool-item tool-active is-share"><a href="javascript:;" data-report-click="{&quot;mod&quot;:&quot;1582594662_002&quot;}"><svg class="icon" aria-hidden="true">
                          <use xlink:href="#icon-csdnc-fenxiang"></use>
                      </svg>分享</a></li>
                      <!--打赏开始-->
                                              <!--打赏结束-->
                                              <li class="tool-item tool-more">
                          <a>
                          <svg t="1575545411852" class="icon" viewBox="0 0 1024 1024" version="1.1" xmlns="http://www.w3.org/2000/svg" p-id="5717" xmlns:xlink="http://www.w3.org/1999/xlink" width="200" height="200"><defs><style type="text/css"></style></defs><path d="M179.176 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5718"></path><path d="M509.684 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5719"></path><path d="M846.175 499.222m-113.245 0a113.245 113.245 0 1 0 226.49 0 113.245 113.245 0 1 0-226.49 0Z" p-id="5720"></path></svg>
                          </a>
                          <ul class="more-box">
                              <li class="item"><a class="article-report">文章举报</a></li>
                          </ul>
                      </li>
                                          </ul>
              </div>
                              <div class="right-toolbox"><a href="https://blog.csdn.net/jaych/article/details/51030939" target="_blank" class="jump-net-article">
              <svg t="1575545252354" class="icon" viewBox="0 0 1024 1024" version="1.1" xmlns="http://www.w3.org/2000/svg" p-id="5597" xmlns:xlink="http://www.w3.org/1999/xlink" width="200" height="200"><defs><style type="text/css"></style></defs><path d="M243.9 1022.2c-62.3 0-124-23.8-171.5-70.8C26.4 905.5 1.5 845 1.5 779.9s24.9-125.6 70.8-171.5l184-184.1c45.9-45.9 106.4-70.8 171.5-70.8s125.6 24.9 171.5 70.8c18.1 18.1 18.1 47 0 65.1s-47 18.1-65.1 0c-28.3-28.3-65.7-43.5-105.9-43.5s-78.1 15.3-105.9 43.7l-184 184c-58.3 58.3-58.3 153.4 0 212.3 28.3 28.3 65.7 43.5 105.9 43.5s78.1-15.3 105.9-43.5l184-184c18.1-18.1 47-18.1 65.1 0 18.1 18.1 18.1 47 0 65.1l-184 184c-46.9 48-109.1 71.2-171.4 71.2z m523.7-423l184-184c94.5-94.5 94.5-248 0-342.5s-248-94.5-342.5 0l-184 184c-18.1 18.1-18.1 47 0 65.1s47 18.1 65.1 0l184-184c28.3-28.3 65.7-43.5 105.9-43.5s78.1 15.3 105.9 43.5c58.3 58.3 58.3 153.4 0 212.3l-184 184c-58.3 58.3-153.4 58.3-212.3 0-18.1-18.1-47-18.1-65.1 0-18.1 18.1-18.1 47 0 65.1 47 47 109.3 70.8 171.5 70.8s123.9-23.2 171.5-70.8z" p-id="5598"></path></svg>
                  站内首发文章</a></div>
                          </div>
          <div class="person-messagebox">
              <div class="left-message"><a href="https://blog.csdn.net/yangpan011">
                  <img src="https://profile.csdnimg.cn/7/F/2/3_yangpan011" class="avatar_pic" username="yangpan011">
                                          <img src="https://g.csdnimg.cn/static/user-reg-year/1x/8.png" class="user-years">
                                  </a></div>
              <div class="middle-message">
                                      <div class="title"><span class="tit"><a href="https://blog.csdn.net/yangpan011" data-report-click="{&quot;mod&quot;:&quot;popu_379&quot;}" target="_blank">yangpan011</a></span>
                                          </div>
                  <div class="text"><span>发布了69 篇原创文章</span> · <span>获赞 106</span> · <span>访问量 31万+</span></div>
              </div>
                              <div class="right-message">
                                          <a href="https://im.csdn.net/im/main.html?userName=yangpan011" target="_blank" class="btn btn-sm btn-red-hollow bt-button personal-letter">私信
                      </a>
                                                          <a class="btn btn-sm attented bt-button personal-watch" data-report-click="{&quot;mod&quot;:&quot;popu_379&quot;}">已关注</a>
                                  </div>
                          </div>
                  </div>