Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Sebastian, Personally I'd subclass from osgTerrain::GeometryTechnique and pass a FloatArray as a vertex attrib array or tex coord array and store the heights there. Robert. On 4 December 2013 13:23, Sebastian Messerschmidt sebastian.messerschm...@gmx.de wrote: Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Sabastian, osgEarth does more or less what Robert describes. We encode each vert's height value (and unit vector) in a vertex attribute. On Dec 6, 2013 5:43 AM, Robert Osfield robert.osfi...@gmail.com wrote: Hi Sebastian, Personally I'd subclass from osgTerrain::GeometryTechnique and pass a FloatArray as a vertex attrib array or tex coord array and store the heights there. Robert. On 4 December 2013 13:23, Sebastian Messerschmidt sebastian.messerschm...@gmx.de wrote: Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Am 06.12.2013 13:14, schrieb Glenn Waldron: Sabastian, osgEarth does more or less what Robert describes. We encode each vert's height value (and unit vector) in a vertex attribute. Thank you Glenn. Are the vertices relative to the Ellipsoid? Or can you point me to the sources where I can check it myself? cheers Sebastian On Dec 6, 2013 5:43 AM, Robert Osfield robert.osfi...@gmail.com mailto:robert.osfi...@gmail.com wrote: Hi Sebastian, Personally I'd subclass from osgTerrain::GeometryTechnique and pass a FloatArray as a vertex attrib array or tex coord array and store the heights there. Robert. On 4 December 2013 13:23, Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de wrote: Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Robert, Hi Sebastian, Personally I'd subclass from osgTerrain::GeometryTechnique and pass a FloatArray as a vertex attrib array or tex coord array and store the heights there. Okay, that would double the data then. I'm a bit nervous about adding more vertex attributes as I already have a lot of data per vertex. Is there any option where I could pass the ModelMatrix (not the ModelViewMatrix) in a draw callback or something? The thing is that I might need more than simply the height, but probably the complete world-space vertex. As I wrote before I guess the precision is lost when using the osg_ViewMatrixInverse * ModelVIewMatrix, so moving this calculation to the CPU would help a lot. Robert. On 4 December 2013 13:23, Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de wrote: Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi, I'm facing the same problem, and I can't pre-compute vertex attributes on CPU because my vertex are moved by the shader him-self (displacement mapping). I've tried to use double matrices, but these leads to a lot of compatibility problem. Since the precision problem is in the temporary computation and not in the final result (result values are small and fit well in a 32 bits float) I precompute (on CPU) lot of values in a 2048 (or 4096) pixels wide Texture1D. In the shader I have only to read in this look-up-table. There are some limitations (linear sampling, max input value in look up table...) but for my problem, it works very well. Cheers, Aurelien -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=57581#57581 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
The heights are relative to the ellipsoid, yes. I can send you a code link when I get back (in on the road atm). On Dec 6, 2013 7:55 AM, Sebastian Messerschmidt sebastian.messerschm...@gmx.de wrote: Am 06.12.2013 13:14, schrieb Glenn Waldron: Sabastian, osgEarth does more or less what Robert describes. We encode each vert's height value (and unit vector) in a vertex attribute. Thank you Glenn. Are the vertices relative to the Ellipsoid? Or can you point me to the sources where I can check it myself? cheers Sebastian On Dec 6, 2013 5:43 AM, Robert Osfield robert.osfi...@gmail.com wrote: Hi Sebastian, Personally I'd subclass from osgTerrain::GeometryTechnique and pass a FloatArray as a vertex attrib array or tex coord array and store the heights there. Robert. On 4 December 2013 13:23, Sebastian Messerschmidt sebastian.messerschm...@gmx.de wrote: Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing listosg-users@lists.openscenegraph.orghttp://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_ latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
HI Wojciech, Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... I've tried. Unfortunately the trigonometric functions are not available for double precision. And I'm clueless with the two-float mechanism and trigonometric functions. Also I guess I would have to push double matrices, and I don't know how OSG is handling them. I can recall Robert explaining that they are indeed pushed as double matrices, so the driver would cut them to float. cheers Sebastian Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Sebastian, Just an extra thought that came to me. Terrain LOD paging may cause the elevation jumps too. When LODs change meshes get denser or less dense - new vertices show up or excess vertices hide and that may also bring lot of surprises. Cheers, Wojtek 2013/12/5 Wojciech Lewandowski w.p.lewandow...@gmail.com Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_ latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users- openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Wojciech, Hi Sebastian, Just an extra thought that came to me. Terrain LOD paging may cause the elevation jumps too. When LODs change meshes get denser or less dense - new vertices show up or excess vertices hide and that may also bring lot of surprises. You're correct. But it seems to happen mostly if get close to the mesh, where the LOD should not change anymore. Right now I'm trying to live with it. But I'm sure there might be clever tricks to make it nicer ;-) I'm trying to get some info from the osgEarth guys, they surely know how to handle those floating point beasts. cheers Sebastian Cheers, Wojtek 2013/12/5 Wojciech Lewandowski w.p.lewandow...@gmail.com mailto:w.p.lewandow...@gmail.com Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org mailto:osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Sebastion, when you deal with geocentric data, the precision is an issue even if it is not on the GPU. The way how TerraVista (you have heard of it, right) was solving this was by introducing a matrix of doubles for the origin or the center of the tile and all the math was then with floats locally to the tile, in other words, all terrain tiles resides in their local floats coordinate system with doubles matrix on top. It is still possible to work with floats this way - that is how I did the geocentric whole earth terrapage - well, Robert pointed me to this trick. My two cents. Nick On Thu, Dec 5, 2013 at 8:28 PM, Sebastian Messerschmidt sebastian.messerschm...@gmx.de wrote: Hi Wojciech, Hi Sebastian, Just an extra thought that came to me. Terrain LOD paging may cause the elevation jumps too. When LODs change meshes get denser or less dense - new vertices show up or excess vertices hide and that may also bring lot of surprises. You're correct. But it seems to happen mostly if get close to the mesh, where the LOD should not change anymore. Right now I'm trying to live with it. But I'm sure there might be clever tricks to make it nicer ;-) I'm trying to get some info from the osgEarth guys, they surely know how to handle those floating point beasts. cheers Sebastian Cheers, Wojtek 2013/12/5 Wojciech Lewandowski w.p.lewandow...@gmail.com Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_latitude*sin_latitude); float height = p/cos(latitude) - N; return vec3(longitude, latitude, height); } void main() { vec3 ws_pos = osg_ViewMatrix * osg_ModelViewMatrix * gl_Vertex; llh = XYZ_to_llh(ws_pos); ... } fragment: void main() { if (llh.z 10.0) { gl_FragColor = vec4(1,0,0,1); } else { gl_FragColor = color; } } Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing listosg-users@lists.openscenegraph.orghttp://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org ___ osg-users mailing list
Re: [osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi Nick, Hi Sebastion, when you deal with geocentric data, the precision is an issue even if it is not on the GPU. The way how TerraVista (you have heard of it, right) was solving this was by introducing a matrix of doubles for the origin or the center of the tile and all the math was then with floats locally to the tile, in other words, all terrain tiles resides in their local floats coordinate system with doubles matrix on top. It is still possible to work with floats this way - that is how I did the geocentric whole earth terrapage - well, Robert pointed me to this trick. My two cents. I'm aware of the concept. In fact this is how OSG renders the geocentric data. The only difference is, that the View and ModelViewMatrices are pulled down to float precision in the shader. So in the ModelViewMatrix the precision is good, as it is kept double precision all the way through the graph. Essentially view and model matrices are canceling out here. My problem is the float precision product (ModelView) and ViewInverse to get the ModelMatrix again. *sigh* Maybe I can install an update callback of some kind to push a float precision ModelMatrix for the terrain tiles and use it directly. Thank you for the hint anyways :-) cheers Sebastain Nick On Thu, Dec 5, 2013 at 8:28 PM, Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de wrote: Hi Wojciech, Hi Sebastian, Just an extra thought that came to me. Terrain LOD paging may cause the elevation jumps too. When LODs change meshes get denser or less dense - new vertices show up or excess vertices hide and that may also bring lot of surprises. You're correct. But it seems to happen mostly if get close to the mesh, where the LOD should not change anymore. Right now I'm trying to live with it. But I'm sure there might be clever tricks to make it nicer ;-) I'm trying to get some info from the osgEarth guys, they surely know how to handle those floating point beasts. cheers Sebastian Cheers, Wojtek 2013/12/5 Wojciech Lewandowski w.p.lewandow...@gmail.com mailto:w.p.lewandow...@gmail.com Hi Sebastian, Perhaps your GPU can use doubles ? Many of them can these days. You may also try to refactor the code to use pairs of floats (as base and offset) which in theory may be almost as precise as double. It may be however very tricky for such non-linear math as geographic projections... Cheers, Wojtek 2013/12/5 Sebastian Messerschmidt sebastian.messerschm...@gmx.de mailto:sebastian.messerschm...@gmx.de Hi, I managed to get the local heights. After realizing, that the osg_ViewMatrix * gl_ModelViewMatrix will bring me into world space I was able to use a XYZ_to_latlonheight function in the vertex shader. There is only one catch with this: precision. It seems that the float matrices will simply cut away to much precision so I get massive flickering. Anyone any idea how to solve, improve this? In the end I simply want to draw a water surface where the height of the geometry is below a certain threshold. My problem here are the fragment where the height is almost equal to the threshold. I also tried to move the LLH-calculation to the fragment shader, but it didn't help. example: vertex-shader: out vec3 lat_lon_height; vec3 XYZ_to_llh(vec3 ws_pos) { //from osgEarth float X = xyz.x; float Y = xyz.y; float Z = xyz.z; float _radiusEquator = 6378137.0; float _radiusPolar = 6356752.3142; float flattening = (_radiusEquator-_radiusPolar)/_radiusEquator; float _eccentricitySquared = 2*flattening - flattening*flattening; float p = sqrt(X*X + Y*Y); float theta = atan(Z*_radiusEquator , (p*_radiusPolar)); float eDashSquared = (_radiusEquator*_radiusEquator - _radiusPolar*_radiusPolar)/(_radiusPolar*_radiusPolar); float sin_theta = sin(theta); float cos_theta = cos(theta); float latitude = atan( (Z + eDashSquared*_radiusPolar*sin_theta*sin_theta*sin_theta), (p - _eccentricitySquared*_radiusEquator*cos_theta*cos_theta*cos_theta) ); float longitude = atan(Y,X); float sin_latitude = sin(latitude); float N = _radiusEquator / sqrt( 1.0 - _eccentricitySquared*sin_latitude*sin_latitude); float height = p/cos(latitude) - N; return
[osg-users] [osgDEM] calculate height above ellipsoid in vertex shader
Hi, I have a osgDEM produced geocentric database. I looked into the geometryTechnique implementation to see if I somehow can access the height of the vertices above the ellipsoid in the vertex shader. Unfortunately I don't have any idea how to calculate this from the give matrices. Is this information somehow available at all? My plan is to overwrite some of the functionality in the geometryTechnique to pass the appropriate matrices to the shader. Anyone having an idea? ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org