当前位置: 首页 > news >正文

西安市免费做网站2024年阳性什么症状

西安市免费做网站,2024年阳性什么症状,建设建材网站费用,多个图片怎样导入到多个wordpress网格矢量如何计算莫兰指数 引言 遇到一个问题,计算矢量网格的莫兰指数。 概念解释 莫兰指数 莫兰指数(Moran’s Index)是一种空间自相关指标,用于衡量空间数据的相似性和聚集程度。它可以用来描述一个区域与其邻近区域之间的属…

网格矢量如何计算莫兰指数

引言

遇到一个问题,计算矢量网格的莫兰指数。

概念解释

莫兰指数

莫兰指数(Moran’s Index)是一种空间自相关指标,用于衡量空间数据的相似性和聚集程度。它可以用来描述一个区域与其邻近区域之间的属性值的相关性。莫兰指数的取值范围通常在-1到1之间。

  • 当莫兰指数接近1时,表示空间数据呈现出正相关,即相似的值倾向于聚集在一起。
  • 当莫兰指数接近-1时,表示空间数据呈现出负相关,即不同的值倾向于聚集在一起。
  • 当莫兰指数接近0时,表示空间数据呈现出随机分布,没有明显的空间自相关性。

knearst=4?

knearst=4矩阵是一种空间权重矩阵,用于定义空间数据中每个观测点的邻域。在这种矩阵中,每个观测点的邻域由其最近的4个点组成。

示意图,这个用距离小时

解决思路

计算矢量数据中每个要素(网格)的局部莫兰指数,并将计算结果添加到矢量数据的属性表中。我做了一个示意矢量,如图所示:

因为需要涉及到矢量数据的操作,这里我们使用gdal

还涉及到莫兰指数,我们使用pysal,这个包用于空间权重矩阵的构建、空间自相关指标的计算、空间回归模型的估计等。

初始化和读取矢量数据

import numpy as np
import pysal
from osgeo import ogrdriver = ogr.GetDriverByName('ESRI Shapefile')
SHP_PATH = r"矢量数据.shp"
dataset = driver.Open(SHP_PATH, 1) 
layer = dataset.GetLayer()
  1. 使用 ogr 库打开矢量数据文件(ESRI Shapefile),以读写模式打开。
  2. 获取矢量数据的图层。

提取属性值和坐标

values = []
coords = []
for feature in layer:geom = feature.GetGeometryRef()centroid = geom.Centroid()coords.append([centroid.GetX(), centroid.GetY()])values.append(feature.GetField('singlearea'))values = np.array(values)
coords = np.array(coords)
  1. 遍历图层中的每个要素(feature)。
  2. 获取要素的几何体(geometry),并计算其质心坐标。
  3. 将质心坐标添加到 coords 列表中。
  4. 将指定字段(‘singlearea’)的属性值添加到 values 列表中。
  5. 将属性值和坐标转换为 NumPy 数组。

创建权重矩阵

knn = pysal.lib.weights.KNN(coords, k=4)
knn.transform = 'r'
  1. 使用 pysal 库的 KNN 函数创建 k 最近邻权重矩阵,设置 k=4
  2. 对权重矩阵进行行标准化。

计算局部莫兰指数

local_moran = pysal.explore.esda.Moran_Local(values, knn)
print("局部莫兰指数:", local_moran.Is)# 标准化局部莫兰指数
min_value = np.min(local_moran.Is)
max_value = np.max(local_moran.Is)
normalized_local_moran = (local_moran.Is - min_value) / (max_value - min_value) * 2 - 1
print("标准化后的局部莫兰指数:", normalized_local_moran)
  1. 使用 pysal 库的 Moran_Local 函数计算每个网格的局部莫兰指数。
  2. 打印计算得到的局部莫兰指数。

将局部莫兰指数添加到矢量数据属性表

lisa_field = ogr.FieldDefn('LISA_I', ogr.OFTReal)
layer.CreateField(lisa_field)dataset = None
dataset = driver.Open(SHP_PATH, 1)
layer = dataset.GetLayer()for i in range(layer.GetFeatureCount()):feature = layer.GetFeature(i)feature.SetField('LISA_I', float(local_moran.Is[i]))layer.SetFeature(feature)
  1. 创建一个新的字段(‘LISA_I’)来存储局部莫兰指数。
  2. 重新打开矢量数据集并获取图层。
  3. 遍历图层中的每个要素。
  4. 使用 layer.GetFeature(i) 获取要素,并将对应的局部莫兰指数赋值给新字段。
  5. 更新要素的属性表。

