PIE-Engine学习笔记03-加载栅格数据
本次学习利用PIE-Engine Studio 加载栅格数据
一、加载Landsat8的影像数据
var img = pie.Image("LC08/01/T1_SR/LC08_119038_20170527")
.select(['B4','B3',"B2"]);
// print(img);
Map.centerObject(img,6);
Map.addLayer(img,{min:0,max:3000},'栅格图')
运行结果:
二、加载指定范围内的影像数据
1.指定范围
顺便来复习之前加载过的矢量数据,结合现在学的加载栅格影像,加载出北京市行政区内的栅格影像数据。
//加载北京市行政区数据
var roi = pie.FeatureCollection('NGCC/CHINA_PROVINCE_BOUNDARY')
.filter(pie.Filter.eq('name','北京市'))
.first()
.geometry();
print(roi);
Map.centerObject(roi, 6);
Map.addLayer(roi,{color:"FF0000",fillColor:"00000000",width:3})
运行结果:
2.根据条件筛选影像
//加载出北京行政区划内的Landsat8的栅格影像数据,根据矢量范围、日期、含云量等筛选
var imgCol = pie.ImageCollection("LC08/01/T1")
.filterBounds(roi)
.filterDate("2017-06-01","2017-10-01")
.filter(pie.Filter.lt('cloudCover',7));
print(imgCol);
// 根据波段、矢量范围进行裁剪
var img = imgCol.select(['B4','B3','B2'])
.mosaic()
.clip(roi);
print(img);
Map.centerObject(img,6);
Map.addLayer(img,{min:0,max:3000})
运行结果
3.举一反三,灵活运用
接下来就是PIE-Engine Studio 最方便的事情了,把省界行政区划一改就获取其他行政区的影像结果,具体代码如下:
//加载河北省行政区数据
var roi = pie.FeatureCollection('NGCC/CHINA_PROVINCE_BOUNDARY')
.filter(pie.Filter.eq('name','河北省'))
.first()
.geometry();
print(roi);
Map.centerObject(roi, 6);
Map.addLayer(roi,{color:"FF0000",fillColor:"00000000",width:3})
//加载出河北省行政区划内的Landsat8的栅格影像数据,根据矢量范围、日期、含云量等筛选
var imgCol = pie.ImageCollection("LC08/01/T1")
.filterBounds(roi)
.filterDate("2017-06-01","2017-10-01")
.filter(pie.Filter.lt('cloudCover',10));
print(imgCol);
// 根据波段、矢量范围进行裁剪
var img = imgCol.select(['B4','B3','B2'])
.mosaic()
.clip(roi);
print(img);
Map.centerObject(img,6);
Map.addLayer(img,{min:0,max:3000})
运行结果:
总结
本节主要学习了栅格影像的加载,同时根据日期、云量等条件去筛选影像,并依据根据行政区划进行裁剪。
欢迎大家批评指正,探讨交流!!!