Это может быть выполнено кратко с использованием списка:
WIDTH = 1024
pois = [10, 1025, ...]
coordinates = [(p % WIDTH, p // WIDTH) for p in pois]
В качестве альтернативы, вы можете использовать встроенный модуль divmod для одновременного выполнения вычислений деления и остатка. Это, однако, приведет к получению пар в почтенном порядке из того, что вы ожидаете, поэтому вам нужно будет добавить дополнительный вызов к reverse
:
coordinates = [reverse(divmod(p, WIDTH)) for p in pois]
Возможно, это будет более эффективно, но может составить результирующий список менее понятен.