【ArcGIS Pro二次开发】(65)：进出平衡SHP转TXT、TXT转SHP

最近一个小伙伴提了这么一个需求，需要把TXT和SHP进行互转。

这种TXT文件其实遇到了好几个版本，都有一点小差异。之前已经做过一个TXT转SHP的工具，但好像不适用。于是针对这个版本，做了互转的2个工具。

---

【SHP转TXT】

一、要实现的功能

 

如上图所示，在【数据处理】组—【TXT相关】面板下，点击【进出平衡@SHP转TXT】工具。

在弹出的工具框中，分别输入参数：

1、输入SHP文件所在的文件夹。

2、输入TXT文件所在的文件夹。

3、这里不用填写，会自动列出【1】中所有的shp要素，不想转换的可以点复选框取消。

4、选择地块名称和地块用途，对应的字段值会写入TXT文件中。

生成结果如下：

---

二、实现流程

核心代码直接贴上，注释已经写得比较清楚了。

需要注意的是，这里获取要素的点集信息，我采取了通过要素的JSON文本来截取的做法，原因是不知道怎么用API来获取，所以用了这么个取巧的方法，以后要是学会了再改吧。这也是无奈之举，要学的东西还很多。

foreach (string fullPath in list_shpPath)
{
    // 初始化写入txt的内容
    string txt_all = "[地块坐标]" + "\r";
 
    string shp_name = fullPath[(fullPath.LastIndexOf(@"\") + 1)..];  // 获取要素名
    string shp_path = fullPath[..(fullPath.LastIndexOf(@"\"))];  // 获取shp名
 
    // 打开shp
    FileSystemConnectionPath fileConnection = new FileSystemConnectionPath(new Uri(shp_path), FileSystemDatastoreType.Shapefile);
    using FileSystemDatastore shapefile = new FileSystemDatastore(fileConnection);
    // 获取FeatureClass
    FeatureClass featureClass = shapefile.OpenDataset(shp_name);
 
    using (RowCursor rowCursor = featureClass.Search())
    {
        int featureIndex = 1;
        while (rowCursor.MoveNext())
        {
            using (Feature feature = rowCursor.Current as Feature)
            {
                // 获取地块名称，地块性质
                Row row = feature as Row;
                string ft_name = "";
                string ft_type = "";
                var areaName = row[field_mc];
                var areaType = row[field_yt];
                if (areaName != null) { ft_name = areaName.ToString(); }
                if (areaType != null) { ft_type = areaType.ToString(); }
 
                // 获取面要素的JSON文字
                Geometry polygon = feature.GetShape();
                string js = polygon.ToJson().ToString();
 
                // 解析JSON文字
                // 取坐标点文字
                string cod = js[(js.IndexOf("[[[") + 3)..js.IndexOf("]]]")];
                // 坐标点列表
                List<string> list_xy = cod.Split("]]").ToList();
                for (int i = 0; i < list_xy.Count; i++)
                {
                    // 坐标行
                    List<string> xy_detils = list_xy[i].Replace(",[[", "").Split("],[").ToList();
 
                    // 加一行title
                    int count = xy_detils.Count;    // 点的个数
                    string title = $"{count},{ft_name},面,{ft_type},@ " + "\r";
                    txt_all += title;
 
                    for (int j = 0; j < xy_detils.Count; j++)
                    {
                        // 点序号
                        int index = j + 1;
                        if (index == xy_detils.Count) { index = 1; }
                        // XY坐标点
                        string x = Math.Round(double.Parse(xy_detils[j].Split(",")[0]), 4).ToString();
                        string y = Math.Round(double.Parse(xy_detils[j].Split(",")[1]), 4).ToString();
                        // 加入文本
                        txt_all += $"J{index},{featureIndex},{x},{y}\r";
                    }
                }
 
            }
            featureIndex++;
        }
    }
    // 写入txt文件
    string txtPath = @$"{folder_txt}\{shp_name.Replace(".shp", "")}.txt";
    if (File.Exists(txtPath))
    {
        File.Delete(txtPath);
    }
    File.WriteAllText(txtPath, txt_all);
}

---

【TXT转SHP】

一、要实现的功能

 

如上图所示，在【数据处理】组—【TXT相关】面板下，点击【进出平衡@TXT转SHP】工具。

在弹出的工具框中，分别输入参数：

1、输入TXT文件所在的文件夹。

2、输入SHP文件所在的文件夹。

3、选择正确的坐标系。

4、这里不用填写，会自动列出【1】中所有的TXT文件，不想转换的可以点复选框取消。

生成结果如下：

---

二、实现流程

TXT转SHP之前已经已经写过一篇文章，可以参看：

【ArcGIS Pro二次开发】(41)：勘测定界txt文件转数据库(批量)https://blog.csdn.net/xcc34452366/article/details/131309938不过文章中用的TXT文件和今天这个有些不同，代码也就不一样，但思路是一致的。

这里只放核心代码，代码中用到一些自定义方法和之前是一样的引用，就不再一一放上：

foreach (string txtPath in list_txtPath)
{
    string shp_name = txtPath[(txtPath.LastIndexOf(@"\") + 1)..].Replace(".txt", "");  // 获取要素名
 
    pw.AddProcessMessage(@$"创建要素：{shp_name}", 10, time_base, Brushes.Black);
 
    // 创建一个空要素
    Arcpy.CreateFeatureclass(shpPath, shp_name, "POLYGON", spatial_reference);
    // 新建字段
    Arcpy.AddField(@$"{shpPath}\{shp_name}.shp", "地块名称", "TEXT");
    Arcpy.AddField(@$"{shpPath}\{shp_name}.shp", "地块性质", "TEXT");
 
    // 打开shp
    FileSystemConnectionPath fileConnection = new FileSystemConnectionPath(new Uri(shpPath), FileSystemDatastoreType.Shapefile);
    using FileSystemDatastore shapefile = new FileSystemDatastore(fileConnection);
    // 获取FeatureClass
    FeatureClass featureClass = shapefile.OpenDataset(shp_name);
 
    // 预设文本内容
    string text = "";
    // 获取txt文件的编码方式
    Encoding encoding = ToolManager.GetEncodingType(txtPath);
    // 读取【ANSI和UTF-8】的不同+++++++（ANSI为0，UTF-8为3）
    // 我也不知道具体原理，只是找出差异点作个判断，以后再来解决这个问题------
    int encoding_index = int.Parse(encoding.Preamble.ToString().Substring(encoding.Preamble.ToString().Length - 2, 1));
 
    if (encoding_index == 0)        // ANSI编码的情况
    {
        Encoding.RegisterProvider(CodePagesEncodingProvider.Instance);
        using (StreamReader sr = new StreamReader(txtPath, Encoding.GetEncoding("GBK"))) { text = sr.ReadToEnd(); }
    }
    else if (encoding_index == 3)               // UTF8编码的情况
    {
        using (StreamReader sr = new StreamReader(txtPath, Encoding.UTF8)) { text = sr.ReadToEnd(); }
    }
 
    // 文本中的【@】符号放前
    string updata_text = ChangeSymbol(text);
 
    // 获取要素txt列表的地块标记
    List<string> Parts = GetParts(updata_text);
 
    for (int i = 0; i < Parts.Count; i++)
    {
        // 地块编号、地块性质
        string dkmc = "";
        string dkxz = "";
 
        // 根据换行符分解坐标点文本
        List<string> lines = Parts[i].Split("@").ToList();
        // 创建空坐标点集合
        var vertices_list = new List>();
 
        for (int j = 1; j < lines.Count; j++)
        {
            var vertices = new List();
            vertices_list.Add(vertices);
        }
        // 构建坐标点集合
        for (int k = 1; k < lines.Count; k++)
        {
            List<string> list_point = lines[k].Split("\r").ToList();
            foreach (var point in list_point)
            {
                if (!point.Contains(","))     // 跳过无坐标部份的文本
                {
                    continue;
                }
                else if (!point.StartsWith("J"))     // 名称、地块编号、功能文本
                {
                    dkmc = point.Split(",")[1];
                    dkxz = point.Split(",")[3];
                }
                else           // 点坐标文本
                {
                    double lat = double.Parse(point.Split(",")[2]);         // 经度
                    double lng = double.Parse(point.Split(",")[3]);         // 纬度
                    vertices_list[k - 1].Add(new Coordinate2D(lat, lng));    // 加入坐标点集合
                }
            }
 
        }
 
        /// 构建面要素
        // 创建编辑操作对象
        EditOperation editOperation = new EditOperation();
        editOperation.Callback(context =>
        {
            // 获取要素定义
            FeatureClassDefinition featureClassDefinition = featureClass.GetDefinition();
            // 创建RowBuffer
            using RowBuffer rowBuffer = featureClass.CreateRowBuffer();
            // 写入字段值
            rowBuffer["地块名称"] = dkmc;
            rowBuffer["地块性质"] = dkxz;
 
            PolygonBuilderEx pb = new PolygonBuilderEx(vertices_list[0]);
            // 如果有空洞，则添加内部Polygon
            if (vertices_list.Count > 1)
            {
                for (int i = 0; i < vertices_list.Count - 1; i++)
                {
                    pb.AddPart(vertices_list[i + 1]);
                }
            }
            // 给新添加的行设置形状
            rowBuffer[featureClassDefinition.GetShapeField()] = pb.ToGeometry();
 
            // 在表中创建新行
            using Feature feature = featureClass.CreateRow(rowBuffer);
            context.Invalidate(feature);      // 标记行为无效状态
        }, featureClass);
 
        // 执行编辑操作
        editOperation.Execute();
    }
    // 保存编辑
    Project.Current.SaveEditsAsync();
}

---

三、工具文件分享

我把工具都集合成工具箱，不再单独放单个工具，可以到这里下载完整工具箱，会不断更新：

【ArcGIS Pro二次开发】：CC工具箱https://blog.csdn.net/xcc34452366/article/details/131506345PS：可以直接点击...bin\Debug\net6.0-windows\下的.esriAddinX文件直接安装。