![](https://img-blog.csdnimg.cn/20201014180756922.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
python
文章平均质量分 59
地学万事屋
读博中......
展开
-
Python Matlab R的Mann-Kendall趋势检验
Python Matlab R的Mann-Kendall趋势检验水文气象中推荐使用Mann-Kendall趋势检验这是一种非参数统计检验方法,在中心趋势不稳定时,关注数据的秩。该方法不需要不需要满足正态分布的假设,因而具有普适性。根据自己需要(图像、并行计算、线趋势图等等)分享python\matlab\R的方法Python进行Mann-Kendall趋势检验代码如下:# -*- coding: utf-8 -*-from __future__原创 2023-02-04 10:54:24 · 2405 阅读 · 1 评论 -
matplotlib你真的明白plt,fig和ax吗
plt系列接口是用来模仿MATLAB的风格的。如果用户熟悉MATLAB,还是建议以plt为主的。对于基本的画图,plt系列足够了。个人觉得,真正不合适的用法,反而是各种混淆plt接口和其它更细节的内部接口函数,混杂着用,成了四不像。和和这就表明如果每次都要手动搜索的话,会非常不确定,因为一种效果有完全不同的实现方式,这次搜到这个,下次搜到那个:笔者在初学的时候,碰见在figure, axes, axis这几个词就头大。原创 2022-11-28 15:54:42 · 1253 阅读 · 1 评论 -
python xarray批量处理nc
xarray是一类非常灵活的结构,被看作新一代nc处理库原创 2022-11-28 15:53:55 · 2936 阅读 · 0 评论 -
jupyterlab教程
总结:Jupyter Lab作为Notebook的上位替代,是未来数据科学的趋势。然而目前的版本还不完善,可以尝鲜使用。4.Jupyter Lab支持两种背景主题,白色和黑色。步骤:settings->jupyterlab theme。还能多窗口操作文件,如一边看文档一边写代码。2.输入以下命令进入jupyter lab。这玩意还是不行,还得MikTex曲线救国。1.安装jupyterlab。8.测试geemap。原创 2022-11-25 15:42:50 · 696 阅读 · 0 评论 -
python和R绘制箱线图+抖动点
通过散点展示数据的直观分布x轴刻度标签如果是“连续”数据,可添加渐变背景。通过渐变来体现升温幅度,美观形象。原创 2022-11-25 15:19:08 · 970 阅读 · 0 评论 -
用python gdal裁剪栅格数据提取添加xy经纬度和栅格值
alfalfa是一种亚洲西南部多年生草本植物,是重要的经济作物。在图中也可以看到,主要分布在热带和南美洲。对于我上述全球影像来说,栅格转点需要6分钟。添加字段和计算几何都需要花费更多的时间。这里我设置nodata为负值,是我本来影像的nodata值,可以在GIS查看。采用python的gdal方法,首先进行影像裁剪。现有一全球经济作物数据alfalfa的产量。上述功能在ArcGIS中是这样实现的。问题:把遥感影像转为一张表。我们想把影像转表,即。原创 2022-11-20 17:20:32 · 1632 阅读 · 0 评论 -
python gdal geopandas basemap cartopy安装
Basemap cartopy geopandas rasterio这几个库存在一定的依赖关系,由于环境的变化很可能哪天不知道就报错了:版本不一致,运行也可能报错,如GDAL3.0的主要变化在于对空间参考系统进行了大规模的重构,在使用上接口发生了变化。首先检查哪些库出了问题。原创 2022-11-20 11:08:42 · 954 阅读 · 1 评论 -
教你自己写Arcpy批处理程序
先上代码,讲解各行代码的意思,从而达到自己写代码的目的为方便说明,逐行讲解:本代码实现【把栅格值提取到点】:输入元素是栅格影像 + 点矢量,输出带像元值属性的矢量表本代码实现【表转excel】:把输出的矢量点属性表转为xls接下来对影像存储的文件夹进行遍历,保证文件夹中都是tif因此,若要批量对数据进行处理,只需修改核心函数就可以了,其次再修改文件名、文件夹、路径和切片就完成了。寻找函数的话,只需要在你想要的工具下面,点击Tool Help即可提供了源代码与数据讲解,我们只需要把核心函数复制过来就行了有了原创 2022-11-19 11:24:04 · 632 阅读 · 1 评论 -
Arcpy批量裁剪mask影像tif
3是输出名字,改变从第4位到第9位,来区别文件(如下图VCF5KYR_2013001_001_201824210156),输出名字截取从第8位到第12位,来区别文件(如下图VCF5KYR_2013001_001_201824210156.tif ,2013是第8-12位)6可以改为任意ArcGIS操作,如near,buffer等。1是需要被裁剪影像的路径。该路径需要存放所有TIF文件。4是输出路径,重新命名为VCF_2013china。2与1完全相同,如下图。5是矢量文件所在位置。原创 2022-11-19 11:16:41 · 411 阅读 · 0 评论 -
ArcGIS pro制图添加经纬网图例
右键 - 属性更改刻度、标签、刻线等等。原创 2022-11-19 10:29:57 · 5939 阅读 · 1 评论 -
解决jupyter TOC勾选了但不显示的问题
有时候TOC(Table of content2)反应很慢,或者勾选了根本就不显示。或者隔三岔五nbextension消失,按以下步骤解决问题:fill:#333;color:#333;color:#333;fill:none;}成功失败重装nbxetension完成降低Jupyter版本。原创 2022-11-18 12:37:56 · 1077 阅读 · 0 评论 -
更新后的python seaborn用于科研图绘
数据科学和科研绘图的人离不开R的ggplot2,尽管python已经有了诸如matplotlib、seaborn等强大绘图包,兼有basemap、cartopy等空间数据可视化能力,但相较于R的ggplot,总感觉欠缺了点东西。此外,整个seaborn的API重写后,操作方式也变化很大,科研人表示要重新开始学习了,不得不感叹开源软件发展的速度,稍微不学习就会落下。最后是分面的功能,这在R中非常容易,现在在python中也很方便了。先以散点图为例,这个add()就像R中的 + geom()…原创 2022-11-18 12:36:34 · 530 阅读 · 0 评论 -
python folium 实现地图平台制作
直接上代码,这次是依赖于python的folium的地图包,folium是基于leaflet的经典js地图库,在python中重写使之能适用于python。原创 2022-11-18 12:34:39 · 593 阅读 · 0 评论 -
python批量读取nc气象数据并转为tif
各类地理数据中,NC格式是很常见的,然而这种格式ArcGIS是无法打开的。一般nc数据中lon和lat是必有的,只不过名字也许不同,可能有的是“longtitude”和“latitude”,我们要的数据是“di_01”。由于不同NC维度不同,数据格式、类型、存储方式都可能不同,因此一定要真正理解代码,才能对其进行修改,使之适合自己的数据。最后是输入路径和输出路径,输出路径是一空文件夹,输入路径记录了全部的nc文件。下图中是结果,红框是变量名,蓝框是数据类型,绿框是数据大小。径记录了全部的nc文件。原创 2022-11-18 12:33:15 · 2504 阅读 · 6 评论 -
python绘制中国降雨分布(添加等高线、比例尺、指南针、南海诸岛、九段线)
画一个比例尺带'N'文字注释主要参数如下:param ax: 要画的坐标区域 Axes实例 plt.gca()获取即可:param labelsize: 显示'N'文字的大小:param loc_x: 以文字下部为中心的占整个ax横向比例:param loc_y: 以文字下部为中心的占整个ax纵向比例:param width: 指南针占ax比例宽度:param height: 指南针占ax比例高度:param pad: 文字符号占ax比例间隙。原创 2022-11-18 12:32:32 · 2723 阅读 · 0 评论 -
python绘制论文常见的散点密度图
加载数据源这里默认是csv或excel的一列,一般数据处理好都是这样,一列Observe,一列Predict。如果数据源不一致(如txt)存进来是list,也强烈推荐存成DataFrame。这里是线性拟合(最小二乘)有k和b,没有计算p值,若需要还要自己算一下。这里是线性拟合(最小二乘)有k和b,没有计算p值,若需要还要自己算一下。经常在论文中见到validation的图,那么这种图应该怎么做呢?RMSE等等是自己写的程序(函数)实现。RMSE等等是自己写的程序(函数)实现。原创 2022-11-18 12:31:50 · 960 阅读 · 0 评论 -
python绘制Bubble气泡图pltscatter
float 或 array,维度为(n, ), 非必需。) 是标记的形状,具体可看matplotlib官方文档。如果没有,则使用颜色数组的相应最小值和最大值。) 颜色条,可以用colormap,参考官方文档。float或 array 维度为(n, )**c **array或list,映射颜色。0 (完全透明) and 1 (不透明).与默认规范一起使用,以将颜色数组。如果为“真”,则使用。.缩放颜色的选项,一般不用。不能与norm同用。原创 2022-11-18 12:29:03 · 2974 阅读 · 0 评论 -
python绘图使用matplotlib色卡
Matplotlib更多的色卡可见:https://matplotlib.org/stable/tutorials/colors/colormaps.html。seaborn在此基础上封装的更厉害,直接使用palette=colorname即可。python绘图一般使用两个库,原创 2022-11-18 12:28:12 · 3030 阅读 · 0 评论 -
python读取降雨格点数据
中国地面降水月值0.5°×0.5°格点数据集文件命名由数据集代码(SURF_CLI_CHN_PRE_MON_GRID_0.5)、年份、月份标识(YYYYMM)组成。中国地面降水月值0.5°×0.5°格点数据集(V2.0)196101-201305。这个数据提供了1961-2013年的长时期站点观测,由气象站插值得到,质量较高。时间分辨率:逐月 3.4.空间属性 3.4.1.地理范围 地理范围描述:中国。接下来求所有降水数据平均值,由于我想求年总降水,数据是月均值。这个数据我也免费提供给大家,见文末。原创 2022-11-18 12:24:08 · 1220 阅读 · 7 评论 -
python 绘制手绘风格图
外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-TKIaA616-1668745282043)(C:/Users/wang/AppData/Roaming/Typora/typora-user-images/image-20220827235703718.png)]这个Matplotlib子函数特别简单,只有三个参数,别看参数少,但功能可不小。原创 2022-11-18 12:22:51 · 283 阅读 · 0 评论 -
cartopy绘制中国降雨地图
Cartopy 是英国气象局开发的地图绘图包,实现了 Basemap 的大部分功能,还可以通过 Matplotlib 的 API 实现丰富的自定义效果。相比于ArcGIS、QGIS和ArcGISpro用鼠标点来点去,编程绘图也是有很大的优点的,方便,可批量,美观。用于添加栅格,x和y对应经纬度,一般来说nc文件自带这个数据,prcp是栅格矩阵,levels是颜色映射的level。似乎Cartopy是目前比较不错的选择,但是实在是太不稳定了,很多功能还不完备,目前的版本是0.21。原创 2022-11-18 12:20:58 · 1183 阅读 · 0 评论