|
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/44442.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/44442.tif ADD: texture/44442.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/44442.tif -o output/out.osgb answer: --geocentric -d dem/n30_w086_1arc_v3.tif ADD: dem/n30_w086_1arc_v3.tif -t texture/44442.tif ADD: texture/44442.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_44442.tif ERROR 1: No PROJ.4 translation for source SRS, coordinate transformation initialization has failed. Failed to reproject texture/44442.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.000138888888884,31.000138888888891) Pixel Size = (0.000277777777778,-0.000277777777778) Metadata: AREA_OR_POINT=Point DTED_CompilationDate=0002 DTED_DataEdition=02 DTED_DigitizingSystem=SRTM DTED_HorizontalAccuracy=0013 DTED_HorizontalDatum=WGS84 DTED_MaintenanceDate=0000 DTED_MaintenanceDescription=0000 DTED_MatchMergeDate=0000 DTED_MatchMergeVersion=A DTED_NimaDesignator=DTED2 DTED_OriginLatitude=0300000N DTED_OriginLongitude=0860000W 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 Left ( -86.0001389, 31.0001389) ( 86d 0' 0.50"W, 31d 0' 0.50"N) Lower Left ( -86.0001389, 29.9998611) ( 86d 0' 0.50"W, 29d59'59.50"N) Upper Right ( -84.9998611, 31.0001389) ( 84d59'59.50"W, 31d 0' 0.50"N) Lower Right ( -84.9998611, 29.9998611) ( 84d59'59.50"W, 29d59'59.50"N) Center ( -85.5000000, 30.5000000) ( 85d30' 0.00"W, 30d30' 0.00"N) Band 1 Block=3601x1 Type=Int16, ColorInterp=Gray NoData Value=-32767 Unit Type: m 3.2. gdalinfo for texture file Driver: GTiff/GeoTIFF Files: 44442.tif Size is 5000, 5000 Coordinate System is: LOCAL_CS["unnamed", UNIT["US survey foot",0.3048006096012192, AUTHORITY["EPSG","9003"]]] Origin = (1655000.330000000100000,559999.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, 559999.510) Lower Left ( 1655000.330, 554999.510) Upper Right ( 1660000.330, 559999.510) Lower Right ( 1660000.330, 554999.510) Center ( 1657500.330, 557499.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 3.3. Conclusion: We need to reproject texture file to match dem file 3.4. Reprojection Since I was pointed that t_srs for texture file is EPSG:2238 and WGS_1984 is EPSG:4326 (at least what i was able to find), then: c:\Terrain\FromUSGS\texture>gdalwarp -s_srs EPSG:2238 -t_srs EPSG:4326 44442.tif 44442_conv.tif Creating output file that is 5386P x 4662L. Processing input file 44442.tif. 0...10...20...30...40...50...60...70...80...90...100 - done. 3.5. gdalinfo for reprojected texture file Driver: GTiff/GeoTIFF Files: 44442_conv.tif Size is 5386, 4662 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 = (-85.495785570955590,30.536183332655433) Pixel Size = (0.000002974135230,-0.000002974135230) 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 ( -85.4957856, 30.5361833) ( 85d29'44.83"W, 30d32'10.26"N) Lower Left ( -85.4957856, 30.5223179) ( 85d29'44.83"W, 30d31'20.34"N) Upper Right ( -85.4797669, 30.5361833) ( 85d28'47.16"W, 30d32'10.26"N) Lower Right ( -85.4797669, 30.5223179) ( 85d28'47.16"W, 30d31'20.34"N) Center ( -85.4877762, 30.5292506) ( 85d29'15.99"W, 30d31'45.30"N) Band 1 Block=5386x1 Type=Byte, ColorInterp=Red Band 2 Block=5386x1 Type=Byte, ColorInterp=Green Band 3 Block=5386x1 Type=Byte, ColorInterp=Blue 3.6. Conclusion: Seems like both files are now in the same projection? 4. Try to use vpbmaster again with the converted texture file: vpbmaster --geocentric -t texture/44442_conv.tif -d dem/n30_w086_1arc_v3.tif -o output/out.osgb Now vpb seems work but the result is: https://drive.google.com/file/d/0ByDDImhSolf6cGFyV21PRXZuZHM/view?usp=sharing https://drive.google.com/file/d/0ByDDImhSolf6N3FvRkNJUW1nMlk/view?usp=sharing Conclusion: vpbmaster did not make texture "elevated" and too much area of elevation compared to small area of texture caused artefacts. Need to crop elevation here? If yes, what tool is the most user friendly for a newbie? 5. Try to use osgdem again - result is the same as №4. Again, what am i doing wrong? Thank you! Cheers, Elias ------------------ Read this topic online here: http://forum.openscenegraph.org/viewtopic.php?p=64048#64048 _______________________________________________ osg-users mailing list [email protected] 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 [email protected] http://lists.openscenegraph.org/listinfo.cgi/osg-users-openscenegraph.org
