osgdem 导入多个DEM文件 重叠问题。

winglemon

比如附件里面是两个用gdal_translate 转的 geotiff 文件， 位置没有错误啊， 可是
我用
osgdem --xx 90 --yy 90 -d a.tiff -d b.tiff -o o.ive
用osgview只看到了一块数据。

如果分开导入
osgdem --xx 90 --yy 90 -d a.tiff  -o a.ive
osgdem --xx 90 --yy 90 -d b.tiff  -o b.ive
用 osgview a.ive b.ive 看到两个dem数据在叠罗汉， 郁闷？？


另外 arcinfo ascii 是不是osgdem不直接支持？ 只能用gdal_translate转成tiff???


多谢了！！

winglemon

2个tiff是冲arcinfo ascii转过来的
文件头分别是

ncols         976
nrows         764
xllcorner     714850
yllcorner     3414600
cellsize      25
NODATA_value  -9999


ncols         978
nrows         766
xllcorner     738725
yllcorner     3415100
cellsize      25
NODATA_value  -9999

FlySky

估计是你的DEM没有正确的坐标转换信息～～～～～～·

winglemon

搞定了，用一个shell搞定了， 原来转完geotiff之后prj文件也得方在一块。
所以干脆救灾arcinfo ascii 目录下面直接转就行了。
shell脚本如下:

#! /bin/bash
cd /media/e/demdir
for file in *.txt
do
   gdal_translate -of "GTiff" $file `basename $file .txt`.tiff
done

afilelist=""
for file in *.tiff
do
   afilelist+=" -d "$file
done
osgdem $afilelist -o chongqing_1_50000/out.ive


有个新问题。没有投影文件的dem怎么转阿？
比如我有文件头， 应该是经纬度的， 用
osgdem 怎么转阿？ 结果总是不对， 两个DEM叠在一块了

ncols         1800

nrows         1200

xllcorner     111

yllcorner     3

cellsize      0.00083333333333346

NODATA_value  -9999

FlySky

。。。。。。。。。。。。。。。。。。投影data一般都应该跟gdal放在环境路径下面～～～～dem都是有投影的，不可能是球形的吧，，经纬度可以直接使用，它的投影方式是Plate Carree～～～～～～～～

winglemon

谢谢 flysky