Never Settle

在GEE中实现时，发现简单的QA60是无法去云的，最近S2地表反射率数据集又进行了更新，原有的属性集也进行了变化，现在的SR数据集名称是“S2_SR_HARMONIZED”。// 云概率低于 20% 保留。SCL在最新的S2_SR_HARMONIZED产品中是存在的，判断云和云阴影。继续用像元云概率进行判断，要求每个像元是云的概率小于10%要想得到研究区无云的图像，可以参考执行以下几个步骤，遥感图像中筛选云占比不超过10%的图像。// 第一次云掩膜：使用QA60。// 去除云(8)和云阴影(9)

叶绿素分子吸收光能(激发能)后，由基态跃迁到激发态(图中的第一单线态和第二单线态)，激发态是不稳定的状态，会重新回到基态，电子由激发态回到基态的过程中，大部分能量转向反应中心推动光化学反应以及接下来的电子传递、光合磷酸化，固定、还原CO。并终将能量贮存在有机物中;一部分能量以热的形式耗散;再有一小部分能量以荧光的形式发出。三者之间是此消彼长相互竞争的关系。因此我们可以用叶绿素荧光来研究光合作用的变化。

Poullain等于2012年提出Poullain模型，用一种朗伯余弦定律Lambertian cosine law和贝克曼散射定律Beckmann scattering law的组合（简称Poullain模型）去表征海岸带表面的后向散射强度分布。Kaasalainen等2018年将Poullain模型用于8波长的多光谱激光雷达植被叶片入射角效应研究中，将粗糙度因子m视为一个广义粗糙度因子，计算出了不同波长下的Poullain模型中参数分布。受Rees著《遥感物理》一书中的启发，

1.微波遥感在植被冠层郁密时信号容易到达饱和，限制了其在高度估算上的应用。2.由于卫星提供能量、大气扰动和激光传感器的传播距离等因素的影响，星载激光雷达系统的激光扫描频率和信噪比较低，利用LM（列文伯格-马夸特）算法进行波形分解容易陷入局部最优解。我评：其实LM算法只是非线性最小二乘法，并不是波形分解算法，分解模型是本质函数决定，LM算法知识求最优解的一种算法，属于求解解析解的优化解问题。...

本实验中给定三棵树和一株草，树高分别为20m，10m，2m（树干1m），草高0.5m。1.需要产生对应不同的叶绿素含量，使用matlab中的随机函数randperm(n,k)p = randperm(n,k) 返回一行从1到n的整数中的k个，而且这k个数也是不相同的。randperm(80,4)ans =41 34 69 80将其记录并赋值给上述四株植物，并执行rapid高光谱激光雷达波形模拟。...

1.RAPID中方格内导入并显示场景时场景文件必须为.map后缀，in不可以。2.RGMDISP.exe中导入场景时需要导入poly.in，后缀为in的文件，方可显示。3.需要运行rapid模型时，.map文件为非必须文件，in文件为必须文件。3.1在软件中手动运行模型时，需要先创建.map文件，然后创建对应的in文件，设定一下太阳和观测角度，然后即可开始运行Go—64bit模型。3.2使用代码创建场景并运行模型时，先使用代码直接创建in文件时，将省略手动创建.map文件的过程，然后代码中设定

RAPID中生化参数分别设置有：1.树干水分含量%2.树枝含水量%3.叶片含水量%4.叶片叶绿素含量mg/m25.比叶面积SLA（m2/kg）比叶面积（specific leaf area, SLA）：叶的单面面积与其干重之比。在同一个体或群落内，一般受光越弱比叶面积越大，所以比叶面积可作为叶遮荫度的指数而使用。叶面积是与产量关系最密切、变化最大，同时又是比较容易控制的一个因素，许多增产措施，包括合理密度和合理肥水技术之所以有显著的增产作用，主要在于适当地扩大了叶面积。叶面积小不好，但

