Создаём свою карту ДТП - История изменений

Максим Дубинин в 10:23, 4 мая 2018

2018-05-04T10:23:43Z

← Предыдущая версия		Версия от 12:23, 4 мая 2018
Строка 266:		Строка 266:

	==Считаем ДТП и смерти==		==Считаем ДТП и смерти==
			QGIS

	==Делаем веб-карту==		==Делаем веб-карту==
			nextgis.com

			Подробнее http://gis-lab.info/qa/nextgiscom.html

Максим Дубинин: /* Преобразовываем в табличный вид */

2018-05-04T10:21:20Z

Преобразовываем в табличный вид

← Предыдущая версия		Версия от 12:21, 4 мая 2018
Строка 11:		Строка 11:

	==Преобразовываем в табличный вид==		==Преобразовываем в табличный вид==
	С XML работать не удобно, поэтому используем следующий скрипт для преобразования XML в CSV.		С XML работать не удобно, поэтому используем следующий скрипт для преобразования XML в CSV. ДОБАВИТЬ СКРИПТ

	Результат: CSV файл со всеми данными по ДТП.		Результат: CSV файл со всеми данными по ДТП.

Максим Дубинин: /* Загружаем исходные данные */

2018-05-04T10:21:10Z

Загружаем исходные данные

← Предыдущая версия		Версия от 12:21, 4 мая 2018
Строка 6:		Строка 6:
	Указываем диапазон дат.		Указываем диапазон дат.

	Запускаем скрипт извлечения данных, по годам.		Запускаем скрипт извлечения данных, по годам. ДОБАВИТЬ СКРИПТ

	Результат: пачка XML с карточками ДТП.		Результат: пачка XML с карточками ДТП.

Trolleway в 10:37, 4 апреля 2018

2018-04-04T10:37:45Z

← Предыдущая версия		Версия от 12:37, 4 апреля 2018
Строка 80:		Строка 80:
	[[Файл:Points near points_2.png]]		[[Файл:Points near points_2.png]]

	Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос.		Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос. Если вы работаете не с Москвой, то вместо 32637 подставьте номер системы координат UTM на нужную местность.

	<syntaxhighlight lang="sql">		<syntaxhighlight lang="sql">
			-- генерируем три слоя буферов с разным расстоянием
			-- Укажите тут нужную СК, если это не Москва.
	CREATE TEMPORARY TABLE buffers100 ON COMMIT DROP AS		CREATE TEMPORARY TABLE buffers100 ON COMMIT DROP AS
	SELECT st_buffer(ST_Transform(wkb_geometry,32637),100) AS wkb_geometry,100 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;		SELECT st_buffer(ST_Transform(wkb_geometry,32637),100) AS wkb_geometry,100 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;
Строка 96:		Строка 98:
	wkb_geometry geometry,		wkb_geometry geometry,
	buffer integer);		buffer integer);

			--для 100 метров
	DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;		DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
	CREATE		CREATE
Строка 113:		Строка 116:
	INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 100 as buffer FROM crossingsarea_disolved_onebuffer;		INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 100 as buffer FROM crossingsarea_disolved_onebuffer;

			--для 50 метров
	DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;		DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
	CREATE		CREATE
Строка 130:		Строка 133:
	INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 50 as buffer FROM crossingsarea_disolved_onebuffer;		INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 50 as buffer FROM crossingsarea_disolved_onebuffer;

			--для 20 метров
	DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;		DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
	CREATE		CREATE
Строка 148:		Строка 151:


			--добавление столбцов и рассчёт числовых значений
	ALTER TABLE crossingsarea_disolved ADD COLUMN count_2014 integer;		ALTER TABLE crossingsarea_disolved ADD COLUMN count_2014 integer;
	ALTER TABLE crossingsarea_disolved ADD COLUMN count_2015 integer;		ALTER TABLE crossingsarea_disolved ADD COLUMN count_2015 integer;
Строка 183:		Строка 187:
	WHERE		WHERE
	ST_Within(dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry);		ST_Within(dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry);



Строка 251:		Строка 254:
	WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;		WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
	DROP TABLE medium2016_cnt;		DROP TABLE medium2016_cnt;


	UPDATE crossingsarea_disolved SET cp_sum = cp_2014+cp_2015+cp_2016;		UPDATE crossingsarea_disolved SET cp_sum = cp_2014+cp_2015+cp_2016;

Trolleway в 10:27, 4 апреля 2018

2018-04-04T10:27:38Z

← Предыдущая версия		Версия от 12:27, 4 апреля 2018
Строка 77:		Строка 77:

	Для перекрёстков строится 3 буфера, они объединяются. Делается spatial self join, и подсчитывается количество и разные суммы атрибутов аварий у перекрёстков. На выходе - полигональный слой с кружками вокруг пеерекрёстков и разными атрибутами.		Для перекрёстков строится 3 буфера, они объединяются. Делается spatial self join, и подсчитывается количество и разные суммы атрибутов аварий у перекрёстков. На выходе - полигональный слой с кружками вокруг пеерекрёстков и разными атрибутами.

	[[Файл:Points near points_2.png]]		[[Файл:Points near points_2.png]]

	Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос.		Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос.

