博文
ArcGIS学习笔记03：全国地级行政区年降水量、年均气温求算
|||
1 准备
1.1 数据源
已插值好的全国年降水量、年均气温数据，ArcGIS的ASCII格式txt文件，每年一个文件。降水数据有45年，气温数据有48年。数据为0.5经纬度分辨率。
1993年全国县级行政区划图，ArcGIS面要素shp文件。
1.2 目的
得到逐年的全国地级行政区年降水量、年均气温序列值。
1.33 软件
ArcGIS 10.0简体中文版
2 步骤
4.1 数据预处理
4.1.1 行政区预处理
从http://www.gov.cn/test/2011-08/22/content_1930111.htm获得2010年全国县级以上行政单位行政编码，用以把1993年全国县级行政区划图整理成2010年地级行政区划图，确保所有的行政编码和对应的行政区名正确。
可以先使用县级行政区划图，把其中的行政编码替换为2010年的行政编码，核对正确。而后根据地级行政区的行政编码规律，使用ArcGIS的合并功能将县级行政区合并为地级行政区。
从水资源公报中获取被水资源公报所承认的地级行政区，与上面对比，整理成符合水资源公报的2010年地级行政区划图，共得到359个地级行政区及县级行政区（有的县级行政区是地级编制，有的是省直管，在水资源公报中认为与地级行政区平级）。
核对shp文件的属性，确保最重要的行政编码无误。
4.1.2 气温、降水文件预处理
打开ArcGIS或ArcCatalog的ArcToolBox，找到转换工具→转为栅格→ASCII转栅格，右键→批处理。
在“输入ASCII栅格文件”中选择txt格式的温度或降水文件，在“输出栅格”中输入对应的输出栅格文件名。可以新建一个地理空间数据库，把所有的文件存在其中。
使用Excel辅助生成批处理表。
点“确定”。
这样就能把txt文件转换成栅格了。但是这些栅格文件都没有投影，需要再定义投影。
打开ArcGIS或ArcCatalog的ArcToolBox，找到数据管理工具→投影和变换→定义投影，右键→批处理。
把刚才生成的栅格文件的绝对路径填写进“输入数据集或要素类”，在“坐标系”中选择需要的投影坐标，最好与行政区图的shp文件的坐标一致，可以导入其坐标。
使用Excel辅助生成批处理表。
点“确定”。
根据数据源，生成的栅格数据的分辨率为0.5经纬度，对一些较小的地级、县级行政区来说还是比较大，可能造成数据空缺，需要把增加栅格分辨率。这有两种方法。一种是数据导出，一种是栅格重采样。数据导出即将要导出的数据由ArcGIS打开，然后右键→导出，格式为grid，分辨率选需要的，例如0.1经纬度，输入要生成的栅格文件名即可。这种方法没有批处理功能，适用于数据量较小的情况下。鉴于数据量较大，使用栅格重采样的方法。
打开ArcGIS或ArcCatalog的ArcToolBox，找到数据管理工具→栅格→栅格数据处理→重采样→右键，重采样。
在“输入栅格”中选择以上生成的0.5经纬度分辨率的温度、降水栅格数据，在“输出栅格数据集”中写入要生成的0.1经纬度分辨率的温度、降水栅格数据文件名，在“输出像元大小”中写上0.1，“重采样技术”使用默认的nearest。
使用Excel辅助生成批处理表。
点“确定”。
4.2 区域分析
打开ArcGIS或ArcCatalog的ArcToolBox，找到Spatial Analyst工具。如果无法运行，找菜单栏中的“自定义”菜单，点“扩展模块”，将其选中即可。
在Spatial Analyst工具→区域分析→以表格显示分区统计，右键，批处理。
在“输入栅格数据或要素区域数据”中选择做好的地级行政区shp文件；“区域字段”选择唯一值，建议使用地级或县级行政区的行政编码；“输入赋值栅格”选择预处理好的0.1经纬度分辨率的栅格数据；“输出表”写上输出的表所在的文件地址；“在计算中忽略NoData”使用默认值true；“统计类型”可以选需要的，也可以全部要，即“all”。
使用Excel辅助生成批处理表。
点“确定”。
4.3 转换为数据表
打开ArcGIS或ArcCatalog的ArcToolBox，找到转换工具→转为dBASE→表转dBASE（批量）。
把4.2中生成的所有表都加入“输入表”框中，选择输出文件夹，点“确定”。
4.4 整理数据
把4.3生成的dbf文件逐个用excel打开，将数据使用excel或access整理即可。
https://blog.sciencenet.cn/blog-485-646387.html
上一篇：ArcGIS学习笔记（02）：按字段属性批量生成shp文件
下一篇：windows10窗口切换问题的可能解决办法
