For pixel to lat/lon, the input vector is, [x, y, 0, 1] transformed. Post multiplying the transform matrix by the vector gives (x, y) coordinates for the lambert transform.
For the reverse, do the inverse lambert transform to get the (x, y) coordinates. Then you you have to solve a 2x2 system of linear equations to get the pixel locations.