（Arcgis）python geopandas库分割shp属性表特定内容，批量导出shp文件

本文主要是介绍（Arcgis）python geopandas库分割shp属性表特定内容，批量导出shp文件，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

一、两个文件。实现目标：从1个shp文件，根据属性表内容提取成200个shp文件，文件名取自txt文本内容

shp文件（要素1-200.shp）：打开属性表前14项相同，后200项不相同。
————任务目标：需要输出200个shp（属性表前14项相同，循环后200项单独提取出1项）

txt文本（numbered.txt）内容：
1 要素1
2 要素2
... ...（总共200行）

import geopandas as gpd# 文件路径
file_path = r"D:\要素\要素1-200.shp"
species_names_file = r"D:\numbered.txt"# 读取要素名称
species_dict = {}
with open(species_names_file, 'r', encoding='utf-8') as file:for line in file:if line.strip():number, name = line.strip().split(maxsplit=1)species_dict[int(number)] = f"{number} {name}"  # 存储编号和名称，包括空格# 读取shapefile
data = gpd.read_file(file_path)# 确定几何列的名称
geometry_column = data.geometry.name# 获取前14个基础要素的列名
base_columns = data.columns[:14].tolist()
if geometry_column not in base_columns:base_columns.append(geometry_column)  # 确保几何列被包括# 循环处理每个要素
for i in range(1, 200):  # 物种编号从1到200species_column = f'要素{i}'  # 构建要素列名if species_column in data.columns:# 选择基础要素和当前要素的列，确保包括几何列selected_columns = base_columns + [species_column]selected_data = data[selected_columns]# 使用包括编号和名称的字符串作为文件名if i in species_dict:filename = species_dict[i]  # 从字典中获取包含编号和名称的字符串new_file_path = fr"D:\要素\{filename}.shp"else:filename = f"{i} 要素_selected"  # 如果编号不存在于字典中，使用默认名称new_file_path = fr"D:\要素\{filename}.shp"# 导出到新的shapefile，使用UTF-8编码selected_data.to_file(new_file_path, encoding='utf-8')print(f"{filename}的数据已成功导出到 {new_file_path}")

二、最后每个要素需要放在单独的文件夹中，文件夹名命名为文件名

文件名                 文件夹名
1 要素1.shp            1 要素1
2 要素2.shp            2 要素2
......

import os
import shutil# 包含 shapefile 文件的目录
source_dir = r"D:\要素\要素1-200"# 遍历源目录中的所有文件
for filename in os.listdir(source_dir):if filename.endswith(".shp"):  # 检查文件是否为 shapefile# 去掉文件扩展名得到基本名称base_name = os.path.splitext(filename)[0]# 使用基本名称创建新的目录路径new_dir = os.path.join(source_dir, base_name)# 如果目录不存在则创建目录if not os.path.exists(new_dir):os.makedirs(new_dir)# 移动 shapefile 及其关联文件for ext in ['.shp', '.shx', '.dbf', '.prj', '.cpg', '.sbn', '.sbx', '.qix', '.fbn', '.fbx', '.ain', '.aih', '.ixs', '.mxs', '.atx']:old_file = os.path.join(source_dir, f"{base_name}{ext}")if os.path.exists(old_file):shutil.move(old_file, new_dir)print(f"已将 {filename} 及其关联文件移动到 {new_dir}")

这篇关于（Arcgis）python geopandas库分割shp属性表特定内容，批量导出shp文件的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

---