我們知道，想發(fā)表SCI期刊論文，除了需要非常好的idea以及非常嚴(yán)謹(jǐn)且優(yōu)美的寫(xiě)作之外，還需要在論文的圖表等方面下功夫。很多時(shí)候，一張圖做好了，審稿人看到了會(huì)心情愉悅，論文發(fā)表后作者本身也會(huì)覺(jué)得是一種享受。

一直以來(lái)，高質(zhì)量的SCI期刊圖表都是在本地做的。為了制作精美的圖片，需要把數(shù)據(jù)都在本地準(zhǔn)備好。例如，我們想實(shí)現(xiàn)下面的圖片，就需要下載DEM數(shù)據(jù)到本地然后再出圖，這種方式還是比較繁瑣的。

既然Google Earth Engine已經(jīng)把各類(lèi)遙感影像數(shù)據(jù)以及一些高級(jí)算法都集成到GEE云平臺(tái)了，那么有沒(méi)有可能把作圖等環(huán)節(jié)也搬到GEE云平臺(tái)上呢？一直以來(lái)，我們都以為GEE出圖比較丑，不適合做這些高質(zhì)量的精美圖片。但是，最近發(fā)現(xiàn)我們錯(cuò)了，使用GEE依然能夠做出非常好看的SCI期刊圖片。

比如，參考上圖右半部分的研究區(qū)圖片，我們可以在GEE上直接利用DEM數(shù)據(jù)做出如下的圖片。兩者圖片雖然在局部地區(qū)存在一些差異，但基本上實(shí)現(xiàn)了最初的功能。最主要的是，整個(gè)作圖過(guò)程僅僅花費(fèi)了10.5秒便完成了。

接下來(lái)，我們就詳細(xì)介紹如何在GEE上面實(shí)現(xiàn)這一精美的圖片。需要指出的是，這個(gè)圖片不是用JavaScript網(wǎng)頁(yè)版本的GEE做出的，而是使用的Python版本。我們使用的IDE是VS Code，配置起來(lái)也很方便。我們最新的寒假課程就講授了如何使用VS Code訪問(wèn)GEE以及有關(guān)配置問(wèn)題（見(jiàn)寒假加油|人工智能與云平臺(tái)大數(shù)據(jù)圖像處理基礎(chǔ)和提升研修班）。

第一步，導(dǎo)入需要的模塊

這里主要是cartopy、geemap和cartoee等模塊，其中cartoee模塊已經(jīng)內(nèi)嵌到geemap，無(wú)需再安裝。本地安裝geemap的時(shí)候，只需使用 pip install geemap -i https://pypi.douban.com/simple/即可完成安裝。

第二步，授權(quán)GEE的訪問(wèn)

主要是讓本地的IDE能夠訪問(wèn)GEE數(shù)據(jù)，使用的是下面語(yǔ)句：

# geemap.ee_initialize()
geemap.set_proxy(port=33210)
ee.Authenticate()
ee.Initialize(project='ee-scistudy')

注意，這里的port后面的端口是每個(gè)人魔法的端口，project后面的項(xiàng)目號(hào)則是自己的project名稱(chēng)。這里我們就不贅敘了，有不懂的請(qǐng)聯(lián)系小編加群咨詢。

第三步，調(diào)用dem數(shù)據(jù)和矢量數(shù)據(jù)，并在地圖中先加載看一下效果，代碼如下：

經(jīng)過(guò)這一步可以得到下面的地圖。在線加載出這個(gè)地圖表明距離最后的成功就不遠(yuǎn)了。

最后一步，將上面的地圖使用plt制作成我們需要的圖片，代碼如下：

bbox = [135.8,39,120.2,55]

import matplotlib.pyplot as plt
import cartopy.io.img_tiles as cimgt

plt.rcParams['font.family'] = 'Times New Roman' 
#  統(tǒng)一設(shè)置
from proplot import rc

# 統(tǒng)一設(shè)置字體

rc['tick.labelsize'] = 20
# 統(tǒng)一設(shè)置xy軸名稱(chēng)的字體大小
rc['axes.labelsize'] = 20
# # 統(tǒng)一設(shè)置軸刻度標(biāo)簽的字體粗細(xì)
rc['axes.labelweight'] = 'bold'
# 統(tǒng)一設(shè)置xy軸名稱(chēng)的字體粗細(xì)
rc['tick.labelweight'] = 'bold'


fig = plt.figure(figsize=(20, 22),facecolor='white')

ax = cartoee.get_map(imgBlend, region=bbox)

cb = cartoee.add_colorbar(ax,vis_params=vis_params,loc='right',orientation='vertical')

# ax.set_title(label='Study Area', fontsize=20)


# add gridlines to the map at a specified interval
cartoee.add_gridlines(ax, interval=[2, 2], linestyle='--')

plt.show()

經(jīng)過(guò)這么幾步，我們就可以在GEE中做出上文中出現(xiàn)的精美圖片。而且大家在做這種圖的時(shí)候，還可以自己加地圖以及換成別的不同的類(lèi)型，比如地形圖DTM可以換成下面的NDVI圖，或者是土地利用的圖片。整個(gè)過(guò)程十分的快捷方便，真的是太舒服了！?。?/span>