如何写一个shp转txt的脚本

以下是一个简单的shp转txt的脚本，使用Python编写：

import shapefile

# 打开shp文件
sf = shapefile.Reader("path/to/shapefile.shp")

# 打开txt文件
with open("path/to/output.txt", "w") as f:
    # 遍历每个记录
    for record in sf.records():
        # 将记录中的每个字段写入txt文件
        for field in record:
            f.write(str(field) + "\t")
        f.write("\n")

注：需要安装shapefile库，可以通过pip install pyshp安装。

相关问题

如何在ArcGIS中编写一个脚本来进行txt转shp

在ArcGIS中，您可以使用Python语言编写一个脚本来批量将TXT文件转换为SHP文件。下面是一个简单的示例脚本，可以帮助您了解如何实现这个功能：

import arcpy
import os

# 设置工作空间
arcpy.env.workspace = r"C:\data"

# 遍历指定目录下的所有TXT文件
for root, dirs, files in os.walk(arcpy.env.workspace):
    for file in files:
        if file.endswith(".txt"):
            # 获取文件路径和名称
            txtFile = os.path.join(root, file)
            fileName = os.path.splitext(file)[0]

            # 创建点要素集
            featureClass = arcpy.CreateFeatureclass_management(
                arcpy.env.workspace, fileName, "POINT")

            # 添加字段
            arcpy.AddField_management(featureClass, "FileName", "TEXT")

            # 打开游标并逐行读取TXT文件中的坐标数据
            cursor = arcpy.da.InsertCursor(featureClass, ["SHAPE@", "FileName"])
            with open(txtFile, "r") as f:
                for line in f:
                    # 解析坐标数据，创建点要素，并写入属性字段
                    x, y = line.strip().split(",")
                    point = arcpy.Point(float(x), float(y))
                    cursor.insertRow((point, fileName))

            # 释放游标
            del cursor

这个脚本的主要思路是遍历指定目录下的所有TXT文件，并逐个读取文件中的坐标数据。对于每个TXT文件，脚本会创建一个点要素集，并添加一个属性字段用于存储文件名。然后，脚本会打开游标并一行一行地读取TXT文件中的坐标数据，解析数据并创建点要素，并将文件名写入属性字段中。最后，脚本会保存点要素集为SHP文件。

需要注意的是，脚本中的代码可能需要根据您的具体情况进行修改。比如，您需要将工作空间设置为自己的数据存储路径，确保TXT文件中的坐标数据格式和分隔符与代码中的匹配，以及根据需要调整点要素的属性字段等。

如何在ArcGIS中编写一个shp转txt的脚本并且可以自定义输出坐标系、生产单位、时间、shp字段信息

在ArcGIS中编写一个shp转txt的脚本可以使用Python语言来实现，具体步骤如下：

1. 打开ArcGIS，点击“工具”菜单，选择“Python”窗口；
2. 在Python窗口中输入以下代码：

import arcpy

# 输入shp文件路径和txt文件路径
in_shp = r"C:\path\to\input.shp"
out_txt = r"C:\path\to\output.txt"

# 自定义输出坐标系
out_coor_system = arcpy.SpatialReference("WGS 1984")

# 自定义生产单位、时间和shp字段信息
production_unit = "xxx"
time = "2021-01-01"
field_info = "field1,field2,field3"

# 创建游标，读取shp中的要素
cursor = arcpy.da.SearchCursor(in_shp, ["SHAPE@XY"] + field_info.split(","))

# 打开输出文件，写入标题行
with open(out_txt, "w") as f:
    f.write("X,Y,%s,ProductionUnit,Time\n" % field_info)

    # 遍历每个要素，写入坐标和属性信息
    for row in cursor:
        x, y = row[0]
        attrs = ",".join(str(a) for a in row[1:])
        f.write("%.6f,%.6f,%s,%s,%s\n" % (x, y, attrs, production_unit, time))

# 完成后打印提示信息
print("转换完成！")

3. 将代码中的`in_shp`和`out_txt`分别替换为要转换的shp文件路径和输出的txt文件路径；
4. 根据需要自定义输出坐标系、生产单位、时间和shp字段信息，分别替换代码中对应的变量；
5. 运行代码，在Python窗口中或者调用脚本文件运行即可。

注意事项：
- 如果要输出其他坐标系的文件，需要在`arcpy.SpatialReference()`中输入对应的坐标系信息；
- 如果要输出其他类型的文件，比如csv或者xlsx，需要修改代码中的文件打开方式和写入方式；
- 如果shp文件中的字段名中有空格或者特殊字符，需要在代码中将其用双引号括起来，比如`"field name"`。