Форумы GIS-Lab.info

Здравствуйте. Есть ли free-ГИС (плагин для free-ГИС) , которая(-ый) содержала (-л) бы наиболее мощный и полный блок расчета полей гидродинамических параметров некой данной характеристики по прямоугольной сеточной области? К примеру: есть растр некой характеристики по горизонтальной XY-регулярной сетке в картографической проекции. Надо посчитать градиенты вдоль осей OX и OY, скорость, дивергенцию скорости, завихренность, лапласиан, slope поля этой некой характеристики и визуализировать поля этих гидродинамических параметров. М.б., есть таковой инструмент в R или Pyhon?

nickleb писал(а):М.б., есть таковой инструмент в R или Pyhon?

тут у вас намешано разного. Сами эти характеристики считаются конечными разностями руками, без проблем, лучше в Питоне, если размер большой, если 100х100, то все равно. А потом полученный растр визуализируйте, сколько хотите. Готовое в R - это навряд ли, в Питоне - более вероятно, но тоже сомнительно, в Матлабе - почти наверняка есть что-то. Но проще в Питоне слепить из numpy+scipy, там вся линейная алгебра есть, например дивергенция тут, правда код не проверял.

gamm писал(а):... проще в Питоне слепить из numpy+scipy, там вся линейная алгебра есть, например дивергенция тут, правда код не проверял.,,,

да, gamm, наверное... потом в matplotlib и нарисовать... спасибо!

С одной стороны, такие базовые вещи наверняка где-то реализованы.

С другой стороны, совсем несложно написать программу самому. Только для больших растров лучше не Питона тянуть за хвост, а кодить на каком-нибудь Си или Фортране, IMHO. Всего три вопроса:

формат растра;
возможность растров с дырами в данных;
конкретный состав выходных величин.

По составу:

"скорость" - модуль градиента, видимо;
"дивергенция скорости" и "Лапласиан" - одно и то же в точности;
"slope" - непонятно, каким боком здесь характеристика рельефа; к тому же она практически взаимозаменяема с модулем градиента, отличие в единицах измерения.

Игорь Белов писал(а):С одной стороны, такие базовые вещи наверняка где-то реализованы.

С другой стороны, совсем несложно написать программу самому. Только для больших растров лучше не Питона тянуть за хвост, а кодить на каком-нибудь Си или Фортране, IMHO. Всего три вопроса:
формат растра;
возможность растров с дырами в данных;
конкретный состав выходных величин.
По составу:
"скорость" - модуль градиента, видимо;
"дивергенция скорости" и "Лапласиан" - одно и то же в точности;
"slope" - непонятно, каким боком здесь характеристика рельефа; к тому же она практически взаимозаменяема с модулем градиента, отличие в единицах измерения.

Игорь, спасибо. Да года четыре назад "закостЫлил" я всё это в basic'e'-Сёрфера по старому навыку - там и функционал вполне подходящий есть, да потом случайно всё за-erase-ил.

.. растры небольшие (60x60)... ASCII-3х-колоночные... в любой др. формат перевести - не проблема... наверное-таки есть, конечно, реализация этого где-то, чтобы там с сетками-схемами различными конечно-разностными "("компасными", "шахматными"...) поиграться... под скоростью подразумевал визуализированный вектор в узлах регулярной сетки по двум составляющим градиента вдоль осей X и Y картографической проекции... Насчёт Лапласиана - понятно... Ещё и ротор, и кривизна есть, и относительная завиренность... Спасибо! Вспомню по шпаргалке-заданию 4-летней давности... Что-"нить" придумаю! Пока для начала актуально визуализировать вектор по двум составляющим градиента...

https://youtu.be/F_phSdrCOA4

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

paleogis писал(а):https://youtu.be/F_phSdrCOA4

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

paleogis, спасибо! Ай, да Дарья! - Как много умеет и успевает! И видел же этот ролик я... - Но, наверное, больше просто Дарьей любовался

И в сторону SAGA также думал...

Чтобы не быть голословным, написал программу-макет на Си, которая выполняет некоторые операции векторного анализа:

g_x - градиент по оси X;
g_y - градиент по оси Y;
grad - модуль градиента;
div - дивергенция градиента.

Ротор градиента не считается, ибо это ноль по определению. Вычисление других операций можно добавлять.

Выбрал формат растров BIL/HDR:

BIL - бинарная матрица без заголовка;
HDR - заголовок в человекочитаемом виде;
PRJ - крайне редкий формат растра, для которого GDAL признаёт PRJ.

Программа не анализирует заголовок глубоко. Она полагается на то, что данные в формате действительных 32-битных чисел под процессор Intel. Всего она читает пять величин: размеры матрицы NROWS и NCOLS, разрешения по осям XDIM и YDIM для вычисления производных и величину NODATA, если она есть.

Программа запускается в командной строке с названием растра (без расширения) в качестве аргумента:

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

nabla Helens2

Программа не копирует файлы HDR и PRJ. Это можно делать вручную или написав командный файл по типу прилагаемого скрипта:

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

. cp_hdr.sh Helens2

Игорь Белов писал(а):Чтобы не быть голословным, написал программу-макет на Си, которая выполняет некоторые операции векторного анализа:...

Игорь, спасибо огромное за Ваше драгоценное время и приложенные Вами усилия, и ,конечно, за саму программку! "Погоняю" на следующей неделе обязательно!

есть такая штуковина, бесплатная, Generic Mapping Tools (GMT), она заточена под геофизику и вообще разную около-геологию, насколько я понял. По крайней мере, наши сейсмики и геофизики были довольны, она создает шикарные карты - они настроили все себе буквально за неделю. Векторные поля она тоже строит, выхлоп в PostScript, который они вставляли в PDF.

gamm писал(а):есть такая штуковина, бесплатная, Generic Mapping Tools (GMT), она заточена под геофизику и вообще разную около-геологию, насколько я понял. По крайней мере, наши сейсмики и геофизики были довольны, она создает шикарные карты - они настроили все себе буквально за неделю. Векторные поля она тоже строит, выхлоп в PostScript, который они вставляли в PDF.

Сам пользуюсь и другим рекомендую.

gamm писал(а):есть такая штуковина, бесплатная, Generic Mapping Tools (GMT), она заточена под геофизику и вообще разную около-геологию, насколько я понял. По крайней мере, наши сейсмики и геофизики были довольны, она создает шикарные карты - они настроили все себе буквально за неделю. Векторные поля она тоже строит, выхлоп в PostScript, который они вставляли в PDF.

если жизнь заставит в OPTIONS его SYNOPSISа разобраться, то - конечно, да... он и с CDF давно "дружит"... и развивается...

... вот интересный обзор свободных программных библиотек:
s.arboreus.com/2008/03/free-software-for-solving-pdes.html

В составе GRASS GIS есть несколько модулей по теме + вот ещё.

Форумы GIS-Lab.info

Free-ГИС с блоком расчета гидродинамических параметров

Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров

Re: Free-ГИС с блоком расчета гидродинамических параметров