GDAL读取Shapefile

最新推荐文章于 2023-05-11 22:27:13 发布

weixin_30797027

最新推荐文章于 2023-05-11 22:27:13 发布

阅读量611

点赞数

原文链接：http://www.cnblogs.com/gispathfinder/p/10890185.html

版权

--------------------------------------------------------------------------------------

#include <gdal_priv.h>
#include <ogrsf_frmts.h>

GDALAllRegister();
OGRDataSource* poDS;
const char* pszDriverName = "ESRI Shapefile";
//
GDALDriverManager* pDriverManager = GetGDALDriverManager(); //new GDALDriverManager();
GDALDriver* pDriver = pDriverManager->GetDriverByName(pszDriverName);
if (pDriver == nullptr)
{
  qDebug() << pszDriverName << " driver not available.";            
  return;
}
//
const char* filePathChar1 = "D:\\data";
OGRSFDriver* pSFDriver = static_cast<OGRSFDriver*>(pDriver);
OGRDataSource* pDataSource = pSFDriver->Open(filePathChar, FALSE);

/*
OGRSFDriverRegistrar* pDriverRegistrar = OGRSFDriverRegistrar::GetRegistrar();
OGRSFDriver* pSFDriver = (OGRSFDriver*)pDriverRegistrar->GetDriverByName(pszDriverName);
OGRDataSource* pDataSource = pSFDriver->Open(filePathChar, FALSE);
*/

------------------------------------------------------------------------------------------------


OGRLayer* pLayer = pDataSource->GetLayerByName(baseNameChar);
/*const char* geometryColumnName = pLayer->GetGeometryColumn();
int fieldIndex = pLayer->FindFieldIndex(geometryColumnName, TRUE);*/
// pEnvelopeDisplay
if (pLayerList.size() == 0) {
  OGREnvelope* pEnvelope = new OGREnvelope();
  pLayer->GetExtent(pEnvelope, TRUE);
  double centerX = (pEnvelope->MinX + pEnvelope->MaxX) / 2;
  double centerY = (pEnvelope->MinY + pEnvelope->MaxY) / 2;
  double distanceX = pEnvelope->MaxX - pEnvelope->MinX;
  double distanceY = pEnvelope->MaxY - pEnvelope->MinY;
  //
}