Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- S-JTSK_Krovak_East_North
- WKID: 5514 Authority: EPSG
- Projection: Krovak
- False_Easting: 0.0
- False_Northing: 0.0
- Pseudo_Standard_Parallel_1: 78.5
- Scale_Factor: 0.9999
- Azimuth: 30.28813975277778
- Longitude_Of_Center: 24.83333333333333
- Latitude_Of_Center: 49.5
- X_Scale: -1.0
- Y_Scale: 1.0
- XY_Plane_Rotation: 90.0
- Linear Unit: Meter (1.0)
- Geographic Coordinate System: GCS_S_JTSK
- Angular Unit: Degree (0.0174532925199433)
- Prime Meridian: Greenwich (0.0)
- Datum: D_S_JTSK
- Spheroid: Bessel_1841
- Semimajor Axis: 6377397.155
- Semiminor Axis: 6356078.962818189
- Inverse Flattening: 299.1528128
- S-JTSK_Krovak_East_North
- WKID: 5514 Authority: EPSG
- Projection: Krovak
- False_Easting: 0.0
- False_Northing: 0.0
- Pseudo_Standard_Parallel_1: 78.5
- Scale_Factor: 0.9999
- Azimuth: 30.28813975277778
- Longitude_Of_Center: 24.83333333333333
- Latitude_Of_Center: 49.5
- X_Scale: -1.0
- Y_Scale: 1.0
- XY_Plane_Rotation: 90.0
- Linear Unit: Meter (1.0)
- Geographic Coordinate System: GCS_S_JTSK
- Angular Unit: Degree (0.0174532925199433)
- Prime Meridian: Greenwich (0.0)
- Datum: D_S_JTSK
- Spheroid: Bessel_1841
- Semimajor Axis: 6377397.155
- Semiminor Axis: 6356078.962818189
- Inverse Flattening: 299.1528128
- library(sp)
- library(rgdal)
- library(rgeos)
- library(raster)
- # import raster
- r1<-raster("F:/.../r1_2013.img")
- #import shp
- d.shp <- readOGR(dsn = "F:/.../dummy_data_arcgis_r",
- layer = "np_by_manag")
- > proj4string(r1) # RASTER
- [1] "+proj=krovak +lat_0=24.83333333333333 +lon_0=49.5 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
- > proj4string(d.shp) # SHAPEFILE
- [1] "+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs"
- >
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement