小琳子要开心呀-CSDN博客

原创 ubuntu操作学习

升级包sudo apt update # 更新apt软件源信息sudo apt upgrade # 更新已安装的包卸载包sudo apt remove package # 删除包sudo apt purge package # 删除包，包括删除配置文件等安装包sudo apt install package # 安装包。

2024-12-26 06:59:32 521

当对染色体未按字典排序的文件使用 -sorted 选项时（例如，对 BED 文件使用 -k1,1 -k2,2n 排序），必须提供定义预期染色体顺序的基因组文件（-g）。使用TASSEL学习GWAS笔记（3/6）：基因型数据可视化：kingship，PCA，MDS。使用TASSEL学习GWAS笔记（4/6）：一般线性模型进行GWAS分析（GLM模型）使用TASSEL学习GWAS笔记（5/6）：混合线性模型进行GWAS分析（MLM模型）-（知乎）使用TASSEL学习GWAS笔记（1-6）完整版。

2024-12-26 06:40:11 824

转载 Zotero+OneDrive多平台在线同步

Zotero+OneDrive多平台在线同步

2024-10-17 04:07:32 304

原创 R语言使用ggplot2包绘图

参考链接：https://blog.csdn.net/weshengxin/article/details/126006349。

2023-05-09 12:48:30 5616 2

原创基于R语言的机器学习学习笔记（多元线性回归, 随机森林，XGBoost）

R语言importance {randomForest}关于变量重要性度量（the variable importance measures）的定义。以下是变量重要性测量的定义。第一个衡量标准是。

2023-05-06 11:34:16 2855

原创 R语言数据处理

R语言中，自定义函数的一般格式为：函数名

2023-03-31 21:10:33 1597 1

原创 python数据处理

参考链接: https://www.yisu.com/zixun/225360.htmldef foo() : #定义函数foo()， m = 3 #定义变量m=3;def bar() : #在foo内定义函数bar() n = 4 #定义局部变量n=4 print(m + n) #m相当于函数bar()的全局变量 return(bar()) #foo()函数内调用函数bar() foo()

2023-03-31 10:21:22 2249 1

原创基于python的ply格式点云数据处理（学习笔记）

open3D直接去看下面推荐的博客和官方文档吧，不再赘述。此博客很多地方参考了下面推荐的博客，后文不在重复指明引用。本文的主要区别是针对我遇到的问题进行学习并整理。官方文档（英文）csdn：梦醒Blue（open3d教程，翻译自官方文档，免费）csdn：点云侠（讲的详细，部分来自官方文档，收费）一般情况读写和处理点云数据使用 open3d 库就够用了，但是我的数据里面字段太多，只有pyntcloud库才能全部读取。

2023-03-31 09:36:27 6859 1

原创 R语言使用ggplot2绘制直方图

R语言使用ggplot2绘制直方图

2022-05-16 22:48:40 2549

原创 R语言使用corrplot绘制相关性图

R语言使用corrplot绘制相关性图

2022-05-16 22:40:25 2150

转载 matlab画图标题中包含下划线的正确显示方法

这里写自定义目录标题转载自：https://blog.csdn.net/lingyunxianhe/article/details/80469347方法：添加 ‘Interpreter’,‘none’代码如下：figureplot((1:10).^2)title('X_1','Interpreter','none')

2021-10-15 09:19:26 2048

原创群体中的人离不开人际关系的烦恼

成年后的烦恼接踵而至，没有良药，只有烟酒。读本科的时候，活的迷茫，因为没有梦想哭过，想来幼稚，却也单纯。后来选了一条路，暂且当做梦想去追求。却接踵而至的烦恼让我怀疑，我的思维和情商配不上我的年龄。25的年纪应该懂不少，却依旧固执。价值观封闭的像一座城，得了疫情的城，进不来也出不去。上一次，用白纸列出七八条心病。终于慢慢的解决了许多。这一次，懒得在写，写来写去，多是人际关系的烦恼。非要什么良方，也只有律己以圣，容人以凡。求圣需修心，苦修的苦，有几人吃得。做凡

2021-08-31 00:18:33 346 3

原创使用R语言和软件计算光响应曲线（表观量子效率，光饱和点，光补偿点，暗呼吸）

