- 博客(47)
- 收藏
- 关注
RSS订阅原创 【ENVI】SEBS模型计算蒸散发
输出的8个文件中包括波段1和2的地表反射率图像、波段4和5的地表温度图像、经纬度图像以及传感器的天顶角、方位角图像，通过单击可用波段列表（ AvailableBandsList ）中的" LoadBand "按钮可以查看相关图像。SEBS程序的目标就是在像元水平上求解陆面特征在瞬时状态下的潜热通量，SEBS中的所有功能都能通过SEBS主菜单来实现，其中包括“提取DEM”，“数据反演”，“地表物理参数”和“蒸发量”4个部分，各功能模块的顺序是根据SEBS在运行过程中对各项参数的要求进行设定的。
2023-05-11 21:06:06
1055
原创 Python和MATLAB和ARCGIS三种不同方法进行SEBAL模型计算蒸散发
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：GEE、ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享SEBAL模型介绍在 1998 年荷兰学者 Bastiaanssen 以能量平衡原理为基础，提出了利用遥感数据估算蒸散的模型 SEBAL。SEBAL（土地表面能量平衡算法）是一个由 25 个计算步骤基于能量平衡的原理通过遥感传感器获得的相关区域内可见光、近红外、热红外光谱辐射等数据反演地表反照率、地表温度等相关地表参数，在与所在地的气象数据及植被信息相结
2022-04-24 00:40:06
2479
3
原创 【Python】SEBS模型计算蒸散发
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：GEE、ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享SEBS模型介绍近年来出现了许多遥感反演蒸散发的模型，SEBS模型是由荷兰籍华人苏中波先生基于 SEBAL模型基础上发展的进行蒸散估算单层模型。该模型参数获取简单， 计算方法明确，且估算精度较高。可用于下垫面不统一且研究区范围大的区域。SEBS模型是一种综合利用遥感反演参数和研究区卫星过境时刻气象数据作为输入参数的模型。具体来说SEBS模型计算蒸散发需
2022-04-16 20:31:35
3609
14
原创 GEE作物时序变化监测
本研究使用哨兵二号数据，对三种不同作物进行计算NDVI，并对三种作物进行长时间NDVI变化监测，最终得到三种作物变化特征。教程照片代码及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：GEE、ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享...
2021-11-14 23:21:06
1652
原创 GEE行政区加载，高程坡度计算和裁剪
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：GEE、ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享01GEE加载地图影像Map.setCenter(116.39, 39.91)”02加载研究区var China_Provinces = ee.FeatureCollection("users");//加入研究区并重命名Map.centerObject(test1,4)//将显示中心调整为test1，缩放级别调整为4Map.addLay
2021-08-22 21:11:59
816
原创 GEE植被覆盖度计算
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享01植被覆盖度植被覆盖度指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比。植被的覆盖度可分为高、中高、中、低四种覆盖类型。覆盖度的检测主要有地表实测法和遥感监测法。植被覆盖度是衡量地表植被状况的一个重要指标，是描述生态系统的重要基础数据，也是区域生态系统环境变化的重要指示，对水文、生态、区域变化等都具有重要意义。使用GEE平台可以快速计算研
2021-08-16 13:15:20
3852
原创 集成学习-数学基础
01Metropolis-Hastings算法import numpy as npimport randomimport matplotlib.pyplot as plt## 设置参数mu = 0.5sigma = 0.1skip = 700 ## 设置收敛步长num = 100000 ##采样点数def Gussian(x,mu,sigma): return 1/(np.sqrt(2*np.pi)*sigma)*np.exp(-np.square((x-mu))/(2*
2021-07-13 00:08:15
289
4
原创 使用ENVI软件中的SARSCAP工具对RADARSAT数据进预处理
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享01SARSCAP设置SARscape 针对不同的数据源、不同的处理，提供了相应的系统默认参数，在处理数据之 前，根据数据选择一套相应的处理参数，参数包括所有的 SAR 数据处理步骤的参数，如：基本参数、滤波、地理编码、配准、去平、干涉、PS 处理，等在处理时，大多数情况都可按照 默认参数进行处理。SARscape 针对不同的数据源提供的系统参数有以下几套：
2021-03-20 21:01:41
2512
1
原创 使用Matlab软件对NDVI影像进行SG滤波处理
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享滤波是将信号中特定波段频率滤除的操作，是抑制和防止干扰的一项重要措施。在图像处理中，滤波是图像预处理的一种。图像处理中滤波将信号中特定的波段频率滤除，从而保留所需要的波段频率信号。根据选择保留的不同频段可以体现这么两个作用，消除图像中混入的噪声，对应的是低通滤波，噪声在图像中一般是高频信号。为图像识别抽取出图像特征，这里的特征一般为边缘纹理的特征，对应的是高通滤波
2021-02-04 01:10:04
2954
10
原创 Python从arxiv下载paper的示例代码
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享Arxiv是一个收集物理学、数学、计算机科学、生物学与数理经济学的论文预印本的网站，始于1991年8月14日。截至2008年10月，Arxiv.org已收集超过50万篇预印本；至2014年底，藏量达到1百万篇。截至2016年10月，提交率已达每月超过10,000篇。那么问题来了，预印版又是什么？预印本也称未定稿本。在学术出版领域，预印本是指尚未在需要同行评审的科学
2021-01-27 19:16:48
391
原创 使用Matlab软件对栅格数据进行相关分析及显著性检验
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享相关分析就是对总体中确实具有联系的标志进行分析，其主体是对总体中具有因果关系标志的分析。它是描述客观事物相互间关系的密切程度并用适当的统计指标表示出来的过程。显著性检验就是事先对总体（随机变量）的参数或总体分布形式做出一个假设，然后利用样本信息来判断这个假设是否合理，即判断总体的真实情况与原假设是否有显著性差异。今天的案例是2000-2015年的NPP数据和2
2021-01-23 00:26:28
6718
1
原创 使用Matlab软件对nc文件的数据产品处理教程
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享NetCDF(network Common Data Form)网络通用数据格式是一种面向数组型并适于网络共享的数据的描述和编码标准。目前，NetCDF广泛应用于大气科学、水文、海洋学、环境模拟、地球物理等诸多领域。用户可以借助多种方式方便地管理和操作 NetCDF 数据集。今天使用MATLAB软件对GIMMS的NDVI产品进行数据处理，得到tif格式的影像。
2021-01-22 00:25:25
3125
1
原创 使用ArcGIS软件进行投影变换和坐标校正
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享01坐标系、地图投影地球表面事物的定位采用二大类坐标：（1） 经纬度坐标，ArcGIS 称地理坐标系（2）二维笛卡尔平面坐标，ArcGIS 称投影坐标系。在实际工作中，经测量得到的空间信息在输入 GIS 数据库之前已经定好了坐标系。不同来源、不同坐标系的空间数据要在一起使用、相互参照时，就要作坐标转换，如果涉及不同的地图投影，要作投影变换。利用 ArcG
2021-01-20 20:04:10
8643
原创 使用Matlab软件进行逐像元Hurst指数分析
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享01理论知识Hurst指数：基于重标极差(R/S)分析方法基础上的赫斯特指数(H)的研究是由英国水文专家Hurst在研究尼罗河水库水流量和贮存能力的关系时，发现用有偏的随机游走(分形布朗运动)能够更好地描述水库的长期存贮能力，并在此基础上提出了用重标极差(R/S)分析方法来建立赫斯特指数(H)。作为判断时间序列数据遵从随机游走还是有偏的随机游走过程的指标。1
2021-01-19 21:06:23
4668
5
原创 使用Python对CMIP6数据批下载与裁剪
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享理论知识：CMIP是国际耦合模式比较计划（Coupled Model Intercomparison Project）的缩写，最早是在1995年由世界气候研究计划（WCRP）下属的耦合模式工作组（WGCM）主持开展的。自CMIP诞生以来，一直致力于促进气候模式的发展和完善，并支持气候变化的评估和预估工作。目前已开展了5次耦合模式比较计划，当前正在进行的是第6次耦
2021-01-19 21:03:41
3134
2
原创 使用Matlab软件对NDVI进行最大值合成
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享理论知识最大合成法：最大合成，就是将来多幅相同的栅格图叠加，每个栅格单元值取多幅中最大的那个，最后合成一幅。归一化植被指数【NDVI=(NIR-R)/(NIR+R)】1、检测植被生长状态、植被覆盖度和消除部分辐射误差。2、-1<=NDVI<=1，负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值，表
2021-01-19 21:00:52
6298
原创 使用python从ftp获取文件并下载到本地
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享文件传输协议（File Transfer Protocol，FTP）是用于在网络上进行文件传输的一套标准协议，它工作在 OSI 模型的第七层， TCP 模型的第四层， 即应用层， 使用 TCP 传输而不是 UDP， 保证客户与服务器之间的连接是可靠的， 而且是面向连接， 为数据传输提供可靠保证。FTP允许用户以文件操作的方式（如文件的增、删、改、查、传送等）与另一
2021-01-19 20:58:18
3604
3
原创 多年NDVI数据时间序列数据的MK检验
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享理论知识Mann-Kendall趋势检验法是世界气象组织推荐并已广泛应用的非参数统计方法，能有效区分某一自然过程是处于自然波动还是存在确定的变化趋势。对于非正态分布的水文气象数据，Mann-Kendall秩次相关检验具有更加突出的适用性。Mann-Kendall也经常用于气候变化影响下的降水、干旱频次趋势检测。Mann-Kendall非参数秩次检验在数据趋势检
2021-01-18 22:47:51
11106
10
原创 ENVI时间序列趋势分析之一元线性回归分析
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享课程作业中需要做一个多年NDVI趋势分析，以下内容仅提供简单的操作指导及大致思路流程。方便我以后再次使用时回忆。大神看到请回避，有错误也可以指出。理论知识：一元线性回归方程反映一个因变量与一个自变量之间的线性关系，当直线方程y=ax+b的a和b确定时，即为一元回归线性方程。经过相关分析后，在直角坐标系中将大量数据绘制成散点图，这些点不在一条直线上，但可以从中找
2021-01-18 22:45:08
9250
26
原创 Python读写HDF5文件
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享HDF5由美国伊利诺伊大学厄巴纳-香槟分校开发，是一种常见的跨平台数据储存文件，可以存储不同类型的图像和数码数据，并且可以在不同类型的机器上传输，同时还有统一处理这种文件格式的函数库。官网提供一个 JAVA 开发的 HDF5 数据可视化工具 HDFView，支持全平台查看数据, 但是注意打开文件的路径中不要包含中文。今天介绍一种python读取HDF5的文件的方
2021-01-18 22:43:06
2491
3
原创 ArcGIS软件气象数据插值教程
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享本次气象数据插值教程，针对的是txt气象数据，CMAD数据中有年份就用CMAD的，没有的就用气象数据的，另外气象数据和太阳辐射下只提取总辐射，一定要先看一下说明，气象数据下8个参数的台站都是一样的。动图跳跃分割线建立文件夹建立文件夹层次，风速文件夹-站点文件夹-该站点的所有时间段该参数值统计，由于参数文档按月份排，每一个文档下有所有站点的数据，所以我是将站
2021-01-18 00:02:51
6941
原创 ArcGIS软件CAD数据通过转换建立GIS数据库
教程照片及其他详细信息请关注微信公众号：夫也的笔记公众号内容包含：ArcGIS、ENVI、MATLAB、Python和R语言教程和实际案例分享。趁着寒假，系统的学习一下ArcGIS，我边学边写，你边看边学。动图跳跃分割线新建Geodatabase某地区有土地使用规划地图，要求输入 ArcGIS 数据库，原始图形是 AutoCAD的 DWG 文件，用 Polyline 和 Line 图形实体绘制地块边界，用 Text 实体注记的各地块编码（图 1）。图 1 用 AutoCAD 显示 landco
2021-01-18 00:01:27
2035
原创 ArcGIS最详细的地图制作教程
教程图片及其他详细介绍，请关注微信公众号：夫也的笔记01进入地图布局视图图 1 进入地图布局视图打开地图文档 Data\ex06\ex06.mxd，激活 Data Frame“人口密度”。选择菜单 View / Layout View，图形的显示由“数据视图（Data View）” 状态进入“布 局视图（Layout View）”状态；或者点击地图窗口右下角的第二个按钮也可，可通过第一和第二个按钮对地图数据视图和布局视图进行切换。此时，空间数据显示在一个新的地图布局里，当前的 Data Fr
2021-01-17 23:57:25
9987
原创 ArcGISModelBuilder创建与应用
实验过程照片及详细教程请关注微信公众号：夫也的笔记ModelBuilder模型构建器是一个用来创建、编辑和管理模型的应用程序。模型是将一系列地理处理工具串联在一起的工作流，它将其中一个工具的输出作为另一个工具的输入。也可以将模型构建器看成是用于构建工作流的可视化编程语言。模型构建器除了有助于构造和执行简单工作流外，还能通过创建模型并将其共享为工具来提供扩展 ArcGIS 功能的高级方法。01打开数据查看模型启动ArcMap不同输入来运行该模型，而无需每次都打开“模型构建器”。打开ArcTool
2021-01-17 23:54:55
1487
原创 基于matlab的栅格数据相关分析及显著性检验
图片详情及更加详细教程，请关注微信公众号：夫也的笔记相关分析在科研过程中常常会碰到，表征了一种两个因素间的相关程度，其值的大小只能说明相关性的强弱。比如当我们有2000-2015年的NPP数据和2000-2015年的降水数据时，我们想查看两者在空间上随时间变化的相关性。本文以产水和NPP为例进行说明具体代码如下：%将两者多年的数据放在三个不同的矩阵中nppsum=zeros(3587*4642,16);for year=2000:2015 filename=strcat('F:\课题项
2021-01-17 14:46:02
9739
20
原创 【Matlab】 气候资料数据集预处理
图片详情及更加详细教程，请关注微信公众号：夫也的笔记不啦不啦，今天介绍一种快速处理整理气象数据的方法，首次学起来很费劲，但是掌握以后对后面的各种任务、论文和项目应该会有帮助吧！先搞起来再说吧！声明一下啦！我不是大神，我只是从网上找到资料和教程以后，自己消化学习一下，然后发出来，当做自己的学习笔记方便回忆。文章收费是因为我也是花钱学来的，所以如果你也需要请为知识付费。本次主要是对中国气象数据网的V3.0地面气象观测数据集进行快速处理，将原先分散的txt各个要素整合到一个列表中，一行代表一个站点。直接从全国
2021-01-17 14:43:43
2152
8
原创 【Matlab】 多年NDVI数据的sen趋势分析
具体内容及详细教程请关注微信公众号：夫也的笔记最近老师在课上安排了一个汇报作业，要求对某一地区的多年NDVI数据进行趋势分析，对于小白的我来说，真是！唉！但是功夫不负有心人，在我充分利用百度搜索、谷歌搜索和Bing搜索后，终于让我找到了相关教程和知识。只是需要用Matlab软件进行一系列操作，废话不多说上代码。趋势分析有很多种，今天介绍的是Theil-Sen median趋势分析，一般还要结合Mann-Kendall建议，这个MK检验后续文章会写。Theil-Sen median趋势分析是一种稳健的非
2021-01-17 14:39:24
15255
19
原创 【ArcGIS】地表模型生成和显示
图片详情及更加详细教程，请关注微信公众号：夫也的笔记01由点状要素产生不规则三角网启动Data\ex13\ex13.mxd，激活 data frame1，可见到二个图层：线状图层“边界”和点状图层“高程点”，打开图层属性表“Attribute of 高程点”，可看到该图层有字段 Height 储存地形高程样本值，关闭属性表，返回 data frame1，选用菜单 View / Data Frame Properties，将 Map Units 和 DisplayUnits 均改为 Meters，返回
2021-01-17 14:35:49
2240
原创 数据分析-建模评价
学习内容模型搭建模型评估学习笔记1、模型搭建处理后的数据就是建模数据，下一步是选择合适模型。模型选择之前需要知道数据集最终是进行监督学习还是无监督学习。模型的选择一方面是通过我们的任务来决定的，另一方面可以根据数据样本量以及特征的稀疏性来决定，开始尝试使用一个基本的模型来作为其baseline，进而再训练其他模型做对比，最终选择泛化能力或性能比较好的模型。2、sklearn的算法选择路径3、模型评估模型评估是为了知道模型的泛化能力。交叉验证（cross-validation）是一种评
2020-08-27 23:07:48
380
原创 数据分析-数据可视化
学习内容Python数据可视化库Matplotlib通过菜鸟教程学习绘制不同图像学习笔记1、显示圆来代表点，使用 ob 作为 plot() 函数中的格式字符串。import numpy as np from matplotlib import pyplot as plt x = np.arange(1,11) y = 2 * x + 5 plt.title("Matplotlib demo") plt.xlabel("x axis caption") plt.ylabel("
2020-08-25 20:20:58
219
原创 数据分析-数据重构
学习内容数据合并数据聚合运算学习笔记1、stack函数是什么？函数原型：np.stack(array,axis,out=None)2、数据合并的方法（1）使用concat方法list_up = [text_left_up,text_right_up]result_up = pd.concat(list_up,axis=1)result_up.head()（2）使用DataFrame自带的方法join和append方法resul_up = text_left_up.join(te
2020-08-22 20:08:49
404
原创 数据分析-数据清洗特征提取
学习内容数据清洗特征提取学习笔记1、dropna和fillna有哪些参数，分别如何使用呢?1、使用dropna函数DataFrame.dropna(axis=0, how='any', thresh=None, subset=None, inplace=False)参数说明：axis:axis=0: 删除包含缺失值的行axis=1: 删除包含缺失值的列how: 与axis配合使用how=‘any’ :只要有缺失值出现，就删除该行货列how=‘all’: 所有的值都缺失，才删除行或
2020-08-20 20:58:51
932
原创 数据分析-基础分析
学习内容数据加载pandas基础探索性数据分析学习笔记1、‘read_csv’和‘read_table’’的区别？解答：read_csv()函数用于读取以’,'分割的文件到DataFrame。read_table()函数用于读取以’/t’分割的文件到DataFrame。2、tsv和csv的区别？TSV 是Tab-separated values的缩写，即制表符分隔值。CSV是Comma-separated values的缩写，即逗号分隔值。（1）从名称上即可知道，TSV是用制表符
2020-08-18 21:08:45
240
原创 Python进阶学习-文件系统
学习内容文件与文件系统文件对象方法简洁的 with 语句OS 模块中关于文件/目录常用的函数序列化与反序列化重点实例打开文件open(file, mode=‘r’, buffering=None,encoding=None,errors=None, newline=None, closefd=True)Open file and return a stream. Raise OSError upon failure.a. file : 必需，文件路径（相对或者绝对路径）。b. m
2020-08-08 18:59:36
123
原创 Python进阶学习-模块
学习内容什么是模块命名空间导入模块包（package）datetime模块date类time类timedelta类重点实例命名空间因为对象的不同，也有所区别，可以分为如下几种：内置命名空间（Built-in Namespaces）：Python 运行起来，它们就存在了。内置函数的命名空间都属于内置命名空间，所以，我们可以在任何程序中直接运行它们，比如id() ,不需要做什么操作，拿过来就直接使用了。全局命名空间（Module：Global Namespaces）：每个模块创建它
2020-08-06 20:37:03
122
原创 Python进阶学习-类对象魔法
学习内容对象 = 属性 + 方法self 是什么？Python 的魔法方法公有和私有继承组合类、类对象和实例对象什么是绑定？一些相关的内置函数基本的魔法方法算术运算符增量赋值运算符属性访问描述符定制序列迭代器重点实例类的方法与普通的函数只有一个特别的区别 —— 它们必须有一个额外的第一个参数名称（对应于该实例，即该对象本身），按照惯例它的名称是 self 。在调用方法时，我们无需明确提供与参数 self 相对应的参数。class Test:def prt(self
2020-08-04 13:44:26
243
原创 Python进阶学习-函数
学习内容函数的定义函数的调用函数文档函数参数函数的返回值变量作用域匿名函数的定义匿名函数的应用重点实例def add(a, b):print(a + b)add(1, 2) # 3add([1, 2, 3], [4, 5, 6]) # [1, 2, 3, 4, 5, 6]def MyFirstFunction(name):"函数定义过程中name是形参"# 因为Ta只是一个形式，表示占据一个参数位置print('传递进来的{0}叫做实参，因为Ta是具体的参数值！'.f
2020-08-02 13:57:57
205
原创 Python入门学习-字典集合序列
学习内容可变类型与不可变类型字典的定义字典的内置方法集合的创建访问集合中的值集合的内置方法集合的转换针对序列的内置函数重点实例那么如何快速判断一个数据类型 X 是不是可变类型的呢？两种方法：麻烦方法：用 id(X) 函数，对 X 进行某种操作，比较操作前后的 id ，如果不一样，则 X 不可变，如果一样，则 X 可变。便捷方法：用 hash(X) ，只要不报错，证明 X 可被哈希，即不可变，反过来不可被哈希，即可变。seq = ('name', 'age', 'sex')
2020-07-30 18:10:05
96
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人