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 或者…...

C++_核心编程_多态案例二-制作饮品
#include <iostream> #include <string> using namespace std;/*制作饮品的大致流程为:煮水 - 冲泡 - 倒入杯中 - 加入辅料 利用多态技术实现本案例,提供抽象制作饮品基类,提供子类制作咖啡和茶叶*//*基类*/ class AbstractDr…...
rknn优化教程(二)
文章目录 1. 前述2. 三方库的封装2.1 xrepo中的库2.2 xrepo之外的库2.2.1 opencv2.2.2 rknnrt2.2.3 spdlog 3. rknn_engine库 1. 前述 OK,开始写第二篇的内容了。这篇博客主要能写一下: 如何给一些三方库按照xmake方式进行封装,供调用如何按…...

以下是对华为 HarmonyOS NETX 5属性动画(ArkTS)文档的结构化整理,通过层级标题、表格和代码块提升可读性:
一、属性动画概述NETX 作用:实现组件通用属性的渐变过渡效果,提升用户体验。支持属性:width、height、backgroundColor、opacity、scale、rotate、translate等。注意事项: 布局类属性(如宽高)变化时&#…...

【第二十一章 SDIO接口(SDIO)】
第二十一章 SDIO接口 目录 第二十一章 SDIO接口(SDIO) 1 SDIO 主要功能 2 SDIO 总线拓扑 3 SDIO 功能描述 3.1 SDIO 适配器 3.2 SDIOAHB 接口 4 卡功能描述 4.1 卡识别模式 4.2 卡复位 4.3 操作电压范围确认 4.4 卡识别过程 4.5 写数据块 4.6 读数据块 4.7 数据流…...

2021-03-15 iview一些问题
1.iview 在使用tree组件时,发现没有set类的方法,只有get,那么要改变tree值,只能遍历treeData,递归修改treeData的checked,发现无法更改,原因在于check模式下,子元素的勾选状态跟父节…...
在鸿蒙HarmonyOS 5中使用DevEco Studio实现录音机应用
1. 项目配置与权限设置 1.1 配置module.json5 {"module": {"requestPermissions": [{"name": "ohos.permission.MICROPHONE","reason": "录音需要麦克风权限"},{"name": "ohos.permission.WRITE…...
什么?连接服务器也能可视化显示界面?:基于X11 Forwarding + CentOS + MobaXterm实战指南
文章目录 什么是X11?环境准备实战步骤1️⃣ 服务器端配置(CentOS)2️⃣ 客户端配置(MobaXterm)3️⃣ 验证X11 Forwarding4️⃣ 运行自定义GUI程序(Python示例)5️⃣ 成功效果的原因分类及对应排查方案
JVM暂停(Stop-The-World,STW)的完整原因分类及对应排查方案,结合JVM运行机制和常见故障场景整理而成: 一、GC相关暂停 1. 安全点(Safepoint)阻塞 现象:JVM暂停但无GC日志,日志显示No GCs detected。原因:JVM等待所有线程进入安全点(如…...
MySQL账号权限管理指南:安全创建账户与精细授权技巧
在MySQL数据库管理中,合理创建用户账号并分配精确权限是保障数据安全的核心环节。直接使用root账号进行所有操作不仅危险且难以审计操作行为。今天我们来全面解析MySQL账号创建与权限分配的专业方法。 一、为何需要创建独立账号? 最小权限原则…...

QT开发技术【ffmpeg + QAudioOutput】音乐播放器
一、 介绍 使用ffmpeg 4.2.2 在数字化浪潮席卷全球的当下,音视频内容犹如璀璨繁星,点亮了人们的生活与工作。从短视频平台上令人捧腹的搞笑视频,到在线课堂中知识渊博的专家授课,再到影视平台上扣人心弦的高清大片,音…...