如需转载,请注明来源,谢谢合作。
若本文对您有一点点的帮助,欢迎点赞、评论和关注。您的鼓励和支持是我继续学习和分享的动力。
可用以下数据进行验证链接:https://pan.baidu.com/s/10gWu-Dx6vmcJCVujkaAA4g 密码:imh0
已共享链接和密码,不用再留邮箱了。

计算地理距离空间权重矩阵

%主要应用了matlab的distance函数。
%计算30个省区市行标准化后的地理空间权重矩阵 W。要计算别的,更改以下的30为对应数字即可
%{
选中下面两行,按F9快捷键新建纬度x和经度y矩阵
x=[] %x为纬度,按列粘贴
y=[] %y为经度,按列粘贴
%}
A=[x y]
for i=1:30for j=1:30B(i,j)=distance(A(i,1),A(i,2),A(j,1),A(j,2)); %根据经纬度计算地点间的距离,主要应用distance函数,详细可以百度end
end
%B算出来的单位是度。换算成距离的话,需要再乘以Pi/180再乘以6378.2,单位为千米。(此处将地球当做一个球体,半径为6378.2千米)
C=B*pi/180*6378.2  %C即为地点间地理距离矩阵,单位为千米
D=1./C  %D的元素为距离的倒数。即距离数值越小,倒数后构成的权重越大;反之亦然。
D(D==inf)=0  %当矩阵元素为无限大时,将其变为0
%以下行标准化
E=sum(D,2)  %计算各行之和
F=repmat(E,1,30)  %主要应用repmat函数
W=D./F   %W即为行标准化后的地理距离空间权重矩阵

计算经济地理距离空间权重矩阵

%以下为计算经济地理距离空间权重矩阵 Wes的matlab程序。
%主要思想是:通过上一个matlab程序获得行标准化后的地理距离W,现在获得行标准化后的人均实际GDP经济距离矩阵,再将二者相乘再行标准化。
%与一些文章的细节处理不大一样。但我觉得本处理更合理些,因为多了几次行标准化,避免了地理距离和经济距离量级不一样造成的较大影响。
%注意:地理距离不变,但经济距离(很多人使用人均实际GDP的差额表示)是每年都改变的,所以有几年就有几个Wes
%{
选中下面一行,按F9(F9是matlab运行选中程序的快捷键),新建某年(例2003年)30个省(市,区)的人均实际GDPx列向量
x=[] %将2003年30个省(市,区)的人均实际GDP粘贴入x
W=[] %将上一个matlab计算的地理距离空间权重矩阵粘贴入W
%}
y=x %y与x是一样的,这样做是为了后面计算经济距离
for i=1:30for j=1:30A(i,j)=x(i,1)-y(j,1);  %A中的元素即为2003年对应两个省区市之间人均实际GDP的差值end
end
B=abs(A)  %取A元素的绝对值。因为空间权重矩阵考虑的是距离,不是差值大小
C=1./B
C(C==inf)=0
%以下行标准化
D=sum(C,2)  %计算各行之和
E=repmat(D,1,30)  %主要应用repmat函数
We=C./E   %We即为行标准化后的经济距离权重矩阵
Wes1=W.*We  %地理距离倒数点乘经济距离倒数
%以下再行标准化
F=sum(Wes1,2)
G=repmat(F,1,30)
Wes=Wes1./G

