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

Vitalts · 2015-01-06 16:46:25

В общем, я уже взялся за это дело, там работы на 10 минут с тестированием, однако, не успел, отвлекли, а теперь вот, домой пора.
ЗЫ: ну и холодрыга на улице, пальцы замерзли сообщение набирать.

Kostik · 2015-01-06 16:57:11

Vort wrote:

Архив последней версии: https://github.com/Vort/scanaerial/arch … er_fix.zip

Что то новенькое, похоже на сыр. Причём результат повторяем.

Vort · 2015-01-06 16:59:08

Kostik wrote:

Что то новенькое, похоже на сыр. Причём результат повторяем.
http://i066.radikal.ru/1501/84/720f87e834a7t.jpg

фотики перечёркнутые чтоли?

Last edited by Vort (2015-01-06 17:02:25)

Kostik · 2015-01-06 17:02:03

Vort wrote:

фотики перечёркнутые чтоли?

Не понял вопроса, какие фотики?
Скачал и обновился, все конфиги дефолтные, ничего не менял. Динамический зум. Это Bing.
Случайное место в Мурманской области, вот OSM файл.

Last edited by Kostik (2015-01-06 17:04:25)

Vort · 2015-01-06 17:02:33

Kostik wrote:

Скачал и обновился, все конфиги дефолтные, ничего не менял. Динамический зум.

это Bing? можно координаты?

Vitalts · 2015-01-06 17:04:03

Перечеркнутые тайлы, по которым ведётся сканирование.

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

По умолчанию же, вроде, не Бинг в архиве?

Kostik · 2015-01-06 17:06:44

Vort wrote:

Kostik wrote:
Скачал и обновился, все конфиги дефолтные, ничего не менял. Динамический зум.
это Bing? можно координаты?

http://www.openstreetmap.org/search?que … 5/35.57111

Vitalts wrote:

не Бинг в архиве?

Хреновый Bing не спорю.

server_api = bing
server_name = Bing
server_url = http://dev.virtualearth.net/REST/v1/Imagery/Metadata/Aerial?include=ImageryProviders&output=xml&key=Arzdiw4nlOJzRwOz__qailc8NiR31Tt51dN2D7cm57NrnceZnCpgOkmJhNpGoppU
projection = EPSG:3857
tile_sizex = 256
tile_sizey = 256

[SCAN]
douglas_peucker_epsilon = 0.60
deactivate_simplifying = 0
colour_str = 30
maxfilter_setting = 3
size_limit = 400

[TAGS]
natural = water

Last edited by Kostik (2015-01-06 17:11:46)

Vort · 2015-01-06 17:12:32

Вобщем, это обведенный пустой тайл.
Решение - уменьшить зум или использовать фиксированный.

Детект пустых тайлов, конечно, не помешал бы.
Но надо подумать как это сделать.

Сразу скажу по тому озеру:
тут нужен colour_str = 60

Last edited by Vort (2015-01-06 17:13:33)

Kostik · 2015-01-06 17:14:48

Vort wrote:

Сразу скажу по тому озеру:
тут нужен colour_str = 60

Конкретно это проще обвести вручную.

Vort · 2015-01-06 17:16:38

Vitalts wrote:

В общем, я уже взялся за это дело, там работы на 10 минут с тестированием, однако, не успел, отвлекли, а теперь вот, домой пора.

тут надо продумать как быть с fixedzoomlevel

Vitalts · 2015-01-06 17:19:35

Фикседзум разве не для случаев параметров без з? Может максзум внедрить?
Сорри, цитату не прочитал. А ничего с ним не делать. Да и вообще, как бы это не показалось странным, я собираюсь передавать конфиг файл первым параметром, не смотря на его необязательность, там ему самое место. Можете ругать, весь во внимании

Last edited by Vitalts (2015-01-06 17:23:12)

Vort · 2015-01-06 17:21:15

Vitalts wrote:

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

а куда ж тогда вставлять опциональное имя конфига?
третья позиция не годится - там зум
четвёртая не годится - если зума нет, она будет третьей

----

Я таки угадал насчёт фотиков.
Вот они: http://ecn.t0.tiles.virtualearth.net/ti … peg?g=3179

