Digi3D
Proyecto con Satélite RPC (Estereoscópico)
Digi3D 2007 presenta un tipos de sensor llamado Satélite RPC (Estereoscópico) que admite diferentes datos de entrada que engloban imágenes de diferentes satélites y datos de proyectos de aerotriangulación de Orima (LPS) en formato Blk.
Los satélites soportados son Ikonos/GeoEye, Quickbird/World View I, II, Alos, Aerotriangulaciones de Orima (LPS) en formato Blk.
Todos los satélites vienen equipados con un sistema de posicionamiento global. Basado en este sistema, se puede determinar la orientación completa de cada tira de imagen capturada.
Las siglas RPC provienen del inglés Rational Polynomial Coefficients. La orientación directa del sensor del satélite permite la determinación de la relación entre la imagen y las coordenadas terreno mediante un polinomio como función de las coordenadas geograficas X, Y, Z dividida por otras funciones. Haciendo uso de polinomios de tercer grado, la orientación puede ser expresada mediante 80 coeficientes. Esta puede ser incluso mejorada mediante puntos de control, en este caso el método se llama método corregido RPC.
En caso de querer abrir uno de los proyectos con los sensores mencionados, deberemos de seleccionar el tipo de sensor Satélite RPC (Estereoscópico) y veremos los siguientes parámetros:
Vamos a explicar el significado de los parámetros usados y la manera de cargarlos mediante un proyecto del satélite Alos (para saber más acerca de este satélite en particular lea el tema Satélite Alos).
Las imágenes de Alos se distribuyen como un conjunto de tiras de imagen separadas que al solaparlas forman una imagen completa.
Existe un archivo de volumen (cuyo nombre comienza por VOL-) que relaciona todas las tiras suministradas.
Cada una de estas tiras, dispone de su propio archivo con los polinomios RPC específicos para trabajar con la tira en cuestión.
Sin embargo, se distribuye un RPC adicional que es el que habría que utilizar si se concatenasen todas las tiras en un único archivo, por lo tanto el número de archivos con RPC suministrados es: 1+ número de tiras.
El satélite dispone de tras sensores: Backward, Nadiral y Frontal.
El cliente no tiene por qué adquirir imágenes de todos los sensores, con tener dos de estas tres, será suficiente para tener la visión estereoscópica.
A su vez, cada sensor dispone de 8 CCDs (cada uno de los cuales genera una tira), pero en las imágenes que distribuye JAXA (empresa propietaria del satélite) no tienen por qué venir todas las tiras, pues es posible que varias de ellas no tengan información (todos sus píxeles serían negros), y por lo tanto no se distribuyen para ahorrar datos.
Ejemplo: Tenemos una imagen que se llama ALPSMN178043865-O1B1___N y los únicos datos disponibles pertenecientes a esta imagen son las tiras generadas por los CCDs 3 al 6, por lo tanto los archivos RPC que tenemos son:
RPC-03-ALPSMN178043865-O1B1___N
RPC-04-ALPSMN178043865-O1B1___N
RPC-05-ALPSMN178043865-O1B1___N
RPC-06-ALPSMN178043865-O1B1___N
RPC-ALPSMN178043865-O1B1___N <----------------( rpc general)
Cada archivo RPC se corresponde con la imagen de cada tira, así el RPC 03 se corresponde con la imagen TIFF 03, sin embargo, el RPC general no está pensado para la concatenación de las fotos 03-04-05-06 sino 04-02-03-04-05-06, aunque no nos hayan distribuido las fotos 01 y 02.
La concatenación no es trivial, pues las imágenes no se concatenan sin más, ya que hay un solape de 32 píxeles entre tiras, de modo que el tamaño de la imagen total para esta imagen será: 6*(ancho de cada imagen) - 5*32
(Nótese cómo no se elimina la franja de 32 píxeles de la última tira).
No todas las imágenes tienen el mismo ancho, el sensor nadiral tiene un tamaño de 4992 píxeles y los backward y forward 4928 píxeles.
Práctica 1: Carga de un modelo formado por dos tiras
Vamos a cargar un par estereoscópico con la información suministrada por JAXA.
1. Examinaremos el directorio de ejemplo y comprobamos que disponemos de dos subdirectorios:
d0901485-001 y d0901485-002
Estos dos directorios contienen tanto los archivos TIFF correspondientes con cada una de las tiras como los RPCs de cada una de estas tiras así como el archivo RPC global.
En el directorio d0901485-001 tenemos la información de las cuatro tiras que forman la imagen nadiral: IMG-03-ALPSMN178043865-O1B1___N.tif,
IMG-04-ALPSMN178043865-O1B1___N.tif,
IMG-05-ALPSMN178043865-O1B1___N.tif ,
IMG-06-ALPSMN178043865-O1B1___N.tif
junto con sus correspondientes archivos RPC:
RPC-03-ALPSMN178043865-O1B1___N,
RPC-04-ALPSMN178043865-O1B1___N,
RPC-05-ALPSMN178043865-O1B1___N,
RPC-06-ALPSMN178043865-O1B1___N.
Y en el directorio d0901485-002 tenemos la información de la imagen backward:
IMG-03-ALPSMB178043920-O1B1___B.tif,
IMG-04-ALPSMB178043920-O1B1___B.tif,
IMG-05-ALPSMB178043920-O1B1___B.tif ,
IMG-06-ALPSMB178043920-O1B1___B.tif
junto con sus correspondientes archivos RPC:
RPC-03-ALPSMB178043920-O1B1___B,
RPC-04-ALPSMB178043920-O1B1___B,
RPC-05-ALPSMB178043920-O1B1___B,
RPC-06-ALPSMB178043920-O1B1___B.
2. Las imágenes suministradas por JAXA están rotadas 90 grados en el sentido de las agujas del reloj.
Como Digi3D 2007 no permite rotar las imágenes 90 grados, el usuario deberá rotar las imágenes 90º en sentido anti- horario en un programa de edición de imágenes como puedes ser Photoshop o GIMP.
Una vez girada la imagen IMG-04-ALPSMN178043865-O1B1___N.tif de la carpeta d0901485-001 la guardaremos con nombre 4.tif
Haremos lo mismo con la imagen IMG-04-ALPSMB178043920-O1B1___B.tif de la carpeta d0901485-002 y también la guardaremos con el mismo nombre 4.tif, ya que se trata de la tira 4.
3. Ya podemos entrar en Digi3D y seleccionar en el cuadro de diálogo Nuevo Proyecto (Ctrl+N) en la pestaña Estereoscópico el tipo de sensor Satélite RPC (Estereoscópico).
En el campo Directorio de proyecto seleccionamos el directorio donde se encuentra nuestro ejemplo.
El resto de los campos los rellenamos de la siguiente manera:
Pulsamos el botón Abrir.
Creación de una imagen completa:
Digi3D puede formar la imagen completa a partir de las tiras separadas.
La imagen creada puede estar formada únicamente por las tiras suminstradas o por todas las tiras que forman parte del sensor.
Para trabajar con el archivo RPC suminstrado por JAXA para la imagen completa, tendremos que unir las imágenes creando una imagen con todos los sensores CCD del sensor, incluidos los no suminstrados, ya que el archivo RPC suminstrado por JAXA está pensado de esta manera.
Práctica 2:
En esta práctica vamos a cargar el proyecto tal y como se distribuye por Jaxa, generaremos una imagen completa y cargaremos los RPCs de la imagen completa.
1. Examinaremos el directorio de ejemplo y comprobamos que disponemos de dos subdirectorios: d0901485-001 y d0901485-002
Estos dos directorios contienen tanto los archivos TIFF correspondientes a cada una de las tiras como los RPCs de cada una de estas tiras así como el archivo RPC global.
En el directorio d0901485-001 tenemos la información de la imagen nadiral:
ALPSMN178043865-O1B1___N
y en el directorio d0901485-002 tenemos la información de la imagen backward:
ALPSMB178043920-O1B1___B
Comprobamos que no existan los archivos ALPSMN178043865-O1B1___N.tif ni
ALPSMB178043920-O1B1___B.tif. Si existen los eliminamos.
2. Ejecutamos el programa Transformador de formato ALOS a Digi3D pulsando el botón Inicio/Todos los programas/Digi21/Digi3D 2007/Herramientas/Transformador de formato ALOS a Digi3D.
Este programa dispone en la parte superior de una ventana sobre la que podemos arrastrar directorios.
3. Abrimos un explorador de archivos, localizamos el directorio del ejemplo y seleccionamos los directorios d0901485-001 y d0901485-002, y los arrastramos sobre la ventana del programa Transformador de formato ALOS a Digi3D. Comprobaremos que el programa muestra los archivos de volumen localizados en esos directorios dentro de la ventana superior.
4. Seleccionamos la opción Generar la imagen completa.
5. Pulsamos el botón Aceptar
El programa se pondrá a procesar las imágenes y cuando finalice mostrará un cuadro de mensaje indicando que ha finalizado la operación satisfactoriamente.
Si se produce algún error, podremos ver en la ventana inferior del programa un informe que nos mostrará la causa del error.
6. Salimos del programa pulsando el botón Salir.
7. Comprobamos que en el directorio d0901485-001 se ha creado el archivo ALPSMN178043865-O1B1___N.tif y que en el directorio d0901485-002 se ha creado el archivo ALPSMB178043920-O1B1___B.tif.
Estos archivos tienen un peso de unos 600 MB aproximadamente ya que tienen muchos píxeles en negro (en la parte inferior, ya que las imágenes han sido rotadas 90 grados en sentido anti horario).
8. Cargamos Digi3D, abrimos el cuadro de diálogo de Nuevo Proyecto (Ctrl+N) en caso de que no se abra automáticamente y en la pestaña titulada Estereoscópico seleccionamos como directorio de proyecto el directorio de ejemplo.
En el campo Tipo de sensor seleccionamos Satélite RPC (Estereoscópico)
En el campo Imagen izquierda localizamos el archivo ALPSMB178043920-O1B1___B.tif del directorio d0901485-002.
En el campo RPCs cámara izquierda localizamos el archivo rpc-ALPSMB178043920-O1B1___B del
directorio d0901485-002.
En el campo Imagen derecha localizamos el archivo ALPSMN178043865-O1B1___N.tif del directorio d0901485-001.
En el campo RPCs cámara izquierda localizamos el archivo RPC-ALPSMN178043865-O1B1___N del directorio d0901485-001.
En el campo Proyección Cartográfica seleccionamos UTM.
En el campo Hemisferio seleccionamos Sur.
En el campo Huso seleccionamos 18.
Pulsamos el botón Abrir.
Comprobaremos cómo Digi3D carga el par estereoscópico y nos permite desplazarnos por el modelo estereoscópicamente.
Aerotriangulación
Es posible que los modelos estén aerotriangulados. Si éstos están aerotriangulados el proceso de carga es muy sencillo, únicamente habrá que seleccionar cómo archivo de RPCs el archivo con el resultado de la aerotriangulación.
Práctica 3: Carga de modelos con orientaciones obtenidas de un programa de cálculo de aerotriangulación
En esta práctica, vamos a cargar el mismo modelo de la práctica 2, pero con las orientaciones generadas mediante una aerotriangulación que se ha realizado con el software Orima de LPS.
1. El software de LPS que genera la imagen completa (concatenación de todas las tiras) no influye en la imagen resultante las tiras de los CCDs no suministrados por JAXA. De modo que debemos crear la imagen completa incluyendo únicamente la información de los sensores suministrados.
Examinaremos el directorio de ejemplo y comprobamos que disponemos de dos subdirectorios: d0901485-001 y d0901485-002.
Estos dos directorios contienen tanto los archivos TIFF correspondientes con cada una de las tiras como los RPCs de cada una de estas tiras así como el archivo RPC global.
En el directorio d0901485-001 tenemos la información de la imagen nadiral: ALPSMN178043865-O1B1___N
y en el directorio d0901485-002 tenemos la información de la imagen backward: ALPSMB178043920-O1B1___B
Comprobamos que no existan los archivos ALPSMN178043865-O1B1___N.tif ni ALPSMB178043920-O1B1___B.tif. Si existen los eliminamos.
2. Ejecutamos el programa Transformador de formato ALOS a Digi3D pulsando el botón Inicio/Todos los programas/Digi21.net/Digi3D 2007/Herramientas/Transformador de formato ALOS a Digi3D.
Este programa dispone en la parte superior de una ventana sobre la que podemos arrastrar directorios.
3. Abrimos un explorador de archivos, localizamos el directorio del ejemplo y seleccionamos los directorios d0901485-001 y d0901485-002, y los arrastramos sobre la ventana del programa Transformador de formato ALOS a Digi3D. Comprobaremos que el programa muestra los archivos de volumen localizados en esos directorios dentro de la ventana superior.
4. Nos aseguramos de que la opción Generar la imagen completa no está seleccionada.
5. Pulsamos el botón Aceptar.
El programa se podrá a procesar las imágenes y cuando finalice mostrará un cuadro de mensaje indicando que ha finalizado la operación satisfactoriamente.
Si se produce algún error, podemos ver en la ventana inferior del programa un informe que nos mostrará la causa del error.
6. Salimos del programa pulsando el botón Salir.
7. Comprobamos que en el directorio d0901485-001 se ha creado el archivo ALPSMN178043865-O1B1___N.tif y que en el directorio d0901485-002 se ha creado el archivo ALPSMB178043920-O1B1___B.tif.
Estos archivos tienen un peso de unos 448 MB aproximadamente.
8. Cargamos Digi3D, abrimos el cuadro de diálogo de Nuevo Proyecto (Ctrl+N) en caso de que no se abra automáticamente y en la pestaña titulada Estereoscópico seleccionamos como directorio de proyecto el directorio de ejemplo.
En el campo Tipo de sensor seleccionamos Satélite RPC (Estereoscópico).
En el campo Imagen izquierda localizamos el archivo ALPSMB178043920-O1B1___B.tif del directorio d0901485-002.
En el campo RPCs la cámara izquierda localizamos el archivo ica_huanca.blk del directorio de ejemplo.
En el campo Imagen derecha localizamos el archivo ALPSMN178043865-O1B1___N.tif del directorio d0901485-001.
En el campo RPCs cámara derecha localizamos el archivo ica_huanca.blk del directorio de ejemplo.
En el campo Proyección cartográfica seleccionamos UTM.
En el campo Hemisferio seleccionamos Sur.
En el campo Huso seleccionamos 18.
Pulsamos el botón Abrir.
Comprobaremos cómo Digi3D carga el par estereoscópico y nos permite desplazarnos por el modelo estereoscópicamente.
Fecha de última modificación: 20 de Diciembre, 2010