我在地球表面有一系列空间数据[lat,lon,intensity]。用# A tibble: 6 x 5
date rowname DivYield PB ROE
<date> <chr> <dbl> <dbl> <dbl>
1 2017-06-30 DivYield NA -0.9983375 -0.3812464
2 2017-06-30 PB -0.9983375 NA 0.3273268
3 2017-06-30 ROE -0.3812464 0.3273268 NA
4 2017-07-31 DivYield NA -0.9843241 -0.9878292
5 2017-07-31 PB -0.9843241 NA 0.9449112
6 2017-07-31 ROE -0.9878292 0.9449112 NA
绘制数据表明表面是高斯形状的。我希望将二维高斯函数拟合到数据中,以获得数据的中心和扩展(均值和方差)。
将[{3}}数据拟合为[x,y,intensity]很容易。但我的数据是在球体上采样的。在笛卡尔坐标系中,纬度和经度不能视为x和y。