osg三维重建的两种方法剖析：三角面片（osgUtil::DelaunayTriangulator）和四角面片（osg::HeightField）
最近项目中需要利用osg重建三维曲面，所以学习了一下。
第一，我先用的狄洛尼三角形的方法，即osgUtil::DelaunayTriangulator，用这种方法的特点是：
1.首先必须给其一个存储三维点集的数组，该方法会对这些杂乱无章的散点自动排序，然后就利用这些排好序的，符合三角网构建规则的散点去构建三角网，需要注意的是经过dt->setInputPointArray(coords);这句话以后，数组coords的值的顺序已经发生改变，不再是原来的coords。
2.再给其贴纹理的时候，必须要首先设置一个颜色数组给它，需要注意的是，纹理的坐标都是0~1范围的，而且是二维的(x,y),所以必须将coords的坐标的x和y值一 一 映射到0~1的范围。
3.必须输出法向量，并用该法向量数组给对应的geometry赋值
具体代码如下：
//创建Delaunay三角网对象
osg::ref_ptr<osgUtil::DelaunayTriangulator> dt = new osgUtil::DelaunayTriangulator();
dt->setInputPointArray(coords);//赋给它三维点集数组
dt->setOutputNormalArray(normals);//输出法向量
//生成三角网
dt->triangulate();
osg::ref_ptr<osg::Geometry> geometry = new osg::Geometry();
//设置坐标
geometry->setVertexArray(coords.get());
//设置描述
geometry->addPrimitiveSet(dt->getTriangles());
//设置法线
geometry->setNormalArray(normals.get());
geometry->setNormalBinding(osg::Geometry::BIND_PER_PRIMITIVE);
//设置纹理坐标（纹理填充）
osg::ref_ptr<osg::Vec2Array> texCoords = ComputerTextureCoords(*(coords.get()));//得到一 一映射后的范围在0~1的二维纹理数组
geometry->setTexCoordArray(0,texCoords.get());
//尝试颜色填充
// osg::ref_ptr<osg::Vec4Array> vextexColorArray = ComputePerVertexColor(*(coords.get()),getOSGColorTable());
// geometry->setColorArray(vextexColorArray );
// geometry->setColorBinding(osg::Geometry::BIND_PER_VERTEX);
//准备纹理图像
//生成一副QImage，给其每个像素用事先定义好的颜色赋值。最后保存成png图像，保存方式就是image.save(strPath)参数是保存全文件名，包括路径和后缀。
QImage image(xCount,yCount,QImage::Format_Indexed8);
QVector<QRgb> colorTable = getColorTable();
image.setColorTable(colorTable);
InterpolateAndDrawImage(vecZs,&image,xCount,yCount,xCount,yCount);
QString strName = ::GetImagePath() + pContourData->GetName() +".png";
image.save(strName);//保存成png图像
//开始用png图像生成纹理
osg::ref_ptr<osg::Image> texImage = osgDB::readImageFile(strName.toStdString());
osg::ref_ptr<osg::Texture2D> tex = new osg::Texture2D;
tex->setImage(texImage.get());
tex->setDataVariance(osg::Object::DYNAMIC);
osg::ref_ptr<osg::StateSet> stateset = new osg::StateSet();
stateset->setTextureAttributeAndModes(0,tex.get(),osg::StateAttribute::ON);
osg::ref_ptr<osg::Geode> geode = new osg::Geode();
geode->addDrawable(geometry.get());
geode->setStateSet(stateset.get());
//设置矩阵变换矩阵
m_pRootSwitch->addChild(geode);
这样利用三角面片重建三维曲面，并且给其贴纹理渲染的效果就已经出来了，但是需要注意的是如果点集是规则网格数据，这种方式的构建就不适合了，应该用四角面片osg::HeightField的方式。
二、四边形面片
这种方式构建的特点如下：
1必须给其分配一个行列号，即构建出的网格有多少行多少列，利用allocate（unsigned int rownum,unsigned int columnnum）函数来分配。
2.必须给其指定一个初始位置，用setOrigin（(const osg::Vec3& origin）来指定，注意三维顶点数组的z轴设为0即可，点集合中最小的x、y即可。
3.必须给其指定x和y方向每行每列之间的间隔用setXInterval(float dx );setYInterval(float dy)来指定。
4.有了初始位置，行列号，和间隔，就可以算出四边形面片每个顶点的位置了，然后在每个位置上设置高程值即可，用setHeight（float height）
具体实现过程如下，纹理渲染过程与三角面片类似，只是少了一个步骤就是：不需要给其设置纹理坐标，跟简单容易。
osg::ref_ptr<osg::HeightField> pHeightField = new osg::HeightField();
pHeightField->allocate(xCount,yCount);
pHeightField->setOrigin(osg::Vec3(xMin,yMin,0));
pHeightField->setXInterval( xdelta );
pHeightField->setYInterval( ydelta );
float x,y;
for(int i = 0; i < yCount; i++)
{
for(int j = 0; j < xCount; j++)
{
x = xMin + j * xdelta;
y = yMin + i * ydelta;
double z = pInterpolater->GetInterpolatedZ(x,y,input.begin(),input.end());
vecZs.push_back(-z+zMin);
coords->push_back(GeoToGeoNormal(osg::Vec3f(x,y,0)));
pHeightField->setHeight(j,yCount-i-1,-z);//循环得到每个顶点，然后为其设置z值
}
}
///绘制纹理图像
QImage image(xCount,yCount,QImage::Format_Indexed8);
QVector<QRgb> colorTable = getColorTable();
image.setColorTable(colorTable);
InterpolateAndDrawImage(vecZs,&image,xCount,yCount,xCount,yCount);
QString strName = ::GetImagePath() + pContourData->GetName() +".png";
image.save(strName);
osg::ref_ptr<osg::Image> texImage = osgDB::readImageFile(strName.toStdString());
osg::ref_ptr<osg::Texture2D> tex = new osg::Texture2D;
tex->setImage(texImage.get());
tex->setDataVariance(osg::Object::DYNAMIC);
osg::ref_ptr<osg::StateSet> stateset = new osg::StateSet();
stateset->setTextureAttributeAndModes(0,tex.get(),osg::StateAttribute::ON);
osg::ref_ptr<osg::Geode> geode = new osg::Geode();
geode->addDrawable(new osg::ShapeDrawable(pHeightField.get()));
geode->setStateSet(stateset.get());
//附录：贴纹理用到的三个函数：
第一个将真实顶点坐标一 一映射到(0~1)的范围
osg::ref_ptr<osg::Vec2Array> COSG3DSurfaceNode::ComputerTextureCoords( const osg::Vec3Array & vp)
{
osg::ref_ptr<osg::Vec2Array> texCoords = new osg::Vec2Array();
int nSize = vp.size();
float maxX = vp[0].x();
float minX = maxX;
float maxY = vp[0].y();
float minY = maxY;
for (int i=0;i<nSize;i++)
{
maxX = maxX<vp[i].x()?vp[i].x():maxX;
minX = minX>vp[i].x()?vp[i].x():minX;
maxY = maxY<vp[i].y()?vp[i].y():maxY;
minY = minY>vp[i].y()?vp[i].y():minY;
}
for(int i = 0;i< nSize; i++)
{
float xValue = 1-(maxX-vp[i].x())/(maxX - minX);
//float yValue = 1-(maxY-vp[i].y())/(maxY - minY);
float yValue = 1-(vp[i].y()-minY)/(maxY - minY);
texCoords->push_back(osg::Vec2(xValue,yValue));
}
return texCoords;
}
第二，由颜色配置文件得到颜色数组
QVector<QRgb> COSG3DSurfaceNode::getColorTable()
{
QVector<QRgb> table;
QString colotpath = ::GetImagePath() + "colorbar.txt";
QFile file(colotpath);
if (!file.open(QIODevice::ReadOnly))
{
assert(false);
}
QTextStream WXStream(&file);
for (int i=0;i<32;i++)
{
float RColor(0.0),GColor(0.0),BColor(0.0),AColor(1.0);
osg::Vec4 Vcolor;
WXStream>>RColor>>GColor>>BColor;
table.push_back(qRgb(RColor*255,GColor*255,BColor*255));
}
return table;
}
第三、生成纹理图像，根据颜色配置文件对图像像素 一 一进行赋值
void COSG3DSurfaceNode::InterpolateAndDrawImage(const std::vector<float>& vecData,QImage* pImage,int xCount,int yCount,int imageSizeX,int imageSizeY)
{
assert(xCount > 0);
assert(yCount > 0);
assert(vecData.size() == xCount * yCount);
float min = vecData[0];
float max = vecData[0];
for(int i = 0; i < vecData.size() ; i++)
{
if(max < vecData[i]) max = vecData[i];
if(min > vecData[i]) min = vecData[i];
}
double colorFactor = getColorTable().size() / (max - min);
/// todo 图片插值
for(int i = 0; i < imageSizeX; i++)
{
for(int j = 0; j < imageSizeY; j++)
{
pImage->setPixel(i,j,(vecData[j*xCount + i] - min) * colorFactor );
}
}
pImage->save(::GetImagePath() + "123.png");
}
m_pRootSwitch->addChild(geode);
最终效果图：
