第13节 实例-彩色转灰度(做假红外)

感谢

感谢网友提出这个问题。大家有问题也可以在评论区提出,有问必有答。

功能描述

网友问怎么做出来一个假红外,其实只需要一个RTT就可以了,但是网友表示用了也不行,估计是有BUG。为了模拟网友的环境,重现问题,这个是用osg341+osgEarth28,以及CompositeViewer 来做的,其它版本的osg+oe或者只使用osg和使用osgViewer应该也是可以的。就不再多进行测试了。

本节代码在网盘中:

注意:务必使用浏览器打开:
链接:https://pan.baidu.com/s/13gwJLwo_LbRnN3Bl2NXXXw
提取码:xrf5

以下是原色和转换的结果

image.png

以上是原色

以下是转换结果

image.png

关键代码

上网一查,彩色转黑白算法竟然有一大摞。灰色的图片其rgb值是一样的,比如(0.5, 0.5, 0.5)就是一张灰度图。因此由彩色转黑白关键就是由彩色的rgb算出灰度gray,然后最终的颜色就是(gray, gray , gray),RGB转gray的算法如下如下:

/*
浮点算法:Gray=R0.3+G0.59+B0.11
整数方法:Gray=(R
30+G59+B11)/100
移位方法:Gray=(R76+G151+B*28)>>8;
平均值法:Gray=(R+G+B)/3;
仅取绿色:Gray=G;
*/
上面是网上查的,下面是代码里实际用的

    "vec3 W = vec3(0.2125, 0.7154, 0.0721);\n"
    "vec4 mask = texture2D(rtt, vTexCoord);\n"
    "float luminance = dot(mask.rgb, W);\n"

其实也就相当于Gray=R0.2125+G0.7154+B*0.0721,如果效果不满意可以自己再查查相关算法。

框架

1.我们先使用了一个RTT相机,把渲染的纹理取出来,这就是正常的纹理。
view->getCamera()->attach(osg::Camera::COLOR_BUFFER, t2d);
上面就给一个相机关联了一个纹理,它的渲染结果会存到这张纹理上。
2.接着我们又创建了一个相机,这个相机会最后渲染,覆盖整个屏幕,在这个相机上画一个和屏幕一样大的方块,贴上这个纹理就可以了,如下代码geom就是这个方块。

    geom->getOrCreateStateSet()->setTextureAttributeAndModes(0, t2d, osg::StateAttribute::ON);

3.最后对这个相机应用shader,在片元shader中做如上算法的处理。

    "#version 120\n"
    "uniform sampler2D rtt;\n"
    "varying vec2 vTexCoord;\n"
    "void main(void)\n"
    "{\n"
    "vec3 W = vec3(0.2125, 0.7154, 0.0721);\n"
    "vec4 mask = texture2D(rtt, vTexCoord);\n"
    "float luminance = dot(mask.rgb, W);\n"
    "gl_FragColor = vec4(vec3(luminance), 1.0);\n"
    "}\n"

以下是全部代码

可以不用osgEarth和CompositeViewer ,用oe只是为了重现客户问题,思想都是一样的,有兴趣的快试试吧。如果修改不定,又需要纯osg的,可以在评论区留言。我再给大家改一个。

#include 
#include 
#include 
#include 
#include 
#include 

using namespace osgEarth;
using namespace osgEarth::Util;

osg::ref_ptr t2d = new osg::Texture2D;
int screenWidth = 1280;
int screenHeight = 1024;

static const char* vs = {
    "#version 120\n"
    "varying vec2 vTexCoord;\n"
    "void main(void)\n"
    "{\n"
    "    gl_Position = gl_ModelViewProjectionMatrix * gl_Vertex;\n"
    "    vTexCoord = gl_MultiTexCoord0.xy;\n"
    "}\n"
};


/*
 1.浮点算法:Gray=R*0.3+G*0.59+B*0.11
 2.整数方法:Gray=(R*30+G*59+B*11)/100
 3.移位方法:Gray=(R*76+G*151+B*28)>>8;
 4.平均值法:Gray=(R+G+B)/3;
 5.仅取绿色:Gray=G;
 */

