Method interpolate

deps/bcd-1.1/src/core/MultiscaleDenoiser.cpp:473–512

	}

	void MultiscaleDenoiser::interpolate(Deepimf& o_rInterpolatedImage, const Deepimf& i_rImage)
	{
		const int width = i_rImage.getWidth();
		const int height = i_rImage.getHeight();
		const int depth = i_rImage.getDepth();
		const int upscaledWidth = o_rInterpolatedImage.getWidth();
		const int upscaledHeight = o_rInterpolatedImage.getHeight();

		assert(width == upscaledWidth / 2);
		assert(height == upscaledHeight / 2);
		assert(depth == o_rInterpolatedImage.getDepth());

		const float mainPixelWeight = 9.f / 16.f;
		const float adjacentPixelWeight = 3.f / 16.f;
		const float diagonalPixelWeight = 1.f / 16.f;

		int upscaledLine, upscaledCol, z, line, col, adjacentLine, adjacentCol;

		for(upscaledLine = 0; upscaledLine < upscaledHeight; ++upscaledLine)
			for(upscaledCol = 0; upscaledCol < upscaledWidth; ++upscaledCol)
			{
				line = upscaledLine / 2;
				col = upscaledCol / 2;
				adjacentLine = clampPositiveInteger(line + ((upscaledLine % 2) * 2 - 1), height);
				adjacentCol = clampPositiveInteger(col + ((upscaledCol % 2) * 2 - 1), width);
				for(z = 0; z < depth; ++z)
				{
					const PixelPosition p1 = i_rImage.clamp(PixelPosition(line, col));
					const PixelPosition p2 = i_rImage.clamp(PixelPosition(line, adjacentCol));
					const PixelPosition p3 = i_rImage.clamp(PixelPosition(adjacentLine, col));
					const PixelPosition p4 = i_rImage.clamp(PixelPosition(adjacentLine, adjacentCol));

					o_rInterpolatedImage.set(upscaledLine, upscaledCol, z,
							mainPixelWeight * i_rImage.get(p1, z) +
							adjacentPixelWeight * (i_rImage.get(p2, z) + i_rImage.get(p3, z)) +
							diagonalPixelWeight * i_rImage.get(p4, z));
				}
			}

	}

	void MultiscaleDenoiser::downscale(Deepimf& o_rDownscaledImage, const Deepimf& i_rImage)
	{