[R-sig-Geo] WriteOGR to GeoJSON loses CRS
Frede Aakmann Tøgersen
frtog at vestas.com
Tue Aug 12 10:08:41 CEST 2014
Hi
I can get your example (SpatialPointsDataFrame) to work. On my windows box I have:
> library(rgdal)
Loading required package: sp
rgdal: version: 0.8-16, (SVN revision 498)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 1.11.0, released 2014/04/16
Path to GDAL shared files: c:/Programmer/R/R-3.1.0/library/rgdal/gdal
GDAL does not use iconv for recoding strings.
Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
Path to PROJ.4 shared files: c:/Programmer/R/R-3.1.0/library/rgdal/proj
However I cannot get this (using SpatialPolygonsDataFrame) to work:
> library(maptools)
Checking rgeos availability: TRUE
> xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))
> summary(xx)
Object of class SpatialPolygonsDataFrame
Coordinates:
min max
x -84.32385 -75.45698
y 33.88199 36.58965
Is projected: FALSE
proj4string : [+proj=longlat +ellps=clrk66]
Data attributes:
AREA PERIMETER CNTY_ CNTY_ID NAME
Min. :0.0420 Min. :0.999 Min. :1825 Min. :1825 Alamance : 1
> writeOGR(xx, "./.geojson", "xx", "GeoJSON")
No crs information at top of file.
Tried this giving no crs information in xx.geojson:
ogr2ogr -f GeoJSON c:\Users\frtog\xx.geojson c:\\Programmer\\R\\R-3.1.0\\library\\maptools\\shapes\\sids.shp
But this gave crs information:
ogr2ogr -f GeoJSON -s_srs EPSG:4008 -t_srs EPSG:4008 c:\Users\frtog\xx.geojson c:\\Programmer\\R\\R-3.1.0\\library\\maptools\\shapes\\sids.shp
Had to include both -s_srs and -t_srs to get it to work. EPSG:4008 specifies "+proj=longlat +ellps=clrk66" as the shape file was read with.
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling
Technology & Service Solutions
T +45 9730 5135
M +45 2547 6050
frtog at vestas.com
http://www.vestas.com
Company reg. name: Vestas Wind Systems A/S
This e-mail is subject to our e-mail disclaimer statement.
Please refer to www.vestas.com/legal/notice
If you have received this e-mail in error please contact the sender.
> -----Original Message-----
> From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-
> project.org] On Behalf Of Roger Bivand
> Sent: 12. august 2014 09:34
> To: Philip Greenwood
> Cc: r-sig-geo at r-project.org
> Subject: Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS
>
> On Tue, 12 Aug 2014, Philip Greenwood wrote:
>
> > Thanks very much for you help. My writeOGR belongs to rgdal version 0.8-
> 5.
>
> Which is from February 2013, current is 0.8-16. More importantly, we don't
> know which version of GDAL is loaded when rgdal loads - messages are
> printed on-screen to tell you - nor do we know how you installed rgdal.
>
> With GDAL 1.11.0, released 2014/04/16, I see:
>
> dsn <- system.file("vectors", package = "rgdal")[1]
> cities <- readOGR(dsn=dsn, layer="cities")
> summary(cities)
> #...
> #[+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0]
> #...
> td <- tempdir()
> writeOGR(cities, paste(td, ".geojson", sep=.Platform$file.sep),
> "OGRGeoJSON", driver="GeoJSON")
> summary(readOGR(paste(td, ".geojson", sep=.Platform$file.sep),
> "OGRGeoJSON"))
> #...
> #[+proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs]
> #...
>
> with the head of:
>
> file.show(paste(td, ".geojson", sep=.Platform$file.sep))
>
> "crs": { "type": "name", "properties": { "name":
> "urn:ogc:def:crs:OGC:1.3:CRS84" } },
>
> where urn:ogc:def:crs:OGC:1.3:CRS84 is WGS 84 longitude-latitude on p. 18
> of https://portal.opengeospatial.org/files/?artifact_id=24045.
>
> Not quite the same CRS, but a CRS is recorded. I think that your problem
> is related to the version of GDAL loaded into rgdal. The driver page on
> www.gdal.org relates to the current release and may include information
> for the development version. Your version seems to be old.
>
> >
> > gdalinfo --formats returns:
>
> ogrinfo --formats, or equivalently ogrDrivers() in R with rgdal loaded.
>
> Hope this clarifies,
>
> Roger
>
> >
> > Supported Formats:
> > VRT (rw+v): Virtual Raster
> > GTiff (rw+v): GeoTIFF
> > NITF (rw+v): National Imagery Transmission Format
> > RPFTOC (ro): Raster Product Format TOC format
> > HFA (rw+v): Erdas Imagine Images (.img)
> > SAR_CEOS (ro): CEOS SAR Image
> > CEOS (ro): CEOS Image
> > JAXAPALSAR (ro): JAXA PALSAR Product Reader (Level 1.1/1.5)
> > GFF (rov): Ground-based SAR Applications Testbed File Format (.gff)
> > ELAS (rw+): ELAS
> > AIG (ro): Arc/Info Binary Grid
> > AAIGrid (rwv): Arc/Info ASCII Grid
> > SDTS (ro): SDTS Raster
> > OGDI (ro): OGDI Bridge
> > DTED (rwv): DTED Elevation Raster
> > PNG (rwv): Portable Network Graphics
> > JPEG (rwv): JPEG JFIF
> > MEM (rw+): In Memory Raster
> > JDEM (ro): Japanese DEM (.mem)
> > GIF (rwv): Graphics Interchange Format (.gif)
> > BIGGIF (rov): Graphics Interchange Format (.gif)
> > ESAT (ro): Envisat Image Format
> > BSB (ro): Maptech BSB Nautical Charts
> > XPM (rw): X11 PixMap Format
> > BMP (rw+v): MS Windows Device Independent Bitmap
> > DIMAP (ro): SPOT DIMAP
> > AirSAR (ro): AirSAR Polarimetric Image
> > RS2 (ro): RadarSat 2 XML Product
> > PCIDSK (rw+v): PCIDSK Database File
> > PCRaster (rw): PCRaster Raster File
> > ILWIS (rw+v): ILWIS Raster Map
> > SGI (rw+): SGI Image File Format 1.0
> > SRTMHGT (rwv): SRTMHGT File Format
> > Leveller (rw+): Leveller heightfield
> > Terragen (rw+): Terragen heightfield
> > GMT (rw): GMT NetCDF Grid Format
> > netCDF (rw): Network Common Data Format
> > HDF4 (ro): Hierarchical Data Format Release 4
> > HDF4Image (rw+): HDF4 Dataset
> > ISIS3 (ro): USGS Astrogeology ISIS cube (Version 3)
> > ISIS2 (ro): USGS Astrogeology ISIS cube (Version 2)
> > PDS (ro): NASA Planetary Data System
> > TIL (ro): EarthWatch .TIL
> > ERS (rw+): ERMapper .ers Labelled
> > JPEG2000 (rwv): JPEG-2000 part 1 (ISO/IEC 15444-1)
> > L1B (ro): NOAA Polar Orbiter Level 1b Data Set
> > FIT (rw): FIT Image
> > GRIB (ro): GRIdded Binary (.grb)
> > RMF (rw+): Raster Matrix Format
> > WCS (ro): OGC Web Coverage Service
> > WMS (ro): OGC Web Map Service
> > MSGN (ro): EUMETSAT Archive native (.nat)
> > RST (rw+): Idrisi Raster A.1
> > INGR (rw+): Intergraph Raster
> > GSAG (rw): Golden Software ASCII Grid (.grd)
> > GSBG (rw+): Golden Software Binary Grid (.grd)
> > GS7BG (ro): Golden Software 7 Binary Grid (.grd)
> > COSAR (ro): COSAR Annotated Binary Matrix (TerraSAR-X)
> > TSX (ro): TerraSAR-X Product
> > COASP (ro): DRDC COASP SAR Processor Raster
> > R (rwv): R Object Data Store
> > PNM (rw+): Portable Pixmap Format (netpbm)
> > DOQ1 (ro): USGS DOQ (Old Style)
> > DOQ2 (ro): USGS DOQ (New Style)
> > ENVI (rw+v): ENVI .hdr Labelled
> > EHdr (rw+v): ESRI .hdr Labelled
> > GenBin (ro): Generic Binary (.hdr Labelled)
> > PAux (rw+): PCI .aux Labelled
> > MFF (rw+): Vexcel MFF Raster
> > MFF2 (rw+): Vexcel MFF2 (HKV) Raster
> > FujiBAS (ro): Fuji BAS Scanner Image
> > GSC (ro): GSC Geogrid
> > FAST (ro): EOSAT FAST Format
> > BT (rw+): VTP .bt (Binary Terrain) 1.3 Format
> > LAN (ro): Erdas .LAN/.GIS
> > CPG (ro): Convair PolGASP
> > IDA (rw+): Image Data and Analysis
> > NDF (ro): NLAPS Data Format
> > EIR (ro): Erdas Imagine Raw
> > DIPEx (ro): DIPEx
> > LCP (rov): FARSITE v.4 Landscape File (.lcp)
> > RIK (ro): Swedish Grid RIK (.rik)
> > USGSDEM (rw): USGS Optional ASCII DEM (and CDED)
> > GXF (ro): GeoSoft Grid Exchange Format
> > DODS (ro): DAP 3.x servers
> > HTTP (ro): HTTP Fetching Wrapper
> > BAG (ro): Bathymetry Attributed Grid
> > HDF5 (ro): Hierarchical Data Format Release 5
> > HDF5Image (ro): HDF5 Dataset
> > NWT_GRD (ro): Northwood Numeric Grid Format .grd/.tab
> > NWT_GRC (ro): Northwood Classified Grid Format .grc/.tab
> > ADRG (rw+v): ARC Digitized Raster Graphics
> > SRP (rov): Standard Raster Product (ASRP/USRP)
> > BLX (rw): Magellan topo (.blx)
> > Rasterlite (rw): Rasterlite
> > EPSILON (rwv): Epsilon wavelets
> > WKTRaster (ro): PostGIS WKT Raster driver
> > SAGA (rw+v): SAGA GIS Binary Grid (.sdat)
> >
> >
> >
> > On 12/08/2014, at 4:32 PM, Frede Aakmann T?gersen
> <frtog at vestas.com<mailto:frtog at vestas.com>>
> > wrote:
> >
> > Sorry. I was too hasty. I suppose that your problem is with the GDAL
> drivers. I see that http://www.gdal.org/drv_geojson.html points to
> http://www.gdal.org/drv_geojson.html.
> >
> > Which package belong writeOGR to. What is the version of the package?
> What is the version of your GDAL installation. What does e.g. 'gdalinfo --
> formats' show?
> >
> > Yours sincerely / Med venlig hilsen
> >
> >
> > Frede Aakmann T?gersen
> > Specialist, M.Sc., Ph.D.
> > Plant Performance & Modeling
> >
> > Technology & Service Solutions
> > T +45 9730 5135
> > M +45 2547 6050
> > frtog at vestas.com<mailto:frtog at vestas.com>
> > http://www.vestas.com
> >
> > Company reg. name: Vestas Wind Systems A/S
> > This e-mail is subject to our e-mail disclaimer statement.
> > Please refer to www.vestas.com/legal/notice
> > If you have received this e-mail in error please contact the sender.
> >
> >
> > -----Original Message-----
> > From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-
> > project.org] On Behalf Of Frede Aakmann T?gersen
> > Sent: 12. august 2014 08:19
> > To: Philip Greenwood; r-sig-geo at r-project.org
> > Subject: Re: [R-sig-Geo] WriteOGR to GeoJSON loses CRS
> >
> > Hi
> >
> > Google points me to http://geojson.org/geojson-spec.html. See Section 3
> > about CRSs.
> >
> > Yours sincerely / Med venlig hilsen
> >
> >
> > Frede Aakmann T?gersen
> > Specialist, M.Sc., Ph.D.
> > Plant Performance & Modeling
> >
> > Technology & Service Solutions
> > T +45 9730 5135
> > M +45 2547 6050
> > frtog at vestas.com
> > http://www.vestas.com
> >
> > Company reg. name: Vestas Wind Systems A/S
> > This e-mail is subject to our e-mail disclaimer statement.
> > Please refer to www.vestas.com/legal/notice
> > If you have received this e-mail in error please contact the sender.
> >
> >
> > -----Original Message-----
> > From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-
> > project.org] On Behalf Of Philip Greenwood
> > Sent: 12. august 2014 08:13
> > To: r-sig-geo at r-project.org
> > Subject: [R-sig-Geo] WriteOGR to GeoJSON loses CRS
> >
> > Hi I am using WriteOGR to write a SpatialPolygonsDataFrame to a
> > temporary
> > GeoJSON file. However the CRS seems to be lost in the resulting file.
> > Does the GeoJSON driver support CRSs?
> >
> > Thanks
> > Phil
> >
> > ---
> > Philip Greenwood
> > Product Manager
> > Australian Urban Research Infrastructure Network (AURIN)
> > Level 2 West, Alice Hoy Building, University of Melbourne
> > T: +61-(0)3-9035-8549
> > E:
> >
> >
> philip.greenwood at unimelb.edu.au<mailto:philip.greenwood at unimelb.edu
> > .au>
> >
> >
> >
> > [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo at r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo at r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
> > ---
> > Philip Greenwood
> > Product Manager
> > Australian Urban Research Infrastructure Network (AURIN)
> > Level 2 West, Alice Hoy Building, University of Melbourne
> > T: +61-(0)3-9035-8549
> > E:
> philip.greenwood at unimelb.edu.au<mailto:philip.greenwood at unimelb.edu
> .au>
> >
> >
> >
> > [[alternative HTML version deleted]]
> >
> >
>
> --
> Roger Bivand
> Department of Economics, Norwegian School of Economics,
> Helleveien 30, N-5045 Bergen, Norway.
> voice: +47 55 95 93 55; fax +47 55 95 91 00
> e-mail: Roger.Bivand at nhh.no
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
More information about the R-sig-Geo
mailing list