static const char* fs = {
    "#version 120\n"
    "uniform sampler2D rtt;\n"
    "varying vec2 vTexCoord;\n"
    "void main(void)\n"
    "{\n"
    "vec3 W = vec3(0.2125, 0.7154, 0.0721);\n"
    "vec4 mask = texture2D(rtt, vTexCoord);\n"
    "float luminance = dot(mask.rgb, W);\n"
    "gl_FragColor = vec4(vec3(luminance), 1.0);\n"
    "}\n"
};

osg::Camera* createHUD()
{
    osg::Camera* camera = new osg::Camera;

    camera->setProjectionMatrix(osg::Matrix::ortho2D(0, screenWidth, 0, screenHeight));
    camera->setReferenceFrame(osg::Transform::ABSOLUTE_RF);
    camera->setViewMatrix(osg::Matrix::identity());
    camera->setClearMask(GL_DEPTH_BUFFER_BIT);
    camera->setRenderOrder(osg::Camera::POST_RENDER);
    camera->setAllowEventFocus(false);

    osg::Geode* geode = new osg::Geode();
    camera->addChild(geode);
    osg::Geometry* geom = new osg::Geometry;
    geode->addDrawable(geom);

    osg::Vec3Array* vertices = new osg::Vec3Array;
    vertices->push_back(osg::Vec3(0, 0, -0.1));
    vertices->push_back(osg::Vec3(screenWidth, 0, -0.1));
    vertices->push_back(osg::Vec3(screenWidth, screenHeight, -0.1));
    vertices->push_back(osg::Vec3(0, screenHeight, -0.1));
    geom->setVertexArray(vertices);

    osg::Vec2Array* coords = new osg::Vec2Array;
    coords->push_back(osg::Vec2(0.0, 0.0));
    coords->push_back(osg::Vec2(1.0, 0.0));
    coords->push_back(osg::Vec2(1.0, 1.0));
    coords->push_back(osg::Vec2(0.0, 1.0));
    geom->setTexCoordArray(0, coords);


    osg::Vec4Array* colors = new osg::Vec4Array;
    colors->push_back(osg::Vec4(1.0f, 1.0, 0.8f, 0.2f));
    geom->setColorArray(colors, osg::Array::BIND_OVERALL);
    geom->addPrimitiveSet(new osg::DrawArrays(GL_QUADS, 0, 4));

    geom->getOrCreateStateSet()->setMode(GL_LIGHTING, osg::StateAttribute::OFF || osg::StateAttribute::PROTECTED);
    geom->getOrCreateStateSet()->setTextureAttributeAndModes(0, t2d, osg::StateAttribute::ON);

    //写个简单shader来对图像变一下色
    osg::Program* program = new osg::Program;
    program->addShader(new osg::Shader(osg::Shader::VERTEX, vs));
    program->addShader(new osg::Shader(osg::Shader::FRAGMENT, fs));
    geom->getOrCreateStateSet()->setAttributeAndModes(program, osg::StateAttribute::ON);
    geom->getOrCreateStateSet()->addUniform(new osg::Uniform("rtt", 0));
    return camera;
}

int main()
{
    osgViewer::CompositeViewer viewer;

    osgViewer::View* view = new osgViewer::View;
    view->setName("Single view");
    view->setUpViewOnSingleScreen(0);
    view->setSceneData(osgDB::readNodeFile("D:\\code\\imp\\osgoe\\osg341oe28\\oe28\\tests\\gdal_tiff.earth"));
    view->setCameraManipulator(new EarthManipulator());
    viewer.addView(view);
    t2d->setInternalFormat(GL_RGBA);
    view->getCamera()->attach(osg::Camera::COLOR_BUFFER, t2d);

    //此时view里已经有个主相机了,加个HUD给它
    osgViewer::Viewer::Windows windows;
    viewer.getWindows(windows);
    if (windows.empty()) return 1;

    osg::Camera* hudCamera = createHUD();
    hudCamera->setGraphicsContext(windows[0]);
    hudCamera->setViewport(0, 0, windows[0]->getTraits()->width, windows[0]->getTraits()->height);

    view->addSlave(hudCamera, false);

    return viewer.run();
}

你可能感兴趣的:(第13节 实例-彩色转灰度(做假红外))