一、什么是光响应曲线二、使用Li-cor6400测量光响应曲线三、拟合光响应曲线3.1 使用软件拟合（叶子飘老师）3.2 使用R语言代码拟合（四种模型）3.3 使用R语言程序包（photosynthesis）拟合（基于非直角双曲线模型）3.4 使用R语言包(FitAQ)拟合（基于非直角双曲线模型）

2021-08-28 14:53:53 16234 6

原创 excel字符串操作（截取，拼接，替换）

函数：mid需求：提取A1单元格字符串的一部分（第四个字符算起，截取2个字符）放在A2单元格。使用：=mid(A1,4,2)

2021-08-14 16:07:18 7665

转载遥感中的数字量化值DN, 辐射亮度Radiance, 反射率Reflectance,发射率Emissive的意义

转载自：http://blog.sina.com.cn/s/blog_764b1e9d0101c7e6.htmlhttps://www.ixxin.cn/2017/05/27/rs_data_jichu/数字量化值（Digital Number ：DN）像素值的通用术语是数字量化值或DN值，它通常被用来描述还没有校准到具有意义单位的像素值。DN值是卫星接收到辐射进行量化的一个数值，大小与量化深度有关，比如说8bit量化，那么Max(DN) = 255，无单位，无实际意义。如果你只是想

2021-08-04 22:37:51 14881 4

转载批量修改文件名（Dos命令导出文件列表+Excel软件）

转载自知乎：https://www.zhihu.com/question/266365274/answer/307140022基本步骤如下：1.读取需要修改的文件名。条件：需要修改的文件都放在了一个文件内。1.1 需要修改的文件夹内，shift+右击空白处，单击“在此处打开Powershell窗口（S）”。1.2 输入下面代码可将文件名全部读取到本文件夹的1.txt文件内。dir *.* >1.txt可自行修改，比如我需要读取文件名带有N0P1的后缀为.asd的文件.

2021-08-03 13:52:37 949

原创经纬度（度分秒）和十进制相互转换

一、度分秒转十进制http://www.ab126.com/shuxue/10801.html一、十进制转度分秒https://www.osgeo.cn/app/s2082

2021-03-23 12:23:37 5218

原创 matlab 读取并拼接已经绘制的图像，去除拼接处的轴刻度标注

问题：已经画好多张图，如何在连接处，去除轴去除轴标注，并拼接在一起。思路：前提：所有图像使用matlab加载tif图像（imread）,读取行列号，根据行列号提取相应范围的tif数据（去除轴刻度标注）。不推荐使用，个人代码备份。代码：path = 'F:\papar_china_GPP-SIF\geotiff\geotiffwrite\多年平均\合并\';file = dir(fullfile(path,'*.tif'));len = length(file)for i =..

2020-12-22 21:38:04 620

原创土壤水分下载链接（soil moisture）

根区土壤水分下载链接（root_zone_soil_wetness）https://disc.gsfc.nasa.gov/datasets?keywords=soil%20moisture%20L4&page=1&source=Models%2FAnalyses%20MERRA-2Shortname: M2TUNXLNDLong name: MERRA-2 tavgU_2d_lnd_Nx: 2d,diurnal,Time-Averaged,Single-Level,Assimil

2020-12-20 14:23:08 14622 9

原创 matlab——计算VPD（vapor pressure defict）

需求：计算计算VPD（vapor pressure defict）。介绍：饱和水汽压差（简称VPD）是指在一定温度下，饱和水汽压与空气中的实际水汽压之间的差值（百度百科）。因此，温室中VPD的理想范围是0.45kPa至1.25 kPa，理想情况下约为0.85 kPa。通常，大多数植物在VPD在0.8到0.95 kPa之间时生长良好（维基百科）。计算方法：一、先计算饱和水汽压二、饱和水汽压减去实际水汽压。世界气象组织推荐的Goff-Gratch公式如下（计算饱和水汽压）。..

2020-12-20 00:54:16 14310 7

原创 arcgis制图汇总

