forked from terryyangwhu/BH_China
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSentinel_SAR_China.txt
95 lines (84 loc) · 3.39 KB
/
Sentinel_SAR_China.txt
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
//collecting Sentinel-1 SAR subset (two polarization) across China in winter season of 2017
//import boundary
var BD = ee.FeatureCollection("users/cyangpkucues/shp_files/CN-sheng-A");
Map.addLayer(BD);
//Map.centerObject(BD,10);
print(BD);
//---------------------------import S1 VV polarization----------------------------
var s1vv1 = ee.ImageCollection('COPERNICUS/S1_GRD')
.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
.filter(ee.Filter.eq('instrumentMode', 'IW'))
//.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'))
//.filter(ee.Filter.eq('relativeOrbitNumber_start','69'))
//.filter(ee.Filter.eq('relativeOrbitNumber_stop','69'))
.select('VV')
.filterBounds(BD)
.filterDate('2017-01-01','2017-03-31');
//print(s1vv1);
var s1vv1_avg = s1vv1.median().clip(BD);
print(s1vv1_avg);
var s1vv2 = ee.ImageCollection('COPERNICUS/S1_GRD')
.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VV'))
.filter(ee.Filter.eq('instrumentMode', 'IW'))
//.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'))
//.filter(ee.Filter.eq('relativeOrbitNumber_start','69'))
//.filter(ee.Filter.eq('relativeOrbitNumber_stop','69'))
.select('VV')
.filterBounds(BD)
.filterDate('2017-11-01','2017-12-31');
//print(s1vv2);
var s1vv2_avg = s1vv2.median().clip(BD);
print(s1vv2_avg);
var s1vv = s1vv1_avg.addBands(s1vv2_avg);
var s1vv_avg = s1vv1.median().clip(BD);
print(s1vv);
print(s1vv_avg);
Map.addLayer(s1vv_avg);
Export.image.toDrive({
image: s1vv_avg.select('VV'),
description: "s1vv_avg_2017_CN",
fileNamePrefix: "s1vv_avg_2017_CN",
scale: 10,
region: BD,
crs: "EPSG:4326",
maxPixels: 1e13
});
//-------------------------import S1 VH polarization--------------------------
var s1vh1 = ee.ImageCollection('COPERNICUS/S1_GRD')
.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH'))
.filter(ee.Filter.eq('instrumentMode', 'IW'))
//.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'))
//.filter(ee.Filter.eq('relativeOrbitNumber_start','69'))
//.filter(ee.Filter.eq('relativeOrbitNumber_stop','69'))
.select('VH')
.filterBounds(BD)
.filterDate('2017-01-01','2017-03-31');
//print(s1vv1);
var s1vh1_avg = s1vh1.median().clip(BD);
print(s1vh1_avg);
var s1vh2 = ee.ImageCollection('COPERNICUS/S1_GRD')
.filter(ee.Filter.listContains('transmitterReceiverPolarisation', 'VH'))
.filter(ee.Filter.eq('instrumentMode', 'IW'))
//.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'))
//.filter(ee.Filter.eq('relativeOrbitNumber_start','69'))
//.filter(ee.Filter.eq('relativeOrbitNumber_stop','69'))
.select('VH')
.filterBounds(BD)
.filterDate('2017-11-01','2017-12-31');
//print(s1vv2);
var s1vh2_avg = s1vh2.median().clip(BD);
print(s1vh2_avg);
var s1vh = s1vh1_avg.addBands(s1vh2_avg);
var s1vh_avg = s1vh1.median().clip(BD);
print(s1vh);
print(s1vh_avg);
Map.addLayer(s1vh_avg);
Export.image.toDrive({
image: s1vh_avg.select('VH'),
description: "s1vh_avg_2017_CN",
fileNamePrefix: "s1vh_avg_2017_CN",
scale: 10,
region: BD,
crs: "EPSG:4326",
maxPixels: 1e13
});