📄 analysis.cpp

📁 图像分割算法
💻 CPP
📖 第 1 页 / 共 4 页
字号:
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return x;
	}



	Vector<ComplexDouble> IFFT(Vector<ComplexDouble>& x)
	{
		int n = x.Length();
		Vector<ComplexDouble> y(n);
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_DOUBLE, DFTI_COMPLEX, 1, n);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Backward scale
		status = DftiSetValue(DescHandle, DFTI_BACKWARD_SCALE, 1.0/(double)n);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Backward transform
		status = DftiComputeBackward( DescHandle, x.Data(), y.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}

		return y;
	}




	Vector<ComplexDouble>& IFFTI(Vector<ComplexDouble>& x)
	{
		int n = x.Length();
		//Vector<ComplexDouble> y(n);
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_DOUBLE, DFTI_COMPLEX, 1, n);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		//status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		//if(status != DFTI_NO_ERROR)
		//{
		//	StatusDisplay(status);
		//	cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
		//	Utility::RunTimeError("Problem with FFT!");
		//}
		// Backward scale
		status = DftiSetValue(DescHandle, DFTI_BACKWARD_SCALE, 1.0/(double)n);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Backward transform
		status = DftiComputeBackward( DescHandle, x.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}

		return x;
	}
















	// /////////////////
	// 2D Transform
	// /////////////////


	Matrix<ComplexFloat> FFT2(Matrix<ComplexFloat>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		Matrix<ComplexFloat> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_SINGLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Forward transform
		status = DftiComputeForward( DescHandle, x.Data(), y.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return y;
	}


	Matrix<ComplexFloat>& FFT2I(Matrix<ComplexFloat>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		//Matrix<ComplexFloat> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_SINGLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		//status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		//if(status != DFTI_NO_ERROR)
		//{
		//	StatusDisplay(status);
		//	cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
		//	Utility::RunTimeError("Problem with FFT!");
		//}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Forward transform
		status = DftiComputeForward( DescHandle, x.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return x;
	}




	Matrix<ComplexDouble> FFT2(Matrix<ComplexDouble>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		Matrix<ComplexDouble> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_DOUBLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Forward transform
		status = DftiComputeForward( DescHandle, x.Data(), y.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return y;
	}


	Matrix<ComplexDouble>& FFT2I(Matrix<ComplexDouble>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		//Matrix<ComplexDouble> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_DOUBLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		//status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		//if(status != DFTI_NO_ERROR)
		//{
		//	StatusDisplay(status);
		//	cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
		//	Utility::RunTimeError("Problem with FFT!");
		//}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Forward transform
		status = DftiComputeForward( DescHandle, x.Data() );
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return x;
	}












	Matrix<ComplexFloat> IFFT2(Matrix<ComplexFloat>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		Matrix<ComplexFloat> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_SINGLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Backward scale
		status = DftiSetValue(DescHandle, DFTI_BACKWARD_SCALE, 1.0/(double)(rows*cols));
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Commit Dfti descriptor
		status = DftiCommitDescriptor(DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}

		// Compute Forward transform
		status = DftiComputeBackward( DescHandle, x.Data(), y.Data());
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		status = DftiFreeDescriptor(&DescHandle);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			Utility::Warning("Problem when trying to free the FFT descriptor handle!");
		}


		return y;
	}


	Matrix<ComplexFloat>& IFFT2I(Matrix<ComplexFloat>& x)
	{
		int rows = x.Rows();
		int cols = x.Columns();
		//Matrix<ComplexFloat> y(rows,cols);
		int dims[2] = {cols,rows};
		int strides[3] = {0,rows,1};
		
		DFTI_DESCRIPTOR_HANDLE DescHandle = 0;
		long status = DftiCreateDescriptor(&DescHandle, DFTI_SINGLE, DFTI_COMPLEX, 2, dims);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Strides
		status = DftiSetValue(DescHandle, DFTI_INPUT_STRIDES, strides);
		if(status != DFTI_NO_ERROR)
		{
			StatusDisplay(status);
			cerr << "Line: " << __LINE__ << " File: " << __FILE__ << endl;
			Utility::RunTimeError("Problem with FFT!");
		}
		// Not Inplace
		//status = DftiSetValue(DescHandle, DFTI_PLACEMENT, DFTI_NOT_INPLACE);
		//if(status != DFTI_NO_ERROR)
		//{
		//	StatusDisplay(status);
💿 文件大小 126 K
👤 上传用户 SLing2008
📂 所属分类人工智能/神经网络
🏷️ 相关标签

#图像分割算法
⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -