Новогоднее обновление scanaerial

BushmanK · 2015-01-07 18:41:35

Чтобы облака на снимке не смущали, по маске, взятой из BQA, соответствующие пиксели вообще стоит обнулить до значения nodata. В ГИС это, правда, имеет больше смысла, так как при сборке мозаики на этом месте может оказаться читаемое изображение из другого снимка.

Last edited by BushmanK (2015-01-07 18:42:02)

Vort · 2015-01-08 05:28:31

BushmanK wrote:

Чтобы облака на снимке не смущали, по маске, взятой из BQA, соответствующие пиксели вообще стоит обнулить до значения nodata.

от того что будет тучка с дыркой посредине - легче не станет
да и тень от тучи никуда не исчезнет
биты в QA под неё есть, но алгоритм их пока что не заполняет
плюс к этому, алгоритм QA не шибко точный - много ложных срабатываний
если б он был хорош, я бы из него напрямую воду повытаскивал

BushmanK wrote:

В ГИС это, правда, имеет больше смысла, так как при сборке мозаики на этом месте может оказаться читаемое изображение из другого снимка.

если размыть маску да наложить данные с другого снимка, то да - в этом может быть смысл

BushmanK · 2015-01-08 05:51:46

Если не превращать плагин еще и в плагин к GRASS, то с этим, конечно, особо ничего не сделаешь. А так - есть и всякие продвинутые алгоритмы детекции облаков, и с тенями они работают (с последними еще неплохо справляются алгоритмы коррекции локального контраста). Но это так, технологическое отступление.

Vort · 2015-01-09 08:36:02

Добавил многопоточности скрипту обработки снимков:

@echo off

if not exist gdal2tiles.py exit

set min_zoom=6
set max_zoom=11
set scale_r_min=6000
set scale_r_max=7600
set scale_g_min=6200
set scale_g_max=7600
set scale_b_min=7000
set scale_b_max=8200
set scale_ir_min=5000
set scale_ir_max=15000
set src_dir=e:\Downloads\ls\
set src_name=LC81070152014263LGN00
set src_path=%src_dir%%src_name%

if "%1" == "" (
	del *.tif >nul 2>&1
	rmdir /s /q %src_name%_B5_merc_8bit >nul 2>&1
	rmdir /s /q %src_name%_RGB_merc_8bit >nul 2>&1
	rmdir /s /q tiles >nul 2>&1
	mkdir tiles >nul 2>&1

	call GDALShell.bat
	start cmd /c go.bat col
	call go.bat ir
)
if "%1" == "ir" (
	title Generating infrared tiles
	gdalwarp.exe -t_srs EPSG:3857 -r bilinear %src_path%_B5.tif %src_name%_B5_merc.tif
	gdal_translate.exe -scale %scale_ir_min% %scale_ir_max% -ot Byte %src_name%_B5_merc.tif %src_name%_B5_merc_8bit.tif
	gdal2tiles.py -z %min_zoom%-%max_zoom% %src_name%_B5_merc_8bit.tif
	move %src_name%_B5_merc_8bit tiles/ir >nul 2>&1
)
if "%1" == "col" (
	title Generating color tiles
	gdal_merge.py -o %src_name%_RGB.tif -co PHOTOMETRIC=RGB -separate %src_path%_B4.tif %src_path%_B3.tif %src_path%_B2.tif
	gdalwarp.exe -t_srs EPSG:3857 -co "PHOTOMETRIC=RGB" -r bilinear %src_name%_RGB.tif %src_name%_RGB_merc.tif
	gdal_translate.exe -scale %scale_r_min% %scale_r_max% -scale %scale_g_min% %scale_g_max% -scale %scale_b_min% %scale_b_max% -ot Byte %src_name%_RGB_merc.tif %src_name%_RGB_merc_8bit.tif
	gdal2tiles.py -z %min_zoom%-%max_zoom% %src_name%_RGB_merc_8bit.tif
	move %src_name%_RGB_merc_8bit tiles/col >nul 2>&1
)

Vitalts · 2015-01-09 10:00:13

А можно еще одно дополнение реализовать? watermark BBOX, при указании коротого для каждого тайла запрашивался бы тайл со сдвигом и с него бы копировался указанный прямоугольник, тем самым очищая первоначальный тайл.

Vort · 2015-01-09 15:55:03

Vitalts wrote:

А можно еще одно дополнение реализовать? watermark BBOX, при указании коротого для каждого тайла запрашивался бы тайл со сдвигом и с него бы копировался указанный прямоугольник, тем самым очищая первоначальный тайл.

можно поподробнее суть и как эта возможность должна применяться?

Vitalts · 2015-01-09 16:06:09

Имеется WMS, на тайлах которого располагается информация о собственнике, она мешает нормальной работе сканэриал. Если взять тайл с некоторым сдвигом, получим тайл, на котором опять таки будет копирайт, однако, то место, где он был на предыдущем снимке, на текущем будет нормальное изображение. В результате сканэриал будет пыхтеть над чистыми изображениями, и как следствие, результат обработки будет лучше.
Пример:
http://forum.openstreetmap.org/viewtopi … 24#p461024

Vort · 2015-01-10 05:26:22

Это слишком частная задача, чтоб её вносить в общую ветку scanaerial.
Тут скорее нужно писать некое подобие прокси-сервера, специально для этого случая.
В котором или действительно запрашивать WMS со сдвигом, или использовать тайлы без водяных знаков:
http://xgis.maaamet.ee/xGIS/XGis/tile/T … /3/3/3/390

Last edited by Vort (2015-01-10 05:27:47)

Vort · 2015-01-10 07:34:39

