做一个类似谷歌地球的东西

TsChus · 发表于 2010-8-10 14:38:58

我现在有几个地区的清晰的卫星地图，想做成类似谷歌地球那样的软件。就是把这些卫星图按照经纬度贴在地球表面。我刚刚接触osg,有这么几个问题。（请大家指教一下，或者告诉我看哪些例子。osg里面的例子实在太多了，虽然很经典。）

1。怎么确定球体上的经纬度，这个经纬度和显示器坐标怎么对应。
2。osg可以在物体表面指定的区域内贴图吗？怎么把卫星图贴在指定的经纬度坐标范围内。
3。因为卫星图非常大，可能导致程序运行很慢。是不是需要使用mipmap？具体怎么操作或者需要学习哪些方面内容请大家给我说一下。

谢谢各位了、·~~~~

FlySky · 发表于 2010-8-11 08:45:53

您提到的这些问题，OSG完全可以实现～～～您的意图似乎是想底层架构一个球，所以，建议您WorldWind，，跟Google Earth差不多的东西，，基于D3D的，，开源的，有C#和Java版的～先明白球的具体实现原理，再用ＯＳＧ实现的难度不大～～

TsChus · 发表于 2010-8-11 09:19:40

本帖最后由 TsChus 于 2010-8-11 11:37 编辑

2# FlySky

谢谢您的回答，另外一个问题。关于mipmap的。

我看了您在回答一个人的回帖中上传的自己写的mipmap的那个代码，我用32*32,16*16,8*8,4*4,2*2,1*1的几张图，从文件读入。已经可以了。

现在问题是，我每一级的图片都太大了，我看您说要分成小图片。我把每一级的图片都分成512*512了，于是每一级图片变成了很多张512*512的小图。但是，我不清楚接下来我怎么办？？

期待您的回答！

我想用TextureRectangle,可是平面上没有贴图，而且会有warning,说没有提供有效的TextureRectangle。

TsChus · 发表于 2010-8-11 09:31:39

本帖最后由 TsChus 于 2010-8-11 11:38 编辑

2# FlySky

谢谢flysky这些天的回答，受益匪浅！！！

FlySky · 发表于 2010-8-11 11:45:09

不太明白您的实现，，给您一个mipmap的实例：

