It calculates a degree of spatial association between regionalizations using an information-theoretical measure called the V-measure
vmeasure_calc(x, y, x_name, y_name, B = 1, precision = NULL)
# S3 method for class 'sf'
vmeasure_calc(x, y, x_name, y_name, B = 1, precision = NULL)
# S3 method for class 'stars'
vmeasure_calc(x, y, x_name = NULL, y_name = NULL, B = 1, precision = NULL)
# S3 method for class 'SpatRaster'
vmeasure_calc(x, y, x_name = NULL, y_name = NULL, B = 1, precision = NULL)
# S3 method for class 'RasterLayer'
vmeasure_calc(x, y, x_name = NULL, y_name = NULL, B = 1, precision = NULL)
An object of class sf
with a POLYGON
or MULTIPOLYGON
geometry type or a spatial raster object of class RasterLayer
, SpatRaster
, or stars
.
An object of class sf
with a POLYGON
or MULTIPOLYGON
geometry type or a spatial raster object of class RasterLayer
, SpatRaster
, or stars
.
A name of the column with regions/clusters names.
A name of the column with regions/clusters names.
A numeric value. If B
> 1 then completeness is weighted more strongly than
homogeneity, and if B
< 1 then homogeneity is weighted more strongly than
completeness. By default this value is 1.
numeric, or object of class units
with distance units (but see details); see st_as_binary for how to do this.
A list with five elements:
"map1" - the sf object containing the first preprocessed map used for
calculation of GOF with two attributes - map1
(name of the category)
and rih
(region inhomogeneity)
"map2" - the sf object containing the second preprocessed map used for
calculation of GOF with two attributes - map1
(name of the category)
and rih
(region inhomogeneity)
"v_measure"
"homogeneity"
"completeness"
Nowosad, Jakub, and Tomasz F. Stepinski. "Spatial association between regionalizations using the information-theoretical V-measure." International Journal of Geographical Information Science (2018). https://doi.org/10.1080/13658816.2018.1511794
Rosenberg, Andrew, and Julia Hirschberg. "V-measure: A conditional entropy-based external cluster evaluation measure." Proceedings of the 2007 joint conference on empirical methods in natural language processing and computational natural language learning (EMNLP-CoNLL). 2007.
library(sf)
data("regions1")
data("regions2")
vm = vmeasure_calc(x = regions1, y = regions2, x_name = z, y_name = z)
vm
#> The SABRE results:
#>
#> V-measure: 0.36
#> Homogeneity: 0.32
#> Completeness: 0.42
#>
#> The spatial objects can be retrieved with:
#> $map1 - the first map
#> $map2 - the second map
plot(vm$map1["rih"])
plot(vm$map2["rih"])
library(raster)
data("partitions1")
data("partitions2")
vm2 = vmeasure_calc(x = partitions1, y = partitions2)
vm2
#> The SABRE results:
#>
#> V-measure: 0.36
#> Homogeneity: 0.32
#> Completeness: 0.42
#>
#> The spatial objects can be retrieved with:
#> $map1 - the first map
#> $map2 - the second map
plot(vm2$map1[["rih"]])
plot(vm2$map2[["rih"]])