|
要做的工作是将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.
不能正常转换,各位高手帮忙看看,怎么解决, |
|