8. 坐标参考系¶
目标： 
理解坐标参考系。 

关键词： 
坐标参考系(CRS)，地图投影，动态投影，纬度，经度，北距、东距 
8.1. 概述¶
**地图投影**试图在一张纸或电脑屏幕上描绘地球的表面，或地球的一部分。用通俗的话来说，地图投影试图将地球从球形(3D)转化为平面(2D)。
然后，坐标参考系统 (CRS)定义地理信息系统中投影的二维地图如何与地球上的真实地点联系起来。使用哪种地图投影和CRS的决定取决于您想要工作的区域范围、您想要进行的分析，通常还取决于数据的可用性。
8.2. 详细地图投影¶
表示地球形状的传统方法是使用地球仪。然而，这种方法有一个问题。尽管地球仪保留了地球的大部分形状，并展示了大陆尺度特征的空间构型，但它们很难装在口袋里。它们也只方便在非常小的尺度下使用(例如1:100万)。
GIS应用中常用的专题地图数据大多比例尺较大。典型的GIS数据集具有1:250 000或更大的比例尺，取决于数据详细程度。这种大小的地球仪生产起来很困难，成本也很高，而且携带起来更困难。因此，制图者开发了一套被称为“地图投影”的技术，旨在以合理的精度显示二维球形地球。
从近处看，地球似乎是相对平坦的。然而，从太空看，我们可以看到地球相对来说是球形的。正如我们将在即将到来的地图制作主题中看到的，地图，是现实的表现。它们不仅被设计用来表现特征，而且也用来表现形状和空间布局。每个地图投影都有**优点**和**缺点**。地图的最佳投影取决于地图的**比例尺**以及使用目的。例如，一个投影用来绘制非洲大陆可能会存在不可接受的失真，但对你们国家的**大比例尺（详细）地图**来说，是一个很好的选择。地图投影的属性也可能会影响到地图的一些设计特征。有些投影在小区域制图很好用，有些适用于绘制大范围东西方向的区域，有些对大范围南北方向制图区域更合适。
8.3. 地图投影的三个系列¶
The process of creating map projections is best illustrated by positioning a light source inside a transparent globe on which opaque earth features are placed. Then project the feature outlines onto a twodimensional flat piece of paper. Different ways of projecting can be produced by surrounding the globe in a cylindrical fashion, as a cone, or even as a flat surface. Each of these methods produces what is called a map projection family. Therefore, there is a family of planar projections, a family of cylindrical projections, and another called conical projections (see 图 8.3)
Today, of course, the process of projecting the spherical earth onto a flat piece of paper is done using the mathematical principles of geometry and trigonometry. This recreates the physical projection of light through the globe.
8.4. Accuracy of map projections¶
Map projections are never absolutely accurate representations of the spherical earth. As a result of the map projection process, every map shows distortions of angular conformity, distance and area. A map projection may combine several of these characteristics, or may be a compromise that distorts all the properties of area, distance and angular conformity, within some acceptable limit. Examples of compromise projections are the Winkel Tripel projection and the Robinson projection (see 图 8.4), which are often used for producing and visualizing world maps.
It is usually impossible to preserve all characteristics at the same time in a map projection. This means that when you want to carry out accurate analytical operations, you need to use a map projection that provides the best characteristics for your analyses. For example, if you need to measure distances on your map, you should try to use a map projection for your data that provides high accuracy for distances.
8.4.1. Map projections with angular conformity¶
When working with a globe, the main directions of the compass rose (North, East, South and West) will always occur at 90 degrees to one another. In other words, East will always occur at a 90 degree angle to North. Maintaining correct angular properties can be preserved on a map projection as well. A map projection that retains this property of angular conformity is called a conformal or orthomorphic projection.
These projections are used when the preservation of angular relationships is important. They are commonly used for navigational or meteorological tasks. It is important to remember that maintaining true angles on a map is difficult for large areas and should be attempted only for small portions of the earth. The conformal type of projection results in distortions of areas, meaning that if area measurements are made on the map, they will be incorrect. The larger the area the less accurate the area measurements will be. Examples are the Mercator projection (as shown in 图 8.5) and the Lambert Conformal Conic projection. The U.S. Geological Survey uses a conformal projection for many of its topographic maps.
8.4.2. Map projections with equal distance¶
If your goal in projecting a map is to accurately measure distances, you should select a projection that is designed to preserve distances well. Such projections, called equidistant projections, require that the scale of the map is kept constant. A map is equidistant when it correctly represents distances from the centre of the projection to any other place on the map. Equidistant projections maintain accurate distances from the centre of the projection or along given lines. These projections are used for radio and seismic mapping, and for navigation. The Plate Carree Equidistant Cylindrical (see 图 8.6) and the Equirectangular projection are two good examples of equidistant projections. The Azimuthal Equidistant projection is the projection used for the emblem of the United Nations (see 图 8.7).
8.4.3. Projections with equal areas¶
When a map portrays areas over the entire map, so that all mapped areas have the same proportional relationship to the areas on the Earth that they represent, the map is an equal area map. In practice, general reference and educational maps most often require the use of equal area projections. As the name implies, these maps are best used when calculations of area are the dominant calculations you will perform. If, for example, you are trying to analyse a particular area in your town to find out whether it is large enough for a new shopping mall, equal area projections are the best choice. On the one hand, the larger the area you are analysing, the more precise your area measures will be, if you use an equal area projection rather than another type. On the other hand, an equal area projection results in distortions of angular conformity when dealing with large areas. Small areas will be far less prone to having their angles distorted when you use an equal area projection. Alber's equal area, Lambert's equal area and Mollweide Equal Area Cylindrical projections (shown in 图 8.8) are types of equal area projections that are often encountered in GIS work.
Keep in mind that map projection is a very complex topic. There are hundreds of different projections available world wide each trying to portray a certain portion of the earth's surface as faithfully as possible on a flat piece of paper. In reality, the choice of which projection to use, will often be made for you. Most countries have commonly used projections and when data is exchanged people will follow the national trend.
8.5. Coordinate Reference System (CRS) in detail¶
With the help of coordinate reference systems (CRS) every place on the earth can be specified by a set of three numbers, called coordinates. In general CRS can be divided into projected coordinate reference systems (also called Cartesian or rectangular coordinate reference systems) and geographic coordinate reference systems.
8.5.1. Geographic Coordinate Systems¶
The use of Geographic Coordinate Reference Systems is very common. They use degrees of latitude and longitude and sometimes also a height value to describe a location on the earth’s surface. The most popular is called WGS 84.
Lines of latitude run parallel to the equator and divide the earth into 180 equally spaced sections from North to South (or South to North). The reference line for latitude is the equator and each hemisphere is divided into ninety sections, each representing one degree of latitude. In the northern hemisphere, degrees of latitude are measured from zero at the equator to ninety at the north pole. In the southern hemisphere, degrees of latitude are measured from zero at the equator to ninety degrees at the south pole. To simplify the digitisation of maps, degrees of latitude in the southern hemisphere are often assigned negative values (0 to 90°). Wherever you are on the earth’s surface, the distance between the lines of latitude is the same (60 nautical miles). See 图 8.9 for a pictorial view.
Lines of longitude, on the other hand, do not stand up so well to the standard of uniformity. Lines of longitude run perpendicular to the equator and converge at the poles. The reference line for longitude (the prime meridian) runs from the North pole to the South pole through Greenwich, England. Subsequent lines of longitude are measured from zero to 180 degrees East or West of the prime meridian. Note that values West of the prime meridian are assigned negative values for use in digital mapping applications. See 图 8.9 for a pictorial view.
At the equator, and only at the equator, the distance represented by one line of longitude is equal to the distance represented by one degree of latitude. As you move towards the poles, the distance between lines of longitude becomes progressively less, until, at the exact location of the pole, all 360° of longitude are represented by a single point that you could put your finger on (you probably would want to wear gloves though). Using the geographic coordinate system, we have a grid of lines dividing the earth into squares that cover approximately 12363.365 square kilometres at the equator  a good start, but not very useful for determining the location of anything within that square.
To be truly useful, a map grid must be divided into small enough sections so that
they can be used to describe (with an acceptable level of accuracy) the location
of a point on the map. To accomplish this, degrees are divided into minutes
('
) and seconds ("
). There are sixty minutes in a degree, and sixty
seconds in a minute (3600 seconds in a degree). So, at the equator, one second
of latitude or longitude = 30.87624 meters.
8.5.2. Projected coordinate reference systems¶
A twodimensional coordinate reference system is commonly defined by two axes. At right angles to each other, they form a so called XYplane (see 图 8.10 on the left side). The horizontal axis is normally labelled X, and the vertical axis is normally labelled Y. In a threedimensional coordinate reference system, another axis, normally labelled Z, is added. It is also at right angles to the X and Y axes. The Z axis provides the third dimension of space (see 图 8.10 on the right side). Every point that is expressed in spherical coordinates can be expressed as an X Y Z coordinate.
A projected coordinate reference system in the southern hemisphere (south of the equator) normally has its origin on the equator at a specific Longitude. This means that the Yvalues increase southwards and the Xvalues increase to the West. In the northern hemisphere (north of the equator) the origin is also the equator at a specific Longitude. However, now the Yvalues increase northwards and the Xvalues increase to the East. In the following section, we describe a projected coordinate reference system, called Universal Transverse Mercator (UTM) often used for South Africa.
8.6. Universal Transverse Mercator (UTM) CRS in detail¶
The Universal Transverse Mercator (UTM) coordinate reference system has its origin on the equator at a specific Longitude. Now the Yvalues increase southwards and the Xvalues increase to the West. The UTM CRS is a global map projection. This means, it is generally used all over the world. But as already described in the section 'accuracy of map projections' above, the larger the area (for example South Africa) the more distortion of angular conformity, distance and area occur. To avoid too much distortion, the world is divided into 60 equal zones that are all 6 degrees wide in longitude from East to West. The UTM zones are numbered 1 to 60, starting at the antimeridian (zone 1 at 180 degrees West longitude) and progressing East back to the antemeridian (zone 60 at 180 degrees East longitude) as shown in 图 8.11.
As you can see in 图 8.11 and 图 8.12, South Africa is covered by four UTM zones to minimize distortion. The zones are called UTM 33S, UTM 34S, UTM 35S and UTM 36S. The S after the zone means that the UTM zones are located south of the equator.
Say, for example, that we want to define a twodimensional coordinate within the Area of Interest (AOI) marked with a red cross in 图 8.12. You can see, that the area is located within the UTM zone 35S. This means, to minimize distortion and to get accurate analysis results, we should use UTM zone 35S as the coordinate reference system.
The position of a coordinate in UTM south of the equator must be indicated with the zone number (35) and with its northing (Y) value and easting (X) value in meters. The northing value is the distance of the position from the equator in meters. The easting value is the distance from the central meridian (longitude) of the used UTM zone. For UTM zone 35S it is 27 degrees East as shown in 图 8.12. Furthermore, because we are south of the equator and negative values are not allowed in the UTM coordinate reference system, we have to add a so called false northing value of 10,000,000 m to the northing (Y) value and a false easting value of 500,000 m to the easting (X) value. This sounds difficult, so, we will do an example that shows you how to find the correct UTM 35S coordinate for the Area of Interest.
8.6.1. The northing (Y) value¶
The place we are looking for is 3,550,000 meters south of the equator, so the northing (Y) value gets a negative sign and is 3,550,000 m. According to the UTM definitions we have to add a false northing value of 10,000,000 m. This means the northing (Y) value of our coordinate is 6,450,000 m (3,550,000 m + 10,000,000 m).
8.6.2. The easting (X) value¶
First we have to find the central meridian (longitude) for the UTM zone 35S. As we can see in 图 8.12 it is 27 degrees East. The place we are looking for is 85,000 meters West from the central meridian. Just like the northing value, the easting (X) value gets a negative sign, giving a result of 85,000 m. According to the UTM definitions we have to add a false easting value of 500,000 m. This means the easting (X) value of our coordinate is 415,000 m (85,000 m + 500,000 m). Finally, we have to add the zone number to the easting value to get the correct value.
As a result, the coordinate for our Point of Interest, projected in UTM zone 35S would be written as: 35 415,000 m E / 6,450,000 m N. In some GIS, when the correct UTM zone 35S is defined and the units are set to meters within the system, the coordinate could also simply appear as 415,000 6,450,000.
8.7. OnTheFly Projection¶
As you can probably imagine, there might be a situation where the data you want to use in a GIS are projected in different coordinate reference systems. For example, you might get a vector layer showing the boundaries of South Africa projected in UTM 35S and another vector layer with point information about rainfall provided in the geographic coordinate system WGS 84. In GIS these two vector layers are placed in totally different areas of the map window, because they have different projections.
To solve this problem, many GIS include a functionality called onthefly projection. It means, that you can define a certain projection when you start the GIS and all layers that you then load, no matter what coordinate reference system they have, will be automatically displayed in the projection you defined. This functionality allows you to overlay layers within the map window of your GIS, even though they may be in different reference systems. In QGIS, this functionality is applied by default.
8.8. Common problems / things to be aware of¶
The topic map projection is very complex and even professionals who have studied geography, geodetics or any other GIS related science, often have problems with the correct definition of map projections and coordinate reference systems. Usually when you work with GIS, you already have projected data to start with. In most cases these data will be projected in a certain CRS, so you don't have to create a new CRS or even re project the data from one CRS to another. That said, it is always useful to have an idea about what map projection and CRS means.
8.9. What have we learned?¶
Let's wrap up what we covered in this worksheet:
Map projections portray the surface of the earth on a twodimensional, flat piece of paper or computer screen.
There are global map projections, but most map projections are created and optimized to project smaller areas of the earth's surface.
Map projections are never absolutely accurate representations of the spherical earth. They show distortions of angular conformity, distance and area. It is impossible to preserve all these characteristics at the same time in a map projection.
A Coordinate reference system (CRS) defines, with the help of coordinates, how the twodimensional, projected map is related to real locations on the earth.
There are two different types of coordinate reference systems: Geographic Coordinate Systems and Projected Coordinate Systems.
On the Fly projection is a functionality in GIS that allows us to overlay layers, even if they are projected in different coordinate reference systems.
8.10. Now you try!¶
Here are some ideas for you to try with your learners:
Start QGIS
In No projection (or unknown/nonEarth projection)
checkLoad two layers of the same area but with different projections
Let your pupils find the coordinates of several places on the two layers. You can show them that it is not possible to overlay the two layers.
Then define the coordinate reference system as Geographic/WGS 84 inside the Project Properties dialog
Load the two layers of the same area again and let your pupils see how setting a CRS for the project (hence, enabling "onthefly" projection) works.
You can open the Project Properties dialog in QGIS and show your pupils the many different Coordinate Reference Systems so they get an idea of the complexity of this topic. You can select different CRSs to display the same layer in different projections.
8.11. Something to think about¶
If you don't have a computer available, you can show your pupils the principles of the three map projection families. Get a globe and paper and demonstrate how cylindrical, conical and planar projections work in general. With the help of a transparency sheet you can draw a twodimensional coordinate reference system showing X axes and Y axes. Then, let your pupils define coordinates (X and Y values) for different places.
8.12. Further reading¶
Books:
Chang, KangTsung (2006). Introduction to Geographic Information Systems. 3rd Edition. McGraw Hill. ISBN: 0070658986
DeMers, Michael N. (2005). Fundamentals of Geographic Information Systems. 3rd Edition. Wiley. ISBN: 9814126195
Galati, Stephen R. (2006): Geographic Information Systems Demystified. Artech House Inc. ISBN: 158053533X
Websites:
https://foote.geography.uconn.edu/gcraft/notes/mapproj/mapproj_f.html
http://geology.isu.edu/wapi/geostac/Field_Exercise/topomaps/index.htm
The QGIS User Guide also has more detailed information on working with map projections in QGIS.
8.13. What's next?¶
In the section that follows we will take a closer look at Map Production.