It is not possible to develop a map projection that does not lead to some disortion of the patterns that appear on the surface of the earth. The extent to which different projections distort the original images/distortion can be studied with Maple. The fundamental idea behind an invesigation of map distortion is to plot a simple shape (X's and O's, for example) and see how that shape changes when plotted in a different coordinate system. The simplest way to look at map distortion is plot simple shapes in rectangular coordinates. Here we plot a number of circles (and semi-circles) on the centers of a 30 degree equirectangular graticule. The above method of looking at map distortion is not the best approach because the original image (the equirectangular projection) already distorts the original outlines on the Earth's surface. A better approach, developed by Tissot, is to project the image of circles drawn on the surface of a sphere. In what follows we show (admittedly without much explanation) one way to draw circles on the surface of a sphere. T[1] := matrix(3,3,[[sin(theta)*cos(phi), sin(theta)*sin(phi), cos(theta)], [cos(theta)*cos(phi), cos(theta)*sin(phi), -sin(theta)], T_1[1] := matrix(3,3,[[sin(theta)*cos(phi), cos(theta)*cos(phi), -sin(phi)], [sin(theta)*sin(phi), cos(theta)*sin(phi), cos(phi)], pp := plots[spacecurve]( subs(temp,r = 0.1, convert(tt1, list)) , lambda = -Pi..Pi, color = black, thickness = 3, axes = boxed): pp := spacecurve( subs(temp,r = 0.1, theta = Pi/4, phi = Pi/4, {convert(tt1, list), convert(tt2, list), convert(tt3, list)}) , lambda = -Pi..Pi, color = black, thickness = 3, axes = boxed,orientation=[-140,60]): dd := evalf(seq(op({seq(subs(temp,r = 0.2, theta = 2*Pi/12*ii, phi = Pi/6*jj, convert(tt3, list)), ii = [2,3,4])}), jj = 1..12)): dd: ee := evalf(seq(op({seq(subs(temp,r = 0.2, theta = 2*Pi/12*ii, phi = Pi/6*jj*2, convert(tt3, list)), ii = [1,5])}), jj = 1..6)): ee: gg := evalf(seq(op({seq(subs(temp,r = 0.2, theta = 2*Pi/12*ii, phi = Pi/6*jj*2, convert(tt3, list)), ii = [0,6])}), jj = 1)): gg: We now wish to project the surface of the above sphere using one of the map coordinate systems discussed elsewhere. We have automated the creation of distortion images/distortion in a procedure included in the Maps package. The procedure must be accessed using the readlib command. It is now a trivial exercise to create distortion images/distortion of any other projection. In what follows we show some representative examples. Note that in an equal area projection the area of all complete circles is the same, although the shape suffers somewhat. The degree of distortion introduced by the Mercator projection is evident from the illustrations below. mp:=[`Cylindrical Equal Area`, `Eckert IV`, `Lambert Conical Equal Area`, Mercator, Mollweide, Sinusoidal, `Transverse Mercator`, `interrupted Boggs`], Please note that much of the Application Center contains content submitted directly from members of our user community. Although we do our best to monitor for objectionable content, it is possible that we occasionally miss something. If there is something objectionable on this page, please Source.


Яндекс.Метрика Рейтинг@Mail.ru Free Web Counter
page counter
Last Modified: April 22, 2016 @ 12:05 pm