Конвертация данных из CSV в SHP и обратно с OGR
по адресу http://gis-lab.info/qa/csv2shp.html
Рассмотрена конвертация данных из формата CSV в формат ESRI shape и обратно.
Введение
Довольно часто возникает задача конвертации файлов формата csv в shape-файлы. Иногда возникнает и обратная задача – конвертация shape-файлов в csv. Две данные операции и будут рассмотрены в этой статье.
Эта статья ориентирована на пользователя знакомого с командной строкой. В случае, если необходимо сконвертировать в shape один-другой файл csv, то возможно проще воспользоваться плагином QGIS «Delimited text», однако обратную конвертацию в QGIS выполнить на данный момент невозможно.
Стоит отметить, что до выхода версии GDAL 1.6.0 возможности конвертировать шейп-файл в формат csv c сохранением информации о геометрических объектах не было, поэтому советуем обновить вашу версию, если это необходимо.
Прежде чем переходить к преобразованию данных из одного формата в другой, рассмотрим, что же представлят из себя непосредственно формат CSV и более детально остановимся на особенностях драйвера OGR CSV. В этом описании используется информация с официального сайта библиотеки GDAL.
Формат CSV и особенности драйвера OGR CSV
OGR поддерживает чтение и запись табулированных данных (данных разделенных специальным символом-разделителем), хранящихся в текстовых файлах формата csv, являющегося распространенным форматом обмена табличными данными между различными программами. Файл данного типа можно легко создать вручную, например, с помошью любого текстового редактора.
Хотя на практике файл с разделителями может иметь любое расширение (csv,txt,…), драйвер OGR CSV поддерживает только файлы с расширением .csv. Данное требование объясняется тем, что при использовании утилит командной строки OGR не указывается формат входного файла и поэтому его формат определяется на основании расширения. Именем источника данных может быть либо единственный csv файл, либо каталог, содержащий несколько файлов. Чтобы каталог распознался как источник csv данных, необходимо, чтобы по крайней мере половина файлов в каталоге имела расширение .csv. Из каждого .csv файла в каталоге генерируется один слой (одна таблица).
Драйвер OGR CSV поддерживает как чтение, так и запись. Поскольку файл формата csv содержит текстовые строки переменной длины, то чтение выполняется последовательно (строка за строкой). Слой OGR CSV никогда не имеет описания системы координат. Драйвер OGR CSV возвращает тип всех полей как строковый, если нет файла описания типов полей (дополнительный файл с расширением .csvt).
Поддерживаемые типы полей - Integer, Real, String, Date (YYYY-MM-DD), Time (HH:MM:SS+nn) и DateTime (YYYY-MM-DD HH:MM:SS+nn). Как уже отмечалось, файл с описанием типов полей должен иметь расширение .csvt и такое же имя, как и файл csv. Типы полей задаются в двойных кавычках через запятую в одну строку (например, "Integer","String"). Также возможно непосредственно указать ширину и точность каждого поля (например, "Integer(5)","Real(10.7)","String(15)").
В файле csv каждая строка описывает отдельный объект в слое (таблице). Значения полей отделены друг от друга запятыми. В каждой строке должно быть минимум 2 поля. Строки отделены друг от друга разделителями в стиле DOS (CR/LF) или в стиле Unix (LF). Каждая строка должна иметь одинаковое количество полей. Начиная с версии GDAL 1.7.0, драйвер также поддерживает в качестве разделителей между полями – точку с запятой или табуляцию. Это автоопределение будет работать только в том случае, если в первой строке нет других потенциально возможных разделителей. В противном случае в качестве разделителя будет использована запятая.
Составные значения полей (например, содержащие запятые, кавычки) следует помещать в двойные кавычки. Если в значении поля содержатся кавычки, то их следует удвоить.
Драйвер пытается обработать первую строку файла как список имен всех полей. Однако, если одно из полей имеет числовое значение – то подразумевается, что в первой строке содержатся данные, описывающие отдельный объект, а в качестве имен полей используются имена field_1, field_2 и т.д.
Пример (employee.csv):
ID,Salary,Name,Comments 132,55000.0,John Walker,"The ""big"" cheese." 133,11000.0,Jane Lake,Cleaning Staff
Встречаются текстовые файлы, часто называемые csv, включающие данные не разделенные запятыми, а с фиксированной шириной поля. Данный драйвер не поддерживает подобный формат. Если у вас такие данные, преобразуйте их в простой .csv фомат, чтобы он стал понятен драйверу OGR CSV.
Конвертация данных из .csv в .shp
Одна из существенных возможностей OGR - извлечение пространственной информации из файла csv . В случае, если csv файл содержит поля в которых хранятся X и Y координаты точек, то конвертировать такой файл csv в shape-файл можно посредством использования промежуточного драйвера VRT. Описание этого драйвера и всех его параметров доступно на официальном сайте.
Рассмотрим следующий файл csv (test.csv):
Latitude,Longitude,Name 48.1,0.25,"First point" 49.2,1.1,"Second point" 47.5,0.75,"Third point"
Создадим соответствуюший файл vrt (test.vrt):
<OGRVRTDataSource>
<OGRVRTLayer name="test">
<LayerSRS>WGS84</LayerSRS>
<SrcDataSource>test.csv</SrcDataSource>
<GeometryType>wkbPoint</GeometryType>
<GeometryField encoding="PointFromColumns" x="Longitude" y="Latitude"/>
</OGRVRTLayer>
</OGRVRTDataSource>
и выполним команду:
ogr2ogr output test.vrt
Будет создана папка с названием output, а в нее добавлен файл test.shp.
Если же видоизменить команду:
ogr2ogr output.shp test.vrt
то shape-файл с именем output.shp будет создан в текущей папке.
В случае, если файл csv (test.csv) содержит поле в котором хранится описание геометрии объекта в формате WKT:
GeomField,Code "POLYGON((1 2,3 4, 2 1))",8 "POLYGON((11 21,31 41, 21 11))",9
то чтобы конвертировать его в shape-файл нужно отредактировать файл vrt (test.vrt):
<OGRVRTDataSource>
<OGRVRTLayer name="test">
<SrcDataSource>test.csv</SrcDataSource>
<GeometryType>wkbUnknown</GeometryType>
<GeometryField encoding="WKT" field="GeomField"/>
</OGRVRTLayer>
</OGRVRTDataSource>
То есть мы заменили encoding="PointFromColumns" на encoding="WKT" и указали поле, в котором содержится описание геометрии в WKT-формате field=" GeomField ". И после выполнения команды ogr2ogr как показано выше, мы получим шейп файл, соответствующей геометрии, в котором помимо остальных в качестве атрибутивного поля будет содержаться поле с WKT-представлением объектов. В связи с ограничением, накладываемым на размер атрибутивного текстового поля (255 символов) shape-файла, информация о WKT-представлении может оказаться обрезанной. Это не повлияет на саму геометрию, а только на ее копию в атрибутивном поле.
Если поле с WKT описанием геометрии в файле csv имеет имя WKT, например такой:
WKT,Code "POLYGON((1 2,3 4, 2 1))",8 "POLYGON((11 21,31 41, 21 11))",9
то можно обойтись без использования драйвера VRT. В этом случае конвертация из csv в shape выполняется командой:
ogr2ogr output.shp test.csv
Таким образом, что бы получить CSV файл, который легко перевести в shape-файл, нужно правильно заполнить одну из колонок csv файла WKT-описанием геометрии.
Конвертация данных из .shp в .csv
Конвертация данных из шейп файла также выполняется утилитой ogr2ogr. Опции создания csv файла:
- LINEFORMAT: По умолчанию, при создании нового csv файла в качестве разделителя строк используется символ в зависимости от текущей платформы (CR/LF под win32 или LF для других). Данный параметр может быть переопределен с использованием ключа LINEFORMAT, который может принимать значение CRLF (DOS формат) или LF (Unix формат).
- GEOMETRY (Начиная с версии GDAL 1.6.0): По умолчанию, представление геометрии не записывается в csv файл. Существует возможность экспортировать геометрию в представление WKT, определив ключ GEOMETRY=AS_WKT. Также возможно экспортировать точечную геометрию в X, Y, Z компоненты (отдельные поля .csv файла), определив ключ GEOMETRY=AS_XYZ, GEOMETRY=AS_XY или GEOMETRY=AS_YX. Поля с описанием геометрии будут располагаться перед атрибутивными полями
- CREATE_'CSVT'=YES/NO (Начиная с версии GDAL 1.7.0): создание связанного .csvt файла для описания типов полей. Значение по умолчанию NO.
- SEPARATOR=COMMA/SEMICOLON/TAB (Начиная с версии GDAL 1.7.0): Символ разделения полей. Значение по умолчанию : COMMA.
Пример конвертации шейп файла в csv:
ogr2ogr -f CSV output input.shp -lco GEOMETRY=AS_WKT -lco CREATE_CSVT=YES
В результате выполнения данной команды в текущей директории будет создан каталог output и в него помещен файл input.csv, содержащий поле с WKT-представлением геометрических объектов, а также файл input.csvt c описанием типов полей csv файла.
Отметим, что ключ -f используется при конвертировании шейп файла в csv в то время как при конвертации csv в шейп данный ключ не используется. Это происходит потому, что выходной формат по умолчанию «ESRI Shapefile», что равносильно строке -f «ESRI Shapefile».
Примеры программных реализаций
Perl - обрабатывает все файлы в папке (csv), создает vrt и генерирует shape-файлы каждый в своей папке. Условие - во всех файлах lat и long столбцы должны называтся одинаково.
Протестировано для точечных слоев на Suse 11.
Python - обрабатывает все файлы в папке (csv), создает vrt и генерирует shape-файлы.