GLC3020002010202030m
# global land 30介绍
GlobeLand30是中国研制的30米空间分辨率全球地表覆盖数据,目前有GlobeLand30 2000、2010版、2020版。
# 1 数学基础
GlobeLand30数据采用WGS-84坐标系。
南纬85°-北纬85°之间的区域,投影方式采用UTM投影, 6度分带,坐标单位为米,坐标不加带号。
南北纬85°-90°之间的区域,投影方式采用极地方位投影,投影面切于地球南北极点。
# 2 分类系统
GlobeLand30数据共包括10个一级类型,分别是:耕地、林地、草地、灌木地、湿地、水体、苔原、人造地表、裸地、冰川和永久积雪。分类系统定义见表1。
类型 | 内容 | 代码 |
---|---|---|
耕地 | 用于种植农作物的土地,包括水田、灌溉旱地、雨养旱地、菜地、牧草种植地、大棚用地、以种植农作物为主间有果树及其他经济乔木的土地,以及茶园、咖啡园等灌木类经济作物种植地。 | 10 |
林地 | 乔木覆盖且树冠盖度超过30%的土地,包括落叶阔叶林、常绿阔叶林、落叶针叶林、常绿针叶林、混交林,以及树冠盖度为10-30%的疏林地。 | 20 |
草地 | 天然草本植被覆盖,且盖度大于10%的土地,包括草原、草甸、稀树草原、荒漠草原,以及城市人工草地等。 | 30 |
灌木地 | 灌木覆盖且灌丛覆盖度高于30%的土地,包括山地灌丛、落叶和常绿灌丛,以及荒漠地区覆盖度高于10%的荒漠灌丛。 | 40 |
湿地 | 位于陆地和水域的交界带,有浅层积水或土壤过湿的土地,多生长有沼生或湿生植物。包括内陆沼泽、湖泊沼泽、河流洪泛湿地、森林/灌木湿地、泥炭沼泽、红树林、盐沼等。 | 50 |
水体 | 陆地范围液态水覆盖的区域,包括江河、湖泊、水库、坑塘等。 | 60 |
苔原 | 寒带及高山环境下由地衣、苔藓、多年生耐寒草本和灌木植被覆盖的土地,包括灌丛苔原、禾本苔原、湿苔原、高寒苔原、裸地苔原等。 | 70 |
人造地表 | 由人工建造活动形成的地表,包括城镇等各类居民地、工矿、交通设施等,不包括建设用地内部连片绿地和水体。 | 80 |
裸地 | 植被覆盖度低于10%的自然覆盖土地,包括荒漠、沙地、砾石地、裸岩、盐碱地等。 | 90 |
冰川和永久积雪 | 由永久积雪、冰川和冰盖覆盖的土地,包括高山地区永久积雪、冰川,以及极地冰盖等。 | 100 |
苔原的四个亚类型灌木苔原、禾本苔原、湿苔原与裸地苔原对应编码为71、72、73、74。
# 3 数据内容
GlobeLand30分幅数据由地表覆盖数据文件、坐标信息文件、分类影像接图表文件、元数据文件4部分组成。其中:
\1) 地表覆盖数据文件是指存储分幅地表覆盖分类信息的文件;
\2) 坐标信息文件是指记录分幅数据坐标信息的文件;
\3) 分类影像接图表文件是指记录分类所用的主要影像范围及获取时间的矢量文件;
\4) 元数据文件指记录分幅数据元数据信息的文件。
为便于数据拷贝和下载传输,将分幅数据进行压缩处理,形成了一个zip格式的压缩文件,压缩文件中包括上述4类文件。
# 4 存储单元
根据所位于的纬度情况,GlobeLand30数据采用3种分幅方式存储,具体如下:
1)在南北纬60°区域内,按照5°(纬度)×6°(经度)的大小进行分幅,其中央经线与其所在6度带一致;
2)在南北纬60°至85°区域内,按照5°(纬度)×12°(经度)的大小进行分幅。由于分幅跨2个6度带的图幅,按照奇数6度带的中央经线进行投影。
3)在南北纬85°以上区域,保存为1幅。
数据名称由16位字符组成,具体规则如下:
南北纬缩写(1位)+6度带号(2位)+“”+起始纬度(2位)+“”+产品年代(4位)+LC(地表覆盖缩写)+分辨率(3位)。其中:
1)南北纬纬度缩写:N表示北纬,S表示南纬。
2)6度代号:图幅所在的6度分带的带号值。对于跨2个6度带的图幅,代号应填写奇数带号。
3)起始纬度:在北半球填写图幅下边缘线的纬度值,在南半球填写图幅上边缘线的纬度值,见图1。
图1 南北半球分幅起始纬度示意图
4)产品年代:表示地表覆盖产品年代版本,如2000年、2010年等。
5)分辨率:LC产品的分辨率,只表示整米,不满3位的在前面加0。如30米用030表示。
对于南北纬85°以上区域,数据命名按:南北纬缩写(1位)+00(2位)+“”+85(2位)+“”+产品年代(4位)+LC(地表覆盖缩写)+分辨率(3位)。
例如:N49_25_2000LC030图幅,表示的是北半球49带(经度范围东经108°-114°)北纬25°至30°矩形区域内的2000年期30米地表覆盖数据。
不同类型文件命名规则如下:
1)地表覆盖数据文件命名:分幅数据名称+“.tif”
2)坐标信息文件命名:分幅数据名称+“.tfw”
3)分类影像接图表文件命名:分幅数据名称+“_IMG.shp”
4)元数据文件命名:分幅数据名称+ “_MAT.xml”(2000和2010版) /“_MAT.xls”(2020版)
图2 GlobeLand30 V2000和V2010数据分幅分布
GlobeLand30 V2000版和V2010版分幅总数均为853幅,其分布如图2所示。GlobeLand30 V2020版分幅总数为966幅,其分布如图3所示。
图3 GlobeLand30 V2020数据分幅分布
# 5 数据源
GlobeLand30数据研制所使用的分类影像主要是30米多光谱影像,包括美国陆地资源卫星(Landsat)的TM5、ETM+、OLI多光谱影像和中国环境减灾卫星(HJ-1)多光谱影像,2020版数据还使用了16米分辨率高分一号(GF-1)多光谱影像。
影像选取原则是:在保障影像无云(少云)前提下,择优选择数据生产基准年或更新年度年±2年内植被生长季的多光谱影像。对于获取困难区域,可放宽影像获取时间,确保影像全球覆盖的完整性。
# 6 精度
GlobeLand30 V2010数据精度评价由同济大学牵头完成。从全球853幅数据中抽取80个图幅,布设超过15万个检验样本,得出GlobeLand30 V2010数据的总体精度为83.50%,Kappa系数0.78。
GlobeLand30 V2020数据精度评价由中国科学院空天信息创新研究院牵头完成。基于景观形状指数抽样模型进行全套数据布点,共布设样本超过23万个。得出GlobeLand30 V2020数据的总体精度为85.72%,Kappa系数0.82。
# 7.数据下载
GlobeLand30是我国首套全球公开的基础地理信息数据产品,具有里程碑意义,数据可以直接从官网下载:
官方网站:http://www.globeland30.org/home.html?type=data
# 8.数据镶嵌
(数据镶嵌参考:公众号 走天涯徐小洋地理数据科学)
先把下载好的文件解压,TIFF文件都放到一个文件夹里
使用Model builder构建处理流,投影,SetNull
需要注意的是,模型构建器里面SetNull需要用SQL语句来写,图形化界面是灰的,不能直接用,SQL语句粘贴下面的代码就行
Value = 0 Or Value = 255
由于数据量很大,需要漫长的等待……
Mosaic to New Raster,镶嵌至新栅格
Extract by mask,按掩膜提取,处理结果如下所示:
GLC30—2020年全国镶嵌结果
GLC30—2010年全国镶嵌结果
GLC30—2000年全国镶嵌结果