Сохранение информации о проекции
Добавлено: 19 янв 2014, 11:37
Предположим, у нас есть слой в EPSG:4326, а мы хотим перевести его в сферический меркатор EPSG:3857. Возьмем к примеру слой admin из набора Geosample из БД PostGIS. Нажимаем правой кнопкой на имени слоя и выбираем Save As и указываем проекцию EPSG:3857. Получаем на выходе Shape файл с двумя файлами с описанием проекции *.prj и *.qpj. Содержимое первого:
Второго:
Не вникая в подробности форматов, предположу, что в первом случае проекция описана как меркатор на эллипсоиде, а во втором (по крайней мере в строке proj) как меркатор на сфере. Почему так происходит? Я столкнулся с ситуацией, когда ogr, читая информацию о проекции такого слоя интерпретирует его как меркатор на эллипсоиде.
Код: Выделить всё
PROJCS["WGS_84_Pseudo_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Mercator"],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],PARAMETER["standard_parallel_1",0.0]]Код: Выделить всё
PROJCS["WGS 84 / Pseudo-Mercator",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],PROJECTION["Mercator_1SP"],PARAMETER["central_meridian",0],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"],AXIS["X",EAST],AXIS["Y",NORTH],AUTHORITY["EPSG","3857"]]