osgDEM无法生成谷歌球体？

sky11811

osgDEM无法生成谷歌球体？
将坐标系统指定为mercator后，osgDEM并不报错，能运行，但是cpu一直占用很高且不生成任何东西。
何故？是否只能在wgs84下运行？

sky11811

顶顶。

FlySky

没啥问题，建议先处理好文件坐标系

sky11811

谢谢回复，我只提供一了张全球卫星影像，大小为2048*2048，谷歌墨卡托投影。
我的命令如下：
osgdem -t L:\tex\qiuMerca_t.tif --geocentric --cs epsg:3785 -l 4 -o L:\tex\sky.ive
运行后不报错，dos窗口的提示信息也没有提示有警告或错误，但是就是不生成任何东西，且CPU占用很高。
命令有问题吗？请指教。

sky11811

下面是运行结果。
请帮忙分析下，我正在谋划vpb教程，需要收集尽可能多的实践经验。
Microsoft Windows XP [版本 5.1.2600]
(C) 版权所有 1985-2001 Microsoft Corp.

C:\Documents and Settings\Administrator>osgdem -t L:\tex\qiuMerca_t.tif --geocen
tric --cs epsg:3785 -l 4 -o L:\tex\sky.ive
--geocentric
-t L:\tex\qiuMerca_t.tif
--cs "epsg:3785" converted to PROJCS["Popular Visualisation CRS / Mercator (depr
ecated)",GEOGCS["Popular Visualisation CRS",DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,AUTHORITY["EPSG","7059"]],TOWG
S84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6055"]],PRIMEM["Greenwich",0,AUTHORITY["EPS
G","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY
["EPSG","4055"]],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARA
METER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",
0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],EXTEN
SION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +
y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"],AUTHORITY["EPSG","3785
"]]
-o L:\tex\sky.ive
Adding terrainTile
DataSet::_run() 0 0
started DataSet::createDestination(4)
AR=1.000000 C1=1 R1=1
createNewDestinationGraph
Time for _destinationGraph->computeMaximumSourceResolution() = 0.006171
Time for createDestinationGraph 0.007761
Time for after_computeNeighbours 0.000189
Time for consolodateRequiredResolutions 0.000002
Time for after_reproject 0.000663
Time for after_sort 0.000187
completed DataSet::createDestination(4)
There are 1 contributing source files:
    L:\tex\qiuMerca_t.tif
mkpath(L:\tex\)
Need to create output task directory = L:\tex\sky_root_L0_X0_Y0
mkpath(L:\tex\sky_root_L0_X0_Y0)
   created directory L:\tex\sky_root_L0_X0_Y0
Task output directory = L:\tex\sky_root_L0_X0_Y0\
started DataSet::writeDestination(L:\tex\sky.ive)
_readRow 1
   reading tile level=0 X=0 Y=0
imageName = sky_L0_X0_Y0.dds
DestinationTile::readFrom(SetName=, FileName=L:\tex\qiuMerca_t.tif)
_equalizeRow 1
   equalizing tile level=0 X=0 Y=0
_writeRow 1


到这里就一直卡住了。

sky11811

我刚又下载了最新的vpb并编译，现在可以生成了，但是结果很诡异。
应该还是坐标系统的问题，请大大们帮助。

sky11811

我刚又下载了最新的vpb并编译，现在可以生成了，但是结果很诡异。命令：osgdem -t L:\tex\qiuMerca_t.tif --geocentric --cs epsg:3785 -l 4 -o L:\tex\sky.ive
应该还是坐标系统的问题，请大大们帮助。

就是这个结果

sky11811

如果命令为 osgdem -t L:\tex\qiuMerca_t.tif --geocentric --cs epsg:900913 -l 4 -o L:\tex\sky.ive
则结果如下图：

sky11811

如果命令为 osgdem -t L:\tex\qiuMerca_t.tif --geocentric --cs epsg:4326 -l 4 -o L:\tex\sky.ive
则结果正确，如下图：

sky11811

但是问题是，如果使用epsg:4326 则我的原始影像会被重新投影，且地球两极会有黑洞，因为墨卡托投影的纬度范围是大约正负85度。

sky11811

如果使用osgEarth，同样一副影像，使用Profile: global-mercator，最终结果则会完全没有问题。如图

osgEarth结果

sky11811

vpb有没有和osgEarth的global-mercator一样效果的坐标系统。

sky11811

冰天雪地裸体跪求。

cooloboy

怎么 图片 都看不见？？？？

xairwolfcn

Google 坐标系是  EPSG900913！ 不光是墨卡托 ！

xairwolfcn

在global mapper里我是这样设置的 ，生成的三维地形精度相当高！

sky11811

global mapper也能用来生成三维地形？好像是重投影工具吧。
EPSG900913我自然知道，不过现在已经是epsg:3785了。而且二者我都试过了，一样的。
虽然没有对我的问题有帮助，但还是感谢回复。