光谱响应函数（Spectral Response Function，SRF），又称光谱响应曲线，是传感器固有的系统参数。对于某束太阳光，其经大气，入射到物体表面后，又沿某一方向反射到传感器入瞳处的能量L是固定的，对于传感器S1如果某一波段a1nm的光谱响应函数为0.9，则入射进传感器的能量则为0.9L，另一传感器S2的光谱响应函数为0.8，则入射进传感器的能量为0.8L。此外，通常计算传感器大气顶层波段平均太阳辐照度时，也用到SRF。由于是波段内平均，所以波长从a到b需要乘以太阳光谱曲线（sun spec

实在是太忙，没时间整理，没办法，抽时间整理一下吧。FY-4A/LPW水汽产品4km经纬度查找表生成代码：# -*- coding: utf-8 -*-"""Created on Tue Apr 23 11:39:02 2019@author: Administrator"""import osimport numpy as np import mathimport netCDF4 as ncimport h5py from osgeo import gdalimport op

pieces of the leaf were immediately cut using a hole punch with a diameter of 6 mm.在叶片表面切割一个直径6毫米孔大小的部分。The pieces were ground in a clear mortar with a pestle, and then the leaf pigment mixture wasplaced in a 25‐ml volumetric flask with 95% ethanol.将其.

分三步：1.定义原始数据的投影，无特殊说明，一般是认为基于WGS-84基准椭球体下进行UTM投影。（怀来试验站中央经线为115度，（115+180）/6 = 50N）注意带号。2.定义新的坐标系。选择地理坐标系，因为输出结果只要经纬度。3.数据框属性设置地图显示以度分秒为单位，而非米。结束！详细操作见视频。有用的话记得回来点个赞！举手之劳，不胜感激！...

亟人之所需，你所要的，这里全都有。欢迎关注，打扰！个人博客：https://blog.youkuaiyun.com/qq_37970770

我们知道简单散射实体一次散射的雷达方程原始公式为：Pr=PtD24πR4βt2σηatmηsysP_{r}=\frac{P_{t}D^{2}}{4\pi R^{4}\beta_{t}^{2}}\sigma\eta_{atm}\eta_{sys}Pr​=4πR4βt2​Pt​D2​σηatm​ηsys​（1）式中，PrP_{r}Pr​是传感器接收功率，DDD为接收孔径，RRR为目标与传感器之间的距离，βt\beta_{t}βt​为波束发散角，是立体角，也称为波束宽；σ\sigmaσ为目标物的后向散射截面。

dIλkλρds=−Iλ\frac{dI_{ \lambda }}{k_{\lambda}\rho ds}=-I_{\lambda}kλ​ρdsdIλ​​=−Iλ​以上，就是在忽略多次散射和发射的增量贡献时的比尔-布格-朗伯定律，kλk_{\lambda}kλ​为质量消光截面，为每单位质量的消光面积，为质量吸收截面和质量散射截面之和，单位是 cm−1cm^{-1}cm−1. ρ\rhoρ为物质密度，二者之和为消光系数。如果单个粒子的消光截面表示为σe\sigma_{e}σe​，单位体积的粒子数密度为NN

不矛盾。这是本科学习知识的一个误区，属于知识没有掌握牢固到位部分，正如很多人现在依然分不清普朗克定律的单位是辐射通量密度，单位是W/m2/um，而不是一个模糊的辐射，模糊的通量，亦或者不知道就模糊说为辐射亮度，学遥感的人，上述物理量的单位必须早已烂熟于心。漫反射体，本科学的时候都知道是各向同性，而我们同学又可能在别的地方看到了：朗伯体遵循余弦定律。那问题来了，对于入射的一束光，经过漫反射体反射，准备出射时，既然已经是漫反射体，出射时各个方向应该是一样的，各项同性的，何来余弦定律之说呢？注意了，此处就是

摘要关于很多matlab矢量裁剪栅格的内容，本文也来简要说说如何使用矢量裁剪栅格。思想（重要）利用矢量裁剪栅格，我没有找到什么好用的函数，其核心思想目前我仍然是取值赋值，先找到矢量边界图形的boundary range，即经度和纬度最大最小值，然后利用栅格图像中的boundary range中的x_min和y_max，以及图像的空间分辨率 ° ，对矢量边界图形的经度和纬度最大最小值求其行列号。以上过程可以得到矢量边界图形对应栅格中的起始行列号和终止行列号，裁剪的图像的本质就是找到其下标范围，对栅格

