I am trying to do some change detection between SFM derived DEM products, but the DEM i would like to use as my base elevation (and one other) is in NAD83, the remainder are WGS84. I would like to reproject the raster values (elevations) to WGS84 datum. Arcmap 10.4 cannot do vertical datum shifts, and VDatum does not seem to be changing the values within the raster when asked to do a datum shift, because it does not 'know' that the values are elevation. Can I use GDAL or rasterio to tackle this problem? I am new to geospatial analysis in python, but I am familiar with the language .