2
import gdal, ogr, 을 사용하고 EPSG (예 : 여기 2154)를 주면 Python의 (거의) 빈 MapInfo TAB 파일 또는 투영의 요약 행 (TAB 형식)을 작성하는 코드는 무엇입니까Python GDAL/OGR : MapInfo TAB 파일을 만드는 방법은 무엇입니까?
CoordSys Earth Projection 3,999,0,0,0,0,7,3,46.5,44.0,49.0,700000,6600000
모든 EPSG 코드는 무엇입니까? 탭용 파일에 대한 몇 가지 상응 있습니까
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"]]
:
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에 대한 글을 참고하세요 :
는 ESRI와 함께, 다음 줄은 .PRJ 파일을 만드는 것이 나란히하려면?
print proj.ExportToMICoordSys()
이 제공 :