demField = np.where(demData["rasterData"] == demData["header"]["nodata_value"], np.nan,
        demData["rasterData"])