maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git


  
1

	
2

	
3

	
4

	
5

	
6

	
7

	
8

	
9

	
10

	
11

	
12

	
13

	
14

	
15

	
16

	
17

	
18

	
19

	
20

	
21

	
22

	
23

	
24

	
25

	
26

	
27

	
28

	
29

	
30

	
31

	
32

	
33

	
34

	
35

	
36

	
37

	
38

	
39

	
40

	
41

	
42

	
43

	
44

	
45

	
46

	
47

	
48

	
49

	
50

	
51

	
52

	
53

	
54

	
55

	
56

	
57

	
58

	
59

	
60

	
61

	
62

	
63

	
64

	
65

	
66

	
67

	
68

	
69

	
70

	
71

	
72

	
73

	
74

	
75

	
76

	
77

	
78

	
79

	
80

	
81

	
82

	
83

	
84

	
85

	
86

	
87

	
88

	
89

	
90

	
91

	
92

	
93

	
94

	
95

	
96

	
97

	
98

	
99

	
100

	
101

	
102

	
103

	
104

	
105

	
106

	
107

	
108

	
109

	
110

	
111

	
112

	
113

	
114

	
115

	
116

	
117

	
118

	
119

	
120

	
121

	
122

	
123

	
124

	
125

	
126

	
127

	
128

	
129

	
130

	
131

	
132

	
133

	
134

	
135

	
136

	
137

	
138

	
139

	
140

	
141

	
142

	
143

	
144

	
145

	
146

	
147

	
148

	
149

	
150

	
151

	
152

	
153

	
154

	
155

	
156

	
157

	
158

	
159

	
160

	
161

	
162

	
163

	
164

	
165

	
166

	
167

	
168

	
169

	
170

	
171

	
172

	
173

	
174

	
175

	
176

	
177

	
178

	
179

	
180

	
181

	
182

	
183

	
184

	
185

	
186

	
187

	
188

	
189

	
190

	
191

	
192

	
193

	
194

	
195

	
196

	
197

	
198

	
199

	
200

	
201

	
202

	
203

	
204

	
205

	
206

	
207

	
208

	
209

	
210

	
211

	
212

	
213

	
214

	
215

	
216

	
217

	
218

	
219

	
220

	
221

	
222

	
223

	
224

	
225

	
226

	
227

	
228

	
229

	
230

	
231

	
232

	
233

	
234

	
235

	
236

	
237

	
238

	
239

	
240

	
241

	
242

	
243

	
244

	
245

	
246

	
247

	
248

	
249

	
250

	
251

	
252

	
253

	
254

	
255

	
256

	
257

	
258

	
259

	
260

	
261

	
262

	
263

	
264

	
265

	
266

	
267

	
268

	
269

	
270

	
271

	
272

	
273

	
274

	
275

	
276

	
277

	
278

	
279

	
280

	
281

	
282

	
283

	
284

	
285

	
286

	
287

	
288

	
289

	
290

	
291

	
292

	
293

	
294

	
295

	
296

	
297

	
298

	
299

	
300

	
301

	
302

	
303

	
304

	
305

	
306

	
307

	
308

	
309

	
310

	
311

	
312

	
313

	
314

	
315

	
316

	
317

	
318

	
319

	
320

	
321

	
322

	
323

	
324

	
325

	
326

	
327

	
328

	
329

	
330

	
331

	
332

	
333

	
334

	
335

	
336

	
            #!/usr/bin/env bash
#
#	Testing GMT gmt mapproject on examples in Snyder.
#
#
#	usage:  snyder.sh
#
#	script will report trouble if we do not match Snyder (+- 360 degrees for longitudes)
#	Slop is 0.11 m for gmt projected values

blabber () {
	gf=($(echo "${sf[*]}" | gmt mapproject $* -F -C))
	gi=($(echo "${gf[*]}" | gmt mapproject $* -F -C -I --FORMAT_FLOAT_OUT=%lg))
	echo "${sf[*]} ${gi[*]} ${si[*]} ${gf[*]} $p" | $AWK -f test.awk
}

