перепроецирование gdalwrap

[XdenisX]

Коллеги, очень нужна ваша помощь. Пытаюсь конвертировать с помощью gdalwrap геотиф файл со следующей мета информацией:

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

gdalinfo k1.tif
Driver: GTiff/GeoTIFF
Files: k1.tif
Size is 1152, 1150
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Polar_Stereographic"],
    PARAMETER["latitude_of_origin",90],
    PARAMETER["central_meridian",45],
    PARAMETER["scale_factor",0.994],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (277823.000000000000000,-1291703.000000000000000)
Pixel Size = (994.000000000000000,-994.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  277823.000,-1291703.000) ( 57d 8'18.28"E, 78d 8'25.01"N)
Lower Left  (  277823.000,-2434803.000) ( 51d30'34.47"E, 68d11'2.53"N)
Upper Right ( 1422911.000,-1291703.000) ( 92d46'1.84"E, 72d48'55.16"N)
Lower Right ( 1422911.000,-2434803.000) ( 75d18'8.02"E, 64d59'17.00"N)
Center      (  850367.000,-1863253.000) ( 69d31'53.06"E, 71d42'11.32"N)
Band 1 Block=1152x2 Type=Byte, ColorInterp=Red
Band 2 Block=1152x2 Type=Byte, ColorInterp=Green
Band 3 Block=1152x2 Type=Byte, ColorInterp=Blue

Конвертирую следующей строкой:

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

gdalwrap -s_srs "+proj=stere +lat_0=90 +lon_0=45 +ellps=sphere +x_0=0 +y_0=0 +units=m +no_defs" -t_srs "EPSG:3576" k1.tif k111.tif

Получается вот что:

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

gdalinfo k111.tif
Driver: GTiff/GeoTIFF
Files: k111.tif
Size is 1660, 1628
Coordinate System is:
PROJCS["WGS 84 / North Pole LAEA Russia",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","3576"]]
Origin = (-1503783.547560292761773,-1108700.852565036620945)
Pixel Size = (961.480047733915285,-961.480047733915285)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-1503783.548,-1108700.853) 
Lower Left  (-1503783.548,-2673990.370) 
Upper Right (   92273.332,-1108700.853) 
Lower Right (   92273.332,-2673990.370) 
Center      ( -705755.108,-1891345.611) 
Band 1 Block=1660x1 Type=Byte, ColorInterp=Red
Band 2 Block=1660x1 Type=Byte, ColorInterp=Green
Band 3 Block=1660x1 Type=Byte, ColorInterp=Blue

Когда гружу в ArcMap он говорит, что в файле потеряна пространственная информация, geoserver тоже не может создать coverage.xml

В чем дело не понимаю

Спасибо.

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

Небольшая поправка: gdalwarp От англ-го warp.

По сути: вы уверены, что параметры в "-s_srs" правильные? Например, почему "+ellps=sphere", хотя во входном геотифе WGS84? Можно попробовать сделать без указания входной проекции, только с выходной:

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

gdalwarp -t_srs "EPSG:3576" k1.tif k111.tif

[Voltron]

А зачем вообще -s_srs узакана, если в исходном файле и так присутсвует привязка?

[XdenisX]

А зачем вообще -s_srs узакана, если в исходном файле и так присутсвует привязка?

Я пробовал не указывать параметры исходной проекции - но результат тот же

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

Интересно: размер входного растра -- 1152x1150, выходного -- 1660x1628..

[gamm]

Я пробовал не указывать параметры исходной проекции - но результат тот же

укажите выходную проекцию в "+" формате, без ссылок на EPSG - видимо, чего-то в базе не хватает

[gamm]

Интересно: размер входного растра -- 1152x1150, выходного -- 1660x1628..

могдо развернуть - в полярной стереографической Россия "на попа" стоит, если на 45 градусов =>sqrt(2), то примерно так и будет, даже без изменения размера пикселя

[XdenisX]

Я пробовал не указывать параметры исходной проекции - но результат тот же

укажите выходную проекцию в "+" формате, без ссылок на EPSG - видимо, чего-то в базе не хватает

Спасибо, именно это помогло
Прописал в -t_srs "+proj=laea +lat_0=90 +lon_0=90 +ellps=WGS84 +datum=WGS84 +units=m"