Строка 254:		Строка 256:
	</syntaxhighlight>		</syntaxhighlight>

			Выгружаем слой в файл
			<syntaxhighlight lang="bash">
			ogr2ogr -progress -overwrite -f GPKG crossings_dtp.gpkg PG:"host=localhost dbname=gis " crossingsarea_disolved -skipfailures
			</syntaxhighlight>

	==Считаем ДТП и смерти==		==Считаем ДТП и смерти==

	==Делаем веб-карту==		==Делаем веб-карту==

Trolleway в 10:22, 4 апреля 2018

2018-04-04T10:22:58Z

← Предыдущая версия		Версия от 12:22, 4 апреля 2018
Строка 77:		Строка 77:

	Для перекрёстков строится 3 буфера, они объединяются. Делается spatial self join, и подсчитывается количество и разные суммы атрибутов аварий у перекрёстков. На выходе - полигональный слой с кружками вокруг пеерекрёстков и разными атрибутами.		Для перекрёстков строится 3 буфера, они объединяются. Делается spatial self join, и подсчитывается количество и разные суммы атрибутов аварий у перекрёстков. На выходе - полигональный слой с кружками вокруг пеерекрёстков и разными атрибутами.
			[[Файл:Points near points_2.png]]
	Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос.		Подключитесь к БД в pgAdmin или Valentina Studio, выполните этот запрос.

Trolleway в 10:16, 4 апреля 2018

2018-04-04T10:16:59Z

Внесённые изменения

Trolleway в 09:59, 4 апреля 2018

2018-04-04T09:59:31Z

← Предыдущая версия		Версия от 11:59, 4 апреля 2018
Строка 48:		Строка 48:

	Скрипт создаст Shapefile с точками, в атрибуте WAYS_CNT будет количество дорог, подходящих к точке.		Скрипт создаст Shapefile с точками, в атрибуте WAYS_CNT будет количество дорог, подходящих к точке.

	[[Файл:Crossing detection result.png]]		[[Файл:Crossing detection result.png]]

Trolleway: /* Скрипт детектирования перекрёстков в данных Openstreetmap */

2018-04-04T09:59:22Z

Скрипт детектирования перекрёстков в данных Openstreetmap

← Предыдущая версия		Версия от 11:59, 4 апреля 2018
Строка 22:		Строка 22:

	===Скрипт детектирования перекрёстков в данных Openstreetmap===		===Скрипт детектирования перекрёстков в данных Openstreetmap===
	В OSM нет тега для перекрёстка, но их можно получить, взяв те точки, в которых пересекается несколько веев (более двух). Для пробы этот скрипт написан на питоновских биндингах GDAL, и Москву он считает за ночь. Вероятно в PostGIS этот же алгоритм можно сделать более быстрым.		В OSM нет тега для перекрёстка, но их можно получить, взяв те точки, в которых пересекается несколько веев (более двух). Для пробы этот скрипт был написан на питоновских биндингах GDAL, и Москву он считает за ночь. Вероятно в PostGIS этот же алгоритм можно сделать более быстрым.

	Если вы работаете на Windows или Mac OS, то установите NextGIS QGIS, найдите и запустите в меню программ "NextGIS command prompt".		Если вы работаете на Windows или Mac OS, то установите NextGIS QGIS, найдите и запустите в меню программ "NextGIS command prompt".

Trolleway в 09:58, 4 апреля 2018

2018-04-04T09:58:19Z

← Предыдущая версия		Версия от 11:58, 4 апреля 2018
Строка 45:		Строка 45:
	python detect_crossings.py -s highway-line.shp -d crossings.shp -f "HIGHWAY IN ('motorway','motorway_link','trunk','trunk_link','primary','primary_link','secondary','secondary_link','tertiary','tertiary_link','unclassified','residential')"		python detect_crossings.py -s highway-line.shp -d crossings.shp -f "HIGHWAY IN ('motorway','motorway_link','trunk','trunk_link','primary','primary_link','secondary','secondary_link','tertiary','tertiary_link','unclassified','residential')"
	</syntaxhighlight>		</syntaxhighlight>
			Фильтр нужен обязательно, потому что без него у вас получатся перекрёстки тротуаров с тропинками в лесу

	[[Файл:Crossing detection result.~~svg~~]]		Скрипт создаст Shapefile с точками, в атрибуте WAYS_CNT будет количество дорог, подходящих к точке.
			[[Файл:Crossing detection result.png]]

	==Создаём перекрестки==		==Создаём перекрестки==