2 # Script to do some testing of various transformation that do not depend
11 echo "Usage: ${0} <path to 'cs2cs' program>"
16 if test -z "${EXE}"; then
20 if test ! -x ${EXE}; then
21 echo "*** ERROR: Can not find '${EXE}' program!"
25 # Would be great to have a universale way of selecting a locale with
26 # a decimal separator that is not '.'
27 if command locale >/dev/null 2>/dev/null; then
28 if test `locale -a | grep fr_FR.utf8`; then
29 echo "Using locale with comma as decimal separator"
30 export LC_ALL=fr_FR.UTF-8
31 export PROJ_USE_ENV_LOCALE=1
35 echo "============================================"
36 echo "Running ${0} using ${EXE}:"
37 echo "============================================"
42 echo "doing tests into file ${OUT}, please wait"
45 echo "##############################################################" >> ${OUT}
46 echo Test raw ellipse to raw ellipse >> ${OUT}
48 $EXE +proj=latlong +ellps=clrk66 \
49 +to +proj=latlong +ellps=bessel \
51 79d58'00.000"W 37d02'00.000"N 0.0
52 79d58'00.000"W 36d58'00.000"N 0.0
55 echo "##############################################################" >> ${OUT}
56 echo Test NAD27 to raw ellipse >> ${OUT}
58 $EXE +proj=latlong +datum=NAD27 \
59 +to +proj=latlong +ellps=bessel \
61 79d00'00.000"W 35d00'00.000"N 0.0
64 echo "##############################################################" >> ${OUT}
65 echo Between two 3parameter approximations on same ellipsoid >> ${OUT}
67 $EXE +proj=latlong +ellps=bessel +towgs84=5,0,0 \
68 +to +proj=latlong +ellps=bessel +towgs84=1,0,0 \
70 0d00'00.000"W 0d00'00.000"N 0.0
71 79d00'00.000"W 45d00'00.000"N 0.0
74 echo "##############################################################" >> ${OUT}
75 echo 3param to raw ellipsoid on same ellipsoid >> ${OUT}
77 $EXE +proj=latlong +ellps=bessel +towgs84=5,0,0 \
78 +to +proj=latlong +ellps=bessel \
80 0d00'00.000"W 0d00'00.000"N 0.0
81 79d00'00.000"W 45d00'00.000"N 0.0
84 echo "##############################################################" >> ${OUT}
85 echo Test simple prime meridian handling. >> ${OUT}
87 $EXE +proj=latlong +datum=WGS84 +pm=greenwich \
88 +to +proj=latlong +datum=WGS84 +pm=1 \
90 0d00'00.000"W 0d00'00.000"N 0.0
91 79d00'00.000"W 45d00'00.000"N 0.0
94 echo "##############################################################" >> ${OUT}
95 echo Test support for the lon_wrap switch. >> ${OUT}
97 $EXE +proj=latlong +datum=WGS84 \
98 +to +proj=latlong +datum=WGS84 +lon_wrap=180 \
100 1d00'00.000"W 10d00'00.000"N 0.0
101 0d00'00.000"W 10d00'00.000"N 0.0
102 0d00'00.000"E 10d00'00.000"N 0.0
103 1d00'00.000"E 45d00'00.000"N 0.0
104 179d00'00.000"E 45d00'00.000"N 0.0
105 181d00'00.000"E 45d00'00.000"N 0.0
106 350d00'00.000"E 45d00'00.000"N 0.0
107 370d00'00.000"E 45d00'00.000"N 0.0
110 echo "##############################################################" >> ${OUT}
111 echo Test simple prime meridian handling within a projection. >> ${OUT}
113 $EXE +proj=utm +zone=11 +datum=WGS84 +pm=3 \
114 +to +proj=latlong +datum=WGS84 +pm=1w \
118 echo "##############################################################" >> ${OUT}
119 echo Test geocentric x/y/z generation. >> ${OUT}
121 $EXE +proj=latlong +datum=WGS84 \
122 +to +proj=geocent +datum=WGS84 \
124 0d00'00.001"W 0d00'00.001"N 0.0
125 0d00'00.001"W 0d00'00.001"N 10.0
126 79d00'00.000"W 45d00'00.000"N 0.0
127 45d00'00.000"W 89d59'59.990"N 0.0
130 echo "##############################################################" >> ${OUT}
131 echo Test geocentric x/y/z consumption. >> ${OUT}
133 $EXE +proj=geocent +datum=WGS84 \
134 +to +proj=latlong +datum=WGS84 \
136 6378137.00 -0.00 0.00
137 6378147.00 -0.00 0.00
138 861996.98 -4434590.01 4487348.41
139 0.00 -0.00 6356752.31
142 #echo "#############################################################" >> ${OUT}
143 #echo Test conversion between geocentric latlong and geodetic latlong >> ${OUT}
145 # The +geoc flag does not currently work with pj_transform() so this is
148 #$EXE +proj=latlong +datum=WGS84 +geoc \
149 # +to +proj=latlong +datum=WGS84 \
151 #0d00'00.000"W 0d00'00.000"N 0.0
152 #79d00'00.000"W 45d00'00.000"N 0.0
153 #12d00'00.000"W 45d00'00.000"N 0.0
154 #0d00'00.000"W 90d00'00.000"N 0.0
157 echo "##############################################################" >> ${OUT}
158 echo "Test stere projection (re: win32 ticket 12)" >> ${OUT}
160 $EXE +proj=latlong +datum=WGS84 \
161 +to +proj=stere +lat_0=90 +lon_0=0 +lat_ts=70 +datum=WGS84 \
166 echo "##############################################################" >> ${OUT}
167 echo "Test stere without lat_ts (#147)" >> ${OUT}
169 $EXE +proj=latlong +datum=WGS84 \
170 +to +proj=stere +lat_0=40 +lon_0=10 +datum=WGS84 \
175 echo "##############################################################" >> ${OUT}
176 echo "Test sts projection (re: ticket 12)" >> ${OUT}
178 $EXE +proj=latlong +ellps=WGS84 \
179 +to +proj=kav5 +ellps=WGS84 +units=m \
183 $EXE +proj=kav5 +ellps=WGS84 +units=m \
184 +to +proj=latlong +ellps=WGS84 \
186 383646.088858 5997047.888175
189 echo "##############################################################" >> ${OUT}
190 echo "Test RSO Borneo projection (re: ticket 62)" >> ${OUT}
192 $EXE +proj=latlong +a=6377298.556 +rf=300.8017 \
193 +to +proj=omerc +a=6377298.556 +rf=300.8017 +lat_0=4 +lonc=115 \
194 +alpha=53d18\'56.9537 +gamma=53d7\'48.3685 +k_0=0.99984 \
195 +x_0=590476.87 +y_0=442857.65 \
197 116d2'11.12630 5d54'19.90183
199 echo "##############################################################" >> ${OUT}
200 echo "Test extended transverse mercator (#97)" >> ${OUT}
202 $EXE +proj=etmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
203 +to +proj=latlong +datum=WGS84 \
212 echo "##############################################################" >> ${OUT}
213 echo "Test extended transverse mercator inverse (#97)" >> ${OUT}
215 $EXE +proj=latlong +datum=WGS84 \
216 +to +proj=etmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
219 15d22'16.108"W 17d52'53.478"N 0.000
220 10d40'55.532"W 17d42'48.526"N 0.000
221 1d32'21.33"W 17d3'47.233"N 0.000
222 15d4'42.357"E 14d48'56.372"N 0.000
225 echo "##############################################################" >> ${OUT}
226 echo "Test transverse mercator (#97)" >> ${OUT}
228 $EXE +proj=tmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
229 +to +proj=latlong +datum=WGS84 \
238 echo "##############################################################" >> ${OUT}
239 echo "Test transverse mercator inverse (#97)" >> ${OUT}
241 $EXE +proj=latlong +datum=WGS84 \
242 +to +proj=tmerc +k=0.998 +lon_0=-20 +datum=WGS84 +x_0=10000 +y_0=20000 \
245 15d22'16.108"W 17d52'53.478"N 0.000
246 10d40'55.532"W 17d42'48.526"N 0.000
247 1d32'21.33"W 17d3'47.233"N 0.000
248 15d4'42.357"E 14d48'56.372"N 0.000
251 echo "##############################################################" >> ${OUT}
252 echo "Test robinson projection (#113)" >> ${OUT}
254 $EXE +proj=latlong +datum=WGS84 \
255 +to +proj=robin +datum=WGS84 \
261 $EXE +proj=robin +datum=WGS84 \
262 +to +proj=latlong +datum=WGS84 \
264 -2612095.95 4276351.58 0.00
265 -2963455.42 4805073.65 0.00
266 1741397.30 4276351.58 0.00
268 echo "##############################################################" >> ${OUT}
269 echo "Test hammer projection (pull request #329)" >> ${OUT}
270 $EXE +proj=latlong +datum=WGS84 \
271 +to +proj=hammer +datum=WGS84 \
277 $EXE +proj=hammer +datum=WGS84 \
278 +to +proj=latlong +datum=WGS84 \
280 -2711575.08 4395506.62 0.00
281 -2964412.70 4929091.33 0.00
282 1811748.54 4377349.50 0.00
284 echo "##############################################################" >> ${OUT}
285 echo "Test healpix forward projection on sphere" >> ${OUT}
286 $EXE +proj=latlong +a=1 +lon_0=0 \
287 +to +proj=healpix +a=1 +lon_0=0 -f '%.'5'f' \
292 $EXE +proj=latlong +R=5 \
293 +to +proj=healpix +R=5 -f '%.'5'f' \
297 0 -41.810314895778596
306 echo "Test healpix forward projection on ellipsoid" >> ${OUT}
307 $EXE +proj=latlong +a=1 +lon_0=0 +ellps=WGS84 \
308 +to +proj=healpix +a=1 +lon_0=0 +ellps=WGS84 -f '%.'5'f' \
313 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
314 +to +proj=healpix +a=5 +e=0.8 +r_a=4.3220011711888882 -f '%.'5'f' \
318 0 -41.810314895778596
327 echo "Test healpix inverse projection on ellipsoid" >> ${OUT}
328 $EXE +proj=latlong +a=1 +lon_0=0 +ellps=WGS84 \
329 +to +proj=healpix +a=1 +lon_0=0 +ellps=WGS84 -f '%.'5'f' -I\
332 -1.5707963267948966 0
334 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
335 +to +proj=healpix +a=5 +e=0.8 +r_a=4.3220011711888882 -f '%.'5'f' -I\
338 0.0 2.0547874222147415
339 0.0 -2.0547874222147415
340 6.788983564106746 0.0
341 -6.788983564106746 0.0
342 -13.577967128213492 0.0
343 -10.183475346160119 6.788983564106746
344 -10.183475346160119 -6.788983564106746
345 0.0 3.351278550178025
346 0.0 -3.351278550178025
349 echo "##############################################################" >> ${OUT}
350 echo "Test rHEALPix forward projection on sphere north=0 south=0" >> ${OUT}
351 $EXE +proj=latlong +R=5 \
352 +to +proj=rhealpix +R=5 +north_square=0 +south_square=0 -f '%.'5'f' \
355 -180 -25.714285714285715
357 60.0 41.809314895778598
360 echo "##############################################################" >> ${OUT}
361 echo "Test rHEALPix forward projection on sphere north=1 south=1" >> ${OUT}
362 $EXE +proj=latlong +R=5 \
363 +to +proj=rhealpix +R=5 +north_square=1 +south_square=1 -f '%.'5'f' \
366 -180 -25.714285714285715
368 60.0 41.809314895778598
371 echo "##############################################################" >> ${OUT}
372 echo "Test rHEALPix inverse projection on sphere north=0 south=0" >> ${OUT}
373 $EXE +proj=latlong +R=5 \
374 +to +proj=rhealpix +R=5 +north_square=0 +south_square=0 -f '%.'5'f' -I\
377 0.0 3.9269908169872414
378 0.0 -3.9269908169872414
379 7.853981633974483 0.0
380 -7.853981633974483 0.0
383 echo "##############################################################" >> ${OUT}
384 echo "Test rHEALPix inverse projection on sphere north=1 south=1" >> ${OUT}
385 $EXE +proj=latlong +R=5 \
386 +to +proj=rhealpix +R=5 +north_square=1 +south_square=1 -f '%.'5'f' -I\
389 0.0 3.9269908169872414
390 0.0 -3.9269908169872414
391 7.853981633974483 0.0
392 -7.853981633974483 0.0
395 echo "##############################################################" >> ${OUT}
396 echo "Test rHEALPix forward projection on ellipsoid north=0 south=0" >> ${OUT}
397 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
398 +to +proj=rhealpix +a=5 +e=0.8 +r_a=4.3220011711888882 +north_square=0 +south_square=0 -f '%.'5'f' \
402 0 -41.810314895778596
407 echo "##############################################################" >> ${OUT}
408 echo "Test rHEALPix forward projection on ellipsoid north=1 south=1" >> ${OUT}
409 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
410 +to +proj=rhealpix +a=5 +e=0.8 +r_a=4.3220011711888882 +north_square=1 +south_square=1 -f '%.'5'f' \
414 0 -41.810314895778596
419 echo "##############################################################" >> ${OUT}
420 echo "Test rHEALPix inverse projection on ellipsoid north=0 south=0" >> ${OUT}
421 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
422 +to +proj=rhealpix +a=5 -I +e=0.8 +r_a=4.3220011711888882 +north_square=0 +south_square=0 -f '%.'5'f'\
425 0.0 2.0547874222147415
426 0.0 -2.0547874222147415
427 6.788983564106746 0.0
428 -6.788983564106746 0.0
431 echo "##############################################################" >> ${OUT}
432 echo "Test rHEALPix inverse projection on ellipsoid north=1 south=1" >> ${OUT}
433 $EXE +proj=latlong +a=5 +e=0.8 +r_a=4.3220011711888882\
434 +to +proj=rhealpix +a=5 -I +e=0.8 +r_a=4.3220011711888882 +north_square=1 +south_square=1 -f '%.'5'f'\
437 0.0 2.0547874222147415
438 0.0 -2.0547874222147415
439 6.788983564106746 0.0
440 -6.788983564106746 0.0
443 echo "##############################################################" >> ${OUT}
444 echo "Test geos projection" >> ${OUT}
445 echo "Test geos on a sphere" >> ${OUT}
446 $EXE +proj=latlong +ellps=sphere \
447 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere -E >>${OUT} <<EOF
449 -43d11'47" -22d54'30"
453 echo "Test geos on a ellipsoid" >> ${OUT}
454 $EXE +proj=latlong +ellps=sphere \
455 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 -E >>${OUT} <<EOF
457 -43d11'47" -22d54'30"
461 echo "Test inv geos on a sphere" >> ${OUT}
462 $EXE +proj=latlong +ellps=sphere \
463 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere -I -E >>${OUT} <<EOF
465 -3780930.93 -2326595.36
466 1608689.65 -3412115.56
467 3825202.59 2885980.79
469 echo "Test inv geos on a ellipsoid" >> ${OUT}
470 $EXE +proj=latlong +ellps=sphere \
471 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 -I -E >>${OUT} <<EOF
473 -3787026.57 -2314765.32
474 1612331.00 -3397031.37
475 3832522.65 2872185.29
477 echo "Test geos on a sphere with alternate sweep" >> ${OUT}
478 $EXE +proj=latlong +ellps=sphere \
479 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere +sweep=x -E >>${OUT} <<EOF
481 -43d11'47" -22d54'30"
485 echo "Test geos on a ellipsoid with alternate sweep" >> ${OUT}
486 $EXE +proj=latlong +ellps=sphere \
487 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 +sweep=x -E >>${OUT} <<EOF
489 -43d11'47" -22d54'30"
493 echo "Test inv geos on a sphere with alternate sweep" >> ${OUT}
494 $EXE +proj=latlong +ellps=sphere \
495 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=sphere +sweep=x -I -E >>${OUT} <<EOF
497 -3772913.22 -2339604.71
498 1601377.77 -3415545.15
499 3812722.89 2902474.62
501 echo "Test inv geos on a ellipsoid with alternate sweep" >> ${OUT}
502 $EXE +proj=latlong +ellps=sphere \
503 +to +proj=geos +h=35785831.0 +lon_0=0 +ellps=WGS84 +sweep=x -I -E >>${OUT} <<EOF
505 -3779077.27 -2327750.87
506 1605067.15 -3400461.47
507 3820138.08 2888664.15
510 echo "##############################################################" >> ${OUT}
511 echo "Test the Natural Earth Projection" >> ${OUT}
512 $EXE +proj=latlong +a=1 +lon_0=0 \
513 +to +proj=natearth +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
516 0.0 22.5 0 0.0 2525419.569383768
517 0.0 45.0 0 0.0 5052537.389973222
518 0.0 67.5 0 0.0 7400065.6562573705
519 0.0 90.0 0 0.0 9062062.394736718
520 45.0 0.0 0 4356790.016612169 0.0
521 45.0 22.5 0 4253309.544984069 2525419.569383768
522 45.0 45.0 0 3924521.5829515466 5052537.389973222
523 45.0 67.5 0 3354937.47115583 7400065.6562573705
524 45.0 90.0 0 2397978.2448443635 9062062.394736718
525 90.0 0.0 0 8713580.033224339 0.0
526 90.0 22.5 0 8506619.089968137 2525419.569383768
527 90.0 45.0 0 7849043.165903093 5052537.389973222
528 90.0 67.5 0 6709874.94231166 7400065.6562573705
529 90.0 90.0 0 4795956.489688727 9062062.394736718
530 135.0 0.0 0 1.3070370049836507E7 0.0
531 135.0 22.5 0 1.2759928634952208E7 2525419.569383768
532 135.0 45.0 0 1.177356474885464E7 5052537.389973222
533 135.0 67.5 0 1.0064812413467491E7 7400065.6562573705
534 135.0 90.0 0 7193934.734533091 9062062.394736718
535 180.0 0.0 0 1.7427160066448677E7 0.0
536 180.0 22.5 0 1.7013238179936275E7 2525419.569383768
537 180.0 45.0 0 1.5698086331806187E7 5052537.389973222
538 180.0 67.5 0 1.341974988462332E7 7400065.6562573705
539 180.0 90.0 0 9591912.979377454 9062062.394736718
541 echo "##############################################################" >> ${OUT}
542 echo "Test the Natural Earth II Projection" >> ${OUT}
543 $EXE +proj=latlong +a=1 +lon_0=0 \
544 +to +proj=natearth2 +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
546 0.0 0.0 0 0.00000000 0.00000000
547 0.0 22.5 0 0.00000000 2531453.57080958
548 0.0 45.0 0 0.00000000 5051471.50086845
549 0.0 67.5 0 0.00000000 7395411.22478983
550 0.0 90.0 0 0.00000000 9073776.52662810
551 45.0 0.0 0 4239151.18200719 0.00000000
552 45.0 22.5 0 4138348.61904244 2531453.57080958
553 45.0 45.0 0 3830621.33880773 5051471.50086845
554 45.0 67.5 0 3158326.32836996 7395411.22478983
555 45.0 90.0 0 957973.37034235 9073776.52662810
556 90.0 0.0 0 8478302.36401439 0.00000000
557 90.0 22.5 0 8276697.23808488 2531453.57080958
558 90.0 45.0 0 7661242.67761547 5051471.50086845
559 90.0 67.5 0 6316652.65673992 7395411.22478983
560 90.0 90.0 0 1915946.74068470 9073776.52662810
561 135.0 0.0 0 12717453.54602160 0.00000000
562 135.0 22.5 0 12415045.85712730 2531453.57080958
563 135.0 45.0 0 11491864.01642320 5051471.50086845
564 135.0 67.5 0 9474978.98510988 7395411.22478983
565 135.0 90.0 0 2873920.11102705 9073776.52662810
566 180.0 0.0 0 16956604.72802880 0.00000000
567 180.0 22.5 0 16553394.47616980 2531453.57080958
568 180.0 45.0 0 15322485.35523090 5051471.50086845
569 180.0 67.5 0 12633305.31347990 7395411.22478983
570 180.0 90.0 0 3831893.48136940 9073776.52662810
572 echo "##############################################################" >> ${OUT}
573 echo "Test the Compact Miller projection" >> ${OUT}
574 $EXE +proj=latlong +a=1 +lon_0=0 \
575 +to +proj=comill +a=6371008.7714 +b=6371008.7714 -f '%.'7'f' \
578 0.0 22.5 0 0.0 2537439.6610749415
579 0.0 45.0 0 0.0 5391682.432264133
580 0.0 67.5 0 0.0 8661480.510260897
581 0.0 90.0 0 0.0 12009484.264916677
582 45.0 0.0 0 5003778.588046594 0.0
583 45.0 22.5 0 5003778.588046594 2537439.6610749415
584 45.0 45.0 0 5003778.588046594 5391682.432264133
585 45.0 67.5 0 5003778.588046594 8661480.510260897
586 45.0 90.0 0 5003778.588046594 12009484.264916677
587 90.0 0.0 0 10007557.176093187 0.0
588 90.0 22.5 0 10007557.176093187 2537439.6610749415
589 90.0 45.0 0 10007557.176093187 5391682.432264133
590 90.0 67.5 0 10007557.176093187 8661480.510260897
591 90.0 90.0 0 10007557.176093187 12009484.264916677
592 135.0 0.0 0 15011335.76413978 0.0
593 135.0 22.5 0 15011335.76413978 2537439.6610749415
594 135.0 45.0 0 15011335.76413978 5391682.432264133
595 135.0 67.5 0 15011335.76413978 8661480.510260897
596 135.0 90.0 0 15011335.76413978 12009484.264916677
597 180.0 0.0 0 20015114.352186374 0.0
598 180.0 22.5 0 20015114.352186374 2537439.6610749415
599 180.0 45.0 0 20015114.352186374 5391682.432264133
600 180.0 67.5 0 20015114.352186374 8661480.510260897
601 180.0 90.0 0 20015114.352186374 12009484.264916677
603 echo "##############################################################" >> ${OUT}
604 echo "Test pconic (#148)" >> ${OUT}
606 $EXE +proj=latlong +datum=WGS84 \
607 +to +proj=pconic +units=m +lat_1=20n +lat_2=60n +lon_0=60W +datum=WGS84 \
611 $EXE +proj=pconic +units=m +lat_1=20n +lat_2=60n +lon_0=60W +datum=WGS84 \
612 +to +proj=latlong +datum=WGS84 \
614 -2240096.40 -6940342.15
616 echo "##############################################################" >> ${OUT}
617 echo "Test laea" >> ${OUT}
620 +proj=laea +lat_0=45 +lon_0=-100 +units=m +datum=WGS84 +no_defs \
621 +to +proj=latlong +datum=WGS84 \
625 echo "##############################################################" >> ${OUT}
626 echo "Test forward calcofi projection" >> ${OUT}
627 $EXE +proj=latlong +ellps=clrk66 \
628 +to +proj=calcofi +ellps=clrk66 \
630 120d40'42.273"W 38d56'50.766"N
632 123d59'56.066"W 30d25'4.617"N
634 echo "Test inverse calcofi projection" >> ${OUT}
635 $EXE +proj=calcofi +ellps=clrk66 \
636 +to +proj=longlat +ellps=clrk66 \
642 echo "##############################################################" >> ${OUT}
643 echo "Check inverse error handling with ob_tran (#225)" >> ${OUT}
645 +o_proj=moll +o_lon_p=LON_POLE +o_lat_p=LAT_POLE +lon_0=180 +ellps=WGS84 \
650 echo "Test inverse handling" >> ${OUT}
651 $EXE -I +proj=ob_tran \
652 +o_proj=moll +o_lon_p=LON_POLE +o_lat_p=LAT_POLE +lon_0=180 +ellps=WGS84 \
656 echo "##############################################################" >> ${OUT}
657 echo "Test MGI datum gives expected results (#207)" >> ${OUT}
660 +proj=latlong +datum=WGS84 \
661 +to +init=epsg:31284 \
665 echo "##############################################################" >> ${OUT}
666 echo "Test omerc sensitivity with locations 90d from origin(#114)" >> ${OUT}
669 +proj=latlong +ellps=WGS84 \
670 +to +proj=omerc +ellps=WGS84 +lon_1=62.581150 +lat_1=74.856102 \
671 +lon_2=53.942810 +lat_2=74.905884 +units=km +no_rot +no_defs \
673 56.958381652832 72.8798
676 echo "##############################################################" >> ${OUT}
677 echo "Test omerc differences between poles (#190)" >> ${OUT}
680 +proj=latlong +ellps=WGS84 \
681 +to +proj=omerc +ellps=WGS84 +datum=WGS84 +no_rot \
682 +lon_1=-27 +lat_1=70 +lon_2=-38 +lat_2=80 +lat_0=70 \
690 # Again, for the south pole.
692 +proj=latlong +ellps=WGS84 \
693 +to +proj=omerc +ellps=WGS84 +datum=WGS84 +no_rot \
694 +lon_1=-27 +lat_1=-70 +lon_2=-38 +lat_2=-80 +lat_0=-70 \
702 echo "##############################################################" >> ${OUT}
703 echo "Test qsc" >> ${OUT}
706 +proj=latlong +datum=WGS84 \
707 +to +proj=qsc +datum=WGS84 +no_defs \
712 +proj=qsc +datum=WGS84 +no_defs \
713 +to +proj=latlong +datum=WGS84 \
715 2073986.94908809568733 -1680858.27222427958623
717 echo "##############################################################" >> ${OUT}
718 echo "Test bug 229" >> ${OUT}
721 +init=epsg:4326 +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0 \
722 +to +proj=latlong +datum=WGS84 +no_defs \
726 echo "##############################################################" >> ${OUT}
727 echo "Test bug 229 (2)" >> ${OUT}
730 +init=epsg:4326 +to +init=epsg:4326 \
734 echo "##############################################################" >> ${OUT}
735 echo "Test bug 244 " >> ${OUT}
739 +to +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \
741 -140.100000 -87.000000
743 echo "##############################################################" >> ${OUT}
744 echo "Test bug 244 (2)" >> ${OUT}
747 +proj=aeqd +lon_0=130.0 +lat_0=40.0 +a=6378137 +b=6378137 +units=m +no_defs \
748 +to +init=epsg:4326 \
750 987122.418330284 -14429896.539530909
752 echo "##############################################################" >> ${OUT}
753 echo "Test bug 245 (use +datum=carthage)" >> ${OUT}
756 +proj=longlat +datum=WGS84 +no_defs +to +proj=utm +zone=32 +datum=carthage +no_defs \
760 echo "##############################################################" >> ${OUT}
761 echo "Test bug 245 (use expansion of +datum=carthage)" >> ${OUT}
764 +proj=longlat +datum=WGS84 +no_defs +to +proj=utm +zone=32 +a=6378249.2 +b=6356515 +towgs84=-263.0,6.0,431.0 +no_defs \
768 echo "##############################################################" >> ${OUT}
769 echo "Test SCH forward projection" >> ${OUT}
772 +proj=latlong +datum=WGS84 +to +proj=sch +datum=WGS84 +plat_0=30.0 +plon_0=45.0 \
773 +phdg_0=-12.0 +nodefs \
782 echo "##############################################################" >> ${OUT}
783 echo "Test SCH inverse projection" >> ${OUT}
786 +proj=sch +datum=WGS84 +plat_0=30.0 +plon_0=45.0 +phdg_0=-12.0 +nodefs +to \
787 +proj=latlong +datum=WGS84 \
794 ##############################################################################
795 echo "##############################################################" >> ${OUT}
796 echo "Test issue #316 (switch utm to use etmerc)" >> ${OUT}
799 +proj=latlong +datum=WGS84 +no_defs +to +proj=utm +zone=35 +datum=WGS84 +no_defs \
803 ##############################################################################
804 echo "##############################################################" >> ${OUT}
805 echo "Test issue #316 (switch utm to use etmerc)" >> ${OUT}
808 +proj=latlong +datum=WGS84 +no_defs +to +proj=etmerc +datum=WGS84 +no_defs +k=0.9996 +lon_0=27 +x_0=500000 -f %.6f \
812 ##############################################################################
813 echo "##############################################################" >> ${OUT}
814 echo "Test nzmg forward projection" >> ${OUT}
817 +proj=latlong +datum=WGS84 +to \
818 +proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=WGS84 +units=m +no_defs \
822 ##############################################################################
823 echo "##############################################################" >> ${OUT}
824 echo "Test nzmg inverse projection" >> ${OUT}
827 +proj=nzmg +lat_0=-41 +lon_0=173 +x_0=2510000 +y_0=6023150 +ellps=WGS84 +units=m +no_defs +to \
828 +proj=latlong +datum=WGS84 \
830 2680778.57267967 6132228.07645127 0.
832 ##############################################################################
833 echo "##############################################################" >> ${OUT}
834 echo "Test misrsom forward projection" >> ${OUT}
837 +proj=latlong +datum=WGS84 +no_defs +to \
838 +proj=misrsom +path=41 +ellps=WGS84 +units=m \
840 48.64966165540372 66.2263195368941 0.
842 ##############################################################################
843 echo "##############################################################" >> ${OUT}
844 echo "Test misrsom inverse projection" >> ${OUT}
847 +proj=misrsom +path=41 +ellps=WGS84 +units=m +no_defs +to \
848 +proj=latlong +datum=WGS84 -f '%0.7f' \
850 7461300.0 528000.0 0.0
852 echo "##############################################################" >> ${OUT}
853 echo "Test patterson forward projection" >> ${OUT}
856 +proj=latlong +datum=WGS84 \
857 +to +proj=patterson +a=6371008.7714 +b=6371008.7714 +units=m \
869 echo "##############################################################" >> ${OUT}
870 echo "Test patterson inverse projection" >> ${OUT}
873 +proj=patterson +a=6371008.7714 +b=6371008.7714 +units=m \
874 +to +proj=latlong +datum=WGS84 \
876 -20015114.352186374 11409566.822831295
877 -15011335.76413978 8729502.054111844
878 -10007557.176093187 5366413.421153781
879 -5003778.588046594 2551415.729669344
881 5003778.588046594 -2551415.729669344
882 10007557.176093187 -5366413.421153781
883 15011335.76413978 -8729502.054111844
884 20015114.352186374 -11409566.822831295
889 # do 'diff' with distribution results
890 echo "diff ${OUT} with ${OUT}.dist"
891 diff -b ${OUT} ${NAD_DIR}/${OUT}.dist
892 if [ $? -ne 0 ] ; then
894 echo "PROBLEMS HAVE OCCURED"
895 echo "test file ${OUT} saved"
900 echo "test file ${OUT} removed"