Hot answers tagged coordinate-system
You should not have to do any particular conversion on your OSM data or Leaflet returned coordinates. There is a high confusion on this subject, but I think it comes down to a simple explanation: OSM data is in WGS84 datum (EPSG 4326). This is how you plot the Earth onto an ellipsoid. Web Mercator (EPSG 3857) is how you project this ellipsoid onto a plane....
sr.ImportFromEPSG() You have the S and P transposed. For the record, I don't see an entry for EPSG::29613 so that code may be invalid or deprecated. http://epsg.io/?q=29613
The integrized grid is sparse, in a sense it is a ragged array. Looks like the sinusoidal one really is a full array, albeit with empty tiles. The integrized sinusoidal grid is still used by the ocean colour group for L3BIN MODIS (and VIIRS, SeaWiFS, CZCS etc), recently converted from HDF4 to NetCDF4 (as compound types within groups) fwiw. The bins are ...
For proj.4 those two mean exactly the same and conversion from one system to another changes only the EPSG code. If you work with geodetic problems which require accurate ellipsoid model you will need some other tools. # ETRS89 <4258> +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs <> # WGS 84 <4326> +proj=longlat +datum=WGS84 ...
Do read the Leaflet tutorial on using WMS - it specifically addresses how to use WMS services in other projections. Also note that you can not have raster layers in different CRSs in Leaflet. It's simply not supported due to the complexity of on-the-fly raster reprojection.
It is a pure Python problem With shapely, your 3D point is represented by pt3D = Point(11.52951677300007,0.7729100360000416, -50000) list(pt3D.coords) [(11.52951677300007, 0.7729100360000416, -50000.0)] Using slicing for example (there are others solutions) pt2D = Point(list(pt3D.coords)[:2]) print pt2D POINT (11.52951677300007 0.7729100360000416) ...
Alternatively to your solution (which seems to work in ARCGIS), you can set up an oblique mercaor projection. This would work in QGIS and GDAL-based software as well: PROJ.4 : +proj=omerc +lat_0=55.102476 +lonc=1.42624 +alpha=-43.22 +k=0 +x_0=0 +y_0=0 +gamma=0 +datum=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs OGC WKT : PROJCS["unnamed", GEOGCS["WGS ...
Try making adding a new field with a long int as the data type. Right click on the new field and select Calculate Geometry Choose the option (length in meters) select a segment and use the measure tool to make sure it is the same length as it shows in the attribute table.
Only top voted, non community-wiki answers of a minimum length are eligible