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

[特殊字符] 第十三讲 | 地统计模拟与空间不确定性评估

news 来源：原创 2025/9/10 11:01:39

📘 专栏：科研统计方法实战分享 | 地学/农学人的数据分析工具箱
✍️ 作者：平常心0715

🎯 关键词：地统计模拟、随机函数、空间不确定性、条件模拟、SGS、R语言

🧠 核心导语

在现实数据有限、空间异质性强的环境中，仅靠单一插值预测并不足以揭示背后的“不确定性”。此时，“地统计模拟”成为了空间研究不可或缺的工具！

💡 什么是地统计模拟？

与Kriging插值输出单一预测值不同，地统计模拟（Geostatistical Simulation）通过多次随机模拟，生成一系列可能的空间结果图，从而刻画不确定性，特别适合：

✅ 土壤污染风险图绘制
✅ 矿产资源储量估算
✅ 空间模型误差传播分析

🔍 模拟方式对比

方法特点
SGS（Sequential Gaussian Simulation）连续变量建模首选，模拟值更平滑
SIS（Sequential Indicator Simulation）离散/分类型变量建模如土地利用等
Turning Bands高频模拟，适合大数据场景

📦 R语言实战简例（以 gstat 包为例）

library(gstat)
library(sp)# 假设已有空间点数据 `dat`
coordinates(dat) <- ~x + y# 拟合变异函数
vgm_mod <- variogram(value ~ 1, dat)
vgm_fit <- fit.variogram(vgm_mod, vgm("Sph"))# 生成模拟网格
grd <- spsample(dat, n = 500, type = "regular")
gridded(grd) <- TRUE# 进行SGS模拟（100 realizations）
sim_result <- krige(value ~ 1, dat, grd, model = vgm_fit, nmax = 20, nsim = 100)

📊 不确定性评估

你可以通过这些方式评估模拟结果：

  • 平均图（Mean Map）：各位置的平均预测值；

  • 标准差图（Std. Dev. Map）：反映空间预测的稳定性；

  • 95%置信区间图：量化置信程度。

✍️ 科研写作句式

“通过顺序高斯模拟方法构建100组空间预测结果，进一步量化研究区土壤污染分布的不确定性，结果显示热点区域的不确定性较高，需进一步实地调查。”

🧩 小结

✅ 插值只给一个答案，模拟提供多个“可能”；
✅ SGS 是连续型变量建模利器；
✅ 不确定性分析帮助我们对空间结果有更科学判断！

📌 下一讲预告：空间异质性检验与地统计局部指标（LISA）应用

如果你觉得这篇文章有帮助，欢迎点赞 + 收藏！
我也会继续带来更多可直接实操的地理/农学统计方法系列分享~

相关文章：

[特殊字符] 第十三讲 | 地统计模拟与空间不确定性评估

&#x1f4d8; 专栏&#xff1a;科研统计方法实战分享 | 地学/农学人的数据分析工具箱 ✍️ 作者&#xff1a;平常心0715 &#x1f3af; 关键词&#xff1a;地统计模拟、随机函数、空间不确定性、条件模拟、SGS、R语言 &#x1f9e0; 核心导语 在现实数据有限、空间异质性强的…...

编程日记 2025/9/10 11:01:39

Vue接口平台学习六——接口列表及部分调试页面

一、实现效果图及界面布局简单梳理 整体布局分左右&#xff0c;左边调试&#xff0c;右边显示接口列表 左侧&#xff1a; 一个输入框按钮&#xff1b;下面展示信息&#xff0c;大部分使用代码编辑器就好了&#xff0c;除了请求体传文件类型需要额外处理。然后再下方显示响应信…...

编程日记 2025/9/11 8:03:00

Spring 中的 @Cacheable 缓存注解

1 什么是缓存 第一个问题&#xff0c;首先要搞明白什么是缓存&#xff0c;缓存的意义是什么。 对于普通业务&#xff0c;如果要查询一个数据&#xff0c;一般直接select数据库进行查找。但是在高流量的情况下&#xff0c;直接查找数据库就会成为性能的瓶颈。因为数据库查找的…...

编程日记 2025/9/11 8:04:56

Context的全面解析：在不同技术应用中的通用作用与差异

