设为首页收藏本站CAD论坛CAD教程CAD下载明经币充值
切换到宽版

明经CAD社区

 找回密码
 注册

QQ登录

只需一步，快速开始

搜索
明经CAD社区»社区 › 开发篇 › ObjectARX/VC++ 编程技术 › 利用gdal库读写shp简单代码示例
返回列表 发新帖
查看: 826|回复: 3

利用gdal库读写shp简单代码示例

  [复制链接]
gzxl
发表于 2025-8-16 10:10:14 | 显示全部楼层 |阅读模式
本帖最后由 gzxl 于 2025-8-16 10:12 编辑

利用 gdal 库读写 shp 简单代码示例

包含目录：添加 gdal 对应所在的 include
库目录：添加 gdal 对应所在的 lib
附加依赖项： gdal.lib

一 简介
  Shape files是ESRI提供的一种矢量数据格式，它没有拓扑信息，一个Shape files由一组文件组成，其中必要的基本文件包括坐标文件（.shp）、索引文件（.shx）和属性文件（.dbf）三个文件。

  1 坐标文件（.shp）的结构说明
    坐标文件(.shp)用于记录空间坐标信息。它由头文件和实体信息两部分构成

  2. 属性文件(.dbf)的结构说明
    属性文件(.dbf)用于记录属性信息。它是一个标准的DBF文件，也是由头文件和实体信息两部分构成

  3. 索引文件（.shx）的结构说明
    索引文件（.shx）主要包含坐标文件的索引信息，文件中每个记录包含对应的坐标文件记录距离坐标文件的文件头的偏移量。

二 《建设工程规划条件核实测量成果数据标准》
  1、分层及文件格式
  
  2、属性表
    属性表分类包含：
      2.1 建筑类规划核实属性表(图形文件)
      2.2 建筑类规划核实属性表(表格文件)
      2.3 市政类规划核实工程属性表（图形文件）
    截图其中一个，其他详细参加《建设工程规划条件核实测量成果数据标准》
      2.3 市政类规划核实工程属性表（图形文件）
        2.3.1 规划核实（交通工程）信息 （HY_JTGCGHHS）
        

