Hi Robert,

 I didn't use any options except:

osgdem -d DEM_1m_VictoriaCrater_1.tif -l 4 -o test.ive -a test.osga

DEM_1m_VictoriaCrater_1.tif was generated using OpenEV from the data set on:
http://webgis.wr.usgs.gov/download/VictoriaCrater_HiRISE/VictoriaCrater_HiRISE.zip

After you unzip it, you will find under the directory VictoriaCrater_HiRISE/GIS_EQUI360_OC:
DEM_1m_VictoriaCrater.tif
DEM_1m_VictoriaCrater.prj
DEM_1m_VictoriaCrater.rrd
DEM_1m_VictoriaCrater.aux
DEM_1m_VictoriaCrater.asc.gz

I then open DEM_1m_VictoriaCrater.tif using the OpenEV tool, and export DEM_1m_VictoriaCrater_1.tif as GTiff. I don't know if going through OpenEV is neccessary but when I ran

osgdem -d DEM_1m_VictoriaCrater.tif -l 4 -o test.ive -a test.osga

on the original DEM_1m_VictoriaCrater.tif file, osgdem died with:

ERROR 6: Unable to load PROJ.4 library (libproj.so), creation of
OGRCoordinateTransformation failed.

Thanks alot Robert,

Linh
Robert Osfield wrote:

Hi Linh,

VPB/osgdem decorates databases with a CoordinateSystemNode which will
record the cooridnate system it uses.  VPB does default to using WGS84
when working geocentric, but you haven't mentioned any of the osgdem
options so I can't really say whether this will have been applied.

I must admit I have yet to work with anything other than earth bound
imagery/dem's so have not been able to test it out on other planetary
data.  Would it be possible to provide me with some?

Robert.

On Thu, Apr 24, 2008 at 3:40 AM, Linh Phan <[EMAIL PROTECTED]> wrote:
Hi,

I have GeoTIFF file with a DEM_1m_VictoriaCrater.prj:

PROJCS["Plate_Carree",GEOGCS["EQUIRECTANGULAR_CYLINDRICAL_MARS",DATUM["D_MARS",
 
SPHEROID["MARS",3396190,0]],PRIMEM["Reference_Meridian",0],UNIT["Degree",0.017453292519943295]],
PROJECTION["Plate_Carree"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",180],
PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

osgdem successfully generate a dem for the DEM_1m_VictoriaCrater.tif.
I  was wondering if osg's getWorldIntersectPoint:

intersectPoint=hlist.front().getWorldIntersectPoint();

will return the intersectPoint in the Mars body frame, ie will the norm
of "intersectPoint._v" be around 3396190  meters, the Mars radius?

Thanks alot,

Linh


_______________________________________________
osg-users mailing list
[email protected]
http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org

_______________________________________________
osg-users mailing list
[email protected]
http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org

_______________________________________________
osg-users mailing list
[email protected]
http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org

Reply via email to