专栏名称: GEE遥感训练营
专注GEE遥感算法,包括遥感影像下载、遥感影像制图、遥感GIS空间分析、遥感生态评价、遥感影像融合、遥感去干扰等多元遥感云计算
目录
相关文章推荐
51好读  ›  专栏  ›  GEE遥感训练营

利用GEE仅仅只要1分钟实现SCI精美制图!

GEE遥感训练营  · 公众号  ·  · 2025-02-10 14:52

正文

00

精华部分


我们利用GEE,实现了1分钟制作任意地区的研究区位置图,当然也可以是其他地图,功能非常炸裂 ,结果如下:

01

概况


俗话说: 好的图表决定了你的SCI是否能被录用。 以前,我们使用GEE进行精美制图的时候,通常需要借助arcgis来进行编辑处理,最后实现SCI的精美制图。 GEE集成了海量的数据,为制图提供了丰富的元素。 但是,它也存在无法进行编辑,只能导出数据以后进行制图。因此, 之前我们进行制图的时候,都是利用GEE生成数据,Arcgis上进行快速制图。
我们之前的做法可以参照这篇推文: GEE 多源数据综合制图 图片

02

GEE快速制图分析


GEE无法编辑窗口进行制图的弊端,能否通过其他方法解决呢? 答案是肯定的。 python集成了海量的第三方应用,可以实现地图编辑和快速制图。 由于 吴秋生老师开发了geemap等库,可以快速实现对GEE接口的调用。 同时,python可以快速实现制图, 这就为GEE+python快速制图提供了便利。 因此,通过python调用geemap,结合python第三方库,就能轻松解决。
首先,要实现GEE快速制图,需要安装以下几个库:
(1)geemap;(2)matplotlib
安装方法很简单,直接使用pip或者conda进行install就行。

03

GEE快速制图流程


(1) 第三方库导入 ,具体如下:
import eeimport geemapfrom geemap import cartoeeimport matplotlib.pyplot as plt

值得注意的是,geemap的cartoee库是吴老师之前的书上提到的内容,不清楚的可以在以下网址学习到:https://book.geemap.org/chapters/08_cartoee.html

(2) 端口设置与初始化 :我们需要设置好端口,并对GEE进行初始化,通过调用token实现对GEE进行连接。这里,我们的端口为10809,每个人的端口有所不同,可以通过计算机属性查询得到,代码如下:
geemap.set_proxy(port=10809)ee.Initialize()
(3) 数据输入与筛选 :本次数据我们以四川省为例,具体使用省-市-县哪一级别,根据自己的需求进行修改,矢量数据我们之前在公众号提供了2024年免费的 审图号:GS(2024)0650号 )数据,具体可以查看推文: GEE矢量详解教程|利用最新官方审图号的中国行政区获取研究区范围!审图号:GS(2024)0650号 栅格数据我们以高层DEM为例,添加了水系(2019年的Nature数据)美观背景 ,代码如下:
china_provinces = ee.FeatureCollection("projects/ee-tilmacatanla/assets/boundry/china_provinces");china_city = ee.FeatureCollection("projects/ee-tilmacatanla/assets/boundry/china_city");china_county = ee.FeatureCollection("projects/ee-tilmacatanla/assets/boundry/china_county");roi = china_provinces.filter(ee.Filter.eq('name','四川省'))# print(roi.getInfo())styling_sichuan = {"color": "FF0000", "width": 2, "lineType": "solid", "fillColor": "00000000"}styling_provinces = {"color": "000000", "width": 1, "fillColor": "00000000"}river = (ee.FeatureCollection('WWF/HydroSHEDS/v1/FreeFlowingRivers')           .filter(ee.Filter.lt('RIV_ORD',5)))styling_river = {"color": "#37c4d6", "width": 0.6, "fillColor": "00000000"}Map=geemap.Map()Map.addLayer(river.style(**styling_river), {}, "river")Map.addLayer(china_provinces.style(**styling_provinces), {}, "Province")Map.addLayer(roi.style(**styling_sichuan), {}, "sichuan")Map.centerObject(roi,6)Map







请到「今天看啥」查看全文


推荐文章
互联网分析沙龙  ·  24 个高效率产品设计工具
8 年前
格上财富  ·  会聊天的人,一开口就赢了
7 年前