1suppressPackageStartupMessages(library(rgdal)) 2getPROJ4VersionInfo() 3getGDALVersionInfo() 4d <- system.file("vectors", package="rgdal") 5#shps <- ogrListLayers(d) 6(o <- OGRSpatialRef(d, "trin_inca_pl03")) 7if (new_proj_and_gdal()) comment(o) 8(o <- OGRSpatialRef(d, "kiritimati_primary_roads")) 9if (new_proj_and_gdal()) comment(o) 10(o <- OGRSpatialRef(d, "cities")) 11if (new_proj_and_gdal()) comment(o) 12(o <- OGRSpatialRef(d, "scot_BNG")) 13if (new_proj_and_gdal()) comment(o) 14(o <- OGRSpatialRef(file.path(d, "airports.gml"), "airports")) 15if (new_proj_and_gdal()) comment(o) 16(o <- OGRSpatialRef(file.path(d, "ps_cant_31.MIF"), "ps_cant_31")) 17if (new_proj_and_gdal()) comment(o) 18(o <- OGRSpatialRef(file.path(d, "Up.tab"), "Up")) 19if (new_proj_and_gdal()) comment(o) 20(o <- OGRSpatialRef(file.path(d, "test_trk2.gpx"), "tracks")) 21if (new_proj_and_gdal()) comment(o) 22(o <- OGRSpatialRef(file.path(d, "PacoursIKA2.TAB"), "PacoursIKA2")) 23if (new_proj_and_gdal()) comment(o) 24d <- system.file("pictures", package="rgdal") 25#f <- c("SP27GTIF.TIF", "cea.tif", "erdas_spnad83.tif", "scaleoffset.vrt", 26# "test_envi_class.envi") 27ds <- GDAL.open(file.path(d, "SP27GTIF.TIF")) 28(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=TRUE)) 29if (new_proj_and_gdal()) comment(o) 30(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=FALSE)) 31if (new_proj_and_gdal()) comment(o) 32GDAL.close(ds) 33ds <- GDAL.open(file.path(d, "cea.tif")) 34(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=TRUE)) 35if (new_proj_and_gdal()) comment(o) 36(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=FALSE)) 37if (new_proj_and_gdal()) comment(o) 38GDAL.close(ds) 39ds <- GDAL.open(file.path(d, "erdas_spnad83.tif")) 40(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=TRUE)) 41if (new_proj_and_gdal()) comment(o) 42(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=FALSE)) 43if (new_proj_and_gdal()) comment(o) 44GDAL.close(ds) 45ds <- GDAL.open(file.path(d, "scaleoffset.vrt")) 46(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=TRUE)) 47if (new_proj_and_gdal()) comment(o) 48(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=FALSE)) 49if (new_proj_and_gdal()) comment(o) 50GDAL.close(ds) 51ds <- GDAL.open(file.path(d, "test_envi_class.envi")) 52(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=TRUE)) 53if (new_proj_and_gdal()) comment(o) 54(o <- getProjectionRef(ds, OVERRIDE_PROJ_DATUM_WITH_TOWGS84=FALSE)) 55if (new_proj_and_gdal()) comment(o) 56GDAL.close(ds) 57 58