Hills and valleys hidden under the ice of Thwaites Glacier have an impact on ice flow and future ice loss, but there are not many three-dimensional observations of their location or size. We apply a mathematical theory to new high-resolution observations of the ice surface to predict the bed topography beneath the ice. There is a good correlation with ice-penetrating radar observations. The method may be useful in areas with few direct observations or as a further constraint for other methods.