Menghitung nilai piksel berturut-turut untuk set raster menggunakan ArcGIS Spatial Analyst?

Saya menggunakan ArcGIS 10 dengan Spatial Analyst.

Saya memiliki satu set raster (total 8) yang hanya setiap berisi 1 atau 0 untuk setiap sel. Setiap raster mewakili nilai data tahun yang berbeda. Demi argumen tahun 1 hingga tahun 8.

Saya dapat menambahkan semua raster yang akan memberi saya kisi terakhir dengan nilai mulai dari 0 hingga 8. Angka 8 menunjukkan bahwa sel secara konstan 1 untuk himpunan raster (semua tahun).

Saya ingin mencari tahu untuk masing-masing sel nomor 1 terpanjang berturut-turut.

Jadi misalnya total grid dapat merekam untuk satu sel nilai katakanlah 5 tetapi lebih dari 8 grid sel itu memiliki jumlah berturut-turut 1 sama dengan 3. Atau cara lain untuk menyatakan ini adalah selama 3 tahun sel itu adalah 1 kemudian mulai berosilasi antara nol dan satu.

Keterampilan pemrosesan raster saya tidak sepanas keterampilan pemrosesan vektor saya dan saya telah melihat dengan baik pada file bantuan ESRI tapi saya tidak tahu bagaimana seseorang akan mencapai ini dengan menggunakan alat geo-processing dari rak?

Ada ide?

Karena ini adalah operasi lokal, mari kita cari tahu bagaimana melakukannya untuk sel tunggal: Peta Aljabar akan mengurus sisanya.

Catatan pertama bahwa urutan raster jelas penting. Karenanya, statistik sel yang dipotret sekali pakai, seperti jumlah sel, tidak akan melakukannya.

Jika kami menemukan urutan seperti 01110101 pada sel tertentu, kami akan memprosesnya dari awal hingga akhir dan

1. Mulailah dengan hitungan nol.

2. Bertambah hitungan setiap kali kita menjumpai 1.

3. Setel ulang penghitungan setiap kali kami menemukan angka 0, setelah menyimpan penghitungan terakhir .

4. Pada akhirnya, ambil jumlah maksimum yang disimpan (termasuk jumlah akhir).

Langkah 1 diimplementasikan dengan kisi nol konstan. Langkah 2 dan 3 bergantung pada apa yang kita temui: karena itu operasi bersyarat . Langkah 4 jelas adalah maksimum lokal. Kami akan mengkodekan ini, kemudian, sedikit lebih formal sebagai:

count = 0
result = 0
For each value:
    If (value==1):
        count=count+1
    else
        result = max(result, count)
        count=0
result = max(result, count)

Itu paling baik dilakukan dengan skrip Python ketika Anda memiliki banyak kotak, tetapi dengan delapan itu tidak sulit untuk membuka gulungan dan menulis langkah-langkahnya dengan tangan. Ini mengungkapkan sedikit masalah: result=max(longest,count)ini adalah semacam "efek samping" yang sulit dikodekan dengan operasi raster. (Tapi itu bisa dilakukan, seperti yang ditunjukkan pada solusi kedua di bawah ini.) Ini juga tidak efisien, karena menambah perhitungan ekstra pada setiap langkah. Karena itu kami memodifikasi pendekatan sedikit, dengan tujuan untuk menunda maxoperasi sampai akhir. Ini akan membutuhkan menyimpan hitungan terpisah di setiap tahap.

Dalam menjalani proses ini saya juga menemukan jalan pintas untuk langkah pertama. Ini mengarah ke solusi berikut, yang walaupun sedikit panjang dan intensif-RAM, sederhana dan melibatkan langkah-langkah yang dijalankan dengan cepat:

result1 = "grid1"
result2 = con("grid2"==1, "result1"+1, 0)
result3 = con("grid3"==1, "result2"+1, 0)
result4 = con("grid4"==1, "result3"+1, 0)
result5 = con("grid5"==1, "result4"+1, 0)
result6 = con("grid6"==1, "result5"+1, 0)
result7 = con("grid7"==1, "result6"+1, 0)
result8 = con("grid8"==1, "result7"+1, 0)
CellStatistics(["result1", "result2", "result3", "result4", "result5", "result6", "result7" "result8"], "max")

