![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
遥感算法
文章平均质量分 53
_养乐多_
记录GEE(GoogleEarthEngine),python处理遥感数据,机器学习算法应用,深度学习模型训练、推理和部署的笔记。
展开
-
GEE:根据DEM绘制全球地形等高线
数据集:The Shuttle Radar Topography Mission (SRTM, see Farr et al. 2007) digital elevation data。本文记录了在GEE平台上绘制地形等高线的代码脚本。等高线的间隔可以设置。原创 2023-03-05 21:24:18 · 288 阅读 · 0 评论 -
GEE:排序影像(根据云量、时间)
本文记录了在影像集合里面选择云量最少的影像。主要用到了.sort(‘CLOUD_COVER’,true).first()函数。原创 2023-03-01 22:51:17 · 898 阅读 · 0 评论 -
GEE:绘制累积降雨量折线图、降雨量均值折线图
本文记录了在GEE(GoogleEarthEngine)平台上,使用UCSB-CHG/CHIRPS/DAILY数据集,绘制一个指定点的一段时间内的降水量累计值折线图的代码,也记录了绘制一个区域一段时间内的降水量均值的折线图代码。其中累积降水量使用每日数据累积,降水量均值折线图使用的是每周的降水量数据。即本文记录了将影像集合按天 ‘day’ 合成和按周 ‘week’ 合成的方法。操作对象:一个点研究区:一个点的缓冲区累计值:一个点一段时期内降水量的累积值二、降水量均值折线图绘制操作对象:整个影像研原创 2022-12-03 17:32:59 · 2065 阅读 · 0 评论 -
GEE:关系、条件和布尔运算
本文记录了在GEE(Google Earth Engine)平台上的关系运算符和布尔运算符,分别应用到了三个不用的场景(筛选低 NDVI 和低 NDWI 区域、根据阈值分类夜光遥感数据、使用一景影像无云像素替换另一景影像云像素),并展示了结果。关系运算函数包括:eq(),gt(),gte(),lt(),lte(),布尔运算函数包括:and(),or(),not()。这些方法可以用来掩膜、绘制分类地图和重新赋值。原创 2022-12-03 16:44:01 · 2446 阅读 · 0 评论 -
GEE:在研究区内随机生成指定数量点
本文记录了在GEE平台上指定实验区域内生成指定数量随机点的方法,展示了结果,记录了代码。var point= ee.FeatureCollection.randomPoints(roi, 17)原创 2022-12-02 16:17:01 · 1471 阅读 · 0 评论 -
GEE:筛选影像集合中的非空影像
【代码】GEE:筛选空影像。原创 2022-12-02 15:39:43 · 621 阅读 · 2 评论 -
GEE:提取一个小区域影像的所有波段像素值到csv
本文记录了一个将一个影像的所有像素的不同波段值和经纬度坐标保存到csv的代码。数据集合:COPERNICUS/S2_SR_HARMONIZED。原创 2022-11-24 21:51:15 · 1623 阅读 · 13 评论 -
GEE:MODIS计算生物量
本文记录在Google Earth Engine(GEE)平台上通过NPP和GPP数据来计算生物量。MODIS数据。原创 2022-11-23 23:18:48 · 2256 阅读 · 7 评论 -
GEE:提取多个点的时间序列数据,保存为csv格式(以NDVI为例)
本文第一个代码记录了提取一个点的LAI值时间序列的代码,并且记录了以csv格式保存结果到本地的方法。第二个代码记录了将一个点的时间序列以折线图的形式打印在控制台中。原创 2022-11-23 22:11:23 · 5404 阅读 · 50 评论 -
GEE:LST地表温度反演函数(针对Landsat8的T1影像和T1_L2影像集合)
LST 反演的代码,使用到了 Landsat-8 的 L1 级别数据和 L2 级别数据。因此,本文分别举例了针对 image 和 imageCollection 做 LST 反演的代码。数据集:LANDSAT/LC08/C02/T1_L2。操作波段:B10。数据集:LANDSAT/LC08/C01/T1。操作波段:ST_B10。原创 2022-11-08 17:50:11 · 5650 阅读 · 22 评论 -
GEE:Landsat-T1级别影像去云和辐射校正
landsatT1级别数据是需要做辐射校正的,看到很多人使用LandsatT1数据,记录一下T1数据做辐射校正的方法和使用的波段。原创 2022-11-08 16:12:08 · 1167 阅读 · 1 评论 -
GEE:将MODIS数据从500米分辨率重采样到30米以后,下载的像元是斜的
将MODIS数据从500米分辨率重采样到30米以后,下载的像元是斜的。本文提供了解决该问题的方法。原创 2022-10-28 15:43:13 · 3078 阅读 · 3 评论 -
GEE:DTW(Dynamic Time Warping)动态时间规整,Sentinel-2 时间序列分类
本文分享了外网上查到的 GEE平台上实现时间加权或时间约束动态时间规整 (TW/TC-DTW) 的模块,以及例子。原创 2022-07-10 21:38:28 · 3209 阅读 · 2 评论 -
GEE:利用libsvm()创建一个支持向量机SVM分类器进行土地利用分类
在GEE平台上,利用libsvm()函数进行土地利用分类原创 2022-07-07 15:02:46 · 2475 阅读 · 0 评论 -
python:使用xlrd、pandas、numpy做MK趋势分析
GEE上的MK趋势检验分析可以参考https://blog.csdn.net/qq_35591253/article/details/121617262一、MK趋势检验出图结果展示:输入格式:读取数据,提取年份和趋势分析数据:import numpy as npimport pandas as pdfrom matplotlib import pyplot as pltplt.rcParams['font.family'] = ['MicroSoft YaHei']plt.rcPa原创 2022-03-24 23:17:36 · 6655 阅读 · 6 评论 -
GEE:计算土地利用覆盖数据中的各个地类的面积/占比
计算2015年土地利用数据各个地类的面积方法一:面积被打印到了字典中。discrete_classification是土地利用类型的标签,sum是土地利用数据各个地类的面积(m^2)。怎么制作成时间序列的折线图,后续会补充。var roi = tablevar dataset = ee.Image("COPERNICUS/Landcover/100m/Proba-V-C3/Global/2015").select('discrete_classification').clip(roi);Ma原创 2021-11-15 22:04:41 · 7400 阅读 · 3 评论 -
python:nc文件转换为tiff格式(仿射变换)
import netCDF4 as ncfrom osgeo import gdal,osr,ogrimport numpy as npimport osfrom netCDF4 import Datasetimport globimport pandas as pd# 单个nc数据ndvi数据读取为多个tif文件,并将ndvi值化为-1-1之间# def NC_to_tiffs(data,Output_folder):# nc_data_obj = nc.Dataset(dat原创 2021-11-10 11:54:26 · 2604 阅读 · 2 评论 -
GEE:快速下载NO2、O3、SO2、 HCHO、CH4、CO和气溶胶数据(Sentinel-5P数据集)
Sentinel-5P数据集: Sentinel-5P是欧空局(Europe Space Agency, ESA)于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪(Tropospheric Monitoring Instrument,TROPOMI),可以有效的观测全球各地大气中痕量气体组分,包括NO2、O3、SO2、 HCHO、CH4、CO和气溶胶等重要的与人类活动密切相关的指标,加强了对气溶胶和云的观测。光谱仪:TROPOMI ,是目前世界上技术最先进、空间分辨率最高的大原创 2021-11-01 22:20:48 · 8572 阅读 · 30 评论 -
GEE:快速下载2020年土地利用/覆盖数据(10米分辨率)
数据集: ESA WorldCover 10m v100数据说明: 该产品提供了基于Sentinel-1和Sentinel-2数据的2020年10米分辨率的全球土地覆盖地图。WorldCover产品有11种土地覆盖类别,是在ESA WorldCover项目框架内生成的,该项目是欧洲航天局第5个地球观测信封计划(EOEP-5)的一部分。源代码链接: https://code.earthengine.google.com/a21637bdf6f94fe68d9eb7928ac91695?noload=tru原创 2021-10-28 00:15:02 · 7542 阅读 · 9 评论 -
GEE:降水量数据下载(每日/每月/每年/每五年)
数据集: CHIRPS Daily: Climate Hazards Group InfraRed Precipitation With Station Data (Version 2.0 Final)数据说明: Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) is a 30+ year quasi-global rainfall dataset. CHIRPS incorporates 0.05° resol原创 2021-10-27 23:40:43 · 45663 阅读 · 29 评论 -
GEE:遥感指数——NDFI(归一化退化指数)
NDFI是什么?NDFI被称为归一化退化指数。可以用来表达森林植被的退化程度和森林的健康程度。NDFI最开始发表的文章是(Combining spectral and spatial information to map canopy damage from selective logging and forest fires)。该指数在国际期刊中使用次数很多,使用该指数的文章多次发表在遥感顶刊RSE上,在国内少有见到相关文章(2021/10/20),因此我希望通过分享该指数可以使国内更多的人关注起来,原创 2021-10-20 12:44:17 · 3517 阅读 · 8 评论 -
GEE:生成20年NDVI动图和下载(彩色,使用MODIS数据集)
原代码链接https://code.earthengine.google.com/2d3796ec7be491f805fac4163acedd0a?noload=truevar roi = ee.Geometry.Polygon( [[112.51963551306746,26.815989026318352],[112.6974766995909,26.815989026318352],[112.6974766995909,26.931749679533954],[112.51原创 2020-10-23 23:37:53 · 6159 阅读 · 11 评论 -
GEE:Landsat影像快速下载
源代码:https://code.earthengine.google.com/df20aa8a78134144cfc5771d80014c1b?noload=truevar ROI = table.geometry();//引入一个变量,筛选时间。var data_collection = L8.filterDate('2020-06-01','2020-09-30');//用clip()函数里面的行政边界切割var ROI_collection = ee.Algorithms.Land.原创 2021-06-30 16:05:06 · 5335 阅读 · 4 评论 -
GEE:哨兵数据下载,将不同景按同一天日期合成
输入:输出:源代码:https://code.earthengine.google.com/dc1c75bdbeb215d4671db79862acbffb?noload=true如图所示,输入的影像集合里面有四幅影像,其中三幅是S2卫星同一天(5月7号)拍摄的,将同一天拍摄的s2数据按中值合成为一副。7月21号的合成为一副。最终输出两幅影像。//# 哨兵数据下载,将不同景按同一天日期合成 // date: 2021-06-03// author: Ruan xiaohao v原创 2021-06-03 21:59:55 · 4124 阅读 · 13 评论 -
GEE:快速下载数字高程DEM数据
DEM数据下载,使用SRTM Digital Elevation Data Version 4,一个2000年 的全球高程数据 。下载以后的DEM不用再做裁剪,可以直接使用。源代码:https://code.earthengine.google.com/f8e39c6c96a935499b09595e2fd83b28?noload=true效果如下图所示,不光可以下载DEM数据,也可以下载坡向。var roi = tablevar roi_geometry = roi.geometry()原创 2021-04-28 09:51:51 · 14887 阅读 · 16 评论 -
GEE:批量下载NDVI数据,按年最大值合成,GIMMS NDVI from AVHRR Sensors (3rd Generation)
NDVI效果如下图所示:所用数据集为:GIMMS NDVI from AVHRR Sensors (3rd Generation)这个数据集有1981年到2013年每个月的NDVI数据。原创 2021-04-23 17:18:46 · 8514 阅读 · 6 评论 -
matlab:自组织映射(SOM)神经网络
自组织映射网络在维基百科上讲的很清楚https://en.wikipedia.org/wiki/Self-organizing_map中文网站这位大佬解释的很清楚:https://blog.csdn.net/wj176623/article/details/52526617matlab实现SOM的代码如下:%% SOM神经网络的数据分类--柴油机故障诊断%% 清空环境变量clcclear%% 录入输入数据% 载入数据load shiyandatda1990.mat;%转置后符合原创 2021-04-08 15:14:56 · 5517 阅读 · 8 评论 -
GEE:mosaic、qualityMosaic、cat、addBands的区别
如图所示,mosaic()是按有值的最后一个值合成qualityMosaic()是按最大值合成操作对象都是影像集合cat()是将几幅影像拼接起来,连接起来区别:cat()是Concatenate;mosiaic()是CompositesaddBands()是给一副影像添加波段,操作对象是一个影像。上面图片展示的两个例子结果相同,但是方法不同。...原创 2021-04-06 23:02:47 · 3903 阅读 · 0 评论 -
《写论文就像做饭》
发表论文=文献调研-论文写作-论文修改-投稿检索——分析——发现——卖点但是首先我们需要确定,我们发表我们期刊,有的期刊侧重理论,有的期刊侧重实际应用,我们需要写出来符合期刊口味的论文,才好发表。发表论文实在就是比较论文,都有人投稿这个期刊,为什么期刊编辑,为什么录用你的文章,不用别人的文章,如果可以回到好这个问题,说明差不多可以发表了。发表论文之前,我自己几个问题1、为什么我这个论文可以发表到这个期刊上面,和已经发表到这个期刊上面的论文也什么区别。2、我做了什么工作3、我为什么要做这个工作原创 2021-04-02 15:49:56 · 666 阅读 · 3 评论 -
遥感数据下载网站汇总
1、GEE(Google Earth Engine)云平台2、美国地质调查局:查影像号https://earthexplorer.usgs.gov/和下载https://espa.cr.usgs.gov/index/?next=%2Fordering%2Fnew%2F%3F3、https://ladsweb.modaps.eosdis.nasa.gov/4、Earth Observation Group (EOG):https://www.ngdc.noaa.gov/eog/5、地理空间数据云:ht原创 2021-04-02 15:45:21 · 1507 阅读 · 0 评论 -
GEE:对矢量数据/行政边界设置填充颜色为空,只显示边界
//设置显示样式:color代表边界颜色;fillcolor代表填充颜色var styling={color:'red',fillColor:'00000000'}//JX就是矢量边界Map.addLayer(JX.style(styling),{},"JX")矢量化效果原创 2020-12-30 20:37:22 · 4969 阅读 · 7 评论 -
GEE:LandTrendr时间序列曲线拟合
毕竟是时间序列分析,LandTrendr的时间序列曲线是怎么拟合的?在一个像素点的位置上,它是怎么输出变化年份、扰动幅度、恢复幅度和时间的?以下代码给出了回答,可以参考LT-GEE官方的说明文档。用LandTrendr选取指定像素点的APP:https://emaprlab.users.earthengine.app/view/lt-gee-pixel-time-seriesLT-GEE的结果包括:每个像素时间序列的观测年份;二维频谱时间空间中的x轴值;(默认)每个像素时间序列的观测值的源值;二原创 2021-03-30 19:46:46 · 6911 阅读 · 18 评论 -
GEE:变化检测算法,干扰检测,LandTrendr代码与说明(变化年份和扰动幅度示例)
LandTrendr是什么?LandTrendr是一组频谱时间分割算法,可用于中等分辨率卫星图像(主要是Landsat)的时间序列中的变化检测以及生成很大程度上不存在年际信号噪声的基于轨迹的频谱时间序列数据。LT最初是用IDL(交互式数据语言)实现的,但是在Google工程师的帮助下,它已移植到GEE平台。LandTrendr官方说明网站:https://emapr.github.io/LT-GEE/index.html该网站介绍了在Google Earth Engine(GEE)中运行LandT原创 2021-03-30 19:11:31 · 13891 阅读 · 45 评论 -
GEE:批量下载数据,按月合成,不同景,全实验区域,一键执行
有时候我们下载的数据可能因为遥感影像采集时间不同,下载的数据需要拼接,这样很麻烦。所以在这里给大家一个GEE代码,可以直接下载按月份拼接的遥感影像。比如,将2021年3月的不同景的数据按最大值合成且拼接在一起。此代码可以下载2014年-2020年每个月的数据,且按照矢量边界裁剪。方便了下载数据的繁琐过程。当然也可能有一些莫名的错误,希望大家指正,欢迎讨论。(此代码没有去云,读者可以根据需要自行加入)结果如图所示源代码:https://code.earthengine.google.com/0b原创 2021-03-29 15:03:39 · 9548 阅读 · 15 评论 -
GEE:如何直接下载指定投影坐标系的遥感影像?(附投影坐标系查询网站)
我们在处理遥感影像时,经常会出现一个问题,就是数据的投影坐标系不是理想的坐标系。然而在ArcGIS中重新转换坐标系,步骤麻烦繁琐,时而伴随莫名其妙的问题。所以,我们可以在下载数据时,直接下载指定坐标系的数据。主要操作就是在toDrive()函数里面加入crs:"EPSG:32649"Export.image.toDrive({ image: crude_object_removal,//指定的数据名称 description: 'RF2020', //下载的数据名 scale: 3原创 2021-01-21 01:38:10 · 5502 阅读 · 2 评论 -
GEE:变化检测算法介绍,和代码(Landtrendr,Bfast,CCDC,CODED,多时间特征叠加的简图分类方法,CUSUM)遥感时间序列
Landtrendr的原理图原创 2020-10-20 18:27:06 · 17350 阅读 · 3 评论 -
遥感算法:遥感算法总结
遥感数据预处理(去云,辐射校正,合并,切割,大气校正,几何校正)遥感影像分类(线性,阈值,机器学习,深度学习)遥感时间序列(Landtrendr,Bfast,VCT,CUSUM,CCDC)原创 2020-10-18 23:23:27 · 4367 阅读 · 2 评论