【ArcGIS风暴】ArcGIS中国地表覆盖数据GlobeLand30预处理(批量投影、拼接、掩膜提取)附成品下载
结果示意图:
GlobeLand30是30米空间分辨率全球地表覆盖数据,目前可供下载的有3年的数据:2000-2010-2020。本文主要讲解在ArcGIS10.6平台下进行GlobeLand30的预处理操作,主要预处理步骤包括:
批量分幅投影转换
、批量分幅无效值处理
、批量图幅拼接
和数据集掩膜提取
。
文章目录
- 一、GlobeLand30 ArcGIS批处理完整流程
- 1. 批量投影转换
- 2. 批量处理无效值
- 3. 批量图幅拼接
- 4. 添加分类字段并赋值
- 5. 掩膜提取
- 二、GlobeLand30成品数据下载
一、GlobeLand30 ArcGIS批处理完整流程
GlobeLand30的预处理步骤主要有:批量分幅投影转换
、批量分幅无效值处理
、批量图幅拼接
和数据集掩膜提取
。
1. 批量投影转换
由于原数据集使用的是6°带UTM投影,图幅从43-53跨了11个带,因此需要先对所有的图幅进行投影转换,然后进行后续操作。为了便于精确计算面积,本文将其转换为Albers等积投影,地理坐标系为WGS1984。
Python源代码:
由于原始数据集分幅是保存在不同的文件夹,因此第一步批量投影变换需要遍历根目录下的所有文件夹,过程稍微有些复杂,从第二步开始,将所有的处理结果保存在一个文件夹内,减少了代码量。
arcpy.ProjectRaster_management(raster,out, "PROJCS['MyAlbers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['central_meridian',105.0],PARAMETER['Standard_Parallel_1',25.0],PARAMETER['Standard_Parallel_2',47.0],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]", "NEAREST", cellsize , "", "", "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")n=n+1print str(n)+" rasters are processed !!!"
查看投影坐标系:
2. 批量处理无效值
第一步投影结束后,影像的黑边依然存在,这将影像后面的拼接处理结果,黑边对应的像元值为0,只需要将其设置为NoData即可,对应的工具是栅格计算器中的SetNull函数,我们继续使用Python批处理。
Python源码:
for inRaster in inRasters:outRaster=outpath+inRaster##print outRasterprint "Processing "+inRaster +" ......"outSetNull = SetNull(inRaster,inRaster,WhereClause)outSetNull.save(outRaster)n=n+1
print str(n)+" rasters are processed!"
处理过程:
处理结果:
3. 批量图幅拼接
在ArcGIS中拼接栅格图幅有两个工具:【镶嵌】和【镶嵌至新栅格】,当然了需要手动添加数据,好在我们的数据在同一个文件夹下,可以直接全选拖动,当然了,也可以写Python代码来完成。
方法一:使用【镶嵌至新栅格】工具
方法二:Python代码
arcpy.MosaicToNewRaster_management(str, outPath, outName, "", "8_BIT_UNSIGNED", "", "1", "LAST", "FIRST")
处理过程:
需要注意的是,原始数据集是tif格式,且每个图幅都是有金字塔的,加载和显示速度都比较快,但是SetNull后会丢失金字塔,加载速度很慢。
4. 添加分类字段并赋值
打开属性表,添加一个文本型字段type,打开编辑器,输入分类名称并保存。
5. 掩膜提取
使用【按掩膜提取】工具,按研究区范围进行提取。
配色参考值:
注意:拼接后的数据集是有金字塔的,但是使用【按掩膜提取】工具后又会丢失金字塔,必须重新构建金字塔,否则加载速度异常慢。
二、GlobeLand30成品数据下载
请私信作者。
【ArcGIS风暴】ArcGIS中国地表覆盖数据GlobeLand30预处理(批量投影、拼接、掩膜提取)附成品下载相关推荐
- 【GIS风暴】30米分辨率地表覆盖数据GlobeLand30原始数据集简介及下载地址
数据集预览: GlobeLand30是30米空间分辨率全球地表覆盖数据,目前可供下载使用的有3年的数据:2000-2010-2020,本文主要讲述GlobeLand30的官网下载地址和数据集简介. 数 ...
- 使用QGIS分区统计工具实现栅格分类数据的分区计算面积——GlobeLand30地表覆盖数据为例
在栅格分析中,常常碰到使用分类后的栅格数据按照特定分区统计面积的需求,今天,我将使用QGIS的分区统计工具,演示地表覆盖数据按照地表分类分区域统计面积的过程,希望能给有这方面需求的朋友提供参考. 0 ...
- GlobeLand30地表覆盖数据下载及预处理(ArcGIS去黑边,镶嵌,裁剪,符号化)
1.下载土地利用数据 GlobeLand30数据共有三个年份(2000,2010,2020),选择自己需要的图幅后,提交下载申请 2.去黑边 由于数据有黑边,需要进行处理后再进行拼接 打开ArcGI ...
- 2020年ESA中国区10m地表覆盖数据下载
开始写博客这么久,但没多少文章,说来惭愧. 从本篇开始,内容会逐步涉及数据获取.处理等实际的过程.也会跟gma库 (详见自建 | 地理与气象数据分析(geographic and meteor ...
- 大数据治理平台建设方案(文末附PDF下载)
这份材料我给满分!分享一份非常好的大数据治理平台解决方案材料,这份PPT将理论与实践相结合,值得仔细阅读,建议收藏. 文档目录主要包含了以下几点: 数据治理概述 某行数据现状及问题 数据治理阶段目标 ...
- GlobeLand30地表覆盖数据下载及预处理(去黑边,镶嵌,裁剪)
一.数据下载 网址:Welcome-GlobeLand30 提供免费下载,包括2000年.2010年.2020年 登录--选择数据--提交申请--审核通过后即可下载 二.栅格数据预处理(去黑边.拼接. ...
- ArcGIS和Global Mapper中按分幅图框进行批量DOM分幅裁剪命名(附练习数据)
工作中按分幅矢量对影像数据进行分幅并按要求命名是我们常做的一个操作,今天小助手就汇总下在ArcGIS和Global Mapper中进行批量分幅. 所需数据: Shp矢量 正射影像 文末有数据下载链接 ...
- Arcgis之国土报备(报盘数据)txt批量转shp格式工具箱
在勘测定界.土地报件.增减挂钩.开发整理等相关项目中经常需要反将标准txt格式转换到SHP矢量格式 1.TXT格式依据标准<勘测定界界址点坐标交换格式> 2.软件功能:工具实现了将土地报 ...
- Arcgis实例操作20--线分割面数据、删除线多余节点、提取点群最小边界几何
1.如何用线要素对面要进行素分割?(更方便的方法,建议使用) 1)加载分割面工具. 打开[菜单栏]|[自定义]|[自定义模式]对话框,选择[命令]|[拓扑]|[分割面工具],将其拖动至任意工具条内后释 ...
最新文章
- 【Gradle】管理库工程release及debug
- 在主机中通过xdebug远程调试Vagrant虚拟机中drush脚本的方法
- 【转】学习asp.net比较完整的流程
- dataguru北京线下沙龙-第二部 《Oracle 索引优化思路--案例分享 -- 刘盛》
- photoshop 插件_什么是Photoshop插件,扩展程序和附加组件?
- 大学生学科竞赛管理系统/竞赛管理系统的设计与实现
- 计算机控制系统信号恢复之零阶保持器
- python生成泊松分布_Python Numpy泊松分布
- caxa自动保存的文件在哪里_CAXA CAD教程之软件临时文件的清理
- 傻瓜式安装卸载office
- 以赛促产 以赛引才 |第六届世界智能大会·中国华录杯数据湖算法大赛正式启动
- Windows UI 测试
- 谢烟客---------Linux之find查找
- PnL Explained FAQ
- vue怎么给pc端浏览器设置一个最小屏幕_图文讲解uni-app的PC宽屏适配方案
- 3-33在图 3-31中,以太网交换机有6个接口,分别接到5台主机和一个路由器。在下面表中的“动作”一栏中,表示先后发送了4个帧。假定在开始时,以太网交换 机的交换表是空的。试把该表中其他的栏目都填写
- 水逆了一整年的王源,2020年年初能靠《大主宰》翻身吗?
- 【学术前沿分析】1 论文数据统计
- word文档的数据损坏如何进行修复
- 第一个OpenDayLight项目:HelloWorld
热门文章
- SQL Server的数据导入MySQL数据库方法简介
- C# WPF中添加调试信息查看窗体
- 《Redis核心技术与实战》学习总结(1)
- WPF开源项目:AIStudio.Wpf.AClient
- Masuit.Tools,一个免费的轮子
- 跟我一起学.NetCore之MVC过滤器,这篇看完走路可以仰着头走
- 使用 C# sdk 连接 高可用的 rabbitmq 镜像集群
- docker registry 镜像同步
- 使用ASP.NET Core 3.x 构建 RESTful API - 4.3 HTTP 方法的安全性和幂等性
- 居然不知道和的区别?