Вызов GRASS GIS из скрипта на Python - История изменений

Максим Дубинин: /* Реализация */

2015-08-26T21:04:47Z

Реализация

← Предыдущая версия		Версия от 23:04, 26 августа 2015
Строка 60:		Строка 60:
	import grass.script as grass		import grass.script as grass
	import grass.script.setup as gsetup		import grass.script.setup as gsetup
	gsetup.init(gisbase,		gsetup.init(gisbase, gisdbase, location, mapset)
	gisdbase, location, mapset)
	print grass.gisenv()		print grass.gisenv()

Строка 67:		Строка 66:

	Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде ''gsetup.init''. Обратите внимание на то, что использование команды ''gsetup.init'' избавляет программиста от создания/редактирования файла путей к проекту (см. переменную GISRC).		Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде ''gsetup.init''. Обратите внимание на то, что использование команды ''gsetup.init'' избавляет программиста от создания/редактирования файла путей к проекту (см. переменную GISRC).

			Область (location) заданная до ''gsetup.init'' должна существовать, иначе инициализации не произойдет.

	В Windows при установке GRASS через OSGeo4W/NextGIS_QGIS набор переменных может выглядеть так:		В Windows при установке GRASS через OSGeo4W/NextGIS_QGIS набор переменных может выглядеть так:

KolesovDmitry: /* Запуск команд */

2015-06-22T07:16:35Z

Запуск команд

← Предыдущая версия		Версия от 09:16, 22 июня 2015
Строка 93:		Строка 93:
	Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными. Обратите внимание, что флаги (например --o перезапись) передаются в flags.		Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными. Обратите внимание, что флаги (например --o перезапись) передаются в flags.

	Для обработки выводов команд GRASS может понадобится работа с пайпами, например для получения списка растров области в виде массива:		Для обработки выводов команд GRASS может понадобится работа с каналами (пайпами), например для получения списка растров области в виде массива:

	<syntaxhighlight lang="python">		<syntaxhighlight lang="python">

Максим Дубинин: /* Запуск команд */

2015-06-15T19:51:23Z

Запуск команд

← Предыдущая версия		Версия от 21:51, 15 июня 2015
Строка 100:		Строка 100:
	for line in p.stdout:		for line in p.stdout:
	result.append(line)		result.append(line)
			</syntaxhighlight>

			Для некоторых команд GRASS параметры передаются не в виде "ключ - значение", а простым аргументом, в этом случае запуск осуществляется не через run_command, а саму команду, например для r.mapcalc:

			<syntaxhighlight lang="python">
			grass.mapcalc('dh3 = std/ave')
	</syntaxhighlight>		</syntaxhighlight>

Максим Дубинин: /* Обработка данных */

2015-06-15T18:24:11Z

Обработка данных

← Предыдущая версия		Версия от 20:24, 15 июня 2015
Строка 77:		Строка 77:
	Запускать python в этом случае нужно после OSGeo4W.bat.		Запускать python в этом случае нужно после OSGeo4W.bat.

	== ~~Обработка данных~~ ==		== Запуск команд ==
	Здесь происходит собственно обработка данных, которая зависит от целей и задач пользователя. Пользователь может вызывать любые команды GRASS или запускать сторонние модули Python или написать свой собственный обработчик.		Здесь происходит собственно обработка данных, которая зависит от целей и задач пользователя. Пользователь может вызывать любые команды GRASS или запускать сторонние модули Python или написать свой собственный обработчик.

Строка 91:		Строка 91:
	</syntaxhighlight>		</syntaxhighlight>

	Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными.		Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными. Обратите внимание, что флаги (например --o перезапись) передаются в flags.

	Для обработки выводов команд GRASS может понадобится работа с пайпами, например для получения списка растров области в виде массива:		Для обработки выводов команд GRASS может понадобится работа с пайпами, например для получения списка растров области в виде массива:

Максим Дубинин: /* Итоговый скрипт */

2015-06-15T16:40:40Z

Итоговый скрипт

