Страница 1 из 1

Расчет фрактальной размерности Минковского Qgis

Добавлено: 08 дек 2018, 23:22
Блоппер
Доброго времени суток!
Мне необходимо выполнить расчет фрактальной размерности Минковского Истринского района Московской области. Я полный новичок и не знаю с чего начать. Прочитала кучу статей, но так и не смогла разобраться. Откуда мне вытащить гидрографическую сеть района? Как запустить этот калькулятор?
Буду благодарна за подробные указания, спасибо

Re: Расчет фрактальной размерности Минковского Qgis

Добавлено: 09 дек 2018, 07:17
gamm
если речь о покрытии объектов ящиками (https://en.wikipedia.org/wiki/Minkowski ... _dimension), то достаточно GDAL и чего-нибудь скриптового (R, Python даже bash). В цикле создаете растровые сетки, растеризуете на них вектор gdal_rasterize с параметром -at -a ID -tr xres yres, -at дает нужное покрытие, ID - поле с ID объектов, -tr res_i res_i - очередной размеры ячеек

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

dal_rasterize -b 1 at -a ID -init 0 -l mask -ot Int32 -tr res_i res_i mask.shp work_i.tif
потом подсчитываете число пикселей с каждым ID, лучше в R, Python, например векторизовать (gdal_polygonize.py), и использовать SQL запрос в GDAL (или просто грузить растры в R, Python, и там считать одной командой).

Накопив статистику по разным разрешениям, посчитать размерность. Желательно сделать несколько вариантов в каждом разрешении со сдвигом границ сетки в пределах res_i (можно подвигать созданные сетки с использованием gdal_translate, а потом растеризовать на них, -init 0 их почистит)

как-то так ...

Re: Расчет фрактальной размерности Минковского Qgis

Добавлено: 09 дек 2018, 15:26
Эдуард Казаков
Для QGIS 2.* есть готовый инструмент:
http://gis-lab.info/qa/minkowski-dimension-qgis.html

Re: Расчет фрактальной размерности Минковского Qgis

Добавлено: 09 дек 2018, 16:52
jerry-maori
1. Гидрографию можно взять в виде вектора из той же OSM и растеризовать.
2. Не забывать о самом вкусном -- о твёрдой уверенности в том, что то наше облако точек будет аппроксимироваться именно степенью в натуральных величинах ( ну или линейно в билограрифмах)