Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hi! Finally, when i croped extra data from the source dem file (i used QGis), then translation works fine and i get needed terrain. Thank you! Cheers, Elias -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=64088#64088 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hy Elias, your procedure seems correct to me, it is more your conclusion that I don't understand... What do you mean by "vpb didn't make texture elevated" ? Ha you fully zoomed on the textured zone in order to be sure that it is really not elevated. When I see the size of your dem, I would say that the elevations look reduced, due to the very large size of represented zone and the fact that there are only small elevations (try to calculate the ratio of highest point of dem and the length or width of your terrain). As far as I know you can define the extend of generated terrain in osgdem (and vpb) with command line option : -e x y w h Extents of the model to generate The artefact is due to the fact that the whole current LOD containing your texture will be black (that's why when you get farer away you see more black texture) I think the result will be correct (or look correct to you, as I think that it is already correct :-) ) once you reduce the extend of your generated database. Regards, Christian Le 11/06/2015 10:52, Elias Tarasov a écrit : Hello! My statement about the problem, that had been solved, looks a little premature. The problem: Terrain is displayed incorrectly. Results: 1. https://drive.google.com/file/d/0ByDDImhSolf6cGFyV21PRXZuZHM/view?usp=sharing 2. https://drive.google.com/file/d/0ByDDImhSolf6N3FvRkNJUW1nMlk/view?usp=sharing Step-by-step actions to make a terrain 0. Gathering required data 0.1. Choose texture: https://drive.google.com/file/d/0ByDDImhSolf6S3E0M1VmR2xDd1U/view?usp=sharing 0.2. Choose dem: https://drive.google.com/file/d/0ByDDImhSolf6cDNtWjlUa3VqdXc/view?usp=sharing 0.3. Downloaded texture: https://drive.google.com/file/d/0ByDDImhSolf6eXRpYlhxMFN3NnM/view?usp=sharing 0.4. Downloaded dem: https://drive.google.com/file/d/0ByDDImhSolf6MXVLUV9kOG1xR2c/view?usp=sharing 1. Straightforward vpbmaster --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/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 2. osgdem osgdem --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/out.osgb Adding terrainTile DataSet::_run() 0 0 Now checking for plug-in osgPlugins-3.3.8/osgdb_nvtt.dll osg::Registry::addImageProcessor(ImageProcessor) Loaded plug-in osgPlugins-3.3.8/osgdb_nvtt.dll and located ImageProcessor started DataSet::createDestination(30) reprojecting to file temporaryfile_2.tif ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed. Failed to reproject texture/2.tif Time for after_reproject 0.003934 And here program crashes. 3. Manual reprojection using gdalwarp 3.1. gdalinfo for dem file: Driver: GTiff/GeoTIFF Files: n30_w086_1arc_v3.tif Size is 3601, 3601 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4326"]] Origin = (-86.0001384,31.0001391) Pixel Size = (0.0002778,-0.0002778) Metadata: AREA_OR_POINT=Point DTED_CompilationDate=0002 DTED_DataEdition=02 DTED_DigitizingSystem=SRTM DTED_HorizontalAccuracy=0013 DTED_HorizontalDatum=WGS84 DTED_MaintenanceDate= DTED_MaintenanceDescription= DTED_MatchMergeDate= DTED_MatchMergeVersion=A DTED_NimaDesignator=DTED2 DTED_OriginLatitude=030N DTED_OriginLongitude=086W DTED_Producer=USCNIMA DTED_RelHorizontalAccuracy=NA DTED_RelVerticalAccuracy=0004 DTED_SecurityCode_DSI=U DTED_SecurityCode_UHL=U DTED_UniqueRef_DSI=H24 084 DTED_UniqueRef_UHL=H24 084 DTED_VerticalAccuracy_ACC=0005 DTED_VerticalAccuracy_UHL=0005 DTED_VerticalDatum=E96 Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hello! My statement about the problem, that had been solved, looks a little premature. The problem: Terrain is displayed incorrectly. Results: 1. https://drive.google.com/file/d/0ByDDImhSolf6cGFyV21PRXZuZHM/view?usp=sharing 2. https://drive.google.com/file/d/0ByDDImhSolf6N3FvRkNJUW1nMlk/view?usp=sharing Step-by-step actions to make a terrain 0. Gathering required data 0.1. Choose texture: https://drive.google.com/file/d/0ByDDImhSolf6S3E0M1VmR2xDd1U/view?usp=sharing 0.2. Choose dem: https://drive.google.com/file/d/0ByDDImhSolf6cDNtWjlUa3VqdXc/view?usp=sharing 0.3. Downloaded texture: https://drive.google.com/file/d/0ByDDImhSolf6eXRpYlhxMFN3NnM/view?usp=sharing 0.4. Downloaded dem: https://drive.google.com/file/d/0ByDDImhSolf6MXVLUV9kOG1xR2c/view?usp=sharing 1. Straightforward vpbmaster --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/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 2. osgdem osgdem --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/out.osgb Adding terrainTile DataSet::_run() 0 0 Now checking for plug-in osgPlugins-3.3.8/osgdb_nvtt.dll osg::Registry::addImageProcessor(ImageProcessor) Loaded plug-in osgPlugins-3.3.8/osgdb_nvtt.dll and located ImageProcessor started DataSet::createDestination(30) reprojecting to file temporaryfile_2.tif ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed. Failed to reproject texture/2.tif Time for after_reproject 0.003934 And here program crashes. 3. Manual reprojection using gdalwarp 3.1. gdalinfo for dem file: Driver: GTiff/GeoTIFF Files: n30_w086_1arc_v3.tif Size is 3601, 3601 Coordinate System is: GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]] Origin = (-86.0001384,31.0001391) Pixel Size = (0.0002778,-0.0002778) Metadata: AREA_OR_POINT=Point DTED_CompilationDate=0002 DTED_DataEdition=02 DTED_DigitizingSystem=SRTM DTED_HorizontalAccuracy=0013 DTED_HorizontalDatum=WGS84 DTED_MaintenanceDate= DTED_MaintenanceDescription= DTED_MatchMergeDate= DTED_MatchMergeVersion=A DTED_NimaDesignator=DTED2 DTED_OriginLatitude=030N DTED_OriginLongitude=086W DTED_Producer=USCNIMA DTED_RelHorizontalAccuracy=NA DTED_RelVerticalAccuracy=0004 DTED_SecurityCode_DSI=U DTED_SecurityCode_UHL=U DTED_UniqueRef_DSI=H24 084 DTED_UniqueRef_UHL=H24 084
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hello! My statement about the problem, that had been solved, looks a little premature. The problem: Terrain is displayed incorrectly. Results: 1. https://drive.google.com/file/d/0ByDDImhSolf6cGFyV21PRXZuZHM/view?usp=sharing 2. https://drive.google.com/file/d/0ByDDImhSolf6N3FvRkNJUW1nMlk/view?usp=sharing Step-by-step actions to make a terrain 0. Gathering required data 0.1. Choose texture: https://drive.google.com/file/d/0ByDDImhSolf6S3E0M1VmR2xDd1U/view?usp=sharing 0.2. Choose dem: https://drive.google.com/file/d/0ByDDImhSolf6cDNtWjlUa3VqdXc/view?usp=sharing 0.3. Downloaded texture: https://drive.google.com/file/d/0ByDDImhSolf6eXRpYlhxMFN3NnM/view?usp=sharing 0.4. Downloaded dem: https://drive.google.com/file/d/0ByDDImhSolf6MXVLUV9kOG1xR2c/view?usp=sharing 1. Straightforward vpbmaster --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/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 2. osgdem osgdem --geocentric -d dem/n30_w086_1arc_v3.tif -t texture/2.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/2.tif ADD: texture/2.tif -o output/out.osgb Adding terrainTile DataSet::_run() 0 0 Now checking for plug-in osgPlugins-3.3.8/osgdb_nvtt.dll osg::Registry::addImageProcessor(ImageProcessor) Loaded plug-in osgPlugins-3.3.8/osgdb_nvtt.dll and located ImageProcessor started DataSet::createDestination(30) reprojecting to file temporaryfile_2.tif ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed. Failed to reproject texture/2.tif Time for after_reproject 0.003934 And here program crashes. 3. Manual reprojection using gdalwarp 3.1. gdalinfo for dem file: Driver: GTiff/GeoTIFF Files: n30_w086_1arc_v3.tif Size is 3601, 3601 Coordinate System is: GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]] Origin = (-86.0001384,31.0001391) Pixel Size = (0.0002778,-0.0002778) Metadata: AREA_OR_POINT=Point DTED_CompilationDate=0002 DTED_DataEdition=02 DTED_DigitizingSystem=SRTM DTED_HorizontalAccuracy=0013 DTED_HorizontalDatum=WGS84 DTED_MaintenanceDate= DTED_MaintenanceDescription= DTED_MatchMergeDate= DTED_MatchMergeVersion=A DTED_NimaDesignator=DTED2 DTED_OriginLatitude=030N DTED_OriginLongitude=086W DTED_Producer=USCNIMA DTED_RelHorizontalAccuracy=NA DTED_RelVerticalAccuracy=0004 DTED_SecurityCode_DSI=U DTED_SecurityCode_UHL=U DTED_UniqueRef_DSI=H24 084 DTED_UniqueRef_UHL=H24 084
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hello, how I found it ? Easy (LOL!), the EPSG:2238 is the official map projection for north of Florida and from the values in the file (Center ( 1657500.330, 572499.510)) you can see (okay, I admit if you are used to look at this kind of data :-) ) that it is a distance in feet from a reference point. That is typical for a Mercantor or Lambert projection. In every case I would advise you to check in a converter what latitude and longitude you will obtain from the centre of the image and double check with google maps or equivalent that you find again the same point. Because you can convert a lot of different projections to latitude and longitude, as the equations are nearly the same only the correction factors or the reference ellipsoid will change. So it will convert but you are not sure to be in the right place in the world... For your information, the EPSG:3857 is a non officially recognized projection as it uses a sphere as earth reference and not a ellipsoid. Furthermore this projection system is centred on Equator/Greenwich so you may have in the United States an error up to 43 km in map projection, (nearly 21 km on the ground). Regards, Christian Le 09/06/2015 23:15, Elias Tarasov a écrit : Hello! Christian Schulte wrote: Hello! I totally agree with Sebastian. These data is surely geo-referenced, as I myself often used in the past data from USGS. Locking at the metadata you provided I would say that you are in an EPSG:2238 Lambert Conformal Conic projection, which is the main projection system for north Florida. You may try to specify to GDAL the source projection using this EPSG code. Regards, Christian Yes, it worked! However, i dont understand how did you figure out that projection is 2238. No such a number is presented into metadata. For a test, i also used EPGS:3857 and it seems to work similar. Thanks a lot! Elias. Post generated by Mail2Forum -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=64020#64020 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org -- SCHULTE Christian Ingénieur Recherche Responsable du Laboratoire de Simulation ONERA - DCSD/PSEV Département Commande des Systèmes et Dynamique du Vol ONERA - Centre de Salon de Provence BA 701 13661 SALON AIR Cedex - France Tel :04.90.17.01.45 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Am 09.06.2015 um 23:05 schrieb Elias Tarasov: Hello! Thanks for quick reply, it helped to solve a problem. SMesserschmidt wrote: Am 09.06.2015 um 11:35 schrieb Elias Tarasov: ... there is a bunch of links to get data. Since --geocentric option allows not to use elevation data, then only textures are needed. Who told you this? Of course you can use elevation data in geocentric mode ... Simply use -t for imagery and -d for digitial elevation data ... Don't doubt the arguments we used this time. Yes, there is no -d option and thus no elevation map specified. As we have already indicated to use the geocentric system to build from the source, VPB will automatically use flat sea-level elevation data and construct the earth geometry according to the given GeoTiff imagery. Yes but you still CAN use elevation data relative to the ellipsoid. And here is a problem: i can't get georeferenced textures from USGS. They are referenced, but in a non-supported coordinate frame. I never had problems of this kind yet. But there are definitively geo-referenced data-sets on USGS, as I'm using them myself. Yes, they are. Just in a different projection, it seems. Use LandSAT imagery, which is WGS84 projected, so there should no problems here. Actually, problems are exist. This is the LandSat image from USGS ... Now datum is WGS84, but vpbmaster shows again Error: vpbmaster can not run without all source data being in the correct destination coordinates system, please reproject them. Maybe because PROJCS[WGS 84 / UTM zone 33N, not just WGS84? That is UTM with reference datum WGS84*. Don't mix them up. One is the projection (UTM at the zone 33North) using the reference ellipsoid WGS84. Usually it will be transformed/reprojected automatically. Can you try osgdem instead of vpbmaster? I think the distributed version (vpbmaster) indeed won't work on data not being adapted. Also make sure you have the gdal/proj libraries accessible for the osgdem application. *(WGS84 defines the sea level by approximating the earth's shape by an ellipsoid) If you need some sample data I can share a small set on googledrive etc. That's very kind, Sebastian, but the problem seems to have been solved already. To use data with vpb i need to reproject them. gdalwarp -s_srs EPSG:2238 -t_srs EPSG:4326 src.tif dst.tif Now produced *.osgb image looks correctly into OSG viewer. However, it also looks quite plain. So elevation data is needed? If you need elevation data or not is totally up to you. Usually you don't need to re-project manually (see above). Cheers Sebastian ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
[osg-users] [vpb] Correct way to get texture data from USGS or any other source?
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\vol001gdalinfo 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.3300010,574999.511) Pixel Size = (1.000,-1.000) 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, 56.510) Upper Right ( 166.330, 574999.510) Lower Right ( 166.330, 56.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\vol001gdalwarp -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\vol001gdalinfo 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.3300010,574999.511) Pixel Size = (1.000,-1.000) 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, 56.510) (Invalid angle,Invalid angle) Upper Right ( 166.330, 574999.510) (Invalid angle,Invalid angle) Lower Right ( 166.330, 56.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:\Terrainvpbmaster --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 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Am 09.06.2015 um 11:35 schrieb 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. Who told you this? Of course you can use elevation data in geocentric mode ... Simply use -t for imagery and -d for digitial elevation data And here is a problem: i can't get georeferenced textures from USGS. They are referenced, but in a non-supported coordinate frame. I never had problems of this kind yet. But there are definitively geo-referenced data-sets on USGS, as I'm using them myself. Use LandSAT imagery, which is WGS84 projected, so there should no problems here. If you need some sample data I can share a small set on googledrive etc. All files i choose from different sets there are without georeferenced information. gdalinfo output: ... Yes it seems there is no valid projection. But still the data is referenced. What am i doing wrong? Without having your data, this is hard to tell Thanks a lot and best regards! Elias -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=63998#63998 ___ 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] [vpb] Correct way to get texture data from USGS or any other source?
Hello! I totally agree with Sebastian. These data is surely geo-referenced, as I myself often used in the past data from USGS. Locking at the metadata you provided I would say that you are in an EPSG:2238 Lambert Conformal Conic projection, which is the main projection system for north Florida. You may try to specify to GDAL the source projection using this EPSG code. Regards, Christian Le 09/06/2015 12:59, Sebastian Messerschmidt a écrit : Am 09.06.2015 um 11:35 schrieb 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. Who told you this? Of course you can use elevation data in geocentric mode ... Simply use -t for imagery and -d for digitial elevation data And here is a problem: i can't get georeferenced textures from USGS. They are referenced, but in a non-supported coordinate frame. I never had problems of this kind yet. But there are definitively geo-referenced data-sets on USGS, as I'm using them myself. Use LandSAT imagery, which is WGS84 projected, so there should no problems here. If you need some sample data I can share a small set on googledrive etc. All files i choose from different sets there are without georeferenced information. gdalinfo output: ... Yes it seems there is no valid projection. But still the data is referenced. What am i doing wrong? Without having your data, this is hard to tell Thanks a lot and best regards! Elias -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=63998#63998 ___ 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 -- SCHULTE Christian Ingénieur Recherche Responsable du Laboratoire de Simulation ONERA - DCSD/PSEV Département Commande des Systèmes et Dynamique du Vol ONERA - Centre de Salon de Provence BA 701 13661 SALON AIR Cedex - France Tel :04.90.17.01.45 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hello! Thanks for quick reply, it helped to solve a problem. SMesserschmidt wrote: Am 09.06.2015 um 11:35 schrieb 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. Who told you this? Of course you can use elevation data in geocentric mode ... Simply use -t for imagery and -d for digitial elevation data The book is called OpenSceneGraph 3 Cookbook Authors: Rui Wang Xuelei Qian Citation: To build databases in geographic coordinate, we can simply use the --geocentric option while executing vbpmaster. The complete command is: # vpbmaster -t data/TrueMarble.4km.10800x5400.tif --geocentric -o output/out.osgb Don't doubt the arguments we used this time. Yes, there is no -d option and thus no elevation map specified. As we have already indicated to use the geocentric system to build from the source, VPB will automatically use flat sea-level elevation data and construct the earth geometry according to the given GeoTiff imagery. And here is a problem: i can't get georeferenced textures from USGS. They are referenced, but in a non-supported coordinate frame. I never had problems of this kind yet. But there are definitively geo-referenced data-sets on USGS, as I'm using them myself. Yes, they are. Just in a different projection, it seems. Use LandSAT imagery, which is WGS84 projected, so there should no problems here. Actually, problems are exist. This is the LandSat image from USGS c:\Terrain\EPP191R025_7F19990915gdalinfo test_texture.tif Driver: GTiff/GeoTIFF Files: test_texture.tif Size is 17884, 16140 Coordinate System is: PROJCS[WGS 84 / UTM zone 33N, GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]], PROJECTION[Transverse_Mercator], PARAMETER[latitude_of_origin,0], PARAMETER[central_meridian,15], PARAMETER[scale_factor,0.9996], PARAMETER[false_easting,50], PARAMETER[false_northing,0], UNIT[metre,1, AUTHORITY[EPSG,9001]], AUTHORITY[EPSG,32633]] Origin = (391518.750,5685165.750) Pixel Size = (14.250,-14.250) Metadata: AREA_OR_POINT=Point TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch) TIFFTAG_XRESOLUTION=72 TIFFTAG_YRESOLUTION=72 Image Structure Metadata: INTERLEAVE=PIXEL Corner Coordinates: Upper Left ( 391518.750, 5685165.750) ( 13d26'37.19E, 51d18'26.83N) Lower Left ( 391518.750, 5455170.750) ( 13d30'34.78E, 49d14'22.90N) Upper Right ( 646365.750, 5685165.750) ( 17d 5'58.23E, 51d17'56.30N) Lower Right ( 646365.750, 5455170.750) ( 17d 0'37.84E, 49d13'54.52N) Center ( 518942.250, 5570168.250) ( 15d15'57.15E, 50d17' 0.02N) Band 1 Block=17884x1 Type=Byte, ColorInterp=Red Band 2 Block=17884x1 Type=Byte, ColorInterp=Green Band 3 Block=17884x1 Type=Byte, ColorInterp=Blue Now datum is WGS84, but vpbmaster shows again Error: vpbmaster can not run without all source data being in the correct destination coordinates system, please reproject them. Maybe because PROJCS[WGS 84 / UTM zone 33N, not just WGS84? If you need some sample data I can share a small set on googledrive etc. That's very kind, Sebastian, but the problem seems to have been solved already. To use data with vpb i need to reproject them. gdalwarp -s_srs EPSG:2238 -t_srs EPSG:4326 src.tif dst.tif Now produced *.osgb image looks correctly into OSG viewer. However, it also looks quite plain. So elevation data is needed? Thanks a lot! Elias -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=63998#63998 ___ Post generated by Mail2Forum -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=64019#64019 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
Re: [osg-users] [vpb] Correct way to get texture data from USGS or any other source?
Hello! Christian Schulte wrote: Hello! I totally agree with Sebastian. These data is surely geo-referenced, as I myself often used in the past data from USGS. Locking at the metadata you provided I would say that you are in an EPSG:2238 Lambert Conformal Conic projection, which is the main projection system for north Florida. You may try to specify to GDAL the source projection using this EPSG code. Regards, Christian Yes, it worked! However, i dont understand how did you figure out that projection is 2238. No such a number is presented into metadata. For a test, i also used EPGS:3857 and it seems to work similar. Thanks a lot! Elias. Post generated by Mail2Forum -- Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=64020#64020 ___ osg-users mailing list osg-users@lists.openscenegraph.org http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
