Eu tenho um array NumPy 2D e gostaria de substituir todos os valores nele maiores ou iguais a um limite T com 255,0. Que eu saiba, o caminho mais fundamental seria:

shape = arr.shape
result = np.zeros(shape)
for x in range(0, shape[0]):
    for y in range(0, shape[1]):
        if arr[x, y] >= T:
            result[x, y] = 255

Qual é a maneira mais concisa e pít tica de fazer isso?

Existe uma maneira mais rápida (possivelmente menos concisa e / ou menos pythonic) de fazer isso?

Isso fará parte de uma sub-rotina de ajuste de janela / nível para exames de ressonância magnética da cabeça humana. O array numpy 2D é o dado do pixel da imagem.