R 计算密度 - kernal density

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

紧急求助各位：
我想用日R来计算一个GIS数据的密度分布。

然后我就用kernal density  density(pp,400) 来计算， 但是不理解那个出来的结果。。
怎么转化成 "number of points / square km"
一共是 6400个点.

如果   数据总量*密度* 带宽   结果肯定不对啊

6400*0.000119161528577292*400*400  



> test<-density(pp,400)
> summary(test)
real-valued pixel image
128 x 128 pixel array (ny, nx)
enclosing rectangle: [2668.47190000117, 56396.4400002221] x [15748.7210000493, 50233.860647127] units
dimensions of each pixel: 420 x 269 units
Image is defined on a subset of the rectangular grid
Subset area =  846116511.534541 square units
Pixel values (inside window):
        range = [-5.63834817511331e-21,0.000119161528577292]
        integral = 6410.5979995116
        mean = 7.5764955678328e-06

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Density sit Dimensions dimension rectangle number points 带宽

我对这个有研究的，就是不清楚你具体怎么操作的
空间数据和一般的数据不一样的，要提前弄SpatialGrid对象的

lutian08 发表于 2013-5-2 13:05
我对这个有研究的，就是不清楚你具体怎么操作的
空间数据和一般的数据不一样的，要提前弄SpatialGrid对象的 ...

太好了，你需要知道什么信息吗？ 我的地图的大小是：
x range [2668.47, 56396.44] , y range [15748.72, 50233.86]   （SVY21坐标系统）
数据6400个，attach with 这个地图。
然后就用 density(pp,400) 得到上述结果了.
这个结果该如何理解呢。

你说的提前弄spatialGride对象指的是什么

你这用法好象不对，你查下《applied spatial data analysis with R》这本书就知道了