gdal2tiles плохое качество

[nd0ut]

Ребят, я новичёк в этом всём, сильно не пинайте

У меня есть данные с envisat asar в ихнем формате N1. Задача - отобразить их на карте.

Единственный более-менее рабочий способ, который я нашел - это использовать gdal2tiles. Но тайлы, полученные с его помощью получаются не качественные.

От чего так происходит? Можно ли пофиксить? Может есть другой способ решения моей задачи?

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

Единственный более-менее рабочий способ, который я нашел - это использовать gdal2tiles. Но тайлы, полученные с его помощью получаются не качественные.

Слово "некачественные" ни о чём не говорит, поясните, что конкретно вам нужно получить в итоге и чем вас не устраивает результат от gdal2tiles. И сформулируйте задачу полностью.

[nd0ut]

Под словом "некачественные" я подразумеваю пикселизацию изображения, а также заметную потерю детализации. И еще чёрные края.

В итоге надо получить результат как на скриншоте. Это в общем-то и есть задача полностью.

[nd0ut]

Вот человек столкнулся с такой же проблемой.

[Boris]

на первый взгляд имеют место следующие ошибки:
1. на тайлы делится исходное изображение, но в другой проекции. с учетом того, что производится большой поворот размер пикселя увеличивается на корень из двух, а полученные пиксели усредняют соседей неоптимальным образом;
2. где-то стоит автоподстройка уровней отображения градаций серого в пределах одного (сильно уменьшившегося) изображения. в результате контраст, присутствовавший в исходном изображении пропадает;
3. правая половина нижнего снимка имеет странный вид, т.к. в ней видно два изображения (смещение? наложение? полупрозрачность на карту?)

[Дмитрий Барышников]

Эх, начнем по порядку.
1) Почему вы уверены что для отображения на карте снимков с ENVISAT нужен именно gdal2tiles?
2) На какой карте какого софта вы хотите отобразить?
3) Версия утилит gdal?
4) Вывод gdalinfo входного файла.
5) Каким именно образом вы смотрите результат?

[ericsson]

4) Вывод gdalinfo входного файла.
5) Каким именно образом вы смотрите результат?

Я бы еще добавил "каким образом смотрите оригинал".
Там, наверняка, либо 16 бит, либо плавающая точка, так что разница в контрасте довольно предсказуемая, т.к. при нарезке на тайлы происходит неявное преобразование типов.

[nd0ut]

Эх, начнем по порядку.
1) Почему вы уверены что для отображения на карте снимков с ENVISAT нужен именно gdal2tiles?
Я в этом не уверен. Я не нашёл другого софта, с помощью которого это можно сделать.

2) На какой карте какого софта вы хотите отобразить?
Openlayers

3) Версия утилит gdal?
1.10.1

4) Вывод входного файла.
pastebin (с -mm ключём)

5) Каким именно образом вы смотрите результат?
В той картинке с двумя результатами:
слева - тифф, сделанный командой

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

gdal_translate ASA_WSM_1PNPDE20120120_091343_000002753110_00424_51731_1962.N1 -of GTiff -ot UInt16 -scale 0 3000 0 53000 translated.tif

справа - скриншот результата работы gdal2tiles

[nd0ut]

О, у меня получилось. Да, там было 16 бит, сконвертировал в 8, теперь всё хорошо.
Теперь еще один вопрос: на основе чего нужно выбирать параметры для scale? должны ли быть эти параметры общие для всех данных со спутника или индивидуальны для каждого конкретного снимка?

[nd0ut]

И ещё вопрос. Как заставить gdal2tiles генерить тайлы для всех тиффов (их много) ?
Когда я пытаюсь с помощью gdalbuildvrt создать виртуальный датасет из всех тиффов, чтобы потом скормить его в gdal2tiles, я получаю ошибку

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

Warning 6: gdalbuildvrt does not support ungeoreferenced image.

[Дмитрий Барышников]

Один вопрос одна тема.

[ Сообщение с мобильного устройства ]

[ericsson]

О, у меня получилось. Да, там было 16 бит, сконвертировал в 8, теперь всё хорошо.
Теперь еще один вопрос: на основе чего нужно выбирать параметры для scale? должны ли быть эти параметры общие для всех данных со спутника или индивидуальны для каждого конкретного снимка?

Ну, ответ простой: нужно найти эти величины для исходного изображения, для чего построить гистограмму. По ней увидите, где у вас там значащие величины, а где - пустота.
Мне для этого удобнее всего использовать ImageJ - там все визуально. В соседней теме про композиты Landsat я только что про это писал.

[nd0ut]

И ещё вопрос. Как заставить gdal2tiles генерить тайлы для всех тиффов (их много) ?
Когда я пытаюсь с помощью gdalbuildvrt создать виртуальный датасет из всех тиффов, чтобы потом скормить его в gdal2tiles, я получаю ошибку

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

Warning 6: gdalbuildvrt does not support ungeoreferenced image.

Оказалось, что перед этим надо сделать gdalwarp и всё ок.

всем спасибо за ответы