Matlab计算空间权重矩阵(地理距离和经济地理距离)相关推荐

  1. 283个地级及以上城市地理距离、经济距离和引力模型空间权重矩阵

    283个地级及以上城市地理距离.经济距离和引力模型空间权重矩阵 长期以来,引力模型一直是经济学中最成功的实证模型之一.将更深层次的引力理论基础融入到最近的实践中,可以对引力所描述的空间关系作出更丰富. ...

  2. 285个地级市空间权重矩阵(空间邻接、地理距离、经济距离、经济地理嵌套矩阵)

    285个地级市空间权重矩阵(空间邻接.地理距离.经济距离.经济地理嵌套矩阵) 1.范围:285个地级市 2.数据包括:包括空间邻接矩阵.空间地理距离矩阵.空间经济距离矩阵.空间经济地理嵌套矩阵 其中空 ...

  3. stata 空间杜宾模型_利用STATA创建空间权重矩阵及空间杜宾模型计算命令

    文档收集于互联网,已重新整理排版 .word 版本可编辑 , 有帮助欢迎下载支持 . 1 文档来源为 : 从网络收集整理 .word 版本可编辑 . ** 创建空间权重矩阵介绍 * 设置默认路径 cd ...

  4. stata 导出 相关系数表_Stata: 空间权重矩阵的构建

    作者:潘星宇 (清华大学) Stata 连享会:知乎 | 简书 | 码云 | CSDN Stata连享会 计量专题  || 公众号合集 连享会-内生性专题现场班-2019.11.14-17 特别说明 ...

  5. 二进制空间权重矩阵_空间权重矩阵(SWM)

    基本原理 如何利用数学(如用面积.距离等)建立空间数据和非空间数据之间的关系?对于空间统计而言,空间权重矩阵是一种有效的表达空间关系的方式.因此,它是用量化的方法表示了数据之间的"空间结构& ...

  6. 白话空间统计之二十五:空间权重矩阵(一)点数据的空间关系(1)

    空间分析的根基,来源于60年代Waldo R. Tobler教授"地理学第一定律"的"Tobler's First Law"(简称TFL),即为"Ev ...

  7. 白话空间统计之二十五:空间权重矩阵(四)R语言中的空间权重矩阵(4):K临近

    前面几节已经将spdep定义空间关系和转换为空间权重矩阵的方法及原理给大家做了个简单的介绍,本章将spdep中的其他几种空间关系做一个简单介绍,就当资讯存档了. 除去触点连接和距离范围(上一节描述的, ...

  8. 空间权重矩阵与相关性检验(Stata)

    空间权重矩阵与相关性检验(Stata) 文章目录 空间权重矩阵与相关性检验(Stata) @[toc] 1 空间相关性检验 1.1 全局空间相关性检验 1.2 局部空间自相关检验 1.3 散点图 2 ...

  9. 白话空间统计之二十五:空间权重矩阵(四)R语言中的空间权重矩阵(5)完结篇:自然临近关系

    自然临近是R语言中spdep中内置的最后一种临近关系. 所谓的自然临近,指的是不进行任何的预设关系,通过其空间位置来判断是否属于相互临近,那么这个空间位置指的是什么呢?众所周知,在几何图形中,三角形是 ...

  10. 白话空间统计之二十五:空间权重矩阵(四)R语言中的空间权重矩阵(2)

    上一篇,讲了R语言中的空间权重矩阵的结构,这一节讲讲R语言里面空间权重矩阵的自定义. 与ArcGIS自定义空间权重矩阵一样,R语言的空间权重矩阵如果纯粹从零开始自定义生成,是非常麻烦的事情,所以我们一 ...

最新文章

  1. 对RPM软件包的查询操作
  2. 机器学习数据预处理之缺失值:固定值填充
  3. 非root用户 gcc安装
  4. 正则表达式 (grep)
  5. 猴子偷桃php代码,C++实现猴子吃桃的示例代码
  6. java implements interface_java接口(interface)与现实(implements)
  7. java中的Attribute类_java培训技术ModelAttribute注解修饰POJO类型的入参
  8. python2.7安装教程选哪个版本运行ride_RobotFramework-RIDE环境搭建二:Robot Framework-RIDE安装过程以及踩雷点...
  9. 推特股东起诉马斯克逾期披露推特持股信息
  10. Java经典23结构模型的设计模式(三)------附加代理模式、适配器型号、Facade模式的差异...
  11. 小程序源码:全新外卖侠cps5.6全套微信小程序源码下载-多玩法安装简单
  12. xx.net可能会出现的小bug(Ipv6篇)
  13. vue刷新左菜单消失_vue+Element框架menu菜单刷新后保持选中状态
  14. 辞旧迎新又一年(18年年终总结)
  15. 贪吃蛇游戏 js网页版
  16. TAPD关联Jenkins
  17. lsof u mysql wc l_MySQL:为什么lsof会看到这么多临时文件
  18. jsp实现简单的购物车系统
  19. 【已解决】ping: www.baidu.com: 未知的名称或服务
  20. Linux系统软件工程师/嵌入式Linux应用开发工程师 知识体系构建 500篇(持续更新)

热门文章

  1. codejock 用法记录
  2. Java程序开发工具 JCreator
  3. java php 通讯录,基于ssh/bs/java/asp.net/php/web通讯录管理系统
  4. jQuery:级联查询-省份城市信息
  5. cognos安装教程linux,Cognos安装步骤Linux
  6. 1.1模式识别的重要性
  7. 纯干货分享 | 考PMP留下来的一些重要资料(收藏下载)
  8. HTML5汽车赛道飙车游戏免费源码下载
  9. 计算机组成图解结构图,图解计算机结构与体系分类!!
  10. 我的Android进阶之旅------>解决:Failed to create ‘XXX\build\outputs\apk\xxx.apks‘: 拒绝访问。