计算CV变异系数，分析数据波动情况

%弧度转化为度分秒function angle=ArcToAngle(arc)arc0=arc*180/pi;deg1=fix(arc0);%取整min=fix((arc0-deg1)*60);sec=((arc0-deg1)*60-min)*60;angle=deg1+min/100+sec/10000;

自己学习编的程序，自己已全部运行，为了防止代码丢失，发个博客记录一下。一下代码比较多，每个代码块是一个函数，相互之间要调用，缺一不可。function M= GetM(B,a,e2)M=a*(1-e2)/sqrt(GetW(B,e2)^3);%计算法线长N，调用GetW函数function N=GetN(B,a,e2)N=a/GetW(B,e2);%计算辅助函数Wfunction...

var ndvi=image.normalizedDifference(["B5","B4"]);normalizedDifference([“B5”,“B4”])就代表(band5 - band4) / (band5 + band4)。也可编出计算公式直接计算。NDVI = (近红外波段 - 红波段) / (近红外波段 + 红波段)不同类型传感器NDVI计算的方式：Landsat8: NDVI

var roi = ee.FeatureCollection("users/FXS/xiangshan");//导入指定区域的矢量文件var image=ee.ImageCollection("LANDSAT/LC08/C01/T1_SR").filterDate("2020-10-10","2020-11-10")//筛选出指定时间内的影像.filterBounds(roi)

基于python的Hurst计算预测未来发展趋势（长时序栅格影像）

GDAL下载及安装

1.打开需要辐射定标及大气校正的图像2.进行辐射定标点击“基本工具——预处理——定标工具——Landsat定标”。之后在弹出的窗口选择带有多波段的图像，如：然后点击OK，在弹出的窗口设置参数，这些数值是图像自带的，尽量不要修改，每幅图的都不一样。记住这些数，下面还会用到。点击OK，定标完成。3.因为大气校正的图像格式是BIL格式，所以要转换该数据存储方式（BSQ转为BIL或BIP...

上时间序列栅格影像MK+Sen趋势检验

基于python的Hurst计算预测未来发展趋势（长时序栅格影像）