气象数据处理方法:spss和Excel
一、 下载原始txt数据中的经纬度处理:将度分处理成度,Excel处
理
首先除以100,处理成小数格式,这里第一个实际是52度58分,
在Excel中用公式:=LEFT(O2,FIND(\
需注意:
当为整数时,值为空,这时需查找出来手动修改,或者将经纬度这一列的小数位改成两位再试试,可能好使(这个我没尝试)
第二步:
将经纬度转换成投影坐标,在arcgis实现
将Excel中的点导入arcgis,给定坐标系为wgs84地理坐标,然后投影转换成自己定义的等面积的albers投影(因为anusplina软件需要投影坐标,这里转换成自己需要的坐标系)
第三步:spss处理
将下载的txt数据导入spss之后,编辑变量属性,删掉不需要的列,然后将最后需要的那些变量进行数据重组
本实验下载的数据是日均温数据,全国800+个站点 2012年366天的数据。相当于有800+ * 366行数据
1.变量
变量属性:变量属性这里的设置决定了在SPLINA这个模块中输入数据的格式,本实验spss处理的气象数据的格式统一用这个:(A5,2F18.6,F8.2,F8.2),一共5列。
即:台站号,字符串,5位;
经纬度:都是浮点型,18位,6个小数位 海拔:浮点型,8位,2个小数位 日均温:浮点型,8位,2个小数位
2.数据重组,将个案重组成变量:
后几步都默认就行:
重组之后结果:变成了800+行,370列,就相当于数据变成了:行代表每个站点,列是代表每一天的数据。
3. 因为anusplin这个软件需要的是投影坐标,在重组完的基础上,将经纬度这两列替换成投影之后的经纬度。
方法1:直接复制粘贴即可
方法二:用合并文件,添加变量功能
运行结果:
下面说明:
本实验中站点59981的位置在南沙那边,不在现有的DEM范围内,固可以删掉这个点不用。 如果在arcgis中插值,删掉这个点
如果在anusplin中,可以通过设置XY范围自动舍去这个点,固可以不删 此处统一将这个点删掉. 所以一共838个站点
变量属性:变量属性这里的设置决定了在SPLINA这个模块中输入数据的格式,本实验spss处理的气象数据的格式统一用这个:(A5,2F18.6,F8.2,F8.2),一共5列。
即:台站号,字符串,5位;
经纬度:都是浮点型,18位,6个小数位 海拔:浮点型,8位,2个小数位 日均温:浮点型,8位,2个小数位
4. 最后:保存成固定ASCII格式
保存成这种格式之后就可以被anusplin软件读取了。
总结:anusplin软件设置的一个易错点:设置xy范围中,splina的范围要比第二步lapgrd中dem的范围大~和dem一样或者比DEM小都会出错
一、软件xy设置,以后可用这个范围(全国的):
第一步- splina 3054296.000 2595703.750
945413.000 6945413.750 或者:-3054296.52594703.5
860413.104751626945413.10475162 海拔:-400到9000
第二步lapgrd:-2954296.500 2494703.500
960413.104751626845413.10475162
二、 本实验spss处理的气象数据的格式统一用这个:
(A5,2F18.6,F8.2,F8.2) 4.37版本的软件中SPLINA的操作
1.界面介绍
2.运行输出
运行完之后会出现很多文件:
三、 4.2版本的软件中Cmd界面中输入命令:
1. 程序和数据介绍
2. 运行过程
然后出现如下DOS界面:
然后在这个界面下输入语句
3. 本实验中输入如下:
1. h:
2. cd h:\\anusplin\\run(定位到本机存放cmd,exe和数据的位置,cmd,exe和数据要放
在同一个文件夹下)
3. lapgrd
4. 解释lapgrd
Dos界面中08011.log,是输出的日志文件的名字,这里边的“名字”自己设置,是用来指定输出日志文件的名字,无影响。