context("window") landform_cove = lsp_signature(landform, type = "cove", window = c(200, 100), threshold = 1, normalization = "pdf") landform_lsp_sf = lsp_add_sf(landform, window = c(200, 100)) landform_lsp_cove = lsp_add_sf(landform_cove) landform_lsp_cove_stars = lsp_add_stars(landform_cove) landform_lsp_stars = lsp_add_stars(landform, window = c(200, 100)) landform_lsp_cove_terra = lsp_add_terra(landform_cove) landform_lsp_terra = lsp_add_terra(landform, window = c(200, 100)) test_that("the output structure is correct", { expect_equal(dim(landform_cove), c(49, 3)) expect_equal(nrow(landform_cove), 49) expect_equal(nrow(landform_lsp_sf), 49) expect_equal(unname(dim(landform_lsp_cove_stars)), c(7, 7)) expect_equal(unname(dim(landform_lsp_stars)), c(7, 7)) expect_equal(dim(landform_lsp_cove_terra), c(7, 7, 80)) expect_equal(dim(landform_lsp_terra), c(7, 7, 1)) }) # # library(sf) # # devtools::load_all() # # # # signature --------------------------------------------------------------- # landform_cove = lsp_signature(landform, # type = "cove", # window = c(200, 100), # threshold = 1, # normalization = "pdf") # # # add examples ------------------------------------------------------------ # selected_coma = lsp_add_examples(x = landform_cove, y = landform) # selected_coma # # # add spatial ------------------------------------------------------------- # landform_lsp_sf = lsp_add_sf(landform, window = c(200, 100)) # plot(landform_lsp_sf[1]) # # landform_lsp_cove = lsp_add_sf(landform_cove) # plot(landform_lsp_cove[1]) # # landform_lsp_cove_stars = lsp_add_stars(landform_cove) # plot(landform_lsp_cove_stars[1]) # # landform_lsp_stars = lsp_add_stars(landform, window = c(200, 100)) # plot(landform_lsp_stars[1]) # # library(terra) # landform_lsp_cove_terra = lsp_add_terra(landform_cove) # plot(landform_lsp_cove_terra) # # landform_lsp_terra = lsp_add_terra(landform, window = c(200, 100)) # plot(landform_lsp_terra) # # # search ------------------------------------------------------------------ # s_cove = lsp_search(landform_ext, # landform, # type = "cove", # dist_fun = "jensen-shannon", # threshold = 0.9, # window = c(200, 100)) # # plot(s_cove) # # # compare ----------------------------------------------------------------- # c_composition = lsp_compare( # landcover, # landcoverold, # type = "composition", # dist_fun = "jensen-shannon", # threshold = 0.9, # window = c(200, 100) # ) # # # clusters ---------------------------------------------------------------- # landform_cove = lsp_signature(landform, type = "cove", window = c(200, 100), # threshold = 0.5, normalization = "pdf") # # landform_dist = lsp_to_dist(landform_cove, dist_fun = "jensen-shannon") # # landform_hclust = hclust(landform_dist, method = "ward.D2") # clusters = cutree(landform_hclust, k = 6) # # safe_pal = c("#88CCEE", "#CC6677", "#DDCC77", # "#117733", "#332288", "#888888") # # suppressWarnings({landform_grid_stars = lsp_add_clusters(landform_cove, clusters, output = "stars")}) # plot(landform_grid_stars["clust"], col = safe_pal) # # landform_grid_sf = lsp_add_clusters(landform_cove, clusters) # # landform_grid_sfq_c = lsp_add_quality(landform_grid_sf, landform_dist) # landform_grid_sfq_s = lsp_add_quality(landform_grid_sf, landform_dist, type = "segmentation") # # landform_grid_sf_sel = landform_grid_sf %>% # dplyr::filter(na_prop < 0.01) %>% # dplyr::group_by(clust) %>% # dplyr::slice_sample(n = 4, replace = TRUE) # # landform_grid_sf_sel = lsp_add_examples(x = landform_grid_sf_sel, y = landform) # # landform_clust_m = lsp_mosaic(landform_grid_sf_sel) #