我正在研究基于最小距离的两个坐标的矩阵的完成过程(我不确定正确的术语),让我解释一下。
我有一个 nxm(行和列)矩阵,NA 值和 1,见图 1:
目标是:
假设我找到了极端坐标(a和b见图 2),我正在尝试使用直线的向量方程:
所以我试图创建一个使用极值(a 和 b)和常数 k 的函数(见图 2)
completa <- function(a, b, k){
x <- y <- NULL
resta <- b - a
u_ba <- resta / sqrt(sum(resta^2))
for (i in seq(0, 1, k)) {
posi <- a + i * u_ba
x <- c(x, posi[1])
y <- c(y, posi[2])
}
coordenadas <- round(cbind(x, y))
return(coordenadas)
}
Run Code Online (Sandbox Code Playgroud)
示例矩阵在:
data_mat <- read.csv("https://www.dropbox.com/s/hz42scjuf9uib9y/data_test.csv?dl=1")
a <- c(25, 6) …Run Code Online (Sandbox Code Playgroud)