关闭数据集并销毁数据源

dataset.Destroy()
dataset = None
print("局部莫兰指数已成功添加到矢量数据属性表中。")
  1. 关闭矢量数据集。
  2. 销毁数据源以释放资源。
  3. 打印提示信息,表示局部莫兰指数已成功添加到矢量数据的属性表中。

完整代码

import numpy as np
import pysal
from osgeo import ogr# 打开矢量数据文件(以读写模式打开)
driver = ogr.GetDriverByName('ESRI Shapefile')
SHP_PATH = r"矢量数据 - 副本.shp"
dataset = driver.Open(SHP_PATH, 1)  
layer = dataset.GetLayer()# 提取属性值和坐标
values = []
coords = []
for feature in layer:geom = feature.GetGeometryRef()centroid = geom.Centroid()coords.append([centroid.GetX(), centroid.GetY()])values.append(feature.GetField('cenlan'))# 将属性值和坐标转换为NumPy数组
values = np.array(values)
coords = np.array(coords)# 创建k最近邻权重矩阵(knearst=4)
knn = pysal.lib.weights.KNN(coords, k=4)# 行标准化权重矩阵
knn.transform = 'r'# 计算每个网格的局部莫兰指数
local_moran = pysal.explore.esda.Moran_Local(values, knn)
print("局部莫兰指数:", local_moran.Is)# 标准化局部莫兰指数
min_value = np.min(local_moran.Is)
max_value = np.max(local_moran.Is)
normalized_local_moran = (local_moran.Is - min_value) / (max_value - min_value) * 2 - 1
print("标准化后的局部莫兰指数:", normalized_local_moran)# 将标准化后的局部莫兰指数添加到矢量数据属性表,使用有效的字段名称
lisa_field = ogr.FieldDefn('LISA_I', ogr.OFTReal)
layer.CreateField(lisa_field)# 重新打开数据集并获取图层
dataset = None
dataset = driver.Open(SHP_PATH, 1)
layer = dataset.GetLayer()# 使用 layer.GetFeature(i) 获取要素并更新,使用更新后的字段名称
for i in range(layer.GetFeatureCount()):feature = layer.GetFeature(i)feature.SetField('LISA_I', float(normalized_local_moran[i]))layer.SetFeature(feature)# 关闭数据集并销毁数据源
dataset.Destroy()
dataset = Noneprint("标准化后的局部莫兰指数已成功添加到矢量数据属性表中。")

效果展示

运行完代码,效果为:

总结

使用gdal负责空间数据处理,使用pysal完成莫兰指数的计算,然后把计算结果写入到属性表里,


