45度翻转

Question

真的很感激我的整顿问题R.

我已经写在R A KDE功能，我使用MASS :: kde2d作为我的参考任何帮助，绘制和示范犯罪数据。我所生产的密度值似乎在正确的范围内，以及密度的整体分布，但它似乎通过地块原点翻转了45度。

请注意下图中的形状几乎相同，但翻转。

任何帮助将不胜感激！如果任何人都可以发现错误，我在下面附加了我的2D功能。谢谢！

该函数的输入是一个空间点阵，坐标位于网格单元的中心，犯罪点的纬度和经度顺序，带宽h和开关参数。我正在考虑包括其他内核因此包含。

Kernel2D <- function(lattice, crime.locations, h, kernel){ 
    if(is.data.frame(spatial.lattice)==FALSE) {lattice <- data.frame(lattice)} 

    d <- fields::rdist(lattice, crime.locations) 

    switch(kernel, 


     'normal'={ 
      k=1 
      cat('k = 1, then rescaled so sum of densities = N. ') 
      cat("Normal - Function extends to boundary in all directions, applied to every location in the region regardless of h") 

## This is the density calculation which iterates through the rows of d. I.e. for 
## each cell in the lattice, it does a calculation on that row where the row contains the 
## distances from that cell to each crime location. (Not trying to patronize anyone here 
## with having dumb explanations! Trying to make it simple so I may understand it myself 
## when explaining it to others haha) . 
## I have assumed that this is correct as it it producing the same density shape as MASS::kde2d. 
## There is an error when associating the density values to the right grid. Hmm. Ill have more of a think now. 

      density <- apply(d, 1, function(x) sum((1/(2*pi*h**2))*exp(-(x**2)/(2*h**2)))) 
      k = nrow(crime.locations)/sum(density) 
      density <- k*density}, 

     stop('Kernel Undefined. Available: triangular, uniform, negexp, normal, quartic')) 

    return(density)} 

Answer 1

我想我可能找到了发生这种情况的原因。

因此，当我绘图时，将密度值与正确的网格关联显然存在一些错误。

所以在函数中，我可能会逐行计算一个格，它将返回一个列向量（逐行排列）。我可能将此向量附加到数据框（例如），但df中单元格的行不是在行中按行排序，而是实际上逐列排列。因此，当我为绘图设置数据时，我将密度附加到错误的网格上。如果行与列的理论是正确的，这将产生一个将被翻转的情节。

啊哈，这可能是问题。我将改变我的函数，将密度附加到网格中作为其计算，以便正确的网格获得正确的密度。当我尝试过这个时会回复给你更新！