/**********************************************************
*Write by FlySky
*zzuxp@163.com http://ww.OsgChina.org
**********************************************************/
#include <osgViewer/Viewer>
#include <osg/Node>
#include <osg/Geode>
#include <osg/Group>
#include <osg/Image>
#include <osg/TexGen>
#include <osg/Texture1D>
#include <osg/TexEnv>
#include <osg/StateSet>
#include <osgDB/ReadFile>
#include <osgDB/WriteFile>
#include <osgUtil/Optimizer>
//创建一个四边形
osg::ref_ptr<osg::Geode> createQuad()
{
osg::ref_ptr<osg::Geode> geode = new osg::Geode() ;
osg::ref_ptr<osg::Geometry> geom = new osg::Geometry() ;
geode->addDrawable(geom.get());
//设置顶点
osg::ref_ptr<osg::Vec3Array> vec = new osg::Vec3Array;
vec->push_back(osg::Vec3(-10.0f,0.0f,-10.0f));
vec->push_back(osg::Vec3(-10.0f,0.0f,10.0f));
vec->push_back(osg::Vec3(10.0f,0.0f,10.0f));
vec->push_back(osg::Vec3(10.0f,0.0f,-10.0f));
geom->setVertexArray(vec.get());
//设置法线
osg::ref_ptr<osg::Vec3Array> nor = new osg::Vec3Array;
nor->push_back( osg::Vec3f( 0.0f, -1.0f, 0.0f ) );
geom->setNormalArray(nor.get());
geom->setNormalBinding(osg::Geometry::BIND_PER_PRIMITIVE_SET);
//设置纹理坐标
osg::ref_ptr<osg::Vec2Array> tex = new osg::Vec2Array;
tex->push_back( osg::Vec2f(0.0f, 0.0f) );
tex->push_back( osg::Vec2f(0.0f, 1.0f) );
tex->push_back( osg::Vec2f(1.0f, 1.0f) );
tex->push_back( osg::Vec2f(1.0f, 0.0f) );
geom->setTexCoordArray(0,tex.get());
geom->addPrimitiveSet( new osg::DrawArrays( osg::PrimitiveSet::QUADS, 0, 4 ) );
return geode.get();
}
static void fillImage(unsigned char* ptr, unsigned int size)
{
//白色
if (size==1)
{
float r = 0.5f;
osg::Vec4 color(1.0f,1.0f,1.0f,1.0f);
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
//蓝色
if(size==2)
{
osg::Vec4 color(0.0f,0.0f,1.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//黄色
if(size==4)
{
osg::Vec4 color(0.0f,1.0f,0.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//红色
if(size==8)
{
osg::Vec4 color(1.0f,0.0f,0.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//粉红色
if(size==16)
{
osg::Vec4 color(1.0f,0.0f,1.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//黄色
if(size==32)
{
osg::Vec4 color(1.0f,1.0f,0.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//蓝绿色
if(size==64)
{
osg::Vec4 color(0.0f,1.0f,1.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//灰白色
if(size==128)
{
osg::Vec4 color(0.5f,0.5f,0.5f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
//黑色
if(size==256)
{
osg::Vec4 color(0.0f,0.0f,0.0f,1.0f);
for(unsigned int r=0;r<size;++r)
{
for(unsigned int c=0;c<size;++c)
{
*ptr++ = (unsigned char)((color[0])*255.0f);
*ptr++ = (unsigned char)((color[1])*255.0f);
*ptr++ = (unsigned char)((color[2])*255.0f);
*ptr++ = (unsigned char)((color[3])*255.0f);
}
}
}
}
int main()
{
osg::ref_ptr<osgViewer::Viewer> viewer = new osgViewer::Viewer();
osg::ref_ptr<osg::Group> root = new osg::Group();
//创建一个平面
osg::ref_ptr<osg::Geode> geode = createQuad() ;
osg::ref_ptr<osg::StateSet> stateset = new osg::StateSet();
osg::ref_ptr<osg::Image> image = new osg::Image();
//创建一个MipmapDataType列表，用来各层图片数据的偏移地址
osg::Image::MipmapDataType mipmapData;
//纹理的尺寸的最大值，必须为2的幂次
unsigned int s = 256;
//计算所需分配的数组的大小
unsigned int totalSize = 0;
for(unsigned int i=0; s>0; s>>=1, ++i)
{
if (i>0)
{
mipmapData.push_back(totalSize);
}
totalSize += s*s*4;
}
//申请一个数据
unsigned char* ptr = new unsigned char[totalSize];
//设置image的尺寸大小，数据及数据格式
image->setImage(256, 256, 256, GL_RGBA, GL_RGBA, GL_UNSIGNED_BYTE, ptr, osg::Image::USE_NEW_DELETE,1);
//将偏移地址传入imge对象
image->setMipmapLevels(mipmapData);
//向image中填充各层数据
s=256 ;
for(unsigned int i=0; s>0; s>>=1, ++i)
{
fillImage(ptr, s);
ptr += s*s*4;
}
//创建一个二维纹理对象
osg::ref_ptr<osg::Texture2D> texture = new osg::Texture2D;
//设置贴图
texture->setImage( 0, image.get() );
//设置边界处理为REPEATE
texture->setWrap(osg::Texture::WRAP_S, osg::Texture::REPEAT);
texture->setWrap(osg::Texture::WRAP_T, osg::Texture::REPEAT);
//设置滤波
texture->setFilter(osg::Texture::MIN_FILTER, osg::Texture::NEAREST_MIPMAP_NEAREST);
texture->setFilter(osg::Texture::MAG_FILTER, osg::Texture::NEAREST);
//启用二维纹理对象
stateset->setTextureAttributeAndModes( 0, texture.get(), osg::StateAttribute::ON );
geode->setStateSet(stateset.get());
root->addChild(geode.get());
//优化场景数据
osgUtil::Optimizer optimizer ;
optimizer.optimize(root.get()) ;
viewer->setSceneData(root.get());
viewer->realize();
viewer->run();
return 0 ;
}

复制代码

TsChus · 发表于 2010-8-19 12:36:06

5# FlySky

谢谢你的回答，看了你的代码mipmap的问题解决了。
可是最近做类似谷歌地球的这个东西还是没进展。world wind的源码太庞大了，我实在没有能力在短期内一个人看完。。。。。。

请教您一些问题：
1. 我用osgsimulation例子里面的createEarth函数生成了一个地球（代码如下）。
我现在有银川地区的高清卫星图片（大小为几十兆），怎么把它按照银川的坐标贴到地球上呢？图片是jpg格式的，一起下下来的还有一个jgw文件。不知道该怎么利用这些东西。

2. 看到以前您回复别人的帖子，对于大图片建议要分成小的切片。于是我把这张图片分成了很多256*256的小图片。可是接下来该怎么往地球上贴呢？

3. osgearth会不会让我的工作简化呢？

osg::Node* createEarth()
{
osg::TessellationHints* hints = new osg::TessellationHints;
hints->setDetailRatio(5.0f);

osg::ShapeDrawable* sd = new osg::ShapeDrawable(new osg::Sphere(osg::Vec3(0.0,0.0,0.0), osg::WGS_84_RADIUS_POLAR), hints);

osg::Geode* geode = new osg::Geode;
geode->addDrawable(sd);

std::string filename = osgDB::findDataFile("Images/land_shallow_topo_2048.jpg");
geode->getOrCreateStateSet()->setTextureAttributeAndModes(0, new osg::Texture2D(osgDB::readImageFile(filename)));

osg::CoordinateSystemNode* csn = new osg::CoordinateSystemNode;
csn->setEllipsoidModel(new osg::EllipsoidModel());
csn->addChild(geode);

return csn;

}

期待您的回复！

萧虓魏 · 发表于 2011-4-2 10:27:29

我也想做这个应该怎么做呢？

lanyucao · 发表于 2011-8-8 13:39:58

正在研究！

做一个类似谷歌地球的东西

联系我们