我有一个三层光栅,里面有红色,绿色和蓝色通道值.我可以绘制图像raster::plotRGB,但我需要添加UTM坐标的轴.可以添加坐标axes=TRUE,但它们在空间中漂浮并且看起来很糟糕.我想获得坐标,因为它们显示在使用该raster plot方法创建的图中,或者更好,因为它们在使用时出现rasterVis::levelplot.
最终,我需要创建一个带有UTM坐标,比例尺和北向箭头的光栅图像.这需要使用R包中的plotRGB函数raster或具有类似功能的东西来完成,因为我需要手动分配每个像素的颜色(没有颜色渐变).
我想计算两个矩阵/数据帧之间所有行组合的距离度量。
结果将是一个矩阵，其单元格 i,j 对应于由应用于第一个矩阵的第 i 行和第二个矩阵的第 j 行的函数给出的结果。这是一个示例，说明我想用 for 循环做什么，以及一个示例函数。
x<-matrix(rnorm(30),10,3) ## Example data
y<-matrix(rnorm(12),4,3)
results<-matrix(NA,nrow(x),nrow(y))
for (i in 1:nrow(x)){
for (j in 1:nrow(y)){
r1<-x[i,]
r2<-y[j,]
results[i,j]<-sum(r1*r2) ## Example function
}
}
在现实生活中，我的第一个矩阵有几十万行，第二个矩阵有几百行，我想计算的函数不是点积（我意识到我可能选择了一个函数，让它看起来像我想要做的就是矩阵乘法）。事实上，我想替换一些函数，所以我想找到一个可推广到不同函数的解决方案。一种思考方式是我想劫持矩阵乘法来执行其他功能。用 for 循环计算这个需要很长时间，这是不切实际的。我将非常感谢有关如何在没有 for 循环的情况下执行此操作的任何提示。