Aberration retrieval for high-NA optical systems using the Extended Nijboer-Zernike theory