GEO2311
Web Mercator
The International Association of Oil & Gas Producers (IOGP)
Geomatics Committee
"If you’re drilling a well in a potentially sensitive area and you’re spending millions of dollars to do it, you want to be certain that well is in the right place.
The IOGP Geomatics Committee helps to provide some of that certainty.
Committee members include specialists in surveying, geodesy, cartography and geospatial data management."
EPSG Dataset
"The EPSG Geodetic Parameter Dataset is a collection of definitions of coordinate reference systems and coordinate transformations. These may be global, regional, national or local in application. To learn more visit www.epsg.org.
The Dataset is maintained by IOGP members through the Geodesy Subcommittee of the Geomatics Committee. It is distributed publicly at no cost as both an online (www.epsg-registry.org) and offline (www.epsg.org) resource, both requiring acceptance of the Terms of Use."
Web mercator
- datum: sfærisk - ikke ellipsoide-basert
- brukes av
Bing Maps
Open Street Map
Google (?)
- EPSG:3857
- se http://epsg.io/3857
- før Web mercator ble lagt til i EPSG-datasette, ble koden 900913 brukt (tallkoden ligner på teksten "google")
PROJCS["WGS 84 / Pseudo-Mercator",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["X",EAST],
AXIS["Y",NORTH],
EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0
+lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"],
AUTHORITY["EPSG","3857"]]EPSG 3857, med ellipsoidisk datum (!)
Opprinnelig EPSG-definisjon for Web Mercator EPSG:3785
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,
AUTHORITY["EPSG","7059"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6055"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4055"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["X",EAST],
AXIS["Y",NORTH],
EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
+k=1.0 +units=m +nadgrids=@null +wktext +no_defs"],
AUTHORITY["EPSG","3785"]]EPSG:4326
- geografiske koordinater
- ikke projisert
- Brukes i
GPS
GPX-filer
KML-filer
Hvordan ser Web mercator ut?
demo i ArcGIS
UTM sone 32, areal i kvadratmeter
Web Mercator, areal i kvadratmeter
Equirektangulær projeksjon (Plate Carrée), areal i kvadratmeter (denne projeksjonen velges automatisk av ArcGIS hvis datasettet har geografiske koordinater, EPSG:4326)
Arealberegning med Google Maps API
<!DOCTYPE html>
<html>
<head>
<title>Arealberegning</title>
<meta charset="utf-8">
<script src="http://ajax.googleapis.com/ajax/libs/jquery/1.10.2/jquery.min.js"></script>
<script src="https://maps.googleapis.com/maps/api/js?v=3.exp&libraries=geometry&sensor=false"></script>
<script>
$(document).ready(function(){
var p1 = new google.maps.LatLng(10.68012091, 60.78962180);
var p2 = new google.maps.LatLng(10.68108467, 60.78905971);
var p3 = new google.maps.LatLng(10.67966384, 60.78851428);
var p4 = new google.maps.LatLng(10.67890237, 60.78898582);
var p5 = new google.maps.LatLng(10.68004903, 60.78956696);
var p6 = new google.maps.LatLng(10.68012091, 60.78962180);
var area = google.maps.geometry.spherical.computeArea([p1,p2,p3,p4,p5]);
$('div').append("Areal: " + area);
});
</script>
</head>
<body>
<div></div>
</body>
</html>Arealberegning med Google Maps API
Påstand fra undertegnede:
Hvis Google opererer med samme lengde på en lengdegrad overalt på kloden,
og samme lengde på breddegrader overalt på kloden,
.. så er arealberegningene i Google Maps Spherical geometry-biblioteket utført i en equirektangulær projeksjon (Plate Carrée) - det samme som standard-projeksjonen i ArcMap. Dette kommer i så fall ikke klart fram fra dokumentasjonen.
https://developers.google.com/maps/documentation/javascript/geometry
Arealberegninger med
geografisk datatype i PostGIS
(koordinater oppgitt som lengde/breddegrad)
select ST_Area(ST_GeogFromText('SRID=4326;
POLYGON((1 1, 2 1, 2 2, 1 2, 1 1))'))
Resultat ved ekvator (breddegrad 1-2): 12304814950.0731
(Dette resultatet er tilnærmet likt Googles resultat)
select ST_Area(ST_GeogFromText('SRID=4326;
POLYGON((10 60, 11 60, 11 61, 10 61, 10 60))'))
Resultat ved 60 grader nord: 3108618015.39229
Arealet på 60 grader nord er ca. 50 % av arealet ved ekvator.
HIG
Verktøy Projeksjon Areal % avvik fra seeiendom.no
ArcMap Web Mercator 29989,0 318
Equirectangular 14635,4 104
UTM Sone 32 7163,3 0
PostGIS Geografisk 7167,5 0
Google Maps Equirectangular (?) 14381,8 101
seeiendom.no UTM Sone 32 7167,0 0
Boligtomt i Østre Toten
Verktøy Projeksjon Areal % avvik fra seeiendom.no
ArcMap Web Mercator 4774,9 316
Equirectangular 2340,8 104
UTM Sone 32 1151,0 0
PostGIS Geografisk 1151,2 0
Google Maps Equirectangular (?) 2297,5 100
seeiendom.no UTM Sone 32 1147,4 0
Oppsummert - resultater fra arealberegninger