Solved

Change MNT vertical datum from gvd28 to gvd2013 question


Badge

Hi,

 

I want to change the vertical datum of a MNT from gvd28 to gvd2013. The input and the output are projected in NAD83 MTM7.

 

I am using the following transformers

MNTThe output tif seems correct except for some artifacts (white horizontal lines) that appear in the image

MNT_2What can be the cause of those artifacts? Should I care?

 

Thanks for your input!

icon

Best answer by kailinatsafe 31 August 2022, 18:37

View original

2 replies

Userlevel 3
Badge +13

Hello @bgeorges​, thanks for posting! Are the points expected to move only vertically in this reprojection? Not horizontally?

 

If so, it's possible that the properties of the raster created by the NumericRasterizer slightly differ from the input raster, which is causing us to end up with some rows with no points, where the points fall just into the surrounding rows. Kinda like this:

 

imageThere are some potential things you can do to workaround this. One option, on the NumericRasterizer, is to change the Ground Extents option to Specify ground extents, and use the values from the input raster (as extracted by the RasterPropertyExtractor). This will only be helpful if the points are horizontally in the same place after the reprojection.

 

Alternatively, instead of a NumericRasterizer, try using a SurfaceModeller to generate a raster instead. This should work regardless of whether the points are horizontally in the same place after the reprojection.

 

On an aside, I've filed a ticket with our development team to have vertical datum support added for rasters (FMEENGINE-74575). I've also found this related idea, in the event you'd like to add some insight! Let me know if you have any related questions! Happy to help, Kailin

Badge

Thank you for this great answer. It is very instructive!

I'll try your suggestions.

 

The surfaceModeller solve my problem!

Reply