[原][osg][osgearth]简单的通过osgDB，读取高程tif，修改高程tif

 ReadResult result;

     osg::ref_ptr<osgDB::ReaderWriter> reader = osgDB::Registry::instance()->getReaderWriterForExtension("tif");

     std::string name("D:\\gd.tif");

     osgDB::ReaderWriter::Options* opt= NULL;

     osgDB::ReaderWriter::ReadResult rr = reader->readImage(name, opt);

     if (rr.validImage())

     {

         result = ReadResult(rr.takeImage());

         result.getImage()->setName("nameNoUse.tif");

     }

     if (result.succeeded())

     {

         result.getObject();

         result.metadata();

         osg::ref_ptr<osg::Image> image = result.getImage();

         osgEarth::ImageToHeightFieldConverter conv;

         osg::HeightField* hf = conv.convert(image.get());

         for (unsigned col = ; col < hf->getNumColumns(); ++col)

         {

             for (unsigned row = ; row < hf->getNumRows(); ++row)

             {

                 float height = hf->getHeight(col, row);

                 if (height < 1.0)

                 {

                     float newh = cos(height*3.141593f);

                     //float rf = rand()% 500;

                     hf->setHeight(col, row, -* newh);

                 }

                 else// if(height > 1)

                 {

                     //height = 100;//下断点看看

                 }

             }

         }

         osg::Image* newimage = conv.convert(hf);

         std::string nameofnew("D:\\gd2.tif");

         reader->writeImage(*newimage, nameofnew);

     }

如题

加SB的“平滑”功能

float fBegin = 0.1;

    //float fEnd = 0.000001;

    float fLowestValue = 1000.0;

    int fWide = 100.0;

    ReadResult result;

    osg::ref_ptr<osgDB::ReaderWriter> reader = osgDB::Registry::instance()->getReaderWriterForExtension("tif");

    std::string name("D:\\gd.tif");

    osgDB::ReaderWriter::Options* opt= NULL;

    osgDB::ReaderWriter::ReadResult rr = reader->readImage(name, opt);

    if (rr.validImage())

    {

        result = ReadResult(rr.takeImage());

        result.getImage()->setName("guandao.tif");

    }

    if (result.succeeded())

    {

        result.getObject();

        result.metadata();

        osg::ref_ptr<osg::Image> image = result.getImage();

        osgEarth::ImageToHeightFieldConverter conv;

        osg::HeightField* hf = conv.convert(image.get());

        int *fFlag = new int[hf->getNumColumns()*hf->getNumRows()];

        for (unsigned col = ; col < hf->getNumColumns(); ++col)

        {

            for (unsigned row = ; row < hf->getNumRows(); ++row)

            {

                fFlag[col*hf->getNumRows() + row] = ;

                float height = hf->getHeight(col, row);

                if (height < fBegin)

                {

                    fFlag[col*hf->getNumRows() + row] = ;

                    hf->setHeight(col, row, -fLowestValue);

                    /*

                    float newh = -1000.0;

                    if(height > 0.00001)

                        newh = 0.1 - (0.1 - height)/ (0.1-0.00001)*1000.0;

                    hf->setHeight(col, row, newh);*/

                }

            }

        }

        for (int i = ; i < hf->getNumColumns()*hf->getNumRows(); i++)

        {

            if (fFlag[i] == )//如果这值在海面以下

            {

                bool isNearSide = false;

                int nowX = i/hf->getNumRows();

                int nowY = i%hf->getNumRows();

                for (int j = ; j <= fWide; j++)

                {

                    //从离此值最近的值开始找附近的岸边，往外延伸

                    //向东南西北四个方向找，没层都遍历一圈

                    for ( int x = ;x <= j;x++ )

                    {

                        //如果找到有岸边

                        int fDifValueX = x;

                        int fDifValueY = j - x;

                        int realX = nowX - fDifValueX;

                        if (realX > )

                        {

                            int realY = nowY - fDifValueY;

                            if (realY > )

                            {

                                if (fFlag[realX*hf->getNumRows() + realY] == )//如果是岸边

                                    isNearSide = true;

                            }

                            realY = nowY + fDifValueY;

                            if (realY < hf->getNumRows())

                            {

                                if (fFlag[realX*hf->getNumRows() + realY] == )//如果是岸边

                                    isNearSide = true;

                            }

                        }

                        realX = nowX + fDifValueX;

                        if (realX < hf->getNumColumns())

                        {

                            int realY = nowY - fDifValueY;

                            if (realY > )

                            {

                                if (fFlag[realX*hf->getNumRows() + realY] == )//如果是岸边

                                    isNearSide = true;

                            }

                            realY = nowY + fDifValueY;

                            if (realY < hf->getNumRows())

                            {

                                if (fFlag[realX*hf->getNumRows() + realY] == )//如果是岸边

                                    isNearSide = true;

                            }

                        }

                    }

                    //查找这个范围内是否有值,如果有值则用此值

                    if (isNearSide)

                    {

                        float fRealHeight = fBegin - j * fLowestValue / fWide;

                        hf->setHeight((i / hf->getNumRows()), (i % hf->getNumRows()), fRealHeight);

                        break;//退出当前寻找的延伸

                    }

                }

            }

        }

        osg::Image* newimage = conv.convert(hf);

        std::string nameofnew("D:\\gd2.tif");

        reader->writeImage(*newimage, nameofnew);

        delete[]fFlag;

    }