Вопрос к сообществу:
Существуют ли TMS сервера с проекцией, отличной от EPSG:3857?
Существуют ли TMS сервера с размерами тайлов, отличающимися от 256x256?
Если существуют, есть у кого-то примеры?
Если существуют, поддерживает ли такое JOSM?

Last edited by Vort (2015-01-10 07:35:27)

chnav · 2015-01-10 11:02:04

Vort wrote:

Существуют ли TMS сервера с проекцией, отличной от EPSG:3857?

Из известных - Яндекс EPSG:3395

Last edited by chnav (2015-01-10 11:05:05)

freeExec · 2015-01-10 11:27:42

Следовательно эту проекцию запретить

OverQuantum · 2015-01-10 12:05:47

Я знаю WMS, который работает только в EPSG:3301

wms:http://kaart.maaamet.ee/wms/kaart?FORMAT=image/png&TRANSPARENT=TRUE&VERSION=1.1.1&SERVICE=WMS&REQUEST=GetMap&LAYERS=MA-KAART&STYLES=&SRS={proj}&WIDTH={width}&HEIGHT={height}&BBOX={bbox}

Vort · 2015-01-10 14:18:15

chnav wrote:

Из известных - Яндекс EPSG:3395

спасибо
проверил - работает как надо
забавно что JOSM этой проекции не знает

OverQuantum wrote:

Я знаю WMS, который работает только в EPSG:3301

вопрос поддержки проекций в TMS меня как раз из-за него и заинтересовал
у них есть какой-то сервер, похоже что TMS. суть в том что он без водяных знаков
если знать что его тайлы (к примеру, по ссылке выше - 3/3/3/xxx) подхватятся scanaerial-ом, можно будет разобраться что значит 4й параметр
---
похоже, там какой-то не такой зум
так что врядли удастся напрямик подключить

Last edited by Vort (2015-01-10 14:55:49)

Vort · 2015-01-10 16:21:04

Добавлена опция {switch:} для TMS серверов:
https://github.com/Vort/scanaerial/arch … switch.zip

Vort · 2015-01-11 17:33:59

Последние изменения уже в основной ветке.
Только есть один ньюанс: мне пришлось самостоятельно править английский текст.
В связи с этим, прошу сообщество помочь найти и исправить ошибки в моих изменениях файла README.org и страницы в вики.

Komяpa · 2015-01-15 09:57:31

На всякий случай: scanaerial, когда писался, предполагалось, что всегда будет цепляться на локально установленный twms с кешем.

Все операции фильтрации-обработки-допривязки-подавления пустых тайлов делались средствами twms :3

Hind · 2015-01-15 10:28:15

Кстати, TWMS умеет раздавать геотиффы в виде тайлов без предварительной нарезки?

BushmanK · 2015-01-15 23:28:02

Судя по тому, что написано тут https://code.google.com/p/twms/wiki/Config - только тайлы.

ad47 · 2015-02-14 12:42:54

Попробовал обновлённую версию - понравилось. Даже в Московской области можно найти крупные водоёмы, обведённые по старым снимкам Landsat, которые неплохо уточняются по Landsat8. Иногда неплохо получается векторизовать что-то не очень крупное и по снимкам высокого разрешения, но после этого требуется ручная проверка и упрощение контуров (туда могут попасть тени от деревьев и окрестные леса, а диагональные линии векторизуются ступеньками в несколько пикселов) и иногда проще сразу обвести вручную.

Опыт показал, что для широты 55° для снимков Landsat8 (яркость из канала 8 + цвета из 7-5-3) z12 не достаточно. Вот пример достаточно хорошего совпадения котура, полученного из z13, со снимком MapBox, в то время как контур из z12 не совпадает даже с исходным Landsat (по размеру озеро вписывается в квадрат 1х1 км). https://www.dropbox.com/s/0v5r07qegztrk … 2.png?dl=0

Кстати, локальный сервер поднимать не пришлось, хватило строки такого вида:

server_url = file:/D:\Tiles\LCxxx\{zoom}\{x}\{-y}.png

Announcement

#101 2015-01-07 18:41:35

Re: Новогоднее обновление scanaerial

#102 2015-01-08 05:28:31

Re: Новогоднее обновление scanaerial

#103 2015-01-08 05:51:46

Re: Новогоднее обновление scanaerial

#104 2015-01-09 08:36:02

Re: Новогоднее обновление scanaerial

#105 2015-01-09 10:00:13

Re: Новогоднее обновление scanaerial

#106 2015-01-09 15:55:03

Re: Новогоднее обновление scanaerial

#107 2015-01-09 16:06:09

Re: Новогоднее обновление scanaerial

#108 2015-01-10 05:26:22

Re: Новогоднее обновление scanaerial

#109 2015-01-10 07:34:39

Re: Новогоднее обновление scanaerial

#110 2015-01-10 11:02:04

Re: Новогоднее обновление scanaerial

#111 2015-01-10 11:27:42

Re: Новогоднее обновление scanaerial

#112 2015-01-10 12:05:47

Re: Новогоднее обновление scanaerial

#113 2015-01-10 14:18:15

Re: Новогоднее обновление scanaerial

#114 2015-01-10 16:21:04

Re: Новогоднее обновление scanaerial

#115 2015-01-11 17:33:59

Re: Новогоднее обновление scanaerial

#116 2015-01-15 09:57:31

Re: Новогоднее обновление scanaerial

#117 2015-01-15 10:28:15

Re: Новогоднее обновление scanaerial

#118 2015-01-15 23:28:02

Re: Новогоднее обновление scanaerial

#119 2015-02-14 12:42:54

Re: Новогоднее обновление scanaerial

Board footer