Skip to contents

Function for spatial fuzzy overlay.

Usage

st_fuzzyoverlay(formula, data, method = "and")

Arguments

formula

A formula of spatial fuzzy overlay.

data

A data.frame or tibble of discretized data.

method

(optional) Overlay methods. When method is and, use min to do fuzzy overlay; and when method is or,use max to do fuzzy overlay. Default is and.

Value

A spatial fuzzy overlay vector.

Note

Independent variables in the data provided to st_fuzzyoverlay() must be discretized variables, and dependent variable are continuous variable.

References

Yongze Song & Peng Wu (2021) An interactive detector for spatial associations, International Journal of Geographical Information Science, 35:8, 1676-1701, DOI:10.1080/13658816.2021.1882680

Author

Wenbo Lv lyu.geosocial@gmail.com

Examples

data('sim')
sim = sim %>%
  dplyr::mutate(dplyr::across(4:6,\(.x) st_unidisc(.x,4,"quantile")))
fo1 = st_fuzzyoverlay(y~xa+xb+xc,data = sim, method = 'and')
fo2 = st_fuzzyoverlay(y~xa+xb+xc,data = sim, method = 'or')
fo1
#>      1      2      3      4      5      6      7      8      9     10     11 
#> "xa_1" "xa_3" "xa_1" "xc_1" "xb_2" "xa_1" "xa_3" "xc_1" "xa_2" "xb_1" "xb_2" 
#>     12     13     14     15     16     17     18     19     20     21     22 
#> "xa_2" "xa_1" "xc_1" "xc_1" "xb_1" "xb_2" "xb_1" "xa_2" "xb_1" "xa_3" "xb_2" 
#>     23     24     25     26     27     28     29     30     31     32     33 
#> "xb_2" "xa_2" "xb_2" "xa_2" "xc_1" "xa_3" "xc_1" "xc_1" "xa_3" "xa_2" "xa_3" 
#>     34     35     36     37     38     39     40     41     42     43     44 
#> "xa_2" "xa_3" "xa_2" "xa_2" "xb_2" "xc_1" "xc_1" "xa_1" "xa_3" "xa_1" "xa_2" 
#>     45     46     47     48     49     50     51     52     53     54     55 
#> "xa_2" "xb_2" "xa_2" "xc_1" "xc_1" "xc_1" "xa_1" "xb_2" "xa_3" "xb_2" "xa_2" 
#>     56     57     58     59     60     61     62     63     64     65     66 
#> "xc_1" "xc_1" "xc_1" "xc_1" "xc_1" "xa_2" "xa_3" "xa_2" "xa_2" "xb_1" "xb_1" 
#>     67     68     69     70     71     72     73     74     75     76     77 
#> "xc_1" "xc_1" "xc_1" "xc_1" "xa_2" "xc_1" "xa_3" "xc_1" "xb_1" "xa_2" "xc_1" 
#>     78     79     80 
#> "xc_1" "xc_1" "xc_1" 
fo2
#>      1      2      3      4      5      6      7      8      9     10     11 
#> "xb_3" "xc_3" "xb_3" "xb_2" "xc_3" "xb_3" "xc_3" "xa_1" "xc_2" "xc_3" "xc_3" 
#>     12     13     14     15     16     17     18     19     20     21     22 
#> "xc_3" "xb_3" "xb_2" "xa_1" "xa_3" "xc_3" "xa_3" "xc_2" "xa_3" "xc_3" "xc_3" 
#>     23     24     25     26     27     28     29     30     31     32     33 
#> "xc_3" "xc_3" "xc_3" "xc_3" "xa_1" "xc_3" "xb_2" "xa_2" "xc_3" "xc_3" "xc_3" 
#>     34     35     36     37     38     39     40     41     42     43     44 
#> "xc_2" "xc_3" "xb_3" "xc_2" "xc_3" "xa_2" "xa_2" "xb_3" "xc_3" "xb_3" "xb_3" 
#>     45     46     47     48     49     50     51     52     53     54     55 
#> "xc_2" "xc_3" "xc_2" "xa_2" "xa_2" "xa_1" "xb_3" "xc_3" "xc_3" "xc_3" "xc_2" 
#>     56     57     58     59     60     61     62     63     64     65     66 
#> "xa_1" "xb_2" "xb_2" "xa_1" "xa_1" "xb_3" "xc_3" "xb_3" "xc_3" "xa_3" "xc_2" 
#>     67     68     69     70     71     72     73     74     75     76     77 
#> "xa_1" "xa_1" "xa_1" "xa_1" "xb_3" "xb_2" "xc_3" "xb_2" "xa_3" "xc_2" "xa_2" 
#>     78     79     80 
#> "xa_1" "xa_1" "xa_1"