Trovare il valore della cella raster più vicino in base al punto vettoriale?

Ho due livelli in ArcGIS: uno è un punto vettoriale e l'altro è un livello raster. Voglio aggiungere il valore della griglia di dati raster al punto vettoriale. Il problema è che il punto vettoriale è lontano da qualsiasi griglia raster esistente con valore, quindi il comando "Campione" non restituisce alcun valore per quel punto vettoriale.

Quindi la mia domanda è: come posso trovare la cella raster più vicina con un punto vettoriale ed estrarre il valore dalla cella?

In realtà ci sono due cose che devo fare:

1. Identifica il outliner all'interno di un ranger (di solito 1-2 celle di distanza dal mio raster esistente)
2. Assegna loro un valore dal raster più vicino, se si trovano in un intervallo

Modificato: ho circa 3000 punti vettoriali e il mio problema è questo:

Devo aggiungere il valore raster in base alle posizioni ai punti vettoriali. Ho usato lo strumento "Esempio" e funziona bene per la maggior parte dei punti.

Tuttavia, c'è una situazione come questa

Il punto giusto può ottenere il valore dal raster ("Esempio" funziona) ma il punto sinistro non può a causa di problemi di allineamento.

Dal raster al poligono non funziona perché ricevo un messaggio di errore che non riesco a trasformarli in vettore perché sono fuori dal dominio. Inoltre ho punti come questo

Quei punti superiori non dovrebbero ottenere alcun valore perché sono troppo lontani dal raster.

Ho pensato in qualche modo di aggirare:

Passaggio 1. Eseguire prima il "campione"

Passaggio 2. Selezionare il valore null dopo le operazioni "campione"

Passaggio 3. costruire zone cuscinetto (entro 1 decimale) in base ai punti null

Passaggio 4. utilizzare le statistiche zonali ?? ? o altri strumenti per eseguire il ciclo tra tutte le celle raster nel buffer, trovare quella più vicina, estrarne il valore e inserirlo nel punto vettoriale.

Sono bloccato al passaggio 4. Non so quale strumento di ArcSDK posso usare per ottenere questa funzione.

Oppure ... Hai un'idea migliore per affrontarlo?

Ho 20 livelli raster come questo e vorrei farlo automaticamente (usando modellista e arcobject).

Un profilo raster occupa un rettangolo delimitato in basso a sinistra dall'origine, con coordinate (diciamo) (x0, y0) e in alto a destra con (x1, y1): sono proprietà facilmente individuabili da qualsiasi raster. È possibile utilizzare queste informazioni per spostare eventuali punti al di fuori dell'estensione del raster sul punto più vicino sul confine del raster ed estrarre i valori del raster nelle nuove posizioni.

Compiti come questo sono spesso prontamente risolti considerando dimensioni inferiori . L'estensione di un raster unidimensionale è solo un intervallo di numeri, diciamo da x0 a x1> x0, scritto in modo convenzionale [x0, x1]. Supponiamo di voler trovare il punto più vicino in questo intervallo a un determinato numero x. Esistono tre casi: x <x0 (x è a sinistra), x0 <= x <= x1 (x è nell'intervallo) e x1 <x (x è a destra). Chiaramente i punti più vicini sono rispettivamente x0, xe x1. Una formula per questo risultato è

x -> min (max (x, x0), x1),

come puoi vedere contemplando ciascuno dei tre casi separatamente.

Poiché l'estensione di un raster è il prodotto cartesiano di due intervalli, [x0, x1] * [y0, y1], la stessa formula funziona in due dimensioni. Basta applicarlo a ciascun intervallo. Quindi le coordinate del punto proiettato possono essere calcolate (in un calcolatore di campo, per esempio) come

(x, y) -> (min (max (x, x0), x1), min (max (y, y0), y1).

Ecco questa formula in azione con 100 punti posizionati casualmente attorno e all'interno di una griglia:

Le linee tratteggiate associano visivamente i punti all'esterno dell'estensione alle posizioni in cui verranno spostati.

Crea un livello punti da queste coordinate calcolate e applica lo Samplestrumento per estrarre i valori raster. Unisci i risultati al livello punto originale.

In pratica , per evitare problemi con arrotondamenti a virgola mobile che potrebbero impedire ai nuovi punti di trovarsi all'interno della griglia, sarebbe saggio rendere x0 e y0 leggermente più grandi della vera origine (è possibile aggiungere metà della dimensione della cella a loro) e, allo stesso modo, rendere x1 e y1 leggermente più piccoli.

In sintesi , il flusso di lavoro consiste nel calcolare due campi (le nuove coordinate) nella tabella del livello punto originale, creare un livello evento punto da quelle nuove coordinate, eseguire lo Samplestrumento ed eseguire un join del database (non un join spaziale).

Quando hai molti più punti vettoriali rispetto alle celle nel raster , dovresti preferire la soluzione data da @celenius (che propone di convertire il raster in punti e usare un join spaziale). Normalmente, tuttavia, i raster hanno un sacco di celle - da milioni a miliardi - e la conversione in punti richiede così tanto tempo e dispendiosa da dover essere affrontata con cautela.

In alternativa, è possibile accelerare la soluzione di Celenius utilizzando un calcolo raster per creare una griglia intera con valori solo lungo il suo confine. La conversione che per un layer di punti è semplice e veloce perché avrà pochi punti. Chiamiamo questo "strato limite". Campionare la griglia nei punti del livello limite, copiando così i valori della griglia nella tabella degli attributi del livello limite. Dopo aver eseguito Samplecon il livello punto originale, rimuovere tutti i punti in cui è stato ottenuto correttamente un campione. Unisci spazialmente lo strato limite ai punti rimanenti per completare il processo di campionamento.

(Un modo per creare una griglia con valori solo lungo il suo confine è applicare un'operazione di vicinato focale che fallisce lungo il confine: esempi includono la ricerca di pendenze, pendenza delle colline, statistiche focali di quartiere 3 x 3. Questo crea una griglia con una cella singola anello spesso di NoData attorno al suo bordo. Usando IsNull e SetNull per rilevare le celle NoData e convertirle tra NoData e celle di dati si produce una griglia piena di valori NoData tranne attorno al limite.)

È possibile convertire il raster in punti ( ref ), quindi eseguire un join spaziale sui punti (dai raster) al punto.

(Ho appena imparato che i punti possono essere uniti spazialmente qui ieri )

È possibile convertire il livello raster in un poligono utilizzando lo strumento Raster in poligono ed eseguire l'unione spaziale (fare clic con il tasto destro del mouse sul livello punto e selezionare Join e relativi - Unisci, al primo menu a discesa selezionare l'opzione di unione spaziale e selezionare il più vicino per l'opzione sul pulsante readio) o utilizzare lo strumento Vicino.