OpenSceneGraph Forum Forum Index OpenSceneGraph Forum
Official forum which mirrors the existent OSG mailing lists. Messages posted here are forwarded to the mailing list and vice versa.
 
   FAQFAQ    SearchSearch    MemberlistMemberlist    RulesRules    UsergroupsUsergroups    RegisterRegister 
 Mail2Forum SettingsMail2Forum Settings  ProfileProfile   Log in to check your private messagesLog in to check your private messages   Log inLog in 
   AlbumAlbum  OpenSceneGraph IRC ChatOpenSceneGraph IRC Chat   SmartFeedSmartFeed 

Correct way to get texture data from USGS or any other source?


 
Post new topic   Reply to topic    OpenSceneGraph Forum Forum Index -> VirtualPlanetBuilder [vpb]
View previous topic :: View next topic  
Author Message
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Tue Jun 09, 2015 9:35 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

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\vol001>gdalinfo 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.330000000100000,574999.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, 574999.510)
Lower Left ( 1655000.330, 569999.510)
Upper Right ( 1660000.330, 574999.510)
Lower Right ( 1660000.330, 569999.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\vol001>gdalwarp -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\vol001>gdalinfo 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.330000000100000,574999.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, 574999.510) (Invalid angle,Invalid angle)
Lower Left ( 1655000.330, 569999.510) (Invalid angle,Invalid angle)
Upper Right ( 1660000.330, 574999.510) (Invalid angle,Invalid angle)
Lower Right ( 1660000.330, 569999.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:\Terrain>vpbmaster --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
Back to top
View user's profile Send private message
SMesserschmidt (Sebastian Messerschmidt)
Forum Moderator


Joined: 10 Sep 2013
Posts: 764

PostPosted: Tue Jun 09, 2015 10:59 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

Am 09.06.2015 um 11:35 schrieb Elias Tarasov:
Quote:
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

Quote:
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.


Quote:
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.
Quote:



What am i doing wrong?
Without having your data, this is hard to tell
Quote:
Thanks a lot and best regards!

Elias

------------------
Read this topic online here:
http://forum.openscenegraph.org/viewtopic.php?p=63998#63998









------------------
Post generated by Mail2Forum
Back to top
View user's profile Send private message
Christian Schulte
Guest





PostPosted: Tue Jun 09, 2015 11:32 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

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 :
Quote:
Am 09.06.2015 um 11:35 schrieb Elias Tarasov:
Quote:
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

Quote:
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.


Quote:
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.
Quote:



What am i doing wrong?
Without having your data, this is hard to tell
Quote:
Thanks a lot and best regards!

Elias

------------------
Read this topic online here:
http://forum.openscenegraph.org/viewtopic.php?p=63998#63998













--
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



------------------
Post generated by Mail2Forum
Back to top
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Tue Jun 09, 2015 9:05 pm    Post subject:
Re: Correct way to get texture data from USGS or any other source?
Reply with quote

Hello!
Thanks for quick reply, it helped to solve a problem.
SMesserschmidt wrote:
Am 09.06.2015 um 11:35 schrieb Elias Tarasov:
Quote:
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.

Quote:

Quote:

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.
Quote:

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_7F19990915>gdalinfo 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",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32633"]]
Origin = (391518.750000000000000,5685165.750000000000000)
Pixel Size = (14.250000000000000,-14.250000000000000)
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.19"E, 51d18'26.83"N)
Lower Left ( 391518.750, 5455170.750) ( 13d30'34.78"E, 49d14'22.90"N)
Upper Right ( 646365.750, 5685165.750) ( 17d 5'58.23"E, 51d17'56.30"N)
Lower Right ( 646365.750, 5455170.750) ( 17d 0'37.84"E, 49d13'54.52"N)
Center ( 518942.250, 5570168.250) ( 15d15'57.15"E, 50d17' 0.02"N)
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?

Quote:

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
Back to top
View user's profile Send private message
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Tue Jun 09, 2015 9:15 pm    Post subject:
Re: Correct way to get texture data from USGS or any other source?
Reply with quote

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
Back to top
View user's profile Send private message
SMesserschmidt (Sebastian Messerschmidt)
Forum Moderator


Joined: 10 Sep 2013
Posts: 764

PostPosted: Wed Jun 10, 2015 6:27 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

Am 09.06.2015 um 23:05 schrieb Elias Tarasov:
Quote:
Hello!
Thanks for quick reply, it helped to solve a problem.

SMesserschmidt wrote:
Quote:
Am 09.06.2015 um 11:35 schrieb Elias Tarasov:

Quote:
...
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.
Quote:

Quote:

Quote:
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.

Quote:
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)
Quote:


Quote:
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


------------------
Post generated by Mail2Forum
Back to top
View user's profile Send private message
Christian Schulte
Guest





PostPosted: Wed Jun 10, 2015 7:26 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

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 Smile ) 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 :
Quote:
Hello!


Christian Schulte wrote:
Quote:
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










--
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



------------------
Post generated by Mail2Forum
Back to top
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Thu Jun 11, 2015 8:50 am    Post subject:
Reply with quote

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


Last edited by Elias Tarasov on Thu Jun 11, 2015 9:49 am; edited 1 time in total
Back to top
View user's profile Send private message
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Thu Jun 11, 2015 8:52 am    Post subject:
Reply with quote

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
Back to top
View user's profile Send private message
Christian Schulte
Guest





PostPosted: Thu Jun 11, 2015 11:00 am    Post subject:
Correct way to get texture data from USGS or any other source?
Reply with quote

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 Smile ) once you reduce the extend of your generated database.

Regards,

Christian



Le 11/06/2015 10:52, Elias Tarasov a écrit :

Quote:
Quote:
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:TerrainFromUSGStexture>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
(
Only registered users can see emails on this board!
Get registred or enter the forums!
)
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


------------------
Post generated by Mail2Forum
Back to top
Elias Tarasov
User


Joined: 20 Apr 2015
Posts: 36

PostPosted: Sun Jun 14, 2015 5:34 pm    Post subject:
Reply with quote

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
Back to top
View user's profile Send private message
Display posts from previous:   
Post new topic   Reply to topic    OpenSceneGraph Forum Forum Index -> VirtualPlanetBuilder [vpb] All times are GMT
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum
You cannot attach files in this forum
You cannot download files in this forum

Similar Topics
Topic Author Forum Replies Posted
No new posts Render to texture with GL3. Raizel General 3 Tue May 23, 2017 10:33 am View latest post
No new posts How to replace the texture data in th... xdfanfan General 1 Wed Apr 05, 2017 12:44 pm View latest post
No new posts Multitexture: projection a texture on... Ekaterina General 2 Mon Feb 06, 2017 1:48 pm View latest post
No new posts How to set texture coordinates Ekaterina General 0 Fri Jan 13, 2017 5:39 pm View latest post
No new posts Camera::resizeAttachments + texture rect James Turner General 1 Wed Jan 04, 2017 5:22 pm View latest post


Board Security Anti Bot Question MOD - phpBB MOD against Spam Bots
Powered by phpBB © 2001, 2005 phpBB Group
Protected by Anti-Spam ACP