NDVI 计算函数：C++ 实现

该函数 CBandMath::NDVICal 用于计算 NDVI，它接受以下参数：

• m_pathname：图像文件路径。
• Bands：图像波段数。
• Width：图像宽度。
• Height：图像高度。
• DataType：数据类型，1 代表 TM 数据。

以下是函数的具体实现步骤：

1. 根据传入的 DataType 确定红波段和近红外波段的通道号。
2. 创建输出文件名为 'NDVI'。
3. 分配内存空间，用于存储计算后的 NDVI 数据 NDVIWr 以及读取红波段和近红外波段的像素值 Red 和 NearRed。
4. 打开图像文件 fReadFile 和输出文件 fWriteFile。
5. 从 fReadFile 中读取红波段和近红外波段的像素值到 Red 和 NearRed 数组中。
6. 计算 NDVI 并将结果存储到 NDVIWr 数组中，计算公式为 (NearRed[k] - Red[k]) / (NearRed[k] + Red[k])。
7. 将 NDVIWr 数组写入 fWriteFile。
8. 释放 NDVIWr, Red 和 NearRed 的内存空间。
9. 关闭 fReadFile 和 fWriteFile。

void CBandMath::NDVICal(CString m_pathname, int Bands, int Width, int Height, int DataType)
{
	int fRed, fNearRed;

	if (DataType == 1) // TM
	{
		fRed = 3;
		fNearRed = 5;
	}

	CString m_OutFile;

	float *NDVIWr = new float[Height * Width];
	unsigned char *Red = new unsigned char[Height * Width];
	unsigned char *NearRed = new unsigned char[Height * Width];

	CFile fReadFile, fWriteFile;
	if (!fReadFile.Open(m_pathname, CFile::modeReadWrite)) // 打开图像文件
	{
		AfxMessageBox('不能打开所读文件!');
		return;
	}
	// 打开写入文件
	if (!fWriteFile.Open('NDVI', CFile::modeCreate | CFile::modeWrite))
	{
		AfxMessageBox('请选择输出文件!');
		return;
	}

	fReadFile.Seek((fRed - 1) * Height * Width, CFile::begin);
	fReadFile.Read(Red, Height * Width);

	fReadFile.Seek((fNearRed - 1) * Height * Width, CFile::begin);
	fReadFile.Read(NearRed, Height * Width);

	for (int k = 0; k < Height * Width; k++)
		NDVIWr[k] = ((float)NearRed[k] - Red[k]) / (NearRed[k] + Red[k]);

	fWriteFile.Write(NDVIWr, sizeof(float) * Height * Width);

	delete[] NDVIWr;
	delete[] Red;
	delete[] NearRed;

	fReadFile.Close();
	fWriteFile.Close();

}

该函数示例代码实现了基本的 NDVI 计算功能，可以在实际应用中根据需求进行修改和完善。例如，可以加入对不同数据类型和波段组合的支持，以及对 NDVI 数据进行更精细的处理和分析。