After doing auto bed leveling, keep the plane equation coefficients. You can use M500/M501 to store/load plane coefficients, and then use M504 to apply it to the rotation matrix.
Following is the G code sample:

// Doing auto bed leveling, and then store the results.
G28
G29
M500

// Home all axis, then apply plane equation coefficients, which was calculated by auto leveling before.
G28
M504