def METHOD_NAME(self, lat, lon) -> Tuple[float, float, float]:
    xs, ys, zs = rasterio.warp.transform("EPSG:4326", self.crs, [lon], [lat], [0])
    x = xs[0] * self.scale - self.offset[0]
    y = ys[0] * self.scale - self.offset[1]
    z = zs[0] * self.scale - self.offset[2]
    y, z = z, -y
    return x, y, z