You can construct it as an array of tuples and use the array builder syntax to see if that works faster. Something like:
im_arr = [(1 - (r[x] / 10813),)*im_elev[x] for x in r_sorted]