Skip to main content
  • EN

Hye, 

I'm trying to create a 3D terrain model (export format : .obj/.fbx/.3ds or any other textured mesh format) from imagery raster (France, IGN) and elevation raster (France, IGN).

Inputs formats are : 

1. For imagery

Multiple *.jp2 files with *.tab files associated (same file name) 

A sample of .tab file : 

!table
!version 300
!charset WindowsLatin1
Definition Table
File "01-2015-0830-6540-LA93-5M00-E100.jp2"
Type "RASTER"
(830000.00,6540000.00) (0,0) Label "Pt 1",
(840000.00,6540000.00) (2000,0) Label "Pt 2",
(840000.00,6530000.00) (2000,2000) Label "Pt 3",
(830000.00,6530000.00) (0,2000) Label "Pt 4"
CoordSys Earth Projection 2003, 33, 7, 3, 46.5, 44, 49.00000000001, 700000, 6600000, -792421, 5278231, 3520778, 9741029
Units "m"

2. For elevation

Multiple *.asc files with *.asc.aux.xml files associated (same file name) 

An extract of *.asc file :

ncols        1000
nrows        1000
xllcorner    74962.500000000000
yllcorner    6750037.500000000000
cellsize     75.000000000000
NODATA_value  -99999.00
 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -9999....]

An extract of *.asc.aux.xml file : 

<PAMDataset>
  <SRS>COMPD_CS0"RGF93 / Lambert-93 + NGF-IGN69 height",PROJCSe"RGF93 / Lambert-93",GEOGCS*"RGF93",DATUM "Reseau_Geodesique_Francais_1993",SPHEROID<"GRS 1980",6378137,298.257222101,AUTHORITYm"EPSG","7019"]],TOWGS84h0,0,0,0,0,0,0],AUTHORITY "EPSG","6171"]],PRIMEMF"Greenwich",0,AUTHORITYs"EPSG","8901"]],UNITS"degree",0.0174532925199433,AUTHORITY2"EPSG","9122"]],AUTHORITY9"EPSG","4171"]],PROJECTION,"Lambert_Conformal_Conic_2SP"],PARAMETERi"standard_parallel_1",49],PARAMETERT"standard_parallel_2",44],PARAMETERI"latitude_of_origin",46.5],PARAMETER""central_meridian",3],PARAMETERf"false_easting",700000],PARAMETERd"false_northing",6600000],UNIT""metre",1,AUTHORITY""EPSG","9001"]],AXISt"X",EAST],AXIS,"Y",NORTH],AUTHORITYt"EPSG","2154"]],VERT_CST"NGF-IGN69 height",VERT_DATUMR"Nivellement General de la France - IGN69",2005,AUTHORITY["EPSG","5119"]],UNIT["metre",1,AUTHORITYN"EPSG","9001"]],AXISG"Up",UP],AUTHORITY["EPSG","5720"]],AUTHORITYD"EPSG","5698"]]</SRS>
  <PAMRasterBand band="1">
    <Metadata>
      <MDI key="STATISTICS_MAXIMUM">102.4</MDI>
      <MDI key="STATISTICS_MEAN">43.77518150174</MDI>
      <MDI key="STATISTICS_MINIMUM">-2.79</MDI>
      <MDI key="STATISTICS_STDDEV">26.137841691288</MDI>
    </Metadata>
  </PAMRasterBand>
</PAMDataset>

I'm discovering FME and I'm working with a demo licence, so I need a complete explaination of the workflow please... 

And... I'm French, so... my english is very bad :/ Sorry !

,

Hye, 

I'm trying to create a 3D terrain model (export format : .obj/.fbx/.3ds or any other textured mesh format) from imagery raster (France, IGN) and elevation raster (France, IGN).

Inputs formats are : 

1. For imagery

Multiple *.jp2 files with *.tab files associated (same file name) 

A sample of .tab file : 

!table
!version 300
!charset WindowsLatin1


