weixin_49232187的博客

原创 【Zotero】【国标csl调教（八）】注意卷次号丢失的问题

即使是CNKI导入（并非百度学术），也会偶现数据字段丢失的问题。，即使笔者已经把csl搞定，也避免不了数据字段的问题。常见如下：年份直接加括号，表示没有卷次号。Zotero中如下显示：没有卷次。需要自行添加，结果如下。

原创 【Zotero】【国标csl调教（七）】导入专著M、学位论文D以及百度学术的问题

百度学术导入的字段（期卷号、页码）等会有严重错误，不建议通过Baidu Scholar导入，但是Google导入的专著、CNKI导入的学位论文并不一定有出版地，需要自行添加。国标对硕博论文【D】和专著【M】要求写上出版地。添加后在参考文献表中refresh即可。建议在文库编目字段自行查找修改。

原创 【ArcGIS】比例尺报错/比例尺不匹配处理

很久很久以前遇到的情况，比例尺数据明显与图片不匹配，且排除了Grid中显示比例尺单位和实际单位的问题。关键在于后面，tab是选项卡的意思，即DataFrame 的属性中General选项卡。

原创 【Arcpy】批量提取栅格数据到点（Extract Multi Values To Points ）

由于总列数超过255，故只能一年31个数据输出一次Excel，然后清空属性表开始下一次提取。2011-2019年3月逐日数据，提取到点属性表后再输出Excel。4.InRasterList变量的简单写法。5.清空字段部分的简单写法`

原创 【ArcPy】批量删除shp属性表中的字段

操作时要关闭在ArcGIS打开的该shp，否则会提示schema lock。先前批量操作添加了很多字段，重新调试想要清空，只留下原有的必要字段。

原创 【Arcpy】不同文件夹下的栅格计算代码构建

上篇文章针对的是同一文件夹下的栅格文件，在env.workspace中设置工作空间，使用ListRasters获取栅格列表然后利用Python切片功能构造相邻文件作差的效果本篇将解决不同文件夹下栅格文件的提取，以及代码构建过程和常见bug处理。

原创 【Arcpy】相邻时相栅格数据批量作差的代码构建

ListRasters的返回值是String，所以后面需要用arcpy.Raster创建栅格对象。文件夹中有一个月的栅格数据，试将相邻两天数据作差（前一天-后一天），并以后一天数据命名。《面向ArcGIS的Python脚本编程》人民邮电出版社。

原创 【ENVItask（二）】IDL批量计算TVDI的代码构建思路及bug

给出了基于EnviAPPStore中的TVDI工具task的批量计算TVDI的IDL代码，本文介绍对于ENVItask的学习过程和代码思路。

原创 【ENVItask（一）】批量计算TVDI的方法及IDL代码

在ENVI安装目录下的F:\ENVI\ENVI53\custom_code文件夹找到TVDI工具对应的task文件，根据文件参数编写程序。视频 【无限可能的ENVITask开发技术】 https://www.bilibili.com/video/BV1K14y1t7R9/利用IDL调用task及自定义task参考以下文章和视频。Hanzt大大编写的ENVI拓展工具可以计算TVDI。

原创 【Zotero】【国标csl调教（六）】author date格式正文引用中只有两个作者 需要表示为A and B的设置

注意，在线工具导出的csl会自动把原来的注释全部删除掉，如图，在notepad++中利用compare插件对比原csl和在线工具导出的csl。一般要求，author date格式中正文引用（citation）如果只有两个作者，要写成A和B（A and B）的形式，但是csl中设置。默认这里and是空，需要显示and（中文自动变成“和”）选text，1.and只修改citation，与bibliography无关。添加and=“text”，注意前后有空格。定位在name form这一行（

原创 【Zotero】【国标csl调教（五）】同一位置引用多个文献的方法

在这里插入图片描述](https://img-blog.csdnimg.cn/direct/5c6ba022a99d47af84a0bdb9481b9499.png。接第四期，经典视图中点击“多重来源”，并且可以调整多个文献的相对顺序。同时，如果选择同一个作者，可以做到一般要求的格式，只显示一次名字。

原创 【Zotero】【国标csl调教（四）】 author-date格式时作者在正文内 括号只显示年份

可以在首选项-引用-文档编辑软件里勾上“使用经典版添加引注对话框”，这样每次添加引注的时候就不要手动搜索或者切换视图了。在author date格式中，有时会出现需要让作者名在正文内，括号只显示年份的情况。添加引注默认的视图不是很好用，为了后面方便其他操作，改为默认采用经典视图。在添加引注时选择经典视图，右下角有个省略作者，勾上就好了。

原创 【Zotero】【国标csl调教（三）】author-date格式 显示et al的作者个数

et-al-use-first代表具体的展现形式，如输入1代表了满足et-al-min的条件之后，展现形式为第一个作者后面加上et al。和系列第一篇一样，注意author date格式里，citation和bibliography两处要修改。在正文引用（citation）或者参考文献表（bibliography）中，想限制显示的作者个数。et-al-min代表满足的条件，如输入3就代表了超过三个作者后使用et al。

原创 【Zotero】【国标csl调教（二）】参考文献表中出现添加日期，以及[J/OL] [D/OL]的修改方法

Zotero参考文献修改需要的步骤太多，站内有的教程就不总结了，为了以后换电脑或者numeric和author-date再弄一遍的时候，可以直接把以前配置好的复制过去。主要还是补充一下站内没有的以及一些写的不是很全面的教程。常见于一些原本是【J/OL 】或者[D/OL]的文献。，但是这个文章没有把引用日期的代码截全。如图，在参考文献表中出现日期。给这一段代码添加注释。

原创 【Zotero】【国标csl调教（一）】使用 GB/T 7714-2005 author-date格式时正文的英文文献“等”改为et al的方法

关于去除URL和DOI号，三人显示“等”或et al，作者名首字母大写等，站内有许多参考文章，如。

原创 【Arcgis】用高分辨率数据计算像元丰度的方法

关于像元丰度计算的方法，网上全部是基于高光谱数据利用端元法之类的，这里参考了一些帖子，总结一下利用高分辨率数据计算像元丰度的方法。此处的Field，工具会把Field相同的栅格转化为同一个面元素，这也是之前一定要制造各不相同栅格的原因。因为第四步的结果，出现了大量大于1089的统计结果，最大值（低分像元全部被高分像元填满）为1156。Convertion Tools，栅格转Point，Point再转栅格，这个栅格是各个点均不同的。首先确定参考栅格，即低分辨率数据，参考。6.计算像元丰度，栅格计算器。

原创 【Arcgis】Cell Statistics工具实现栅格像元累计识别

具体的需求描述是：有一系列栅格数据，经重分类后，0表示Nodata，1表示未识别，2表示识别，需统计这一系列栅格中哪些像元点被识别过至少一次。

原创 【ArcGis】巧用栅格计算器进行土地利用类型图的筛选&离散数据的重采样方法

有一幅30m分辨率的土地利用类型分类数据，现在要把研究数据（500m）按照这个土地利用类型裁剪，只保留耕地区域，如果采用栅格转矢量再裁剪，内存占用过大，容易死机。上述式子的意思是：栅格“clcd“中等于1的地方保持dtvdi的值，缺失的假语句对应的地区（clcd其他值）为NoData。函数的基本语句：con（条件（语句），真（语句），假（语句））

原创 【MATLAB】批量统计多幅栅格图像中像元数等于某一值的点的个数

一般用在重采样之后，批量重采样使用arcgis模型构建器。

原创 【arcgis】多幅栅格求和

原创 【arcgis】记录一下时间序列栅格图像动图的tips

2.日期格式注意，手动填写把0写上、如2016/03/03。1.添加栅格图像不要一次加过多，可以少量多次不然会报错。

原创 【GEE】使用expression之后导出数据由于范围过大导致分幅的问题及解决方案

关于gee输出是自动分割的问题，有很多文章介绍，主要就是研究区太大，这里我推断是expression的问题导致周围本来nodata的数据参与了运算和输出，不过与其刨根问底找expression的毛病，这里直接给一个治标不治本的方法。下载后在arcgis查看，发现两幅图像是一整副图像被分割而成，并且在研究区以外的nodata也都显示。在gee中使用expression进行栅格计算，输出到drive后发现一幅图像变成了两幅，如图所示。输出的时候再裁一遍就好了。

原创 【RS】MCTK工具处理MOD09GA数据--单波段导出和投影

MCTK工具对解决预处理MOD09GA数据

原创 【RS】使用wget批量下载MODIS数据的补充说明

对wget工具批量下载MODIS数据方法的补充

原创 【Origin】Excel连接时的日期格式问题

为了表示历年梅雨期起止日及长度的变化，可以采用origin中的浮动柱形图绘制，但是在把Excel文件连接到origin时出现问题用复制粘贴的方法没有办法将Excel的日期格式导入origin，origin工作簿中现实的是数字格式（类似43200这种），这时在将origin单元格格式改为日期，则直接显示–，意为格式错误如果先粘贴进数字格式，再选中列→右键属性→单元格格式→日期，这样则直接提示该列与格式不匹配综上可以得出结论，origin没有像Excel一样将数字格式直接转化为日期的功能当我采用连接到E

原创 【MATLAB】高塔观测数据波文比法计算感热及潜热通量

本文作为【MATLAB】高塔观测数据通量梯度法计算感热及潜热通量的后续，利用波文比法计算感热和潜热通量并作图数据导入参考上文此处我们忽略土壤热通量Qs代码Bowen.mT1 =Ta_20m_Avg;T2 =Ta_60m_Avg;%tempurature(C)Tmean=mean([T1,T2],2);theta1 =T1 +273;theta2 =T2 +273;%potential temperature(K)Rn=Rn_Avg;Qs=0;p =amb_press_Avg

原创 【MATLAB】高塔观测数据通量梯度法计算感热及潜热通量

数据导入利用MATLAB菜单栏的“导入数据”功能，将高塔观测数据导入要点Loggernet将原始二进制数据转为ASCII【data】【card convert】按文件夹选择MATLAB导入数据 输出类型中选择【列向量】导入65080 30min数据和9345flux的amb_press_Avg（气压）因为二者观测时间并不完全一致，因此在计算中应将时间取二者并集运行代码main.m%flux gradient method %author:Iris LeeM=96;%num.

原创 【MATLAB】分时段分类汇总代码

在微气象观测中，经常会需要对短时间内相同时段的数据做分类汇总。但是由于Excel并不支持隔项分类汇总功能。笔者写了一段小代码来实现这个功能。待处理数据5月到9月每天相同时段（0-23时）做汇总平均值代码S=table2array(data);导入名为data的数据，table转array是所有通过“导入数据”功能进行数据处理的必要步骤for t=0:23hour=S(:,3)';O3=S(:,4)';hour-数据中第三列为小时（分类项），并转置为行向量O3-数据中第四列为