Last edited by Vort (2015-01-06 17:21:50)

Vitalts · 2015-01-06 17:29:21

Можно же определять является ли параметр конфигом, например по наличию расширения .cfg, и в зависимости от этого делать оффсет положения правых параметров.

Vort · 2015-01-06 17:38:05

Vitalts wrote:

Можно же определять является ли параметр конфигом, например по наличию расширения .cfg, и в зависимости от этого делать оффсет положения правых параметров.

или убрать опциональность третьего параметра и для fixed_zoom`а передавать туда 0

Vitalts · 2015-01-06 17:41:04

Но тогда пользователям сканэриал без зума придётся перенастраивать свои параметры запуска, не кошерно.

Vort · 2015-01-06 17:59:10

С "сыром" покончено:
https://github.com/Vort/scanaerial/arch … _tiles.zip
(если я нигде не налажал)

Vitalts wrote:

Но тогда пользователям сканэриал без зума придётся перенастраивать свои параметры запуска, не кошерно.

делать параметр-мутант, который то ли зум, то ли конфиг - тоже не супер
<s>реализовуй, наверно, в своей локальной ветке, пока что-нибудь годное не придумается</s>
готово. см ниже.

Last edited by Vort (2015-01-06 18:30:30)

Vitalts · 2015-01-06 18:13:48

Вообще-то никто и никогда не ограничивал кол-во необязятельных параметров, важно лишь соблюдать оговоренный порядок следования, дабы можно было определён что есть что. Это норма, общепринято лишь, что бы необязательные параметры следовали после обязательных. В следствии чего, для общественности готов перенести конфиг и в конец параметров, но мне этот вариант не по душе, ну не место ему там.
Наиболее однозначный вариант, передавать параметры парами ключ-значение, тут на порядок вообще наплевать. Однако, и это требует изменений параметров запуска настроенных скриптов, чего хотелось бы избежать.

Vort · 2015-01-06 18:15:45

Кажется, знаю что можно сделать.

Фиксированный зум:
scanaerial.py {lat} {lon}
scanaerial.py {lat} {lon} 0
scanaerial.py {lat} {lon} 0 config

Динамический зум:
scanaerial.py {lat} {lon} {TZoom}
scanaerial.py {lat} {lon} {TZoom} config

Сейчас попробую накодить

Last edited by Vort (2015-01-06 18:19:47)

Vort · 2015-01-06 18:29:03

Имя конфига 4м параметром готово:
https://github.com/Vort/scanaerial/arch … g_name.zip

---
Ох... кто ж это всё будет в вики документировать )

Last edited by Vort (2015-01-06 18:32:58)

BushmanK · 2015-01-07 08:25:01

К вопросу о том, работать ли с композитом или с одним каналом.
Вспомним простую векторную геометрию: модуль вектора (его длина) равна квадратному корню из суммы квадратов разностей соответствующих координат концов вектора.
Имея один восьмибитный канал с подправленным контрастом, мы получаем типичное значение расстояния между значениями, соответствующими воде и суше, равное некоему числу A=sqrt(dX^2)
Имея в композите два канала - B=sqrt(dX^2+dY^2)
Если предположить, что порядок dX и dY - один и тот же, получим B=sqrt(2*dX^2). То есть расстояние будет больше в sqrt(2) рааза, то есть где-то 1,4 раза. Для трех каналов с приблизительно одинаковым контрастом получим уже sqrt(3), то есть около 1,7 раза. Даже если разница в каком-то канале будет совсем небольшой, под корнем стоит сумма, так что даже это дает небольшое увеличение различимости двух любых классов, которые хоть как-то различаются на использованных каналах.

Так что композит выгоду дает без всяких сомнений.

Vort · 2015-01-07 10:07:50

Какая же ужасная вещь эти GDAL скрипты.
Неинтуитивные, корявые, глюченые.

Но несмотря на это, подкрасить тучки мне всё-таки удалось

@echo off

set min_zoom=6
set max_zoom=11
set src_level_min=4000
set src_level_max=14000
set src_dir=e:\Downloads\ls\
set src_name=LC81100132014252LGN00
set src_path=%src_dir%%src_name%

del %src_name%_clouds.tif >nul 2>&1
del %src_name%_clouds_merc.tif >nul 2>&1
del %src_name%_clouds_merc_8bit.tif >nul 2>&1
del %src_name%_B5_merc.tif >nul 2>&1
del %src_name%_B5_merc_8bit.tif >nul 2>&1
del %src_name%_B5_merc_clouds_gb_8bit.tif >nul 2>&1
del %src_name%_B5_merc_clouds_rgb_8bit.tif >nul 2>&1

call GDALShell.bat
gdal_calc.py -A %src_path%_BQA.TIF --outfile=%src_name%_clouds.tif --calc="(A&0xC000)*(A&0xC000>0x4000)"
gdalwarp.exe -r bilinear -t_srs EPSG:3857 %src_name%_clouds.tif %src_name%_clouds_merc.tif
gdal_translate.exe -scale 0 60000 -ot Byte %src_name%_clouds_merc.tif %src_name%_clouds_merc_8bit.tif
gdalwarp.exe -r bilinear -t_srs EPSG:3857 %src_path%_B5.TIF %src_name%_B5_merc.tif
gdal_translate.exe -scale %src_level_min% %src_level_max% -ot Byte %src_name%_B5_merc.tif %src_name%_B5_merc_8bit.tif
gdal_calc.py -A %src_name%_B5_merc_8bit.tif -B %src_name%_clouds_merc_8bit.tif --outfile=%src_name%_B5_merc_clouds_gb_8bit.tif --calc="A ^ B"
gdal_merge.py -co PHOTOMETRIC=RGB -separate %src_name%_B5_merc_8bit.tif %src_name%_B5_merc_clouds_gb_8bit.tif %src_name%_B5_merc_clouds_gb_8bit.tif
ren out.tif %src_name%_B5_merc_clouds_rgb_8bit.tif
gdal2tiles.py -z %min_zoom%-%max_zoom% %src_name%_B5_merc_clouds_rgb_8bit.tif

Last edited by Vort (2015-01-09 08:36:34)

Vort · 2015-01-07 11:53:38

Хоть результат и получился интересным, пользоваться им оказалось неудобно.

Буду думать дальше.

Vort · 2015-01-07 14:56:17

Красные тучи мне сильно резали глаза.
Поэтому, сделал два набора тайлов - из цветного снимка и из инфракрасного.
Если что-то не видно на инфракрасном снимке - включаю слой с цветным и рассматриваю его.

Вот текущий скрипт:

@echo off

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=LC81100132014252LGN00
set src_path=%src_dir%%src_name%

if not exist gdal2tiles.py exit

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

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

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

Last edited by Vort (2015-01-09 08:36:22)

Vort · 2015-01-07 16:48:11

к примеру: белое пятно на ИК снимке по центру озера может быть: островом, растениями или облаком.
на цветном снимке это пятно можно рассмотреть и узнать его тип по его цвету.

Last edited by Vort (2015-01-07 16:49:09)

Announcement

#76 2015-01-06 16:46:25

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

#77 2015-01-06 16:57:11

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

#78 2015-01-06 16:59:08

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

#79 2015-01-06 17:02:03

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

#80 2015-01-06 17:02:33

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

#81 2015-01-06 17:04:03

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

#82 2015-01-06 17:06:09

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

#83 2015-01-06 17:06:44

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

#84 2015-01-06 17:12:32

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

#85 2015-01-06 17:14:48

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

#86 2015-01-06 17:16:38

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

#87 2015-01-06 17:19:35

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

#88 2015-01-06 17:21:15

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

#89 2015-01-06 17:29:21

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

#90 2015-01-06 17:38:05

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

#91 2015-01-06 17:41:04

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

#92 2015-01-06 17:59:10

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

#93 2015-01-06 18:13:48

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

#94 2015-01-06 18:15:45

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

#95 2015-01-06 18:29:03

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

#96 2015-01-07 08:25:01

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

#97 2015-01-07 10:07:50

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

#98 2015-01-07 11:53:38

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

#99 2015-01-07 14:56:17

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

#100 2015-01-07 16:48:11

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

Board footer