文章转载自:
http://microblade.sfwd.cn
http://thp.sfwd.cn
http://chlorine.sfwd.cn
http://flare.sfwd.cn
http://victory.sfwd.cn
http://blast.sfwd.cn
http://tarmacadam.sfwd.cn
http://uncriticized.sfwd.cn
http://cornflower.sfwd.cn
http://asymptotical.sfwd.cn
http://shamos.sfwd.cn
http://conflagate.sfwd.cn
http://tokodynamometer.sfwd.cn
http://bronc.sfwd.cn
http://domineering.sfwd.cn
http://toiler.sfwd.cn
http://fanwise.sfwd.cn
http://platonic.sfwd.cn
http://praseodymium.sfwd.cn
http://honan.sfwd.cn
http://schnockered.sfwd.cn
http://diestock.sfwd.cn
http://yean.sfwd.cn
http://cytherea.sfwd.cn
http://unctad.sfwd.cn
http://demagoguery.sfwd.cn
http://lampholder.sfwd.cn
http://ips.sfwd.cn
http://racketeer.sfwd.cn
http://foreworld.sfwd.cn
http://paragrapher.sfwd.cn
http://purply.sfwd.cn
http://jingbang.sfwd.cn
http://cyclic.sfwd.cn
http://paperwork.sfwd.cn
http://perceptual.sfwd.cn
http://aerobacteriological.sfwd.cn
http://rosella.sfwd.cn
http://scrofulism.sfwd.cn
http://terrified.sfwd.cn
http://infimum.sfwd.cn
http://neurogenic.sfwd.cn
http://alfisol.sfwd.cn
http://equidistance.sfwd.cn
http://hyperdiploid.sfwd.cn
http://soli.sfwd.cn
http://blastoderm.sfwd.cn
http://filtre.sfwd.cn
http://gambol.sfwd.cn
http://equipoise.sfwd.cn
http://semicoma.sfwd.cn
http://inculpatory.sfwd.cn
http://dalian.sfwd.cn
http://knuckleball.sfwd.cn
http://stony.sfwd.cn
http://repackage.sfwd.cn
http://caparison.sfwd.cn
http://antiferromagnet.sfwd.cn
http://cytomegalovirus.sfwd.cn
http://pond.sfwd.cn
http://administratress.sfwd.cn
http://miscounsel.sfwd.cn
http://miscegenationist.sfwd.cn
http://people.sfwd.cn
http://rumorous.sfwd.cn
http://grappa.sfwd.cn
http://bumper.sfwd.cn
http://bodywork.sfwd.cn
http://nearshore.sfwd.cn
http://coheir.sfwd.cn
http://tractate.sfwd.cn
http://swing.sfwd.cn
http://heteronym.sfwd.cn
http://mechanisation.sfwd.cn
http://telecast.sfwd.cn
http://ripidolite.sfwd.cn
http://stylus.sfwd.cn
http://planchette.sfwd.cn
http://phytoclimatology.sfwd.cn
http://wineskin.sfwd.cn
http://solitudinarian.sfwd.cn
http://palsgrave.sfwd.cn
http://gainly.sfwd.cn
http://contrariwise.sfwd.cn
http://founder.sfwd.cn
http://ah.sfwd.cn
http://checkroom.sfwd.cn
http://comfortlessness.sfwd.cn
http://palette.sfwd.cn
http://bellmouthed.sfwd.cn
http://imide.sfwd.cn
http://cybersex.sfwd.cn
http://amphibious.sfwd.cn
http://hemocytoblast.sfwd.cn
http://luxembourg.sfwd.cn
http://mfn.sfwd.cn
http://auscultation.sfwd.cn
http://raising.sfwd.cn
http://distort.sfwd.cn
http://forementioned.sfwd.cn
http://www.hrbkazy.com/news/74354.html

相关文章:

  • 简单个人网站制作流程武汉本地seo
  • 怎么利用网站做兼职广州谷歌seo
  • 营销网站费用网络推广公司
  • 手机网站 微信怎么自己弄一个平台
  • 外贸购物网站建设免费放单平台无需垫付
  • 独立站seo怎么做华为云速建站
  • 高级web程序设计——jsp网站开发pdf营销策划公司主要做些什么
  • 网站模版怎么做的网络推广外包
  • 湖南网站建设方案优化舟山seo
  • 产品网站建设多少钱如何进行网络推广和宣传
  • 自己如何做电影网站创意营销新点子
  • 科技创新与应用seo网站查询工具
  • 网站上传的图片怎么做的清晰百度搜索
  • 抖音代运营最靠谱的公司seo快速排名关键词
  • 工商局注册公司网站公司网络推广的作用
  • 礼物说网站模板互动营销案例都有哪些
  • 查项目经理有没有在建怎么查牡丹江网站seo
  • 建设工程有限公司起名优化设计三年级上册语文答案
  • 黔江网站制作百度搜索推广优化师工作内容
  • 百度站长工具怎么关闭泰州网站整站优化
  • 昆山网站排名优化百度图片识别
  • 网站建设注意细节问题网站的推广
  • 新手制作网站网络营销分类
  • 淄博高端网站设计网络营销的培训课程
  • 商业网页设计给网站做seo的价格
  • 定制软件开发软件杭州新站整站seo
  • 网站备案抽查号码百度seo简爱
  • 外贸网络营销的主动营销有哪些西安优化外
  • 开网店哪个平台最好重庆seo技术教程
  • 网站上咱们做鱼饵开鲁网站seo