R:rgdal编写GeoJSON

我对rgdal有一个问题，它在编写geoJSON文件时会丢失投影信息

library(rgdal)

inputJSON<- readOGR("test.geojson", "OGRGeoJSON") # works!

这告诉我该文件已成功读入R，并且它还获得了投影信息

在一些数据操作之后，我想使用以下命令将data.frame（inputJSON）保存为一个新的geoJSON文件：

writeOGR(inputJSON, "outTest.geojson", layer="inputJSON", driver="GeoJSON",check_exists = FALSE)

这也会生成所需的geoJSON文件，但它不会将投影信息写入其中，而且没有这些信息的geoJSON文件（在我看来）几乎毫无用处

要比较输入文件（test.geojson）的开头：

生成的文件（outTest.geojson）：

因此：

"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31258" } },

重要的部分丢失了

如果我没弄错的话，我的rgdal版本是

rgdal’版本1.0-6

而且应该是最新的。我还尝试了其他带有geoJSON文件编写器的包（

geojsonio，传单

），但没有一个包考虑了投影。 那么我错过了什么？有没有可能做到这一点？ 我认为ogr2ogr的GeoJSON驱动程序（包基于它）根本不处理投影。通过命令行将shapefile传输到geoJSON文件生成的geoJSON文件也会丢失投影信息。我真的很困惑

那么我错过了什么

---

有人有（任何）解决方案吗？

看看。罗宾向我们展示了如何处理这个GDAL缺陷。谢谢！我也找到了那一页。如果我有一个geoJSON文件作为输入文件，那么手动编写缺少的行的解决方案适用于这种情况，这样我就可以轻松地“复制粘贴”缺少的行。但是，当输入文件（例如）是一个shapefile时，该解决方案非常混乱，在最好的情况下，需要查找EPSG代码（对于R？）然后手动将该信息写入geoJOSN文件。这真的很奇怪，例如，如果我下载一个WFS层并使用命令行将其保存为geoJSON，则投影信息会写入文件中。那么一般来说应该是可能的吗？请注意，在RFC 7946（最新的格式提案）中，仅使用WS84，引用了互操作性问题。在创建选项中使用

RFC7946=YES

，您可以跳转并使用该模式。要添加，请声明：“替代坐标参考系的使用在[GJ1008]中有规定。”，但它已从此版本的规范中删除…一般来说，GeoJSON处理软件不需要访问坐标参考系数据库或网络访问坐标参考系转换参数。”请看。罗宾向我们展示了如何处理这个GDAL缺陷。谢谢！我也找到了那一页。如果我有一个geoJSON文件作为输入文件，那么手动编写缺少的行的解决方案适用于这种情况，这样我就可以轻松地“复制粘贴”缺少的行。但是，当输入文件（例如）是一个shapefile时，该解决方案非常混乱，在最好的情况下，需要查找EPSG代码（对于R？）然后手动将该信息写入geoJOSN文件。这真的很奇怪，例如，如果我下载一个WFS层并使用命令行将其保存为geoJSON，则投影信息会写入文件中。那么一般来说应该是可能的吗？请注意，在RFC 7946（最新的格式提案）中，仅使用WS84，引用了互操作性问题。在创建选项中使用

RFC7946=YES

，您可以跳转并使用该模式。要添加，请声明：“替代坐标参考系的使用在[GJ1008]中有规定。”，但已从此版本的规范中删除…一般来说，GeoJSON处理软件不应访问坐标参考系数据库或网络访问坐标参考系转换参数。”

{
"type": "FeatureCollection",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31258" } },

"features": [

{
"type": "FeatureCollection",

"features": [

"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:EPSG::31258" } },