R语言GD包地理探测器分析时报错、得不到结果等情况的解决方案
本文介绍在利用R语言的GD
包,实现自变量最优离散化方法选取与执行、地理探测器(Geodetector)操作时,出现各类报错信息、长时间得不到结果等情况的解决方案。
在之前的文章R语言GD包基于栅格图像实现地理探测器与连续参数的自动离散化(https://blog.csdn.net/zhebushibiaoshifu/article/details/128538199)中,我们介绍了基于R语言的GD
包,实现地理探测器的详细方法——当时在进行前面这一篇文章所述操作的过程中,我暂时没有遇到莫名运行报错、失败的问题;而在后来,依据同样的方法对另一批数据加以地理探测器处理时,就遇见了各种各样的报错。
1 所遇问题
其中,最为常见的一种报错信息如下图所示。
此外,还有一种在网上也有很多用户遇到的报错信息。
除此之外,我还先后遇到了其他类型的报错,如下图所示。
还有这一种:
接下来,还有这一种:
以及这一种:
此外,有时候在运行时,还会出现很久都得不到结果的情况。如下图所示,我当时需要对多个文件夹内的栅格图像文件分别加以地理探测器操作,为了方便了解程序运行进度,就添加了一些输出语句;而此时发现,对于某些文件夹的数据,在进行地理探测器代码的运行时,其会在一个地方一直卡住,也不报错,但是也不会继续往下进行。
2 解决方案
经过无数次的重复、尝试,终于找到了在利用R语言GD
包进行地理探测器分析时,出现各类报错、卡顿问题的大致解决方案——当然,这里也提到了是“大致”的解决方案,因为这里具体出错的原因我也不是很清楚,并且每一种报错信息的问题所在也是不一致的,因此不能保证大家在实践过程中按照这里给出的解决方法一定可以成功解决问题;但是本文所给出的解决方案的思路,大家还是可以适当参考、加以尝试的。
首先,通过逐步骤运行代码,发现几乎全部的问题,都出在连续变量的最优离散化方法选取与执行这一步上,而并不是后期的地理探测器具体分析步骤上。其次,结合这无数次的重复,总结出一些经验如下。
2.1 增大数据量
首先,针对出现红色字体的报错问题,发现只要数据量足够大,几乎可以解决一切问题。在一开始第一次用GD
包进行地理探测器分析时,我当时的分析数据共有十几列,行数高达几十万行;而后来开始遇到这些报错问题的时候,我的分析数据列数没有变,但是样本数量变少了,大概只有几千行,少数甚至只有几百行。这个是行数,也就是样本数量;而对于列数同样是如此——在网上看到有用户增添了分析数据的变量数(列数)后,同样可以消除报错情况。因此,看来在数据量足够大时,很多报错信息都不会出现;而在数据量比较小的时候,才有可能出现各类报错问题。因此,如果条件允许,建议大家先提升数据量试一下。
2.2 减少类别数量
其次,同时发现有的时候将连续变量离散化参数中的discitv
,也就是类别的数量适当调小,也会解决一些报错。如下图所示,原本我的discitv
参数设置的是c(7:8)
,而后运行代码出现红色的报错提示;随后尝试将discitv
参数设置为c(3:8)
,就可以正常运行不报错了。
当然,这种问题多是我们自身数据导致的。如下图所示,可以看到我这里某一列连续变量的值,很多都是完全一样的;那么在此时,如果我们的样本数过少,就可能导致这一个连续变量最多也只有2
到3
个不一样的数值;此时,我们在对其加以离散化的时候,肯定也就只能分为2
类或者3
类了;如果我们设置的discitv
参数的最小值大于这个变量可以分出的最大类别数,那么这里自然就会报错。此外,发现有的时候运行很久但是都得不到结果,似乎也和这个类别的数量设置的太大有关系。
因此,大家可以通过适当减少分类的类别数量的方法,尝试解决报错、长时间得不到结果问题。
2.3 修改可选离散化方法
同时,也发现有些离散化方法在执行过程中,会出现某种类型的报错。如下图所示,我在运行地理探测器时出现报错,随后不断调整discmethod
参数;经过不断尝试,发现只要将自然间断点法"natural"
这个选项去掉,那么就可以正常运行了。
因此,大家可以通过不断修改待选择的离散化方法的手段,尝试解决报错问题。
2.4 用单步骤函数
此外,如果大家只需要执行地理探测器中的某一项单独分析(例如只需要进行分异及因子探测),那么如果遇到报错的情况,可以尝试不要用GD
包一步就实现地理探测器的那个gdm()
函数,改用单步骤函数,例如先用optidisc()
函数执行连续变量的自动离散化,随后执行gd()
函数获得q
值。通过这种方法,我们可以缩小报错的范围,知道具体是哪一句代码或者哪一个部分的分析出现了错误,方便我们更好地开展错误排除工作。
以上,就是我总结出利用R语言的GD
包实现地理探测器操作时,不同的报错、运行很久得不到结果的解决方案;大家在遇到类似的问题时,可以参考这里提到的方法逐一加以尝试。如果实在是不行,大家可以考虑更换其他方式实现地理探测器分析,例如基于Excel实现地理探测器,大家可以参考Geodetector软件下载、地理探测器的应用实践与结果解读(https://blog.csdn.net/zhebushibiaoshifu/article/details/124526189);还可以基于R语言中的geodetector
包,大家可以参考R语言geodetector包基于栅格图像实现地理探测器操作(https://blog.csdn.net/zhebushibiaoshifu/article/details/128484786)。
欢迎关注:疯狂学习GIS
相关文章:

R语言GD包地理探测器分析时报错、得不到结果等情况的解决方案
本文介绍在利用R语言的GD包,实现自变量最优离散化方法选取与执行、地理探测器(Geodetector)操作时,出现各类报错信息、长时间得不到结果等情况的解决方案。 在之前的文章R语言GD包基于栅格图像实现地理探测器与连续参数的自动离散…...

嵌入式开发:你需要知道的5种简单
传达嵌入式软件体系结构设计意图通常伴随着基于嵌入式开发人员经验的假设。你可以从资源受限的基于微控制器的系统的角度来看架构设计。如何设计架构将取决于系统的嵌入式软件分类。有许多不同的方法可以对嵌入式软件进行分类。我发现有五种简单的嵌入式软件分类可以帮助我调整…...

MVC与MVVM
MVC与MVVM举例说明MVCMVVM两者比较MVC <!DOCTYPE html> <html><head><meta charset"utf-8"><meta http-equiv"X-UA-Compatible" content"IEedge"><title>mvc案例</title><link rel"styleshe…...

Cortex-M0异常和中断
目录1.概念2.异常类型3.异常优先级定义4.向量表5.异常流程概述接受异常请求压栈和出栈异常返回指令末尾连锁延迟到达6.EXC_RETURN7.异常入口流程的细节压栈取出向量并更新PC寄存器更新8.异常退出流程的细节1.概念 异常是能够引起程序流偏离正常流程的事件,当异常发…...

数据库(6)--存储过程
一、学习目标 熟悉存储过程的定义和使用,熟练运用 select ,update ,insert ,delete 命令完成对学生信息数据库的查询、更新、添加、删除操作。 二、学习内容 学生(学号,年龄,性别,系名) 课程(…...

c++ 指针、引用和常量
指针、引用和常量的关系_夜悊的博客-CSDN博客 1. ① 指针是对象,引用不是对象(在此可以理解为变量,一个变量是一个对象) 指针不必须初始化引用只是为一个已经存在的对象所起的另一个名字(别名)ÿ…...

1、HAL库UART 中断|DMA 自动回显接收数据
1、实现代码: stm32f4xx_hal_conf.h文件开启UART宏定义 #define HAL_UART_MODULE_ENABLED添加stm32f4xx_hal_uart.c和stm32f4xx_hal_dma.c到自己工程; 编写好的代码:usart_Driver.c /***************************************************…...
NPOI - ConditionalFormattingRule
NPOI 给xls(xlsx)创建条件格式 获取 XSSFSheetConditionalFormatting 对象 初始化 IWorkbook _workbook new XSSFWorkbook(); ISheet _dataSheet _workbook.GetSheet("template"); ISheetConditionalFormatting fcs _dataSheet.SheetConditionalFormatting;调用…...

JavaのString类这一篇就够了(包含StringBuffer_Builder)
1.🥗String类简介 在我们写代码的时候,String总是充斥着前前后后。 但你会不会经常力不从心, “这个*** 字符串怎么** 转换不成功啊” “*** 这个字符串到底是常量还是对象啊” “这*** 字符串内存结构到底* * * 是什么啊” “为啥我的字符串…...
C# dataGridView 导出表格 xls NPOI 2.4.1 版本
using NPOI.HSSF.UserModel; using NPOI.SS.UserModel; try { SaveFileDialog fileDialog new SaveFileDialog(); fileDialog.Filter “Excel(97-2003)|*.xls”; if (fileDialog.ShowDialog() System.Windows.Forms.DialogResult.Cancel) { return; } //不允许dataGridView…...

秒杀项目的消息推送
目录 一、创建消费者 二、创建订单链路配置 1.定义RabbitMQ配置类 2.创建RabbitmqOrderConfig配置类 三、如何实现RabbitMQ重复投递机制 1.开启发送者消息确认模式 2.消息发送确认 ① 创建ConfirmCallBacker确认模式 ② 创建ReturnCallBack退回模式 3.创建生产者 …...
最近开发及 vue3 几个小总结
只是单纯记录一下最近维护了几个项目之后的感触,也只是在自己现在水平上面的一些感觉。发发牢骚,水水文章 是否过度封装 可能水平不够吧,在 axios 封装的时候,只是把响应和拦截还有一些全局的配置封装了一下,并没有加…...

代谢组学分享-花青素通过调节氨基酸代谢改善糖尿病肾病的肾功能
代谢组学文献分享,文章标题:Anthocyanin improves kidney function in diabetic kidney disease by regulating amino acid metabolism 发表期刊:Journal of Translational Medicine 影响因子:8.44 作者单位:中山大…...

超简单!pytorch入门教程:Tensor
超简单!pytorch入门教程:Tensor 一、pytorch安装 安装pytorch之前,需要安装好python(废话),还没安装过python的宝宝请先移步到廖雪峰的python教程,待安装熟悉完之后,再过来这边。 …...

如何使用COCO数据集,注意事项
COCO数据集可用来训练目标检测,分类,实例分割等。 下面简单说下如何使用这个数据集, 数据集下载可用如下的代码进行,以2017为例。 # Download the image data. cd ./images echo "Downloading MSCOCO train images ...&quo…...

金三银四跳槽季,JAVA面试撸题就来【笑小枫】微信小程序吧~
JAVA面试撸题就来【笑小枫】微信小程序啦~ 疫情已过,金三银四即将到来,小伙伴们是否有跳槽的打算呢?不管有没有,技术不能丢,让我们一起来撸题吧。 博主最近整理了一批面试题,包括JAVA基础、多线程与锁、Red…...

分享115个HTML电子商务模板,总有一款适合您
分享115个HTML电子商务模板,总有一款适合您 115个HTML电子商务模板下载链接:https://pan.baidu.com/s/158y3jP0tv7ZikxNOBMKsSg?pwdt970 提取码:t970 Python采集代码下载链接:采集代码.zip - 蓝奏云 import os import shuti…...
Python 字符串
字符串是 Python 中最常用的数据类型。我们可以使用引号 ( 或 " ) 来创建字符串。创建字符串很简单,只要为变量分配一个值即可。例如:var1 Hello World!var2 "Python Runoob"Python 访问字符串中的值Python 不支持单字符类型ÿ…...
总线定义,车载总线:车载etherNet or CAN
总线的定义总线是连接多个设备或者接入点的数据传输通路。这里面的关键词是多个设备或者接入点,所以不要过于局限的看待总线。根据互联的设备/接入点不同,传输的数据带宽,速率,距离不同和应用场景的不同都可能有不同的总线。不同的…...

Python(for和while)循环嵌套及用法
Python 不仅支持 if 语句相互嵌套,while 和 for 循环结构也支持嵌套。 所谓嵌套(Nest),就是一条语句里面还有另一条语句,例如 for 里面还有 for,while 里面还有 while,甚至 while 中有 for 或者…...
进程地址空间(比特课总结)
一、进程地址空间 1. 环境变量 1 )⽤户级环境变量与系统级环境变量 全局属性:环境变量具有全局属性,会被⼦进程继承。例如当bash启动⼦进程时,环 境变量会⾃动传递给⼦进程。 本地变量限制:本地变量只在当前进程(ba…...
GitHub 趋势日报 (2025年06月08日)
📊 由 TrendForge 系统生成 | 🌐 https://trendforge.devlive.org/ 🌐 本日报中的项目描述已自动翻译为中文 📈 今日获星趋势图 今日获星趋势图 884 cognee 566 dify 414 HumanSystemOptimization 414 omni-tools 321 note-gen …...

毫米波雷达基础理论(3D+4D)
3D、4D毫米波雷达基础知识及厂商选型 PreView : https://mp.weixin.qq.com/s/bQkju4r6med7I3TBGJI_bQ 1. FMCW毫米波雷达基础知识 主要参考博文: 一文入门汽车毫米波雷达基本原理 :https://mp.weixin.qq.com/s/_EN7A5lKcz2Eh8dLnjE19w 毫米波雷达基础…...
flow_controllers
关键点: 流控制器类型: 同步(Sync):发布操作会阻塞,直到数据被确认发送。异步(Async):发布操作非阻塞,数据发送由后台线程处理。纯同步(PureSync…...

MeshGPT 笔记
[2311.15475] MeshGPT: Generating Triangle Meshes with Decoder-Only Transformers https://library.scholarcy.com/try 真正意义上的AI生成三维模型MESHGPT来袭!_哔哩哔哩_bilibili GitHub - lucidrains/meshgpt-pytorch: Implementation of MeshGPT, SOTA Me…...

如何做好一份技术文档?从规划到实践的完整指南
如何做好一份技术文档?从规划到实践的完整指南 🌟 嗨,我是IRpickstars! 🌌 总有一行代码,能点亮万千星辰。 🔍 在技术的宇宙中,我愿做永不停歇的探索者。 ✨ 用代码丈量世界&…...
shell脚本质数判断
shell脚本质数判断 shell输入一个正整数,判断是否为质数(素数)shell求1-100内的质数shell求给定数组输出其中的质数 shell输入一个正整数,判断是否为质数(素数) 思路: 1:1 2:1 2 3:1 2 3 4:1 2 3 4 5:1 2 3 4 5-------> 3:2 4:2 3 5:2 3…...
虚幻基础:角色旋转
能帮到你的话,就给个赞吧 😘 文章目录 移动组件使用控制器所需旋转:组件 使用 控制器旋转将旋转朝向运动:组件 使用 移动方向旋转 控制器旋转和移动旋转 缺点移动旋转:必须移动才能旋转,不移动不旋转控制器…...
StarRocks 全面向量化执行引擎深度解析
StarRocks 全面向量化执行引擎深度解析 StarRocks 的向量化执行引擎是其高性能的核心设计,相比传统行式处理引擎(如MySQL),性能可提升 5-10倍。以下是分层拆解: 1. 向量化 vs 传统行式处理 维度行式处理向量化处理数…...

2025-06-01-Hive 技术及应用介绍
Hive 技术及应用介绍 参考资料 Hive 技术原理Hive 架构及应用介绍Hive - 小海哥哥 de - 博客园https://cwiki.apache.org/confluence/display/Hive/Home(官方文档) Apache Hive 是基于 Hadoop 构建的数据仓库工具,它为海量结构化数据提供类 SQL 的查询能力…...