Compare height values between the surface model and Height model download with FME

The example describes how to load a Surface Model File (LAZ) and Elevation Data from Elevation Model Download and compare the elevation values. The difference is written to a GeoTIFF.

Use fmw files to solve the task

To solve the task, there are two fmw files you can use. The first file brings up the bounding box for the area to be compared and then calls the second file with max / min coordinates as input.

This is what the contents of file number one look like:

File number one: CompareforDMFranYtmodellMedWCSFranHojdmodellDownloadWS1.fmw

Excerpt from FMW file number one
Contents of FMW file number one. Click on the image to see it larger.

This is what the contents of file number two look like:

File number two: CompareforDMFranYtmodelMedWCSFranHojdmodellDownloadWS2.fmw

Excerpt from FMW file number two
Contents of FMW file number two. Click on the image to see it larger.

The reason why there are two files is because we did not come up with a better way to extract WCS heights from Elevation model download with boundingbox coordinates from the surface model file to be checked.

It should be able to fix with for example a FeatureReader but it did not work when we tested it.

Description of the contents of the files

In file number one reads a surface model file (laz) and from it its bounding box is removed.

The workspace parameters in the WorkspaceRunner transformer in file number one should be those needed to run file number two, see example in the image below:

Excerpts from settings in file number one
Parameters in file number one needed to run file number two. Click on the image to see it larger.

In file number two, a raster DEM is created from the surface model file pss that has been described in document Convert Surface Model FileMedFME.docx. Also read altitude data from Altitude data download pss that you can see examples in the Download height data guide with FME .

The height grid from the WCS is then sampled to get the same resolution as the DEM from the surface model. Then the NoData value is set to - 999 for both data streams, after which the difference between the height data values is calculated by taking the difference between their pixel values.

Example of output

The image below shows examples of output. In the difference picture you can see forest areas and buildings because they are very different from the ground level model from Höjddata download:

Example of output with dark and light areas on a map
Example of output. Click on the image to see it larger.

Contents of this page may be automatically translated, we take no responsibility for the accuracy of the translation. Feel free to contact our customer support centre if you have any questions.

Read more about our website