Find spiral 2D matrix (List of Lists) of ints in a clockwise spiral
Just reviewing solution 1.


The solution is clever (repeatedly extract the first row and rotate the remainder by a quarter-turn) and the code terse.
However, the code destructively modifies the input matrix by popping the first row. This would be surprising to the caller and might make it more difficult to use or to test.
The code does not work in Python 3: it fails with the exception:

TypeError: 'zip' object is not subscriptable


In Python 3 you would have to write list(zip(...)). Since Python 2.7 is close to its end of life it is a good idea to write code that would be easy to port to Python 3. In this case you could use

from future_builtins import zip


to get a version of the built-in zip function that's compatible with Python 3.
If the input is an \$n×n\$ matrix, then rotating the remainder takes \$O(n^2)\$ and this has to be done \$O(n)\$ times, leading to a total runtime of \$O(n^3)\$. But intuitively it ought to be possible to output the spiral in \$O(n^2)\$.

We can achieve this by implementing the same approach (repeatedly extract the first row and rotate the remainder by a quarter-turn) but instead of modifying the matrix, we instead modify our system of coordinates. This is probably easiest to see if I give the code first:

def find_spiral_2(matrix):
    "Return list of elements of matrix in clockwise spiral order."
    h, w = len(matrix), len(matrix[0]) # Height and width of remainder
    i, j = 0, -1                # Current position in spiral
    di, dj = 0, 1               # Current direction of movement
    spiral = []
    while h &gt; 0:
        # Output top row
        for _ in range(w):
            i += di
            j += dj
            spiral.append(matrix[i][j])
        # Remove top row and rotate quarter-turn
        h, w = w, h - 1
        di, dj = dj, -di
    return spiral


The way this works is that we maintain the state of the algorithm in the form of the height \$h\$ and width \$w\$ of the remainder of the matrix, the coordinates \$(i, j)\$ of our current position along the spiral, and a vector \$(di, dj)\$ giving the current direction of movement along the spiral path.

To output the top row, we follow the current direction for \$w\$ steps, leaving us in the top-right corner of (the remainder of) the matrix. Then, as in your solution, we remove the top row and rotate the remainder though a quarter-turn—or rather, we modify our coordinates to have this effect. Removing the top row subtracts 1 from \$h\$. Rotating the remainder through a quarter-turn causes \$w\$ and \$h\$ to swap (height becomes width and vice versa) and the current direction of movement to rotate through a quarter-turn. You should check that the line di, dj = dj, -di achieves this.

It should be clear that this does not modify matrix, and runs in \$O(n^2)\$ as required.