Hi,
I also tried with dem and orthoimage downloaded at:
http://www.geobase.ca/geobase/en/index.html
I get the same problem, the projection doesn't work.
With "gdal_info data.dem", I get :
Driver: USGSDEM/USGS Optional ASCII DEM (and CDED)
Files: 013c07_0100_deme.dem
Size is 1201, 1201
Coordinate System is:
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.257222101,
AUTHORITY["EPSG","7019"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9108"]],
AXIS["Lat",NORTH],
AXIS["Long",EAST],
AUTHORITY["EPSG","4269"]]
Origin = (-60.750104166666667,52.500104166666667)
Pixel Size = (0.000208333333333,-0.000208333333333)
Metadata:
AREA_OR_POINT=Point
Corner Coordinates:
Upper Left ( -60.7501042, 52.5001042) ( 60d45'0.38"W, 52d30'0.38"N)
Lower Left ( -60.7501042, 52.2498958) ( 60d45'0.38"W, 52d14'59.63"N)
Upper Right ( -60.4998958, 52.5001042) ( 60d29'59.63"W, 52d30'0.38"N)
Lower Right ( -60.4998958, 52.2498958) ( 60d29'59.63"W, 52d14'59.63"N)
Center ( -60.6250000, 52.3750000) ( 60d37'30.00"W, 52d22'30.00"N)
Band 1 Block=1201x1201 Type=Int16, ColorInterp=Undefined
NoData Value=-32767
Unit Type: m
and with "gdal_info data.tif", I get :
Driver: GTiff/GeoTIFF
Files: S4_06032_5215_20061018_m20_1_utm20.tif
Size is 3826, 3622
Coordinate System is:
PROJCS["unnamed",
GEOGCS["NAD83",
DATUM["North_American_Datum_1983",
SPHEROID["GRS 1980",6378137,298.2572221010002,
AUTHORITY["EPSG","7019"]],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-63],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","26920"]]
Origin = (630040.000000000000000,5828260.000000000000000)
Pixel Size = (20.000000000000000,-20.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 630040.000, 5828260.000) ( 61d 4'49.40"W, 52d35'19.68"N)
Lower Left ( 630040.000, 5755820.000) ( 61d 6'29.89"W, 51d56'16.38"N)
Upper Right ( 706560.000, 5828260.000) ( 59d57'8.07"W, 52d33'54.44"N)
Lower Right ( 706560.000, 5755820.000) ( 59d59'47.44"W, 51d54'53.11"N)
Center ( 668300.000, 5792040.000) ( 60d32'3.63"W, 52d15'10.72"N)
Band 1 Block=3826x2 Type=Byte, ColorInterp=Gray
I don't understand why it doesn't work.
I tried with --geocentric and without.
I'm sending you what I obtain with "osgdem --geocentric -d data.dem -t data.tif
-l 3 -o data.osg" in attachment.
Thank you!
Cheers,
lucie
------------------
Read this topic online here:
http://forum.openscenegraph.org/viewtopic.php?p=29838#29838
Attachments:
http://forum.openscenegraph.org//files/terrain_328.jpg
_______________________________________________
osg-users mailing list
[email protected]
http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
