Hey Linh, Robert, I am doing a very similar thing with Moon data, and am running into problems with the WGS84 coordinate system. That is, I can't change it. Osgdem is working great for a flat surface, but if I try to make a spherical geometry I can't seem to change the radius. I have tried the radius-equator and radius-polar options, but they don't seem to change anything. I am currently going through code trying to find the problem. One thing I notice is that the --wkt and --wkt-file options in osgdem don't do anything. It looks like are just stubs in the code. I'll keep looking to see if I can figure out why my changes to the EllipsoidModel aren't working. My best guess right now is that the command line is changing one EllipsoidModel, but someplace else another EllipsoidModel is being created but is only using default values.
And Linh, I'm trying to do intersections on the Moon too like this, perhaps we could share notes when we get something working... PaulB 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
