查看: 1592|回复: 0

osgdem转换后影像被切割掉一部分

[复制链接]

该用户从未签到

发表于 2011-8-18 15:58:57 | 显示全部楼层 |阅读模式
要做的工作是将tif格式用vbp转换为ive格式,读取到osgearth中去,现在碰到的问题是:直接使用osgdem转换(代码为osgdem -t jilinshi_1.tif --geocentric --terrain -o jilin.ive),生成的jilin.ive加载到osgearth中,图片被切割了右边一部分。。。
        如果在osgdem中没有加入--geocentric选项,则生成的jilin.ive正常,没有被切割,但无法载入osgeargth中(应该是坐标系问题)。
        然后我用gdalinfo查看了jinlinshi_1.tif,显示如下:
E:\osgearthlearning>gdalinfo jilinshi_1.tif
Driver: GTiff/GeoTIFF
Files: jilinshi_1.tif
Size is 21993, 21082
Coordinate System is:
PROJCS["IMAGINE GeoTIFF Support
Copyright 1991 - 2006 by Leica Geosystems Geospatial Imaging, LLC. All Rights Re
served
@(#)$RCSfile: egtf.c $ IMAGINE 9.1 $Revision: 22.0 $ $Date: 2006/05/24 17:16:00
EST $
Projection = Gauss Kruger",
    GEOGCS["xian_80",
        DATUM["xian_80",
            SPHEROID["IAG 75",6378140,298.2570005923655]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",126],
    PARAMETER["scale_factor",1],
    PARAMETER["false_easting",42500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (42533371.500000000000000,4874746.000000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
  TIFFTAG_SOFTWARE=IMAGINE TIFF Support
Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved
@(#)$RCSfile: etif.c $ $Revision: 1.10.1.9.1.9.2.11 $ $Date: 2004/09/15 18:42:01
EDT $
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
  TIFFTAG_RESOLUTIONUNIT=1 (unitless)
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (42533371.500, 4874746.000) (126d24'58.05"E, 44d 0'27.12"N)
Lower Left  (42533371.500, 4864205.000) (126d24'55.66"E, 43d54'45.61"N)
Upper Right (42544368.000, 4874746.000) (126d33'11.66"E, 44d 0'25.03"N)
Lower Right (42544368.000, 4864205.000) (126d33'8.49"E, 43d54'43.52"N)
Center      (42538869.750, 4869475.500) (126d29'3.46"E, 43d57'35.39"N)
Band 1 Block=64x64 Type=Byte, ColorInterp=Gray
Band 2 Block=64x64 Type=Byte, ColorInterp=Undefined
Band 3 Block=64x64 Type=Byte, ColorInterp=Undefined

使用的应该是Gauss Kruger投影坐标系,应该要转换成地球地心坐标系(地理坐标系)才能集成到osgearth中。
然后我考虑使用gdalwrap进行转换,
E:\osgearthlearning>gdalwarp -t_srs WGS84 jilinshi_1.tif jilin.tif
Processing input file jilinshi_1.tif.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.

不能正常转换,各位高手帮忙看看,怎么解决,
您需要登录后才可以回帖 登录 | 注册

本版积分规则

OSG中国官方论坛-有您OSG在中国才更好

网站简介:osgChina是国内首个三维相关技术开源社区,旨在为国内更多的技术开发人员提供最前沿的技术资讯,为更多的三维从业者提供一个学习、交流的技术平台。

联系我们

  • 工作时间:09:00--18:00
  • 反馈邮箱:1315785073@qq.com
快速回复 返回顶部 返回列表