GIS之arcgis系列08:arcpy实现批量excel转矢量点

文件夹内有很多excel表,每个excel表内有很多数据,爬取每条数据的日期、精度、纬度,根据经纬度将数据转换成矢量点,并指定坐标系。


完整版代码如下:

代码均已经过测试,可直接copy至arcgis工具箱中使用。

(copy后注意改成适用于自己数据的代码,微调相关信息)

python 复制代码
# -*- coding: utf-8 -*-
import pandas as pd
import arcpy
import os

input_folder = arcpy.GetParameterAsText(0)
output_folder = arcpy.GetParameterAsText(1)
if not os.path.exists(output_folder):
    os.makedirs(output_folder)
spatial_ref = arcpy.SpatialReference(4326)
temp_folder = output_folder

def process_excel_file(file_path, output_folder, temp_folder):
    df = pd.read_excel(file_path, header=0)

    columns_to_delete = [3, 4, 5, 6, 7, 10, 12, 13, 15, 16, 18, 21, 22]
    df.drop(df.columns[columns_to_delete], axis=1, inplace=True)

    new_column_names = {
        df.columns[2]: 'temmean',
        df.columns[3]: 'pre',
        df.columns[4]: 'temmax',
        df.columns[5]: 'windmax',
        df.columns[6]: 'temmin',
        df.columns[7]: 'rainmax'
    }
    df.rename(columns=new_column_names, inplace=True)

    date_col = 'Date'
    df[date_col] = pd.to_datetime(df[date_col])

    unique_dates = df[date_col].dt.date.unique()

    for date in unique_dates:
        date_df = df[df[date_col].dt.date == date]
        date_str = date.strftime('%Y%m%d')
        doy = date.timetuple().tm_yday
        temp_csv = os.path.join(temp_folder, "{}.csv".format(date_str))
        date_df.to_csv(temp_csv, index=False, encoding='utf-8')

        temp_layer = "temp_layer"
        if arcpy.Exists(temp_layer):
            arcpy.Delete_management(temp_layer)

        arcpy.MakeXYEventLayer_management(
            temp_csv, 'LON', 'LAT', temp_layer, spatial_ref
        )

        year_folder = file_path.split("\\")[-1].split(".")[0]
        print (year_folder)
        year_folder_path = os.path.join(output_folder, year_folder)
        if not os.path.exists(year_folder_path):
            os.makedirs(year_folder_path)

        output_shapefile = os.path.join(year_folder_path, "{}_{}.shp".format(year_folder, doy))
        arcpy.FeatureClassToFeatureClass_conversion(temp_layer, year_folder_path, os.path.basename(output_shapefile))


        os.remove(temp_csv)

for filename in os.listdir(input_folder):
    if filename.endswith(".xlsx"):
        file_path = os.path.join(input_folder, filename)
        process_excel_file(file_path, output_folder, temp_folder)

---------------完------------

相关推荐
南草徽9 天前
arcgis分割 (Split)
arcgis
张太行_9 天前
MySQL与Excel比较
数据库·mysql·excel
cwtlw9 天前
Excel学习03
笔记·学习·其他·excel
郭优秀的笔记9 天前
计算本地Excel某两列的差异值
excel
涔溪10 天前
实现 el-table 中键盘方向键导航功能vue2+vue3(类似 Excel)
前端·vue.js·excel
云博客-资源宝10 天前
Excel函数大全
机器学习·excel·概率论
BillKu10 天前
【前后前】导入Excel文件闭环模型:Vue3前端上传Excel文件,【Java后端接收、解析、返回数据】,Vue3前端接收展示数据
java·前端·excel
野生工程师11 天前
【ArcGIS】在线影像底图调用
arcgis
焚 城11 天前
Excel数据导出小记
.net·excel