← Предыдущая версия		Версия от 18:40, 15 июня 2015
Строка 109:		Строка 109:

	<syntaxhighlight lang="python">		<syntaxhighlight lang="python">
	#!/usr/bin/env python		#!/usr/bin/env python
	# -- coding: utf-8 --		# -- coding: utf-8 --

	import os		import os
	import sys		import sys
	gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"		gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"
	gisdbase = os.path.join(os.environ['HOME'], "grassdata")		gisdbase = os.path.join(os.environ['HOME'], "grassdata")
	location = "geosample"		location = "geosample"
	mapset = "PERMANENT"		mapset = "PERMANENT"
	sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))		sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))
	import grass.script as grass		import grass.script as grass
	import grass.script.setup as gsetup		import grass.script.setup as gsetup
	gsetup.init(gisbase,		gsetup.init(gisbase,
	gisdbase, location, mapset)		gisdbase, location, mapset)
	print grass.gisenv()		print grass.gisenv()

	# Обработка		# Обработка
	grass.message('Raster maps:')		grass.message('Raster maps:')
	for rast in grass.list_strings(type = 'rast'):		for rast in grass.list_strings(type = 'rast'):
	print rast		print rast
	grass.run_command('r.stats', input=rast, flags='a')		grass.run_command('r.stats', input=rast, flags='a')
	</syntaxhighlight>		</syntaxhighlight>

Максим Дубинин: /* Обработка данных */

2015-06-15T16:40:14Z

Обработка данных

← Предыдущая версия		Версия от 18:40, 15 июня 2015
Строка 85:		Строка 85:

	<syntaxhighlight lang="python">		<syntaxhighlight lang="python">
	grass.message('Raster maps:')		grass.message('Raster maps:')
	for rast in grass.list_strings(type = 'rast'):		for rast in grass.list_strings(type = 'rast'):
	print rast		print rast
	grass.run_command('r.stats', input=rast, flags='a')		grass.run_command('r.stats', input=rast, flags='a')
	</syntaxhighlight>		</syntaxhighlight>

	Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными.		Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными.

			Для обработки выводов команд GRASS может понадобится работа с пайпами, например для получения списка растров области в виде массива:

			<syntaxhighlight lang="python">
			result = []
			p = grass.pipe_command('g.mlist', type = 'rast', quiet=True)
			for line in p.stdout:
			result.append(line)
			</syntaxhighlight>

	== Очистка системы ==		== Очистка системы ==

Максим Дубинин: /* Реализация */

2015-06-12T21:23:58Z

Реализация

← Предыдущая версия		Версия от 23:23, 12 июня 2015
Строка 65:		Строка 65:

	</syntaxhighlight>		</syntaxhighlight>

			Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде ''gsetup.init''. Обратите внимание на то, что использование команды ''gsetup.init'' избавляет программиста от создания/редактирования файла путей к проекту (см. переменную GISRC).

	В Windows при установке GRASS через OSGeo4W/NextGIS_QGIS набор переменных может выглядеть так:		В Windows при установке GRASS через OSGeo4W/NextGIS_QGIS набор переменных может выглядеть так:
Строка 73:		Строка 75:
	</syntaxhighlight>		</syntaxhighlight>

	Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде ''gsetup.init''. Обратите внимание на то, что использование команды ''gsetup.init'' избавляет программиста от создания/редактирования файла путей к проекту (см. ~~переменную GISRC)~~.		Запускать python в этом случае нужно после OSGeo4W.bat.

	== Обработка данных ==		== Обработка данных ==

Максим Дубинин: /* Реализация */

2015-06-12T21:22:12Z

Реализация

← Предыдущая версия		Версия от 23:22, 12 июня 2015
Строка 51:		Строка 51:

	<syntaxhighlight lang="python">		<syntaxhighlight lang="python">
	import os		import os
	import sys		import sys
	gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"		gisbase = os.environ['GISBASE'] = "/usr/lib/grass64"
	gisdbase = os.path.join(os.environ['HOME'], "grassdata")		gisdbase = os.path.join(os.environ['HOME'], "grassdata")
	location = "geosample"		location = "geosample"
	mapset = "PERMANENT"		mapset = "PERMANENT"
	sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))		sys.path.append(os.path.join(os.environ['GISBASE'], "etc", "python"))
	import grass.script as grass		import grass.script as grass
	import grass.script.setup as gsetup		import grass.script.setup as gsetup
	gsetup.init(gisbase,		gsetup.init(gisbase,
	gisdbase, location, mapset)		gisdbase, location, mapset)
	print grass.gisenv()		print grass.gisenv()

			</syntaxhighlight>

			В Windows при установке GRASS через OSGeo4W/NextGIS_QGIS набор переменных может выглядеть так:

			<syntaxhighlight lang="python">
			gisbase = os.environ['GISBASE'] = "c:/tools/NextGIS_QGIS/apps/grass/grass-6.4.4/"
			gisdbase = os.environ['GISDBASE'] = "e:/users/maxim/thematic/dhi/"
	</syntaxhighlight>		</syntaxhighlight>