三 代码示例：依据前面的数据标准输出ARC/INFO SHP格式
  1. TCHAR szKword[10];
  2. acedInitGet(0, _T("A B"));
  3. int nReturn = acedGetKword(_T("\n输入分层及文件格式选项[交通工程(A)/管线工程(B)]:"), szKword);
  4. if (nReturn == RTNORM)
  5. {
  6.     struct resbuf *rbuf;  // 结果缓冲区链表
  7.     ads_name ssname;
  8.     int res;
  9.     if (_tcscmp(szKword, _T("A")) == 0) // 交通工程
  10.     {
  11.         rbuf = acutBuildList(RTDXF0, _T("LWPOLYLINE"), RTNONE);
  12.         res = acedSSGet(TEXT("X"), NULL, NULL, rbuf, ssname);
  13.         acutRelRb(rbuf);
  14.         if (res != RTNORM)
  15.         {
  16.             acutPrintf(_T("\n图面没有可用的交通道路线(多段线)!\n"));
  17.             return;
  18.         }
  19.     }
  20.     else if (_tcscmp(szKword, _T("B")) == 0) // 管线工程
  21.     {
  22.         rbuf = acutBuildList(RTDXF0, _T("LINE"), RTNONE);
  23.         res = acedSSGet(TEXT("X"), NULL, NULL, rbuf, ssname);
  24.         acutRelRb(rbuf);
  25.         if (res != RTNORM)
  26.         {
  27.             acutPrintf(_T("\n图面没有可用的管线(直线)!\n"));
  28.             return;
  29.         }
  30.     }
  31.     Adesk::Int32 length;
  32.     acedSSLength(ssname, &length);
  33.     ads_name ename;
  34.     AcDbObjectId endId;

  36.     vector<vector<OGRPoint>> vertexPoint;
  37.     int noCloseNum = 0;
  38.     double dMaxArea = 0.0, dMaxLength = 0.0;
  39.     for (int i = 0; i < length; i++)
  40.     {
  41.         acedSSName(ssname, i, ename);
  42.         acdbGetObjectId(endId, ename);
  43.         AcDbEntity* pEnt = NULL;
  44.         acdbOpenObject(pEnt, endId, AcDb::kForRead);
  45.         if (_tcscmp(szKword, _T("A")) == 0) // 交通工程
  46.         {
  47.             if (pEnt->isKindOf(AcDbPolyline::desc()))
  48.             {
  49.                 AcDbPolyline *pPolyline = AcDbPolyline::cast(pEnt);
  50.                 if (!pPolyline->isClosed())
  51.                 {
  52.                     pPolyline->close();
  53.                     noCloseNum++;
  54.                     continue;
  55.                 }
  56.                 double dArea = 0.0;
  57.                 pPolyline->getArea(dArea);
  58.                 if (dArea > dMaxArea)
  59.                     dMaxArea = dArea;
  60.                 vector<OGRPoint> ogrPts;
  61.                 for (int j = 0; j < pPolyline->numVerts(); j++)
  62.                 {
  63.                     AcGePoint3d verPoint;
  64.                     if (pPolyline->getPointAt(j, verPoint) == Acad::eOk)
  65.                     {
  66.                         OGRPoint point;
  67.                         point.setX(verPoint.x);
  68.                         point.setY(verPoint.y);
  69.                         ogrPts.push_back(point);
  70.                     }
  71.                 }
  72.                 pPolyline->close();
  73.                 vertexPoint.push_back(ogrPts);
  74.                 ogrPts.clear();
  75.             }
  76.         }
  77.         else if (_tcscmp(szKword, _T("B")) == 0) // 管线工程
  78.         {
  79.             if (pEnt->isKindOf(AcDbLine::desc()))
  80.             {
  81.                 AcDbLine *pLine = AcDbLine::cast(pEnt);
  82.                 AcGePoint3d startPt, endPt;
  83.                 pLine->getStartPoint(startPt);
  84.                 pLine->getEndPoint(endPt);
  85.                 vector<OGRPoint> ogrPts;
  86.                 OGRPoint ogrStartPt, ogrEndPt;
  87.                 ogrStartPt.setX(startPt.x);
  88.                 ogrStartPt.setY(startPt.y);
  89.                 ogrEndPt.setX(endPt.x);
  90.                 ogrEndPt.setY(endPt.y);
  91.                 ogrPts.push_back(ogrStartPt);
  92.                 ogrPts.push_back(ogrEndPt);
  93.                 vertexPoint.push_back(ogrPts);
  94.                 ogrPts.clear();
  95.                 double dist = (AcGePoint2d(startPt.x, startPt.y)).distanceTo(AcGePoint2d(endPt.x, endPt.y));
  96.                 dMaxLength += dist;
  97.                 pLine->close();
  98.             }
  99.         }
  100.         pEnt->close();
  101.     }
  102.     acedSSFree(ssname);

  104.     int nFormatWidth = 0;
  105.     if (dMaxArea != 0.0)
  106.     {
  107.         CString sArea = ToString(dMaxArea, 2);
  108.         for (int i = 0; i < sArea.GetLength(); i++)
  109.         {
  110.             if ((sArea >= '0' && sArea <= '9') || sArea == '.')
  111.                 nFormatWidth++;
  112.         }
  113.     }
  114.     if (dMaxLength != 0.0)
  115.     {
  116.         CString sMaxLength = ToString(dMaxLength, 2);
  117.         for (int i = 0; i < sMaxLength.GetLength(); i++)
  118.         {
  119.             if ((sMaxLength >= '0' && sMaxLength <= '9') || sMaxLength == '.')
  120.                 nFormatWidth++;
  121.         }
  122.     }
  123.     if (nFormatWidth == 0)
  124.         nFormatWidth = 8;

  126.     if (vertexPoint.size() > 0)
  127.     {
  128.         CString ReturnPach;
  129.         Getbrowse(ReturnPach); // 获得文件夹路径
  130.         if (_tcscmp(ReturnPach, _T("")) == 1)
  131.         {
  132.             CString sPrefixName;
  133.             if (_tcscmp(szKword, _T("A")) == 0) // 交通工程
  134.                 sPrefixName = _T("HY_JTGCGHHS");
  135.             else if (_tcscmp(szKword, _T("B")) == 0) // 管线工程
  136.                 sPrefixName = _T("HY_GXGCGHHS");

  138.             CString sFileName = ReturnPach + _T("\\") + sPrefixName + _T(".shp");

  140.             const char *pszDriverName = "ESRI Shapefile";
  141.             // 注册所有的文件格式驱动
  142.             GDALDriver *poDriver;
  143.             GDALAllRegister();
  144.             // 得到shp文件的处理器
  145.             poDriver = GetGDALDriverManager()->GetDriverByName(pszDriverName);
  146.             if (poDriver == NULL)
  147.             {
  148.                 acutPrintf(_T("\n%s 获取驱动程序 ESRI Shapefile 错误!\n"), pszDriverName);
  149.                 return;
  150.             }
  151.             // 创建shp文件
  152.             GDALDataset *poDS;
  153.             const char* fileName = CStringToCharArray(sFileName);
  154.             poDS = poDriver->Create(fileName, 0, 0, 0, GDT_Unknown, NULL); //创建shp文件
  155.             if (poDS == NULL)
  156.             {
  157.                 acutPrintf(_T("\n创建输出文件失败!\n"));
  158.                 return;
  159.             }
  160.             // 创建图层
  161.             OGRLayer *poLayer;
  162.             if (_tcscmp(szKword, _T("A")) == 0) // 交通工程
  163.             {
  164.                 poLayer = poDS->CreateLayer("HY_JTGCGHHS", NULL, wkbPolygon, NULL);
  165.                 OGRFieldDefn KJFWXXIDField("KJFWXXID", OFTString); // 空间范围信息ID
  166.                 OGRFieldDefn YSDMField("YSDM", OFTString);         // 要素代码
  167.                 OGRFieldDefn DZJGHField("DZJGH", OFTString);       // 电子监管号
  168.                 OGRFieldDefn YDMJField("YDMJ", OFTReal);           // 用地面积
  169.                 OGRFieldDefn JTGCXMLXField("JTGCXMLX", OFTString); // 交通工程项目类型
  170.                 OGRFieldDefn DLDJField("DLDJ", OFTString);         // 道路等级
  171.                 OGRFieldDefn DLCDField("DLCD", OFTReal);           // 道路长度
  172.                 OGRFieldDefn BZHDMKDField("BZHDMKD", OFTReal);     // 标准横断面宽度
  173.                 OGRFieldDefn GDJTLXField("GDJTLX", OFTString);     // 轨道交通类型
  174.                 OGRFieldDefn GDCDField("GDCD", OFTReal);           // 轨道长度
  175.                 OGRFieldDefn KZXKDField("KZXKD", OFTReal);         // 轨道控制线宽度
  176.                 KJFWXXIDField.SetWidth(36);
  177.                 YSDMField.SetWidth(10);
  178.                 DZJGHField.SetWidth(19);
  179.                 YDMJField.Set("YDMJ", OFTReal, nFormatWidth, 2, OJRight);
  180.                 JTGCXMLXField.SetWidth(2);
  181.                 DLDJField.SetWidth(10);
  182.                 DLCDField.Set("DLCD", OFTReal, nFormatWidth, 2, OJRight);
  183.                 BZHDMKDField.Set("BZHDMKD", OFTReal, nFormatWidth, 2, OJRight);
  184.                 GDJTLXField.SetWidth(10);
  185.                 GDCDField.Set("GDCD", OFTReal, nFormatWidth, 2, OJRight);
  186.                 KZXKDField.Set("KZXKD", OFTReal, nFormatWidth, 2, OJRight);
  187.                 poLayer->CreateField(&KJFWXXIDField);
  188.                 poLayer->CreateField(&YSDMField);
  189.                 poLayer->CreateField(&DZJGHField);
  190.                 poLayer->CreateField(&YDMJField);
  191.                 poLayer->CreateField(&JTGCXMLXField);
  192.                 poLayer->CreateField(&DLDJField);
  193.                 poLayer->CreateField(&DLCDField);
  194.                 poLayer->CreateField(&BZHDMKDField);
  195.                 poLayer->CreateField(&GDJTLXField);
  196.                 poLayer->CreateField(&GDCDField);
  197.                 poLayer->CreateField(&KZXKDField);
  198.             }
  199.             else if (_tcscmp(szKword, _T("B")) == 0) // 管线工程
  200.             {
  201.                 poLayer = poDS->CreateLayer("HY_GXGCGHHS", NULL, wkbPolygon, NULL);
  202.                 OGRFieldDefn KJFWXXIDField("KJFWXXID", OFTString); // 空间范围信息ID
  203.                 OGRFieldDefn YSDMField("YSDM", OFTString);         // 要素代码
  204.                 OGRFieldDefn DZJGHField("DZJGH", OFTString);       // 电子监管号
  205.                 OGRFieldDefn SPGXLXField("SPGXLX", OFTString);     // 审批管线类型
  206.                 OGRFieldDefn GXFJField("GXFJ", OFTString);         // 审批管线分级
  207.                 OGRFieldDefn SPGXCDField("SPGXCD", OFTReal);       // 审批管线长度
  208.                 OGRFieldDefn JGGXCDField("JGGXCD", OFTReal);       // 竣工管线长度
  209.                 OGRFieldDefn SPGXGJField("SPGXGJ", OFTString);     // 审批管线管径
  210.                 OGRFieldDefn JGGXGJField("JGGXGJ", OFTString);     // 竣工管线管径
  211.                 KJFWXXIDField.SetWidth(36);
  212.                 YSDMField.SetWidth(10);
  213.                 DZJGHField.SetWidth(19);
  214.                 SPGXLXField.SetWidth(10);
  215.                 GXFJField.SetWidth(10);
  216.                 SPGXCDField.Set("SPGXCD", OFTReal, nFormatWidth, 2, OJRight);
  217.                 JGGXCDField.Set("JGGXCD", OFTReal, nFormatWidth, 2, OJRight);
  218.                 SPGXGJField.SetWidth(200);
  219.                 JGGXGJField.SetWidth(200);
  220.                 poLayer->CreateField(&KJFWXXIDField);
  221.                 poLayer->CreateField(&YSDMField);
  222.                 poLayer->CreateField(&DZJGHField);
  223.                 poLayer->CreateField(&SPGXLXField);
  224.                 poLayer->CreateField(&GXFJField);
  225.                 poLayer->CreateField(&SPGXCDField);
  226.                 poLayer->CreateField(&JGGXCDField);
  227.                 poLayer->CreateField(&SPGXGJField);
  228.                 poLayer->CreateField(&JGGXGJField);
  229.             }
  230.             if (poLayer == NULL)
  231.             {
  232.                 acutPrintf(_T("\n图层创建失败!\n"));
  233.                 GDALClose(poDS);
  234.                 poDS = NULL;
  235.                 vertexPoint.clear();
  236.                 return;
  237.             }
  238.             for (int i = 0; i < vertexPoint.size(); i++)
  239.             {
  240.                 vector<OGRPoint> ogrPts = vertexPoint.at(i);

  242.                 OGRFeature *poFeature;
  243.                 poFeature = OGRFeature::CreateFeature(poLayer->GetLayerDefn());

  245.                 // 创建几何和 Feature
  246.                 if (_tcscmp(szKword, _T("A")) == 0) // 交通工程
  247.                 {
  248.                     OGRLinearRing ogrring;
  249.                     int pNum = (int)ogrPts.size();
  250.                     ogrring.setNumPoints(pNum);
  251.                     for (int j = 0; j < ogrPts.size(); j++)
  252.                     {
  253.                         ogrring.setPoint(j, ogrPts[j].getX(), ogrPts[j].getY(), ogrPts[j].getZ());
  254.                     }
  255.                     OGRPolygon polygon; // 多边形
  256.                     polygon.addRing(&ogrring);
  257.                     poFeature->SetGeometry(&polygon);
  258.                     // 创建字段
  259.                     poFeature->SetField("KJFWXXID", "7CF7F6D9-9210-4BC7-AFAD-F288F2357A99");
  260.                     poFeature->SetField("YSDM", "8008010240");
  261.                     poFeature->SetField("DZJGH", "/");
  262.                     poFeature->SetField("YDMJ", polygon.get_Area());
  263.                     poFeature->SetField("JTGCXMLX", "1");
  264.                     poFeature->SetField("DLDJ", "03");
  265.                     poFeature->SetField("DLCD", 0.0);
  266.                     poFeature->SetField("BZHDMKD", 0.0);
  267.                     poFeature->SetField("GDJTLX", "");
  268.                     poFeature->SetField("GDCD", 0.0);
  269.                     poFeature->SetField("KZXKD", 0.0);
  270.                 }
  271.                 else if (_tcscmp(szKword, _T("B")) == 0) // 管线工程
  272.                 {
  273.                     OGRLineString pLine; // 线
  274.                     int pNum = (int)ogrPts.size();
  275.                     pLine.setNumPoints(pNum);
  276.                     for (int j = 0; j < ogrPts.size(); j++)
  277.                     {
  278.                         pLine.setPoint(j, ogrPts[j].getX(), ogrPts[j].getY(), ogrPts[j].getZ());
  279.                     }
  280.                     poFeature->SetGeometry((OGRGeometry *)(&pLine));
  281.                     // 创建字段
  282.                     poFeature->SetField("KJFWXXID", "7CF7F6D9-9210-4BC7-AFAD-F288F2357A99");
  283.                     poFeature->SetField("YSDM", "8008010250");
  284.                     poFeature->SetField("DZJGH", "/");
  285.                     // GS 给水   PS 排水   TX 通信   DL 电力
  286.                     // RQ 燃气   GR 供热   XF 消防   QT 其他
  287.                     poFeature->SetField("SPGXLX", "GS");
  288.                     poFeature->SetField("GXFJ", "");
  289.                     poFeature->SetField("SPGXCD", 0.0);
  290.                     poFeature->SetField("JGGXCD", dMaxLength);
  291.                     poFeature->SetField("SPGXGJ", "");
  292.                     poFeature->SetField("JGGXGJ", "");
  293.                 }

  295.                 if (poLayer->CreateFeature(poFeature) != OGRERR_NONE)
  296.                 {
  297.                     acutPrintf(_T("\n在 shapefile 中创建特性失败!\n"));
  298.                     OGRFeature::DestroyFeature(poFeature);
  299.                     GDALClose(poDS);
  300.                     poDS = NULL;
  301.                     vertexPoint.clear();
  302.                     return;
  303.                 }
  304.                 // 资源清理
  305.                 OGRFeature::DestroyFeature(poFeature);
  306.             }
  307.             // 资源清理
  308.             GDALClose(poDS);
  309.             poDS = NULL;
  310.             vertexPoint.clear();
  311.             acutPrintf(_T("\n已输出shp文件格式!\n"));
  312.         }
  313.     }
  314. }


