Setiap perangkat lunak yang dapat memproyeksikan koordinat secara akurat dapat menghitung pendeteksian Tissot yang akurat .
Sumber yang baik untuk formula adalah Snyder, John, Proyeksi Peta - Manual Kerja , terutama di hlm. 20-26. (Saya tidak akan mereproduksi mereka di sini karena situs ini tidak memiliki alat yang tepat untuk mengkomunikasikan rumus matematika.) Mereka memerlukan keempat turunan pertama dari koordinat yang diproyeksikan (x, y) sehubungan dengan koordinat bola (lat, lon) = (phi, lambda):
dx / d(phi), dx / d(lambda);
dy / d(phi), dy / d(lambda).
Segala sesuatu yang lain tentang TI dihitung dalam hal ini (menggunakan beberapa fungsi aritmatika dan trigonometri: cosinus, sinus invers utama, dan tangen invers utama). Perhitungannya membutuhkan deskripsi bentuk bumi. Untuk akurasi terbesar, gunakan datum ellipsoidal dengan sumbu semimajor a dan eksentrisitas e. (Ini akan diketahui oleh perangkat lunak.)
Buku Snyder memiliki instruksi tentang cara menghitung semuanya kecuali turunan ini. Lakukan secara numerik. Saya telah mendapatkan hasil yang sangat baik dengan menggunakan estimasi selisih hingga pusat orde pertama pada jarak h = 10 ^ (- 5.2) radian (biasanya sekitar 50 meter): ini adalah kompromi yang baik antara mencoba untuk menjadi sangat dekat dan kehilangan terlalu banyak presisi dari roundoff floating point (dengan asumsi presisi ganda), karena kesalahan yang dibuat sebanding dengan (10 ^ (- 5.2)) ^ 2 = 10 ^ (- 10.4) dan 10 ^ (- 5.2) sama dengan 10 ^ 10.4 kali akurasi presisi ganda IEEE dari 10 ^ (- 15.6) dan masih jauh lebih besar dari ketepatan tipikal dalam proyeksi, yang biasanya berjalan dari 10 ^ (- 10) menjadi sekitar 10 ^ (- 14).
Jadi, bagaimana Anda menghitung perkiraan perbedaan hingga? Bagian ini sangat mudah. Untuk mendapatkan dx / d (phi) pada suatu titik (phi, lambda), minta GIS Anda untuk memproyeksikan poinnya
(phi - h/2, lambda) --> (x0,y0),
(phi + h/2, lambda) --> (x1,y1).
Gunakan taksiran
dx / d(phi) = (x1 - x0)/h,
dy / d(phi) = (y1 - y0)/h.
Demikian pula, proyeksikan poinnya
(phi, lambda - h/2) --> (x2,y2),
(phi, lambda + h/2) --> (x3,y3)
dan gunakan estimasi
dx / d(lambda) = (x3 - x2)/h,
dy / d(lambda) = (y3 - y2)/h.
Itu membutuhkan empat proyeksi dan sedikit aritmatika. (Anda dapat menguranginya menjadi tiga dengan menggunakan perbedaan non-sentral, tetapi akurasinya sedikit turun. Ini adalah tujuan bijak untuk akurasi tinggi, tanpa membiarkan h menjadi terlalu kecil, kecuali Anda yakin GIS Anda menggunakan survey-grade (milimeter) keakuratan dalam rumus proyeksi.)
Dari turunan ini, bersama dengan rumus Snyder (memperhatikan modifikasi yang dijelaskan pada 4-19 dan 4-21), Anda dapat memperoleh panjang sumbu sumbu Tissot Indicatrix di (phi, lambda) dan orientasinya. Pada peta skala dunia TI akan sangat kecil sehingga tidak terlihat, jadi hal terakhir yang harus dilakukan adalah memutuskan seberapa banyak Anda ingin mengubah skala setiap TI. Saya menentukan faktor skala dengan mencari tahu seberapa besar peta akan, menemukan ukuran tipikal TI di seluruh peta, dan menskalakan sehingga TI tersebut akan sekitar 6% selebar peta. Ini awal yang baik; Saya membiarkan pengguna menyesuaikan ukuran TI dari sana. Tentu saja Anda akan menskala ulang semua TI dengan jumlah yang sama, sehingga mereka dapat dibandingkan, dan masing-masing akan ditata ulang di sekitar pusatnya sendiri (yang diperoleh dengan proyeksi kelima, (phi, lambda) -> (x, y) ).
Tambahan bagus untuk penggambaran elips TI adalah untuk menunjukkan arah meridian lokal dan paralel: kemudian, sekilas, Anda dapat menilai konvergensi kisi . Saya juga menunjukkan konsentris lingkaran standar (mewakili tidak ada distorsi) dengan masing-masing TI karena meningkatkan kemampuan pembaca untuk mengukur jumlah distorsi yang diwakili oleh setiap elips.
Yang perlu diperhatikan dalam proyeksi Mollweide ini adalah TI ekstrim di dekat kutub selatan. Itu masih elips yang sempurna dan secara akurat menggambarkan distorsi peta di sana.
