Vectorized Diagonal of any Multi-Dimensional Matrix

Suppose you have a matrix which was generated by the following MATLAB command:
>> ND = repmat(diag(1:3),[1 1 2 2 1 2]);
The size of this matrix is found to be: [3 3 2 2 1 2]

This routine can return the diagonal terms of this matrix simply by specifying the row dimension and the column dimension of the 2-D plane of interest:
>> dND = diagND(ND,1,2);

The resulting Matrix, dND, has the following size [1 3 2 2 1 2] and the expected solution of:
[1 2 3] across dimensions 3 through 6.