They are not the same. EPSG:4326 refers to WGS 84 whereas EPSG:900913 refers to WGS84 Web Mercator. EPSG:4326 treats the earth as an ellipsoid while EPSG:900913 treats it as a sphere. This affects calculations done based on treating the map as a flat plane which is why your features got plotted on the wrong places.
900913 = GOOGLE (spelled with numbers). There is no official EPSG code 900913, because for some time this reference system wasn't added into the EPSG list of spatial reference systems. After some time, it eventually was under the code 3857 (and there was even some other code for some time). Looking at their properties, they are the same, so you could use one ...
If an Esri well-known ID is below 32767, it corresponds to the EPSG ID. WKIDs that are 32767 or above are Esri-defined. Either the object isn't in the EPSG Geodetic Parameter Dataset yet, or it probably won't be added. If an object is later added to the EPSG Dataset, Esri will update the WKID to match the EPSG one, but the previous value will still work.
There are three different ways to change/manipulate projections in QGIS:
1. Project Properties > CRS (or you can click on the icon where EPSG is written in the bottom right corner)
This defines the proprieties of the canvas (the background layer). For exemple, if I want to display a background map from openlayers plugin (google street etc.) I will need to ...
In ol.proj.transform you need to specify the fromProjection before the toProjection, then it should work.
As Michael Gentry explains in his answer, another Problem is that you have to specify the longitude (west-east thus x) first and then the latitude (south-north thus y).
And a better way to set the center is to get the current view and set the center ...
Assuming your areas of interest are rather small compared to the globe, you could set up a custom transverse mercator projection.
You have to know the geographic coordinates lat_0 and lon_0 of the origin of your CRS, and the direction of x and y axis:
+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Your solution is a little out of date (look at Fiona - Preffered method for defining a schema). For a better solution look at that proposed by Sean Gillies in gistfile1.py to parse a delimited text file data and create a new shapefile
1) How do I define the spatial reference: Fiona crs Module
from fiona.crs import from_epsg
It is a projection of a spheroid on a flat surface. Every projection has strengths and weaknesses and will preserve some elements of direction, distance or area better or worse than others (which is why careful selection of a suitable, that is local where possible, projection is so important).
So, while the unit of measurement in EPSG:3857 is indeed ...
There is no EPSG SRID 102743. Note that EPSG is the authority, and 102743 is the SRID. If you look up SRID 102743 on spatialreference.org, the listing is for ESRI:102743, meaning that ESRI (the publishers of ArcGIS) is the authority, not EPSG (European Petroleum Survey Group, now absorbed by the International Association of Oil & Gas Producers). The ...
QGIS handles SRS much the same as ESRI. Dive into QGIS and you'll find all your favorite SRS waiting for you.
EPSG is the European Petroleum Survey Group. It provides 'official' spatial reference systems and spheroids and its codes are ubiquitous. ESRI
takes most of its SRS from the EPSG dataset. So you are already using the ...
The EPSG:32719 is a UTM Zone 19 south setup with:
$ grep -B1 '<32719>' /usr/share/proj/epsg
# WGS 84 / UTM zone 19S
<32719> +proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs <>
which is a shortcut defined in /usr/share/proj/epsg. For the background you can refer to the PROJ.4 description. If you use the full parameter set for ...
9820 is the EPSG code for the laea projection definition. You can find it at https://www.epsg-registry.org/ under retrieve by code. This code numer has nothing to do with the commonly used EPSG codes like 4326 and 3857. These include projection parameters specified for a certain country or region.
For laea centered on Europe, the code is EPSG:3035, with ...
I found the following workaround. I am unsure if it is the most efficient, but it does work for me.
path = r"C:\temp\test2.tif"
d = gdal.Open(path)
proj = osr.SpatialReference(wkt=d.GetProjection())
EDIT: Slightly more condensed
This appears to be documented in the SDK help, e.g. IGeometryServer.FindSRByWKID:
AuthorityName is usually "EPSG" or "ESRI", but can also be an arbitrary string. It can also be the empty string if you want the default authority name associated with the new spatial reference. Clients can associate their own authority names with factory codes that are ...
I have a new stack exchange account and don't have a reputation high enough to comment on the "Uncaught TypeError: Failed to execute 'putImageData' on 'CanvasRenderingContext2D': float parameter 3 is non-finite." error. This occurs because you have the lat and the long inputs backwards.
map.getView().setCenter(ol.proj.transform([lat, long], 'EPSG:4326', '...
You did nothing wrong. It's just the fact that a Mercator projection can not display the south pole.
The dataset you have contains a whole line at the bottom on 90°S, which can mathematically not be reprojected to Mercator. The points would be in inifinity.
Easy solution is to clip the Antarctic at 89°S, i.e. remove the bottom line vertices before ...
Simply use (Converting elevations into correct units with pyproj?, Proj4 String for NAD83(2011) / Louisiana South (ftUS), ...)
preserve_units=True (as you say, pyproj assumes that your coordinates are in meters, therefore)
for northing, easting, up in targets:
print "Converting", easting, northing
# Transform using pyproj (gives wrong answer)
KML files are always WGS84/lat lon (EPSG:4326). Save the KML file as an EPSG:27700 shapefile instead.
Alternatively, you could download the Ordnance Survey OpenData Boundary-Line dataset, which is already in EPSG:27700. See:
The official OGC “Well-known Text Representation of Spatial Reference Systems” for EPSG 4326 (http://spatialreference.org/ref/epsg/4326/ogcwkt/) is (your second projection):
You can also do this in one line using gdal.Info as so:
epsg = int(gdal.Info(input, format='json')['coordinateSystem']['wkt'].rsplit('"EPSG","', 1)[-1].split('"'))
This does essentially the same thing as the answer provided by Lennert. It reads the WKT representation of the file's spatial reference, then parses the string to extract the EPSG.
EPSG:31463 seems to be deprecated and EPSG:31467 is the correct one.
http://spatialreference.org/ref/epsg/dhdn-3-degree-gauss-zone-3/ is stating:
ProjCRS related to inappropriate coordinate system resulting in incorrect axes and name. Use EPSG:31467 instead.
you can see see this, if you compare the AXIS definitions:
PROJCS["DHDN / 3-degree Gauss zone ...
As you found, EUREF-FIN is a realization of ETRS89 in Finland.
Some countries, when displaying coordinates or transmitting data, use the convention northing, easting order while other countries use easting, northing order. Most GIS software wants easting, northing so the EPSG registry (http://www.epsg-registry.org NOT epsg.io) often has both versions ...
I am not sure anyy open source supports this geomagnetic to geographic.
But if you are having few coordinates, pls. try this..
NASA has published the algorithm, you can try
The first archive ("Roads Data - tables and layers") doesn't contain the projection information, which is pretty bad form. However, the second archive ("Roads ArcReader Projects") has the projection information in its .prj files:
This is a feature, not a bug. Most Leaflet users are not aware that the map's display CRS is different from the map's data CRS, and that's fine; it makes life simpler for most people.
If your question would have been...
I have data in EPSG:3857 but Leaflet will only accept EPSG:4326 LatLngs, how can I add my data?
Then the answer would have been: