Apply geometric transformation to an EBSD map - equivalent to imwrap

[Remsdel]

Using porosity recognizable in an EBSD indexed map and a corresponding SEM-BSE (taken at 0 tilt) as reference points I can define an affine geometric transformation that link the two maps. Distortions is small but not negligible.
Knowing the transformation matrix, the geometric transformation can be applied using the imwarp MATLAB function (image toolbox) to transform the SEM-BSE map to match the EBSD map. The function interpolates the pixels from the SEM-BSE map so that the transformed map has the same resolution as the EBSD map.
Now I would like to do the opposite: transform the EBSD map to match the SEM-BSE map. However the imwarp function does not work for EBSD data. I assume that with a combination of affinetrans(ebsd, A) and interp(ebsd,xNew,yNew), there might be a possibility to replicate what imwrap is doing but it is not straighforward (to me). Would somebody have the skills to adapt imwrap for ebsd data or give the code that could replace it (using affinetrans and interp probably), it would be fantastic and a useful tool to correct distortions in an EBSD map if it can be correlated to a less-distorted image (e.g. a SEM-BSE image at zero tilt) using fiducial markers (e.g. pores).