pi = 3.1415926535898普朗克常数h=6.62607015* 10.0^(-34) J·sc = 2.99792458 * 10.0^8 m/slamda: um波尔茨曼常数 k=1.380649 * 10^(-23) J/K逆函数中：按照辐射通量密度反算亮温时：c1 = 2 * pi * h * c ^2= 3.841 * 10.0^8 W/m2...

题目：静止卫星遥感图像太阳及卫星天顶、方位角(SOZ/SOA/SAZ/SAA)计算方法作者：白艺亭日期：2020.02.25摘要遥感图像角度数据（太阳天顶角、太阳方位角、卫星天顶角、卫星方位角）是遥感图像基础数据的重要组成部分，在遥感图像预处理和遥感应用如气溶胶光学厚度反演中是必不可少的参数之一。国产遥感数据经常缺失角度数据支持，大多数产品都不给出成像时对应的逐像元角度数据，需要用户手动计...

数据：GLASS12B02.V02.A2009001.2019093.hdf语言：IDL任务：读取1982-2017文件夹内的所有8天时间分辨率的GPP遥感图像数据，并找出地面站点经纬度处的像元值，用txt存储，第0-2列为地面站点txt原来的年份，纬度和经度，第3-48列存储遥感图像值。如图：站点txt文件内容：代码如下pro read_GLASS_GPP time_0 = ...

作业一请利用下列数据统计各类土地利用类型的NDVI均值。NDVI值通常位于-1至1之间，请将该地区NDVI异常值发生的像元位置和具体数值，生成一个点shape文件。数据：(1) Landsat8地面反射率文件LC08_L1TP_121041_121042_20190312_T1_sr_scale_mosaic_clip.tifLandsat8的NDVI公式：NDVI = (Band5-B...

envi/idl 批处理 释放 内存

最近很忙，每天都特别忙碌，也很少有时间去打球了。想念南京的老铁们，你们还好吗~过得可好，我在北京，想念你们的第53天。祝愿大家一切都好，北京一切还好，就是北京物价很贵，没有工资补助的我露出了卑微的笑容~今天和大家分享的是IDL怎么读取MOD021KM MYD021KM中的角度数据，将其重采样到其他的反射率光谱数据的维数大小，即ns = 1354，nl = 2030.重要的是如何避免最容易错的...

摘要：笔者主要花了半个小时研究了一下AERONET中国站点分布图。思路是先从AERONET官网，点击亚洲中国区域，将所有站点及经纬度复制下来，然后制作成点图层，最后检索确定中国区域的2级站点图层。一、AERONET亚洲中国区域站点经纬度打开AERONET网站：https://aeronet.gsfc.nasa.gov/ ，点击左侧Aerosol Optical Depth， 结果显示，如下图：...

本文研究了MODIS遥感图像的4种下载方式，该方式包括单景下载或批量下载方式，并同时比较了不同方式的下载速度，研究表明：单景下载时，选择LAADS DAAC一个一个点击的方式下载，还算比较快，但也不算太快，网络也不太稳定。批量下载时通过IDM方式下载最优，速度最快，网络极其稳定。需要指出的是，对于不同学校不同通信商所辖的校园网宽带及其他网络连接情况下的下载速度可能有不同，有待进一步研究。

自从数据发布以来局部中国区域文件官方只给出了图像文件，并未给出对应的经纬度查找表，而对于只研究中国区域的同学来说，下载全圆盘图像文件显得费时费力，然后处理起来又占据内存。为了对中国区域图像文件进行几何校正，将标称投影转为等经纬度投影，博主和同学老王进行了交流，终成此文。下面将利用Python实现区域图像文件的经纬度查找表建立！~将中国区域行列号转为经纬度查找表，并且存为tiff文件，band1...