Формат ASCIIGRID с примерами

Эта страница является черновиком статьи.

Формат ASCIIGRID это форма представления регулярно-сетевых данных (растров). Исходно этот формат использовался как обменный в ГИС ARC/INFO (файлы в этом формате создавались командой GRIDASCII), сейчас его основная часть поддерживается всеми основными ГИС (QGIS, ArcGIS). Файл в формате ASCIIGRID представляет собой текстовый файл в котором содержимое разделяется на две части: 1) заголовок, где записывается общая информация о растре, его размеры, географическая привязка и т.д., 2) собственно данные, где записываются данные.

Данные в этом формате удобны прежде всего удобством быстрого визуального анализа содержимого: много-мало данных, какие значения, есть ли nodata и т.д.. Разумеется, в этом формате не рекомендуется хранить большие растры.

Поддерживается чтение и запись в этот формат с помощью GDAL.

Заголовок

ncols - количество колонок (целочисленное);
nrows - количество рядов (целочисленное);
xllcorner - х-координата левого нижнего угла, долгота (десятичные градусы DD.D(D), DD:MM:SS не поддерживается) или смещение на восток (метры) (целочисленное или с плавающей точкой)
yllcorner - y-координата левого нижнего угла, широта или смещение на север
cellsize - линейный размер квадратной ячейки (целочисленное или с плавающей точкой)
nodata_value - значение "нет данных", этот параметр не обязателен, он нужен для того, чтобы программное обеспечение автоматически делало все ячейки с этим значением прозрачными (целочисленное или с плавающей точкой). Часто используемое значение -9999, нецифровые значения, например "NA" использовать нельзя

Данные

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

Примеры

Все примеры можно копировать и вставлять в текстовые файлы с расширением *.asc, они должны работать. Чтобы получить такие же иллюстрации в QGIS, нужно: 1. добавить растр, 2. зайти в свойства\Стиль и изменить Цветовая карта на Псевдоцвета, 3. зайти в свойства\Прозрачность и удалить прозрачный пиксел со значением -9999, а также заменить в значении no data -9999 на 0.

Стандартный заголовок

Это же изображение в табличном виде:

-9999	-9999	-2080	-1909	-9999	-601	-9999	-1362
-510	-387	-350	-182	-9999	-9999	-9999	27
13	-92	-9999	-33	54	28	-182	-750
59	-165	18	167	-182	-9999	-608	666
-62	-9999	-805	-9999	-9999	-9999	-9999	-497
-9999	-9999	332	1240	-9999	-9999	-9999	-9999

И непосредственно в формате ASCIIGRID

ncols   8
nrows        6
xllcorner    -40.5
yllcorner    25.25
cellsize     10
NODATA_value  -9999
-9999 -9999 -2080 -1909 -9999 -601 -9999 -1362
-510 -387 350 -182 -9999 -9999 -9999 27 
13 -92 -9999 -33 54 28 -182 -750
59 -165 18 167 -182 -9999 -608 666 
-62 -9999 -805 -9999 -9999 -9999 -9999 -497 
-9999 -9999 332 1240 -9999 -9999 -9999 -9999

Неквадратная ячейка

В примере приведенном выше, ячейка растра квадратная (т.е. высота равна широте). Если необходима прямоугольная ячейка, где высота не равна широте, вместо параметра cellsize можно использовать DX и DY, например:

ncols        8
nrows        6
xllcorner    -40.5
yllcorner    25.25
dx           11.6
dy           5.025
NODATA_value  -9999
-9999 -9999 -2080 -1909 -9999 -601 -9999 -1362
-510 -387 350 -182 -9999 -9999 -9999 27 
13 -92 -9999 -33 54 28 -182 -750
59 -165 18 167 -182 -9999 -608 666 
-62 -9999 -805 -9999 -9999 -9999 -9999 -497 
-9999 -9999 332 1240 -9999 -9999 -9999 -9999

ASCIIGRID файлы с параметрами DX и DY вместо CELLSIZE не поддерживаются многими пакетами GIS (в т.ч. ArcGIS), но поддерживаются QGIS.

Вот как будет выглядеть файл выше, где cellsize заменен на DX, DY:

Привязка по центру ячейки, а не углу

Альтернативой привязки по углу, т.е. использования xllcorner и yllcorner является привязка по центру ячейки, в этом случае нужно использовать xllcenter и yllcenter. Результатом, при неизменных координатах будет сдвиг растра на половину ячейки вниз и влево (по диагонали).

ncols        8
nrows        6
xllcenter    -40.5
yllcenter    25.25
cellsize     10
NODATA_value  -9999
-9999 -9999 -2080 -1909 -9999 -601 -9999 -1362
-510 -387 350 -182 -9999 -9999 -9999 27 
13 -92 -9999 -33 54 28 -182 -750
59 -165 18 167 -182 -9999 -608 666 
-62 -9999 -805 -9999 -9999 -9999 -9999 -497 
-9999 -9999 332 1240 -9999 -9999 -9999 -9999

Если взять один и тот же растр и заменить xllcenter,yllcenter на xllcorner,yllcorner, а сверху наложить сетку границ пикселей с первого растра, то это будет выглядеть так:

Практика

Таким образом, чтобы создать растр размером m*n не используя никаких инструментов кроме текстового редактора и табличного процессора, нужно сделать следующее.

1. Создать таблицу нужного размера в табличном процессоре (OpenOffice Calc, MS Excel и т.д.) или любом редакторе. Скопировать и вставить в текстовый документ. Заменить табуляции на пробелы. 2. Прописать заголовок файла в описанном выше формате. Еще один пример для матрицы 360*180, представляющий массив с разрешением 1 градус на весь мир:

ncols        360
nrows        180
xllcorner    -180.000000000000
yllcorner    -90.000000000000
cellsize     1.000000000000
nodata_value 0
12 321 1 12 3 4 1 23 4 1 234 1234 1 2 3 4 2 1 2 3 ...

3. Поменять/задать расширение текстового файла на asc. 4. Полученный файл можно открыть в QGIS и экспортировать в любой другой растровый формат или работать прямо с ним.