Sentinel-1的一些数据介绍https://developers.google.com/earth-engine/tutorials/community/sar-basics
单独下载一张影像,影像信息如代码里的
// 导出HH HV 和angle
//下载指定文件名称的S1数据,导出为北极极投影,40m,双极化数据,波段为HH HV 入射角
var image = ee.Image.load('COPERNICUS/S1_GRD/S1A_EW_GRDM_1SDH_20160201T161634_20160201T161734_009755_00E419_7921');
// image = image.select("VH");
//转成float32位
image = image.float();
Map.centerObject(image,7);
Map.addLayer(image, {min: -40, max: 5},"HH", true);
// 对上述极化影像进行打包
var composite = ee.Image.cat([
image.select('HH'),
image.select('HV'),
image.select('angle')
]);
// region=ee.Geometry.BBox(image.getInfo())
//导出文件到Drive,40米分辨率
Export.image.toDrive({
image: composite,
description: "S1A_EW_GRDM_1SDH_20160201T161634_20160201T161734_009755_00E419_7921",
fileNamePrefix: "S1A_EW_GRDM_1SDH_20160201T161634_20160201T161734_009755_00E419_7921",
region: image.geometry().getInfo(),
crs: "EPSG:3413", //坐标系参数,极投影
scale: 40,
maxPixels: 1e13
});
print("over")
对影像进行一些bandmath计算,注意可以用上面的代码显示目前的波段命名,再进行重命名
var adddiff3 = function(image) {
return image.addBands(image.expression(
'HV - HH', {
'HV': image.select(['HV']),
'HH': image.select(['HH'])
}
));
};
PS:不知道为啥早上的下载网速很快,晚上很难下载下来,因此最好白天下载。