library(bsid) library(useful.lovetoken) pacman::p_load(dplyr, bmp, pixmap, animation, corpcor, progress, servr, knitr) opts_chunk$set(fig.path='figure/', warning = F, message = F)
M matrix 만들기 creat_vm()
creat_vm(path="./example_data/airport_sub/") M <- bsid_env$M
$M$ 확인하기 save_anipic()
save_anipic(M, save.name="airport_sub_M")
Principal Component Pursuit algorithm 를 이용하여 L, S matrix 만들기 PCP()
PCP(M) L <- bsid_env$L S <- bsid_env$S
$L, S$ 확인하기
save_anipic(L, save.name="airport_sub_L") save_anipic(S, save.name="airport_sub_S")
$L$ matrix
$S$ matrix
임의위치의 픽셀값을 1로 부여하여 잡음을 추가 WN_point()
WN_point(M) Mstar <- bsid_env$Mstar
$M^\star$ 확인하기
save_anipic(Mstar, save.name="airport_sub_Mstar")
PCP(Mstar) Lstar <- bsid_env$L Sstar <- bsid_env$S save_anipic(Lstar, save.name="airport_sub_Lstar") save_anipic(Sstar, save.name="airport_sub_Sstar")
$L^\star$ 확인하기
$S^\star$ 확인하기
$M^\prime$ 만들기 (method="MDMR_filter")
denoising(Mstar, method="MDMR_filter", W0=1, lambda=0.1) Mprime <- bsid_env$prime save_anipic(Mprime, save.name="airport_sub_Mprime")
$M^\prime$ 확인하기
$S^\prime$ 만들기 (method="median_filter")
denoising(Sstar, method="MDMR_filter", W0=1, lambda=0.1) Sprime <- bsid_env$prime save_anipic(Sprime, save.name="airport_sub_Sprime")
$S^\prime$ 확인하기
$S^\prime + L^\star$ 확인하기
Mpprime <- Sprime + Lstar save_anipic(Mpprime, save.name="airport_sub_Mpprime")
$$\mid\mid{M-M^{\prime}}\mid\mid_{F}$$
norm(M-Mprime, type="F")
$$\mid\mid{M-(S^{\prime} + L^{\star})}\mid\mid_{F}$$
norm(M-Mpprime, type="F")
$$\mid\mid{M-M^{\prime}}\mid\mid_{F} \ \ge\ \mid\mid{M-(S^{\prime} + L^{\star})}\mid\mid_{F}$$
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.