文章/答案/技术大牛

发布

问从C到Neon的双线性插值
EN

Stack Overflow用户

提问于 2013-03-19 21:54:18

回答 2查看 3.5K关注 0票数 0

我正在尝试使用霓虹灯对图像进行下采样。因此，我尝试通过编写一个函数来练习霓虹灯，该函数使用霓虹灯来减去两个图像，并且我成功了。现在我回来用霓虹灯内插法写双线性插值。现在我有两个问题，从一行和一列得到4个像素，以及从4个像素计算插值(灰度)，或者从一行和一列的8个像素计算插值(灰度)。我试着去考虑这个问题，但是我认为算法应该被重写？

void resizeBilinearNeon( uint8_t *src, uint8_t *dest,  float srcWidth,  float srcHeight,  float destWidth,  float destHeight)
{

    int A, B, C, D, x, y, index;

       float x_ratio = ((float)(srcWidth-1))/destWidth ;
       float y_ratio = ((float)(srcHeight-1))/destHeight ;
       float x_diff, y_diff;

       for (int i=0;i<destHeight;i++) {
          for (int j=0;j<destWidth;j++) {
               x = (int)(x_ratio * j) ;
               y = (int)(y_ratio * i) ;
               x_diff = (x_ratio * j) - x ;
               y_diff = (y_ratio * i) - y ;
               index = y*srcWidth+x ;

               uint8x8_t pixels_r = vld1_u8 (src[index]);
               uint8x8_t pixels_c = vld1_u8 (src[index+srcWidth]);

               // Y = A(1-w)(1-h) + B(w)(1-h) + C(h)(1-w) + Dwh
               gray = (int)(
                           pixels_r[0]*(1-x_diff)*(1-y_diff) +  pixels_r[1]*(x_diff)*(1-y_diff) +
                           pixels_c[0]*(y_diff)*(1-x_diff)   +  pixels_c[1]*(x_diff*y_diff)
                           ) ;

               dest[i*w2 + j] = gray ;
           }
  }

c++

optimization

arm

neon

回答 2

Stack Overflow用户

回答已采纳

发布于 2013-03-19 22:46:34

@MarkRansom关于最近邻居与2x2双线性插值的关系不正确；使用4个像素的双线性将产生比最近邻居更好的输出。他是正确的，平均适当的像素数(如果按>2:1缩放，则大于4)将产生更好的输出。然而，霓虹灯不会对图像下采样有所帮助，除非按整数比例进行缩放。

NEON和其他SIMD指令集的最大好处是能够使用相同的操作一次处理8或16个像素。通过以您的方式访问单个元素，您将失去SIMD的所有好处。另一个问题是将数据从霓虹灯寄存器移动到ARM寄存器是一个很慢的操作。下采样图像最好由GPU或优化的ARM指令完成。

票数 1

Stack Overflow用户

发布于 2013-03-20 00:10:06

霓虹灯肯定会有助于使用双线性滤波以任意比例进行下采样。关键是vtbl.8指令的巧妙使用，它能够对来自预加载阵列的8个连续目标像素执行并行查找表：

 d0 = a [b] c [d] e [f]  g  h, d1 =  i  j  k  l  m  n  o  p 
 d2 = q  r  s  t  u  v  [w] x, d3 = [y] z [A] B [C][D] E  F ...
 d4 = G  H  I  J  K  L   M  N, d5 =  O  P  Q  R  S  T  U  V ...

可以很容易地计算括号中像素的小数位置：

 [b] [d] [f] [w] [y] [A] [C] [D],  accessed with vtbl.8 d6, {d0,d1,d2,d3}
 The row below would be accessed with            vtbl.8 d7, {d2,d3,d4,d5}

递增vadd.8 d6，d30；如果d30 =1 1 1 ... 1，则给出原点右侧像素的查找索引等。

没有理由从两行中获取像素，只是为了说明这是可能的，并且如果需要，该方法也可以用于实现轻微的失真。

在实时应用中，使用例如lanzcos可能有点夸张，但使用霓虹灯仍然是可行的。较大因子的下采样当然需要(重)滤波，但可以通过2:1的迭代平均和抽取轻松实现，并且仅在结束时使用分数采样。

对于要写入的任意8个连续像素，可以计算向量

  x_positions = (X + [0 1 2 3 4 5 6 7]) * source_width / target_width;
  y_positions = (Y + [0 0 0 0 0 0 0 0]) * source_height / target_height;

  ptr = to_int(x_positions) + y_positions * stride;
  x_position += (ptr & 7); // this pointer arithmetic goes only for 8-bit planar
  ptr &= ~7;               // this is to adjust read pointer to qword alignment

  vld1.8 {d0,d1}, [r0]
  vld1.8 {d2,d3], [r0], r2 // wasn't this possible? (use r2==stride)

  d4 = int_part_of (x_positions);
  d5 = d4 + 1;
  d6 = fract_part_of (x_positions);
  d7 = fract_part_of (y_positions);

  vtbl.8 d8,d4,{d0,d1}  // read top row
  vtbl.8 d9,d5,{d0,d1}  // read top row +1
  MIX(d8,d9,d6)             // horizontal mix of ptr[] & ptr[1]
  vtbl.8 d10,d4,{d2,d3} // read bottom row
  vtbl.8 d11,d5,{d2,d3} // read bottom row
  MIX(d10,d11,d6)           // horizontal mix of ptr[1024] & ptr[1025]
  MIX(d8,d10,d7)

  // MIX (dst, src, fract) is a macro that somehow does linear blending
  // should be doable with ~3-4 instructions

要计算整数部分，使用8.8bit分辨率就足够了(实际上不必计算666+0 1 2 3 ..7)并将所有中间结果保存在simd寄存器中。

免责声明--这是概念性的伪c/向量代码。在SIMD中，有两个并行任务需要优化:所需的最小算术运算量是多少，以及如何最小化不必要的数据混洗/复制。在这方面，与SSE相比，具有三寄存器方法的too NEON更适合于高性能DSP。第二个方面是乘法指令的数量，第三个优点是交错指令。

票数 3

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/15501429

复制

相似问题

问从C到Neon的双线性插值
EN

回答 2

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问从C到Neon的双线性插值EN

回答 2

Stack Overflow用户

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问从C到Neon的双线性插值
EN