Перепривязка растров и нестандартные СК

[Alt]

Здравствуйте.

Есть несколько геопривязанных растров с world-файлами и файлом проекции:

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

PROJCS["Custom",
       GEOGCS["GCS_Sphere",
                    DATUM["D_Sphere",SPHEROID["Sphere",6371000,0]],
                    PRIMEM["Greenwich",0],
                    UNIT["Degree",0.0174532925199433]],
       PROJECTION["Mercator"],
       PARAMETER["False_Easting",0],
       PARAMETER["False_Northing",0],
       PARAMETER["Central_Meridian",0],
       PARAMETER["Standard_Parallel_1",0],
       UNIT["Meter",1]] 

В списках EPSG не нашел аналогов.

Что можно предпринять чтобы перепроецировать эти растры в проекции Гаусса-Крюгера ?
В примерно такую:

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

PROJCS["Pulkovo 1942 / 3-degree Gauss-Kruger zone 10",
            GEOGCS["Pulkovo 1942",
                         DATUM["Pulkovo_1942",SPHEROID["Krassowsky 1940",6378245,298.3,AUTHORITY["EPSG","7024"]], AUTHORITY["EPSG","6284"]],
PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4284"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",30],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",10500000],
PARAMETER["false_northing",0],
UNIT["metre",1,AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","2526"]]"

Благодарю.

[Максим Дубинин]

Не надо искать обязательно аналоги в EPSG, она же не перебирает все возможные сочетания
Если у вас есть prj, то просто делаете примерно так:

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

gdalwarp -s_srs ESRI::Input.prj ESRI::output.prj -t_srs input.tif output.tif

[Alt]

Благодарю.
А как поступить если понадобится указать эту проекцию для растрового слоя Mapserver ?

[Максим Дубинин]

1. создаете описание СК в формате proj4
2. прописываете его в epsg, например как рассказано ниже по тексту вот здесь
http://gis-lab.info/qa/ms-gm.html
3. используете в map-файле