GEE、USGS、地理空间数据云上下载武汉地区的影像数据
目錄
一、GEE
二、USGS
三、地理空間數(shù)據(jù)云
總結(jié)
這次的主要任務(wù)是學(xué)會(huì)下載哨兵二號的遙感影像數(shù)據(jù)，學(xué)長推薦了幾個(gè)常用的網(wǎng)站如GEE、EOS、USGS、地理空間數(shù)據(jù)云，我計(jì)劃把幾個(gè)網(wǎng)站都嘗試一下，下面就按順序進(jìn)行嘗試：
一、GEE
網(wǎng)址：https://code.earthengine.google.com
GEE(Google Earth Engine)，即谷歌地球引擎,據(jù)說是世界上使用量最大的云平臺(tái)（如果真是這樣的話，那說明遙感和地信的分量確實(shí)很重了）。
登錄這個(gè)網(wǎng)站需要科學(xué)上網(wǎng)，首先需要用一個(gè)谷歌賬號登錄(關(guān)于谷歌賬號的注冊門道也有點(diǎn)多，這里就不細(xì)說了），然后用該賬號提出申請，申請界面如下：
雖然這里的郵箱使用的是qq郵箱，但還是很不建議使用，最好使用學(xué)校的官方郵箱，申請成果的可能性會(huì)大一些。
機(jī)構(gòu)就選擇所在的大學(xué)，類型選擇acadamic即可，國家選擇China，至于最后的“申請賬號來干什么”，我當(dāng)時(shí)是先用百度翻譯把我的申請理由給翻譯成英文后再復(fù)制上去的，具體理由就是說一下要做研究之類的。我當(dāng)時(shí)這樣申請之后可以說是秒通過了（郵箱里會(huì)收到通知），這也是在之前失敗了n次之后的結(jié)果。
來到主界面，是這個(gè)樣子的：
首先導(dǎo)入武漢的矢量文件，注意是壓縮包（zip）格式：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
通過如下代碼可實(shí)現(xiàn)導(dǎo)入：
? ? ? ??
導(dǎo)入結(jié)果：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
在資源搜索欄中尋找哨兵二號的數(shù)據(jù)集，并選擇導(dǎo)入（IMPORT）：
? ? ? ? ?
可以看到代碼編輯區(qū)為我們準(zhǔn)備好了代碼：
接下來使用如下代碼下載影像：
var district = table;var dsize = district.size(); var district_geometry = district.geometry(); print(dsize); Map.centerObject(district_geometry,7);//去云 function maskS2clouds(image) {var qa = image.select('QA60');// Bits 10和11分別是云和卷云var cloudBitMask = 1 << 10;var cirrusBitMask = 1 << 11;// Both flags should be set to zero, indicating clear conditions.var mask = qa.bitwiseAnd(cloudBitMask).eq(0).and(qa.bitwiseAnd(cirrusBitMask).eq(0));return image.updateMask(mask).divide(10000); }//處理數(shù)據(jù)，生成數(shù)據(jù)集 var dataset = ee.ImageCollection('COPERNICUS/S2').filterBounds(district_geometry)//區(qū)域篩選.filterDate('2020-01-01', '2020-12-30')//時(shí)間范圍// 選擇云量.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 1)).map(maskS2clouds).select(['B8', 'B4', 'B3']);//可視化參數(shù) var rgbVis = {min: 0.0,max: 0.3,bands: ['B8', 'B4', 'B3'],//按照8波段、4波段、3波段放入通道合成 };var IMAGE = dataset.mosaic();//下載數(shù)據(jù) Export.image.toDrive({image:IMAGE,description:'image_really_none',//設(shè)置下載任務(wù)tasks的名稱scale:10,//分辨率是10mmaxPixels: 1e13,//單幅影像輸出的最大像元數(shù)region:district_geometry,//下載范圍fileFormat: 'GeoTIFF',//設(shè)置影像導(dǎo)出格式，注意GeoTIFF格式需要所有波段存儲(chǔ)類型一致（如不能同時(shí)存Int16和Int32）formatOptions: {cloudOptimized: true}});可以在右上角的Tasks欄中看到加載的任務(wù)：
? ? ? ? ? ? ? ? ??
點(diǎn)擊“Open in Drive”,可以從云端下載影像：
將下載后的影像加載到ENVI中進(jìn)行查看：
二、USGS
網(wǎng)址：EarthExplorer (usgs.gov)
USGS是美國地質(zhì)調(diào)查局的官方網(wǎng)站，同樣需要科學(xué)上網(wǎng)，我這里使用的是學(xué)校郵箱注冊的官網(wǎng)賬號，登陸后就可以進(jìn)行數(shù)據(jù)下載了。
首先，我還是通過想通過矢量數(shù)據(jù)直接定位武漢，所以就選擇"KML/Shapefile Upload",再選擇Shapefile格式，加載文件（注意這里加載的文件必須是zip壓縮文件形式）：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
但是很遺憾的是出現(xiàn)了如下狀況：
這里顯示我的武漢市范圍矢量文件不只一個(gè)多邊形，我有些奇怪，就把文件加載到arcgis里看了下，果然在右上角發(fā)現(xiàn)了一個(gè)小范圍的多邊形：
所以我選擇從行政區(qū)劃的角度來選擇，在左側(cè)欄中輸入關(guān)鍵字“wuhan”即可得到結(jié)果：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
為了找出下載的具體范圍，我們需要進(jìn)行選定，這里方法比較多，比如說下圖中的畫圓法，先點(diǎn)擊一下可得圓心，再點(diǎn)擊一下可得半徑：? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?
也可以選擇多邊形法，這樣畫的范圍可能更為準(zhǔn)確:
下面還可以選擇影像時(shí)間，默認(rèn)為全部時(shí)間段的影像：
接下來選擇數(shù)據(jù)源，這里是哨兵二號：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
接下來可以選擇一些具體的參數(shù)，如云量：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
接下來就可以檢索到符合條件的一些影像數(shù)據(jù)了，可以選擇第一個(gè)像腳印一樣的符號來把影像加載到地圖上來查看范圍（自動(dòng)會(huì)用不同顏色標(biāo)識）：
選擇第五個(gè)圖標(biāo)即可進(jìn)行下載：
? ? ? ??
三、地理空間數(shù)據(jù)云
網(wǎng)址：地理空間數(shù)據(jù)云 (gscloud.cn)
在這個(gè)上面下數(shù)據(jù)是比較簡單的了，之前做遙感實(shí)習(xí)時(shí)的數(shù)據(jù)也主要是從這上面下的。打開官網(wǎng)后選擇“高級檢索”，可以看到界面是非常簡潔明了的：
選擇數(shù)據(jù)集為哨兵二號，選擇用矢量文件來劃定區(qū)域：
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ??
比較奇怪的是，這里用矢量數(shù)據(jù)劃定區(qū)域后找不到合適的影像（可能和之前在USGS下載時(shí)的問題一樣，由于矢量文件不只一個(gè)多邊形所導(dǎo)致）：
所以這里就采取另一種方法：先根據(jù)行政區(qū)劃找到覆蓋武漢市影像，然后再用武漢市范圍的矢量數(shù)據(jù)進(jìn)行剪裁，這就需要找到能夠覆蓋武漢市全部范圍的圖像，總共是3張（之前做實(shí)習(xí)下的是Landsat8的數(shù)據(jù)，也是三張即可覆蓋武漢全境），注意要檢索能夠下載的數(shù)據(jù)，因?yàn)楹芏鄶?shù)據(jù)能被檢索到但我們卻不能下載：
總結(jié)
哨兵二號的數(shù)據(jù)可以從很多網(wǎng)站上下載得到，除了上面介紹的這些還有很多，擔(dān)上面這三種應(yīng)當(dāng)可以滿足大部分時(shí)候的需求了。
總結(jié)
以上是生活随笔為你收集整理的GEE、USGS、地理空间数据云上下载武汉地区的影像数据的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 01大数据之路
- 下一篇： 批量修改图片 *jfif 格式