假设与地球相比,您感兴趣的区域很小,则可以设置自定义横向墨卡托投影。
您必须知道CRS原点的地理坐标lat_0和lon_0,以及x和y轴的方向:
+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
要与Arcgis一起使用,.prj文件应包含:
PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",51.4],PARAMETER["central_meridian",7],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
如果它们与北,东方向不平行,则需要使用倾斜的墨卡托投影进行旋转:
+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
此投影的.prj文件为:
PROJCS["Hotine_Oblique_Mercator_Azimuth_Center",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",51.4],PARAMETER["longitude_of_center",7],PARAMETER["azimuth",-10],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
这两个都在一起(蓝色的tmerc和红色的omerc):
lat_0和lon_0(或lonc)是您本地CRS原点的坐标。您必须使用GPS或其他地图来探索这些。您必须根据本地X轴上两个已知点的坐标来计算角度。
唯一无法正常工作的是单位毫米。它应该工作,但在当前的QGIS中不起作用:-(
http://hub.qgis.org/issues/9414
更新资料
当前的QGIS版本现在可以使用以mm为单位的局部投影。