Использование настроенной системы координат в ArcGIS Desktop?

16

Я не так много знаю о системах координат ...

В моем офисе мы используем пространственные данные, поступающие с археологических памятников. Каждый сайт имеет свою собственную систему координат XYZ. Три простые ортогональные декартовы оси . В последние годы мы управляли этими пространственными данными с помощью программного обеспечения ГИС (ArcGIS), без использования конкретной системы координат (просто оставьте ее как «неопределенная»)

Я хотел бы знать, существует ли какая-либо GCS, предназначенная для работы с такими наборами данных с использованием простой декартовой ортогональной оси, без искажений сетки типичной GCS. Кроме того, я хотел бы знать, подходит ли эта система для использования ее в онлайн-картографическом приложении.

Кстати, мы управляем 2D (ArcMap) и 3D (ArcScene) средами и работаем с "мм" в качестве базовой единицы длины.

Если такой вещи не существует, возможно, кто-то знает, как ее создать.

Maik
источник

Ответы:

17

Предполагая, что ваши области интереса довольно малы по сравнению с земным шаром, вы можете настроить собственную поперечную проекцию Меркатора.

Вы должны знать географические координаты lat_0 и lon_0 происхождения вашего CRS и направление осей 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 теперь могут работать с локальными проекциями, используя мм в качестве единиц измерения.

Andrej
источник