文章/答案/技术大牛

发布

社区首页 >问答首页 >Python GDAL/OGR :如何创建MapInfo TAB文件？

问Python GDAL/OGR :如何创建MapInfo TAB文件？
EN

Stack Overflow用户

提问于 2017-01-25 17:38:05

回答 1查看 1.5K关注 0票数 2

使用import gdal、ogr和给定的EPSG (例如，这里是2154)，用python编写(几乎)空的MapInfo标签文件的代码是什么，或者仅仅是投影的摘要行(标签格式)：

CoordSys Earth Projection 3,999,0,0,0,0,7,3,46.5,44.0,49.0,700000,6600000

任何EPSG代码？

为了将其与ESRI并列，以下行创建一个.prj文件：

from osgeo import osr

def make_prj(epsg, path_out):
    proj=osr.SpatialReference()
    proj.ImportFromEPSG(int(epsg))
    with open(path_out, "w") as out:
        print >>out, proj.ExportToWkt()

对于epsg=2154，这是这样写的：

PROJCS["RGF93 / Lambert-93",GEOGCS["RGF93",DATUM["Reseau_Geodesique_Francais_1993",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6171"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4171"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["standard_parallel_1",49],PARAMETER["standard_parallel_2",44],PARAMETER["latitude_of_origin",46.5],PARAMETER["central_meridian",3],PARAMETER["false_easting",700000],PARAMETER["false_northing",6600000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","2154"]]

TAB文件有没有等价物？

mapinfo

epsg

python

gdal

ogr

回答 1

Stack Overflow用户

发布于 2017-01-25 20:11:57

找到解决方案：

print proj.ExportToMICoordSys()

它提供了：

Earth Projection 3, 104, "m", 3, 46.5, 49, 44, 700000, 6600000

票数 3

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/41848065

复制

相似问题

问Python GDAL/OGR :如何创建MapInfo TAB文件？
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问Python GDAL/OGR :如何创建MapInfo TAB文件？EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问Python GDAL/OGR :如何创建MapInfo TAB文件？
EN