Sintaks aktual bervariasi dengan versi ArcMap Anda. (Misalnya, saya CellStatisticsbaru mengenal versi 10, tetapi operasi maksimum lokal selalu tersedia.)

Dalam contoh dengan input 01110101, urutan kisi "hasil *" akan berisi nilai 0, 1, 2, 3, 0, 1, 0, 1, jadi pada akhirnya CellStatisticsakan mengembalikan 3, panjang string terpanjang dari 1

Jika RAM langka, solusinya dapat dimodifikasi untuk menggunakan kembali hasil antara, dengan biaya sekitar dua kali lipat waktu eksekusi:

result = "grid1"
temp = con("grid2"==1, "result"+1, 0)
result = CellStatistics[["temp", "result"], "max"]
temp = con("grid3"==1, "temp"+1, 0)
result = CellStatistics[["temp", "result"], "max"]
...
temp = con("grid8"==1, "temp"+1, 0)
CellStatistics[["temp", "result"], "max"]

Dalam contoh dengan input 01110101, nilai ("temp", "result") akan menjadi (NoData, 0) setelah baris pertama dan setelah setiap pasangan ("Con", "CellStatistics") operasi nilainya adalah (1 , 1), (2, 2), (3, 3), (0, 3), (1, 3), (0, 3), (1, 3). Sekali lagi nilai akhir adalah 3.

Pola reguler ekspresi Peta Aljabar di salah satu solusi menunjukkan bagaimana kode algoritma sebagai loop dalam skrip, mengubah indeks yang sesuai dengan setiap iterasi.

Hanya mengobrol tentang hal ini dan bertanya-tanya apakah Anda dapat mendekati masalah dengan memperlakukan kisi masukan sebagai aliran biner. Ini akan memungkinkan Anda untuk menggabungkan mereka untuk memberikan integer ringkasan unik untuk urutan - yaitu 01110101 = 117. Nilai ini kemudian dapat direklasifikasi untuk memberikan jumlah maksimum 1s berturut-turut.

Berikut ini contoh yang menunjukkan satu cara untuk menggabungkan delapan kisi:

2*(2*(2*(2*(2*(2*(2*"g8" + "g7") + "g6") + "g5") + "g4") + "g3") + "g2") + "g1"

Operasi bitwise juga dapat ditekan ke dalam layanan untuk langkah ini. Atau, Anda dapat menggunakan kombinasi diikuti oleh perhitungan bidang. (Perhitungan lapangan akan memiliki ekspresi yang mirip dengan yang sebelumnya.)

The meja reklasifikasi harus menyediakan panjang run max untuk semua nilai antara 00000000B = 0 dan 11111111B = 255. Dalam rangka, di sini mereka adalah:

0, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 3, 4, 1, 1, 1, 2, 1, 1, 2, 3, 2, 2, 2, 2, 3, 3, 4, 5, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 3, 4, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 5, 6, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 3, 4, 1, 1, 1, 2, 1, 1, 2, 3, 2, 2, 2, 2, 3, 3, 4, 5, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 3, 4, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 6, 7, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 3, 4, 1, 1, 1, 2, 1, 1, 2, 3, 2, 2, 2, 2, 3, 3, 4, 5, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 2, 2, 3, 4, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 4, 4, 5, 6, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 3, 4, 2, 2, 2, 2, 2, 2, 2, 3, 2, 2, 2, 2, 3, 3, 4, 5, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 7, 8

Pendekatan ini terbatas pada sekitar 20 kisi di ArcGIS: menggunakan lebih dari ini dapat membuat tabel atribut yang sulit. ( Combinesecara khusus terbatas pada 20 kisi.)

Pernahkah Anda berpikir untuk mengubah nilai dari 0 dan 1 menjadi nilai dengan kekuatan 2 (1,2,4,8,16,32). Ketika Anda menggabungkan 8 kisi Anda akan mendapatkan nilai unik untuk setiap sel yang akan memberi Anda info berturut-turut (yaitu: nilai 3 berarti tahun 1 dan 2, di mana nilai 54 akan menjadi tahun 6 hingga 8).

Hanya pemikiran saja