ошибки gdalwarp

Ответить
progressive
Новоприбывший
Сообщения: 3
Зарегистрирован: 11 апр 2012, 15:13
Репутация: 0

ошибки gdalwarp

Сообщение progressive » 11 апр 2012, 15:33

добрый день.
я сверх новичок). в деле буквально 2 дня.

передо мной стоит задача сделать красивую карту. точно такую как тут http://mapbox.com/tilemill/docs/guides/terrain-data/

но вот беда.. в ввиду практически полного отсутствия знаний в этой области, я завис буквально на первых шагах. вот что выдает.
ERROR 1: Too many points (2116 out of 2116) failed to transform, unable to compute output bounds.
почитав форум чуть чуть, я понял что в обнаружении проблемы может помочь gdalinfo. вот мой лог-

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

Driver: GTiff/GeoTIFF
Files: mos.tif
Size is 2582, 2242
Coordinate System is:
PROJCS["LAEA        D000",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",48],
    PARAMETER["longitude_of_center",48],
    PARAMETER["false_easting",48],
    PARAMETER["false_northing",48],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (-1668800.000000000000000,2724700.000000000000000)
Pixel Size = (2000.000000000000000,-2000.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (-1668800.000, 2724700.000) (  5d46' 4.89"E, 67d56'12.88"N)
Lower Left  (-1668800.000,-1759300.000) ( 30d39'47.49"E, 30d30'31.93"N)
Upper Right ( 3495200.000, 2724700.000) (115d41' 4.79"E, 56d21'28.68"N)
Lower Right ( 3495200.000,-1759300.000) ( 83d12'53.93"E, 25d16'46.23"N)
Center      (  913200.000,  482700.000) ( 61d15'33.88"E, 51d37'22.81"N)
Band 1 Block=2582x3 Type=Byte, ColorInterp=Red
  Minimum=0.000, Maximum=255.000, Mean=82.814, StdDev=53.050
  Metadata:
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=82.814132493465
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=53.049521862569
Band 2 Block=2582x3 Type=Byte, ColorInterp=Green
  Minimum=0.000, Maximum=255.000, Mean=84.158, StdDev=64.858
  Metadata:
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=84.157859669392
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=64.858199711693
Band 3 Block=2582x3 Type=Byte, ColorInterp=Blue
  Minimum=0.000, Maximum=255.000, Mean=89.951, StdDev=34.481
  Metadata:
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=89.951196992007
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=34.481184592024
ребят, у кого есть время. помогите, на пальцах объясните что к чему.
Последний раз редактировалось Voltron 11 апр 2012, 15:40, всего редактировалось 1 раз.
Причина: не забываем использовать теги [code] и [quote]

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: ошибки gdalwarp

Сообщение Александр Мурый » 11 апр 2012, 15:37

Какая опер.система, какая версия GDAL? Какие конкретно шаги из статьи надо выполнить? В ответ на какую команду выдаётся ошибка?

Желательно прикрепить сам файл.
Редактор материалов, модератор форума

Voltron
Гуру
Сообщения: 2627
Зарегистрирован: 29 мар 2007, 14:12
Репутация: 34
Откуда: Ukraine

Re: ошибки gdalwarp

Сообщение Voltron » 11 апр 2012, 15:40

Вывод gdalinfo не имеет смысла без информации о том, какая команда у вас выдает ошибку

progressive
Новоприбывший
Сообщения: 3
Зарегистрирован: 11 апр 2012, 15:13
Репутация: 0

Re: ошибки gdalwarp

Сообщение progressive » 11 апр 2012, 16:04

я на mac os 10.7
GDAL 1.9.0

кстати да, насчет файлов. небольшое лирическое отступление: я работаю в винодельческой отрасли, в Молдове, и мне нужен конкретно южный ее регион.
я написал в гугле moldova geotiff и скачал воот отсюда http://www.eea.europa.eu/data-and-maps/ ... del-etopo5

ошибку выдает тут:

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

gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3785 -r bilinear mos.tif  mos-3785.tif
ERROR 1: Too many points (2116 out of 2116) failed to transform, unable to compute output bounds.
там на сайте откуда я качал, есть строка Coordinate reference system: EPSG:4008
я менял в команде с EPSG:4326 на EPSG:4008 но результат тот же.

мне нужно добиться такой картинки, которая получается в конце примера.
Последний раз редактировалось Voltron 11 апр 2012, 16:53, всего редактировалось 1 раз.
Причина: для выделения команд используем тег [code]

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: ошибки gdalwarp

Сообщение Александр Мурый » 11 апр 2012, 16:26

По ссылке ETOPO5 -- глобальная (!) модель рельефа, к тому же это слишком грубый рельеф для ваших целей, мне кажется.

Во-первых, рекомендую скачать разграфку для данных SRTM (разрешение 90м), выложенных на гис-лабе. Затем скачать нужные данные отсюда.

Если нужен ASTER GDEM с разрешением 30м, то советую пользоваться сервисом Reverb, где есть ASTER GDEM 2-й (улучшенной) версии. На этом сайте надо зарегистрироваться, далее нужно выделить нужный участок и указать нужные данные для поиска --> заказать данные.
Reverb_Moldova.png
Reverb_Moldova.png (303.84 КБ) 8271 просмотр
И SRTM, и ASTER GDEM v2 идут в формате GeoTIFF, так что дальше можно действовать по статье :)
Редактор материалов, модератор форума

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Репутация: 793
Ваше звание: званий не имею
Откуда: Москва

Re: ошибки gdalwarp

Сообщение Александр Мурый » 11 апр 2012, 16:37

progressive, написал в личку.
Редактор материалов, модератор форума

progressive
Новоприбывший
Сообщения: 3
Зарегистрирован: 11 апр 2012, 15:13
Репутация: 0

Re: ошибки gdalwarp

Сообщение progressive » 11 апр 2012, 16:55

спасибо большое за такую оперативность.

Ответить

Вернуться в «GDAL/OGR»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость