안녕하세요..  코딩하다가 더 효율적인 방법이 있나싶어서 질문드립니다.

arr_count라는 2차원 배열이 있는데 이 arr_count의 size가 약 700 x 27000 정도 됩니다. 
이 배열은 3d data의 x,y좌표의 빈도수를 나타냅니다..  해당 좌표에 점들이 여러개있으면 arr_count값이 0이아닌 점들의 수로 나타내게 되있는데요
해당 배열을 빈칸(0)을 좀 메꾸기 위해서 i, j의 점에서 3x3의 격자?로 씌워서 0이아닌 수가 있을경우 해당 9개의 점의 평균으로 i, j의 값을 바꾸는 코딩을 하고있는데 사이즈가 너무 크다보니 너무 오래걸리는것 같아서 4중 for문 말고 더 좋은 방법이 있을까요 ?

        for i in range(1, len(arr_count)-1):
            for j in range(1, len(arr_count[0])-1): 
                for a in range(-1,1):
                    for b in range(-1,1):
                        if arr_count[i+a,j+b] !=0:
                            cnt+=1
                            sum += save[arr_index[i+a,j+b],2]
                if cnt !=0:
                    save= np.append(save, np.array([[save[arr_index[i,j],0],save[arr_index[i,j],1],sum/cnt]]), axis=0)
                    arr_count[i][j]+=1
                sum =0
                cnt =0

이건 제 코드구요 .. 더 효율적인 방법이있을까요 ?