cat << EOF > test.awk
#!/bin/awk -f
{
	true_lon = \$1
	true_lat = \$2
	GMT_lon = \$3
	GMT_lat = \$4
	true_x = \$5
	true_y = \$6
	GMT_x = \$7
	GMT_y = \$8

#	Check longitudes & latitudes

	dx = GMT_lon - true_lon
	dy = GMT_lat - true_lat
	if (dy < 0.0) dy = -dy;
	if (dx < 0.0) dx = -dx;
	if (dx > 350) dx -= 360.0;
	if (! (dx == 0.0 || dx == -360.0 || dx = 360.0)) {
		printf "%s: Bad longitude conversion, d = %lg\n", \$9, dx
	}
	if (! (dy == 0.0)) {
		printf "%s: Bad latitude conversion, d = %lg\n", \$9, dy
	}

#	Check x & y

	dx = GMT_x - true_x
	dy = GMT_y - true_y
	if (dy < 0.0) dy = -dy;
	if (dx < 0.0) dx = -dx;
	if (! (dx <= 0.11)) {
		printf "%s: Bad x conversion, d = %lg\n", \$9, dx
	}
	if (! (dy <= 0.11)) {
		printf "%s: Bad y conversion, d = %lg\n", \$9, dy
	}
}
EOF

gmt set PROJ_SCALE_FACTOR 1

# Make ellipsoids (spheres) of radius 1 and 3:
unit=1.0
three=3.0
# Make special Clarke-1866 ellipsoid using Snyders value for e^2=0.00676866
# which gives 1/f = 294.97861076:
snyder=6378206.4,294.97861076
# International 1924 ellipsoid
int=International-1924

#----------------------------------------------------------------------------
# Mercator, Spherical

sf=(-75 35)
si=(1.8325957 0.6528366)
p="Mercator(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-50/50 -Jm180/0/1:1 > fail
#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Mercator, Ellipsoidal

sf=(-75 35)
si=(11688673.7 4139145.6)
p="Mercator(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --FORMAT_FLOAT_OUT=%9.1lf -R0/360/-50/50 -Jm180/0/1:1 >> fail
#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Transverse Mercator, Spherical

sf=(-73.5 40.5)
si=(0.0199077 0.7070276)
p="TM(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-77/-72/35/45 -Jt-75/0/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Oblique Mercator, Spherical
# Commented out since GMT and Snyder use different origins...
#
#sf=(120 -30)
#si=(-2.4201335 -0.0474026)
#p="Oblique-Mercator(Spherical)"
#blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-80/80 -Job0/45/-90/0/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Cylindrical Equal-Area, Spherical

sf=(80 35)
si=(2.3428242 0.6623090)
p="Cylindrical-Equal-Area(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-255/105/-80/80 -Jy-75/30/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Cylindrical Equal-Area, Ellipsoidal

sf=(-78 5)
si=(-332699.8 554248.5)
p="Cylindrical-Equal-Area(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --PROJ_SCALE_FACTOR=0.9996 --FORMAT_FLOAT_OUT=%9.1lf -R-255/105/-80/80 -Jy-75/5/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Miller, Spherical

sf=(-75 50)
si=(-1.3089969 0.9536371)
p="Miller(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-180/180/-80/80 -Jj0/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Cassini, Spherical

sf=(-90 25)
si=(-0.2367759 0.7988243)
p="Cassini(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-77/-73/-23/-17 -Jc-75/-20/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Cassini, Ellipsoidal

sf=(-73 43)
si=(163071.1 335127.6)
p="Cassini(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --FORMAT_FLOAT_OUT=%9.1lf -R-76/-72/38/42 -Jc-75/40/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Albers, Spherical

sf=(-75 35)
si=(0.2952720 0.2416774)
p="Albers(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-100/-90/30/40 -Jb-96/23/29.5/45.5/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Albers, Ellipsoidal

