1.目标:近期需要对气象站点数据进行空间插值,前提是将含有气候数据的站点数据在空间中展示出来。
2. 问题:GIS中可以通过添加表数据———为表数据添加xy坐标——数据导出为shp文件
但是涉及年份较多的时候批量转换可以节省很多时间
3.方法(依旧找到了两种方法。R中应该也能实现,只是我一时摸索不出来)
3.1 Python接口
# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
import glob
import os
import arcgisscripting
gp=arcgisscripting.create()
spatial_ref=r'F:\Shujuchuli\Export_Output_feature_to_shp.prj' #转换点文件需要的投影信息数据prj
env.workspace =r'F:\Shujuchuli\Pre_huanglongxian' #工作环境设置、输入路径
pathout=r'F:\Shujuchuli\Pre_csv_to_shp' #输出文件路径
# 点文件信息
x_corrods='Longitude' #经度字段
y_corrods='Latitude' #纬度字段
Z_filed='MEAN_year' # 转换点文件需要保留的字段
try:
for file1 in arcpy.ListFiles("*.csv"): #列出当前工作环境下csv文件
print file1
info=os.path.basename(file1).split('.')[0] #获取csv文件的基本文件名,并去掉后缀
intable=file1 #转换