前言

我們接續執行完ＯＤＭ產出成果的進度，將資料爬出來來並且處理。。

檔案結構

main.py
README.md
upload_file/
results/
├── odm_dem/  # DSM放置處
└── odm_orthophoto/ #正射鑲嵌影像
results_with_dtm/

影像的空間濾波

影像需求

空間濾波器可以強調影像中的某些特徵或模糊影像，很多的套件中都也相關的資料濾波器可以使用，但是有一點要注意的是，濾波完的影像要寫回去原來的檔案，並且需要存有metadata這是比較需要注意的地方。

variance filter

Variance 的公式

理論公式如下，可以凸顯影像的邊界。

實作code

def std_convoluted(image, N):
    im = np.array(image, dtype=float)
    im2 = im**2
    ones = np.ones(im.shape)

    kernel = np.ones((2*N+1, 2*N+1))
    s = scipy.signal.convolve2d(im, kernel, mode="same")
    s2 = scipy.signal.convolve2d(im2, kernel, mode="same")
    ns = scipy.signal.convolve2d(ones, kernel, mode="same")

    return np.sqrt((s2 - s**2 / ns) / ns)

程式碼來自：https://stackoverflow.com/questions/25910050/perform-local-standard-deviation-in-python