- 博客(52)
- 资源 (2)
- 问答 (1)
- 收藏
- 关注
RSS订阅原创 python修改抖动图facecolor和箱型图alpha
抖动图和箱型图在数据分析和可视化中具有重要的作用，可以帮助人们更好地理解数据的特征和变化趋势。当我们想使用Python绘制抖动图的时候，可以使用seaborn的库来绘制，但其自身的参数限制了我们更改facecolor，绘制箱型图使更改透明度也难以更改，因接下来将针对这个问题进行解决。然而当你想要设置facecolor时会发现这个参数不起作用，实际上Python绘图本质是一系列的patch及collections等绘制，因此可以访问带内部的collections来修改颜色。
2024-04-07 21:43:51
343
原创 基于Python的简单颜色替换
当我们临时需要改变一个照片的颜色，使其符合我们想要的主题色时，对于不会PS的我就只能使用一下Python来实现这个简单的过程。
2024-04-04 00:00:00
183
原创 python实现裂区试验方差分析
方差分析（Analysis of Variance，ANOVA）是一种统计方法，用于比较三个或三个以上组别的平均值是否存在显著差异。它通过比较组内变异和组间变异的大小来判断组别间的平均值是否有显著差异。方差分析的基本假设包括对各组数据的正态性、方差齐性和独立性的要求。通过方差分析，可以得出组别间的平均值是否存在显著差异，以及各组之间的比较结果。方差分析是一种重要的统计方法，可以帮助研究者进行多组别的平均值比较和因素影响分析。数据要求的格式为长表，内容包括主区因子、幅区因子、区组和各个变量的数据这些列。
2024-02-24 08:47:15
990
原创 基于python的遥感影像灰色关联矩阵纹理特征计算
目前熵值特征计算结果与ENVI有差异，不过只是数值差异，使用影像打开的结果显示一致，说明只是值的范围差异，不影响分析，其他特征均与ENVI的计算结果一致。
2024-02-18 21:03:41
662
4
原创 作物模型的区域模拟实现过程初探
作物生长模型能够根据作物的生长环境模拟作物的生长发育，为农业生产及科学研究提供有价值的参考。但模型多数是基于单点或者几个站点的模拟结果进行差值来表征区域的模拟结果。这里来探讨一下模型的网格化模拟实现。一般的作物生长模型都需要土壤、气象、管理、初始条件等几个必要的因素。在区域模拟时还需要种植区域的栅格图。基本流程假设我们现在有一副小麦区域分布图 wheat.tif以下以Python代码为例提供一些过程的核心代码种植区域获取首先需要获取种植区域的坐标将栅格的种植区域读取为arrayfrom .
2022-05-28 21:41:41
1037
12
原创 R语言机器学习
简单介绍一下R语言的几个分类算法便于以后查找# install.packages("DMwR")# install.packages("kernlab")# install.packages("Hmisc")# install.packages("corrplot")# install.packages("mlogit")# install.packages("randomForest")# install.packages("e1071")# 决策树方法set.seed(42)#加载
2021-12-28 15:36:53
2288
原创 R语言判别分析
判别分析又称“分辨法”，是按照一定的判别准则，建立一个或多个判别函数，用研究对象的大量资料确定判别函数中的待定系数，并计算判别指标。据此即可确定某一样本属于何类。判别分析目的——对已知分类的数据建立由数值指标构成的分类规则，然后把这样的规则应用到未知分类的样本中去分类，以识别未知样本所属的类别。library('MASS')data(iris)head(iris)colnames(iris) <- tolower(colnames(iris))head(iris)ld <- lda
2021-12-28 14:35:52
1237
2
原创 R语言聚类分析
基本过程-读数据-数据标准化-计算距离-分类library('flexclust') # 导入数据集data(nutrient) # 加载数据row.names(nutrient) <- tolower(row.names(nutrient)) # 把行索引小写head(nutrient)#由于能量变化范围比其他变量更大，缩放数据有利于均衡各变量的影响。nutrient_scaled <- scale(nutrient) # 数据标准化head(nutrient_sca
2021-12-28 14:35:16
7708
原创 R语言主成分分析
主成分分析介绍可以看此视频library(FactoMineR)data("iris")iris_pca <- PCA(iris[,-5], graph =F)summary(iris_pca)library(factoextra)eig_val <- get_eigenvalue(iris_pca) # 获得特征值，其实summary里也有eig_valscree_plot <-fviz_eig(iris_pca)scree_plotfviz_pca_var(iri
2021-12-28 14:34:23
411
原创 R语言回归笔记
主要内容一元回归多元回归存在共线性的回归岭回归偏最小二乘一元回归library(rio)lm <- import('lm1.xlsx')fitlm <- lm(weight~N,data = lm)summary((fitlm))# 模型R20.911，截距和系数均极显著方程为weight=3N-87.51667summary((fitlm))Call:lm(formula = weight ~ N, data = lm)Residuals: Min
2021-12-28 14:32:58
484
1
原创 R语言统计分析之方差分析
R语言是一个很好的做统计分析的语言，本文对R的方差分析做个总结单因素方差分析数据sssdf <- read.csv('sss.csv')# 正态分布检验shapiro.test(df[which(df[, 2] == "A"), ] $x)shapiro.test(df[which(df[, 2] == "B"), ] $x)shapiro.test(df[which(df[, 2] == "C"), ] $x)shapiro.test(df[which(df[, 2] == "D"
2021-12-26 09:41:25
4860
1
原创 pandas实践之表格整理
试验数据整理目的一个试验有不同的处理，每个处理又有不同小区，每个小区不同指标有不同重复类似下图这个利用Excel的透视表功能和vlookup函数结合就能实现对小区的统计分析，不过稍显麻烦这里pandas涉及的主要操作有数据的分组统计（涉及多元索引）行列的计算空值填充（需要注意的是python计算均值的时候会忽略空值）分组后索引更改透视表的制作（包括透视表的制作与逆向生成）不同表之间的融合（类似VLOOKUP）日期类型的转换数据框多因素排序多s
2021-12-23 09:17:21
1793
原创 基于python的多光谱影像植被指数计算
利用GDAL库对多光谱影像进行分析，光谱5波段，顺序import numpy as npfrom osgeo import gdal, osrdef export_tif(out_tif_name, var_lat, var_lon, NDVI): # 分辨率计算 N_Lat = len(var_lat) N_Lon = len(var_lon) LonMin, LatMax, LonMax, LatMin = [min(var_lon), max(var_lat)
2021-12-09 19:18:30
3620
5
原创 大疆精灵4多光谱辐射信息分析
大疆多光谱无人机顶部设置了辐射传感器，可捕捉太阳辐照度并记录于影像文件中，当进行数据后期处理时，太阳辐照度数据将可用于对影像进行光照补偿，排除环境光对数据采集的干扰，有助于使用者获得更准确的 NDVI 结果，提高不同时段采集到的数据的准确度与一致性。但很多人说可能不准，因此本文对大疆的辐射数据进行提取查看结果首先对获取的...
2021-11-24 21:47:01
3235
原创 Python分类文件（大疆精灵4多光谱版PM4影象分类）
任务描述大疆精灵4多光谱版有六个通道，每次拍摄都会生成六张图片。图片的顺序为 0. 可见光 1. 蓝光 2. 绿光 3. 红光 4. 红边 5. 近红外。在用其他软件进行拼接时需要对六个通道的照片分别进行拼接，而不同通道的图片仅是最后一位数字不同，因此要将相同通道的图片放到一起，易于利用软件拼接。基于Python的os库和shutil库进行文件转移操作，代码如下`import shutilimport ospath = r’E:\博士期间记录\照片\站内\高光谱\20210811’files
2021-11-03 22:11:22
1057
2
原创 R语言基础学习笔记
最近在看一些R语言，顺便记录一下R语言的基础操作，以后查找方便R语言更多的用于数据分析中，算是一个比较有特色的一个语言学习的主要内容为数据类型数据结构函数控制流（判断循环）文件读取数据类型R数据类型主要为字符character，数字numeric，和逻辑值logicalis.numeric()可以判断类型as.numeric()可转换类型，这与一般语言基本类似，不过一般数值型分为浮点、整型等数据结构结构主要为向量 数组 矩阵 数据框 列表等，向量是R语言的基础，单个数字也是单个向量，
2021-08-05 22:53:12
555
原创 python 批量下载指定时间段的邮箱附件
邮箱附件下载在收到很多人发的邮件时，一个一个的接收附件比较麻烦，考虑使用python实现指定日期邮件附件的下载首先需要开启邮箱的pop协议，一般在设置里，以QQ邮箱为例，设置后会给个密码，这个就是程序访问邮箱的登录密码需要改的东西有：用户名密码文件储存路径时间import emailimport osimport poplibimport timefrom email.header import decode_headerfrom email.parser import Pars
2021-05-16 11:30:58
2422
3
原创 ImportError:DLL load failed:找不到指定模块的解决方法
今天换个电脑安装nlopt优化包后一直出现找不到C++编译的_nlopt模块，包的内容如下图所示。网上说的原因可能是路径不对、版本不对然后在非官方python下载库里下载64位的nlopt, 通过本地安装完美解决问题64位里面长这样说明以前的版本和64位不太兼容...
2021-04-18 08:29:50
262
原创 非线性优化算法求解线性规划（python）
Nlopt是一种求解非线性模型最优解的一种集合算法，尝试一下线性模型中的应用问题：x1+x2+x3<=956x1+5x2+2x3<=4005x1+2x2<=20012x1+10x2+16x3<1200x1,x2,x3>=0使得6x1+4x2+3x3最大这个其实可以用Excel或者单纯搜索法手动计算，这里测试下nlopt中的SUBPLEX算法在上面的应用import nlopt # 导入模块# 定义代价函数class ObjectiveFunction(
2021-04-05 10:58:18
1602
3
原创 python 提取图片中绿色植被，计算冠层覆盖度
任务描述提取绿色植被部分对应的红外光谱（即红绿蓝三波段提取绿色，掩膜到红外波段）import cv2import numpy as np# 第一种办法--------------------------------------------------------------im=cv2.imread('FLIR06.jpg')#提取图像的三个通道B, G, R = cv2.split(im)#计算植被指数cive = 0.441*R-0.811*G+0.385*B+18.78745g
2021-02-20 19:58:23
4997
7
原创 python 绘图，绘制模型模拟与实测值对比图
模型模拟值一般为连续的值，实测值为点，且一般带有误差# 设置绘图地点sites = ['内黄', '许昌', '洛宁', '罗山']sites_e = ['NH', 'XC', 'LN', 'LS']annote = ['a', 'b', 'c', 'd']data_dir = r'C:\Users\Administrator\Desktop\model calibration' # 模型校准文件夹weather_dir = r'C:\Users\Administrator\Desktop\2
2021-02-17 12:27:56
5709
3
原创 S-G滤波之包络滤波
由于遥感观测总是受到云层，大气影响，总会造成植被指数骤降问题，包络滤波的就是根据植被的特性，将骤降值认为是错误值，选择最大值进行滤波SG滤波代码参考&后来&博客定义滤波函数import pandas as pdimport matplotlib.pyplot as plt#%% SG_filterimport numpy as npfrom pandas.plotting import register_matplotlib_convertersregister_matplo
2021-01-23 19:29:03
3123
3
原创 PCSE.WOFOST模型调参
根据实测值对模型进行调参先利用土壤物理参数对土壤的水分特征进行估计，由于气象文件分辨率为0.5度，定义函数使坐标转换为对于气象文件的坐标。# 参数定义def soil_data(s, c, om, b=1.5): s = s # sand % c = c # clay % om = om # organic matter % b = b # bulk g/cm3 # 永久萎蔫点 PWP pwp = 1.14*(-0.00024*s+0.004
2021-01-02 16:06:02
4330
7
原创 python的gdal库，读取tif影像为DataFrame,或将DataFrame转为tif
将tif转为DataFrame存为csv此数据为modis的LAI产品，仅有LAI一个波段from osgeo import gdalimport pandas as pdgdal.AllRegister()filePath = r'C:\Users\Administrator\遥感文件\MOD15A2H.006_Lai_500m_doy2020097_aid0001.tif'dataset = gdal.Open(filePath)adfGeoTransform = dataset.Ge
2021-01-01 20:04:39
2484
原创 python：有序的列表中模糊匹配给定的数值
给定一个值，在列表中查找近似的对应的值列表需要有序，正序逆序均可from bisect import bisect_leftdef takeClosest(myList, myNumber): if myNumber >= max(myList): return max(myList) elif myNumber <= min(myList): return min(myList) if myList[0] < myList[1
2021-01-01 17:17:02
878
原创 用python对Excel每列不确定个数的元素转换为行
目的如上图，每列有很多相同的数值，想要每个编号的病级转成一行但手动复制太慢先读取数据获取每个编号的个数import pandas as pddata = pd.read_excel(r'C:\Users\Desktop\data.xlsx')count = data.groupby(by='编号').sum()创建一个字典储存结果，对data进行循环，先循环编号，再循环编号内数据result = {}value = []h=0for i in range(359): fo
2020-11-04 19:00:05
375
1
原创 python将nc文件转为tiff
import numpy as npimport netCDF4 as ncfrom osgeo import gdal,osrvar = 'SA'data = r'C:\Users\13290\Desktop\soil data\{}.nc'.format(var)f = nc.Dataset(data)var_lon = f['lon'][:]var_lat = f['lat'][:]data = f[var][0, :]data_arr = np.asarray(data)dat.
2020-08-06 11:24:49
7812
21
原创 netCDF4(.nc)文件读取转为tif或csv（python）
.nc文件（network Common Data Format）文件是气象上常用的数据格式，python上读取.nc使用较多的库为netCDF4这个库。基本操作参考fangzuliang的博客"""author: shuaijieintro: date: 08/03/2020 18:47"""import netCDF4 as ncimport matplotlib.pyplot as pltimport pandas as pdimport seaborn as snsimport
2020-08-05 17:34:40
6559
8
原创 对有重复的结果筛选
结果每三个为一个重复，选择其中两个差距小的平均数据如下"""author: shuaijieintro: 在三个测氮的值中选择两个相近的平均date: 07/27/2020 11:28"""import pandas as pddef main(): fp = pd.read_excel(r'C:\Users\admire\Desktop\测氮结果示例.xlsx') # 读取数据 result = [] identify = [] for i in r
2020-07-27 16:22:09
207
原创 python对有重复值的列求和之后放到第一次出现的位置
对第二列求和，但是重复的要把数据放到第一次重复的地方，其他的空着。整体思路是计数，之后建立一个列表储存列信息，之后相乘。"""author: shuaijieintro: 分类求和date: 07/22/2020 19:17"""import pandas as pddef main(): data = pd.read_excel(r'C:\Users\13290\Desktop\施肥量.xlsx') sum_data = data.groupby(by='姓名', s.
2020-07-22 21:18:05
803
原创 prosail模型敏感性分析
敏感性分析一般借助simlab软件，具体敏感性分析的方法可以参考https://www.jianshu.com/p/2656af9f8d81首先通过simlab软件进行样本采样，生成.SAM文件（我定义了7个参数），然后读取文件中的参数值，代入PROSAIL模型不了解PROSAIL的可以参考https://blog.csdn.net/weixin_45452300/article/detai...
2020-04-06 22:07:55
6859
39
原创 PROSAIL模型python版简单说明
PROSAIL模型是prospect 和sail模型的融合，其本来是用fortran语言，现在也有python版的，可以直接pip导入PROSAIL库。其参数ParameterDescription of parameterUnitsTypical minTypical maxNLeaf structure parameterN/A0.82.5cabC...
2020-04-05 22:50:15
12120
47
原创 对asd光谱数据的多个TXT文件进行求平均
对文本文件进行操作，将多个文本文件求平均合并为一个文件。def sv(f, f1, waves, wave_number, block): f1.write('waves\t') for i in range(waves): i += 1 f1.write(str(i)) f1.write('\t') f1.write('\...
2020-04-05 20:51:44
648
GDAL-2.2.4-cp37-cp37m-win_amd64.whl
2020-08-06
pyodm创建任务时无法获得node信息
2021-08-28
TA创建的收藏夹 TA关注的收藏夹
TA关注的人