2017-01-25 13 views
2

Использование импорта gdal, ogr, и заданного EPSG (например, здесь 2154), что такое код в python для записи (почти) пустого файла TAB MapInfo или просто сводной строки проекции (формат TAB):Python GDAL/OGR: как создать файл TAB 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?

ответ

2

Найдено решение:

print proj.ExportToMICoordSys() 

это дает:

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

 Смежные вопросы

  • Нет связанных вопросов^_^