四 代码示例：读shp文件
  1. // 打开文件对话框
  2. CFileDialog dlg(TRUE, NULL, NULL, OFN_PATHMUSTEXIST | OFN_FILEMUSTEXIST, TEXT("shp文件(*.shp)|*.shp||"), NULL);
  3. dlg.GetOFN().lpstrTitle = _T("读取shp数据文件");
  4. if (dlg.DoModal() == IDOK)
  5. {
  6.     CString strPath = dlg.GetPathName();
  7.     USES_CONVERSION;  //需要这一步才能使用W2A
  8.     const char* pszShapeFile = string(W2A(strPath)).c_str();

  10.     GDALAllRegister();
  11.     GDALDataset *poDS;
  12.     CPLSetConfigOption("SHAPE_ENCODING", "");  //解决中文乱码问题
  13.     // 读取shp文件
  14.     poDS = (GDALDataset*)GDALOpenEx(pszShapeFile, GDAL_OF_VECTOR, NULL, NULL, NULL);
  15.     if (poDS == NULL)
  16.     {
  17.         acutPrintf(_T("Open failed.\n%s"));
  18.         return;
  19.     }
  20.     OGRLayer *poLayer;
  21.     poLayer = poDS->GetLayer(0); // 读取层
  22.     OGRFeature *poFeature;
  23.     poLayer->ResetReading();
  24.     int i = 0;
  25.     while ((poFeature = poLayer->GetNextFeature()) != NULL)
  26.     {
  27.         if (poFeature->GetFieldAsDouble("AREA") < 1)
  28.             continue; //去掉面积过小的 polygon
  29.         i = i++;
  30.         OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn();
  31.         int iField;
  32.         // 获得字段的数目，不包括前两个字段（FID,Shape);
  33.         int n = poFDefn->GetFieldCount();
  34.         for (iField = 0; iField < n; iField++)
  35.         {
  36.             // 输出每个字段的值
  37.             acutPrintf(_T("%s\n"), poFeature->GetFieldAsString(iField));
  38.         }
  39.         OGRFeature::DestroyFeature(poFeature);
  40.     }
  41.     GDALClose(poDS);
  42. }






本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

x
回复

使用道具 举报

ymgupl
发表于 2025-8-16 14:46:47 | 显示全部楼层
G佬NB              
回复 支持 反对

使用道具 举报

yonjay
发表于 2025-8-18 20:33:38 | 显示全部楼层
G大，SHP字段、长度、类型可以写到配置中，一个接口直接完成
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|手机版|CAD论坛|CAD教程|CAD下载|联系我们|关于明经|明经通道 ( 粤ICP备05003914号 )  
©2000-2023 明经通道 版权所有 本站代码，在未取得本站及作者授权的情况下，不得用于商业用途

GMT+8, 2025-10-6 01:40 , Processed in 0.188497 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表