正射投影坐标如何转为地心坐标？？？？？？

white_water125

想问下，如果想把模型叠加到地球表面上时，假如模型与地球坐标系不一致，如何转换两者的坐标，如地形模型是用VPB生成的，坐标系为投影坐标系，如何转换为地球的地心坐标系呢？

white_water125

对已生成的模型节点直接进行坐标转换，不知道可不可以，那样会不会太复杂，繁琐呢？

white_water125

对已生成的模型节点直接进行坐标转换，不知道可不可以，那样会不会太复杂，繁琐呢？

lucyli

要做的工作是将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.

不能正常转换，各位高手帮忙看看，怎么解决，

lucyli

跟楼主的问题一样，等待高手帮忙解决一下

sky11811

简单的矩阵变换是不行的，平面变为曲面，一个4*4的矩阵做不到，所以只有重新生成地形，让二者坐标系统一致。