Size: a a a

спбгеотех

2020 July 31

EK

Eduard Kazakov in спбгеотех
White Whale
Всем привет! Я слышала тут сидят мало того что умные люди, так ещё есть и картографы (Эдик привет из нансена :))
Надеюсь мне смогут ответить на интересующий меня вопрос и что-нибудь подсказать. У меня есть пачка спутниковых данных в проекции espg:4326 в виде numpy 2d массивов и у меня возникла проблема с расчетом реальных площадей  (в км²). Есть ли какие-нибудь непринуждённые и изящные пути посчитать интересующие меня площади некоторых областей без репроецирования всех данных с помощью питона?
Привет! В epsg:4326 фактическая площадь пикселей претерпевает заметные изменения от строчки к строчке (при перемещении в меридиональном направлении), и просто и изящно вряд ли получится, без некоторой цепочки процессинга. У тебя же там объекты - сложные контуры по матрицам-маскам, да?
Данные в матрицах хранятся какие, кстати, по всему земному шару, -180,180,-90,90 ?
источник

WW

White Whale in спбгеотех
Eduard Kazakov
Привет! В epsg:4326 фактическая площадь пикселей претерпевает заметные изменения от строчки к строчке (при перемещении в меридиональном направлении), и просто и изящно вряд ли получится, без некоторой цепочки процессинга. У тебя же там объекты - сложные контуры по матрицам-маскам, да?
Данные в матрицах хранятся какие, кстати, по всему земному шару, -180,180,-90,90 ?
Да, по всему земному шару. Маски сложные :)
Ну вроде сейчас пересчитали исходя из заданных координат (из исходных данных) к каждому пикселю в маске. Площади стали чуть больше похожи на реальность. Но насколько это правильно (просто пересчитывая площадь каждого пикселя принимая что он все ещё квадратный на всех широтах) и хотелось бы знать существуют ли функции попроще, которые учитывают особенности проекции
источник

f

fr1 in спбгеотех
White Whale
Завезли но я совсем не картограф к сожалению. Даже близко. Ничего не поняла в нем :/
трансформировать координаты из wgs84 в соответствующую зону UTM (если у вас там не размеры регионов) и считать по старинке в обычной декартовой системе
источник

EK

Eduard Kazakov in спбгеотех
Просто могло бы быть, если бы территория была локальной, одно море например: тогда с приемлемой точностью можно было одной командой перепроецировать в локальную проекцию, и брать декартову площадь в ней, умножив количество пикселей попавших в маску на площадь одного пикселя.

Из того, что будет работать, но придётся потратить некоторое время и сделать себе код:
1. numpy 2d массив -> GDAL датасет (не проблема, т.к. пространственные параметры матрицы известны)
2. gdal_polygonize на результат п.1, получается векторный слой, где объекты под маской имеют значение атрибута, равное значениям пикселей в этой маске.
3. через python-ogr открыть вектор из п.2, удалить все объекты которые по значению атрибута не попали в маску
4. В цикле пробегать для всех исходных nympy-матриц пункты 1-3, итоговые векторные объекты закидывая в новый (единый) векторный слой, помечая в атрибутах принадлежность объекта матрице (дата, или идентфикатор какой).
5. Открыть получившийся в п.4 векторный слой в QGIS, настроить вычисление площадей на эллипсоиде, рассчитать для всех объектов площадь через калькулятор полей ($area).

Может кто подскажет, умеет ли pyproj или другая библиотека в сфероидическую геодезию над объектами-полигонами, или геометриями доставаемыми OGR-ом. Если есть что-то - то 4,5 не нужны, сразу вектор из 3 в соответствующую функцию и вот она площадь.
источник

AR

Alex Rybakov in спбгеотех
О..не знал, спасибо 👍
источник
2020 August 02

A

Asiia in спбгеотех
Alex Rybakov
ESRI по подписке дает доступ к данным AIRBUS....
Добрый день, а что значит по подписке? Если есть наградная лицензия на pro - даёт?
источник

AR

Alex Rybakov in спбгеотех
Asiia
Добрый день, а что значит по подписке? Если есть наградная лицензия на pro - даёт?
Не в курсе про наградные... Надо пробовать. Описание в блоге Esri-Cis ищите
источник

A

Asiia in спбгеотех
Ну в смысле тем у кого лицензия. Спасибо
источник
2020 August 03

LK

Leonid Kolesnichenko in спбгеотех
Здравствуйте! Подскажите, как правильно добавить МСК-координаты в Geoserver, чтобы не было конфликтов с другими уже существующими СК?
источник

d

dmiteco [E16] in спбгеотех
Leonid Kolesnichenko
Здравствуйте! Подскажите, как правильно добавить МСК-координаты в Geoserver, чтобы не было конфликтов с другими уже существующими СК?
Пропиши в файл конфигурации ск, со своим epsg, которого нет в списке общепринятых.
источник

АТ

Алексей Трошин... in спбгеотех
А в EPSG нет какого-либо диапазона, зарезервированного под такие случаи?
источник

n

nikikiki in спбгеотех
По негласному соглашению вроде задают  srid > 100 000 у пользовательских СК.
источник

LK

Leonid Kolesnichenko in спбгеотех
nikikiki
По негласному соглашению вроде задают  srid > 100 000 у пользовательских СК.
интересно почитать про негласные соглашения
источник

AS

Alexander Semenov in спбгеотех
Кто в курсе, где границы финского залива скачать можно? Датасет с геофабрика на северо-запад че-то не очень.
источник

AS

Alexander Semenov in спбгеотех
Ответ “купи у некстгиса” не подходит
источник

DK

Dmitry Kolesov in спбгеотех
Нарисовать по снимку самому?
источник

AS

Alexander Semenov in спбгеотех
Что-нибудь готовое хочется
источник

IL

Ilya Levashev in спбгеотех
Тут очевидно косяк в подложке)
источник

IL

Ilya Levashev in спбгеотех
или это вектор с геофабрика?
источник

AS

Alexander Semenov in спбгеотех
вектор с геофабрика
источник