【GEE学习日记】栅格提取到点（根据经纬度提取栅格值）

一、背景
前面完成了MODIS的LST影像的加载和裁剪（按掩膜提取），接下来需要提取一些点上的LST值。
使用数据：
（1）MOD11A1
（2）研究区范围面矢量
（3）需要提取的点（可以现场写，也可以导入，后面会写出两种情况的代码）
二、参考文献
《Google Earth Engine（像元值提取）》
《GEE案例：根据点位提取栅格值（acrgis的多值提取到点）》
《二. GEE加载矢量数据》
三、代码
先import研究区矢量图

3.1 现场写2个点

//按掩膜提取函数
function ExtractByMask(Image)
{return Image.clip(roi)
}
//乘以参数函数
function Multiply(Image)
{return Image.multiply(0.02)
}var roi=table.geometry() //获得研究区
Map.centerObject(roi,5) //设置地图的显示中心
//选择数据集
var LST_dataset=ee.ImageCollection("MODIS/061/MOD11A1").filterDate('2002-01-01','2003-01-01').select('LST_Day_1km').filterBounds(roi);
print(LST_dataset)
var ft=ee.FeatureCollection(ee.List([])) //创建一个要素集合var p1 = ee.Geometry.Point([93.36, 32.23]) //第1个点的位置
var p2 = ee.Geometry.Point([93.23, 33.14]) //第2个点的位置
//把两个点合在一起组成要素合集
var pts = ee.FeatureCollection(ee.List([ee.Feature(p1).set('name','p1'),ee.Feature(p2).set('name','p2')]))
//设置栅格值提取到点函数
var fill = function(img, ini) {var inift = ee.FeatureCollection(ini)var ft2 = img.sampleRegions({collection:pts,properties:ee.List(['name']),//scale:30});var date = img.date().format()var ft3 = ft2.map(function(f){return f.set("date", date)})return inift.merge(ft3)
}
//先乘以参数再裁剪
var LST_Multiply=LST_dataset.map(Multiply)
var LST_ExtractByMask=LST_Multiply.map(ExtractByMask)
//在ImageCollection上迭代
var newft = ee.FeatureCollection(LST_dataset.iterate(fill, ft))
print(newft)
//导出表格
Export.table.toDrive({collection: newft,description: 'sample_get3',fileFormat: 'CSV'
});
//设置显示属性
var color={min:200,max:350,palette:['FFFFFF','CE7E45','DF923D','F1B555','FCD163','99B718','74A901','66A000','529400','3E8601','207401','056201','004C00','023B01','012E01','011D01','011301'],
};
Map.addLayer(LST_ExtractByMask,color,'LST') //添加LST图层
Map.addLayer(pts) //添加点
print('OK')

结果如下图所示：

打开云盘

3.2 导入点文件
我的研究区和所需要的点位置如下图所示：

先上传点矢量文件：

四、总结

原始数据集是365张影像，但由于只有193张影像在该点处有值，因此得到的表格只有193行，其余天数都是没有nodata

【GEE学习日记】栅格提取到点（根据经纬度提取栅格值）相关推荐

GEE学习记录~~~
目录 20220510--学习GEE平台上遥感信息的处理 20220511--GEE中上传tif文件 20221203--GEE计算导出辐射日均值csv 20230114--GEE计算多年辐射平均值导 ...
FireBird Embedded学习日记
FireBird Embedded学习日记 FireBird简介 Firebird是一个真正的关系数据库,支持存储过程.视图.触发器.事务等大型关系数据库的所有特性:不用考虑授权费用(免费),不用 ...
python学习日记2-3周
python学习日记2-3周怠惰怠惰 improt 导入 import phone.apple.iphone6 import phone.apple.iphone7 import phone.sam ...
10.16 my学习日记 (XPath的基础语法，lxml库的应用）
10.16 my学习日记(XPath的基础语法,lxml库的应用) XPath的基础语法 XPath查找标签 XPath谓语 lxml库在爬虫中的应用 etree库etree_Element对象使用 ...
gee学习笔记第二天
gee学习笔记第二天文章目录 gee学习笔记第二天前言一.继续熟悉jupyter lab notebook 二.获取和找到数据 1.选定区域和设置过滤器前言第二天学习如何准备水煮鱼的材 ...
【前端学习日记】利用reveal.js把实验报告做成一个简单的幻灯片
一.整体效果把电磁场的实验报告做成网页PPT,原文是这里:<[电磁场实验作业]有限差分法(FDM)求解静电场电位分布_轩辕衍的博客-CSDN博客> 二.核心代码讲解 0.创建页面第一篇 ...
java的圆周率_java学习日记，圆周率的打印
前段时间看到听说学习java每天写技术贴会对自己提升很大,我现在学习java也就2个周,算不上技术贴,就写写学习日记吧. 昨天师傅给我出了一道题,说是试试用java打印圆周率. 刚开始我的思路是,如果 ...
GPU（CUDA）学习日记（十一）------ 深入理解CUDA线程层次以及关于设置线程数的思考
GPU(CUDA)学习日记(十一)------ 深入理解CUDA线程层次以及关于设置线程数的思考标签: cuda存储线程结构网格 2012-12-07 16:30 6298人阅读评论(4)收藏举 ...
GPU（CUDA）学习日记（十三）------ CUDA内存简介
GPU(CUDA)学习日记(十三)------ CUDA内存简介标签: cuda存储线程结构 2012-12-07 16:53 2902人阅读评论(0)收藏举报分类: GPU(16) CUDA ...

【GEE学习日记】栅格提取到点（根据经纬度提取栅格值）

【GEE学习日记】栅格提取到点（根据经纬度提取栅格值）相关推荐

最新文章

热门文章