# Generates a surface
**URL:** https://discourse.vtk.org/t/generates-a-surface/11238
**Category:** Support
**Created:** [April 18, 2023, 5:38am UTC](https://discourse.vtk.org/t/generates-a-surface/11238 "2023-04-18T05:38:40Z")
**Posts on this page:** 1
**Page:** 1
### Author:  [@Frankenstein](https://discourse.vtk.org/u/Frankenstein)
#### Post date: [April 18, 2023, 5:38am UTC](https://discourse.vtk.org/t/generates-a-surface/11238/1 "2023-04-18T05:38:40Z")
I’m creating a surface using a point cloud, and I want to control the color of each point.Just like this picture，
I have used Class vtkSurfaceReconstructionFilte to create surface,but I can’t control the color.I want more advice, thanks
here is my code
#ifndef INITIAL\_OPENGL
#define INITIAL\_OPENGL
#include \
VTK\_MODULE\_INIT(vtkRenderingOpenGL2)
VTK\_MODULE\_INIT(vtkRenderingVolumeOpenGL2)
VTK\_MODULE\_INIT(vtkInteractionStyle)
VTK\_MODULE\_INIT(vtkRenderingFreeType)
#endif
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include \
#include
#define pi 3.1415926535
#include \
#include \
#include \
#include \ //点数据
#include \
#include \
int main()
{
vtkOutputWindow::SetGlobalWarningDisplay(0);
```
std::vector>XX;
std::vector X;
std::vector>YY;
std::vector Y;
std::vector>ZZ;
std::vector Z;
double wzj=180
;double r=100;
int m=73;int n=37;
```
/\* for(int i=0;i\<73;i++)
{
X.push\_back(cos(((180-wzj)/2+180+wzj_i/(m-1))pi/180)r);
Y.push\_back(sin(((180-wzj)/2+180+wzji/(m-1))pi/180)r);
}
for(int i=0;i\<37;i++)
{
XX.push\_back(X);
YY.push\_back(Y);
}/
std::ifstream ifs\_p1;
ifs\_p1.open(“C:/Rbjse\_work/untitled/module0/result\_p1.txt”);
std::vector\pp1;
std::vector pp1\_tmp;
double tmp;
for(int i=0;i\<72+1;i++)
{
for(int j=0;j\<36+1;j++)
{
ifs\_p1\>\>tmp;
pp1\_tmp.push\_back(tmp500);
}
pp1.push\_back(pp1\_tmp);
pp1\_tmp.clear();
}
for(int i=0;i\<37;i++)
{
for(int j=0;j\<73;j++)
{
X.push\_back(cos(((180-wzj)/2+wzjj/(m-1))pi/180)(r+pp1[j][i]));
Y.push\_back(sin(((180-wzj)/2+wzj_j/(m-1))_pi/180)_(r+pp1[j][i]));
}
XX.push\_back(X);
YY.push\_back(Y);
X.clear();
Y.clear();
}
```
for(int i=0;i<37;i++)
{
for(int j=0;j<73;j++)
{
Z.push_back(i*100/37);
}
ZZ.push_back(Z);
Z.clear();
}
vtkSmartPointer pts =
vtkSmartPointer::New();
for(int i=0;i<37;i++)
{
for(int j=0;j<73;j++)
{
pts->InsertNextPoint(XX[i][j],ZZ[i][j],YY[i][j]);
}
}
vtkSmartPointer scalars_surf = vtkSmartPointer::New(); // 创建存储顶点属性的float数组
scalars_surf->SetNumberOfValues(n*m);
for (int i = 0; i < n*m; i++) // 为属性数组中的每个元素设置标量值（这个标量值可以当作颜色值）
{
scalars_surf->SetValue(i, i);
}
vtkSmartPointer hueLut = vtkSmartPointer::New();
hueLut->SetNumberOfColors(n*m); // 指定颜色查找表中有多少种颜色
hueLut->SetHueRange(0.6667, 0.0); //蓝到红渐变
hueLut->Build();
vtkSmartPointer points = vtkSmartPointer::New();
points->SetPoints(pts); //获得网格模型中的几何数据：点集
points->GetPointData()->SetScalars(scalars_surf);
vtkSmartPointer surf =
vtkSmartPointer::New();
surf->SetInputData(points);
surf->SetNeighborhoodSize(10);
surf->SetSampleSpacing(2);
surf->Update();
vtkSmartPointer contour =
vtkSmartPointer::New();
contour->SetInputConnection(surf->GetOutputPort());
//contour->SetValue(0, 0.0);
contour->Update();
//
vtkSmartPointer vertexGlyphFilter =
vtkSmartPointer::New();
vertexGlyphFilter->AddInputData(points);
vertexGlyphFilter->Update();
vtkSmartPointer contourMapper =
vtkSmartPointer::New();
contourMapper->SetInputData(contour->GetOutput());
contourMapper->SetScalarRange(0, n*m); // 设置标量值的范围（0-100）
contourMapper->ScalarVisibilityOn();
contourMapper->SetColorModeToDefault(); // 默认的映射器行为，即把unsigned char类型的标量属性数据当作颜色值，不执行隐式。对于其他类型的标量数据，将通过查询表映射。
contourMapper->SetLookupTable(hueLut);
vtkSmartPointer contourActor =
vtkSmartPointer::New();
contourActor->SetMapper(contourMapper);
/// double pointView[4] = { 0, 0, 0.5, 1 };
double contourView[4] = { 0.5, 0, 1, 1 };
vtkSmartPointer contourRender =
vtkSmartPointer::New();
contourRender->AddActor(contourActor);
contourRender->SetViewport(contourView);
contourRender->SetBackground(1, 1, 1);
vtkSmartPointer rw =
vtkSmartPointer::New();
rw->AddRenderer(contourRender);
rw->SetSize(640, 320);
rw->SetWindowName("3D Surface Reconstruction ");
rw->Render();
vtkSmartPointer rwi =
vtkSmartPointer::New();
rwi->SetRenderWindow(rw);
rwi->Initialize();
rwi->Start();
return 0;
```
}
here is my result
