преобразования системы координат при помощи gdalwarp

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

Офиц.код Google Mercator EPSG:3857. Попробуйте с ним вместо EPSG:900913

[malo-lilo]

Команды я запускаю через FWTools.
На EPSG:3857 заменила, но тоже не помогло.
Может я изначально что-то не так делаю? Хотя вроде следую инструкции, вот этой http://ru.topoaxis.com/blog/111/

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

Приведите (если это не секрет) вывод команды gdalinfo some.jpg

[malo-lilo]

Нет, не секрет.

Driver: JPEG/JPEG JFIF
Files: data\need\some.jpg
data\need\some.jpg.aux.xml
Size is 9856, 8400
Coordinate System is:
PROJCS["Popular Visualisation CRS / Mercator",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,
AUTHORITY["EPSG","7059"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6055"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4055"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","3785"],
AXIS["X",EAST],
AXIS["Y",NORTH]]
Origin = (5802478.136555763900000,7686822.867098996400000)
Pixel Size = (19.109257071294099,-19.109257071294099)
Metadata:
PyramidResamplingType=NEAREST
Image Structure Metadata:
SOURCE_COLOR_SPACE=YCbCr
INTERLEAVE=PIXEL
COMPRESSION=JPEG
Corner Coordinates:
Upper Left ( 5802478.137, 7686822.867) ( 52d 7'28.37"E, 56d38'22.79"N)
Lower Left ( 5802478.137, 7526305.108) ( 52d 7'28.37"E, 55d50'18.10"N)
Upper Right ( 5990818.974, 7686822.867) ( 53d48'59.19"E, 56d38'22.79"N)
Lower Right ( 5990818.974, 7526305.108) ( 53d48'59.19"E, 55d50'18.10"N)
Center ( 5896648.555, 7606563.987) ( 52d58'13.78"E, 56d14'27.99"N)
Band 1 Block=9856x1 Type=Byte, ColorInterp=Red
Image Structure Metadata:
COMPRESSION=JPEG
Band 2 Block=9856x1 Type=Byte, ColorInterp=Green
Image Structure Metadata:
COMPRESSION=JPEG
Band 3 Block=9856x1 Type=Byte, ColorInterp=Blue
Image Structure Metadata:
COMPRESSION=JPEG

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

Выглядит вроде нормально.. А если попробовать EPSG:3785 ? Такой код написан для проекции Popular Visualisation CRS / Mercator

Или, на худой конец, можно попробовать так (без указания исходной проекции):
gdalwarp -t_srs EPSG:4326 some.jpg foo1.tif

[malo-lilo]

Спасибо вам за помощь, amuriy!
Файл я нашла, он оказался в какой-то богом забытой папке с непонятным адресом.. Правда этот tif oказался сильно растянутым.
Теперь у меня другая проблема: listgeo не хочет мне создавать world файл для foo1.tif.
Программа пишет что команда должна выглядеть так:
Usage: listgeo [-norm] [-t tabledir] [filename]

А в статье в которой описан алгоритм получения геопривязанного изображения с Google Maps написана вот такая команда:
listgeo –tfw foo.tif

[Voltron]

Зачем пользоваться левой утилитой, если gdalwarp может создавать world-файл при выполнении перепроецирования?

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

gdalwarp -s_srs EPSG:900913 -t_srs EPSG:4326 -co TFW=YES C:\data\need\some.jpg C:\data\need\foo1.tif

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

Voltron, listgeo отнюдь не "левая", а, наоборот, очень родная геотиффовская утилита. Другое дело, что статья корявая, лучше читать родную документацию, там же на чистом языке Шекспира всё написано В статье не сказано, что для получения tfw-файла вывод listgeo надо перенаправить в файл:

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

listgeo –tfw foo.tif > foo.tfw

[malo-lilo]

Ура! И еще раз ура, с вашей помощью amuriy и Voltron я наконец-то справилась с поставленной задачей, получила геопривязанную картинку и теперь могу продолжать работу дальше.
Спасибо вам большое:)
А статья действительно корявая, там даже вместо listgeo написано lsitgeo...