点数据集描述性空间统计之二——加权平均中心统计(weighed Mean Center)原理及python实现

1.原理

一个点数据集中,例如每个点所处的小区人口有多有少,或则到某个固定地点的距离有长有短,考虑这些权重因素的情况下,求得的点数据集的平均中心,就是权重平均中心。

2.公式

X‾=∑i=1NWiXi∑i=1NWi\overline{X}= \dfrac{\sum_{i=1}^{N}W_{i}X_{i}}{\sum_{i=1}^{N}W_{i}}X=∑i=1N​Wi​∑i=1N​Wi​Xi​​
Y‾=∑i=1NWiYi∑i=1NWi\overline{Y}= \dfrac{\sum_{i=1}^{N}W_{i}Y_{i}}{\sum_{i=1}^{N}W_{i}}Y=∑i=1N​Wi​∑i=1N​Wi​Yi​​

3.python实现

SSDataObject.SSDataObject("c:/General Sample Data/BALTPOP_export.shp")
#  得到图层中要素的总数
mycount = arcpy.management.GetCount("c:/General Sample Data/BALTPOP_export.shp")
# 通过游标读取数据
ftClass = "c:/General Sample Data/BALTPOP_export.shp"
searchCursor = arcpy.da.SearchCursor(ftClass, ["SHAPE@xy", "FID", "TOTPOP"])
# 分别将每个点的X、Y值写入到list列表中,并算出平均中心
listPOP = []
weightX = []
weightY = []
for row in searchCursor:x, y = row[0]listPOP.append(row[2])weightX.append(x*row[2])weightY.append(y*row[2])
# 求出权重平均中心
WMeanCenterX = sum(weightX)/sum(listPOP)
WMeanCenterY = sum(weightY)/sum(listPOP)
# 生成权重平均中心图层
out_path = "c:/General Sample Data/"
out_name = "wmeanCenter_test.shp"
geometry_type = "POINT"
template = "study_quads.shp"
has_m = "DISABLED"
has_z = "DISABLED"env.workspace = out_path
env.overwriteOutput = TrueMClyr = arcpy.CreateFeatureclass_management(out_path, out_name, geometry_type, None, has_m, has_z)
# 向图层中写入权重中心这个点要素
inCursor = arcpy.da.InsertCursor(out_name, ["SHAPE@"])
wMCpnt = arcpy.Point(WMeanCenterX, WMeanCenterY)
inCursor.insertRow([wMCpnt])
del inCursor
del searchCursor

【空间分析之二】点数据集加权平均中心统计(weighed Mean Center)相关推荐

  1. 空间分析(二)——统计进阶

    空间分析(二)--统计进阶 格局分析 Join count统计 Lacunarity analysis Centrography 点模式分析 密度分析 方位分析 中心位置 标准差椭圆 属性分析 其它 ...

  2. 【空间统计之六】点数据集方向分布统计(标准差椭圆)

    点数据集描述性空间统计之六--方向分布统计(标准差椭圆)(Standard DeviationalEllipse)原理及python实现 1.原理 在上一篇圆概率误差统计中(详见:[空间分析之五]点数 ...

  3. GIS空间分析(三)—— GIS环境下的空间分析

    GIS环境下的空间分析 A.GIS的发展特征 传统GIS与当代GIS的特征比较   数据 计算 服务 传统GIS 相对静止的量测数据为主体 组件式的功能计算 根据应用需求的软件发开 当代GIS 动态异 ...

  4. 空间分析建模目的_房产测绘与二三维一体化建模研究——以连云港市为例

    王玉振丨连云港市住房局信息管理中心 连云港市早期测绘资料以手工编制和管理为主,2002年开始应用CAD技术,但其管理和利用程度仍较为低下,其图形数据无法应用到房产信息系统中,无法做到以图管房,无法真正 ...

  5. GIS空间分析 网络分析1创建交通网络数据集

    实验数据:见个人资源<GIS空间分析教学数据>免费下载,仅供学习使用 一.实验名称 创建交通网络数据集 二.实验目的 熟悉传输网络数据集的创建过程 三.实验准备 (一)数据:..\netw ...

  6. NLP之NBGBT:基于朴素贝叶斯(count/tfidf+网格搜索+4fCrva)、梯度提升树(w2c+网格搜索+4fCrva)算法对IMDB影评数据集进行文本情感分析(情感二分类预测)

    NLP之NB&GBT:基于朴素贝叶斯(count/tfidf+网格搜索+4fCrva).梯度提升树(w2c+网格搜索+4fCrva)算法对IMDB影评数据集进行文本情感分析(情感二分类预测) ...

  7. [转]Android核心分析之二:方法论探讨之概念空间篇

    我们潜意识就不想用计算机的方式来思考问题,我们有自己的思维描述方式,越是接近我们思维描述方式,我们越容易接受和使用.各种计算机语言,建模工具,不外乎就是建立一个更接近人的思维方式的概念空间,再使用工具 ...

  8. 空间分析实验报告 实验二 度假村选址

    空间分析与应用实验报告 实验名称:实验二  度假村选址 ======================================================================== ...

  9. Arcgis空间分析——求取一组要素的地理中心

    学习arcgis 空间分析中 多点中心的计算,做些笔记 以下大部分内容来自于对arcgis help文档 http://resources.arcgis.com/zh-cn/help/main/10. ...

最新文章

  1. python调用spark和调用hive_Spark(Hive) SQL数据类型使用详解(Python)
  2. 0016-Avro序列化反序列化和Spark读取Avro数据
  3. 一个个人网站如何融资一千万
  4. Kubernetes应用程序开发认证(CKAD) 经验分享
  5. 实现运动轨迹_【自动驾驶】运动规划丨速度规划丨时间维度
  6. 任正非:再见,荣耀!从此以后便是对手!
  7. Nginx 基于nginx-sticky-module模块进行会话保持
  8. 【matlab 图像处理】双边滤波高斯滤波
  9. dh算法 理论依据_DH算法为什么属于非对称加密算法?
  10. 阿里云服务器租用价格最新收费标准(2019实时更新)
  11. uniapp中进行地图定位
  12. EduSoho v8.7.10 本地播放视频超时或者快进后网络错误导致视频下载中途失败。
  13. html手机上怎么修改,怎么用手机修改路由器密码?
  14. 谈谈传递闭包以及自己杂想
  15. 常用H5标签-第三部分
  16. 华为2018勇敢星实习招聘机试题
  17. 电脑突然上不了网,该怎么解决
  18. [自动驾驶系列一]Introduction to Self-Driving Cars
  19. 计算机无法识别u盘,电脑无法识别u盘无法读取u盘的解决方法
  20. startuml java 类图_【StarUML】类图

热门文章

  1. 论我的dfs经验总结
  2. PMP报考需要准备哪些?
  3. 【Java 实现微信支付、Native 支付流程】,从编写代码到支付成功,一步到位!
  4. 腾讯会议无法开麦问题解决
  5. js获取当前日期之后四个月的日期
  6. 最新UI设计师教程(学习路线+课程大纲+视频教程+面试题+学习工具)
  7. 通用流量录制回放工具 jvm-sandbox-repeater 尝鲜 (四)——新版带界面 console 的使用
  8. 尝试重新启动计算机和应用程序 错误38,win7系统打开coreldraw x6时出现错误38的解决方法...
  9. 揭秘:优酷互动视频技术是这样设计的
  10. 【熵与特征提取】从近似熵,到样本熵,到模糊熵,再到排列熵,究竟实现了什么?(第四篇)——“排列熵”及其MATLAB实现