Context的全面解析&#xff1a;在不同技术应用中的通用作用与差异 引言&#xff1a; 在软件开发中&#xff0c;“Context”这个概念被广泛使用。它不仅限于某个特定的技术或编程语言&#xff0c;实际上&#xff0c;Context 作为一种抽象的设计模式&#xff0c;贯穿在许多开发领…...

编程日记 2025/9/11 8:05:30

机器学习（2）——逻辑回归

文章目录 1. 什么是逻辑回归?2. 核心思想3. 逻辑回归模型的训练&#xff1a;4. 参数估计&#xff08;损失函数与优化&#xff09;4.1. **损失函数&#xff1a;**4.2. 极大似然估计&#xff08;MLE&#xff09;4.3. 优化方法 5. 决策边界6. 模型评估指标7 . 假设与适用条件8. 逻…...

编程日记 2025/9/11 8:01:40

Sentinel核心算法解析の滑动窗口算法

文章目录 前言一、回顾&#xff1a;快速失败二、固定窗口算法三、滑动窗口算法三、源码体现3.1、ArrayMetric的初始化3.2、addPass3.2.1、currentWindow3.2.2、wrap.value().addPass 总结 前言 在Sentinel中&#xff0c;流控效果有快速失败、预热和排队等待。其中快速失败的统计…...

编程日记 2025/9/7 7:36:51

ida 使用记录

文章目录 伪代码-汇编hexstring快捷键 伪代码-汇编 流程图界面——F5——伪代码界面——再点Tab——流程图界面——再按空格——汇编界面流程图界面——空格——汇编界面 hex view - open subviews - hex dump string view - open subviews - string快捷键&#xff1a; sh…...

编程日记 2025/9/10 9:53:43

数字统计：

1.题意&#xff1a; 在1~N之间寻找d出现的个数&#xff0c;然后输出即可&#xff1b;例如&#xff1a;d2,N23&#xff0c;那么满足条件的有2,12,21,23&#xff0c;所以是4个 2.思路&#xff1a; 1.暴力枚举&#xff08;不可能&#xff09;&#xff1a;可以先写出来去找规律 …...

编程日记 2025/9/11 8:02:58

【架构师从入门到进阶】第五章：DNSCDN网关优化思路——第八节：网关-注入攻击与预防

【架构师从入门到进阶】第五章&#xff1a;DNS&CDN&网关优化思路——第八节&#xff1a;网关-注入攻击与预防 SQL注入攻击的原理攻击者获取数据库表结构预防SQL注入的方法 这篇文章我们来看SQL注入。 SQL注入攻击的原理 SQL注入攻击的原理呢&#xff1f;我们来简单说…...

编程日记 2025/9/10 13:06:59

波束形成（BF）从算法仿真到工程源码实现-第五节-线性约束最小方差波束形成算法(LCMV)

一、概述 本节我们讨论线性约束最小方差波束形成算法(Linearly constrained minimum variance,LCMV)波束形成算法&#xff0c;包括原理分析及代码实现。 更多资料和代码可以进入https://t.zsxq.com/qgmoN &#xff0c;同时欢迎大家提出宝贵的建议&#xff0c;以共同探讨学习。 …...

编程日记 2025/8/31 21:03:16

Java类加载机制原理与应用

前言 Java 中的类加载机制&#xff08;Class Loading Mechanism&#xff09;是 JVM 架构中的核心组成部分&#xff0c;它控制着类从编译后的 .class 文件被加载到内存、并最终变成可以被程序使用的对象的全过程。涉及类加载器、双亲委派模型及加载过程。下面我们从原理到实际应…...

编程日记 2025/9/7 15:13:48

android display 笔记（十三）surfcaeflinger的DEQUEUED、QUEUED

BufferQueue 的核心作用 BufferQueue 是 生产者-消费者模型 的核心组件&#xff0c;协调应用&#xff08;生产者&#xff09;和 SurfaceFlinger&#xff08;消费者&#xff09;之间的图形缓冲区&#xff08;GraphicBuffer&#xff09;传递。 生产者&#xff1a;应用&#xff0…...

编程日记 2025/9/10 1:25:51

数据库预热

