R: RGDAL написание GeoJSON
У меня проблема с rgdal, который сбрасывает информацию о проекции при записи файлов geoJSON.
library(rgdal)
inputJSON<- readOGR("test.geojson", "OGRGeoJSON") # works!
Если я наберу summary(inputJSON) я получаю следующий вывод:
Object of class SpatialPointsDataFrame
Coordinates:
min max
coords.x1 0 499690.8
coords.x2 0 321771.2
Is projected: TRUE
proj4string :
[+proj=tmerc +lat_0=0 +lon_0=13.33333333333333 +k=1 +x_0=450000 +y_0=-5000000
+ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs]
Это говорит мне о том, что файл был успешно прочитан в R и что он также получает информацию о проекции.
После некоторых манипуляций с данными я хочу сохранить свой data.frame (inputJSON) как новый файл geoJSON с помощью следующей команды:
writeOGR(inputJSON, "outTest.geojson", layer="inputJSON", driver="GeoJSON",check_exists = FALSE)
при этом также создается нужный файл geoJSON, но он не записывает в него информацию о проекции, а файл geoJSON без этой информации в значительной степени (для моих целей) бесполезен.
Для сравнения начало входного файла (test.geojson):
{
"type": "FeatureCollection",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31258" } },
"features": [
Созданный файл (outTest.geojson):
{
"type": "FeatureCollection",
"features": [
так вот:
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31258" } },
важная часть отсутствует.
Если я не ошибаюсь, моя версия rgdal
ргдаль версия 1.0-6
и должен быть самым новым. Ну, я также пробовал другие пакеты с доступной записи файлов geoJSON (geojsonio,leafletR) но ни один из них не рассматривает проекцию. Так чего мне не хватает? Есть ли возможность сделать эту работу? Я думаю, что драйвер GeoJSON ogr2ogr (на котором основан пакет) не обрабатывает проекцию вообще. При переносе шейп-файла в файл geoJSON через командную строку полученный файл geoJSON также пропускает информацию о проекции. Я действительно смущен!
Так чего мне не хватает?
У кого-нибудь есть (любое) решение?