OGR库读写mif/tab文件lonefox使用经验小结

来源：互联网发布：数据库架构师北京编辑：程序博客网时间：2024/05/16 14:00

===================================

本人log中所有未注明转载的文章和blog一般为本人原创或整理加工，
原创文章版权本人（lonefox）所有；转载文章版权归原作者所有；

http://blog.csdn.net/boythl

欢迎转载，但请注明出处，保留作者和版权信息。

===================================

最近项目中使用GDAL/OGR库读写mif和tab文件, 经过几天的折腾,积累了点点经验,且记之.

1. OGRSFDriverRegistrar::Open(Path, FALSE);只能读文件

对于OGRSFDriverRegistrar::Open方法, ogr内置的mif/TAB读写库MITAB是不支持该方法的更新的,看下面的MITAB源码:

OGRDataSource *OGRTABDriver::Open( const char * pszFilename, int bUpdate ){ OGRTABDataSource *poDS; if( bUpdate )//请注意这儿,如果bUpdate为TRUE,则直接返回NULL { return NULL; } //other code}

再深入源码，你会发现，通过Open方法传给IMapInfoFile类访问文件的方式是TABRead，而非TABWrite

也许你会说"这是OGRTABDriver的方法, 可以看看OGRMIFDriver的方法?", 遗憾, 我们只能用 RegisterOGRTAB(); 注册mif和tab通用的驱动,RegisterOGRMIF();是一个只有定义而没有函数实现的美丽谎言, 不知道OGR库留着这个定义的目的何在。也不存在OGRMIFDriver，OGRTABDriver通吃mif和tab文件。这也是我们的第二点。

2.只使用RegisterOGRTAB();注册tab/MIF文件读写驱动

原因上面已经提过。单凭RegisterOGRTAB();已足够。

mif和tab文件的读取都是通过Open方法，MITAB库内部会调用SmartOpen方法区分是mif/mid后缀或是tab后缀。

3.写mif/tab文件的方法

OGRSFDriver *pDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName("MapInfo File");
OGRDataSource *pDSTAB = pDriver->CreateDataSource(strTabPath);//创建一个文件，根据扩展名自动判断是创建mif文件或tab文件

lonefox的解释： MITAB库在实现CreateDataSource方法时会调用OGRTABDataSource的Create方法，而非Open方法

Create方法内部会调用IMapInfoFile类的Open方法，但区别在于访问模式参数传递的是"wb"！

而该模式将标记 m_eAccessMode = TABWrite; 不是OGRDataSource::Open方法对应的TABRead模式；

顺带提句，IMapInfoFile内只有TABWrite和TABRead, 是不能又读又写的

typedef enum{ TABRead, TABWrite, TABReadWrite /* ReadWrite not implemented yet *///注意这个注释} TABAccess;

4.读写mif/tab的方法

如果我一个文件已有部分数据，但我想打开后继续加数据，怎么办？

我采取一种曲线救国的方法解决该问题： i. 打开该文件; ii.按上面第三条的方法新建一个临时文件供写入； iii. 使用OGRDataSource类的CopyLayer方法将 i 中打开的旧数据贴进新layer，iv. 调用OGRLayer::SyncToDisk方法保存新文件；v. 删除旧文件将临时文件更名

若有大侠有更好的方法，一定要不吝指教才好！！！

5. mif/tab格式互转

利用OGRDataSource::CopyLayer，mif和tab格式的互转就再简单不过了。直接源码吧

BOOL MIF2TAB( const CString& MIFPath, const CString& TabPath ){RegisterOGRTAB();OGRSFDriver *pDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName("MapInfo File");if( pDriver == NULL ){return FALSE;}//打开原始mif//CNILayer::CreateEmptyMID(MIFPath);//将空mid改非空,避免ogr加载失败OGRDataSource* pDSMIF = OGRSFDriverRegistrar::Open(MIFPath);if (pDSMIF == NULL){TRACE1("warning:MIF file open failed %s/n", MIFPath);return FALSE;}OGRLayer* pLayerMIF = pDSMIF->GetLayer(0);//创建结果TABOGRDataSource *pDSTAB = pDriver->CreateDataSource(TabPath);if (pDSTAB == NULL){TRACE1("error:TAB file created failed %s/n", TabPath);OGRDataSource::DestroyDataSource( pDSMIF );return FALSE;}//拷贝数据OGRLayer* pLayerTAB = pDSTAB->CopyLayer(pLayerMIF, "test", NULL);//关闭MIF文件OGRDataSource::DestroyDataSource(pDSMIF);if( pLayerTAB == NULL ){TRACE2("error:copied layer to TAB failed %s->%s/n", MIFPath, TabPath);OGRDataSource::DestroyDataSource( pDSTAB );return FALSE;}//保存结果OGRErr err = pLayerTAB->SyncToDisk();//关闭TAB文件OGRDataSource::DestroyDataSource( pDSTAB );if (err != OGRERR_NONE){TRACE1("error:saved TAB file failed %s/n", TabPath);return FALSE;}return TRUE;}

tab转mif和上面代码雷同，略之。

6. ogr库的bug

6.1. 读取空mif文件失败的限制

如果一个mif文件是空文件（当然，文件格式是符合mif格式的，只是图层中没有feature），即mid文件尺寸为0；那么，ogr的open方法会返回NULL，意味着读取失败。

这是因为MITAB库在打开mid文件时会尝试读一行

/* Put the MID file at the correct location, on the first feature */ if (m_eAccessMode == TABRead && (m_poMIDFile->GetLine() == NULL))//空mid时GetLine会返回NULL { Close(); if (bTestOpenNoError) CPLErrorReset(); return -1; }

这段代码位于int MIFFile::Open(const char *pszFname, const char *pszAccess, GBool bTestOpenNoError /*=FALSE*/ )方法内部。

有两种方法解决该问题： 1. 修改ogr库源码，先前我们是这么干的；2. 打开前向空mid文件插入一个换行符，该字符的存在能“欺骗”上文GetLine方法不返回NULL

6.2. 字符集的限制

ogr写出的mif/tab文件，字符集全部都是Neutral，根本不管先前的charset，鄙人驽钝，未发现修改字符集的方法，有大虾知道的话请千万要不吝指教啊！

但在IMapInfoFile类中存在一个SetCharset接口，该方法由IMapInfoFile默认实现，未被override；且该方法在MITAB库内部未被调用。

问题是OGR库未提供类似接口。

解决方法： 1. 修改ogr库，增加一个设置charset的接口。2. 抛弃OGR库，直接使用IMapInfoFile类操纵mif文件

6.3. 只支持少量sql查询，不支持空间查询

sql查询的相关内容可以看swq.h文件；空间查询在OGR1.60版还未看到接口，期待。

7. 尽量使用GetNextFeature代替GetFeatureCount

我们看看OGRLayer::GetFeatureCount方法的实现代码

int OGRLayer::GetFeatureCount( int bForce ){ OGRFeature *poFeature; int nFeatureCount = 0; if( !bForce ) return -1; ResetReading(); while( (poFeature = GetNextFeature()) != NULL ) { nFeatureCount++; delete poFeature; } ResetReading(); return nFeatureCount;}

由上可以看到，大部分情况下，GetFeatureCount的参数取False时将直接得到-1；如果采用TRUE，不好意思，就是一次整图层所有feature的遍历。

很遗憾的是，MITAB库并未override该方法，也就没有更高效的算法了。

因此，如果想判断图层是否有feature，例如sql查询的查询结果是否为空，不明确要求得到要素个数时，请先使用ResetReading再判断GetNextFeature是否为NULL，避免使用GetFeatureCount。

切记，OGR的GetFeatureCount不同于ado库的GetRecordCount。

时间匆忙，lonefox上班间隙随手写就，再加就几天的使用，经验有限，错误之处恳请斧正；另，GIS同行欢迎加好友共同探讨！