介绍 Database Warm-up &#x1f9e0; 一句话理解 数据库是在应用启动阶段&#xff0c;提前建立数据库连接 或 执行轻量 SQL 操作&#xff0c;从而 加快首个请求的响应速度 的一种优化手段 &#x1f3af; 为什么需要数据库预热&#xff1f; 当 FastAPI 或其他 Web 服务刚启…...

编程日记 2025/9/4 18:01:50

C语言—程序的编译和链接

1. 翻译环境和运行环境 在ANSI S的任何一种实现中&#xff0c;存在两个不同的环境 第一种是翻译环境&#xff0c;在这个环境中源代码被转换为可执行的机器指令&#xff08;二进制指令&#xff09; 第二种是执行环境&#xff0c;它用于实际执行代码 2. 翻译环境 翻译环境是由…...

编程日记 2025/9/9 3:38:55

Neo4j GDS-10-neo4j GDS 库中相似度算法介绍

neo4j apoc 系列 Neo4j APOC-01-图数据库 apoc 插件介绍 Neo4j GDS-01-graph-data-science 图数据科学插件库概览 Neo4j GDS-02-graph-data-science 插件库安装实战笔记 Neo4j GDS-03-graph-data-science 简单聊一聊图数据科学插件库 Neo4j GDS-04-图的中心性分析介绍 Ne…...

编程日记 2025/9/6 4:54:52

Unity 动画

Apply Root Motion 勾选的话就会使用动画片段自带的位移 Update Mode &#xff08;动画重新计算骨骼位置转向缩放的数值&#xff09;&#xff1a; Normal &#xff1a; 随Update走&#xff0c;每次Update都计算Animate Physics &#xff1a;与 fixed Update() 同步&#xff0…...

编程日记 2025/9/10 11:31:55

【位运算】只出现一次的数字 II

文章目录 137. 只出现一次的数字 II解题思路一&#xff1a;借用数组的位运算解法二&#xff1a;不使用数组的位运算 137. 只出现一次的数字 II 137. 只出现一次的数字 II ​ 给你一个整数数组 nums &#xff0c;除某个元素仅出现 一次 外&#xff0c;其余每个元素都恰出现 **…...

编程日记 2025/9/6 15:21:12

模型开发中的微调是干什么

在模型开发中&#xff0c;微调&#xff08;Fine-tuning&#xff09; 是指利用预训练模型&#xff08;Pre-trained Model&#xff09;的参数作为初始值&#xff0c;在特定任务或数据集上进一步调整模型参数的过程。它是迁移学习&#xff08;Transfer Learning&#xff09;的核心…...

编程日记 2025/9/7 16:09:17

leetcode 204. Count Primes

题目描述 这是道纯数学类问题。 先回忆一下&#xff0c;素数的定义。 质数&#xff08;英文名&#xff1a;Prime number&#xff09;又称素数&#xff0c;是指在大于1的自然数中&#xff0c;除了1和它本身以外不再有其他因数的自然数。 质数又称素数。一个大于1的自然数&…...

编程日记 2025/9/7 11:45:16

fastadmin后端添加页面,自主控制弹出框关闭,关闭父页面弹框

