C++实现四步相移法：光栅投影与三维重建

#define PI  3.1415926
#define step 4   //步数
#define res_x 1024  //分辨率
#define res_y 768
#define res_z 1
#define gray_max 220  //最大/最小灰度值
#define gray_min 50

void OnNewGrating()
{
// TODO: 在此添加命令处理程序代码
//this->m_ViewController->RemoveSelectedActor();
//vtkRenderer* ren = this->m_ViewController->GetRenderer();
CString Name = NULL;
CString FilePathName;
CString m_ImageFormat;

vtkSmartPointer<vtkBMPReader> reader = vtkSmartPointer<vtkBMPReader>::New();
reader->SetFileName("D:\\Documents\\Pictures\\Grating graph\\demo.bmp");  //读取一幅1024*768的BMP图片
reader->Update();

int dims[3];
reader->GetOutput()->GetDimensions(dims);

int nbOfComp;
nbOfComp = reader->GetOutput()->GetNumberOfScalarComponents();
for (int m = 0; m < 3; m++) {

for (int n = 0; n < step; n++) {

for (int k = 0; k < dims[2]; k++)
{
for (int j = 0; j < dims[1]; j++)
{
for (int i = 0; i < dims[0]; i++)
{
int fre = pow(8, m);
unsigned char * pixel =
(unsigned char *)(reader->GetOutput()->GetScalarPointer(i, j, k)); //返回第(x,y,z)个像素值的地址
*pixel = (sin(i * 2 * fre* PI / res_x + n * PI / 2) + 1)*(gray_max - gray_min) / 2 + 50;  //伽马校正，灰度范围取50~220
*(pixel + 1) = (sin(i * fre * 2 * PI / res_x + n * PI / 2) + 1)*(gray_max - gray_min) / 2 + 50;
*(pixel + 2) = (sin(i * fre * 2 * PI / res_x + n * PI / 2) + 1)*(gray_max - gray_min) / 2 + 50;
Name = "频率为";
char buf1[10];   //字符后接整形变量
itoa(fre, buf1, 10);
Name.Append(buf1);
Name.Append("，初相位为");
char buf2[10];
itoa(90 * n, buf2, 10);
Name.Append(buf2);
Name.Append("°");
}
}
}
FilePathName = ("D:\\Documents\\Pictures\\Grating graph\\");  //文件保存路径
FilePathName.Append(Name);
FilePathName.Append(".bmp");

vtkSmartPointer<vtkImageViewer2> imageViewer =
vtkSmartPointer<vtkImageViewer2>::New();