Форумы GIS-Lab.info

Уважаемые эксперты!

В ходе одного исследования возникла необходимость использовать ГИС для визуализации данных. Выбрал бесплатную QGis. Но реализовать свои ожидания не смог. Прошу помощи в следующих вопросах:

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

2) Вырузка координат сетки
Необходимо выгрузить коордитаны центров ячеек построенной сетки с заданным шагом (п. а предыдущего вопроса). Выгрузить данные лучше в Excel.

3) Возможно ли определение расстояния от заданных точек, до, например сети дорог (кратчайшее). Например. Необходимо определить, на каком расстоянии от дороги находится заданная точка. При этом будет здорово использовать карты OpenStreetMaps.

Буду рад любой помощи (готового решения, ссылки на литературу или просто сочувствия )

Вопрос обширный Попробуйте почитать в разделе Статьи. Там Вам могут оказаться полезными статьи о Sextante, Построении сеток, Создания карты из текстового файла, Плстроение графа дорог.
Поиск работает.
На всякий случай: Если мои советы не помогут- сочувствие я тоже прилогаю: )

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

И еще. Обозначенная тема не информативна и обширна.

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

triangulator писал(а):а) построить сетку с определенным шагом

Вектор\Выборка\Регулярная сетка

Назначить уникальный ID каждой ячейке сетки (не помню делает ли это инструмент выше).

triangulator писал(а):б) пространственным запросом выбрать данные, попадающие в соты сетки

Вектор\Управление данными\Объединение атрибутов по районам - объединить точки с полигонами сетки.
(каждой точке присоединится код полигона сетки).

triangulator писал(а):в) вычислить среднее (лучше медианное) значение внутри каждой соты

Это вообще лучше делать не в QGIS, а скажем в R, банальная функция aggregate(). Дальше выгрузить в DBF/CSV в стиле: и присоединить эту таблицу обратно к сгенерированной сетке через Join

Максим Дубинин писал(а):

triangulator писал(а): Это вообще лучше делать не в QGIS, а скажем в R, банальная функция aggregate(). Дальше выгрузить в DBF/CSV в стиле:

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

CELLID MEDIANVAL
1 1000
2 1100

и присоединить эту таблицу обратно к сгенерированной сетке через Join

Правильно ли я понял насчет R. Эта программа позволит мне вычислить различные статистические параметры в разрезе заданных областей? Поясню. Имеется массив данных с координатами точки и определенным её значением. Необходимо сгруппировать данные значения по заданным областям, и затем посчитать статистику по каждой области (среднее, медиана, коэф. вариации и пр.).
Спсибо за ответы.

правильно, но вы можете вычислять где угодно еще, если вам так удобнее

Максим Дубинин писал(а): и присоединить эту таблицу обратно к сгенерированной сетке через Join

В QGIS, при присоединении вычисленных данных, они принимают тип "string" и соответственно классификацию уже не провести. В чем может быть проблема?

Нашел способ. В модуле SEXTANTE выбираем инструмент QGIS..-Vector general tools - Join Attributes table

все это можно самому легко закодить за неделю-две. Единственная проблема здесь - вывести полученную в результате привязанною картинку. Также, если вы используете, например, Bitmap то не надо забывать, что у него ограничение на 4 миллиарда по длине, поэтому надо будет делать несколько картинок, если не париться.

DevGis, сообщение совершенно не в тему.

У самого была подобная задача как у автора темы, учитывая совсем поверхностные знания в этих вопросах и короткие ответы в этой ветке, расскажу как это сделал - может кому пригодится, кто только узнал о QGIS и R.
I. QGIS
1. Загружаем шейпфайл топоосновы
2. Загружаем файл с данными отбора проб, формат .csv, (например, water.csv - названия полей латиницей)
3. Создаем регулярную сетку с заданным шагом (вектор/выборка/регулярная сетка: границы как у слоя water, создать сетку как полигоны) и сохраняем в шейп-файл, например grid.shp
4. Проставляем каждой станции ИД квадрата регулярной сетки (вектор/управление данными/объединение атрибутов по районам: целевой слой-станции, объединяем со слоем grid, свойства атрибутов - сохраняем атрибуты первого обнар. объекта. Сохраняем в water_grid.shp)
5. сохраняем water_rgid.shp в water_grid.csv (слой/сохранить как)
II.Блокнот и EXCEL
В Блокноте заменяем "," на ";", далее удаляем пробелы. Удостоверяемся в отсутствии не численных значений.
В EXCEL пустым ячейкам присваиваем значение "NA". Удаляем поля не численного формата, это важно, иначе функция aggregate в R будет выдавать ошибку (оставляем только ИД, зачения показателей и долготу с широтой - нужны только для того чтоб эту таблицу зугрузить в QGIS Сохраняем с разделителем ";"(MS-DOS)
III. R
выполняем функцию aggregate с помощью функции write.table сохраняем результаты вычислений например в aggdata.csv
IV. блокнот
открываем aggdata.csv в блокноте, удаляем все NA
V. QGIS - классификация данных
1. добавляем aggdata.csv
2. в SEXTANTE делаем связь таб. GRID и aggdata и сохраняем в новую таблицу и делаем по ней классификацию

по повросу 3, автора темы, можно здесь посмотреть.

losdimas, не хотите ли оформить вашу заметку на вики и опубликовать в виде статьи на ГИС-Лаб? Я уверен, кому-то бы пригодилось

Хорошо, при наличии времени - попробую, правда я не писатель

losdimas писал(а):Хорошо, при наличии времени - попробую, правда я не писатель

Здесь никто не писатель, если что. Заранее спасибо!

losdimas спасибо огромное! Сейчас как раз во второй раз возникла такая же необходимость. Заглянул в эту ветку - а тут готовое решение

PS в прошлый раз просто раскрасил точки в зависимости от их значения (через градуированный знак)