+  int aCnt = 0;
+  float* anImageData = reinterpret_cast<float*> (theImage.ChangeData());
+
+  for (int aRow = 0; aRow < aH; aRow += 2)
+  {
+    for (int aCol = 0; aCol < aW; aCol += 3)
+    {
+      float aInvNbSamples = 1.f / values[aRow * aW + aCol + aW];
+      anImageData[aCnt++] = values[aRow * aW + aCol] * aInvNbSamples;
+      anImageData[aCnt++] = values[aRow * aW + aCol + 1] * aInvNbSamples;
+      anImageData[aCnt++] = values[aRow * aW + aCol + 1 + aW] * aInvNbSamples;
+    }
+  }