Definition Table
File "01-2015-0830-6540-LA93-5M00-E100.jp2"
Type "RASTER"
(830000.00,6540000.00) (0,0) Label "Pt 1",
(840000.00,6540000.00) (2000,0) Label "Pt 2",
(840000.00,6530000.00) (2000,2000) Label "Pt 3",
(830000.00,6530000.00) (0,2000) Label "Pt 4"
CoordSys Earth Projection 2003, 33, 7, 3, 46.5, 44, 49.00000000001, 700000, 6600000, -792421, 5278231, 3520778, 9741029
Units "m"

2. For elevation

Multiple *.asc files with *.asc.aux.xml files associated (same file name) 

An extract of *.asc file :

ncols        1000
nrows        1000
xllcorner    74962.500000000000
yllcorner    6750037.500000000000
cellsize     75.000000000000
NODATA_value  -99999.00
 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -99999.00 -9999Â...]

An extract of *.asc.aux.xml file : 

<PAMDataset>
  <SRS>COMPD_CS0"RGF93 / Lambert-93 + NGF-IGN69 height",PROJCS0"RGF93 / Lambert-93",GEOGCS9"RGF93",DATUM9"Reseau_Geodesique_Francais_1993",SPHEROID9"GRS 1980",6378137,298.257222101,AUTHORITY-"EPSG","7019"]],TOWGS8490,0,0,0,0,0,0],AUTHORITY9"EPSG","6171"]],PRIMEM."Greenwich",0,AUTHORITY0"EPSG","8901"]],UNITc"degree",0.0174532925199433,AUTHORITYu"EPSG","9122"]],AUTHORITY<"EPSG","4171"]],PROJECTION "Lambert_Conformal_Conic_2SP"],PARAMETERb"standard_parallel_1",49],PARAMETERS"standard_parallel_2",44],PARAMETER9"latitude_of_origin",46.5],PARAMETER_"central_meridian",3],PARAMETER1"false_easting",700000],PARAMETER,"false_northing",6600000],UNIT]"metre",1,AUTHORITY6"EPSG","9001"]],AXISw"X",EAST],AXIST"Y",NORTH],AUTHORITYI"EPSG","2154"]],VERT_CS1"NGF-IGN69 height",VERT_DATUM]"Nivellement General de la France - IGN69",2005,AUTHORITY_"EPSG","5119"]],UNITR"metre",1,AUTHORITY_"EPSG","9001"]],AXISa"Up",UP],AUTHORITY4"EPSG","5720"]],AUTHORITY_"EPSG","5698"]]</SRS>
  <PAMRasterBand band="1">
    <Metadata>
      <MDI key="STATISTICS_MAXIMUM">102.4</MDI>
      <MDI key="STATISTICS_MEAN">43.77518150174</MDI>
      <MDI key="STATISTICS_MINIMUM">-2.79</MDI>
      <MDI key="STATISTICS_STDDEV">26.137841691288</MDI>
    </Metadata>
  </PAMRasterBand>
</PAMDataset>

I'm discovering FME and I'm working with a demo licence, so I need a complete explaination of the workflow please... 

And... I'm French, so... my english is very bad :/ Sorry !

Hi @loicroybonjob,

The basic process would be to use CsmapReprojectors to reproject the data into a common coordinate system (if they are not already in one), then use a TINGenerator to convert the elevation raster to a surface. Next, use an AppearanceSetter to texture the surface, with Texture Coordinate Generation Parameters -> Texture Mapping Type set to Top Down Georeferenced.

The resulting textured surface can then be written to your desired 3d output format.

texturedem.fmw

Reply


Community Stats

31,789
Posts
120,751
Replies
39,405
Members

Latest FME

Terms & ConditionsCookie settings
A product of
Safe Software respectfully acknowledges that we live, learn and work on the traditional and unceded territories of the Kwantlen, Katzie, and Semiahmoo First Nations.