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

Вопросы по нескольким пакетам сразу, или вопросы, которые непонятно к какой ГИС отнести
Ответить
nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 27 апр 2017, 07:42

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

gamm
Гуру
Сообщения: 2536
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 517
Ваше звание: программист
Откуда: Казань

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

Сообщение gamm » 27 апр 2017, 15:55

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

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 27 апр 2017, 16:08

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

Игорь Белов
Гуру
Сообщения: 1412
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 868
Откуда: Казань

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

Сообщение Игорь Белов » 27 апр 2017, 16:18

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

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

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 27 апр 2017, 16:54

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

С другой стороны, совсем несложно написать программу самому. Только для больших растров лучше не Питона тянуть за хвост, а кодить на каком-нибудь Си или Фортране, IMHO. Всего три вопроса:
  1. формат растра;
  2. возможность растров с дырами в данных;
  3. конкретный состав выходных величин.
По составу:
  • "скорость" - модуль градиента, видимо;
  • "дивергенция скорости" и "Лапласиан" - одно и то же в точности;
  • "slope" - непонятно, каким боком здесь характеристика рельефа; к тому же она практически взаимозаменяема с модулем градиента, отличие в единицах измерения.
Игорь, спасибо. Да года четыре назад "закостЫлил" я всё это в basic'e'-Сёрфера по старому навыку - там и функционал вполне подходящий есть, да потом случайно всё за-erase-ил. :oops: .. растры небольшие (60x60)... ASCII-3х-колоночные... в любой др. формат перевести - не проблема... наверное-таки есть, конечно, реализация этого где-то, чтобы там с сетками-схемами различными конечно-разностными "("компасными", "шахматными"...) поиграться... под скоростью подразумевал визуализированный вектор в узлах регулярной сетки по двум составляющим градиента вдоль осей X и Y картографической проекции... Насчёт Лапласиана - понятно... Ещё и ротор, и кривизна есть, и относительная завиренность... Спасибо! Вспомню по шпаргалке-заданию 4-летней давности... Что-"нить" придумаю! Пока для начала актуально визуализировать вектор по двум составляющим градиента...

Аватара пользователя
paleogis
Модератор
Сообщения: 1112
Зарегистрирован: 22 мар 2009, 08:54
Статьи: 2
Репутация: 200
Ваше звание: Модератоо

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

Сообщение paleogis » 28 апр 2017, 06:41

https://youtu.be/F_phSdrCOA4

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

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 28 апр 2017, 07:11

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

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

Игорь Белов
Гуру
Сообщения: 1412
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 868
Откуда: Казань

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

Сообщение Игорь Белов » 29 апр 2017, 19:05

Чтобы не быть голословным, написал программу-макет на Си, которая выполняет некоторые операции векторного анализа:
  • 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
Вложения
nabla.7z
(132.75 КБ) 23 скачивания

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 29 апр 2017, 19:49

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

gamm
Гуру
Сообщения: 2536
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 517
Ваше звание: программист
Откуда: Казань

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

Сообщение gamm » 30 апр 2017, 07:51

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

Игорь Белов
Гуру
Сообщения: 1412
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 868
Откуда: Казань

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

Сообщение Игорь Белов » 30 апр 2017, 09:23

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

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 30 апр 2017, 09:33

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

nickleb
Гуру
Сообщения: 733
Зарегистрирован: 22 май 2010, 20:20
Репутация: 65

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

Сообщение nickleb » 01 май 2017, 17:20

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

Александр Мурый
Гуру
Сообщения: 5015
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 707
Ваше звание: званий не имею
Откуда: Москва

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

Сообщение Александр Мурый » 01 май 2017, 18:21

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

Ответить

Вернуться в «Общий - ПО»