sf=(-75 35)
si=(1885472.7 1535925.0)
p="Albers(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --FORMAT_FLOAT_OUT=%9.1lf -R-100/-90/30/40 -Jb-96/23/29.5/45.5/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Lambert Conic, Spherical

sf=(-75 35)
si=(0.2966785 0.2462112)
p="Lambert-Conic(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-100/-90/30/40 -Jl-96/23/33/45/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Lambert Conic, Ellipsoidal

sf=(-75 35)
si=(1894410.9 1564649.5)
p="Lambert-Conic(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --FORMAT_FLOAT_OUT=%9.1lf -R-100/-90/30/40 -Jl-96/23/33/45/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Equidistant Conic, Spherical

sf=(-75 35)
si=(0.2952057 0.2424021)
p="Equidistant-Conic(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R-100/-90/30/40 -Jd-96/23/29.5/45.5/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Stereographic, Spherical

sf=(-75 30)
si=(0.3807224 -0.1263802)
p="Stereographic(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-90/90 -Js-100/40/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Stereographic, Ellipsoidal

sf=(-90 30)
si=(971630.8 -1063049.3)
p="Stereographic-General(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --PROJ_SCALE_FACTOR=0.9999 --FORMAT_FLOAT_OUT=%9.1lf -R0/360/-90/90 -Js-100/40/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Stereographic, Ellipsoidal

sf=(150 -75)
si=(-1573645.4 -572760.1)
p="Stereographic-Polar(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$int --PROJ_SCALE_FACTOR=0.994 --FORMAT_FLOAT_OUT=%9.1lf -R0/360/-90/0 -Js-100/-90/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Gnomonic, Spherical

sf=(-110 30)
si=(-0.1542826 -0.1694739)
p="Gnomonic(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-90/90 -Jf-100/40/60/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Lambert Azim, Spherical

sf=(100 -20)
si=(-4.2339303 4.0257775)
p="Lambert-Azimuthal(Spherical)"
blabber --PROJ_ELLIPSOID=$three --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-90/90 -Ja-100/40/180/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Lambert Azim, Ellipsoidal Polar

sf=(5 80)
si=(1077459.7 288704.5)
p="Lambert-Azimuthal-Polar(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$int --FORMAT_FLOAT_OUT=%9.1lf -R0/360/0/90 -Ja-100/90/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Lambert Azim, Ellipsoidal

sf=(-110 30)
si=(-965932.1 -1056814.9)
p="Lambert-Azimuthal-General(Ellipsoidal)"
blabber --PROJ_ELLIPSOID=$snyder --FORMAT_FLOAT_OUT=%9.1lf -R0/360/-90/90 -Ja-100/40/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Azim Equidistant, Spherical

sf=(100 -20)
si=(-5.8311398 5.5444634)
p="Azimuthal-Equidistant(Spherical)"
blabber --PROJ_ELLIPSOID=$three --FORMAT_FLOAT_OUT=%9.7lf -R0/360/-90/90 -Je-100/40/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Van der Grinten, Spherical

sf=(-160 -50)
si=(-1.1954154 -0.9960733)
p="Van-der-Grinten(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R95/455/-90/90 -Jv-85/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Sinusoidal, Spherical

sf=(-75 -50)
si=(0.1682814 -0.8726646)
p="Sinusoidal(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R90/450/-90/90 -Ji-90/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Mollweide, Spherical

sf=(-75 -50)
si=(0.1788845 -0.9208758)
p="Mollweide(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R90/450/-90/90 -Jw-90/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Eckert IV, Spherical

sf=(-75 -50)
si=(0.1875270 -0.9519210)
p="Eckert-IV(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R90/450/-90/90 -Jkf-90/1:1 >> fail

#----------------------------------------------------------------------------

#----------------------------------------------------------------------------
# Eckert VI, Spherical

sf=(-75 -50)
si=(0.1693623 -0.9570223)
p="Eckert-VI(Spherical)"
blabber --PROJ_ELLIPSOID=$unit --FORMAT_FLOAT_OUT=%9.7lf -R90/450/-90/90 -Jks-90/1:1 >> fail

#----------------------------------------------------------------------------