Saya menggunakan ArcGIS Pro untuk mencoba mengidentifikasi kawah yang dibayangi secara permanen di Kutub Utara Bulan. Gambar di bawah adalah tangkapan layar DEM yang saya gunakan yang dikumpulkan dari laser altimeter pada Lunar Reconnaissance Orbiter. Ini tersedia untuk diunduh dari USGS di sini :
Saya dapat menggunakan alat Hillshade dengan sudut ketinggian rendah dan model_shadows dihidupkan untuk membuat raster seperti yang di bawah ini yang menunjukkan daerah gelap kawah untuk kondisi pencahayaan tertentu. Gambar di bawah ini menunjukkan hasil Hillshade dengan Azimuth 315 deg dan Altitude 5 deg.
Dalam dunia yang disederhanakan, saya bisa menjalankan hillshade ini untuk berbagai input Azimuth dan menjumlahkan raster yang dihasilkan. Pixel yang masih memiliki nilai 0 akan menjadi daerah yang dibayangi secara permanen. Namun ada satu masalah dengan rencana ini. Bulan tidak rata, terutama di atas wilayah yang begitu luas. DEM ini mencakup seluruh permukaan Bulan utara dari 75 derajat lintang dan bahkan lebih rendah di sudut-sudut. Karena itu, ketinggian matahari akan sangat bervariasi untuk berbagai piksel dalam gambar.
Untuk meringkas: Apakah ada cara untuk menghitung bayangan yang disebabkan oleh medan pada permukaan bola yang besar?
Alat Hillshade hanya beroperasi pada peta planar.
EDIT: Esri memiliki halaman bantuan di mana mereka pada dasarnya memberi Anda rumus hillshade dan menjelaskan semua persyaratan. Saya dapat memodifikasi rumus ini dengan menghitung sudut azimuth dan zenith yang unik untuk setiap piksel. Hasilnya di bawah ini hampir apa yang saya cari. Dalam gambar ini, matahari berada di bagian atas halaman dan mudah untuk melihat bagaimana kelengkungan bulan membuat sisi lain lebih gelap karena sudut zenith yang lebih tinggi. Satu-satunya hal yang hilang dari gambar ini adalah bayangan, itulah sebabnya beberapa dinding kawah pada gambar bawah masih menyala meskipun berada di sisi yang jauh dari cahaya. Sayangnya, saya tidak dapat menemukan dokumentasi tentang bagaimana bayangan dihitung.