Скрипты для подсчёта площадей полигонов в квадратных метрах в проекции Geographic.
Довольно часто покрытия хранятся в проекции Geographic, в геодезических координатах (единицы измерения — десятичные градусы). Если требуется подсчитать площадь полигонов в таком покрытии, чтобы она была не в квадратных градусах, а в квадратных метрах, можно использовать пакет, который можно скачать здесь (985 байт). В командной aml'ке (!control.aml) нужно связать переменную prjaml с файлом описания проекции, которую вы желаете использовать для расчёта площади. Можно использовать уже имеющиеся файлы — ddtoeqr.аml для расчётов в пределах России (коническая равновеликая проекция Альберса) или ddtoeqc.aml для расчётов в любой части мира (нормальная цилиндрическая равновеликая проекция Ламберта).
|
NB |
Обе aml'ки для расчётов в российской системе координат предназначены (система координат 1942 года). Само собой разумеется, если ваши покрытия находятся на другом эллипсоиде с другим датумом, их нельзя использовать будет. |
Затем следует связать переменную cover с названием покрытия, в котором нужно посчитать площади. Затем — запустить arc &run !control.aml. После выполнения пакета в атрибутивной таблице исходного покрытия появляется поле REAL_AREA, в котором будет занесена площадь в квадратных метрах.
Пакет работает неспеша, особенно если исходные покрытия большие. Зато с его помощью можно обрабатывать сколь угодно большое количество слоёв в пакетном режиме, встраивая его в более сложные пакеты обработки данных. Без участия оператора.
Следует обратить внимание, что в исходном покрытии изначально не должно быть поля REAL_AREA. Оно образуется в процессе работы пакета.
Для тех, кто работает в ArcView 3.x, полезным может оказаться скрипт
View.CalculateFeatureGeometry, (1 668 байт), разработанный компанией ESRI (http://www.esri.com/) для быстрого подсчёта геометрических характеристик объектов — площади, периметра и длины, в зависимости от проекции и характера локализации. Скрипт удобен для операций с единичными покрытиями. Работает быстро. Добавляет поля REAL_AREA и R_PERIMETER для полигональных объектов, и R_LENGTH для линейных объектов, если они ещё не присутствуют. Если поля с такими именами уже существуют, значения в них будут пересчитаны. В поля заносятся значения, рассчитанные в той проекции, в которой находится View.