![](https://img-blog.csdnimg.cn/20201014180756918.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
python
Soul_taker
这个作者很懒,什么都没留下…
展开
-
用降水、比湿和温度计算相对湿度(nc版、python)
有时收集到的气象数据为比湿,但模型中需要的输入为相对湿度,我们可以利用降水、比湿和温度进行换算。还可以批量计算nc格式的哦原创 2022-08-25 23:42:48 · 4068 阅读 · 6 评论 -
合并年nc数据并统计生长季数据
利用python做lai的基础处理,可以进行数据的生长季尺度提取,年尺度合并原创 2022-06-23 11:56:15 · 583 阅读 · 0 评论 -
针对QGIS二次开发时环境搭建的坑(解决no module named qgis、no module named processing)
qgis的python环境搭建,解决no module named qgis/no module named processing问题原创 2021-12-15 20:44:55 · 3323 阅读 · 2 评论 -
气象专业文件nc的读取、裁剪与输出(python)
气象专业文件nc的读取、裁剪与输出如今的大数据时代,很多气象数据都是nc格式的,尤其外网的大尺度数据,怎么进行nc文件的处理呢?当然arcgis也可以,不过我推荐python中的xarray包。xarray包相当灵活便捷,我写了一段nc文件读取、裁剪和输出的函数,供大家参考。import xarray as xrimport matplotlibdef read_clip_output_ncdata(input_data_path,output_data_path,min_x,min_y,max_x原创 2021-11-11 10:49:59 · 12408 阅读 · 8 评论 -
gee数据批量合成及下载(python)
gee数据批量合成及下载写了一个在Gee中最大合成为年、月数据,然后批量下载的程序。python客户端链接import eeimport geemapimport osimport geetools# we need to set this before initialize, we can get this from 工具-internet选项-链接-局域网设置os.environ ['HTTP_PROXY'] ='http://127.0.0.1:1080'os.environ [原创 2021-09-13 15:29:48 · 872 阅读 · 0 评论 -
python图像签名背景消除
python图像签名背景消除利用python的PIL库可以进行图像的二值化,从而将图像的背景删除,可以参考https://blog.csdn.net/qq_16669583/article/details/106318056?ivk_sa=1024320u博主的文章,我又加了一个自动筛选、批量处理的功能。#coding:utf-8from PIL import Imagefrom PIL import ImageOpsimport osimport redef initTable(thres原创 2021-08-20 16:03:53 · 537 阅读 · 0 评论 -
基于GEE的年最大合成、sen+mk检验算法
基于GEE的年最大合成、sen+mk检验算法1、实验平台与数据实验平台为google earth engine、python3.6,数据为gee中自带的NDVI数据。2、思路对长序列NDVI数据(8天合成数据)进行NDVI年最大合成, 然后进行sen+mk检验,检测NDVI的变化趋势及显著情况。3、代码年最大合成def img_max_year(start_year,end_year,data_path,data_name): ''' this function can get原创 2021-06-23 10:03:01 · 4103 阅读 · 9 评论 -
python数据清洗小计(气象)
需求数据格式如下:1-7分别代表7个大气模式对站点数据多年的相关系数,根据显著性检验,0.231为临界值,我们需要1)找出小于临界值的数据,赋值为0;2)将不是0的数据从大到小排序,并赋权重,各权重相加为1;3)清洗后的数据,模式编号和站点不能变化。代码实现导入包,获得数据import pandas as pdimport numpy as npdata = pd.read_excel("ttt.xlsx",sheet_name="权重")将index赋值为站点,并给所有数据做绝对值,避原创 2021-03-25 10:55:13 · 553 阅读 · 0 评论 -
knn python 手动实现
为了更好的理解knn算法,用python进行编程实现。原创 2021-03-24 17:21:55 · 291 阅读 · 0 评论 -
利用python将月数据统计为季节数据
利用python将月数据统计为季节数据在时间序列分析中,我们经常利用不同的时间尺度对数据进行描述和分析。下面利用python实现月度数据转换为季度数据的方法。数据形式输入数据为Excel,数据形式如图,第一行为标签,时间、年、月、数据。代码代码比较简单,pandas实在好用。import pandas as pddef season_count(path,output): ''' Get the seasonal data of each year :param pa原创 2021-03-15 20:59:42 · 5769 阅读 · 3 评论 -
时间序列分解随笔
时间序列分解对于时间序列来说,传统的时间序列预测模型对于不平稳的数据是不能进行处理的,需要对时间序列进行分解,使得时间序列区域平稳,才可以进行时间序列的模拟预测。本次采用python的seasonal_decompose分解来进行时间序列的分解。工作环境python3.6、statsmodels库、pandas库、matplotlib库。时间序列分解导入包from statsmodels.tsa.seasonal import seasonal_decomposeimport pandas原创 2020-10-28 19:58:49 · 587 阅读 · 0 评论 -
批量下载科研数据
批量下载科研数据(纪念我的第一个爬虫)最近需要下载NASA的一些全球数据集,由于时间序列较长,一个个点链接下太费劲,就搞了一个简单的爬虫,以供大家方便下载数据。环境python3.6、requests、pandas、time需求分析我要下载的是GLDAS数据集,这个集合中包含不同尺度的数据,首先在https://disc.gsfc.nasa.gov/information/data-release?title=New%20and%20Reprocessed%20GLDAS%20Version%20原创 2020-10-02 09:09:33 · 704 阅读 · 0 评论 -
GEE 提取长系列NDVI像元值
GEE 提取长系列NDVI像元值跟gee的python本地环境战斗很久,发现自己的网上的不太科学,所以只能用gee的线上Colab版本,对gee 的配置不熟悉的伙伴可以跳转大神的网站看看https://blog.csdn.net/shi_weihappy/article/details/93849370,相关介绍啥的比较详细了,向大神致敬。今天主要介绍利用GEE进行长系列NDVI像元值的提取。工作环境gee线上Colab版本即可。初始化GEE环境import ee# Trigger the au原创 2020-09-15 19:27:19 · 3160 阅读 · 2 评论 -
MK趋势分析
MK趋势分析Mann-Kendall(MK)趋势分析在气象、水文学中已得到广泛应用,该方法不要求变量遵从正太分布,故可以检验水文、气象序列。计算原理对于一个时间序列Xt = (x1,x2,…,xn),MK趋势检验的统计量S计算公式如下:其中,sgn为检验函数:然后计算方差Var以及Z统计量:利用双边检验对数据进行检验,当-Z1-a<=Zs<=Z1-a时,即表明序列具有显著变化趋势,反之趋势不显著。当Z1-a给定的显著水平a=0.05时,临界值为1.96。工作环境Python3原创 2020-07-02 18:48:27 · 13205 阅读 · 4 评论 -
arcgis批量出gis图
arcgis批量出gis图接上一篇,把批量导出的图层利用统一底图,批量导出为gis文件。工作空间arcgis10.2,Python2.7以及arcpy。导入包如下:# -*- coding:utf-8 -*-import sysreload(sys)sys.setdefaultencoding('utf-8')import arcpy as arcimport osimport re读取文件夹下文件在path中写入原mxd文件路径,以便读取文件名(有点儿多余,懒得改了)。new_p原创 2020-06-30 17:50:43 · 1744 阅读 · 4 评论 -
arcgis批量对比,导出图层
arcgis批量对比,导出图层最近工作中遇到了这么一个问题,把套好底图的mxd中的关键数据导入另一个底图当中,也就是需要换底图。一个个的换实在费劲,开发了新的程序进行批量导出。首先就是对关键数据的图层导出,我们都知道,arcgis中只导出数据的话,调好的颜色、分类的等会变,需要在新图上重新调整,所以我们需要把调整好的关键数据导出为图层形式,这样就不再需要重新对图层进行调整啦。今天我们就来看看如何通过对比批量导出图层数据。1.工作环境工作环境依然是arcgis10.2和python2.7,所需的库有a原创 2020-06-23 13:20:32 · 3956 阅读 · 0 评论 -
arcgis批量处理栅格数据
arcgis批量处理栅格数据在ArcGIS中,栅格数据的处理比较繁琐,自带的工具箱中可做的统计有限。本次借助arcpy对栅格数据进行批量处理,并可以实现栅格到numpy数组的转变,从而可以对其进行中位数、平均数、众数等更多元化的统计分析。1. 工作环境选用的arcgis版本为10.2,Python版本为2.7,还需要numpy、arcpy等库的支持。2.导入所用的库# -*- coding:utf-8 -*-import arcpyfrom arcpy import envimport nu原创 2020-06-15 19:22:31 · 3524 阅读 · 0 评论 -
利用python进行png图像的读写操作
利用python进行png图像的读写操作最近遇到了对png灰度图像进行读取,并统计图像中众数的需求,现将代码记录如下。ps:适用于单波段图,形如(X,Y)。1.导入需要的包import matplotlib.pyplot as plt # 绘图import matplotlib.image as mpimg # 显示图像import numpy as np # 处理数据from sci...原创 2019-10-13 15:36:11 · 14955 阅读 · 0 评论 -
关于数据的批量回归插补缺失值
关于数据的批量回归插补缺失值缺失值是数据分析工作中常见的数据问题,我们通常通过中位数、众数、平均数等进行缺失值插补,对于需要高精度插补的水文气象等数据,我们经常使用回归的方式插补,但是,利用excel回归插补需要耗费大量的经历,尤其站点较多的情况下,更是让人头疼,利用python,可以快速方便的进行数据的回归插补。数据形式数据的形式大体为上图,站点有100多个,年份是51到18年,可见数据...原创 2019-09-23 10:04:07 · 3183 阅读 · 2 评论 -
关于pandas里gruoupby函数的理解
关于pandas里gruoupby函数的理解最近在学习pandas库,pandas里面有很多好用的函数,今天来记录以下groupby函数。数据准备首先建立演示数据。import pandas as pddf = pd.DataFrame({'Animal': ['Falcon', 'Falcon', 'Parrot', 'Pa...原创 2019-09-19 20:58:15 · 450 阅读 · 0 评论 -
批量txt数据转换为excel
批量txt数据转换为excel日常工作中经常会遇到大量的数据处理,尤其是大量的txt文件需要转换为excel,只能一个个粘贴复制,有了如下代码,分分钟即可完成繁琐的数据处理工作。环境准备需要python和python里的pandas库和openpyxl库,如果没有可以到管网下载python,并利用pip install + 库名进行pandas库和openpyxl库的安装。代码首先,我们...原创 2019-09-18 17:48:29 · 16039 阅读 · 5 评论