双立方插值实现

双立方插值实现，里面还有一个说明文档，大家可以拿去参考

// Bicubic interpolation.cpp : 定义控制台应用程序的入口点。
//

#include “stdafx.h“
#include 
#include 

double cubicInterpolate（double p[4] double x）
{
	return p[1] + 0.5 * x*（p[2] - p[0] + x*（2.0*p[0] - 5.0*p[1] + 4.0*p[2] - p[3] + x*（3.0*（p[1] - p[2]） + p[3] - p[0]）））;
}

double bicubicInterpolate（double p[4][4] double x double y）
{
	double arr[4];
	arr[0] = cubicInterpolate（p[0] y）;
	arr[1] = cubicInterpolate（p[1] y）;
	arr[2] = cubicInterpolate（p[2] y）;
	arr[3] = cubicInterpolate（p[3] y）;
	return cubicInterpolate（arr x）;
}

double tricubicInterpolate（double p[4][4][4] double x double y double z）
{
	double arr[4];
	arr[0] = bicubicInterpolate（p[0] y z）;
	arr[1] = bicubicInterpolate（p[1] y z）;
	arr[2] = bicubicInterpolate（p[2] y z）;
	arr[3] = bicubicInterpolate（p[3] y z）;
	return cubicInterpolate（arr x）;
}

double nCubicInterpolate（int n double* p double coordinates[]）
{
	assert（n > 0）;
	if （n == 1）
	{
		return cubicInterpolate（p *coordinates）;
	}
	else
	{
		double arr[4];
		int skip = 1 << （n - 1） * 2;
		arr[0] = nCubicInterpolate（n - 1 p coordinates + 1）;
		arr[1] = nCubicInterpolate（n - 1 p + skip coordinates + 1）;
		arr[2] = nCubicInterpolate（n - 1 p + 2 * skip coordinates + 1）;
		arr[3] = nCubicInterpolate（n - 1 p + 3 * skip coordinates + 1）;
		return cubicInterpolate（arr *coordinates）;
	}
}

int main（）
{
	// Create array
	double p[4][4] = { { 1334 }{ 7234 }{ 1636 }{ 2572 } };

	// Interpolate
	std::cout << bicubicInterpolate（p 0.1 0.2） << ‘\n‘;

	// Or use the nCubicInterpolate function
	double co[2] = { 0.1 0.2 };
	std::cout << nCubicInterpolate（2 （double*）p co） << ‘\n‘;
	
	system（“pause“）;
	return 0;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1780  2017-06-08 14:09  双立方插值实现\Bicubic interpolation.cpp

     文件     372474  2017-06-09 14:03  双立方插值实现\Bicubic interpolation学习总结.pdf

     目录          0  2018-07-11 09:49  双立方插值实现

----------- ---------  ---------- -----  ----

               374254                    3