ESRI WKID 102100

[ericsson]

Столкнулся с одним сервисом, где СК определена через ESRI WKID 102100, пытаюсь понять, что это, вроде как должен быть аналог гугловскому меркатору на сфере. На spatialreference.org нашел вот такой PRJ, соответствующий этой СК:

Код: Выделить всё

PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",
GEOGCS["GCS_WGS_1984",
DATUM["D_WGS_1984",
SPHEROID["WGS_1984",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Mercator_Auxiliary_Sphere"],
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",0.0],
PARAMETER["Standard_Parallel_1",0.0],
PARAMETER["Auxiliary_Sphere_Type",0.0],
UNIT["Meter",1.0]]

Похоже на правду?

[Александр Мурый]

Пруф:

. . .
ESRI:102113 - First used by ESRI, but then replaced by 102100 (equivalent, deprecated)
ESRI:102100 - Built differently, but equivalent to 102113 (equivalent, deprecated)
. . .
?? - ESRI ArcGIS 10 changed the WKID for this projection from ESRI:102100 to EPSG:3857 [8]

По-моему, сходство налицо
Отсюда:

Код: Выделить всё

PROJCS["unnamed",
    GEOGCS["unnamed ellipse",
        DATUM["unknown",
            SPHEROID["unnamed",6378137,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Mercator_2SP"],
    PARAMETER["standard_parallel_1",0],
    PARAMETER["central_meridian",0],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["Meter",1]

Единственное, чётко не определён WGS84, но это не мешает PROJ4 работать.

[ericsson]

Тем мне менее, запрошенный кусок куда-то на юг почти на 20 км улетел. Но X-координата вполне себе в порядке.

[Александр Мурый]

"Улёты" случаются, подтверждаю

Полная строка 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

Если я верно понимаю, на желание данных улететь влияет параметр "+nadgrids=@null"

[ericsson]

Поскольку параметр вводить мне некуда (я банальным REST'ом дергаю растры), обошел проблему используя bboxSR и imageSR равные 32637, то есть WGS84-UTM37N - там все совпало.