On Tue, 12 Aug 2014, Frede Aakmann Tøgersen wrote:

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:

The definition is clearly not being recognised by the OGR driver, possibly needing a +towgs84= key or +datum= defintion:

proj4string(xx) <- CRS("+proj=longlat +ellps=clrk66 +datum=NAD27")

gives:

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

which is: +proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs

Without a +towgs84= or +datum=, it appears that the driver is choosing to omit an incomplete description.

Roger


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
fr...@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-boun...@r-project.org [mailto:r-sig-geo-bounces@r-
project.org] On Behalf Of Roger Bivand
Sent: 12. august 2014 09:34
To: Philip Greenwood
Cc: r-sig-geo@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
<fr...@vestas.com<mailto:fr...@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
fr...@vestas.com<mailto:fr...@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-boun...@r-project.org [mailto:r-sig-geo-bounces@r-
project.org] On Behalf Of Frede Aakmann T?gersen
Sent: 12. august 2014 08:19
To: Philip Greenwood; r-sig-geo@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
fr...@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-boun...@r-project.org [mailto:r-sig-geo-bounces@r-
project.org] On Behalf Of Philip Greenwood
Sent: 12. august 2014 08:13
To: r-sig-geo@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.greenw...@unimelb.edu.au<mailto:philip.greenw...@unimelb.edu
.au>



[[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@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.greenw...@unimelb.edu.au<mailto:philip.greenw...@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.biv...@nhh.no

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


--
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.biv...@nhh.no
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to