[osg-users] [vpb] Correct way to get texture data from USGS or any other source?

[Elias Tarasov]

Hello!
I try to build map using vpb in ECEF.
According to manuals i've read it needs to start: 
vpbmaster --geocentric -t texture_file -o output_file
So, clearly i need georeferenced texture file. 
On that page:
http://www.osgvisual.org/projects/osgvisual/wiki/OsgTerrainData
there is a bunch of links to get data. Since --geocentric option allows not to use elevation data, then only textures are needed.
And here is a problem: i can't get georeferenced textures from USGS. 
All files i choose from different sets there are without georeferenced information.

gdalinfo output:

c:\Terrain\1_42822\FL\2007\200701_bay_county_fl_1ft_sp_clr\vol001>gdalinfo 42822.tif
Driver: GTiff/GeoTIFF
Files: 42822.tif
       42822.tfw
Size is 5000, 5000
Coordinate System is:
LOCAL_CS["unnamed",
    UNIT["US survey foot",0.3048006096012192,
        AUTHORITY["EPSG","9003"]]]
Origin = (1655000.330000000100000,574999.510000000010000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 1655000.330,  574999.510)
Lower Left  ( 1655000.330,  569999.510)
Upper Right ( 1660000.330,  574999.510)
Lower Right ( 1660000.330,  569999.510)
Center      ( 1657500.330,  572499.510)
Band 1 Block=5000x1 Type=Byte, ColorInterp=Red
  Overviews: 1024x1024
Band 2 Block=5000x1 Type=Byte, ColorInterp=Green
  Overviews: 1024x1024
Band 3 Block=5000x1 Type=Byte, ColorInterp=Blue
  Overviews: 1024x1024

gdalwarp output:

c:\Terrain\1_42822\FL\2007\200701_bay_county_fl_1ft_sp_clr\vol001>gdalwarp -t_srs "+proj=longlat +ellps=WGS84" 42822.tif 42822_conv.tif
ERROR 1: No PROJ.4 translation for source SRS, coordinate
transformation initialization has failed.
Creating output file that is 5000P x 5000L.
Processing input file 42822.tif.
0...10...20...30...40...50...60...70...80...90...100 - done.

gdalwarp converted file output:

c:\Terrain\1_42822\FL\2007\200701_bay_county_fl_1ft_sp_clr\vol001>gdalinfo 42822_conv.tif
Driver: GTiff/GeoTIFF
Files: 42822_conv.tif
Size is 5000, 5000
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["unknown",
        SPHEROID["WGS84",6378137,298.257223563]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433]]
Origin = (1655000.330000000100000,574999.510000000010000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=1
  TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 1655000.330,  574999.510) (Invalid angle,Invalid angle)
Lower Left  ( 1655000.330,  569999.510) (Invalid angle,Invalid angle)
Upper Right ( 1660000.330,  574999.510) (Invalid angle,Invalid angle)
Lower Right ( 1660000.330,  569999.510) (Invalid angle,Invalid angle)
Center      ( 1657500.330,  572499.510) (Invalid angle,Invalid angle)
Band 1 Block=5000x1 Type=Byte, ColorInterp=Red
Band 2 Block=5000x1 Type=Byte, ColorInterp=Green
Band 3 Block=5000x1 Type=Byte, ColorInterp=Blue

Reprojection doesn't work and it seems very reasonably, because gdal doesn't know the source projection.

vpbmaster output:

c:\Terrain>vpbmaster --geocentric -t c:/terrain/42822_conv.tif -o c:/terrain/out.osgb
--geocentric
-t c:/terrain/42822_conv.tif
ADD: c:/terrain/42822_conv.tif
-o c:/terrain/out.osgb
Adding terrainTile
Error: vpbmaster can not run without all source data being in the correct destination coordinates system, please reproject them.
Recieved signal 15, doing TERMINATE_RUNNING_TASKS_THEN_EXIT.
Setting up MachinePool to use all 8 cores on this machine.
MachinePool::signal(15)
Machine::signal(15)
Machine::cancelThreads() hostname=, threads=8
  Cancel thread
  Cancel thread
  Cancel thread
  Cancel thread
  Cancel thread
  Cancel thread
  Cancel thread
  Cancel thread
Completed Machine::cancelThreads() hostname=, threads=8

BlueMarble data works, but since the best accuracy is about 250 meters, i can't create a proper detailed terrain.

What am i doing wrong?
Thanks a lot and best regards!

Elias

------------------
Read this topic online here:
http://forum.openscenegraph.org/viewtopic.php?p=63998#63998