Александр Мурый в 18:00, 2 декабря 2013

2013-12-02T18:00:36Z

← Предыдущая версия		Версия от 20:00, 2 декабря 2013
Строка 1:		Строка 1:
			{{Статья\|Опубликована\|grass-external-scripting}}

	Очень часто работа с GRASS GIS происходит в интерактивном режиме, когда пользователь вводит команды в консоли или вызывает их через графический интерфейс пользователя. При этом часто повторяемые действия пользователь может автоматизировать, создав скрипт, который будет доступен для повторного использования. GRASS GIS поддерживает большое количество языков программирования, на которых можно создавать пользовательские скрипты, наиболее распространенными из них являются скриптовый язык командной оболочки (shell, bash) и Python. Пример создания скрипта на языке Python описывается в статье [http://gis-lab.info/qa/grass64-python.html Создание скрипта на Python для ГИС GRASS]. После создания скрипта пользователь может запускать его из консоли, в которой запущена ГИС GRASS.		Очень часто работа с GRASS GIS происходит в интерактивном режиме, когда пользователь вводит команды в консоли или вызывает их через графический интерфейс пользователя. При этом часто повторяемые действия пользователь может автоматизировать, создав скрипт, который будет доступен для повторного использования. GRASS GIS поддерживает большое количество языков программирования, на которых можно создавать пользовательские скрипты, наиболее распространенными из них являются скриптовый язык командной оболочки (shell, bash) и Python. Пример создания скрипта на языке Python описывается в статье [http://gis-lab.info/qa/grass64-python.html Создание скрипта на Python для ГИС GRASS]. После создания скрипта пользователь может запускать его из консоли, в которой запущена ГИС GRASS.

Александр Мурый в 13:56, 2 декабря 2013

2013-12-02T13:56:02Z

← Предыдущая версия		Версия от 15:56, 2 декабря 2013
Строка 1:		Строка 1:
	Очень часто работа с GRASS GIS происходит в интерактивном режиме, когда пользователь вводит команды в консоли или вызывает их через графический интерфейс пользователя. При этом часто повторяемые действия пользователь может автоматизировать, создав скрипт, который будет доступен для повторного использования. GRASS GIS поддерживает большое количество языков программирования, на которых можно создавать пользовательские скрипты, наиболее распространенными из них являются bash и Python. Пример создания скрипта на языке ~~python~~ описывается в статье [http://gis-lab.info/qa/grass64-python.html Создание скрипта на Python для ГИС GRASS]. После создания скрипта пользователь может запускать его из консоли, в которой запущена ГИС GRASS.		Очень часто работа с GRASS GIS происходит в интерактивном режиме, когда пользователь вводит команды в консоли или вызывает их через графический интерфейс пользователя. При этом часто повторяемые действия пользователь может автоматизировать, создав скрипт, который будет доступен для повторного использования. GRASS GIS поддерживает большое количество языков программирования, на которых можно создавать пользовательские скрипты, наиболее распространенными из них являются скриптовый язык командной оболочки (shell, bash) и Python. Пример создания скрипта на языке Python описывается в статье [http://gis-lab.info/qa/grass64-python.html Создание скрипта на Python для ГИС GRASS]. После создания скрипта пользователь может запускать его из консоли, в которой запущена ГИС GRASS.

	Однако, иногда может потребоваться полностью автономный скрипт, который может работать, даже если консоль GRASS GIS не запущена, например, при вызове скрипта из какой-либо своей программы. Это может быть ~~к примеру,~~ полезным при автоматической обработке геоданных, когда нужно по расписанию обновлять, генерировать или анализировать данные. Возможно также создание таких скриптов, которые автоматически создают проект GRASS, выполняют определенную работу (импортируют данные, анализируют их) и, по завершении обработки, удаляют проекты. Методика создания подобных скриптов подробно описана в статье [http://geoinformatics.fsv.cvut.cz/gwiki/On_scripting_GRASS_GIS:_Building_location-independent_command_line_tools On scripting GRASS GIS: Building location-independent command line tools]. В данной статье рассматривается лишь часть методов, описанных вышеуказанной работе -- создание скрипта, способного использовать готовые проекты GRASS.		Однако, иногда может потребоваться полностью автономный скрипт, который может работать, даже если консоль GRASS GIS не запущена, например, при вызове скрипта из какой-либо своей программы. Это может быть полезным при автоматической обработке геоданных, когда нужно по расписанию обновлять, генерировать или анализировать данные. Возможно также создание таких скриптов, которые автоматически создают проект GRASS, выполняют определенную работу (импортируют данные, анализируют их) и, по завершении обработки, удаляют проекты. Методика создания подобных скриптов подробно описана в статье [http://geoinformatics.fsv.cvut.cz/gwiki/On_scripting_GRASS_GIS:_Building_location-independent_command_line_tools On scripting GRASS GIS: Building location-independent command line tools]. В данной статье рассматривается лишь часть методов, описанных вышеуказанной работе -- создание скрипта, способного использовать готовые проекты GRASS.

	Пример подобного скрипта, написанного на shell, приводится в [http://gis-lab.info/forum/viewtopic.php?f=40&t=9303 рецепте], данная статья описывает решение данной задачи для языка Python. Нужно отметить, что описываемый подход можно применить к любому языку программирования, а не только к shell или Python.		Пример подобного скрипта, написанного на shell, приводится в [http://gis-lab.info/forum/viewtopic.php?f=40&t=9303 рецепте], данная статья описывает решение данной задачи для языка Python. Нужно отметить, что описываемый подход можно применить к любому языку программирования, а не только к shell или Python.
Строка 48:		Строка 48:
	Предположим, что пользователь работает с набором данных [[http://gis-lab.info/qa/geosample.html geosample]], и пусть $HOME/grassdata -- каталог, в котором хранится данный набор. Тогда для установки переменных можно воспользоваться следующими командами на языке Python:		Предположим, что пользователь работает с набором данных [[http://gis-lab.info/qa/geosample.html geosample]], и пусть $HOME/grassdata -- каталог, в котором хранится данный набор. Тогда для установки переменных можно воспользоваться следующими командами на языке Python:

			<syntaxhighlight lang="python">
	import os		import os
	import sys		import sys
Строка 61:		Строка 62:
	print grass.gisenv()		print grass.gisenv()

	Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде gsetup.init. Обратите внимание на то, что использование команды gsetup.init избавляет программиста от создания/редактирования файла путей к проекту (см. переменную GISRC).		</syntaxhighlight>

			Третья строка задает переменную окружения GISBASE. Четвертая, пятая и шестая строки задают переменные путей к проекту GRASS, которые будут использоваться в команде ''gsetup.init''. Обратите внимание на то, что использование команды ''gsetup.init'' избавляет программиста от создания/редактирования файла путей к проекту (см. переменную GISRC).

	== Обработка данных ==		== Обработка данных ==
Строка 69:		Строка 72:
	# Выводит называние растра на экран		# Выводит называние растра на экран
	# Выводит статистику по данному растру (значение ячейки растра и общая площадь ячеек с данным значением)		# Выводит статистику по данному растру (значение ячейки растра и общая площадь ячеек с данным значением)

			<syntaxhighlight lang="python">
	grass.message('Raster maps:')		grass.message('Raster maps:')
	for rast in grass.list_strings(type = 'rast'):		for rast in grass.list_strings(type = 'rast'):
	print rast		print rast
	grass.run_command('r.stats', input=rast, flags='a')		grass.run_command('r.stats', input=rast, flags='a')
			</syntaxhighlight>

	Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными.		Читатель легко может заменить эти команды на свои, отражающие суть необходимых ему манипуляций с данными.
Строка 81:		Строка 87:
	== Итоговый скрипт ==		== Итоговый скрипт ==
	Окончательный вид скрипта, готового для использования, приводится ниже:		Окончательный вид скрипта, готового для использования, приводится ниже:

			<syntaxhighlight lang="python">
	#!/usr/bin/env python		#!/usr/bin/env python
	# -- coding: utf-8 --		# -- coding: utf-8 --
Строка 102:		Строка 110:
	print rast		print rast
	grass.run_command('r.stats', input=rast, flags='a')		grass.run_command('r.stats', input=rast, flags='a')
			</syntaxhighlight>

	== Ссылки ==		== Ссылки ==
	* [http://grasswiki.osgeo.org/wiki/Working_with_GRASS_without_starting_it_explicitly Working with GRASS without starting it explicitly (англ.)]		* [http://grasswiki.osgeo.org/wiki/Working_with_GRASS_without_starting_it_explicitly Working with GRASS without starting it explicitly (англ.)]
	* [http://gis-lab.info/forum/viewtopic.php?f=40&t=9303 Рецепт создания внешнего для ГИС GRASS скрипта на shell]		* [http://gis-lab.info/forum/viewtopic.php?f=40&t=9303 Рецепт создания внешнего для ГИС GRASS скрипта на shell]