Форумы GIS-Lab.info

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Warning 6: gdalbuildvrt does not support ungeoreferenced image.

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

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