質問編集履歴
1
エラーの出る前までの全文を載せました。
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -36,6 +36,164 @@
|
|
|
36
36
|
|
|
37
37
|
|
|
38
38
|
|
|
39
|
+
追記です。10/7 18:58
|
|
40
|
+
|
|
41
|
+
エラーが出る最後の行まで全文を載せます。これでわかるでしょうか。
|
|
42
|
+
|
|
43
|
+
|
|
44
|
+
|
|
45
|
+
min_lat =35.357416
|
|
46
|
+
|
|
47
|
+
min_lon =138.724022
|
|
48
|
+
|
|
49
|
+
max_lat =35.369735
|
|
50
|
+
|
|
51
|
+
max_lon =138.737926
|
|
52
|
+
|
|
53
|
+
|
|
54
|
+
|
|
55
|
+
# AVNIR-2のURL\n",
|
|
56
|
+
|
|
57
|
+
url_AVNIR2 = 'https://gisapi.tellusxdp.com/api/v1/av2ori/scene'
|
|
58
|
+
|
|
59
|
+
# LandsatのURL\n",
|
|
60
|
+
|
|
61
|
+
url_Landsat8 = 'https://gisapi.tellusxdp.com/api/v1/landsat8/scene'
|
|
62
|
+
|
|
63
|
+
headers = { 'Authorization': 'Bearer ' + API_TOKEN }
|
|
64
|
+
|
|
65
|
+
# パラメータ(緯度、経度の最大値、最小値)を指定\n",
|
|
66
|
+
|
|
67
|
+
parameter = { 'min_lat': min_lat, 'min_lon': min_lon,
|
|
68
|
+
|
|
69
|
+
'max_lat': max_lat, 'max_lon': max_lon }
|
|
70
|
+
|
|
71
|
+
|
|
72
|
+
|
|
73
|
+
# AVNIR-2の検索\n",
|
|
74
|
+
|
|
75
|
+
r_AVNIR2 = requests.get(url_AVNIR2,
|
|
76
|
+
|
|
77
|
+
params=parameter,
|
|
78
|
+
|
|
79
|
+
headers=headers)
|
|
80
|
+
|
|
81
|
+
|
|
82
|
+
|
|
83
|
+
# Landsat-8の検索\n",
|
|
84
|
+
|
|
85
|
+
r_Landsat8 = requests.get(url_Landsat8,
|
|
86
|
+
|
|
87
|
+
params=parameter,
|
|
88
|
+
|
|
89
|
+
headers=headers)
|
|
90
|
+
|
|
91
|
+
|
|
92
|
+
|
|
93
|
+
# AVNIR-2の検索結果を表示\n",
|
|
94
|
+
|
|
95
|
+
print("AVNIR-2のリスト [撮影日時、rspId、productId]")
|
|
96
|
+
|
|
97
|
+
for res in r_AVNIR2.json():
|
|
98
|
+
|
|
99
|
+
print(",".join([res['acquisitionDate'],
|
|
100
|
+
|
|
101
|
+
res['rspId'], res["productId"]]))
|
|
102
|
+
|
|
103
|
+
|
|
104
|
+
|
|
105
|
+
# Landsat-8の検索結果を表示\n",
|
|
106
|
+
|
|
107
|
+
print("Landsat-8のリスト [撮影日時、path、row、productId]")
|
|
108
|
+
|
|
109
|
+
for res in r_Landsat8.json():
|
|
110
|
+
|
|
111
|
+
print(", ".join([res['acquisitionDate'],
|
|
112
|
+
|
|
113
|
+
str(res["path"]), str(res["row"]),res["productId"]]))
|
|
114
|
+
|
|
115
|
+
|
|
116
|
+
|
|
117
|
+
# 選択したAVNIR-2衛星画像の productId、rspIdを格納\n",
|
|
118
|
+
|
|
119
|
+
AVNIR2_rspId = 'D070P3'
|
|
120
|
+
|
|
121
|
+
AVNIR2_productId = 'ALAV2A097002890'
|
|
122
|
+
|
|
123
|
+
|
|
124
|
+
|
|
125
|
+
# 選択したLandsat-8衛星画像の productId、path、rowを格納\n",
|
|
126
|
+
|
|
127
|
+
Landsat8_path = '108'
|
|
128
|
+
|
|
129
|
+
Landsat8_row = '35'
|
|
130
|
+
|
|
131
|
+
Landsat8_productId = 'LC08_L1TP_108035_20181125_20181210_01_T1'
|
|
132
|
+
|
|
133
|
+
|
|
134
|
+
|
|
135
|
+
# 緯度経度をタイル座標に変換する関数を定義\n",
|
|
136
|
+
|
|
137
|
+
def LatLon2Tile(lat,lon, z):
|
|
138
|
+
|
|
139
|
+
L= 85.05112878
|
|
140
|
+
|
|
141
|
+
tileX = int( (2**(z+7)) * ((lon/180) + 1) )
|
|
142
|
+
|
|
143
|
+
|
|
144
|
+
|
|
145
|
+
tileY = int( (2**(z+7)/math.pi) *
|
|
146
|
+
|
|
147
|
+
(-1 * np.arctanh(math.sin(lat * math.pi/180)) +
|
|
148
|
+
|
|
149
|
+
np.arctanh(math.sin(L * math.pi/180))) )
|
|
150
|
+
|
|
151
|
+
|
|
152
|
+
|
|
153
|
+
return int(tileX/256), int(tileY/256)
|
|
154
|
+
|
|
155
|
+
|
|
156
|
+
|
|
157
|
+
# 今回の指定エリアの中心座標(緯度、経度)\n",
|
|
158
|
+
|
|
159
|
+
center_lat = ( min_lat + max_lat ) / 2
|
|
160
|
+
|
|
161
|
+
center_lon = ( min_lon + max_lon ) / 2
|
|
162
|
+
|
|
163
|
+
|
|
164
|
+
|
|
165
|
+
# 今回はズームレベル zを 12とします。\n",
|
|
166
|
+
|
|
167
|
+
z = 12
|
|
168
|
+
|
|
169
|
+
# タイル座標を取得\n",
|
|
170
|
+
|
|
171
|
+
x, y = LatLon2Tile(center_lat,center_lon, z)
|
|
172
|
+
|
|
173
|
+
print("タイル座標(%d, %d, %d)" % (x, y, z))
|
|
174
|
+
|
|
175
|
+
|
|
176
|
+
|
|
177
|
+
# AVNIR-2のURL\n",
|
|
178
|
+
|
|
179
|
+
url_AVNIR2 = 'https://gisapi.tellusxdp.com/blend/av2ori/'
|
|
180
|
+
|
|
181
|
+
# Landsat-8のURL
|
|
182
|
+
|
|
183
|
+
url_Landsat8 = 'https://gisapi.tellusxdp.com/blend/landsat8/'
|
|
184
|
+
|
|
185
|
+
# ヘッダ\n",
|
|
186
|
+
|
|
187
|
+
headers={"Authorization": "Bearer " + API_TOKEN }
|
|
188
|
+
|
|
189
|
+
|
|
190
|
+
|
|
191
|
+
# AVNIR-2のバンド割り当て（赤：バンド3、緑：バンド2、青：バンド1）\n",
|
|
192
|
+
|
|
193
|
+
query_AVNIR2 = __”r=3&g=2&b=1”__
|
|
194
|
+
|
|
195
|
+
|
|
196
|
+
|
|
39
197
|
|
|
40
198
|
|
|
41
199
|
```ここに言語名を入力
|
|
@@ -44,7 +202,7 @@
|
|
|
44
202
|
|
|
45
203
|
|
|
46
204
|
|
|
47
|
-
### 試したこと
|
|
205
|
+
### 試したこと__イタリックテキスト__
|
|
48
206
|
|
|
49
207
|
|
|
50
208
|
|