#include?“iostream“

double**?Jconv（double**?featurePixeldouble**?featurePixel2int?featureSizeYint?featureSizeX）
{
	int?Heigth?=?2?*?featureSizeY?-?1?;
	int?Width??=?2?*?featureSizeX?-?1?;
	double?**tempPixel;
	tempPixel?=?new?double*[Heigth];
	for?（int?a=0?;?a?	{
		tempPixel[a]?=?new?double[Width]?;
	}
	double?**finalPixel;
	finalPixel?=?new?double*[featureSizeY];
	for?（int?a=0?;?a?	{
		finalPixel[a]?=?new?double[featureSizeX]?;
	}
	//初始化數組
	for（int?i?=?0?;?i?	{
		for?（int?j?=?0?;?j?		{
			tempPixel[i][j]?=?0;
		}
	}
	for（int?i?=?0?;?i?	{
		for?（int?j?=?0?;?j?		{
			finalPixel[i][j]?=?0;
		}
	}

	int?k1?=?（Heigth-1）/2;
	int?k2?=?（Width-1）/2;

	for?（int?m?=?0;?m?<=?k1;?m++）
	{
		for（int?n?=?0?;?n<=?k2;?n++）
		{
			for?（int?i?=?0;?i<=?m;i++）
			{
				for（int?j?=?0;?j<=n;?j++）
				{
					tempPixel[m][n]?=?tempPixel[m][n]?+?featurePixel[i][j]?*?featurePixel2[m-i][n-j]?;
				}
			}
		}
		for（int?n?=?k2+1?;?n		{
			for?（int?i?=?0;?i<=?m;i++）
			{
				for（int?j?=?n;?j				{
					tempPixel[m][n]?=?tempPixel[m][n]?+?featurePixel[i][j-k2]?*?featurePixel2[m-i][n-j+k2]?;
				}
			}
		}
	}

	for?（int?m?=?k1+1;?m?	{
		for（int?n?=?0?;?n<=?k2;?n++）
		{
			for?（int?i?=?m;?i			{
				for（int?j?=?0;?j<=n;?j++）
				{
					tempPixel[m][n]?=?tempPixel[m][n]?+?featurePixel[i-k1][j]?*?featurePixel2[m-i+k1][n-j]?;
				}
			}
		}
		for（int?n?=?k2+1?;?n		{
			for?（int?i?=?m;?i			{
				for（int?j?=?n;?j				{
					tempPixel[m][n]?=?tempPixel[m][n]?+?featurePixel[i-k1][j-k2]?*?featurePixel2[m-i+k1][n-j+k2]?;
				}
			}
		}
	}
	std::cout?<	for（int?i?=?0?;?i?	{
		for?（int?j?=?0?;?j?		{
			std::cout?<		}
		std::cout?<	}

	for（int?i?=?0?;?i?	{
		for?（int?j?=?0?;?j?		{
			finalPixel[i][j]?=?tempPixel[i+（k1+1）/2][j+（k2+1）/2];
		}
	}

	return?finalPixel;
}

int?main（）
{
	std::cout?<	int?serachWidth?=?1;
	int?searchHeight?=?1;
	int?featureSizeX?=?2*serachWidth+1;
	int?featureSizeY?=?2*searchHeight+1?;
	//提取局部圖像數組
	double?**featurePixel;
	featurePixel?=?new?double*[featureSizeY];
	for?（int?a=0?;?a?	{
		featurePixel[a]?=?new?double[featureSizeX]?;
	}
	double?**featurePixel2;
	featurePixel2?=?new?double*[featureSizeY];
	for?（int?a=0?;?a?	{
		featurePixel2[a]?=?new?double[featureSizeX]?;
	}
	//初始化數組
	for（int?i?=?0?;?i?	{
		for?（int?j?=?0?;?j?		{
			featurePixel[i][j]?=?3*i+j+1;
			std::cout?<		}
		std::cout?<	}
	for（int?i?=?0?;?i?