目录镶嵌小图（中国地图右下角绘制南海区域小图）将arcgis的图例以色条形式横放制作经纬格网地图将格网放在图层下面arcgis官方文档https://desktop.arcgis.com/zh-cn/arcmap/10.3/main/map/mapping-and-visualization-in-arcgis-for-desktop.htm镶嵌小图（中国地图右下角绘制南海区域小图）看起来复杂，其实简单。1、菜单栏-插入-数据框。2、新插入的数据框缩小放在中国地图的右

2020-12-18 16:55:01 6404 1

转载 origin画误差阴影

描述：origin画误差阴影参考：https://www.51xxziyuan.com/58/957.html有资源的可以直接使用origin2020版（bilibili有官方账号，以及origin试用版安装教程），这个版本软件可以直接绘制误差阴影。

2020-12-17 10:27:09 5510

转载 ArcGIS横向（水平）图例（色带）

https://blog.csdn.net/qq_33356563/article/details/83998963?utm_medium=distribute.pc_relevant_t0.none-task-blog-searchFromBaidu-1.control&depth_1-utm_source=distribute.pc_relevant_t0.none-task-blog-searchFromBaidu-1.control

2020-12-16 17:12:55 7932

原创 arcgis10.6打开属性、选项等对话框出现卡顿（官方补丁链接）

问题:arcgis10.6打开属性、选项等对话框出现卡顿，十分卡卡卡卡卡卡卡卡卡卡卡卡卡卡，哭死解决：软件官方的锅，有bug，自己别瞎整了，官方提供了补丁，打上就好了。下载地址：（官方补丁下载地址）https://support.esri.com/zh-cn/Products/Desktop/arcgis-desktop/arcmap/10-6-1#downloads?id=7678（知乎大佬提供的）https://www.zhihu.com/question/384752775以

2020-12-16 10:56:58 4492

转载 matlab截取字符串

