Calcolare un nuovo attributo basato sulle modifiche in un altro attributo utilizzando ArcGIS Desktop con Python?

Sto tentando di classificare un insieme di dati punto codificati nel tempo gps in comportamenti basati su attributi diversi.

Ho creato un attributo 0 per la casa e 1 per la trasferta in base alla posizione e ora voglio numerare i viaggi lontano da casa (una serie di punti 01111111111110sarebbe un viaggio perché è iniziato e terminato a casa). Ho aggiunto il campo attributo che avrà i numeri di viaggio, ma non so come calcolare il campo in modo che sia basato sul campo casa / fuori.

Ecco un esempio dei dati GPS (utilizzando "*" per indicare informazioni irrilevanti e semplicemente indicizzare i tempi come 1, 2, ecc.), L'indicatore "Casa / Ospite" sopra descritto e l'indicatore di viaggio desiderato, "Viaggio", che devo calcolare:

Time Lat Lon Home/Away Trip
   1   *   *         0    0
   2   *   *         1    1
   3   *   *         1    1
....
  12   *   *         1    1
  13   *   *         0    0
  14   *   *         0    0
  15   *   *         1    2
  16   *   *         1    2
.... 
  34   *   *         1    2
  35   *   *         0    0
  36   *   *         0    0
  37   *   *         1    3
....

Il mio set di dati è troppo grande per passare manualmente e numerare ogni viaggio nella tabella degli attributi, quindi esiste un modo per calcolare il campo in base a come viene ordinato l'attributo home / away e ogni "gruppo" di punti away è designato come viaggio?

Queste sono le ossa nude di come potrebbe apparire il codice Python (non ho esperienza con il codice).

Espressione:

trip = Reclass(!home!)

codeblock:

def Reclass(home):  
  if (home = 0):  
    return 0   
  elif (home = 1 and lastValue = 0):  
    return _(incremental numbering?)_  
  elif (home = 1 and lastValue = 1):  
    return lastValue  

Dopo aver utilizzato lo script consigliato da Matt Wilson ho apportato alcune modifiche in modo che il mio primo viaggio sia il numero 1, il secondo sia il 2, ecc. Ecc.

Ecco il codice modificato da Matt:

import arcpy
rows = arcpy.UpdateCursor("test2")

trip = 0
for row in rows:
    if row.home == 0:
        prev = row.home
        row.TRIP = trip
        rows.updateRow(row)

    elif row.home == 1 and prev == 0:
        trip += 1
        prev = row.home
        row.TRIP = trip
        rows.updateRow(row)
        rows.next()

    elif row.home == 1 and prev == 1:
        prev = row.home
        row.TRIP = trip
        rows.updateRow(row)
        rows.next()

    row.TRIP = trip
    rows.updateRow(row)


del row, rows

Quindi seleziono solo home = 0 e calcolo il campo del mio viaggio a 0. Viaggi ordinatamente ordinati.

Per questo è possibile utilizzare UpdateCursor , che apre la classe di funzionalità o la tabella e passa attraverso ciascun record (riga) in modo incrementale.

Lo script seguente funziona su questi dati di test

+-----------------------+
| Time| Home_Away|Trip  |
+-----|----------|------+
|  1  |  0       | <nul>|
|  2  |  1       | <nul>|
|  4  |  1       | <nul>|
|  5  |  0       | <nul>|
|  6  |  0       | <nul>|
|  7  |  1       | <nul>|
|  9  |  1       | <nul>|
| 12  |  1       | <nul>|
| 13  |  0       | <nul>|
+-----------------------+

.

import arcpy
fc = r'D:\s\py\pyscratch.gdb\gps_points'

# open the feature class and create the cursor
rows = arcpy.UpdateCursor(fc)

trip = 0
for row in rows:
    if row.HOME_AWAY == 0:
        trip += 1           # start of new trip, increment counter
        row.TRIP = trip     # calc the TRIP field to be current trip#
        rows.updateRow(row) # save
        print "Trip %s started at %s" % (trip, row.TIME)

    # keep cycling through records until HOME_AWAY is not 1
    while row.HOME_AWAY == 1:
        row.TRIP = trip
        rows.updateRow(row)
        rows.next() # move to next record

    # this is for the trailing end of a trip, the second 0
    # print "     %s ended at %s" % (trip, row.TIME)
    row.TRIP = trip
    rows.updateRow(row)

# remove programming objects and data locks
# the data itself is left alone
del row, rows

La fine finale del blocco di viaggio viene effettivamente eseguita anche per l'inizio di un viaggio, ma poiché il contatore di viaggi è corretto, il doppio calcolo sulla riga di inizio-viaggio non ha importanza. Scomponi l'istruzione print in quel blocco per vedere cosa intendo.

Python aggiunge automaticamente un implicito rows.next()alla fine per il for row in rowsblocco.

Ciò presuppone l'integrità dei dati. Si incasinerà se ci sarà mai un numero dispari di zero record Casa / Ospite in una riga ( 000o 00000). Un viaggio che consiste solo in start e stop dovrebbe andare bene, ad esempio una sequenza di 3 trip di 01..10 00 01..10, in cui gli spazi indicano gli spazi tra i viaggi. In altre parole, convalida i risultati!

La guida di ArcGIS 10 in "calcola esempi di campi" mostra come "Calcolare il valore cumulativo di un campo numerico". Questo funzionerà, a condizione che i dati siano fisicamente nell'ordine temporale previsto.

Per applicarlo direttamente, inverti l'indicatore [Home / Away] (sottrai da 1) in modo che "0" significhi "away" e "1" significhi "home". Lo chiamo [Fuori casa] nell'esempio di seguito.

Calcola il suo valore cumulativo - [Cumulativo] nell'esempio.

Aggiungi uno e dividi per due - [Trip] nell'esempio (quasi).

Infine, imposta [Trip] su zero per tutti i record "home". Ora i risultati concordano con l'esempio:

Time Lat Lon Home/Away Trip Away/Home Cumulative 
   1   *   *         0    0         1          1
   2   *   *         1    1         0          1
   3   *   *         1    1         0          1
.... 
  12   *   *         1    1         0          1
  13   *   *         0    0         1          2
  14   *   *         0    0         1          3
  15   *   *         1    2         0          3
  16   *   *         1    2         0          3
.... 
  34   *   *         1    2         0          3
  35   *   *         0    0         1          4
  36   *   *         0    0         1          5
  37   *   *         1    3         0          5
....

Per la cronaca, ecco il codice preso dalla guida di ArcGIS 10. L'ho modificato leggermente in modo da fare tutti i passaggi contemporaneamente: ora devi solo eseguirlo. Dovrebbe essere chiaro dove [Casa / Ospite] viene invertito e dove si verifica il passaggio "aggiungi 1, dividi per 2".

Espressione:

acc(!Home/Away!)

Tipo di espressione:

PYTHON_9.3

Blocco di codice:

t=0
def acc(i):
  global t
  if t:
    t += (1-i)
  else:
    t = 1
  if i:
    return (t+1)/2
  else:
    return 0