Form.api.bindevent($(“form[roleform]”), (data, ret) > { 重写绑定事件,返回false即可 注意:只有返回code1才能拦截,其他值不进行拦截 add: function () {//获取当前search里面的type值var type location.search.split(type)[1];Form.api.bindevent($("form[role…...

编程日记 2025/9/2 12:10:15

LeetCode 255 超通俗讲解：Swift 验证前序是否 BST

文章目录 摘要描述题解答案题解代码分析核心点解释&#xff1a; 示例测试及结果时间复杂度空间复杂度总结未来展望 摘要 在做算法题的时候&#xff0c;树相关的题总是“神神叨叨”的&#xff0c;但其实抓住核心规则&#xff0c;它们也挺有逻辑的。今天这题——LeetCode 255&am…...

编程日记 2025/9/10 20:39:48

Win32++ 使用初探

文章目录 1. 环境要求2. Win32安装3. 项目创建3.1 项目创建&#xff08;1&#xff09;直接使用Win32里的示例Sample&#xff08;2&#xff09;自行创建项目 最近想用 VC写些 UI&#xff0c;但又不太想用 MFC&#xff0c;正好对界面要求不太高&#xff0c;就使用了一下 Win3…...

编程日记 2025/8/29 9:05:51

求解时间复杂度

1.设 t 法 当求解出现while循环时&#xff0c;设t求解 void fun(int n) {int i 1;while(i < n)i i * 2; } 解法&#xff1a; 1.设循环次数为t&#xff1b; 2.将while循环中的语句展开到循环t次 1 2 3 …… t 2 2^2 2^3 …… 2^t 3.跳出循环 2^t > n …...

编程日记 2025/9/10 0:27:24

深度解析：如何高效识别并定位问题关键词

什么是问题关键词&#xff1f; 问题关键词是人们在搜索引擎中输入以查找信息、答案或解决方案的问题。这些查询以问题指示符开头&#xff0c;例如&#xff1a; who、what、where、when、why、how、which、will、would、should、can、could、is、are、was、were、do、does 或 d…...

编程日记 2025/9/8 13:15:38

c++小做——完全数

今天&#xff0c;我们来写一个完全数的代码 首先是 long long n; cin>>n; &#xff08;you~输入的数&#xff09; 然后是 long long b0;//因数的和 long long cnt0;//计数器 接着是 for(long long i2;i<n-1;i) {} 在里面插入 bb-i;再写一个for for(int a1;a&…...

编程日记 2025/9/5 15:10:36

GGML源码逐行调试（下）

目录 前言1. 简述2. 预分配计算图内存2.1 创建图内存分配器2.2 构建最坏情况的计算图2.3 预留计算图内存 3. 分词4. 模型推理与生成4.1 模型推理4.2 采样 结语下载链接参考 前言 学习 UP 主 比飞鸟贵重的多_HKL 的 GGML源码逐行调试 视频&#xff0c;记录下个人学习笔记&#x…...

编程日记 2025/9/8 1:28:13

JavaScript学习教程，从入门到精通， JavaScript 函数全面解析与案例实践（11）

JavaScript 函数全面解析与案例实践 项目导读 JavaScript 函数是编程中的核心概念&#xff0c;是执行特定任务的代码块。本教程将全面讲解函数的定义、参数、返回值及调用方式&#xff0c;并通过实际案例加深理解。 学习目标 掌握 JavaScript 函数的定义与调用方法理解函数…...

编程日记 2025/9/10 10:37:00

音视频之H.265/HEVC编码框架及编码视频格式

一、编码框架&#xff1a; H.265/HEVC采用混合编码框架&#xff0c;包括变换、量化、熵编码、帧内预测、帧预测以及环路滤波等模块。但是&#xff0c;H.265/HEVC几乎在每个模块都引入了新的编码技术。 1、帧内预测&#xff1a; 该模块主要用于去除图像的空间相关性。通过编码后…...

编程日记 2025/9/8 18:25:55

栈与队列：两种经典线性数据结构的深度解析

一、栈&#xff1a;LIFO 特性的完美诠释 &#xff08;一&#xff09;核心概念与抽象模型 定义与特性 栈是一种严格遵循后进先出&#xff08;LIFO&#xff09;原则的线性数据结构&#xff0c;其操作被限制在栈顶&#xff08;Top&#xff09;进行。形象化理解&#xff1a;如同堆…...

编程日记 2025/9/5 20:28:02

0x01、Redis 主从复制的实现原理是什么？

Redis 主从复制概述 Redis 的主从复制是一种机制&#xff0c;允许一个主节点&#xff08;主实例&#xff09;将数据复制到一个或多个从节点&#xff08;从实例&#xff09;。通过这一机制&#xff0c;从节点可以获取主节点的数据并与之保持同步。 复制流程 开始同步&#xf…...

编程日记 2025/9/3 12:05:56

Python实现贪吃蛇一

贪吃蛇是一款经典的小游戏&#xff0c;最近尝试用Python实现它。先做一个基础版本实现以下目标&#xff1a; 1、做一个按钮&#xff0c;控制游戏开始 2、按Q键退出游戏 3、右上角显示一个记分牌 4、随机生成一个食物&#xff0c;蛇吃到食物后长度加一&#xff0c;得10分 5、蛇碰…...

编程日记 2025/9/10 9:48:20

01-libVLC的视频播放器：环境搭建以及介绍

项目展示项目播放器 VLC简介VLC媒体播放器（VideoLAN Client）是一款开源、跨平台的自由多媒体播放器，由VideoLAN项目开发。它支持众多音频与视频格式（如MPEG-2、MPEG-4、H.264、MKV、WebM、WMV、MP3等），以及DVD、VCD和各种流媒体协议。 VLC的特点跨平台支持：Windows、mac…...

编程日记 2025/9/10 18:30:27

linux内核升级

这里介绍一下linux内核升级 因为需要搭建k8s集群内核内核版本过低会导致集群出现问题&#xff0c;为了避免问题发生我们对集群内核进行升级处理 这个是我目前本身的内核版本 用了很多的镜像站去进行更新发现更新不了&#xff08;阿里云不能用了&#xff0c;貌似是删除了&…...

编程日记 2025/9/9 19:33:53

电感详解：定义、作用、分类与使用要点

一、电感的基本定义 电感&#xff08;Inductor&#xff09; 是由导线绕制而成的储能元件&#xff0c;其核心特性是阻碍电流变化&#xff0c;将电能转化为磁能存储。 基本公式&#xff1a; 自感电动势&#xff1a; E -L * (di/dt) &#xff08;L&#xff1a;电感值&#xff0c…...

编程日记 2025/9/8 6:51:51

扩散模型简介

扩散模型简介 基本原理 扩散模型是一种基于概率扩散过程的生成模型&#xff0c;其核心思想是通过正向扩散过程和反向去噪过程生成数据&#xff1a; 正向扩散过程&#xff1a;从真实数据&#xff08;如图像&#xff09;开始&#xff0c;逐步添加高斯噪声&#xff0c;最终将数据…...

编程日记 2025/9/8 6:33:31

MySQL安装实战分享

一、在 Windows 上安装 MySQL 1. 下载 MySQL 安装包 访问 MySQL 官方下载页面。选择适合你操作系统的版本。一般推荐下载 MySQL Installer。 2. 运行安装程序 双击下载的安装文件&#xff08;例如 mysql-installer-community-<version>.msi&#xff09;。如果出现安全…...

编程日记 2025/9/7 1:24:06

掌握 Git 的十大基础命令

李升伟 编译 在 IT 领域&#xff0c;很少有技术能像 Git 一样占据绝对主导地位&#xff0c;几乎无人能及。Git 在软件开发中扮演着核心角色&#xff0c;其影响力之大甚至让其他版本控制系统&#xff08;如 SVN 和 Mercurial&#xff09;几乎被淘汰。如今&#xff0c;我们已难以…...

编程日记 2025/9/10 3:17:41

58-使用wordpress快速创建个人网站

直接找台可以联网的linux&#xff08;我的环境是rocky8.9&#xff09;一顿运行&#xff0c;思路就是安装docker&#xff0c;然后启动一个数据库&#xff0c;然后启动一个wordpress&#xff0c;然后就是把端口暴露出来。 227 yum remove podman 228 yum install -y yum-utils…...

编程日记 2025/9/3 8:24:23

若依前后端分离版运行教程、打包教程、部署教程

后端打包教程 注意&#xff1a;需要先运行redis 2、前端运行教程 2.1安装依赖 2.2运行 打开浏览器查看,地址&#xff1a;http://localhost:80 3、前端打包教程 3.1打包 3.2运行打包好的文件&#xff0c;先找到打包好的文件 这是nginx的文件结构 将打包好的文件放到html目录下…...

编程日记 2025/9/9 14:57:34

【Python3教程】Python3基础篇之数据结构

博主介绍：✌全网粉丝22W+，CSDN博客专家、Java领域优质创作者，掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域✌ 技术范围：SpringBoot、SpringCloud、Vue、SSM、HTML、Nodejs、Python、MySQL、PostgreSQL、大数据、物联网、机器学习等设计与开发。 感兴趣的可…...

编程日记 2025/9/9 21:15:33

transformers的 pipeline是什么:将模型加载、数据预处理、推理等步骤进行了封装

transformers的 pipeline是什么:将模型加载、数据预处理、推理等步骤进行了封装 pipe = pipeline("text-generation", model=model, tokenizer=tokenizer, max_new_tokens=50 )pipeline ：这是 transformers 库中一个非常实用的工具函数。它可以基于预训练模型快速构…...

编程日记 2025/8/17 19:17:13

十七、TCP编程

TCP 编程是网络通信的核心&#xff0c;其 API 围绕面向连接的特性设计&#xff0c;涵盖服务端和客户端的交互流程。以下是基于 ​C 语言的 TCP 编程核心 API 及使用流程的详细解析&#xff1a; 核心 API 概览 ​函数​角色​描述socket()通用创建套接字&#xff0c;指定协议族…...

编程日记 2025/9/10 0:15:25

Obsidian 技巧篇

Obsidian 技巧篇 本篇文章主要汇总分享几个 Ob 中好用的小技巧&#xff0c;包括嵌入视频播放、文本颜色设置、插入大纲、Mermaid 绘制图形。原文见于&#xff1a;Obsidian技巧篇。 嵌入视频播放 <iframe width"860" height"700" src"https://ww…...

编程日记 2025/9/5 19:30:42

使用Fortran读取HDF5数据

使用Fortran读取HDF5数据 下面我将介绍如何在Fortran中读取HDF5文件中的各种类型数组数据&#xff0c;包括一维数组、二维数组、元数组和变长数组。 准备工作 首先需要确保系统安装了HDF5库&#xff0c;并且在编译时链接了HDF5库。例如使用gfortran编译时&#xff1a; gfor…...

编程日记 2025/9/7 0:46:11

L36.【LeetCode题解】查找总价格为目标值的两个商品(剑指offer:和为s的两个数字) (双指针思想,内含详细的优化过程)

目录 1.LeetCode题目 2.分析 方法1:暴力枚举(未优化的双指针) 方法2:双指针优化:利用有序数组的单调性 版本1代码 提问:版本1代码有可以优化的空间吗? 版本2代码 提问:版本2代码有可以优化的空间吗? 版本3代码(★推荐★) 3.牛客网题目:和为s的数字 1.LeetCode题目 …...

编程日记 2025/9/4 19:18:16

mysql 商城商品属性开发的动态解决方案

终极方案&#xff1a;动态属性解决方案 推荐使用 JSON 字段 虚拟列索引 的组合方案 结合灵活存储与查询优化&#xff0c;平衡扩展性与性能 完整实现步骤 步骤 1&#xff1a;创建基础表结构 CREATE TABLE products (id INT PRIMARY KEY AUTO_INCREMENT,name VARCHAR(100) NO…...

编程日记 2025/8/31 10:59:39

Java递归练习----猴子偷桃

问题&#xff1a; 有一堆桃子&#xff0c;猴子第一天吃灵其中的一般&#xff0c;并在多吃了一个&#xff01;以后每天猴子都吃其中的一半&#xff0c;然后多吃一个。当到第十天时&#xff0c;想再吃时&#xff08;即还没吃&#xff09;&#xff0c;发现只有1个桃子了&#xff…...

编程日记 2025/9/9 4:08:59

[干货]PHM学习软件|PHM预测性维护系统

使用步骤教程如下 1、登录 用户名&#xff1a;52phm 密码&#xff1a;xxx &#xff08;区别在于不同用户密钥不一样&#xff09; 2、上传需要分析的数据集 支持数据集格式&#xff1a;csv、xlsx、xls、mat、json 3、主题1&#xff1a;机械参数计算 计算轴承、齿轮、皮带的…...

编程日记 2025/9/8 5:08:00

详解正则表达式中的?:、?= 、 ?! 、?＜=、?＜!

1、?: - 非捕获组 语法: (?:pattern) 作用: 创建一个分组但不捕获匹配结果&#xff0c;不会将匹配的文本存储到内存中供后续使用。 优势: 提高性能和效率 不占用编号&#xff08;不会影响后续捕获组的编号&#xff09; 减少内存使用 // 使用捕获组 let regex1 /(hell…...

编程日记 2025/9/4 13:50:45

Java常见面试问题

一.Liunx 二.Java基础 1.final 2.static 3.与equals 三.Collection 1.LIst 2.Map 3.Stream 四、多线程 1.实现方法 2.线程池核心参数 3.应用场景 五、JVM 1.堆 2.栈 六、Spring 1.面向对象 2.IOC 3.AOP 七、Springboot 1.自动装配 八、SpringCloud 1.Nacos 2.seata 3.ga…...

编程日记 2025/9/10 1:25:53