基于python的站点数据Kriging插值绘图

相关文章推荐

虚心的薯片 · 臺灣大學資訊系統訓練班· 1 周前 ·

好帅的野马 · nodejs调用脚本(python/shel ...· 1 周前 ·

纯真的丝瓜 · 使用aiohttp实现Python异步HTT ...· 3 天前 ·

乖乖的橙子 · python利用requests库模拟pos ...· 3 天前 ·

失落的木瓜 · Python报错：PermissionErr ...· 12 小时前 ·

聪明的猕猴桃 · 获取GCS blob的文件名· 1 年前 ·

长情的回锅肉 · mysql错误代码1304-掘金· 1 年前 ·

自信的橡皮擦 · 在excel里处理数据的VBA代码 - 哔哩哔哩· 2 年前 ·

谈吐大方的伏特加 · Python学习，使用Python实现多表分 ...· 2 年前 ·

鬼畜的面包 · Go-Excelize ...· 2 年前 ·

科研中常常会将站点数据进行插值，绘制成图。常用的二维插值方法有 最近邻法、线性内插法、三次样条内插法 ，此外还有一些基于地理的插值方法，如克里金插值法、IDW反距离加权法。今天我们就克里金插值法介绍一下使用python进行站点数据插值绘图的方法。

绘图模块


    cartopy

、


    shpfile

、


    matplotlib

插值模块：

对于简单的二维插值方法（最近邻法、线性内插法、三次样条内插法），可以用scipy.interpolate中的griddata方法实现。


    scipy.interpolate.griddata(points, values, xi, method='linear')

对于克里金/IDW插值方法，可以自行编写计算函数，也可以调用现成的模块，如 pykrige 。

指北针的添加

导入 import matplotlib.patches as mpatches # 比例尺

加载地图、矢量数据

kriging插值

其他（ colorbar/contour/刻度标签 ...）

绘图用 cartopy ，插值用 scipy / pykrige ，掩膜用 shp2clip 。

站点插值 绘图与利用 tiff绘图 稍有不同，主要就区别在插值这块，tiff绘图是 直接读取tiif边界坐标，插值成xsize/ysize个点位 。

常用shp需要收集，如中国城市、全国国家矢量数据...

推荐文章

虚心的薯片 · 臺灣大學資訊系統訓練班

1 周前

好帅的野马 · nodejs调用脚本(python/shell)和系统命令 - 单曲荨环

1 周前

纯真的丝瓜 · 使用aiohttp实现Python异步HTTP POST请求_python 异步发post请求,不等返回结果的异步

3 天前

乖乖的橙子 · python利用requests库模拟post请求时json的使用 - NewJune

3 天前

失落的木瓜 · Python报错：PermissionError: [Errno 13] Permission denied解决方案详解

12 小时前

聪明的猕猴桃 · 获取GCS blob的文件名

1 年前

长情的回锅肉 · mysql错误代码1304-掘金

1 年前

自信的橡皮擦 · 在excel里处理数据的VBA代码 - 哔哩哔哩

2 年前

谈吐大方的伏特加 · Python学习，使用Python实现多表分批次从数据库导出到Excel_python如何把excel分批输出_sehun_sx的博客-CSDN博客

2 年前

鬼畜的面包 · Go-Excelize API源码阅读（四）——Save()、SaveAs(name string) - 掘金

2 年前