描述:使用matlab截取固定位置字符串.函数:方法:1.使用冒号A = 'abcdabcd'B = A(2:3)B = A(4:end)B = [A(2:4), A(6)]2.使用函数 strfind查找下标位置A = 'abcdabcd'idx = strfind(A,'c')B = A(inx)参考:1.百度经验(https://jingyan.baidu.com/article/f3ad7d0f5e23b409c3345beb.html)...

2020-12-15 11:15:04 13418

转载如何在Matlab循环中动态更改变量、文件、函数名

转载:如何在Matlab循环中动态更改变量、文件、函数名https://jingyan.baidu.com/article/fc07f9895fdb2712fee51978.htmlfor i = 1:10 b = ['a',num2str(i)]; eval([b,'=i'])end生成十个已经赋值的变量: a1到 a10

2020-12-13 16:49:21 2957 1

原创 matlab 画折线图并美化

使用函数:plot (绘图)；set（设置坐标属性）；print（输出图片）代码如下：X = 1:9;Y = 10:10:90;f2 = figure(2);set(gcf,'unit','centimeters','position',[0,0,8,6])plot(X,Y,'-ok','LineWidth',1,'MarkerSize',5,'MarkerFaceColor','k')set(gca,'FontName', 'Times New Roman', 'FontSize'

2020-12-11 16:17:44 7637

转载所有能反演SIF的卫星以及下载链接

转载自科学网张园博客：http://blog.sciencenet.cn/blog-3430406-1223659.html目录SCIAMACHYGOME-2GOSATOCO系列卫星轨道碳观察者TanSatTROPOMI调研了一下叶绿素荧光（SIF）卫星的相关情况，看看哪种数据产品更适合对森林的研究。六个目前能获取数据的传感器，部分参数如下：SCIAMACHY2002年3月1日发射，2012年4月8日失联。回归周期6天...

2020-12-10 19:28:59 6282 5

原创 PPT另存为图片分辨率太低，设置高分辨率（300dpi）

问题：PPT另存为图片分辨率太低解决方案1、下载设置脚本（推荐）参考：如何提高PowerPoint另存为图像的分辨率？ - 木雷的文章 - 知乎 https://zhuanlan.zhihu.com/p/23555929下载链接：PPT_DPI设置脚本.bat2、使用cmd程序windows搜索CMD，打开命令提示符，输入下面代码（设置另存为图片300dpi）。reg add HKCU\Software\Microsoft\Office\16.0\PowerPoint\Op

2020-12-09 10:03:38 32391 3

原创 arcgis 投影变换

目录1、目的：2、背景知识：2.1、地理坐标系统与投影坐标系统的区别（看单位）地理坐标系统投影坐标系统2.2、arcgis中“定义投影”与“投影”或者“投影栅格”的区别定义投影投影2.3、中国常用的投影坐标系统地理坐标投影坐标投影转换方法1、目的：WCS_WGS_1984 转成Krasovsky_1940_Albers2、背景知识：参考：坐标系统与投影变换及在ARCGIS中的应用（https://blog.csdn.net/jax_lee/a.

2020-12-08 16:59:32 20392 9

原创 arcgis 删除创建的自定义地理变换文件

地理变换文件默认存储路径为：C:\Users\Administrator\AppData\Roaming\ESRI\Desktop10.3\ArcToolbox\CustomTransformations找到删除即可。参考：https://www.cnblogs.com/yiliangmi/p/12402710.html

2020-12-08 16:06:17 2453

原创 arcgis 去除影像黑色边框（nodata）

目的：加载进arcmap的tif数据发现有黑边。想要去除黑边。问题：1、在属性中设置“将NoData显示为绿色”，发现对黑色区域无效。2、使用“识别”功能发现黑色区域为“nodata”，而绿色区域没有数值，为空值。原理：需要知道没有数据与显示“NoData”不一样。黑色区域，“识别”显示“NoData”，表示此处没有有效数据，使用“NoData”存储。绿色区域，“识别”显示“为空值，表示此处没有数据。解决办法：1、使用裁剪工具（数据管理-栅格-栅格处理-裁

2020-12-08 15:25:28 13957 6

原创 MOD15A2(LAI,FPAR)

数据下载网址：https://ladsweb.modaps.eosdis.nasa.gov/missions-and-measurements/products/MOD15A2H数据介绍：MODIS 4级全球叶面积指数（LAI）和光合有效辐射分数（FPAR）产品是正弦网格上8天，500米分辨率的产品。MOD15A2H中提供的科学数据集包括LAI，FPAR，质量等级和每个变量的标准偏差。LAI变量定义相对于单位地面面积的等效叶层数，而FPAR则测量被树冠吸收的光合有效波长中可用辐射的比例。

2020-11-24 16:23:29 10587 3

原创气象数据集（MERRA-2再分析数据集）

数据下载链接：https://disc.gsfc.nasa.gov/datasets?project=MERRA-2数据介绍链接：https://gmao.gsfc.nasa.gov/pubs/docs/Bosilovich785.pdf文章链接：https://journals.ametsoc.org/jcli/article/30/14/5419/97087网页截图

2020-11-11 17:26:05 12251 2

原创全球气象数据下载链接 Climatic Research Unit (CRU) 和MERRA-2

链接1 Climatic Research Unit (CRU)：https://catalogue.ceda.ac.uk/uuid/10d3e3640f004c578403419aac167d82链接2European Climate Assessment & Dataset：http://climexp.knmi.nl/selectfield_obs2.cgi?id=someone@somewhere网页1：网页2：...

2020-11-07 00:00:05 8414

原创 matlab代码——画全球图并添加海岸线与颜色条

要点imagesc画彩色影像图 imagesc去除背景添加全球海岸线使用view镜像垂直翻转图像添加colorbar颜色条，并设置长宽，位置调整坐标轴刻度，标签，字体，大小输出为tiff格式（设置尺寸，比例，单位：cm）代码如下：自用备份%% 画图(全球逐项元相关图)load('G:\DATA\5MOD12Q1_UMD\out_matlab\LC_05_50.mat')figure(1)X = [-180 180];Y = [90 -90];f = image..

2020-10-22 10:29:13 4677

原创 matlab num2str怎么设置小数点位数

str = {[num2str(year(y)), 'ing: ', num2str(100*i/1200,'%4.2f'), '%', ' Time: ', num2str(toc/60, '%4.2f'), 'm'], [' pre: ',num2str(toc*20/i,'%4.2f'), 'm', ' rem: ', num2str((toc*20/i-toc/60),'%4.2f'),'m']};1. num2str 设置小数点位数num2str(1.12345...

2020-10-07 16:08:42 13625