GIS-Lab - Вклад [ru]

Свободные данные по границам субъектов РФ

2026-01-09T12:29:02Z

Максим Дубинин: новые ссылки

{{Статья|Опубликована|rusbounds-rosreestr}}
{{Аннотация|Результаты коллективно проекта по созданию слоя границ субъектов РФ.}}

Эта статья предоставляет краткое описание и ссылки для загрузки открытого слоя границ субъектов Российской Федерации, созданных в результате открытого коллективного проекта на базе GIS-Lab. С подробной информацией о проекте в стиле вопрос-ответ можно ознакомиться на [http://gis-lab.info/qa/rusbounds-rosreestr-faq.html отдельной странице]. Координационную страницу проекта можно посмотреть [http://gis-lab.info/projects/rusbounds-rosreestr.html здесь].

Этот проект был бы невозможен без его участников (сортировка по алфавиту, в скобках номера сделанных участниками фрагментов и их дополнительная деятельность по проекту):

'''_DR_''' (130, 172, 118, 119, 120, 122, 101, 125, 099, 145, 150, ,153, 156, 160, 162, 124, 073, 013, 014, 030, 035, 036, 203, 063, 102, 074, 075, 076, 077, 088, 098, 045, 200, 178, 181, 179, 191, 202, 190), '''AkMeR''' (064, 169, 183, 182), '''Anna''' (подготовка слоёв с поворотными точками),''' alexandr cherepanov''' (142, 143, 144), '''amuriy''' (010, 011, 078, 009, 079, 001, 008, 201, 012, 163, 171, 146, 164, 027, 067, 154, 155), '''bim2010''' (157, 189), '''esasha_''' (165, 166), '''Igg''' (087), '''kabanin''' (193, 196, 187, 066, 188, 186, 040, 039, 192, 199, 167, 168), '''KolesovDmitry''' (041, 061, 065), '''Lileo''' (005, 002, 131, 135, 134, 129, 003, 007, 006, 004, 132), '''lzhl''' (140), '''Minos''' (141, 016, 017, 018, 022, 020, 059, 058, 021, 123, 180), '''Nikolaich''' (060, 070, 151, 170), '''old_Bibigon''' (033, 062, 034, 112, 069, 043, 137, 038, 096, 095, 094, 037, 136, 138, 139), '''outmaker''' (028, 029, 032), '''prial78''' (100, 107, 104, 108, 105, 106, 103), '''saxa-tortile''' (115, 113, 117, 116, 158, 152, 159), '''sim''' (047, 048, 042, 046, 044, 049, 024, 097, 015, 109, 110, 111, 176, 177, 025, 026, 175, 174, 023), '''SS_Rebelious''' (031, 147, 148, 149, 114, общая координация и редакция), '''taimyr''' (081, 133, 128, 093, 092, 091, 090, 089, 086, 084, 085, 082, 198, 080, 055, 050, 051, 052, 083, 054, 194, 056, 057, 197, 072, 195, 053, сведение в слои), '''toxa''' (068, 126, 127), '''Trace''' (071), '''vveco''' (173, 185, 184, 121).

==О данных==

До недавнего момента полностью открытый слой границ субъектов РФ отсутствовал. Ряд существующих источников не отвечал одному или нескольким из условий открытости и качества:

* Неясное правообладание или нечетко описанные права использования данных (Data+, ESRI);
* Неизвестный источник (GADM);
* Запрет на коммерческое использование (GADM);
* Недостаточно хорошее качество (OSM);

Этот слой распространяется по лицензии [http://creativecommons.org/licenses/by/3.0/ CC-BY], разрешающей любое использование, включая передачу, изменение, использование в коммерческих проектах, с условием упоминания авторства и наличием гиперссылки. Предлагаемое цитирование данных:

<pre>Росреестр, GIS-Lab.info</pre>

==Использованные источники==

* Границы субъектов- публичные данные [https://rosreestr.ru/wps/portal/cc_ib_list_borders_cadastral_districts?param_infoblock_document_path=infoblock-root/cc_ib_list_borders_cadastral_districts/index.htm Росреестра](старая ссылка: http://www.rosreestr.ru/cadastre/ground_cadastre/border);
* Сухопутная граница России (кроме Калининградской области) - [http://gis-lab.info/data/osm/russia_full/russia_full-admin-l.7z OSM] от 31.03.2010
* Граница г. Москва - [http://gis-lab.info/data/osm/mosobl/mosobl-admin-l.7z OSM] от 04.04.2010 г., выборка Москва-регион
* Береговая линия - [http://gis-lab.info/qa/vmap0.html VMAP0] (bnd-coast-l.shp, "закрыты" устья рек)
* Границы Калиниградской области - [http://gis-lab.info/qa/vmap0.html VMAP0]

==Получение данных==

===Полигональные слои===

Полигональный слой субъектов РФ по состоянию на 2010 г. - современное состояние административно-территориального деления РФ, только сухопутные участки. Чукотский АО разрезан по 180 градусу.

{| class="wikitable" align="center"
|
| Система координат/проекция
| Shape
| MIF
| TAB
| KML
|-
| regions2010_wgs
| Lat/Lon WGS 1984
|
[https://disk.yandex.ru/d/haz4nXE9R8xXaQ скачать]
|
[https://disk.yandex.ru/d/aEp-ymVlvgGyxg скачать]
|
[https://disk.yandex.ru/d/Mv0e8DjxYa0Lrw скачать]
|
[https://disk.yandex.ru/d/FQVCl3tZQsDvEw скачать]
|-
| regions2010_pul
| Lat/Lon Pulkovo 1942
|
[https://disk.yandex.ru/d/XfcfyUow9nuziQ скачать]
|
[https://disk.yandex.ru/d/1sByatcosAc7CA скачать]
|
[https://disk.yandex.ru/d/nrAYIiYoqCL4kA скачать]
|
|-
| regions2010_alb
|
Albers-Siberia ([http://gis-lab.info/qa/gis-lab-projections.html параметры])
|
[https://disk.yandex.ru/d/H_bZhvMK8LV_Jw скачать]
|
[https://disk.yandex.ru/d/pb-c93XtO2J_zQ скачать]
|
[https://disk.yandex.ru/d/AnxSkRL7amQTxg скачать]
|
|}

Полигональный слой субъектов РФ по состоянию на 2000 г. - до объединения регионов, включает также акватории Архангельской области, Республики Карелия и г. Санкт-Петербурга:

{| class="wikitable" align="center"
|
| Система координат/проекция
| Shape
| MIF
| TAB
| KML
|-
| regions2000_wgs
| Lat/Lon WGS 1984
|
[https://disk.yandex.ru/d/zIoOS4KvN6w-rg скачать]
|
[https://disk.yandex.ru/d/0z2gzWz3i97FWA скачать]
|
[https://disk.yandex.ru/d/ouIiBGjF6u8FpQ скачать]
|
[https://disk.yandex.ru/d/hgvVdAalsz70iw скачать]
|-
| regions2000_pul
| Lat/Lon Pulkovo 1942
|
[https://disk.yandex.ru/d/IkRFWjKc4hbxMg скачать]
|
[https://disk.yandex.ru/d/515X74BL8Clumg скачать]
|
[https://disk.yandex.ru/d/AxXDFzq3N2Lawg скачать]
|
|-
| regions2000_alb
|
Albers-Siberia ([http://gis-lab.info/qa/gis-lab-projections.html параметры])
|
[https://disk.yandex.ru/d/1p8oRzFia4Ahqg скачать]
|
[https://disk.yandex.ru/d/Awt7CONDTDC5Zw скачать]
|
[https://disk.yandex.ru/d/IKvlc74KEbo8XQ скачать]
|
|}

===Линейные слои===

Все границы одним слоем. Отсутствуют границы г. Москва и Чеченской и Ингушской республик по причине отсутствия по ним данных Росреестра.

{| class="wikitable" align="center"
|
| Система координат/проекция
| Shape
| OSM
| MIF
| TAB
| KML
|-
| all_rubounds_pul
| Lat/Lon Pulkovo 1942
|
[https://disk.yandex.ru/d/Q-jRzrCQZrUesA скачать]
|
[https://disk.yandex.ru/d/on0jTQpFtJ-kzw скачать]
|
[https://disk.yandex.ru/d/R2hovdXfWvWHZQ скачать]
|
[https://disk.yandex.ru/d/7l3X-kdbuACMwA скачать]
|
[https://disk.yandex.ru/d/mWEm-zDl6RnU7Q скачать]
|-
| all_rubounds_wgs
| Lat/Lon WGS 1984
|
[https://disk.yandex.ru/d/tHnCneVDxuu2PQ скачать]
|
[https://disk.yandex.ru/d/xdrPfUTf6DsLcg скачать]
|
|
|
|}

Атрибутивные поля:

* '''№''' - №№ границы
* '''Region1''' - название 1-го региона
* '''Code_r1''' - код 1-го региона
* '''Region2 '''- название 2-го региона
* '''Code_r2''' - код 2-го региона
* '''Status''' - статус границы (границы регионов, объединенных в период 2004-2009 гг., имеют статус munic, остальные - subj)
* '''Rastr''' - гиперрсылка на первичные растровые схемы
* '''Descript''' - гиперссылка на описание границ
* '''Rastr_georef''' - гиперрсылка на геопривязанные растры (СК Lat/Lon Pulkovo 1942)
* '''Points''' - гиперрсылка на шейп-файл поворотных точек границы
* '''Vect_bound''' - гиперссылка на shape-файл границы
* '''Autor''' - автор
* '''Autor_pr''' - гиперссылка на профиль автора на форуме gis-lab
* '''Corrected_by '''- отмечены границы, допривязанные и откорректированные по дополнительным источникам (мозаики ЕТМ+ и др.)

===PostGIS===

Загрузить слой субъектов можно из базы данных PostGIS: сервер: gis-lab.info, БД: gen, имя/пароль: guest/guest, слой: oblasts ([http://gis-lab.info/qa/openbase.html подробнее]).

==Ссылки по теме==
*[http://gis-lab.info/qa/rusbounds-rosreestr.html Свободные данные по границам субъектов РФ - вопросы и ответы]
*[http://gis-lab.info/qa/vmap0-settl-rus.html Русскоязычная база данных по названиям населенных пунктов]

Перевод QGIS на русский язык: сообщество и руководство для участника

2025-12-17T08:29:08Z

Максим Дубинин: /* Хочу участовать, что мне понадобится и как подготовиться? */

{{Статья|Опубликована|qgis-trans}}
{{Аннотация|Материал посвящен процессу организованного перевода открытой ГИС QGIS на русский язык. Здесь вы найдёте общую информацию о том, как организован перевод, все необходимые ссылки, а также руководство по тому, как стать участником и полезные советы}}

Популярность [https://qgis.org/ru/site/ QGIS] на сегодняшний день впечатляет: эта открытая ГИС используется в индустриальных, научных и образовательных организациях по всему миру. Если вы как-то связаны с миром геоинформационных технологий, то, скорее всего, вам доводилось если не работать с QGIS, то по крайней мере сталкиваться с проектами и инфраструктурами пространственных данных на её основе.

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

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

Перед тем, как приступить к переводу, просим вас полностью прочитать эту статью.

== Как вообще организована локализация QGIS? ==

Перевод на разные языки осуществляется добровольцами (иногда спонсируемыми государственными или частными организациями). Стать таким добровольцем может любой человек, например вы. Для удобства процесса перевода разработчики предлагают использовать специализированное веб-приложение [https://transifex.com transifex.com], которое позволяет редактировать перевод непосредственно в интернет-браузере, без утомительных настроек окружения и установки дополнительных программ. Единственное, что вам нужно - подать заявку на перевод проекта QGIS в transifex, и дождаться её одобрения. Мы предлагаем вам не просто сразу начать перевод, но сделать это более организовано, так, как мы предлагаем в этой статье.

Вы можете познакомиться с [https://qgis.org/ru/site/getinvolved/translate.html описанием процесса локализации] на сайте QGIS.

== Кто организует систематизированный перевод на русский язык? ==

На данный момент перевод координируют:
* Эдуард Казаков (ee.kazakov@gmail.com, t.me/kazakov_e)

Если у вас есть предложения по процессу локализации, свяжитесь с нами по указанным контактам.

== Хочу участовать, что мне понадобится и как подготовиться? ==

Если вы хотите просто сообщить об обнаруженном плохом/несуществующем переводе, но не хотите полноценно участовать в процессе, просто напишите на адрес электронной почты ee.kazakov@gmail.com или в Telegram @kazakov_e / https://t.me/qgis_translate_ru

Для того, чтобы начать работать с нами, вам потребуется:
* Ночная или последняя (не LTR) версия QGIS на компьютере (дистрибутивы и инструкции для всех операционных систем [https://qgis.org/ru/site/forusers/download.html здесь])
* Telegram аккаунт
* Аккаунт в системе [https://transifex.com transifex.com] (регистрация бесплатная)

Если всё это у вас есть, действуйте по следующему алгоритму:
# Подайте заявку на перевод Russian на странице проекта QGIS в transifex: https://app.transifex.com/qgis/QGIS/ (при этом нужно быть авторизованным в transifex).
# Вступите в телеграм-чат https://t.me/qgis_translate_ru и представтесь.
# Напишите в телеграм-чате, что подали заявку в transifex.
# Вас оповестят, когда заявка будет одобрена.
# Переходите на страницу проекта https://app.transifex.com/qgis/QGIS/dashboard/, выбирайте ветку "qgis-application", и можно приступать к работе!

== Как получить актуальную версию перевода и увидеть её в QGIS? ==

Во-первых, скачайте и установите последнюю (не LTR) версию QGIS, доступную на [https://qgis.org/ru/site/forusers/download.html# официальном сайте].

Затем на transifex проделайте следующие действия:
# Авторизуйтесь и перейдите на страницу проекта https://www.transifex.com/qgis/QGIS/dashboard/
# В списке языков найдите Russian, наведите на него мышь и нажмите "Обзор ресурсов"
# В появившемся списке нажимайте на qgis-application
# Во всплывшем окне выбирайте пункт "Скачать файл для использования" - загрузится файл формата .ts

[[Файл:Qgis_translate_resource_download_options.png|600px|thumb|center]]

Далее загруженный файл нужно скомпилировать в формат .qm с помощью программ Qt Linguist. Она устанавливается для Win и Mac отсюда: https://github.com/lelegard/qtlinguist-installers/releases, в Linux устанавливается вместе с QGIS. Запустив Linguist, откройте в ней загруженный файл .ts и скомпилируйте через меню File - Release As: достаточно просто указать путь до выходного файла.

Выходной файл нужно назвать qgis_ru.qm и сохранить в папку i18n в ресурсах QGIS, например C:\Program Files\QGIS\apps\qgis\i18n

Перезапустив QGIS, актуальный перевод подхватится.

== Процесс работы переводчика ==

Перед тем, как начать перевод, выберите часть интерфейса, которую хотите перевести. Мы не рекомендуем хаотично браться за перевод отдельно взятых фраз, гораздо продуктивнее целиком переводить раздел за разделом. Уведомите команду в телеграм-чате о том, что вы начали работать над определенной частью интерфейса. Когда начинается спринт, список разделов предлагается организаторами.

Зайдя на страницу проекта в transifex, начинайте перевод ветки "qgis-application": находите фразы, которые хотите перевести или исправить, в строке поиска, и вводите свою версию перевода (некоторые советы по поиску фраз приведены далее в разделе "Практические советы по работе с transifex и важные замечания").

Если вы хотите привлечь внимание проверяющих к определенным строкам (например, они давно хорошо переведены, но не помечены как "проверенные"), то напишите в чат номера таких строк. Номера копируются из текущего адреса в браузере. Например, у вас выделена хорошо переведенная фраза, и в адресной строке вы видите https://www.transifex.com/qgis/QGIS/translate/#ru/qgis-application/52740834, значит нужно скопировать число 52740834.

Строки, которые вы сами изменяли, копировать не нужно - их мы сможем найти по фильтрам с ником и датами. Числа нужны только чтобы пометить хорошие, но не помеченные "проверенными" переводы, которые вы не трогали (не изменяли).

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

== Практические советы по работе с transifex и важные замечания ==

Основное время вы будете проводить в интерфейсе transifex, он предоставляет инструменты по поиску фраз и редактированию/созданию переводов для них. Фраза - это отдельный текст в интерфейсе, он может быть как простейшим (из одного слова, например layer), так и огромным (например, текст помощи для инструмента в панели анализа). Вы можете искать фразы по:
* Их оригинальному написанию (или фрагменту)
* По их существующему переводу (или фрагменту)
* По файлу с кодом, где они записаны
* По контексту (именованной области программного интерфейса QGIS)
* По множеству других метаданных

За поиск отвечает вот эта панель:

[[Файл:Qgis_translate_search_panel.png|600px|thumb|center]]

В выпадающих меню, видных на изображении выше, вы можете управлять, каким фильтром хотите воспользоваться сейчас. При этом имя фильтра попадёт в сам текст поиска, вот, к примеру, как будет выглядеть поиск по переводу:

[[Файл:Qgis_translate_translation_search.png|600px|thumb|center]]

По мере формулирования условий поиска будет формироваться набор фраз, доступных для перевода, в панели прямо под панелью условий поиска:

[[Файл:Qgis_translate_search_result.png|600px|thumb|center]]

Выбрав фразу кликом левой кнопки мыши, в правой части интерфейса вы увидите инструменты работы с фразой, в частности поле дял ввода перевода, а также раздел "подробности", в котором можно увидеть контекст текста и некоторые другие детали. По этим деталям обычно вы можете понять, является ли эта фраза той, которую вы ищете. Например в примере на скриншоте виден контекст AlgorithmClassification - по нему понятно, откуда фраза. Это важно, так как похожие или одинаковые фразы могут встречаться в разных частях интерфейса программы. Введя перевод, нажимайте кнопку "сохранить перевод".

Приступая к переводу конкретного раздела интерфейса, удобно, увидев название контекста на примере одной из фраз, настроить фильтр поиска на этот контекст: так вы будете видеть все связанные с этой частью интерфейса фразы (это работает не всегда - иногда контекст очень абстрактен). Для поиска по контексту вводите в строку поиска текст вида

context:<имя контекста>, например context:AlgorithmClassification

[[Файл:Qgis_translate_context_search.png|600px|thumb|center]]

Среди фраз вы можете встретить такие, которые содержат знак & между буквами - это ссылка на горячую клавишу. Если вы хотите, чтобы в переводном варианте также работала горячая клавиша, то можете поставить & перед нужной буквой.

Также среди фраз возможно встретить такую, в которой используются html-теги (выглядят, как номера на цветных бирках). Эти теги необходимо сохранять в переводной версии фразы. Самый простой способ это сделать, это скопировать оригинальный текст фразы с помощью кнопки "копировать исходную строку", а затем в копии, не трогая теги, перевести только текст.

[[Файл:Qgis_translate_tags.png|400px|thumb|center]]

При возникновении любых вопросов пишите в телеграм-чат команды переводчиков.

== Общие рекомендации ==
* '''Буква ё'''. Мы стараемся использовать букву ё во всех словах русского языка, в которые она входит.

== Перевод списков ключевых слов (тегов) ==

У инструментов панели анализа и редактора выражений есть теги - перечни ключевых слов, по которым их можно быстро искать. Выглядят, например, так:

export,add,information,measurements,areas,lengths,perimeters,
latitudes,longitudes,x,y,z,extract,points,lines,polygons,sinuosity,fields

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

* В поле перевода копируются все исходные англоязычные теги
* К ним через запятую добавляются переводные версии всех тегов. Для примера выше перевод будет следующим:

''export,add,information,measurements,areas,lengths,perimeters,
latitudes,longitudes,x,y,z,extract,points,lines,polygons,sinuosity,
fields,экспорт,добавить,информация,измерение,площадь,длина,периметр,
широта,долгота,точки,линии,полигоны,поля''

== Переводы некоторых спорных ключевых терминов, которые необходимо соблюдать ==

Здесь будет пополняться список '''дискуссионных''' терминов, которые встречаются повсеместно, и нужна единая стратегия для их перевода.

Если вы не согласны, будем рады обсудить в Telegram-чате https://t.me/qgis_translate_ru и выработать общее решение.

* '''Authentication - Аутентификация''' (мотивация: термин '''аутентификация''' закреплен в ГОСТ Р ИСО/МЭК 9594-8-98)
* '''Extent - Охват'''. Отброшенные альтернативы: граница, экстент
* '''Folder - Каталог'''. Не "папка" и не "директория".
* '''Plugin - Модуль'''. Не "Плагин".
* '''px''' (сокращенное название единицы измерения "пиксель") '''- пикс'''.
* '''Script - Сценарий'''. Не "Скрипт".
* '''Symbol - Символ''' / Symbology - символизация (мотивация: более общая терминология. Привычна ГИС-пользователям в РФ, т.к. используется в ArcGIS много лет). Отброшенные альтернативы: знак.
* '''Tile - Тайл''' (мотивация: термин вошёл в полноценный обиход, даже в гос. контрактах значатся термины вида "тайловый кэш", "векторные тайлы" и подобные им). Отброшенные альтернативы: блок мозаики.
* '''Snapping - Прилипание''' (мотивация: так всегда было в QGIS, и вся документация давно так описывает процесс. Привязка больше ассоциируется с привязкой растров/векторов). Отброшенные альтернативы: привязка.
* '''URL - URL''' (мотивация: аббревиатура закреплена в ГОСТ Р 7.0.5-2008. Понятие URL шире чем адрес). Отброшенные альтернативы: адрес.
* '''Editing mode - Режим редактирования''', '''Edits - Изменения'''. Отброшенные альтернативы: режим правок, правки.
* '''Merge - Слияние''', не объединение
* '''Union - Объединение'''
* '''Data provider - провайдер данных''', не источник данных, иначе возникает конфликт с также распространенным понятием data source
* '''nodata - NoData''', не нет данных, НетДанных и прочие вариации.

== Разделы и термины, требующие уборки/внимания в данный момент ==

== Все нужные ссылки ==

* Телеграм чат команды переводчиков: https://t.me/qgis_translate_ru
* Проект в transifex: https://www.transifex.com/qgis/QGIS/dashboard/
* Сайт QGIS: https://qgis.org
* Дистрибутивы QGIS: https://qgis.org/ru/site/forusers/download.html
* QT Linguist для Windows: https://github.com/lelegard/qtlinguist-installers/releases
* Страница для переводчиков на сайте QGIS: https://qgis.org/ru/site/getinvolved/translate.html

Перевод QGIS на русский язык: сообщество и руководство для участника

2025-12-16T13:35:00Z

Максим Дубинин:

{{Статья|Опубликована|qgis-trans}}
{{Аннотация|Материал посвящен процессу организованного перевода открытой ГИС QGIS на русский язык. Здесь вы найдёте общую информацию о том, как организован перевод, все необходимые ссылки, а также руководство по тому, как стать участником и полезные советы}}

Популярность [https://qgis.org/ru/site/ QGIS] на сегодняшний день впечатляет: эта открытая ГИС используется в индустриальных, научных и образовательных организациях по всему миру. Если вы как-то связаны с миром геоинформационных технологий, то, скорее всего, вам доводилось если не работать с QGIS, то по крайней мере сталкиваться с проектами и инфраструктурами пространственных данных на её основе.

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

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

Перед тем, как приступить к переводу, просим вас полностью прочитать эту статью.

== Как вообще организована локализация QGIS? ==

Перевод на разные языки осуществляется добровольцами (иногда спонсируемыми государственными или частными организациями). Стать таким добровольцем может любой человек, например вы. Для удобства процесса перевода разработчики предлагают использовать специализированное веб-приложение [https://transifex.com transifex.com], которое позволяет редактировать перевод непосредственно в интернет-браузере, без утомительных настроек окружения и установки дополнительных программ. Единственное, что вам нужно - подать заявку на перевод проекта QGIS в transifex, и дождаться её одобрения. Мы предлагаем вам не просто сразу начать перевод, но сделать это более организовано, так, как мы предлагаем в этой статье.

Вы можете познакомиться с [https://qgis.org/ru/site/getinvolved/translate.html описанием процесса локализации] на сайте QGIS.

== Кто организует систематизированный перевод на русский язык? ==

На данный момент перевод координируют:
* Эдуард Казаков (ee.kazakov@gmail.com, t.me/kazakov_e)

Если у вас есть предложения по процессу локализации, свяжитесь с нами по указанным контактам.

== Хочу участовать, что мне понадобится и как подготовиться? ==

Если вы хотите просто сообщить об обнаруженном плохом/несуществующем переводе, но не хотите полноценно участовать в процессе, просто напишите на адрес электронной почты ee.kazakov@gmail.com или в Telegram @kazakov_e / https://t.me/qgis_translate_ru

Для того, чтобы начать работать с нами, вам потребуется:
* Ночная или последняя (не LTR) версия QGIS на компьютере (дистрибутивы и инструкции для всех операционных систем [https://qgis.org/ru/site/forusers/download.html здесь])
* Telegram аккаунт
* Аккаунт в системе [https://transifex.com transifex.com] (регистрация бесплатная)

Если всё это у вас есть, действуйте по следующему алгоритму:
# Подайте заявку на перевод Russian на странице проекта QGIS в transifex: https://www.transifex.com/qgis/QGIS/ (при этом нужно быть авторизованным в transifex).
# Вступите в телеграм-чат https://t.me/qgis_translate_ru и представтесь.
# Напишите в телеграм-чате, что подали заявку в transifex.
# Вас оповестят, когда заявка будет одобрена.
# Переходите на страницу проекта https://www.transifex.com/qgis/QGIS/dashboard/ , выбирайте ветку "qgis-application", и можно приступать к работе!

Сейчас работа организована в формате спринтов, координируемых через Telegram-чат.

== Как получить актуальную версию перевода и увидеть её в QGIS? ==

Во-первых, скачайте и установите последнюю (не LTR) версию QGIS, доступную на [https://qgis.org/ru/site/forusers/download.html# официальном сайте].

Затем на transifex проделайте следующие действия:
# Авторизуйтесь и перейдите на страницу проекта https://www.transifex.com/qgis/QGIS/dashboard/
# В списке языков найдите Russian, наведите на него мышь и нажмите "Обзор ресурсов"
# В появившемся списке нажимайте на qgis-application
# Во всплывшем окне выбирайте пункт "Скачать файл для использования" - загрузится файл формата .ts

[[Файл:Qgis_translate_resource_download_options.png|600px|thumb|center]]

Далее загруженный файл нужно скомпилировать в формат .qm с помощью программ Qt Linguist. Она устанавливается для Win и Mac отсюда: https://github.com/lelegard/qtlinguist-installers/releases, в Linux устанавливается вместе с QGIS. Запустив Linguist, откройте в ней загруженный файл .ts и скомпилируйте через меню File - Release As: достаточно просто указать путь до выходного файла.

Выходной файл нужно назвать qgis_ru.qm и сохранить в папку i18n в ресурсах QGIS, например C:\Program Files\QGIS\apps\qgis\i18n

Перезапустив QGIS, актуальный перевод подхватится.

== Процесс работы переводчика ==

Перед тем, как начать перевод, выберите часть интерфейса, которую хотите перевести. Мы не рекомендуем хаотично браться за перевод отдельно взятых фраз, гораздо продуктивнее целиком переводить раздел за разделом. Уведомите команду в телеграм-чате о том, что вы начали работать над определенной частью интерфейса. Когда начинается спринт, список разделов предлагается организаторами.

Зайдя на страницу проекта в transifex, начинайте перевод ветки "qgis-application": находите фразы, которые хотите перевести или исправить, в строке поиска, и вводите свою версию перевода (некоторые советы по поиску фраз приведены далее в разделе "Практические советы по работе с transifex и важные замечания").

Если вы хотите привлечь внимание проверяющих к определенным строкам (например, они давно хорошо переведены, но не помечены как "проверенные"), то напишите в чат номера таких строк. Номера копируются из текущего адреса в браузере. Например, у вас выделена хорошо переведенная фраза, и в адресной строке вы видите https://www.transifex.com/qgis/QGIS/translate/#ru/qgis-application/52740834, значит нужно скопировать число 52740834.

Строки, которые вы сами изменяли, копировать не нужно - их мы сможем найти по фильтрам с ником и датами. Числа нужны только чтобы пометить хорошие, но не помеченные "проверенными" переводы, которые вы не трогали (не изменяли).

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

== Практические советы по работе с transifex и важные замечания ==

Основное время вы будете проводить в интерфейсе transifex, он предоставляет инструменты по поиску фраз и редактированию/созданию переводов для них. Фраза - это отдельный текст в интерфейсе, он может быть как простейшим (из одного слова, например layer), так и огромным (например, текст помощи для инструмента в панели анализа). Вы можете искать фразы по:
* Их оригинальному написанию (или фрагменту)
* По их существующему переводу (или фрагменту)
* По файлу с кодом, где они записаны
* По контексту (именованной области программного интерфейса QGIS)
* По множеству других метаданных

За поиск отвечает вот эта панель:

[[Файл:Qgis_translate_search_panel.png|600px|thumb|center]]

В выпадающих меню, видных на изображении выше, вы можете управлять, каким фильтром хотите воспользоваться сейчас. При этом имя фильтра попадёт в сам текст поиска, вот, к примеру, как будет выглядеть поиск по переводу:

[[Файл:Qgis_translate_translation_search.png|600px|thumb|center]]

По мере формулирования условий поиска будет формироваться набор фраз, доступных для перевода, в панели прямо под панелью условий поиска:

[[Файл:Qgis_translate_search_result.png|600px|thumb|center]]

Выбрав фразу кликом левой кнопки мыши, в правой части интерфейса вы увидите инструменты работы с фразой, в частности поле дял ввода перевода, а также раздел "подробности", в котором можно увидеть контекст текста и некоторые другие детали. По этим деталям обычно вы можете понять, является ли эта фраза той, которую вы ищете. Например в примере на скриншоте виден контекст AlgorithmClassification - по нему понятно, откуда фраза. Это важно, так как похожие или одинаковые фразы могут встречаться в разных частях интерфейса программы. Введя перевод, нажимайте кнопку "сохранить перевод".

Приступая к переводу конкретного раздела интерфейса, удобно, увидев название контекста на примере одной из фраз, настроить фильтр поиска на этот контекст: так вы будете видеть все связанные с этой частью интерфейса фразы (это работает не всегда - иногда контекст очень абстрактен). Для поиска по контексту вводите в строку поиска текст вида

context:<имя контекста>, например context:AlgorithmClassification

[[Файл:Qgis_translate_context_search.png|600px|thumb|center]]

Среди фраз вы можете встретить такие, которые содержат знак & между буквами - это ссылка на горячую клавишу. Если вы хотите, чтобы в переводном варианте также работала горячая клавиша, то можете поставить & перед нужной буквой.

Также среди фраз возможно встретить такую, в которой используются html-теги (выглядят, как номера на цветных бирках). Эти теги необходимо сохранять в переводной версии фразы. Самый простой способ это сделать, это скопировать оригинальный текст фразы с помощью кнопки "копировать исходную строку", а затем в копии, не трогая теги, перевести только текст.

[[Файл:Qgis_translate_tags.png|400px|thumb|center]]

При возникновении любых вопросов пишите в телеграм-чат команды переводчиков.

== Общие рекомендации ==
* '''Буква ё'''. Мы стараемся использовать букву ё во всех словах русского языка, в которые она входит.

== Перевод списков ключевых слов (тегов) ==

У инструментов панели анализа и редактора выражений есть теги - перечни ключевых слов, по которым их можно быстро искать. Выглядят, например, так:

export,add,information,measurements,areas,lengths,perimeters,
latitudes,longitudes,x,y,z,extract,points,lines,polygons,sinuosity,fields

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

* В поле перевода копируются все исходные англоязычные теги
* К ним через запятую добавляются переводные версии всех тегов. Для примера выше перевод будет следующим:

''export,add,information,measurements,areas,lengths,perimeters,
latitudes,longitudes,x,y,z,extract,points,lines,polygons,sinuosity,
fields,экспорт,добавить,информация,измерение,площадь,длина,периметр,
широта,долгота,точки,линии,полигоны,поля''

== Переводы некоторых спорных ключевых терминов, которые необходимо соблюдать ==

Здесь будет пополняться список '''дискуссионных''' терминов, которые встречаются повсеместно, и нужна единая стратегия для их перевода.

Если вы не согласны, будем рады обсудить в Telegram-чате https://t.me/qgis_translate_ru и выработать общее решение.

* '''Authentication - Аутентификация''' (мотивация: термин '''аутентификация''' закреплен в ГОСТ Р ИСО/МЭК 9594-8-98)
* '''Extent - Охват'''. Отброшенные альтернативы: граница, экстент
* '''Folder - Каталог'''. Не "папка" и не "директория".
* '''Plugin - Модуль'''. Не "Плагин".
* '''px''' (сокращенное название единицы измерения "пиксель") '''- пикс'''.
* '''Script - Сценарий'''. Не "Скрипт".
* '''Symbol - Символ''' / Symbology - символизация (мотивация: более общая терминология. Привычна ГИС-пользователям в РФ, т.к. используется в ArcGIS много лет). Отброшенные альтернативы: знак.
* '''Tile - Тайл''' (мотивация: термин вошёл в полноценный обиход, даже в гос. контрактах значатся термины вида "тайловый кэш", "векторные тайлы" и подобные им). Отброшенные альтернативы: блок мозаики.
* '''Snapping - Прилипание''' (мотивация: так всегда было в QGIS, и вся документация давно так описывает процесс. Привязка больше ассоциируется с привязкой растров/векторов). Отброшенные альтернативы: привязка.
* '''URL - URL''' (мотивация: аббревиатура закреплена в ГОСТ Р 7.0.5-2008. Понятие URL шире чем адрес). Отброшенные альтернативы: адрес.
* '''Editing mode - Режим редактирования''', '''Edits - Изменения'''. Отброшенные альтернативы: режим правок, правки.
* '''Merge - Слияние''', не объединение
* '''Union - Объединение'''
* '''Data provider - провайдер данных''', не источник данных, иначе возникает конфликт с также распространенным понятием data source
* '''nodata - NoData''', не нет данных, НетДанных и прочие вариации.

== Разделы и термины, требующие уборки/внимания в данный момент ==

== Все нужные ссылки ==

* Телеграм чат команды переводчиков: https://t.me/qgis_translate_ru
* Проект в transifex: https://www.transifex.com/qgis/QGIS/dashboard/
* Сайт QGIS: https://qgis.org
* Дистрибутивы QGIS: https://qgis.org/ru/site/forusers/download.html
* QT Linguist для Windows: https://github.com/lelegard/qtlinguist-installers/releases
* Страница для переводчиков на сайте QGIS: https://qgis.org/ru/site/getinvolved/translate.html

Использование языка PERL для автоматизации работы с пространственными данными в среде Windows

2024-02-22T19:53:11Z

Максим Дубинин:

{{Статья|Опубликована|perl}}
{{Аннотация|Как использовать язык PERL для автоматизации рутинных операций с данными с большим количеством примеров.}}

==Зачем это нужно==

Данная статья безусловно не является пособием по языку PERL, однако она может помочь специалистам, активно использующим пространственные данные , в том числе данные дистанционного зондирования в своей работе. Увеличение количества доступных пространственных данных часто ставит перед исследователями задачу автоматизации процедур подготовки данных и их анализа. В большинстве случаев, при работе крупных информационных центров, эти задачи решаются написанием специальных программ и реализуются на платформах UNIX-LINUX. В тоже время, для большинства менее масштабных работ стандартных ГИС приложений и среды Windows обычно бывает достаточно для организации небольших потоков операций. Одним из средств организации таких потоков операций в Win32 является использование языка PERL. Приведенные в статье примеры скриптов PERL могут не являться оптимальными с точки зрения программирования. Однако они позволяют решать поставленные задачи и просты для понимания и изменения пользователями.

==Краткие основы языка PERL==

===О языке PERL===

Язык PERL - Practical Extraction and Report Language (практический язык извлечений и отчетов) – был создан американским программистом Larry Wall для автоматизации работы с текстом и выполнения рутинных задач в операционной системе UNIX. В настоящее время (сентябрь 2006 года) стандартом является 5-ая версия интерпретатора языка, распространяемая свободно и доступная для платформ UNIX, LINUX, WIN32 и MAC OS. Примеры в этой статье используют версию ActivePerl для Windows ([http://www.activestate.com/Products/ActivePerl/?_x=1 скачать дистрибутив ActivePerl]). Основным удобством использования языка PERL в среде Windows является отсутствие необходимости в компилировании программ – Ваши скрипты будут автоматически запускать интерпретатор языка для выполнения, и изменение программы будет не сложнее внесение правок в текстовый файл! При установке интерпретатора проверьте, чтобы он был установлен в директорию по умолчанию - C:\Perl.

Для получения дополнительной информации и справки по языку можно использовать встроенную справочную систему (на английском), а также русскоязычные пособия, доступные в Интернет: Изучаем Perl (http://perl.find-info.ru/perl/016/index.htm), Введение в Perl (http://perl.org.ru/documentation/docs/perl/index.htm) и другие. Следует иметь ввиду, что как и во всяком OpenSource проекте, в PERL существует масса вариантов синтаксиса одних и тех же команд. В этой статье не всегда приводятся самые удобные и короткие выражения – Вы сможете заменить их на более простые при написании своих скриптов.

===Типы данных, загрузка и вывод данных===

В Perl используется три типа данных: переменные, массивы и хеши. В простых скриптах мы будем использовать только переменные и массивы. Имена переменных всегда начинаются с символа $, например: $filename PERL различает тип переменных по контексту данных, и нет необходимости специально декларировать список переменных в начале программы. В любом месте программы вы можете определить новую переменную и присвоить ей значение:

<source lang="perl">
$filename=”mytext.txt”;
$number=20;
$newnumber=$oldnumber;
</source>

Массив представляет собой последовательность переменных любого типа, элементы массива нумеруются последовательно, начиная с 0. Массив обозначается символом @:

<source lang="perl">
@alldata
</source>

Элемент массива определяется его номером:

<source lang="perl">
#первый элемент массива
@alldata[0]
</source>

Вы можете присвоить переменной значение элемента массива:

<source lang="perl">
$data=@alldata[0];
</source>

или определить переменную как размер (число элементов) в массиве:

<source lang="perl">
$number=@alldata;
</source>

Поскольку мы будем рассматривать скрипты, оперирующие с текстами, то нам необходимо представлять себе как прочитать или записать данные в текстовый файл. Встроенной функцией open можно открыть файл для чтения:

<source lang="perl">
open (DRR, "dirlist.txt");
</source>

или для записи:

<source lang="perl">
open (OUT, ">batch.bat");
</source>

Если файл с таким именем существует, он будет стерт и записан заново. После записи данных в файл его необходимо закрыть:

<source lang="perl">
close (OUT);
</source>

===Проверка значений, циклы===

Основные функции языка программирования, которые нам понадобятся для написания скриптов, будут проверка значений и организация циклов процедур. Эти операции в PERL имеют схожий синтаксис:

<source lang="perl">
#фигурная скобка обозначает начало блока операций
условие
{
операции
}
#конец блока операций
</source>

Проверка значений по типу «если да-то, если нет-то» выглядит так:

<source lang="perl">
if ($p1<=$p2)
{
$U=$p1;
}
else
{ $U=$p2;
}
</source>

В данном примере переменной $U будет присвоено минимальное значение из двух вариантов.

Определение цикла операций выглядит следующим образом:

<source lang="perl">
$num=10;
$inc=0;
while ($inc < $num)
{
#набор операций
++$inc;
}
</source>

В данном случае цикл будет выполнен 10 раз, в конце каждого цикла мы увеличиваем значение переменной $inc на единицу - ++$inc.

===Операторы, встроенные функции===

Вот некоторые из операторов языка, которые будут нами использованы в примерах скриптов:
Присвоение значения: «=» или «my … =»

<source lang="perl">
#определение значения переменной
$inc=0;
#чтение содержимого файла, открытого под именем
#DRR в массив @alldir
my @alldir=<DRR>;
</source>

Отношение: < (меньше), > (больше), <= (меньше или равно), == (равно – для чисел), != (не равно – для чисел), eq (равно – для строк), ne (не равно – для строк).

Арифметические операторы: +,-,*,/.

Мы также будем использовать встроенные функции, отвечающие за ввод и вывод данных и выполняющие простые операции над текстом:

open/close: открытие текстового файла для чтения или записи – open, и закрытие файла – close (синтаксис этих функция рассмотрен выше).

print: вывод на печать.

<source lang="perl">
#значение переменной $line будет выведено на экран
print "$line";
#значение будет напечатано в текстовый файл, #открытый под именем OUT
print OUT "$line";
</source>

Следует отметить наиболее распространенные спецсимволы, которые можно выводить на печать:
\n – конец строки,
\t – символ табуляции,
\ - символ пробела (пробел стоит после \),
\символ - любой символ, если он зарезервирован во внутренней структуре PERL (например, для вывода символов “ ; \ необходимо набрать \“ \; \\ соответственно, иначе программа не будет работать!).
split: разделяет строку по заданному символу-разделителю:

<source lang="perl">
($filenam,$pnum,$zone)=split(',',$fileout);
</source>

Например: если $fileout = «image1,s_11,23» (строка таблицы с информацией о имени файла, номере точки и зоны проекции), то переменные будут иметь следующие значения:

<source lang="perl">
$filenam=«image1»
$pnum=«s_11»
$zone=«23» (последняя переменная будет так же содержать символ переноса строки).
</source>

s///: оператор замены. В общем случае синтаксис оператора выгладит так:

переменная =~ s/символ-шаблон/символ для замены шаблона/;

Например, полученная выше переменная $zone содержит, кроме необходимой информации, также символ переноса строки. Что бы очистить его, используем замену:

$zone =~ s/\n//; - переменная $zone равна переменной $zone, в которой символ конца строки \n заменен на отсутствие символа.

system: выполняет системную функцию. Например, для получения списка файлов в директории, используем функцию MS-DOS dir:

<source lang="perl">
$com="C:\\Geomatica_V91\\exe\\easi.exe r script";
system($com);
</source>

PERL запустит интерпретатор MS-DOS и выполнит указанную в переменной строку.

===Общие правила написания скриптов в PERL и их отладка===

Скрипт PERL в среде Windows представляет собой простой текстовый файл с расширением .pl. Если в Вашей системе установлен интерпретатор языка, то файлы *.pl будут иметь значок с желтым шариком, и будут выполняться автоматически по двойному щелчку на имя файла (подобно файлам .bat). Первая строка файла должна содержать путь к интерпретатору языка – в нашем случае строка эта будет выглядеть так:

<source lang="perl">
#!/usr/local/bin/perl
</source>

Далее вводиться текст программы. Все строки, в которых используются операторы или функции, должны заканчиваться точкой с запятой - ;. Исключение составляют строки, определяющие условия или операторы циклов:

<source lang="perl">
while ($inc < $num)
{
print "привет";
++$inc;
}
</source>

Примечания обозначаются символом # - весь текст от символа до конца строки будет считаться примечанием.

Никаких специальных символов конца скрипта в PERL не используется.
Для запуска скриптов, рассчитанных на работу со всеми файлами или поддиректориями в определенной директории, требуется скопировать файл скрипта в эту папку и запустить его. В противном случае Вам придется использовать не относительные, а абсолютные пути к папкам и файлам. Для отладки скриптов полезно бывает использовать интерпретатор MS-DOS. Для этого запустите Command Prompt, перейдите в директорию со скриптами, и наберите имя файла скрипта с расширением. Файл будет запущен, и информация об ошибках появиться на экране.

==Примеры скриптов==

Внимание: для работы с примерами, которые можно загрузить с данной страницы, необходимо сменить их расширение с *.txt на *.pl.

===Автоматизация работы с файлами: Операции с файлами и директориями===

Приведем для примера два простых скрипта. Первый скрипт создает поддиректории во всех поддиректориях рабочей директории (рабочей мы будем называть директорию, из которой запущен скрипт) ([http://gis-lab.info/programs/perl/newfolder.txt скачать скрипт]):

<source lang="perl">
#Стандартное начало скрипта:
#!/usr/local/bin/perl
#Определяем строковую переменную, содержащую команду DOS для вывода списка поддиректорий в файл dirlist.txt (используем ключи команды dir: /A:D – только директории, /O:N – сортировка по имени и /B – краткий формат) :
$com="dir /A:D /O:N /B >dirlist.txt";
#Выполняем команду в операционной системе:
system($com);
#Открываем для чтения файл со списком поддиректорий:
open (DRR, "dirlist.txt");
#Читаем содержимое файла в массив @alldir:
my @alldir=<DRR>;
#Определяем длину массива:
$num=@alldir;
#Указываем номер первого элемента (напомним, что элементы в массиве нумеруются с нуля, и последний элемент будет иметь номер на единицу меньший числа элементов):
$inc=0;
#Запускаем цикл операций, который будет продолжаться до тех пор, пока значение переменной $inc (номер элемента массива) меньше переменной $num (число элементов в массиве). В ходе этого цикла мы прочитаем последовательно все строки файла, цикл закончиться после прочтения последней строки:
while ($inc < $num)
#Начало операций цикла:
{
#Присваиваем переменной $line значение первой строки массива (то есть первой строки текстового файла):
$line = @alldir[$inc];
#Удаляем ненужный символ конца строки:
$line =~ s/\n//;
#Определяем команду для создания поддиректории: команда DOS md (makedir), имя поддиректории составляем из имени директории, сохраненного в переменной $line и нужного повторяющегося имени, для вывода на печать символы \ и “ используем вспомогательный символ \, показывающий, что стоящий за ним символ не является командой (поскольку эти символы зарезервированы PERL). Имя директории указано в кавычках чтобы избежать проблем с пробелами в именах директорий:
$com="md \"$line\\data\"\n";
#Выполняем команду в операционной системе:
system($com);
#Увеличиваем значение переменной $inc на единицу:
++$inc;
#Конец операций цикла:
}
#После выполнения команд удалим файл dirlist.txt. Для этого сначала закроем его:
close (DRR);
#Затем определим команду DOS:
$com="del dirlist.txt\n";
#И выполним ее:
system($com);
</source>

Примечание – хотя данный скрипт может использовать директории с пробелами, многие приведенные ниже скрипты при наличии пробелов в именах директорий и файлов не будут работать. Впрочем, тоже относиться и к большинству распространенных ГИС-приложений.

Второй скрипт скопирует файлы с именем time1.tif из всех поддиректорий рабочей директории в директорию output, заменив имя time1 на имя поддиректории, из которой он был скопирован (то есть файл time1.tif из директории 123_056 получит имя 123_056.tif и так далее). Мы могли бы использовать такой же подход, как в предыдущем случае, но сознательно применили другой путь – теперь скрипт PERL не будет самостоятельно выполнять команды ОС, а создаст исполняемый файл, который можно будет запустить позже. Хотя такой подход менее продуктивен, он позволяет избежать необходимости восстанавливать данные при ошибках в программе – понятно, что исполняемый batch файл легче проверить и убедиться, что он будет работать правильно ([http://gis-lab.info/programs/perl/copyfile.txt скачать скрипт]).

<source lang="perl">
#!/usr/local/bin/perl
#Создаем список директорий и читаем его:
$com="dir /A:D /O:N /B >dirlist.txt";
system ($com);
open (DRR, "dirlist.txt");
#Открываем файл для записи (если файл batch.bat существует, он будет стерт и перезаписан):
open (OUT, ">batch.bat");
#Печатаем в batch команду создания директории output:
print OUT "md output\n";
#Запускаем цикл операций:
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
#Читаем имена директорий:
$line = @alldir[$inc];
$line =~ s/\n//;
#Печатаем команду копирования и переименования файла:
print OUT "copy \"$line\\time1.tif\" \"output\\$line.tif\"\n";
++$inc;
}
#Закрываем файлы:
close (OUT);
close (DRR);
#Удаляем файл dirlist.txt:
$com="del dirlist.txt\n";
system($com);
</source>

===Автоматизация работы с файлами: Распаковка архивов===

Рассмотрим примеры, более близкие к практике. Скрипт, приведенный ниже, может помочь в автоматической распаковке снимков Landsat-7, скаченных с сайта GLCF (http://glcfapp.umiacs.umd.edu:8080/esdi/index.jsp). Мы скачали целиком директории со снимками и теперь хотим удалить все технические файлы и ненужные каналы, и распаковать файлы tif из архивов GZ с помощью архиватора WinRAR (установленного в C:\WinRAR) ([http://gis-lab.info/programs/perl/landsat_unpack_rar.txt скачать скрипт]):

<source lang="perl">
#!/usr/local/bin/perl
#Создаем список директорий и читаем его:
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
open (DRR, "dirlist.txt");
#Открываем файл для записи:
open(OUT, ">batch.bat");
#Запускаем цикл операций:
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
#Читаем имена директорий:
$line = @alldir[$inc];
$line =~ s/\n//;
#Печатаем команду перехода в поддиректорию:
print OUT "cd $line\n";
#Печатаем команды удаления ненужных файлов (если файлы отсутствуют, команды будут пропущены):
print OUT "del /Q *_nn6*\n";
print OUT "del /Q *_nn8*\n";
print OUT "del /Q *_nn1*\n";
print OUT "del /Q *.jpg\n";
print OUT "del /Q *.met\n";
print OUT "del /Q *.hdr\n";
print OUT "del /Q *.i*\n";
#Печатаем команду разархивирования:
print OUT "C:\\WinRAR\\winrar.exe e -ibck *.gz\n";
#Печатаем команду удаления архивов:
print OUT "del /Q *.gz\n";
#Печатаем команду возврата в рабочую директорию:
print OUT "cd ..\n";
++$inc;
}
close (DRR);
close (OUT);
$com="del dirlist.txt";
system($com);
</source>

===Автоматизация работы с файлами: Упаковка файлов===

Если вы получили снимки Landsat из архивов EROS Data Center, то, скорее всего, они будут в формате tif без сжатия. Мы хотели бы сжать эти файлы в архивы формата GZ с помощью архиватора GZip (установленного в C:\Gzip), причем каждый файл tif в отдельный архивный файл. Для этого мы используем [http://gis-lab.info/programs/perl/landsat_pack_gz.txt вот такой скрипт]:

<source lang="perl">
#!/usr/local/bin/perl
#Создаем список директорий и читаем его:
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
open (DRR, "dirlist.txt");
#Запускаем цикл операций:
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
#Читаем имена директорий:
$line = @alldir[$inc];
$line =~ s/\n//;
#Создаем список файлов в поддиректории и читаем его:
$com="dir $line\\*.tif /B >dirlist$inc.txt\n";
system($com);
$newname="dirlist$inc.txt";
open (FRR, $newname);
#Запускаем второй цикл операций внутри первого цикла:
my @files=<FRR>;
$num1=@files;
$inc1=0;
while ($inc1 < $num1)
{
#Читаем имена файлов:
$line1 = @files[$inc1];
$line1 =~ s/.tif\n//;
#Определяем и выполняем команду упаковки файлов:
$com = "C:\\Gzip\\gzip $line\\$line1.tif\n";
system($com);
++$inc1;
#Завершаем второй цикл:
}
++$inc;
Завершаем первый цикл:
}
close (DRR);
close (FRR);
#Удаляем все временные текстовые файлы:
$com="del dirlist*.txt\n";
system($com);
</source>

===Автоматизация работы в PCI Geomatica: Запуск скриптов EASI===

Одним из наиболее удобных приложений PCI Geomatica является язык EASI – простой язык для написания скриптов, позволяющих выполнить большинство распространенных операций подготовки и анализа данных ДЗ. Скрипты EASI имеют расширение .eas и могут быть запущены из командной строки DOS – имя скрипта указывается в качестве атрибута при вызове программы. Простейший скрипт PERL, приведенный ниже, позволяет запускать один и тот же скрипт EASI последовательно в каждой поддиректории рабочей директории – то есть позволяет выполнить однотипные операции над всеми наборами данных, используемых в Вашем проекте. Для удобства использования все необходимые скрипты EASI заранее собраны в директории C:\scripts\, так что перед выполнением скрипта необходимо только указать имя нужного скрипта для исполнения ([http://gis-lab.info/programs/perl/execute_easi.txt скачать скрипт]):

<source lang="perl">
#!/usr/local/bin/perl
$com1="dir /A:D /O:N /B >dirlist.txt";
system($com1);
open (DRR, "dirlist.txt");
open(OUT, ">batch.bat");
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
$line = @alldir[$inc];
$line =~ s/\n//;
#Печатаем комманды в исполняемый файл DOS: переходим в поддиректорию...
print OUT "cd $line\n";
#...копируем скрипт, переименовывая его в script.eas...
print OUT "copy C:\\scripts\\myprocess.eas script.eas /y\n";
#...выполняем скрипт...
print OUT "C:\\Geomatica_V91\\exe\\easi.exe r script\n\n";
#...и удаляем скрипт и файл параметров EASI:
print OUT "del script.eas\n\n";
print OUT "del prm.prm\n\n";
#Возвращаемся в рабочую директорию:
print OUT "cd ..\n";
++$inc;
}
close (OUT);
close (DRR);
$com="del dirlist.txt";
system($com);
</source>

===Автоматизация работы в PCI Geomatica: Генерирование скриптов EASI===

Простота языка EASI делает возможным автоматическое создание скриптов с помощью программы PERL по заданному шаблону. Поскольку задачи, для которых может потребоваться подобная процедура, достаточно индивидуальны, мы не будем приводить примера готового скрипта, а покажем только часть кода. В данной задаче исходными данными служили файлы с таблицами координат углов прямоугольника, спроецированного из синусоидальной проекции в UTM. Требовалось построить пустую растровую базу данных PCIDSK с размером ячейки 28.5 метров таким образом, что бы она полностью включала в себя полученный прямоугольник. Данный кусок кода представляет собой операцию внутри цикла:

<source lang="perl">
#Вне цикла мы отрыли список файлов (в формате: «имя файла, номер полигона, номер зоны UTM») и прочитали его в массив @filelist. В ходе цикла мы открываем строку файла списка, читаем ее, выделяем имя файла и открываем этот файл:
$fileout = @filelist[$inc];
($filenam,$pnum,$zone)=split(',',$fileout);
$zone =~ s/\n//;
open (OUT, "$filenam");
#Читаем содержимое файла – четыре строки с парами координат:
my @filetext = <OUT>;
$line1 = @filetext[0];
$line2 = @filetext[1];
$line3 = @filetext[2];
$line4 = @filetext[3];
#Разбиваем строки с получением значений координат:
($p1_x,$p1_y)=split(' ',$line1);
($p2_x,$p2_y)=split(' ',$line2);
($p3_x,$p3_y)=split(' ',$line3);
($p4_x,$p4_y)=split(' ',$line4);
#Сравниваем координаты для нахождения углов прямоугольника, полностью включающего в себя исходный прямоугольник (extent):
if ($p1_x<=$p3_x)
{
$UL_x=$p1_x;
}
else
{
$UL_x=$p3_x;
}
if ($p3_y<=$p2_y)
{
$UL_y=$p2_y;
}
else
{
$UL_y=$p3_y;
}
if ($p2_x<=$p4_x)
{
$LR_x=$p4_x;
}
else
{
$LR_x=$p2_x;
}
if ($p1_y<=$p4_y)
{
$LR_y=$p1_y;
}
else
{
$LR_y=$p4_y;
}
#Создаем новую директорию с использованием функции PERL:
mkdir "s_$pnum";
#Расширяем размер создаваемой базы данных на 100 метров для избегания краевых эффектов при последующих операциях:
$ULx2=$UL_x-100;
$ULy2=$UL_y+100;
$LRx2=$LR_x+100;
$LRy2=$LR_y-100;
#Определяем размер создаваемой БД в пикселях:
$sizeX2=int(($LRx2-$ULx2)/28.5);
$sizeY2=int(($ULy2-$LRy2)/28.5);
#Рассчитываем новые координаты правого нижнего угола исходя из размера пикселя:
$LRx2n=$ULx2+($sizeX2*28.5);
$LRy2n=$ULy2-($sizeY2*28.5);
#Печатаем команды создания базы данных (файл с именем EASI заранее был открыт для записи):
print EASI "FILE=\"s_$pnum\\clip\n";
print EASI "DBSZ=$sizeX2,$sizeY2\n";
print EASI "DBNC=15,0,0,0\n";
print EASI "DBLAYOUT=\"FILE\n";
print EASI "RUN CIM\n";
print EASI "\n";
#Печатаем команды присвоения новой БД параметров проекции. Обратите внимание на использование переменной $zone, которая была прочитана в самом начале цикла.
print EASI "UPLEFT=$ULx2,$ULy2\n";
print EASI "LORIGHT=$LRx2n,$LRy2n\n";
print EASI "MAPUNITS=\"UTM $zone E012\n";
print EASI "RUN GEOSET\n";
print EASI "\n";
</source>

Примечание – для включения скриптов EASI в непрерывный процесс Вы можете запускать созданные скрипты автоматически из PERL, используя команду:

<source lang="perl">
$com="C:\\Geomatica_V91\\exe\\easi.exe r script";
system($com);
</source>

После чего использовать результаты выполнения скрипта для дальнейших действий.

===Автоматизация работы в ERDAS Imagine: Генерирование пакетов команд===

Скрипт, приведенный в первом примере, создаст пакетный файл команд ERDAS Imagine для относительно простой процедуры Layer merge. В параметрах процедуры требуется указать абсолютные имена файлов. Для этого нам необходимо знать или задать имя рабочей директории. Хотя определить имя директории, из которой запущен скрипт, можно автоматически, мы, в данном случае, использовали переменную, значение которой необходимо изменить на абсолютный путь к рабочей директории перед запуском скрипта. Так же обратите внимание на использование в ERDAS путей в стандарте UNIX. Данный скрипт рассчитан на определенный шаблон имени файлов (стандарт Ortho Landsat для GLCF): p169r057_7t20000127_z37_nn30.tif, и будет объединять только 3,4 и 5 каналы снимков ([http://gis-lab.info/programs/perl/layermerge_erdas.txt скачать скрипт]):

<source lang="perl">
#!/usr/local/bin/perl
#Абсолютный путь к рабочей директории (обратите внимание на
#использование стандарта UNIX в путях!):
$path="D:/glcf";
#Создание списка папок и запуск цикла:
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
open (DRR, "dirlist.txt");
open(OUT, ">merge.bcf");
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
$line = @alldir[$inc];
$line =~ s/\n//;
#Печатаем в файл имя одного из файлов tif из папки для получения шаблона
#имени:
$com="dir $line\\*nn20.tif /B >dirlist$inc.txt\n";
system($com);
$newname="dirlist$inc.txt";
open (FRR, $newname);
my @files=<FRR>;
$line1 = @files[0];
#Разбиваем имя – получаем позицию сцены, дату и номер зоны UTM:
($wrs,$date,$temp1,$temp2)=split('\_',$line1);
#Печатаем длинную строку команды layermerge в пакетный файл ERDAS bcf.
#Обратите внимание, что это одна строка, переносы здесь недопустимы!
print OUT "modeler -nq \$IMAGINE_HOME/etc/models/layermerge.pmdl -meter -state Union None \'FLOAT RASTER n1 FILE OLD NEAREST NEIGHBOR AOI NONE EDGE FILL\"$path/$line/$wrs\_$date\_$temp1\_nn50.tif\"\; FLOAT RASTER n2 FILE OLD NEAREST NEIGHBOR AOI NONE EDGE FILL\"$path/$line/$wrs\_$date\_$temp1\_nn40.tif\"\; FLOAT RASTER n3 FILE OLD NEAREST NEIGHBOR AOI NONE EDGE FILL\"$path/$line/$wrs\_$date\_$temp1\_nn30.tif\"\;' ignore Unsigned_8_bit \"$path/$wrs\_$date\_utm.img\" \'n1(1)\,n2(1)\,n3(1)\'\n";
++$inc;
}
close (OUT);
close (DRR);
$com="del dirlist*.txt";
system($com);
</source>

Второй пример содержит скрипт для расчета среднего и стандартного отклонения в плавающем окне для всех файлов stat.tif из всех поддиректорий рабочей директории. Его основным отличием от предыдущего скрипта является необходимость внесения путей к файлам в отдельные переменные. Смена абсолютного пути к рабочей директории здесь не продумана, так что при необходимости изменения, придется изменить все четыре пути к файлам (более правильным было бы введение дополнительной переменной, как в предыдущем варианте).

<source lang="perl">
#!/usr/local/bin/perl
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
open (DRR, "dirlist.txt");
open(OUT, ">batch.bcf");
my @alldir=<DRR>;
$num=@alldir;
#Поскольку мы будем создавать два набора переменных для каждого файла,
#увеличиваем число шагов вдвое:
$num=$num*2;
$inc=0;
while ($inc < $num)
{
$n=$inc/2;
$line = @alldir[$n];
$line =~ s/\n//;
#Печатаем первую команду – вычисление среднего:
print OUT "variable Inp$inc Auto \"d:/data/$line/stat.tif\"\;\n";
print OUT "variable Out$inc Auto \"d:/data/$line/statmean.img\" Delete_Before\;\n";
<nowiki> print OUT "modeler -nq \$IMAGINE_HOME/etc/models/focalanal.pmdl -meter -state \"\$(Inp$inc)\" Integer \"\$(Out$inc)\" 'Unsigned 8 bit' Float 'matrix (3,3: 1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000)' Mean '\$(Inp$inc.Ulx)' '\$(Inp$inc.Uly)' '\$(Inp$inc.Lrx)' '\$(Inp$inc.Lry)' Map useall \"None\" '' '' '' ''\n";</nowiki>
++$inc;
#Вторая команда – вычисление стандартного отклонения:
print OUT "variable Inp$inc Auto \"d:/data/$line/stat.tif\"\;\n";
print OUT "variable Out$inc Auto \"d:/data/$line/statstd.img\" Delete_Before\;\n";
<nowiki> print OUT "modeler -nq \$IMAGINE_HOME/etc/models/focalanal.pmdl -meter -state \"\$(Inp$inc)\" Integer \"\$(Out$inc)\" 'Float Single' Float 'matrix (3,3: 1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000,1.000)' SD '\$(Inp$inc.Ulx)' '\$(Inp$inc.Uly)' '\$(Inp$inc.Lrx)' '\$(Inp$inc.Lry)' Map useall \"None\" '' '' '' ''\n";</nowiki>
++$inc;
}
close (OUT);
close (DRR);
$com="del dirlist.txt";
system($com);
</source>

===Автоматизация работы в ERDAS Imagine: Генерирование моделей ===

Однако иногда встроенных команд ERDAS бывает недостаточно и возникает необходимость запуска собственных моделей. В таком случае мы использовали следующий подход: сначала создается и тестируется модель, потом она записывается в файл в виде скрипта, затем команды скрипта переносятся в PERL в качестве атрибутов функции print. В таком случае, для каждого набора файлов PERL будет создавать собственную модель, используя абсолютные пути к файлам, а затем все созданные модели будут включены в пакетный файл для автоматического запуска. Ниже приведен пример скрипта для создания и выполнения моделей для расчета стандартного отклонения для 5-канальных файлов stat.tif из всех поддиректорий ([http://gis-lab.info/programs/perl/model_erdas.txt скачать скрипт]):

<source lang="perl">
#!/usr/local/bin/perl
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
open (DRR, "dirlist.txt");
#Открываем для записи пакетный файл:
open(OUT, ">batch.bcf");
my @alldir=<DRR>;
$num=@alldir;
$inc=0;
while ($inc < $num)
{
$line = @alldir[$inc];
$line =~ s/\n//;
#Создаем и открываем файл модели:
open(MOD, ">std$line.mdl");
#Печатаем команды модели для данного набора файлов (нет переносов!):
print MOD "COMMENT \"no comment\"\;\n";
print MOD "SET CELLSIZE MIN\;\n";
print MOD "SET WINDOW UNION\;\n";
print MOD "SET AOI NONE\;\n";
print MOD "Integer RASTER n1_temp FILE OLD NEAREST NEIGHBOR AOI NONE \"D\:/data/$line/stat.tif\"\;\n";
print MOD "Integer RASTER n11_std FILE NEW USEALL ATHEMATIC 8 BIT UNSIGNED INTEGER \"D\:/data/$line/statstd1.img\"\;\n";
print MOD "Integer RASTER n12_std FILE NEW USEALL ATHEMATIC 8 BIT UNSIGNED INTEGER \"D\:/data/$line/statstd2.img\"\;\n";
print MOD "Integer RASTER n13_std FILE NEW USEALL ATHEMATIC 8 BIT UNSIGNED INTEGER \"D\:/data/$line/statstd3.img\"\;\n";
print MOD "Integer RASTER n14_std FILE NEW USEALL ATHEMATIC 8 BIT UNSIGNED INTEGER \"D\:/data/$line/statstd4.img\"\;\n";
print MOD "Integer RASTER n15_std FILE NEW USEALL ATHEMATIC 8 BIT UNSIGNED INTEGER \"D\:/data/$line/statstd5.img\"\;\n";
print MOD "INTEGER MATRIX n6_Low_Pass\;\n";
print MOD "n6_Low_Pass = MATRIX(3, 3\:\n";
print MOD "\t1, 1, 1,\n";
print MOD "\t1, 1, 1, \n";
print MOD "\t1, 1, 1)\;\n";
print MOD "n11_std = FOCAL STANDARD DEVIATION ( \$n1_temp(1) , \$n6_Low_Pass ) \;\n";
print MOD "n12_std = FOCAL STANDARD DEVIATION ( \$n1_temp(2) , \$n6_Low_Pass ) \;\n";
print MOD "n13_std = FOCAL STANDARD DEVIATION ( \$n1_temp(3) , \$n6_Low_Pass ) \;\n";
print MOD "n14_std = FOCAL STANDARD DEVIATION ( \$n1_temp(4) , \$n6_Low_Pass ) \;\n";
print MOD "n15_std = FOCAL STANDARD DEVIATION ( \$n1_temp(5) , \$n6_Low_Pass ) \;\n";
print MOD "QUIT\;\n";
#Закрываем файл модели:
close (MOD);
#Печатаем команду выполнения модели в пакетный файл:
print OUT "modeler -nq D:/data/std$line.mdl -meter -state\n";
++$inc;
}
close (OUT);
close (DRR);
</source>

<nowiki>*Примечание – данный скрипт является примером организации работы с пользовательскими моделями в ERDAS Imagine. Возможны более продуктивные подходы – например, написание моделей с пользовательскими переменными и пр.</nowiki>

===Автоматизация работы в ArcInfo: Генерирование скриптов AML===

Разумеется, PERL имеет огромный потенциал для автоматического создания скриптов AML. К сожалению, автор практически не использовал в своей работе ArcInfo. В качестве примера можно использовать [http://gis-lab.info/programs/perl/aml_example.txt небольшой скрипт PERL] для создания огромного файла для выполнения команды Generate.

===Автоматизация импорта данных MODIS: Использование MRT===

Использование пакетных файлов для запуска MRT уже реализовано в виде, значительно более дружелюбном к пользователю, чем использование скриптов PERL – например, существует специальная программа, [modisimport-tool.html представленная на сайте GIS-Lab]. Тем не менее, мы приведем здесь один из вариантов скриптов, который может быть модифицирован и является общим (хотя и не оптимальным) решением проблемы автоматического импорта, склеивания и перепроецирования продуктов MODIS. Для использования скрипта Вам потребуется изменить путь к программе MRT на тот, который Вы использовали при установке, и разместить свои данные в директории C:\modwork\input. В случае данного скрипта используются 6 ячеек сетки для каждого композита – если число ячеек другое, его надо изменить, также надо изменить параметры проекции и набор каналов. Программа создаст список файлов HDF в Вашей директории, создаст папки для мозаик (C:\modwork\mosaic) и для выгрузки данных (C:\modwork\output), файлы с параметрами перепроецирования для каждого набора мозаик и исполняемый файл для запуска мозаики и перепроецирования ([http://gis-lab.info/programs/perl/mrt_convert.txt скачать скрипт]).

===Автоматизация импорта данных MODIS: Использование MRT-Swath===

Более важная для нас операция – автоматический импорт и перепроецирование продуктов MOD02QKM (Swath продуктов MODIS) с помощью программы MRT-Swath еще не решен в виде отдельного приложения, и поэтому использование скрипта PERL здесь оправдано. Для запуска скрипта надо изменить (или переместить данные):
#Имя папки с данными – по умолчанию файлы MOD02QKM находятся в D:\MODISWORK\MOD02QKM, а геолокационные поля к ним (MOD03) – в папке D:\\MODISWORK\MOD03\.
#Путь к программе – в данной версии это C:\MRTSwath\.
#Параметры проекции и набор каналов. Мы предполагаем, что для всех сцен скачены оба продукта (MOD02QKM и MOD03), и их наборы полностью соотвествуют друг другу. Скрипт создаст папку для выгрузки результатов - D:\MODISWORK\output\, исполняемый файл batch.bat для запуска перепроецирования и набор файлов параметров для каждого снимка. Запускать скрипт надо из директории D:\MODISWORK\ (или той, на которую Вы измените это название!):

<source lang="perl">
#!/usr/local/bin/perl
#Создаем список продуктов MOD02QKM:
$com="dir MOD02QKM\\*.hdf /B >dir1.txt";
system($com);
#Создаем список продуктов MOD03 (мы предполагаем, что для всех снимков
#скачены геолокационные поля):
$com="dir MOD03\\*.hdf /B >dir2.txt";
system($com);
#Создаем директорию для выгрузки данных
$com="md output";
system($com);
#Читаем файлы и создаем исполняемый файл:
open (DR1, "dir1.txt");
open (DR2, "dir2.txt");
open(OUT, ">batch.bat");
my @allmodis=<DR1>;
my @allgeo=<DR2>;
#Запускаем цикл для всех снимков:
$num=@allmodis;
$inc=0;
while ($inc < $num)
{
#Читаем имя файла
$line1 = @allmodis[$inc];
$line2 = @allgeo[$inc];
#Разбиваем имя по шаблону и создаем имя выходного файла:
($temp1,$temp2,$temp3,$temp4,$temp5,$temp6)=split('\.',$line1);
$newfile="$temp1\_$temp2\_$temp3.tif";
#Создаем файл параметров:
open(PRM, ">mod$inc.prm");
#Печатаем имена файлов для импорта:
print PRM "INPUT_FILENAME = D:\\MODISWORK\\MOD02QKM\\$line1\n";
print PRM "GEOLOCATION_FILENAME = D:\\MODISWORK\\MOD03\\$line2\n";
#Печатаем список каналов:
print PRM "INPUT_SDS_NAME = EV_250_RefSB, 0, 1\n";
#Печатаем имя выходного файла:
print PRM "OUTPUT_FILENAME = D:\\MODISWORK\\output\\$newfile\n";
#Печатаем формат и параметры проекции выходного файла:
print PRM "OUTPUT_FILE_FORMAT = GEOTIFF_FMT\n";
print PRM "KERNEL_TYPE (CC/BI/NN) = NN\n";
print PRM "OUTPUT_PROJECTION_NUMBER = ALBERS\n";
print PRM "OUTPUT_PROJECTION_PARAMETER = 0.0 0.0 52.0 64.0 45.0 0.0 8500000.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0\n";
print PRM "OUTPUT_PROJECTION_SPHERE = 15\n";
print PRM "OUTPUT_PIXEL_SIZE = 250\n";
#Закрываем файл параметров:
close (PRM);
#Печатаем в командный файл команду импорта и перепроецирования:
print OUT "C:\\MRTSwath\\bin\\swath2grid -pf=mod$inc.prm\n";
++$inc;
}
close (OUT);
close (DRR);
close (FRR);
$com="del *.txt\n";
system($com);
</source>

===Автоматизация обработки текстовых отчетов===

Обработка текстов и генерирование отчетов является основной задачей PERL, однако при использовании в ГИС такие задачи не являются первостепенными. Тем не менее, при получении большого количества результатов, для их сравнения иногда необходимо применять автоматизацию. Например, после классификации большого количества снимков требуется создание суммарной таблицы соотношения классов. В нашем случае каждый снимок классифицировался на классы «лес» и «не лес», кроме того, ряд пикселей был отмечен как «нет данных». Таким образом, результирующий растровый слой изображения содержит три класса данных:
1 – нет данных
2 – лес
3 - не лес
Мы можем выгрузить гистограмму значений этого слоя в текстовый файл (например, командой EASI histdump) – таким образом, в каждой директории с данными мы создадим файл blockforest.txt, который представляет собой колонку чисел:
10 – число пикселей первого класса
128779 – число пикселей второго класса и т.д....
6576
0
0
...
Создадим скрипт для автоматического сбора данных из всех файлов blockforest.txt и создания таблицы:

<source lang="perl">
#!/usr/local/bin/perl
#Создадим список поддиректорий для обработки:
$com="dir /A:D /O:N /B >dirlist.txt";
system($com);
#Откроем список директорий и создадим файл с результирующей таблицей:
open (DRR, "dirlist.txt");
open(OUT, ">output.txt");
#В результирующую таблицу добавим «шапку»:
print OUT "name\t nodata\t forest\t noforest\n";
my @alldir=<DRR>;
#Запустим цикл:
$num=@alldir;
$inc=0;
while ($inc < $num)
{
#Читаем имя директории:
$line = @alldir[$inc];
$line =~ s/\n//;
#Открываем текстовый фал из нужной директории:
$fileforest="$line\\blockforest.txt";
open (bforest, "$fileforest");
#Читаем значения гистограммы:
my @allforest=<bforest>;
$bforest1 = @allforest[0];
$bforest1 =~ s/\n//;
$bforest2 = @allforest[1];
$bforest2 =~ s/\n//;
$bforest3 = @allforest[2];
$bforest3 =~ s/\n//;
#Печатаем значения переменных в результирующую таблицу:
print OUT "$line\t$bforest1\t$bforest2\t$bforest3\n";
++$inc;
}
close (OUT);
close (DRR);
$com="del dirlist.txt";
system($com);
</source>

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

==Ссылки по теме==
*[Perl.org]

Примеры использования инструментов GDAL

2023-11-04T13:42:50Z

Максим Дубинин: /* Создание растров уменьшенного разрешения (т.н. quicklook, preview) */

{{Статья|Опубликована|gdal-examples}}
{{Аннотация|Перечь примеров для справки}}

[http://www.gdal.org/index.html GDAL/OGR] - библиотека для работы с географическими форматами данных. GDAL представляет собой набор утилит для обработки растровых данных, в то время, как OGR предназначена для работы с векторными форматами. В статье рассматриваются некоторые практические примеры применения утилит этой библиотеки для работы с растровыми данными.

С библиотекой GDAL так же поставляется утилита ogr2ogr, предназначенная для работы с векторными данными. Примеры использования этой утилиты приводятся в [http://gis-lab.info/qa/ogr2ogr-examples.html другой статье].

Если у вас есть свои часто используемые примеры - присылайте автору или дописывайте прямо здесь.

__TOC__

== Общие сведения ==

Утилиты GDAL предназначены для конвертации растровых данных из одного формата в другой и выполнения над ними различных операций. Установить утилиты GDAL для Windows можно с [http://gis-lab.info/qa/qgis-osgeo4w.html помощью OSGeo4W]. Удобный визуальный интерфейс для утилит имеется в [http://gis-lab.info/qa/gdaltools.html GDALTools для QGIS]. В комплект GDAL входят следующие утилиты:

*[http://www.gdal.org/gdal_calc.html gdal_calc.py] - растровый калькулятор, арифметические операции с растрами;
*[http://www.gdal.org/gdal-config.html gdal-config] - получить опции необходимые для создания ПО использующего GDAL;
*[http://www.gdal.org/gdal_contour.html gdal_contour] - получение изолиний по цифровым моделям рельефа (ЦМР);
*[http://www.gdal.org/gdal_edit.html gdal_edit.py ] - редактировать информацию о растре;
*[http://www.gdal.org/gdal_fillnodata.html gdal_fillnodata.py] - заполнение областей имеющих значение NODATA;
*[http://www.gdal.org/gdal_merge.html gdal_merge.py] - создание мозаик и композитных изображений;
*[http://www.gdal.org/gdal_polygonize.html gdal_polygonize.py] - векторизовать растр с получением полигонального слоя;
*[http://www.gdal.org/gdal_proximity.html gdal_proximity.py] - расчитать растр близости;
*[http://www.gdal.org/gdal_rasterize.html gdal_rasterize] - растеризация векторных данных;
*[http://www.gdal.org/gdal_retile.html gdal_retile.py] - создать новый набор тайлов и/или перестроить пирамидные слои;
*[http://www.gdal.org/gdal_grid.html gdal_grid] - создание ЦМР из векторных данных;
*[http://www.gdal.org/gdal_sieve.html gdal_sieve.py] - фильтрация осколочных объектов растра;
*[http://www.gdal.org/gdal_translate.html gdal_translate] - конвертация растров из формата в формат;
*[http://www.gdal.org/gdal2tiles.html gdal2tiles.py] - создание тайловой структуры, KML и простого просмотровщика;
*[http://www.gdal.org/gdaladdo.html gdaladdo] - добавление пирамидных слоёв (overview);
*[http://www.gdal.org/gdalbuildvrt.html gdalbuildvrt] - создание виртуального растра (VRT) из набора;
*[http://www.gdal.org/gdalcompare.html gdalcompare.py ] - сравнение двух изображений;
*[http://www.gdal.org/gdaldem.html gdaldem] - набор инструментов для анализа и визуализации ЦМР;
*[http://www.gdal.org/gdalinfo.html gdalinfo] - информация о растре;
*[http://www.gdal.org/gdallocationinfo.html gdallocationinfo] - запросы информации к растровыми файлам;
*[http://www.gdal.org/gdalmanage.html gdalmanage] - управление растровыми файлами (копирование, переименование, удаление и т.д.);
*[http://www.gdal.org/gdalmove.html gdalmove.py] - трансформирование системы координат растра без ресэмплирования;
*[http://www.gdal.org/gdalsrsinfo.html gdalsrsinfo] - показывается информацию по системе координат в разных форматах (WKT, PROJ.4 и др.);
*[http://www.gdal.org/gdaltindex.html gdaltindex] - построить индекс фрагментов (тайлов) MapServer;
*[http://www.gdal.org/gdaltransform.html gdaltransform]- трансформация координат;
*[http://www.gdal.org/gdalwarp.html gdalwarp] - трансформация изображения в новую систему координат;
*[http://www.gdal.org/nearblack.html nearblack] - конвертация черных/белых границ в нужное значение;
*[http://www.gdal.org/pct2rgb.html pct2rgb.py] - конвертация 8-битных изображений с палитрой в 24-битные RGB изображений;
*[http://www.gdal.org/rgb2pct.html rgb2pct.py] - конвертация 24-битных RGB изображений в 8-битные с палитрой;

Поддерживаемые форматы и используемые ключи можно узнать просто набрав в командной строке имя одной из утилит.

<pre>gdalinfo</pre>

В результате будет получена справка по использованию этой программы:

<pre>Usage: gdalinfo [--help-general] [-mm] [-stats] [-nogcp] [-nomd]
[-noct] [-checksum] [-mdd domain]* datasetname</pre>

Версию GDAL можно посмотреть командой:

<pre>gdalinfo --version</pre>

Список форматов поддерживаемых утилитами GDAL можно посмотреть следующим образом:

<pre>gdalinfo --formats</pre>

Cписок поддерживаемых форматов (список может отличаться как в большую, так и в меньшую сторону, поскольку зависит от того, были ли подключены/отключены соответствующие модули при компиляции программы):

*GRASS (ro): GRASS Database Rasters (5.7+)
*VRT (rw+): Virtual Raster
*GTiff (rw+): GeoTIFF
*HFA (rw+): Erdas Imagine Images (.img)
*AIG (ro): Arc/Info Binary Grid
*AAIGrid (rw): Arc/Info ASCII Grid
*JPEG (rw): JPEG JFIF
*MEM (rw+): In Memory Raster
*GIF (rw): Graphics Interchange Format (.gif)
*BMP (rw+): MS Windows Device Independent Bitmap
*DIMAP (ro): SPOT DIMAP
*PCIDSK (rw+): PCIDSK Database File
*SRTMHGT (rw): SRTMHGT File Format
*GMT (rw): GMT NetCDF Grid Format
*HDF4 (ro): Hierarchical Data Format Release 4
*HDF4Image (rw+): HDF4 Dataset
*ENVI (rw+): ENVI .hdr Labelled
*EHdr (rw+): ESRI .hdr Labelled

==Конвертация==

Извлечь три канала с номерами 1, 2, 3 в новый файл из исходного с перекомбинацией, в котором каналов может быть больше.

<pre>gdal_translate -b 3 -b 2 -b 1 output.tif input.tif</pre>

В результате в текущем каталоге появится результат 3-х канальный файл output.tif. Или в цикле для например 46-канального (win):

<pre>for /L %i in (1,1,46) DO gdal_translate -b %i input.tif output_%i.tif</pre>

Конвертация с обрезкой по заданным координатам:

<pre>gdal_translate -of GTiff -projwin 75.081940 57.250275 89.869980 49.083084 input.tif output.tiff</pre>

Конвертация с компрессией и созданием [http://gis-lab.info/qa/tfw.html world-файла]:

<pre>gdal_translate -co "COMPRESS=LZW" -co "worldfile=yes" input.tif output.tiff</pre>

Конвертация 16 битного одноканального растра в 8 битный:

<pre>gdal_translate -scale -ot Byte input_16bit.tif output.tif</pre>

Пакетная конвертация всех JPG в TIF (Windows):
<pre>for %i in (*.jpg) do gdal_translate %i %~ni.tif </pre>

==Стэки, композиты, мозаики==

Создание композитного изображения из серии отдельных растров, каждый из которых в своем файле TIF. Разрешение выходного файла устанавливается по первому их растров. Таким образом, если первый канал 15 м, а остальные 30 м, то последние будут пересчитаны на 15 м. Чтобы указать, что каждый исходный растр должен попасть в отдельный слой (layer-stacking), а не мозаицирование, используется ключ -separate:

<pre>gdal_merge.py -o output.tif band1.tif band2.tif band3.tif band4.tif band5.tif -separate</pre>

Для мозаицирования (объединения растров располагающихся в пространстве рядом друг с другом), этот ключ нужно убрать. Например чтобы склеить соседние фрагменты (тайлы) рельефа в единое поле:

<pre>gdal_merge -o altay.tif srtm_53_02.tif srtm_53_03.tif srtm_54_02.tif srtm_54_03.tif</pre>

А так можно указать разрешение выходного растра и то, что использовать определенное значение - не нужно (чтобы поля не перекрывали значащие части):

<pre>gdal_merge -n 0 -ps 0.00416 0.00416 -o output.tif in1.tif in2.tif in3.tif in4.tif</pre>

==Работа с NODATA==
Конвертация с заменой одного значения на другое (обычно используется для NODATA):

<pre>gdalwarp -srcnodata -999 -dstnodata 0 input.tif output.tif</pre>

Если в исходном растре nodata записано как nan <code>NoData Value=nan</code>, то конвертировать лучше так:

<pre>gdalwarp -dstnodata nan input.tif output.tif</pre>

Как снять флаг NODATA? Допустим есть растр в котором часть пикселей имеют значение 0 и на них установлен флаг NODATA и нужно его убрать. Делается это так:

<pre>gdal_translate -a_nodata none input_with_NoData.tif output_without_NoData.tif</pre>

==Работа с системами координат==
===Переназначение системы координат===
Если вам нужно просто перепрописать систему координат, без пересчета самого растра:

<pre>gdal_edit -a_srs "EPSG:4326" input.tif</pre>

Если код неизвестен, можно указать описание системы координат в формате Proj.4:

<pre>gdal_edit -a_srs "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs" input.tif</pre>

Узнать строку описания системы координат можно следующим образом:

<pre>gdalsrsinfo template.tif</pre>

===Конвертация с перепроецированием===

gdalwarp позволяет не только конвертировать данные из одного формата в другой, но и одновременно произвести перепроецирование данных из одной системы координат в другую. Для этого используются параметры:

*-a_srs используется для указания системы координат (СК) для данных
*-s_srs используется для перезаписи информации о системе координат
*-t_srs перепроецирования данных в требуемую систему координат

В самом простейшем случае это делается следующим образом (ключ -tr указывает разрешение целевого растра):

<pre>gdalwarp.exe -tr 0.0083 0.0083 -t_srs "EPSG:4326" in.tif out.tif</pre>

Если EPSG-кода у СК конечного растра нет, то указать целевую СК можно так:

<pre>gdalwarp.exe -t_srs "+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=45 +x_0=8500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs" in.tif out.tif</pre>

==Обрезка==
Обрезка по векторному контуру c уменьшением охвата растра (реальная обрезка, а не просто заполнение ненужных областей значениями NODATA):
<pre>gdalwarp -cutline aoi.shp -crop_to_cutline input.tif output.tif</pre>

==Работа с рельефом==
Теневая отмывка рельефа:

<pre>gdaldem hillshade altay.tif altayhill.tif -z 5 -s 111120</pre>
Ключ <code>-s 111120</code> используется для пересчета для растров сделанных в EPSG:4326 в метровые СК. Если исходник уже находится в проекции, то он не нужен.

Цветовая отмывка рельефа:

<pre>gdaldem color-relief altay.tif ramp.txt altay-color.tif</pre>

Пример файла ramp.txt:
<pre>5000 255 255 255
1000 168 112 0
650 198 165 48
400 229 218 97
200 218 229 97
0 112 168 0</pre>

==Работа с цветовой таблицей==
Если на входе есть растр с цветами в шкале серого (grayscale), а на выходе нужно получить изображение с индексированными цветами (палитрой),о можно воспользоваться подходом через VRT.

=== С использованием VRT ===

1. Создаем файл VRT:
<pre>gdal_translate -of VRT input.tif input.vrt</pre>

2. Открываем файл VRT в любом текстовом редакторе, находим:

<pre><ColorInterp>Gray</ColorInterp></pre>

и меняем на:

<pre><ColorInterp>Palette</ColorInterp>
<ColorTable>
<Entry c1="0" c2="0" c3="128"/>
<Entry c1="0" c2="128" c3="0"/>
<Entry c1="0" c2="255" c3="0"/>
<Entry c1="153" c2="204" c3="0"/>
<Entry c1="131" c2="66" c3="37"/>
<Entry c1="51" c2="153" c3="102"/>
</ColorTable>
</pre>

Каждая Entry - это цвет, кодированный в RGB и соответствующий по порядку значениям от 0 до 255. Если цветов меньше 255, то остальные будут добавлены автоматически с RGB = 0,0,0 (черный).

После того как таблица отредактирована, нужно сохранить ее обратно в растр:
<pre>gdal_translate input.vrt result.tif</pre>

==Построение изолиний==
Утилита gdal_contour используется для получения изолиний - линий равных значений по растровым данным.
Полученные линии пересекают все пиксели с одинаковым значением, очерчивая при этом некоторую область.
Чаще всего применяется для построения горизонталей рельефа из ЦМР.

Построение контуров с интервалом в 5 единиц (Единица указывается в единицах измерения исходного растра):

<pre>gdal_contour -i 5 mydem.tif contour.shp</pre>

Построение контуров из первого канала растра, с интервалом в 100 единиц начиная с 1200 и записью значения в поле elev:

<pre>gdal_contour -b 1 -a elev -i 100 -off 1200 mydem.tif contour.shp</pre>

Построение только контуров с фиксированными значениями 1000, 1100 и 1120 и выводом результата в таблицу PostGIS

<pre>gdal_contour -a elev -f PostgreSQL -fl 1000 1100 1120 -nln cont mydem.tif "PG:host=localhost user=iampg password=iampgpass dbname=iamgis"</pre>

==netCDF==

Конвертирование в GeoTIFF:

<pre>gdal_translate -of GTiff -b 1 NETCDF:precip.mon.mean.nc:precip b1.tif</pre>

Конвертирование в GeoTIFF с обрезкой по исходным координатам и созданием TFW (world-)файла:

<pre>gdal_translate -of GTiff -srcwin 0 0 72 72 -co TFW=YES -b 1 NETCDF:precip.mon.mean.nc:precip b1.tif</pre>

Разбиение поднаборов данных ("SUBDATASET") на отдельные файлы netCDF (на выходе — файлы типа "example1", "example2" и т.д.):
<pre>gdal_translate -sds example.nc example</pre>

==HDF4==

В HDF4 распространяется множество данных дистанционного зондирования, например MODIS и ASTER.

Импорт данных ASTER L1A:

<pre>gdalwarp -overwrite -of GTiff HDF4_EOS:EOS_SWATH:"110601_081441.hdf":VNIR_Band1:ImageData b1.tif </pre>

Импорт данных MODIS Level 3,4:

<pre>gdal_translate HDF4_EOS:EOS_GRID:"MCD12Q1.A2001001.h00v08.051.2014287161513.hdf":MOD12Q1:Land_Cover_Type_1 output.tif</pre>

Если одна из частей названия SDS (массива данных внутри HDF) имеет пробелы, ее нужно взять в кавычки. Использовать -geoloc для перепроецирования не нужно.

==Расчеты==
Для разнообразных пересчетов значений пикселей можно использовать gdal_calc. Например, такая команда сбросит в NODATA все пиксели чьё значение больше 16000:

<pre>gdal_calc.bat -A input.tif --outfile=output.tif --calc="A*(A<16000)" --NoDataValue=0</pre>

Читается такое выражение следующим образом: для каждого пикселя растра input.tif, если его значение меньше 16000 - оставить его таким же, иначе - сбросить в 0.

В выражениях могут использоваться логические "и" и "или", т.е. например в примере ниже: если значение пикселя больше или равно 249, но меньше 255 - сделать его единицей, остальные значения (включая 255) установить в 0.

<pre>gdal_calc.bat -A input.tif --outfile=output.tif --calc="1*(logical_and(A>=249,A<255)) " --NoDataValue=0'</pre>

Следует обратить внимание, что часть ELSE выражения по умолчанию всегда будет сбрасывать остальные значения в 0. При этом, часть --NoDataValue=0 только говорит к какому значению нужно "приклеить" ярлык NODATA. Таким образом, если в выражении выше сказать, например, --NoDataValue=255, то, несмотря на то, что NODATA установится на 255, значения не удовлетворяющие условию в calc вовсе не станут равны 255, а останутся равны 0. Что бы установить остальные значения в число отличное от 0, нужно сделать это явным образом, например:

<pre>gdal_calc --overwrite -A input.tif --outfile=input.tif --calc="1*(A==0) + 255*(A>0)" --NoDataValue=255</pre>

Присвоить пикселям со значением 0 значение 1, а остальным присвоить 255, установив на это значение тег NODATA.

==Создание растров уменьшенного разрешения (т.н. quicklook, preview)==
Создать для данных ДЗЗ высокого разрешения так называемый quicklook, т.е. привязанный растр для быстрого предварительного просмотра, можно с помощью gdal_translate, ему можно просто указать конечный размер процентах:

<pre>gdal_translate -of "JPEG" -outsize 20% 20% ALOS_example.tif ALOS_example_preview.jpg -co "WORLDFILE=YES"</pre>

Или с помощью gdalwarp, ему можно передать нужный размер пикселя в единицах системы координат, так же можно указать метод, которым будет происходить объединение значений (например усреднение):

<pre>gdalwarp -tr 8000 8000 -r average input.tif output.tif</pre>

или, если на выходе нужен растр определенного размера:

<pre>gdalwarp -ts 8000 8000 -r average input.tif output.tif</pre>

==Ссылки по теме==
*[http://www.gdal.org/index.html Библиотеки GDAL/OGR]
*[http://gis-lab.info/qa/ogr2ogr-examples.html Примеры использования ogr2ogr]

Обсуждение участника:Jackpit8888

2023-07-15T07:41:04Z

Максим Дубинин: Полностью удалено содержимое страницы

2018-05-04T10:23:43Z

Максим Дубинин:

==Загружаем исходные данные==
Источник данных: http://stat.gibdd.ru

Находим код региона.

Указываем диапазон дат.

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

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

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

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

==Находим пересечения==
Источник данных: http://data.nextgis.com
Совсем свежие данные стоят 200-300р, старые можно скачать бесплатно. Так или иначе получаем свой регион. Нам нужен только один слой: highway-lines.shp (дорожная сеть).

Находим все пересечения и считаем в каждом количество примыканий, определяем таким образом Х-образные пересечения (у них это значение >= 4). Но можете оставить и другие, если нужно.

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

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

Проверяем, что биндинги работают. Выполняем эту команду - она не должна выводить ничего.
<syntaxhighlight lang="python">
python -c "import gdal"
</syntaxhighlight>

Клонируем скрипт

<syntaxhighlight lang="bash">
git clone https://github.com/nextgis/data_processing_scripts.git
cd data_processing_scripts/nextgis_extracts_detect_crossings
</syntaxhighlight>

При запуске скрипта нужно только указать путь к highway-line.shp, который распаковали на предыдущем шаге и задать фильтр по типам дорог - это то же самое, если бы вы отфильтровали highway-lines.shp в ogr2ogr. Этот фильтр вы задаёте в зависимости от конкретно вашего алгоритма анализа, в зависимости от его настройки можно получить аварии только на перекрёстках крупных улиц, или только на перекрёстках мелких улиц.
<syntaxhighlight lang="bash">
python detect_crossings.py -s highway-line.shp -d crossings.shp -f "HIGHWAY IN ('primary')"
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>
Фильтр нужен обязательно, потому что без него у вас получатся перекрёстки тротуаров с тропинками в лесу

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

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

==Создаём перекрестки==
Перекресток - зона определенного радиуса вокруг пересечения (буферная зона).

Перекрестки получаются путем построения буферных зон вокруг пересечений. Соединенные при этом зоны - объединяются в единый объект.

На входе - точки (crossings), точки событий (dtp2014).

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

Расчёт делается в PostGIS. Хорошая инструкция по установке есть на https://www.digitalocean.com/community/tutorials/how-to-install-and-configure-postgis-on-ubuntu-14-04

Загружаем слои из файлов в PostGIS. Подставьте в эти команды логин и пароль к БД.
<syntaxhighlight lang="bash">
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" ../moscow_crossings.gpkg

ogr2ogr -progress -overwrite -f GPKG dtp2014.gpkg --config SHAPE_ENCODING cp1251 dtp_post_14_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2014.gpkg -skipfailures -nln dtp2014

ogr2ogr -progress -overwrite -f GPKG dtp2015.gpkg --config SHAPE_ENCODING cp1251 dtp_post_15_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2015.gpkg -skipfailures -nln dtp2015

ogr2ogr -progress -overwrite -f GPKG dtp2016.gpkg --config SHAPE_ENCODING cp1251 dtp_post_16_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2016.gpkg -skipfailures -nln dtp2016
</syntaxhighlight>

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

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

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

<syntaxhighlight lang="sql">
-- генерируем три слоя буферов с разным расстоянием
-- Укажите тут нужную СК, если это не Москва.
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;
CREATE TEMPORARY TABLE buffers50 ON COMMIT DROP AS
SELECT st_buffer(ST_Transform(wkb_geometry,32637),50) AS wkb_geometry,50 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;
CREATE TEMPORARY TABLE buffers20 ON COMMIT DROP AS
SELECT st_buffer(ST_Transform(wkb_geometry,32637),20) AS wkb_geometry,20 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;

DROP TABLE IF EXISTS crossingsarea_disolved;
CREATE TABLE crossingsarea_disolved (
wkb_geometry geometry,
buffer integer);

--для 100 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers100),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 100 as buffer FROM crossingsarea_disolved_onebuffer;

--для 50 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers50),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 50 as buffer FROM crossingsarea_disolved_onebuffer;

--для 20 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers20),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 20 as buffer FROM crossingsarea_disolved_onebuffer;

--добавление столбцов и рассчёт числовых значений
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_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2014 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2015 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2014 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2015 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cp_sum integer;
COMMENT ON COLUMN crossingsarea_disolved.cp_sum IS 'КОЛИЧЕСТВО АВАРИЙ С ТРУПАМИ';
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_sum integer;
COMMENT ON COLUMN crossingsarea_disolved.cp_sum IS 'КОЛИЧЕСТВО ТРУПОВ';

UPDATE crossingsarea_disolved SET wkb_geometry=ST_SetSRID(wkb_geometry,32637);
UPDATE crossingsarea_disolved SET count_2014=0;
UPDATE crossingsarea_disolved SET count_2015=0;
UPDATE crossingsarea_disolved SET count_2016=0;
UPDATE crossingsarea_disolved SET cp_2014=0;
UPDATE crossingsarea_disolved SET cp_2015=0;
UPDATE crossingsarea_disolved SET cp_2016=0;
UPDATE crossingsarea_disolved SET cps_2014=0;
UPDATE crossingsarea_disolved SET cps_2015=0;
UPDATE crossingsarea_disolved SET cps_2016=0;

CREATE TEMPORARY TABLE medium2014 ON COMMIT DROP AS
SELECT
dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2014 , crossingsarea_disolved
WHERE
ST_Within(dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2014_cnt ON COMMIT DROP AS
SELECT medium2014.crossing_geometry, COUNT(*) AS cnt FROM medium2014 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2014 = medium2014_cnt.cnt
FROM medium2014_cnt
WHERE medium2014_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2014_cnt;

CREATE TEMPORARY TABLE medium2014_cnt ON COMMIT DROP AS
SELECT medium2014.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2014 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2014 = medium2014_cnt.cnt, cps_2014 = medium2014_cnt.corpses
FROM medium2014_cnt
WHERE medium2014_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2014_cnt;

CREATE TEMPORARY TABLE medium2015 ON COMMIT DROP AS
SELECT
dtp2015.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2015 , crossingsarea_disolved
WHERE
ST_Within(dtp2015.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2015_cnt ON COMMIT DROP AS
SELECT medium2015.crossing_geometry, COUNT(*) AS cnt FROM medium2015 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2015 = medium2015_cnt.cnt
FROM medium2015_cnt
WHERE medium2015_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2015_cnt;

CREATE TEMPORARY TABLE medium2015_cnt ON COMMIT DROP AS
SELECT medium2015.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2015 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2015 = medium2015_cnt.cnt, cps_2015 = medium2015_cnt.corpses
FROM medium2015_cnt
WHERE medium2015_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2015_cnt;

CREATE TEMPORARY TABLE medium2016 ON COMMIT DROP AS
SELECT
dtp2016.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2016 , crossingsarea_disolved
WHERE
ST_Within(dtp2016.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2016_cnt ON COMMIT DROP AS
SELECT medium2016.crossing_geometry, COUNT(*) AS cnt FROM medium2016 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2016 = medium2016_cnt.cnt
FROM medium2016_cnt
WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2016_cnt;

CREATE TEMPORARY TABLE medium2016_cnt ON COMMIT DROP AS
SELECT medium2016.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2016 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2016 = medium2016_cnt.cnt, cps_2016 = medium2016_cnt.corpses
FROM medium2016_cnt
WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2016_cnt;

UPDATE crossingsarea_disolved SET cp_sum = cp_2014+cp_2015+cp_2016;
UPDATE crossingsarea_disolved SET cps_sum = cps_2014+cps_2015+cps_2016;
</syntaxhighlight>

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

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

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

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

Создаём свою карту ДТП

2018-05-04T10:21:20Z

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

Создаём свою карту ДТП

2018-05-04T10:21:10Z

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

==Загружаем исходные данные==
Источник данных: http://stat.gibdd.ru

Находим код региона.

Указываем диапазон дат.

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

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

==Преобразовываем в табличный вид==
С XML работать не удобно, поэтому используем следующий скрипт для преобразования XML в CSV.

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

==Находим пересечения==
Источник данных: http://data.nextgis.com
Совсем свежие данные стоят 200-300р, старые можно скачать бесплатно. Так или иначе получаем свой регион. Нам нужен только один слой: highway-lines.shp (дорожная сеть).

Находим все пересечения и считаем в каждом количество примыканий, определяем таким образом Х-образные пересечения (у них это значение >= 4). Но можете оставить и другие, если нужно.

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

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

Проверяем, что биндинги работают. Выполняем эту команду - она не должна выводить ничего.
<syntaxhighlight lang="python">
python -c "import gdal"
</syntaxhighlight>

Клонируем скрипт

<syntaxhighlight lang="bash">
git clone https://github.com/nextgis/data_processing_scripts.git
cd data_processing_scripts/nextgis_extracts_detect_crossings
</syntaxhighlight>

При запуске скрипта нужно только указать путь к highway-line.shp, который распаковали на предыдущем шаге и задать фильтр по типам дорог - это то же самое, если бы вы отфильтровали highway-lines.shp в ogr2ogr. Этот фильтр вы задаёте в зависимости от конкретно вашего алгоритма анализа, в зависимости от его настройки можно получить аварии только на перекрёстках крупных улиц, или только на перекрёстках мелких улиц.
<syntaxhighlight lang="bash">
python detect_crossings.py -s highway-line.shp -d crossings.shp -f "HIGHWAY IN ('primary')"
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>
Фильтр нужен обязательно, потому что без него у вас получатся перекрёстки тротуаров с тропинками в лесу

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

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

==Создаём перекрестки==
Перекресток - зона определенного радиуса вокруг пересечения (буферная зона).

Перекрестки получаются путем построения буферных зон вокруг пересечений. Соединенные при этом зоны - объединяются в единый объект.

На входе - точки (crossings), точки событий (dtp2014).

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

Расчёт делается в PostGIS. Хорошая инструкция по установке есть на https://www.digitalocean.com/community/tutorials/how-to-install-and-configure-postgis-on-ubuntu-14-04

Загружаем слои из файлов в PostGIS. Подставьте в эти команды логин и пароль к БД.
<syntaxhighlight lang="bash">
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" ../moscow_crossings.gpkg

ogr2ogr -progress -overwrite -f GPKG dtp2014.gpkg --config SHAPE_ENCODING cp1251 dtp_post_14_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2014.gpkg -skipfailures -nln dtp2014

ogr2ogr -progress -overwrite -f GPKG dtp2015.gpkg --config SHAPE_ENCODING cp1251 dtp_post_15_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2015.gpkg -skipfailures -nln dtp2015

ogr2ogr -progress -overwrite -f GPKG dtp2016.gpkg --config SHAPE_ENCODING cp1251 dtp_post_16_all_Mow.shp -skipfailures
ogr2ogr -progress -overwrite PG:"host=localhost dbname=gis user= password=" dtp2016.gpkg -skipfailures -nln dtp2016
</syntaxhighlight>

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

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

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

<syntaxhighlight lang="sql">
-- генерируем три слоя буферов с разным расстоянием
-- Укажите тут нужную СК, если это не Москва.
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;
CREATE TEMPORARY TABLE buffers50 ON COMMIT DROP AS
SELECT st_buffer(ST_Transform(wkb_geometry,32637),50) AS wkb_geometry,50 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;
CREATE TEMPORARY TABLE buffers20 ON COMMIT DROP AS
SELECT st_buffer(ST_Transform(wkb_geometry,32637),20) AS wkb_geometry,20 AS buffer_size FROM moscow_crossings WHERE ways_cnt = 4;

DROP TABLE IF EXISTS crossingsarea_disolved;
CREATE TABLE crossingsarea_disolved (
wkb_geometry geometry,
buffer integer);

--для 100 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers100),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 100 as buffer FROM crossingsarea_disolved_onebuffer;

--для 50 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers50),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 50 as buffer FROM crossingsarea_disolved_onebuffer;

--для 20 метров
DROP TABLE IF EXISTS crossingsarea_disolved_onebuffer ;
CREATE
TEMPORARY
TABLE crossingsarea_disolved_onebuffer
ON COMMIT DROP
AS
WITH
clusters(wkb_geometry) AS
(SELECT ST_CollectionExtract(unnest(ST_ClusterIntersecting(wkb_geometry)), 3)
FROM buffers20),
multis(id, wkb_geometry) AS
(SELECT row_number() over() as id, wkb_geometry FROM clusters)
SELECT ST_UNION(wkb_geometry) AS wkb_geometry FROM
(SELECT id, (ST_DUMP(wkb_geometry)).geom AS wkb_geometry FROM multis) d GROUP BY id;
INSERT INTO crossingsarea_disolved SELECT wkb_geometry, 20 as buffer FROM crossingsarea_disolved_onebuffer;

--добавление столбцов и рассчёт числовых значений
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_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2014 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2015 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cp_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2014 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2015 integer;
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_2016 integer;

ALTER TABLE crossingsarea_disolved ADD COLUMN cp_sum integer;
COMMENT ON COLUMN crossingsarea_disolved.cp_sum IS 'КОЛИЧЕСТВО АВАРИЙ С ТРУПАМИ';
ALTER TABLE crossingsarea_disolved ADD COLUMN cps_sum integer;
COMMENT ON COLUMN crossingsarea_disolved.cp_sum IS 'КОЛИЧЕСТВО ТРУПОВ';

UPDATE crossingsarea_disolved SET wkb_geometry=ST_SetSRID(wkb_geometry,32637);
UPDATE crossingsarea_disolved SET count_2014=0;
UPDATE crossingsarea_disolved SET count_2015=0;
UPDATE crossingsarea_disolved SET count_2016=0;
UPDATE crossingsarea_disolved SET cp_2014=0;
UPDATE crossingsarea_disolved SET cp_2015=0;
UPDATE crossingsarea_disolved SET cp_2016=0;
UPDATE crossingsarea_disolved SET cps_2014=0;
UPDATE crossingsarea_disolved SET cps_2015=0;
UPDATE crossingsarea_disolved SET cps_2016=0;

CREATE TEMPORARY TABLE medium2014 ON COMMIT DROP AS
SELECT
dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2014 , crossingsarea_disolved
WHERE
ST_Within(dtp2014.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2014_cnt ON COMMIT DROP AS
SELECT medium2014.crossing_geometry, COUNT(*) AS cnt FROM medium2014 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2014 = medium2014_cnt.cnt
FROM medium2014_cnt
WHERE medium2014_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2014_cnt;

CREATE TEMPORARY TABLE medium2014_cnt ON COMMIT DROP AS
SELECT medium2014.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2014 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2014 = medium2014_cnt.cnt, cps_2014 = medium2014_cnt.corpses
FROM medium2014_cnt
WHERE medium2014_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2014_cnt;

CREATE TEMPORARY TABLE medium2015 ON COMMIT DROP AS
SELECT
dtp2015.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2015 , crossingsarea_disolved
WHERE
ST_Within(dtp2015.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2015_cnt ON COMMIT DROP AS
SELECT medium2015.crossing_geometry, COUNT(*) AS cnt FROM medium2015 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2015 = medium2015_cnt.cnt
FROM medium2015_cnt
WHERE medium2015_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2015_cnt;

CREATE TEMPORARY TABLE medium2015_cnt ON COMMIT DROP AS
SELECT medium2015.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2015 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2015 = medium2015_cnt.cnt, cps_2015 = medium2015_cnt.corpses
FROM medium2015_cnt
WHERE medium2015_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2015_cnt;

CREATE TEMPORARY TABLE medium2016 ON COMMIT DROP AS
SELECT
dtp2016.wkb_geometry, crossingsarea_disolved.wkb_geometry AS crossing_geometry,corpse
FROM
dtp2016 , crossingsarea_disolved
WHERE
ST_Within(dtp2016.wkb_geometry, crossingsarea_disolved.wkb_geometry);

CREATE TEMPORARY TABLE medium2016_cnt ON COMMIT DROP AS
SELECT medium2016.crossing_geometry, COUNT(*) AS cnt FROM medium2016 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET count_2016 = medium2016_cnt.cnt
FROM medium2016_cnt
WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2016_cnt;

CREATE TEMPORARY TABLE medium2016_cnt ON COMMIT DROP AS
SELECT medium2016.crossing_geometry, COUNT(*) AS cnt, SUM(corpse) AS corpses FROM medium2016 WHERE corpse > 0 GROUP BY crossing_geometry;

UPDATE crossingsarea_disolved SET cp_2016 = medium2016_cnt.cnt, cps_2016 = medium2016_cnt.corpses
FROM medium2016_cnt
WHERE medium2016_cnt.crossing_geometry = crossingsarea_disolved.wkb_geometry;
DROP TABLE medium2016_cnt;

UPDATE crossingsarea_disolved SET cp_sum = cp_2014+cp_2015+cp_2016;
UPDATE crossingsarea_disolved SET cps_sum = cps_2014+cps_2015+cps_2016;
</syntaxhighlight>

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

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

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

Donate

2018-05-03T09:37:33Z

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

{{Организационное|donate}}
{{Аннотация|Пожертвования на работу GIS-Lab}}

Наш ресурс существует уже 14 лет, но мы, за исключением 2 случаев специальных сборов на апгрейд сервера, никогда не собирали деньги с наших пользователей, обходились своими ресурсами. Тем временем на нашем ресурсе накопились задачи, требующие денег. Пожертвования являются простым способом поддержать GIS-Lab. Если финансы этого не позволяют, но есть желание поддержать еще как-то, то есть и [http://gis-lab.info/help.html другие способы].

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

На что идут собранные деньги:

* Системное администрирование (периодически)
* Исправление инструментов сайта (периодически)
* Оплата услуг дизайнеров, разработчиков и т.д. (периодически)
* Организация мероприятий (периодически)

Вы можете воспользоваться этой формой, чтобы сделать своё пожертвование, сумму можно изменить на любую:

{{#widget:Iframe
|url=http://gis-lab.info/qa/donate-form.html
|width=740
|border=0
}}

Внимание: Если вы не видите форму, то ее скорее всего скрыло одно из расширений браузера для блокировки баннеров и т.п. вещей. Если вы хотите сделать пожертвование, временно отключите эти расширения.

==Регулярные пожертвования==
Самым эффективным способом поддержки являются регулярные пожертвования.
===Инструкция для Сбербанк-онлайн===
#Зайти в приложение.
#Выбрать карту для оплаты
#Платежи и переводы
#В поиске вбить Яндекс , В списке появится одна строка "Яндекс.Деньги прочие услуги", выбираем его.
#Выбрать "по номеру кошелька"
#Вбиваем номер кошелька 410011586105267 в поле Номер счета
#Сумма - введите сумму
#Продолжить, подтвердить.

==Отчет о расходах==
*03.05.2018: 3000 продление на год gisconf.ru, gis-lab.ru, open-gis.ru, остаток: 20906
*19.04.2018: 825 CleanTalk на год, остаток: 22632
*13.08.2017: 20000 перевод IgorZ за успешное администрирование в течение многих лет, остаток: 9692
*13.08.2017: 300 СМС информирование для Яндекс.Кошелька до 16.11.2017, остаток: 25656
*20.06.2017: 1520 продление домена gis-lab.ru и gisconf.ru (nic.ru)
*27.04.2017: 4500 хостинг (hc.ru)
*25.03.2017: 4500 хостинг (hc.ru)
*25.01.2017: 4500 хостинг (hc.ru)
*23.11.2016: 4000 хостинг (hc.ru)
*29.10.2016: 4500 хостинг (hc.ru)
*23.09.2016: 4500 хостинг (hc.ru)
*28.08.2016: 5000 хостинг (hc.ru)
*16.08.2016: 240 смс-информирование
*25.07.2016: 4000 хостинг (hc.ru)
*24.05.2016: 5000 хостинг (hc.ru)
*23.03.2016: 5000 хостинг (hc.ru)
*25.02.2016: 4000 хостинг (hc.ru)
*04.02.2016: 2000 год хостинга 1 Тб на яндекс.диске для папки aster
*04.02.2016: 2500 оплата апгрейда форума, первый этап
*24.01.2016: 5000 хостинг (hc.ru)
*22.12.2015: 4000 хостинг (hc.ru)
*26.11.2015: 5000 хостинг (hc.ru)
*27.10.2015: 5000 хостинг (hc.ru)
*24.09.2015: 5000 хостинг (hc.ru)
*29.08.2015: 2000 починка блога (сломалась поддержка кириллицы при обновлении)
*27.08.2015: 5000 хостинг (hc.ru)
*23.08.2015: 3000 дизайн баннеров
*25.07.2015: 4000 хостинг (hc.ru)
*20.07.2015: 240 смс-информирование

[[Категория:Служебные]]

Zonal statistics with ZonalStats for QGIS

2018-04-24T11:31:44Z

Максим Дубинин:

{{Статья|Опубликована|zonalstats-qgis-eng}}
{{Аннотация|Description of the new tool for zonal statistics calculation.}}

Thematic classification usually followed by an analysis of the results. One of the types of such analysis is zonal statistics calculation.

QGIS laready contains Zonal statistic plugin, that calculates several values (sum, mean value, total count) for pixels by polygonal vector layer. But this plugin don't allow group features by some field, that is necessary when one object is represented by several features. Calculated statistics added to shapefile attributes and this complicates it's usage for further analysis. And also for singleband binary rasters there is no way to get value of used by pixels area.

For this task we developed ZonalStats plugin for opensource GIS QGIS, that allows to calculate zonal statistics and show short report.

This development is part of the project to facilitate monitoring of FSC (Forest Stewardship Council) certified forestry enterprises in Russia.

{|
|[[Файл:Nextgis.png|link=http://nextgis.ru]]
|Open source geospatial solutions
|}

=== Installation notes ===

Plugin is available via [http://plugins.qgis.org/plugins/plugins.xml official repository].

ZonalStats is designed for QGIS 1.7.2 or higher and have no additional dependencies.

You can obtain sources from our Git:

<pre>git clone https://github.com/nextgis/zonalstats</pre>

=== How does it work ===

After plugin installation and startup by clicking

[[Файл:Zonalstats-icon.png|center]]

main window will appear

[[Файл:zonalstats-qgis-02-eng.png|center]]

At the "Settings" tab you can specify all necessary parameters.

First, select raster layer for analysis from "Raster layer" combobox. Usually this is singleband raster created as result of the thematic classification.

'''IMPORTANT'''! Plugin mainly developed to analyze classification results from [http://gis-lab.info/qa/dtclassifier.html DTclassifier] plugin (0 — feature absence; 1 — feature presence). To use ZonalStats with rasters that have other "encoding" you need to convert them.

Combobox "Polygon layer (zones)" used to select input polygonal vector layer with zones boundaries. Field that will be used for zones labeling selected in "Get zones names from" combobox. If checkbox "Group zones by name" is checked, than zones with identical value of the selected field will be counted as single zone. By default this checkbox is unchecked and each feature treated as single zone.

Field "Save report as" used to select filename for analysis report in HTML format.Additionaly plugin creates CSV file with analysis results that later can be easy loaded in statistical software for further analysis. By default this data file created near the report file, filename is same as for report but with suffix "_data". In you want to save data into another directory tick checkbox "Save data as" and select desired destination.

To start analysis press "OK" button.

ZonalStats counts number of pixels, that corresponds to some class and calculates area used by this class.

You can view results directly in plugin in the "Report" tab.

[[Файл:zonalstats-qgis-03-eng.png|center]]

Report generated by plugin contains three columns:

* firts column contains zones names (from shape field selected by user)
* "Object count" — number of features in zone. Makes sense, only when groupping is enabled. This column displayed number of features that treated as single zone (number of features with same name)
* "Area" — area used by class. Measurement units depends on layer units (meters, kilometers, feets...)

=== Contacts ===

If you want report a bug or want make suggestion — use [https://github.com/nextgis/zonalstats/issues bugtracker]. Also if you have a question about the script or plugin usage, you can [http://www.nextgis.ru/contact/ contact us] directly.

{|
|[[Файл:Nextgis.png|link=http://nextgis.ru]]
|Open source geospatial solutions
|}

=== Links ===

* [http://gis-lab.info/qa/dtclassifier-eng.html Raster classification with DTclassifier for QGIS]
* [http://gis-lab.info/qa/rastercalc-eng.html Raster algebra in QGIS with RasterCalc]
* [http://gis-lab.info/qa/qgis-install-plugin.html Установка модулей расширения в QGIS]
* [http://gis-lab.info/docs/qgis/ Документация по QGIS]

Данные OpenStreetMap в формате shape-файлов

2018-04-23T12:43:03Z

Максим Дубинин: /* Преимущества */

{{Статья|Опубликована|osmshp}}
{{Аннотация|Здесь можно скачать данные OpenStreetMap на любой регион в разных форматах.}}

== Общая информация ==

В рамках данного проекта, из данных OpenStreetMap и других открытых источников создаются обновляемые наборы слоев по любой точке мира, включая страны бывшего СССР и все регионы РФ. Данные наборы слоев доступны в форматах ESRI Shapefile, Mapinfo TAB, GeoJSON, что позволяет использовать их практически в любой ГИС.
{{NextGIS}}

Если вам нужны данные в формате XML или PBF, то их можно найти на [http://gis-lab.info/projects/osm_dump/ другой] странице.

После 12 сентября 2012 года, данные публикуются под лицензией [http://www.opendatacommons.org/licenses/odbl ODBL]. Подробнее о ваших правах и обязанностях можно почитать в [http://wiki.openstreetmap.org/wiki/Legal_FAQ Legal FAQ].

Данные для загрузки: http://data.nextgis.com/osmshp/ 
Cписок слоев: http://data.nextgis.com/osmshp/layer/

Система координат данных: WGS84 
Проекция: широта/долгота

== Преимущества ==
Эти данные обладают рядом особенностей и преимуществ относительно исходных данных в формате OSM XML:

* Данные представлены в распространенных ГИС форматах
* Добавлены готовые проекты для ГИС, благодаря которым данные можно немедленно начинать использовать в ГИС
* Данные разбиты по слоям, у каждого слоя своя четкая структура атрибутов ([https://docs.google.com/spreadsheets/d/1F83dtRH8c7O83E55ox3Kfh8Ibbh2TFL70nF5Iw_33d0/edit?usp=sharing посмотреть])
* Объекты слоёв административно-территориального деления и населенных пунктов имеют коды ОКТМО (источник - Росстат)
* Добавлены стили, условные обозначения

Дополнительные слои, отсутствующие в таком виде в исходных данных:
* Административно-территориальное деление в форме обрезанной по границе суши (исходный вариант так же сохранён)
* Береговые линии, границы суши и водной поверхности

== Условия распространения ==
Данные на июль 2017 г. можно скачать бесплатно (кроме РФ целиком). Для скачивания необходимо получить [http://my.nextgis.com NextGIS ID].

Можно так же заказать:
* разовую платную выгрузку, стоимость от 200р.
* регулярную выгрузку с периодичностью раз в день, каждую неделю или месяц (до полугода).
* выгрузку по произвольному контуру, стоимость от 300р.

Актуальность исходной базы данных - текущий день.

Оплата осуществляется банковской картой или Яндекс.Деньги. Если заказ сделан в выходной день, выгрузка производится в ближайший рабочий. Если заказ сделан в конце рабочего дня, результат могут прислать на следующий. Если вам не прислали результат в течение 2 рабочих дней, оплата будет возвращена.

Выгрузка данных выполняется сразу после заказа. Готовность данных - в среднем 10-15 минут.

== Замечания по форматам выгрузки ==

* ESRI Shape ([https://nextgis.com/data/examples/shapefile-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs) или проект для ArcGIS (файл *.mxd)
* GeoJSON ([https://nextgis.com/data/examples/geojson-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs)
* ESRI Geodatabase ([https://nextgis.com/data/examples/geodatabase-arcgis-example.7z пример]))
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для ArcGIS (файл *.mxd)
* Mapinfo TAB ([https://nextgis.com/data/examples/mapinfo-example.7z пример])
** Кодировка атрибутивной информации - Windows 1251 (CP1251)
** В комплект включены проект для Mapinfо (файл *.wor)
* CSV ([https://nextgis.com/data/examples/csv-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
* PDF ([https://nextgis.com/data/examples/pdf-example.pdf пример])
** Максимальный размер листа - 3 на 3 метра.
** Для большей детализации рекомендуется заказывать не целый регион, а конкретную область с помощью инструмента [http://data.nextgis.com/osmshp/custom_region указания области].

Возможна выгрузка и в другие форматы по запросу.

== Проекты для QGIS ==

В архивы с данными в формате ESRI Shape и GeoJSON включаются несколько вариантов проектов для QGIS.

{|class="wikitable" width="100%"
| width="0%" | [[Файл:OSM shp qgis mapnik.png|300px]]
| width="100%" style="vertical-align:top" |
==== qgis_mapnik ====

Этот вариант проектов был разработан [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8933 old_Bibigon] в рамках [http://gis-lab.info/qa/osm-qgis-projects.html отдельного проекта].

|-
| [[Файл:OSM shp qgis mapnik lc.png|300px]]
| width="100%" style="vertical-align:top" |

==== qgis_mapnik_lowcontrast ====

Вариант оформления предоставленный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8650 sanyasi] и дополненный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=11960 gisuser] на форуме. Аналогично предыдущему варианту оформления - тоже вариации на тему слоя Mapnik в OpenStreetMap, но с пониженной контрастностью, что может быть удобно в случае использования в качестве базовой карты.
|}

Если у вас есть предложения по корректировке этих вариантов, или вы хотите предложить свой вариант, то сделать это можно на форуме в теме-обсуждении данной статьи.

== Ссылки ==
* [http://nextgis.ru/contact Обратная связь]
* [http://gis-lab.info/projects/osm_dump/ Данные OpenStreetMap в форматах XML и PBF]

Donate

2018-04-19T12:00:15Z

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

{{Организационное|donate}}
{{Аннотация|Пожертвования на работу GIS-Lab}}

Наш ресурс существует уже 14 лет, но мы, за исключением 2 случаев специальных сборов на апгрейд сервера, никогда не собирали деньги с наших пользователей, обходились своими ресурсами. Тем временем на нашем ресурсе накопились задачи, требующие денег. Пожертвования являются простым способом поддержать GIS-Lab. Если финансы этого не позволяют, но есть желание поддержать еще как-то, то есть и [http://gis-lab.info/help.html другие способы].

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

На что идут собранные деньги:

* Системное администрирование (периодически)
* Исправление инструментов сайта (периодически)
* Оплата услуг дизайнеров, разработчиков и т.д. (периодически)
* Организация мероприятий (периодически)

Вы можете воспользоваться этой формой, чтобы сделать своё пожертвование, сумму можно изменить на любую:

{{#widget:Iframe
|url=http://gis-lab.info/qa/donate-form.html
|width=740
|border=0
}}

Внимание: Если вы не видите форму, то ее скорее всего скрыло одно из расширений браузера для блокировки баннеров и т.п. вещей. Если вы хотите сделать пожертвование, временно отключите эти расширения.

==Регулярные пожертвования==
Самым эффективным способом поддержки являются регулярные пожертвования.
===Инструкция для Сбербанк-онлайн===
#Зайти в приложение.
#Выбрать карту для оплаты
#Платежи и переводы
#В поиске вбить Яндекс , В списке появится одна строка "Яндекс.Деньги прочие услуги", выбираем его.
#Выбрать "по номеру кошелька"
#Вбиваем номер кошелька 410011586105267 в поле Номер счета
#Сумма - введите сумму
#Продолжить, подтвердить.

==Отчет о расходах==
*19.04.2018: 825 CleanTalk на год, остаток: 22632
*13.08.2017: 20000 перевод IgorZ за успешное администрирование в течение многих лет, остаток: 9692
*13.08.2017: 300 СМС информирование для Яндекс.Кошелька до 16.11.2017, остаток: 25656
*20.06.2017: 1520 продление домена gis-lab.ru и gisconf.ru (nic.ru)
*27.04.2017: 4500 хостинг (hc.ru)
*25.03.2017: 4500 хостинг (hc.ru)
*25.01.2017: 4500 хостинг (hc.ru)
*23.11.2016: 4000 хостинг (hc.ru)
*29.10.2016: 4500 хостинг (hc.ru)
*23.09.2016: 4500 хостинг (hc.ru)
*28.08.2016: 5000 хостинг (hc.ru)
*16.08.2016: 240 смс-информирование
*25.07.2016: 4000 хостинг (hc.ru)
*24.05.2016: 5000 хостинг (hc.ru)
*23.03.2016: 5000 хостинг (hc.ru)
*25.02.2016: 4000 хостинг (hc.ru)
*04.02.2016: 2000 год хостинга 1 Тб на яндекс.диске для папки aster
*04.02.2016: 2500 оплата апгрейда форума, первый этап
*24.01.2016: 5000 хостинг (hc.ru)
*22.12.2015: 4000 хостинг (hc.ru)
*26.11.2015: 5000 хостинг (hc.ru)
*27.10.2015: 5000 хостинг (hc.ru)
*24.09.2015: 5000 хостинг (hc.ru)
*29.08.2015: 2000 починка блога (сломалась поддержка кириллицы при обновлении)
*27.08.2015: 5000 хостинг (hc.ru)
*23.08.2015: 3000 дизайн баннеров
*25.07.2015: 4000 хостинг (hc.ru)
*20.07.2015: 240 смс-информирование

[[Категория:Служебные]]

Данные OpenStreetMap в формате shape-файлов

2018-04-18T22:01:42Z

Максим Дубинин: /* Преимущества */

{{Статья|Опубликована|osmshp}}
{{Аннотация|Здесь можно скачать данные OpenStreetMap на любой регион в разных форматах.}}

== Общая информация ==

В рамках данного проекта, из данных OpenStreetMap и других открытых источников создаются обновляемые наборы слоев по любой точке мира, включая страны бывшего СССР и все регионы РФ. Данные наборы слоев доступны в форматах ESRI Shapefile, Mapinfo TAB, GeoJSON, что позволяет использовать их практически в любой ГИС.
{{NextGIS}}

Если вам нужны данные в формате XML или PBF, то их можно найти на [http://gis-lab.info/projects/osm_dump/ другой] странице.

После 12 сентября 2012 года, данные публикуются под лицензией [http://www.opendatacommons.org/licenses/odbl ODBL]. Подробнее о ваших правах и обязанностях можно почитать в [http://wiki.openstreetmap.org/wiki/Legal_FAQ Legal FAQ].

Данные для загрузки: http://data.nextgis.com/osmshp/ 
Cписок слоев: http://data.nextgis.com/osmshp/layer/

Система координат данных: WGS84 
Проекция: широта/долгота

== Преимущества ==
Эти данные обладают рядом особенностей и преимуществ относительно исходных данных в формате OSM XML:

* Данные представлены в распространенных ГИС форматах
* Добавлены готовые проекты для ГИС, благодаря которым данные можно немедленно начинать использовать в ГИС
* Данные разбиты по слоям, у каждого слоя своя четкая структура атрибутов ([https://docs.google.com/spreadsheets/d/1F83dtRH8c7O83E55ox3Kfh8Ibbh2TFL70nF5Iw_33d0/edit#gid=1897387324 посмотреть])
* Объекты слоёв административно-территориального деления и населенных пунктов имеют коды ОКТМО (источник - Росстат)
* Добавлены стили, условные обозначения

Дополнительные слои, отсутствующие в таком виде в исходных данных:
* Административно-территориальное деление в форме обрезанной по границе суши (исходный вариант так же сохранён)
* Береговые линии, границы суши и водной поверхности

== Условия распространения ==
Данные на июль 2017 г. можно скачать бесплатно (кроме РФ целиком). Для скачивания необходимо получить [http://my.nextgis.com NextGIS ID].

Можно так же заказать:
* разовую платную выгрузку, стоимость от 200р.
* регулярную выгрузку с периодичностью раз в день, каждую неделю или месяц (до полугода).
* выгрузку по произвольному контуру, стоимость от 300р.

Актуальность исходной базы данных - текущий день.

Оплата осуществляется банковской картой или Яндекс.Деньги. Если заказ сделан в выходной день, выгрузка производится в ближайший рабочий. Если заказ сделан в конце рабочего дня, результат могут прислать на следующий. Если вам не прислали результат в течение 2 рабочих дней, оплата будет возвращена.

Выгрузка данных выполняется сразу после заказа. Готовность данных - в среднем 10-15 минут.

== Замечания по форматам выгрузки ==

* ESRI Shape ([https://nextgis.com/data/examples/shapefile-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs) или проект для ArcGIS (файл *.mxd)
* GeoJSON ([https://nextgis.com/data/examples/geojson-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs)
* ESRI Geodatabase ([https://nextgis.com/data/examples/geodatabase-arcgis-example.7z пример]))
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для ArcGIS (файл *.mxd)
* Mapinfo TAB ([https://nextgis.com/data/examples/mapinfo-example.7z пример])
** Кодировка атрибутивной информации - Windows 1251 (CP1251)
** В комплект включены проект для Mapinfо (файл *.wor)
* CSV ([https://nextgis.com/data/examples/csv-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
* PDF ([https://nextgis.com/data/examples/pdf-example.pdf пример])
** Максимальный размер листа - 3 на 3 метра.
** Для большей детализации рекомендуется заказывать не целый регион, а конкретную область с помощью инструмента [http://data.nextgis.com/osmshp/custom_region указания области].

Возможна выгрузка и в другие форматы по запросу.

== Проекты для QGIS ==

В архивы с данными в формате ESRI Shape и GeoJSON включаются несколько вариантов проектов для QGIS.

{|class="wikitable" width="100%"
| width="0%" | [[Файл:OSM shp qgis mapnik.png|300px]]
| width="100%" style="vertical-align:top" |
==== qgis_mapnik ====

Этот вариант проектов был разработан [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8933 old_Bibigon] в рамках [http://gis-lab.info/qa/osm-qgis-projects.html отдельного проекта].

|-
| [[Файл:OSM shp qgis mapnik lc.png|300px]]
| width="100%" style="vertical-align:top" |

==== qgis_mapnik_lowcontrast ====

Вариант оформления предоставленный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8650 sanyasi] и дополненный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=11960 gisuser] на форуме. Аналогично предыдущему варианту оформления - тоже вариации на тему слоя Mapnik в OpenStreetMap, но с пониженной контрастностью, что может быть удобно в случае использования в качестве базовой карты.
|}

Если у вас есть предложения по корректировке этих вариантов, или вы хотите предложить свой вариант, то сделать это можно на форуме в теме-обсуждении данной статьи.

== Ссылки ==
* [http://nextgis.ru/contact Обратная связь]
* [http://gis-lab.info/projects/osm_dump/ Данные OpenStreetMap в форматах XML и PBF]

Данные OpenStreetMap в формате shape-файлов

2018-04-18T21:59:46Z

Максим Дубинин: /* Преимущества */

{{Статья|Опубликована|osmshp}}
{{Аннотация|Здесь можно скачать данные OpenStreetMap на любой регион в разных форматах.}}

== Общая информация ==

В рамках данного проекта, из данных OpenStreetMap и других открытых источников создаются обновляемые наборы слоев по любой точке мира, включая страны бывшего СССР и все регионы РФ. Данные наборы слоев доступны в форматах ESRI Shapefile, Mapinfo TAB, GeoJSON, что позволяет использовать их практически в любой ГИС.
{{NextGIS}}

Если вам нужны данные в формате XML или PBF, то их можно найти на [http://gis-lab.info/projects/osm_dump/ другой] странице.

После 12 сентября 2012 года, данные публикуются под лицензией [http://www.opendatacommons.org/licenses/odbl ODBL]. Подробнее о ваших правах и обязанностях можно почитать в [http://wiki.openstreetmap.org/wiki/Legal_FAQ Legal FAQ].

Данные для загрузки: http://data.nextgis.com/osmshp/ 
Cписок слоев: http://data.nextgis.com/osmshp/layer/

Система координат данных: WGS84 
Проекция: широта/долгота

== Преимущества ==
Эти данные обладают рядом особенностей и преимуществ относительно исходных данных в формате OSM XML:

* Данные представлены в распространенных ГИС форматах
* Добавлены готовые проекты для ГИС, благодаря которым данные можно немедленно начинать использовать в ГИС
* Данные разбиты по слоям, у каждого слоя своя четкая структура атрибутов ([https://docs.google.com/spreadsheets/d/1F83dtRH8c7O83E55ox3Kfh8Ibbh2TFL70nF5Iw_33d0/edit#gid=1897387324 посмотреть])
* Объекты слоёв административно-территориального деления и населенных пунктов имеют коды ОКТМО (источник - Росстат)
* Административно-территориальное деление также представлено в форме обрезанной по границе суши (исходный вариант так же сохранён)
* Добавлены стили, условные обозначения

== Условия распространения ==
Данные на июль 2017 г. можно скачать бесплатно (кроме РФ целиком). Для скачивания необходимо получить [http://my.nextgis.com NextGIS ID].

Можно так же заказать:
* разовую платную выгрузку, стоимость от 200р.
* регулярную выгрузку с периодичностью раз в день, каждую неделю или месяц (до полугода).
* выгрузку по произвольному контуру, стоимость от 300р.

Актуальность исходной базы данных - текущий день.

Оплата осуществляется банковской картой или Яндекс.Деньги. Если заказ сделан в выходной день, выгрузка производится в ближайший рабочий. Если заказ сделан в конце рабочего дня, результат могут прислать на следующий. Если вам не прислали результат в течение 2 рабочих дней, оплата будет возвращена.

Выгрузка данных выполняется сразу после заказа. Готовность данных - в среднем 10-15 минут.

== Замечания по форматам выгрузки ==

* ESRI Shape ([https://nextgis.com/data/examples/shapefile-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs) или проект для ArcGIS (файл *.mxd)
* GeoJSON ([https://nextgis.com/data/examples/geojson-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs)
* ESRI Geodatabase ([https://nextgis.com/data/examples/geodatabase-arcgis-example.7z пример]))
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для ArcGIS (файл *.mxd)
* Mapinfo TAB ([https://nextgis.com/data/examples/mapinfo-example.7z пример])
** Кодировка атрибутивной информации - Windows 1251 (CP1251)
** В комплект включены проект для Mapinfо (файл *.wor)
* CSV ([https://nextgis.com/data/examples/csv-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
* PDF ([https://nextgis.com/data/examples/pdf-example.pdf пример])
** Максимальный размер листа - 3 на 3 метра.
** Для большей детализации рекомендуется заказывать не целый регион, а конкретную область с помощью инструмента [http://data.nextgis.com/osmshp/custom_region указания области].

Возможна выгрузка и в другие форматы по запросу.

== Проекты для QGIS ==

В архивы с данными в формате ESRI Shape и GeoJSON включаются несколько вариантов проектов для QGIS.

{|class="wikitable" width="100%"
| width="0%" | [[Файл:OSM shp qgis mapnik.png|300px]]
| width="100%" style="vertical-align:top" |
==== qgis_mapnik ====

Этот вариант проектов был разработан [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8933 old_Bibigon] в рамках [http://gis-lab.info/qa/osm-qgis-projects.html отдельного проекта].

|-
| [[Файл:OSM shp qgis mapnik lc.png|300px]]
| width="100%" style="vertical-align:top" |

==== qgis_mapnik_lowcontrast ====

Вариант оформления предоставленный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8650 sanyasi] и дополненный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=11960 gisuser] на форуме. Аналогично предыдущему варианту оформления - тоже вариации на тему слоя Mapnik в OpenStreetMap, но с пониженной контрастностью, что может быть удобно в случае использования в качестве базовой карты.
|}

Если у вас есть предложения по корректировке этих вариантов, или вы хотите предложить свой вариант, то сделать это можно на форуме в теме-обсуждении данной статьи.

== Ссылки ==
* [http://nextgis.ru/contact Обратная связь]
* [http://gis-lab.info/projects/osm_dump/ Данные OpenStreetMap в форматах XML и PBF]

Данные OpenStreetMap в формате shape-файлов

2018-04-18T13:39:18Z

Максим Дубинин:

{{Статья|Опубликована|osmshp}}
{{Аннотация|Здесь можно скачать данные OpenStreetMap на любой регион в разных форматах.}}

== Общая информация ==

В рамках данного проекта, из данных OpenStreetMap и других открытых источников создаются обновляемые наборы слоев по любой точке мира, включая страны бывшего СССР и все регионы РФ. Данные наборы слоев доступны в форматах ESRI Shapefile, Mapinfo TAB, GeoJSON, что позволяет использовать их практически в любой ГИС.
{{NextGIS}}

Если вам нужны данные в формате XML или PBF, то их можно найти на [http://gis-lab.info/projects/osm_dump/ другой] странице.

После 12 сентября 2012 года, данные публикуются под лицензией [http://www.opendatacommons.org/licenses/odbl ODBL]. Подробнее о ваших правах и обязанностях можно почитать в [http://wiki.openstreetmap.org/wiki/Legal_FAQ Legal FAQ].

Данные для загрузки: http://data.nextgis.com/osmshp/ 
Cписок слоев: http://data.nextgis.com/osmshp/layer/

Система координат данных: WGS84 
Проекция: широта/долгота

== Преимущества ==
Эти данные обладают рядом особенностей и преимуществ относительно исходных данных в формате OSM XML:

* Это данные в распространенных ГИС форматах
* Данные разбиты по слоям
* У каждого слоя свои атрибуты
* Объекты слоёв административно-территориального деления и населенных пунктов имеют коды ОКТМО (источник - Росстат)
* Добавлены стили, условные обозначения
* Добавлены готовые проекты для ГИС

== Условия распространения ==
Данные на июль 2017 г. можно скачать бесплатно (кроме РФ целиком). Для скачивания необходимо получить [http://my.nextgis.com NextGIS ID].

Можно так же заказать:
* разовую платную выгрузку, стоимость от 200р.
* регулярную выгрузку с периодичностью раз в день, каждую неделю или месяц (до полугода).
* выгрузку по произвольному контуру, стоимость от 300р.

Актуальность исходной базы данных - текущий день.

Оплата осуществляется банковской картой или Яндекс.Деньги. Если заказ сделан в выходной день, выгрузка производится в ближайший рабочий. Если заказ сделан в конце рабочего дня, результат могут прислать на следующий. Если вам не прислали результат в течение 2 рабочих дней, оплата будет возвращена.

Выгрузка данных выполняется сразу после заказа. Готовность данных - в среднем 10-15 минут.

== Замечания по форматам выгрузки ==

* ESRI Shape ([https://nextgis.com/data/examples/shapefile-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs) или проект для ArcGIS (файл *.mxd)
* GeoJSON ([https://nextgis.com/data/examples/geojson-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для QGIS (файл *.qgs)
* ESRI Geodatabase ([https://nextgis.com/data/examples/geodatabase-arcgis-example.7z пример]))
** Кодировка атрибутивной информации - UTF-8
** В комплект включен проект для ArcGIS (файл *.mxd)
* Mapinfo TAB ([https://nextgis.com/data/examples/mapinfo-example.7z пример])
** Кодировка атрибутивной информации - Windows 1251 (CP1251)
** В комплект включены проект для Mapinfо (файл *.wor)
* CSV ([https://nextgis.com/data/examples/csv-example.7z пример])
** Кодировка атрибутивной информации - UTF-8
* PDF ([https://nextgis.com/data/examples/pdf-example.pdf пример])
** Максимальный размер листа - 3 на 3 метра.
** Для большей детализации рекомендуется заказывать не целый регион, а конкретную область с помощью инструмента [http://data.nextgis.com/osmshp/custom_region указания области].

Возможна выгрузка и в другие форматы по запросу.

== Проекты для QGIS ==

В архивы с данными в формате ESRI Shape и GeoJSON включаются несколько вариантов проектов для QGIS.

{|class="wikitable" width="100%"
| width="0%" | [[Файл:OSM shp qgis mapnik.png|300px]]
| width="100%" style="vertical-align:top" |
==== qgis_mapnik ====

Этот вариант проектов был разработан [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8933 old_Bibigon] в рамках [http://gis-lab.info/qa/osm-qgis-projects.html отдельного проекта].

|-
| [[Файл:OSM shp qgis mapnik lc.png|300px]]
| width="100%" style="vertical-align:top" |

==== qgis_mapnik_lowcontrast ====

Вариант оформления предоставленный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=8650 sanyasi] и дополненный пользователем [http://gis-lab.info/forum/memberlist.php?mode=viewprofile&u=11960 gisuser] на форуме. Аналогично предыдущему варианту оформления - тоже вариации на тему слоя Mapnik в OpenStreetMap, но с пониженной контрастностью, что может быть удобно в случае использования в качестве базовой карты.
|}

Если у вас есть предложения по корректировке этих вариантов, или вы хотите предложить свой вариант, то сделать это можно на форуме в теме-обсуждении данной статьи.

== Ссылки ==
* [http://nextgis.ru/contact Обратная связь]
* [http://gis-lab.info/projects/osm_dump/ Данные OpenStreetMap в форматах XML и PBF]

Расширение функциональности MapInfo на основе SmartPanels

2018-04-07T08:16:00Z

Максим Дубинин:

{{Статья|Опубликована|mapinfo-smartpanels}}
{{Аннотация|В данной статье речь пойдет об утилите SmartPanels для MapInfo, позволяющей достаточно просто интегрировать в MapInfo пользовательские специализированные приложения или приложения, расширяющие функциональность MapInfo.}}

==Описание==
В данной статье речь пойдет об утилите SmartPanels для MapInfo, позволяющей достаточно просто интегрировать в MapInfo пользовательские специализированные приложения или приложения, расширяющие функциональность MapInfo.

MapInfo имеет штатный язык программирования MapBasic, который позволяет автоматизировать практически все операции MapInfo. Всем хорош MapBasic, пока дело не доходит до создания с его помощью достаточно развитого пользовательского интерфейса. Тут возникает сразу несколько проблем:
* большая трудоемкость при создании и размещении элементов управления.<ref>В качестве примера можно привести ситуацию, когда в уже созданный диалог с большим количеством элементов управления нужно внести новые или изменить расположение существующих элементов управления.</ref>
* невозможность создавать немодальные окна;
* ограниченный набор элементов управления состоящий всего из 15 элементов.


С другой стороны, существует достаточно широкий класс задач, где развитой пользовательский интерфейс просто необходим. В качестве характерного примера можно привести тот же кадастровый учет, в котором необходимо отображать большое количество разнообразной информации, привязанной к одному участку.

Разработчиками предпринимались различные попытки обойти ограничения MapBasic. К ним можно отнести создание визуальных редакторов диалогов (своеобразные IDE), реализация функционала MapInfo в среде MSAccess с использованием MapX, динамическое связывание MapInfo и MSAccess посредством механизма DDE и т.п. Но какого-то универсального инструмента не существовало. Точнее он был. Однако для его использования необходимо было быть специалистом в C++, а точнее хорошо разбираться в тонкостях создания библиотек DLL (разработка, подключение их в MapBasic, вызов из них внешних подпрограммы и пр.) что, согласитесь, для прикладного программиста MapBasic не является характерным.

Ситуация изменилась коренным образом с выходом утилиты SmartPanels, разработанной специалистами компании ЭСТИ МАП.<ref> Официальный дистрибьютер MapInfo в России и СНГ</ref> Разработка данной утилиты стала возможной, благодаря переходу MapInfo начиная с версии 9.5 на платформу Microsoft .NET. С помощью этой утилиты стало возможным быстро подключать (регистрировать) к MapInfo пользовательские приложения, созданные на любом из языков семейства .NET. Утилита является, как бы, посредником между MapInfo и внешним приложением. При этом для разработки прикладного решения не требуется установленного компилятора MapBasic. Все что нужно–это установленная утилита SmartPanel и среда разработки MSVisualStudio.
====Примечания====
<references/>

==Примеры реализации==
Для наглядного представления о том, что представляет из себя SmartPanels, приведем примеры четырех плагинов. Первые два поставляются в качестве примеров вместе со SmartPanels и представляет собой плагины общего назначения, т.е. расширяющие стандартные возможности MapInfo. Первый из них - «Вкладки» вносит полезное дополнение к пользовательскому интерфейсу MapInfo. Так при развернутых окнах, в верхней части окон появляются вкладки с именами документов и переключения между ними становится намного удобнее.<ref>в MapInfo 16 версии это уже штатный функционал, но и режим отображения окон документов только один – "во все окно".</ref> Вторым дополнением к пользовательскому интерфейсу <ref>Данная возможность обеспечивается непосредственно SmartPanels</ref> является возможность «прикреплять» панели к одной из сторон окна MapInfo с помощью соответствующих инструментов, отображающихся при перемещении панелей. Если ко одной стороне будет прикреплено несколько панелей, то они также будут отображаться в виде вкладок.

[[Изображение:SmartPanels_img1.png|800px|center|frame|Вкладки документов и прикрепляемые панели.]]

Второй плагин – «Информация». Его задачей является отображение и редактирование информации по выбранному объекту в настраиваемой пользователем форме. При первом запуске он ничем не отличается от стандартного окна информации MapInfo.

[[Изображение:SmartPanels_img2.png|800px|center|frame|Вид окна "Информация" SmaprtPanels.]]

Но вызвав с помощью правой клавиши мыши контекстное меню, можно открыть визуальный редактор экранной формы и настроить представление информации в соответствии со своими потребностями.
Кроме этого он позволяет определять для полей списки предопределенных значений, т.е. определять словари. Причем эти списки можно создавать как непосредственно в дизайнере форм, так и в обычных таблицах MapInfo или таблицах базы данных, указывая при этом ссылку на них в дизайнере форм. Это значительно ускоряет ввод и редактирование данных, помогая избежать многих ошибок.

[[Изображение:SmartPanels_img3.png|800px|center|frame|Настройка окна "Информация"]]

Остальные два плагина были созданы под конкретные проекты. Первый из них предназначен для отображения данных переписи населения.
Задача состояла в представлении данных переписи населения в различных разрезах. В MapInfo имелся слой с населенными пунктами и имелась база данных с результатами переписи. Плагин выполняет две основные задачи:
*получает от MapInfo уникальный код населенного пункта (код ОКАТО), осуществлял выборку данных из БД их интерпретацию и представление.
*передает MapInfo команду на позиционирование окна карты в центре с выбранным пользователем в иерархическом списке населенным пунктом.

[[Изображение:SmartPanels_img4.png|800px|center|frame|Выбранный населенный пункт (ОКАТО) в центре карты.]]

При использовании специального инструмента «Информация» при выборе населенного пункта на карте происходит отображение информации в различных разрезах

[[Изображение:SmartPanels_img5.png|800px|center|frame|Представление результатов переписи населения в иерархическом списке.]]

[[Изображение:SmartPanels_img6.png|800px|center|frame|Представление результатов переписи населения в виде круговой диаграммы.]]

[[Изображение:SmartPanels_img7.png|800px|center|frame|Представление результатов переписи населения в виде столбчатой диаграммы.]]

И еще один пример. В данном случае был разработан плагин для отображения панорам Yandex и Google при выборе точки на карте и использования поискового сервиса Yandex.

[[Изображение:SmartPanels_img8.png|800px|center|frame|Использование панорамы Yandex.]]

[[Изображение:SmartPanels_img9.png|800px|center|frame|Использование панорамы Google.]]

===Примечания===
<references/>

==Пример создания плагина для SmartPanel==
Особый интерес, конечно же, представляет возможность создания собственных плагинов. В качестве примера рассмотрим создание приложения, которое будет не просто приветствовать этот мир, но реализовать полезную функцию – клонировать окно отчета. В стандартном MapIinfo такая возможность отсутствует, в отличие от окна карты.
Для начала необходимо скачать и установить SmartPanel. Утилита является абсолютно бесплатной и скачать ее можно по ссылке [[Медиа:SPInstaller32.zip‎]] для 32 разрядной версии и [[Медиа:SPInstaller64.zip]] для 64 разрядной. В терминах SmartPanel все подключаемые пользовательские приложения носят названия плагинов. В качестве примеров вместе с самой утилитой распространяются также бесплатные плагины – «Вкладки» и «Расширенная панель информации». Подробное описание этих плагинов можно найти по той же ссылке.
После установки в основном меню MapInfo появится новый пункт – SmartPanel. Данный пункт содержит два элемента: Настройки и Панели. Пункт «Настройки» предназначен для основных действий с плагинами по их настройке, удалению и добавлению в SmartPanels. Второй пункт отвечает за отображение списка всех зарегистрированных панелей, позволяет их показывать или скрывать. 

[[Image:SmartPanels_img10a.png|450px|center|Меню SmartPanels в Mapinfo до 16 версии]]
[[Image:SmartPanels_img10b.png|450px|center|Меню SmartPanels в Mapinfo в 16 версии]]

Теперь все готово для добавления своего собственного плагина. Для его создания необходимо установить MS Visual Studio. Вполне подойдет и бесплатная версия этой IDE – Community [https://www.visualstudio.com/ru/downloads]. После запуска загрузчика рекомендуется выбрать пункт «Разработка классических приложений .NET»

[[Image:SmartPanels_img11.png|center|frame|Настройка загрузчика MSVisualStudio]]

После установки MS VisualStudio можно приступать непосредственно к созданию нашего проекта:'' Файл – Проект – Cоздать – Библиотека классов (.NET Framework)''

[[Image:SmartPanels_img12.png|center|frame|Общий вид меню «Создание проекта».]]

В созданном проекте нужно указать ссылку на библиотеку SmarPanels.dll. Данная библиотека находится в корневом каталоге установки SmartPanel. Если при установке утилиты были использованы настройки по умолчанию, то эта библиотека будет находится в каталоге C:\Program Files\ESTI MAP\SmartPanels. Для создания ссылки необходимо в панели «Обозреватель решения» вызвать диалог добавления ссылки путем нажатия правой клавиши мыши на пункте «Ссылки» и далее «Добавить ссылку…». После инициализации проекта был создан класс по умолчанию Class1. Название класса можно оставить без изменений, а можно переименовать, но в любом случае этот класс должен реализовывать интерфейс SPNet.IMiPlugin. Это и есть предопределенный шаблон плагина. В контексте этого интерфейса наш класс должен реализовывать пять свойств и три метода.
Свойства:
*PluginName - название плагина, которое будет отображаться в списке установленных плагинов.
*PluginDescription - описание плагина.
*PluginUID – глобально уникальный идентификатор плагина.
*PluginVersion - версия плагина.
*RestoredPanels - список панелей, местоположение которых SmartPanels будет автоматически сохранять и восстанавливать при запуске.
Методы:
*Init - производит инициализацию плагина.
*LoadPluginSettings - загружает состояние плагина.
*Dispose - освобождение используемых ресурсов при завершении приложения.

Для краткости не будем здесь приводить подробное описание этих свойств и методов. В состав дистрибутива SamrtPanel входит полная документация разработчика, где все подробно описано. Для генерации PluginUID можно воспользоваться соответствующим пунктом меню VisualStudio:'' Средства - Создать GUID ''. В конечном итоге проект должен выглядеть аналогично снимку ниже.

[[Image:SmartPanels_img13.png|center|frame|Класс CloneLayout.]]

На данном этапе плагин создан и готов к добавлению в SmartPanels, но он еще не имеет визуального представления. Добавим в проект панель для размещения на ней элементов управления:'' Проект – Добавить пользовательский элемент управления… ''. В качестве имени задайте LayoutControl.cs. На созданном элементе управления разместите группирующую панель (Panel в разделе Контейнеры панели Элементы управления), элемент управления «ListBox» и одну кнопку. (раздел Стандартные элементы управления). Новый элемент управления должен выглядеть примерно следующим образом:

[[Image:SmartPanels_img15.png|center|frame|Пользовательские элементы управления.]]

Экранная форма плагина создана. Теперь надо сообщить об этом SmartPanels. В классе CloneLayout.cs необходимо добавить следующий код:

<syntaxhighlight lang="java">
//Создаем новый объект – панель SmartPanels и его идентификатор
// Guid необходимо сформировать с помощью меню Средства-Создать GUID
MiDockPanel pnl = null;
Guid pnlID = new Guid("{77BCA491-65DA-4C40-A1AE-B792C77DFDAC}");
//Заменияем public List<Guid> RestoredPanels{get{return null;}} на
public List<Guid> RestoredPanels
{
get
{List<Guid> l = new List<Guid>();
l.Add(pnlID);
return l;}
}
//Реализовываем метод Init(), который вызывается при инициализации плагина
public void Init()
{
//Событие MiHost.Inited вызывается после инициализации всех плагинов.
MiHost.Inited += MiHost_Inited;
//Добавляем панель
pnl = MiHost.AddDockPanel("Клонирование отчета", pnlID);
//Создаем экземпляр нашего элемента управления
LayoutControl ctrl = new LayoutControl();
ctrl.Dock = DockStyle.Fill;
//Созданная панель должна иметь вид нашего элемента управления
pnl.Controls.Add(ctrl);
}
//Добавляем метод Init(), который будет выполняться при инициализации плагина
//Подробное описание команд см. Документацию
void MiHost_Inited(object sender, EventArgs e)
{try
{
MiHost.AddPanelToMenu("Клонирование отчета",SPNet.MbCommands.MenuItemShortcut
.CreateShortcut(SPNet.MbCommands.MenuItemShortcut.ModifierKey.Ctrl,(byte)'T'), "Клонирование отчета",pnl);
}
catch (Exception except)
{ MessageBox.Show(MiHost.MainHandle,
"Произошла ошибка при инициализации плагина:" + except.Message,
"Клонирование отчета",
MessageBoxButtons.OK,
MessageBoxIcon.Error);
}
}
</syntaxhighlight>

После сделанных изменения класс CloneLayout.cs должен выглядеть следующим образом:

[[Image:SmartPanels_img16.png|center|frame|Измененный класс CloneLayout.]]

На этом основная работа по созданию и регистрации нового плагина выполнена. На этом этапе возможно посмотреть, как будет выглядеть плагин в MapInfo. Для этого необходимо собрать решение:''Сборка – Собрать решение F6'' и дождаться окончания сборки. По умолчанию готовая сборка в виде файла CloneLayout.dll будет находится в каталоге \bin\Debug корневого каталога проекта. Далее необходимо его загрузить в SmartPanels в Mapinfo: меню ''SmartPanels – Настройки – вкладка Плагины - кнопка Добавить''

[[Image:SmartPanels_img17.png|750px|center|Добавление созданного плагина]]

К сожалению, существует особенность, при которой необходимо перезагружать MapInfo (или SamrtPanels) после установки нового плагина. После перезагрузки в списке установленных плагинов появится новый плагин.

[[Image:SmartPanels_img18.png|center|frame|Новый плагин в MapInfo]]

Теперь реализуем взаимодействие с MapInfo. В класс LayoutControl.cs добавим реализацию метода, который будет опрашивать MapInfo на предмет открытых окон и их типа. Если тип окна «Отчет», то имя окна и его ID будет добавляться в список формы. Для этого будем использовать метод <code>MiHost.MB.GetWindows</code>, возвращающий коллекцию всех открытых в MapInfo окон. Этот метод может получать в качестве параметра фильтр <code>MbCommands.MiWindowInfo.WindowTypeFilter</code>, который будет указывать SmartPanels возвращать окна только определенного типа.
<syntaxhighlight lang="java">
private void fillList() {
listLayout.Items.Clear();
foreach (var lw in MiHost.MB.GetWindows(MbCommands.MiWindowInfo.WindowTypeFilter.WIN_LAYOUT)) {
listLayout.Items.Add(string.Format("{0}#{1}", lw.WindowID, lw.WindowName));
}
//Пусть последний элемент списка всегда будет выделенным
if (listLayout.Items.Count>0) {
listLayout.SelectedIndex = listLayout.Items.Count - 1;
}
}
</syntaxhighlight>
Хотя SmartPanels и имеет развитый программный интерфейс (API), в большинстве случаев, можно обойтись всего двумя командами - <code>MiHost.MB.EvalCommand</code> и <code>MiHost.MB.EvalCommandX</code>, где X обозначает тип возвращаемого значения (EvalCommandS – для String, EvalCommandI – для Integer и т.п.), первая команда просто выполняет команду, а вторая еще и возвращает полученное значение. На следующем шаге создадим метод, вызываемый при нажатии на кнопку «Клонировать». Данный метод будет определять активный элемент списка и засылать в MapInfo команду на дублирование этого окна.
<syntaxhighlight lang="java">
private void btnClone_Click(object sender, EventArgs e)
{
if (listLayout.SelectedItem != null)
{
string item = listLayout.SelectedItem.ToString();
string idWin = item.Substring(item.IndexOf("#")+1);
//Засылаем в MI команду на клонирование окна
MiHost.MB.EvalCommand("Run Command WindowInfo("+idWin+", 15)");
}
}
</syntaxhighlight>
Осталось только определить, как будет обновляться список отчетов. Для этого попросим MapInfo присылать сообщение при открытии или закрытии окон. При получении такого сообщения будем обновлять список отчетов.
<syntaxhighlight lang="java">
public LayoutControl()
{
InitializeComponent();
//Будем следить за новыми окнами в MI
MiHost.AddDocumentWindow += MI_NewWin;
//и за закрытыми окнами в MI
MiHost.RemoveDocumentWindow+= MI_RemoveWin;
}
void MI_NewWin(object sender, EventArgs e) {
fillList();
}
void MI_ RemoveWin(object sender, EventArgs e)
{
fillList();
}
</syntaxhighlight>
После внесенных изменений класс должен выглядеть следующим образом:

[[Image:SmartPanels_img19.png|center|frame|Измененный класс LayoutControl]]

После сборки решения открываем MapInfo и проверяем результат.

[[Image:SmartPanels_img20.png|center|frame|Плагин в действии]]

==Заключение==

SmartPanels существенно облегчает жизнь разработчикам приложений MapInfo и позволяет осуществлять проекты любой сложности. Это могут быть как проекты, реализующие конкретные прикладные задачи, так и проекты, расширяющие функциональные возможности MapInfo. В качестве примера последнего можно привести разрабатываемый в настоящий момент плагин, который является аналогом Редактора моделей анализа в QGIS и ModelBuidler в ArcGIS.