Bagaimana cara membuat Tissot Indicatrix yang akurat?

31

Sebuah Tissot Indicatrix adalah metode yang berguna untuk berkomunikasi sekilas jenis distorsi proyeksi yang diberikan rentan terhadap (pada gambar di bawah, masing-masing lingkaran merah menempati daerah yang sama). Saya telah diberitahu bahwa metode populer untuk menghasilkan TI memiliki masalah mereka sendiri, sampai-sampai kadang-kadang sangat tidak akurat.

Apa masalah dengan metode yang populer, dan apa cara yang paling benar untuk menghasilkan TI yang dapat diakses oleh rata-rata GIS Anda (ette)?

Mercator dan bola dunia dengan tissot

matt wilkie
sumber
2
pertanyaan bagus Saya juga ingin tahu itu.
George Silva

Jawaban:

23

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.

teks alternatif

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.

whuber
sumber
2
whuber adalah yang terbaik: P.
George Silva
1
Saya perhatikan bahwa ESRI memposting artikel tentang membuat tissot indicatrix dengan buffer, apakah metode buffer benar karena indicatrix dan buffer 'tidak sama sekali'? blogs.esri.com/Support/blogs/mappingcenter/archive/2011/03/28/…
SaultDon
1
@Sault Buffer adalah pengganti untuk TI yang benar. Mereka akan bekerja dengan cukup baik secara umum, asalkan mereka tidak meluas lebih dari beberapa derajat, kecuali jika distorsi menjadi ekstrim, di mana mereka tidak akan bekerja sama sekali (karena buffer tumpah melintasi wilayah distorsi tak terbatas, seperti yang ditunjukkan pada ESRI map ini ) atau hasilnya non-elips. Sedikit perubahan akan membuat pendekatan ini bekerja lebih baik: hitung buffer kecil, seperti buffer 50 m, dan rentangkan masing-masing secara seragam (di peta) di sekitar pusatnya agar terlihat.
Whuber
1
Bagaimana Anda menghitung orientasi elips?
Jason Davies
3
@Jason The Tissot indicatrix tidak lain adalah representasi bergambar dari turunan dari proyeksi yang diterapkan pada lingkaran di bidang tangen pada suatu titik. Maka, cara paling sederhana untuk menggambarnya adalah dengan menghitung turunan ini (ini adalah matriks 2 dengan 2 A ), membuat parameter lingkaran, dan menerapkan turunan ke simpul parameter. Voila, kamu memiliki elips. (Persis seperti bagaimana elips pada gambar saya digambar.) Kapaknya dapat dihitung sebagai vektor eigen dari A'A ; ini adalah perhitungan sederhana yang melibatkan tidak ada yang lebih rumit daripada rumus kuadratik.
whuber