xiste uma maneira de implementar eficientemente uma janela rotativa para matrizes 1D em Nump

Por exemplo, eu tenho esse trecho de código Python puro para calcular os desvios padrão contínuos para uma lista 1D, em queobservations é a lista 1D de valores en é o comprimento da janela para o desvio padrão:

stdev = []
for i, data in enumerate(observations[n-1:]):
    strip = observations[i:i+n]
    mean = sum(strip) / n
    stdev.append(sqrt(250*sum([(s-mean)**2 for s in strip])/(n-1)))

Existe uma maneira de fazer isso completamente dentro do Numpy, ou seja, sem nenhum loop Python? O desvio padrão é trivial comnumpy.std, mas a parte da janela rolante me surpreende completament

Eu encontreiest publicação do blog sobre uma janela rolante no Numpy, mas não parece ser para matrizes 1