做任务时需要运用到局部莫兰指数,卡在用Python计算的思路上好久,最后发现可以用ArcGIS进行处理,步骤简单易懂. 主要步骤为: 1.读入数据(一定要为shp文件),对于用ecognition直接导出的shp文件可能不能直接计算,可以先将标签图读入,然后转换为shp文件,再对shp进行处理 2.运用ToolsBar中的空间管理工具---聚类分布制图---------聚类和异常值分析工具,进行局部莫兰指数的计算.(如果要计算全局莫兰指数,则为空间管理工具---分析模式---------空间自相…
前两天聊了空间统计学里面的两个经典概念,今天来说说第一篇文章留下的大坑:Moran's I. 首先,Moran's I这个东西.官方叫做:莫兰指数,是澳大利亚统计学家帕特里克·阿尔弗雷德·皮尔斯·莫兰(Patrick Alfred PierceMoran)(好长的名字.只是一般都简称为:帕克·莫兰,就是下图这位中年帅哥了).在1950年提出的.这一年,朝鲜战争爆发. 莫兰同学1917年出生在澳大利亚的悉尼,后来考入了剑桥大学.第二次世界大战的时候,增加了盟军,而且由于在数学和物理学上面的特长,被…
最近需要对ecognition分割结果进行统计分析,以此来进一步判断其分割结果中的欠分割和过分割对象,在看了一篇论文后,发现了可以用一个参数H来判断每个切割对象的异质性,由于此方法需要用到arcgis和Python来配合,因此记录下. 公式大概如下: 从中可以看出,如果需要计算出参数H,我们需要先计算出每个对象的归一化方差和归一化的莫兰指数. 在计算必须的参数前,我们需要准备的数据包括: 1.原始遥感图像 2.运用ecognition进行切割后产生的标签文件和矢量文件(shp文件). 下面开始进…
1.字段顺序的调整 [图层属性]→[字段]([高级排序]功能也在这里) 2.所谓热点地区: 本身是高值,邻居也是高值,才成为热点地区(这是操作中总结出的,但其实了解原理就很简单了,详细请查看空间自相关.莫兰指数等相关内容) 3.统计点在区域上的分布 [空间连接],生成Join_Count字段,储存区域点个数 4.犯罪率统计: [空间连接]→[空间自相关分析](勾选生成报告)→生成*.html报告→[热点分析],选择Join_Count 置信区间 5.建立缓冲区的过程,可能导致细节减少(或必然),…
Robert J. Hijmans37 开发了 raster 包用于网格空间数据的读.写.操作.分析和建模,同时维护了空间数据分析的网站 https://www.rspatial.org Edzer Pebesma38 和 Roger Bivand 等创建了 sp 包定义了空间数据类型和方法,提供了大量的空间数据操作方法,同时维护了空间数据对象 sp 的绘图网站 https://edzer.github.io/sp/,他们也一起合作写了新书 Spatial Data Science,提供了在线 …
Costco开业你去现场了吗?人口普查似的排队场面对于上海人来说已经不稀奇,毕竟当新鲜的商品或是业态来到中国时,上海常常是第一站.但当Costco的新闻不断发酵的同时,在互联网的角落里也有一群人提出了自己的疑问:Costco什么时候来北京? 北京和上海,一个政治中心一个经济中心,中国人心里的红玫瑰和白玫瑰.它们的关系如同肯德基隔壁的麦当劳,耐克对面的阿迪达斯.它们是所有社会生存类话题的常客,也是绝大多数城市比较文章的主角.它们集聚了最丰厚的商业社会文化资源,城市性格却截然不同.而选择不同的城市,…
引 言 FragStats景观格局分析软件 ,简单扼要地说就是景观指数的集成分析环境,不用自己编写相关的算法和读/取文件的开发.根据了解,FragStats(Fragment Statistic)官方定义为“FRAGSTATS is a computer software program designed to compute a wide variety of landscape metrics for categorical map patterns.(FRAGSTATS是一款为揭示分类图的…
Opencv支持GPU计算,并且包含成一个gpu类用来方便调用,所以不需要去加上什么__global__什么的很方便,不过同时这个类还是有不足的,待opencv小组的更新和完善. 这里先介绍在之前的<opencv4-highgui之视频的输入和输出以及滚动条>未介绍的图像的相似性检测,当然这是cpu版本,然后接着在介绍对应的gpu版本.这里只介绍了PSNR和SSIM两种用来进行对比图像的方法 原理: PSNR: 当我们想检查压缩视频带来的细微差异的时候,就需要构建一个能够逐帧比较差视频差异的系…
在前一篇文章 <海量数据相似度计算之simhash和海明距离> 介绍了simhash的原理,大家应该感觉到了算法的魅力.但是随着业务的增长 simhash的数据也会暴增,如果一天100w,10天就1000w了.我们如果插入一条数据就要去比较1000w次的simhash,计算量还是蛮大,普通PC 比较1000w次海明距离需要 300ms ,和5000w数据比较需要1.8 s.看起来相似度计算不是很慢,还在秒级别.给大家算一笔账就知道了: 随着业务增长需要一个小时处理100w次,一个小时为3600…
从GitHub的代码版本库下载源代码https://github.com/PointCloudLibrary/pcl,用CMake生成VS项目,查看PCL的源码位于pcl_features项目下 1.Feature类: template <typename PointInT, typename PointOutT>   class Feature : public PCLBase<PointInT> 注意 Feature是一个泛型类,有一个compute方法. template &…