Estoy trabajando en un estudio de detección de cambios para observar los efectos de los herbicidas en el crecimiento de la vegetación. Por lo tanto, tengo imágenes de Worldview 2 (4 bandas + banda de pan) para dos años posteriores que me gustaría utilizar para un análisis de detección de cambios, pero necesito algunos consejos sobre los pasos adecuados de preprocesamiento.
Las técnicas que propongo para evaluar el cambio son, en primer lugar, una simple comparación del NDVI entre los dos años y, en segundo lugar, ejecutar una clasificación forestal aleatoria en R. He tratado de seguir las orientaciones indicadas en el Digital Globe Documento sobre el uso radiométrico de las imágenes Worldview 2 pero estoy teniendo un poco de problemas para descifrar exactamente qué fórmulas tengo que utilizar para que mis imágenes sean comparables entre años.
Inicialmente, mi flujo de trabajo era:
1) Convierta el DN en radiancia espectral de la parte superior de la atmósfera utilizando la segunda fórmula de la sección 5.1 del documento (véase la ecuación 1 a continuación) utilizando la calculadora de ráster de ArcGIS.
2) Convertir los resultados del paso 1 en reflectancia media de banda de la parte superior de la atmósfera utilizando la última ecuación de la sección 7 del documento técnico (véase la ecuación 2 a continuación) utilizando la calculadora de ráster de ArcGIS.
3) Ejecutar la clasificación de NDVI y R forest utilizando las capas de reflectancia de banda media de la parte superior de la atmósfera adecuadas.
Sin embargo, después de volver a revisar el documento técnico de Digital Globe, me pregunto si también tengo que completar el equilibrio radiométrico para múltiples escenas que se indica en la sección 6 (véanse las ecuaciones 3 o 4 más abajo) o si este paso ya se tiene en cuenta en la fórmula de conversión a reflectancia media de la banda superior de la atmósfera (ecuación 2).
Ecuaciones del documento técnico propuesto: