El Proyecto REGENTE

Con el fin de establecer una cartografía europea unificada, se hace indispensable la conversión de las coordenadas de los Marcos de los Sistemas Geodésicos Nacionales al Marco ETRF89, lo cual sólo es posible a través de la determinación de transformaciones y superficies de ajuste desde uno a otro marco. Tal determinación requiere el conocimiento de ambas clases de coordenadas en un número muy elevado de estaciones uniformemente distribuidas, debiendo este número ser tanto mayor cuantas más irregularidades presente el Marco local. En el caso de la Península y Archipiélagos, el IGN decidió resolver el problema mediante el Proyecto REGENTE (Red Geodésica Nacional por Técnicas Espaciales), consistente en el establecimiento de una densa red GPS de alta precisión con estaciones coincidentes con vértices de ROI y con clavos de las líneas NAP. La densidad media quedó fijada en una estación por Hoja del MTN 1:50.000, es decir, de una estación por cada 300 km2. REGENTE quedará perfectamente enlazada con la red de referencia europea ETRF89 por medio de las redes ibéricas IBERIA95 y BALEAR98. REGENTE Canarias se apoya, como estación de referencia, en la estación GPS de Maspalomas, incluida en la red ITRF93.

Red Nacional de Estaciones de Referencia GPS (ERGPS)

El Instituto Geográfico Nacional, por medio del Área de Geodesia, está llevando a cabo el establecimiento de una Red de Estaciones Permanentes GPS que permitan obtener coordenadas muy precisas, así como sus campos de velocidades en un Sistema de Referencia Global (ITRFxx). Dichas estaciones pertenecen a la Red de Estaciones Permanentes de EUREF (EUropean REference Frame) y constituyen el orden cero de la Geodesia Española.

Efectos dinámicos laterales en vehículos y puentes ferroviarios sometidos a la acción de vientos transversales

En las últimas décadas el aumento de la velocidad y la disminución del peso de los vehículos ferroviarios de alta velocidad ha provocado que aumente su riesgo de vuelco. Además, las exigencias de los trazados de las líneas exige en ocasiones la construcción de viaductos muy altos situados en zonas expuestas a fuertes vientos. Esta combinación puede poner en peligro la seguridad de la circulación. En esta tesis doctoral se estudian los efectos dinámicos que aparecen en los vehículos ferroviarios cuando circulan sobre viaductos en presencia de vientos transversales. Para ello se han desarrollado e implementado una serie de modelos numéricos que permiten estudiar estos efectos de una forma realista y general. Los modelos desarrollados permiten analizar la interacción dinámica tridimensional tren-estructura, formulada mediante coordenadas absolutas en un sistema de referencia inercial, en un contexto de elementos _nitos no lineales. Mediante estos modelos se pueden estudiar de forma realista casos extremos como el vuelco o descarrilamiento de los vehículos. Han sido implementados en Abaqus, utilizando sus capacidades para resolver sistemas multi-cuerpo para el vehículo y elementos finitos para la estructura. La interacción entre el vehículo y la estructura se establece a través del contacto entre rueda y carril. Para ello, se han desarrollado una restricción, que permite establecer la relación cinemática entre el eje ferroviario y la vía, teniendo en cuenta los posibles defectos geométricos de la vía; y un modelo de contacto rueda-carril para establecer la interacción entre el vehículo y la estructura. Las principales características del modelo de contacto son: considera la geometría real de ambos cuerpos de forma tridimensional; permite resolver situaciones en las que el contacto entre rueda y carril se da en más de una zona a la vez; y permite utilizar distintas formulaciones para el cálculo de la tensión tangencial entre ambos cuerpos. Además, se ha desarrollado una metodología para determinar, a partir de formulaciones estocásticas, las historias temporales de cargas aerodinámicas debidas al viento turbulento en estructuras grandes y con pilas altas y flexibles. Esta metodología tiene cuenta la variabilidad espacial de la velocidad de viento, considerando la correlación entre los distintos puntos; considera las componentes de la velocidad del viento en tres dimensiones; y permite el cálculo de la velocidad de viento incidente sobre los vehículos que atraviesan la estructura. La metodología desarrollada en este trabajo ha sido implementada, validada y se ha aplicado a un caso concreto en el que se ha estudiado la respuesta de un tren de alta velocidad, similar al Siemens Velaro, circulando sobre el viaducto del río Ulla en presencia viento cruzado. En este estudio se ha analizado la seguridad y el confort de la circulación y la respuesta dinámica de la estructura cuando el tren cruza el viaducto. During the last decades the increase of the speed and the reduction of the weight of high-speed railway vehicles has led to a rise of the overturn risk. In addition, the design requests of the railway lines require some times the construction of very tall viaducts in strong wind areas. This combination may endanger the traffic safety. In this doctoral thesis the dynamic effects that appear in the railway vehicles when crossing viaducts under strong winds are studied. For this purpose it has been developed and implemented numerical models for studying these effects in a realistic and general way. The developed models allow to analyze the train-structure three-dimensional dynamic interaction, that is formulated by using absolute coordinates in an inertial reference frame within a non-linear finite element framework. By means of these models it is possible to study in a realistic way extreme situations such vehicle overturn or derailment. They have been implemented for Abaqus, by using its capabilities for solving multi-body systems for the vehicle and finite elements for the structure. The interaction between the vehicle and the structure is established through the wheel-rail contact. For this purpose, a constraint has been developed. It allows to establish the kinematic relationship between the railway wheelset and the track, taking into account the track irregularities. In addition, a wheel-rail contact model for establishing the interaction of the vehicle and the structure has been developed. The main features of the contact model are: it considers the real geometry During the last decades the increase of the speed and the reduction of the weight of high-peed railway vehicles has led to a rise of the overturn risk. In addition, the design requests of the railway lines require some times the construction of very tall viaducts in strong wind areas. This combination may endanger the traffic safety. In this doctoral thesis the dynamic effects that appear in the railway vehicles when crossing viaducts under strong winds are studied. For this purpose it has been developed and implemented numerical models for studying these effects in a realistic and general way. The developed models allow to analyze the train-structure three-dimensional dynamic interaction, that is formulated by using absolute coordinates in an inertial reference frame within a non-linear finite element framework. By means of these models it is possible to study in a realistic way extreme situations such vehicle overturn or derailment. They have been implemented for Abaqus, by using its capabilities for solving multi-body systems for the vehicle and finite elements for the structure. The interaction between the vehicle and the structure is established through the wheel-rail contact. For this purpose, a constraint has been developed. It allows to establish the kinematic relationship between the railway wheelset and the track, taking into account the track irregularities. In addition, a wheel-rail contact model for establishing the interaction of the vehicle and the structure has been developed. The main features of the contact model are: it considers the real geometry

Comparación de algoritmos detectores de puntos singulares para reconocimiento de objetos en vídeo quirúrgico.

El análisis de vídeo laparoscópico ofrece nuevas posibilidades a la navegación quirúrgica al garantizar una incorporación mínima de tecnología en quirófano, evitando así alterar la ergonomía y los flujos de trabajo de las intervenciones. Una de sus principales ventajas es que puede servir como fuente de datos para reconstruir tridimensionalmente la escena laparoscópica, lo que permite dotar al cirujano de la sensación de profundidad perdida en este tipo de cirugía. En el presente trabajo de investigación se comparan dos detectores de puntos singulares, SIFT y SURF, para estimar cuál de los dos podría integrarse en un algoritmo de cálculo de coordenadas 3D, MonoSLAM, basado en la detección y el seguimiento de estos puntos singulares en los fotogramas del vídeo. Los resultados obtenidos posicionan a SURF como la mejor opción gracias a su rapidez y a su mayor capacidad de discriminación entre estructuras anatómicas e instrumental quirúrgico.

Estudio y aplicación de redes neuronales artificiales al análisis y diseño de antenas reflectarrays

En los últimos años, estamos siendo testigos de la alta implantación en la sociedad de dispositivos de comunicación. Lo que hace años estaba reservado a un público reducido, con claras necesidades en comunicación, se ha trasladado al público general, dado la amplia variedad de servicios que sobre los nuevos medios de comunicación se han desarrollado. De hecho, el mayor tráfico de datos en la actualidad no se produce al hilo de necesidades de máxima importancia, sino como producto de nuevos hábitos cotidianos. En este contexto de renovación tecnológica constante en busca de la eficiencia, las antenas reflectoras reflectarray (o, simplemente, los reflectarrays, RAs, [1]) se presentan como una opción competitiva contra los reflectores parabólicos metálicos. En su versión más simple, una antena reflectarray se trata de una estructura compuesta de un elemento alimentador radiante, como puede ser una bocina, y de una superficie plana, consistente en multitud de elementos individuales dispuestos en una rejilla periódica. Sobre esta superficie plana, los frentes de onda provenientes del alimentador son reflejados formando frentes de ondas planas, de una manera análoga a como lo hace un reflector parabólico. A partir de la configuración inicial, y centrándose en el principio de funcionamiento, se ha ido modificando el tipo de elemento RA empleado, consiguiendo RA cada vez más operativos. Es, sobre todo, con el nacimiento de la tecnología impresa cuando las antenas RAs vuelven a cobrar interés. Aunque el uso de tecnología impresa supuso un gran impulso en los RAs, también abrió otros desafíos en lo que al diseño de ellos se refiere. Desde el punto de vista del análisis, es común suponer que el elemento RA se encuentra en un ambiente infinitamente periódico, de forma que se puedan aplicar las condiciones de contorno de Floquet (suposición de periodicidad local). Desde un punto de vista funcional, en general, los elementos RA de tecnología impresa presentan un ancho de banda reducido, que condiciona el ancho de banda del RA completo. Entre las soluciones aportadas, es comúnmente aceptado que las estructuras multicapa, con resonadores a distintas frecuencias cercanas, pueden mitigar en parte el problema del ancho de banda. Por ello, en la actualidad, los elementos RA más comunes están compuestos por varios elementos resonadores, cuyas dimensiones constituyen los parámetros de diseño libres. Es decir, en función de dichas dimensiones, el elemento RA tendrá un valor del coeficiente de reflexión u otro. Esto supone un aumento en la complejidad a la hora de analizar dicho elemento por los métodos numéricos conocidos, como el Método de los Momentos (MoM) o el Método de Elementos Finitos (FEM, por las siglas de su traducción inglesa Finite Element Method), que redundará en un mayor tiempo de cómputo en el análisis. Por otra parte, como se muestra en la Figura R.1, el diseño de un RA conlleva analizar multitud de veces el elemento RA considerado. En efecto, se trata de un método de diseño indirecto, en donde las dimensiones de los parámetros geométricos libres de cada elemento RA se obtienen de manera iterativa usando un optimizador. Se ve claro, entonces, que el aumento en tiempo de análisis del elemento RA repercute en gran medida en el tiempo de diseño total, por lo que una reducción en el tiempo de análisis del elemento RA podría ser muy beneficioso. Uno de los métodos para conseguir reducir el tiempo de diseño de un RA, que se puede encontrar en la literatura, es emplear un modelo de la respuesta del elemento RA en función de los parámetros libres. La cuestión que aflora es cuál es la herramienta idónea para modelar la respuesta del elemento RA. En los últimos años se han propuestos varias formas. La primera de ellas consistía en encontrar un equivalente circuital. Esta aproximación está bien extendida para otras estructuras EM, donde los equivalentes circuitales con componentes LC ofrecen respuestas muy precisas con respecto a las que ofrecen las estructuras EM en sí. A raíz del carácter no lineal de la respuesta, hay autores que han propuesto para el diseño de RAs la creación de tablas de datos (look up tables) que, para cada parámetro de diseño de interés (suele ser el desfase introducido por el elemento) guardan las dimensiones de los parámetros geométricos libres asociados. De esta forma, consiguen un diseño rápido, pero poco versátil, ya que la tabla ofrece un único valor para cada entrada, por lo que es difícil jugar con más de una restricción de diseño. Más recientemente, se está comenzando a utilizar, para la caracterización de estructuras EM, unos sistemas llamados Redes Neuronales Artificiales (ANN, por sus siglas en inglés Artificial Neural Network). El uso fundamental de los mismos en EM es el de servir como interpoladores no lineales. Se trata de sistemas que admiten múltiples parámetros de entradas y múltiples parámetros de salida. Antes de poder ser usados como interpoladores, deben ser entrenados. Para ello, necesitan de un conjunto de pares de los parámetros de entrada a la red, con los valores de las salidas asociados. Algunos usos en electromagnetismo de las ANNs que se pueden encontrar en la literatura son: el modelado de filtros; la caracterización de dispositivos activos; la obtención de modelos que aceleran los algoritmos que calculan la dirección de llegada en antenas de radar; o el diseño de arrays de antenas. Volviendo al modelado de elementos RA, en este trabajo haremos uso de las ANNs para caracterizar distintos tipos de elementos RA. A lo largo de estos últimos años, se ha considerado esta posibilidad como una de las más prometedoras. De hecho, podemos encontrar algunas pocas referencias al respecto, varias de las cuales han sido publicadas por distintos autores durante la elaboración del trabajo recogido en esta Tesis. Como veremos, los resultados que vamos a presentar aportan novedades con respecto a la citada literatura. Particularmente, en este trabajo se ha realizado la caracterización de un elemento RA de tres capas, considerando hasta 9 parámetros de entrada (seis parámetros geométricos, las dos coordenadas del ángulo de incidencia, y la frecuencia) y 4 parámetros de salida complejos (los coeficientes de reflexión para dos polarizaciones ortogonales lineales). Haciendo uso de esta caracterización en el flujo de diseño de RAs, se ha realizado el análisis y el diseño de varias antenas RA con restricciones de diseño de comunicaciones espaciales. Los resultados fueron exitosos comparados con los resultados obtenidos por los métodos tradicionales. De manera puntualizada, podríamos resumir las aportaciones que se verán en esta Tesis como: Caracterización de distintos elementos RA mediante ANNs basadas en el Perceptrón Multicapa (MLP). En concreto, se ha realizado con éxito la caracterización de un elemento RA de parche acoplado a línea de retardo a través de apertura; la caracterización de un elemento RA basado en dipolos sobre substratos de distintas características eléctricas en el rango de centenas de GHz; y la caracterización de un elemento RA basado en 3 parches apilados, con 9 parámetros libres en su caracterización. Uso del FEM, de la técnica de segmentación en subdominios y de la generación y manipulación de accesos MAM para el análisis y la caracterización de elementos RA mediante ANNs. Desarrollo de una nueva técnica de obtención de muestras, para el caso de estructura multicapa cuyo estudio EM se pueda dividir en dos pasos: estudio de cada capa y conexión de capas. De esta forma, se ha podido reducir en varios órdenes de magnitud el tiempo necesario para obtener el set de entrenamiento de las ANNs. Valoración del uso de distintos métodos de entrenamiento de segundo orden para el entrenamiento de redes ANN MLP, en la caracterización de elementos RA. Desarrollo de una nueva técnica para realizar el entrenamiento de redes ANNs basadas en el MLP, denominada como Entrenamiento en Cascada. Dado el alto número de parámetros a caracterizar, era difícil conseguir una red que, partiendo del número de entradas deseado, proporcionara convergencia con precisión suficiente. Con el algoritmo propuesto y probado en esta Tesis, se consiguió entrenar redes de 8 parámetros de entradas (el noveno parámetro, la frecuencia, correspondía a redes diferentes para cada valor) con gran precisión. Desarrollo de un método adaptativo para mejorar la precisión de las ANNs en el análisis de antenas RA. Este método, basado en re-entrenar las ANNs para sub rangos de los parámetros de entrada en donde el error es mayor, aporta una precisión mayor, al mejorar el entrenamiento global de las ANNs, en un tiempo aceptable, ya que solo se incluyen nuevas muestras en torno a los valores donde el error es mayor. Análisis de antena RA completa, con cobertura según especificaciones de la misión AMAZONAS (haz conformado, banda Ku), usando las caracterización el elemento RA obtenida mediante ANNs. La mejora en tiempo de análisis conseguida con respecto al uso del MoM está en un factor 102, con precisiones comparables. Diseño de antenas RA completas, con especificaciones de haz pincel y EuTELSAT (banda Ku). De nuevo, la mejora en tiempo de diseño conseguida están en torno a 102. De todos los puntos anteriores, son de destacar los dos últimos, que forman el objetivo principal de esta Tesis. Esto es, el uso de modelos rápidos de elementos RA mediante ANNs para el análisis y el diseño de antenas para comunicaciones por satélite.

Desarrollo de una herramienta software para el manejo de un teléfono móvil adaptada a personas con discapacidad física severa

El presente proyecto desarrolla una aplicación residente en un terminal móvil, que pretende proporcionar un valor añadido al actual proyecto Localiza, sistema de localización bajo demanda para personas con discapacidad severa. Mediante el desarrollo de este proyecto se pretende facilitar el acceso al teléfono móvil y al ordenador a las personas con discapacidad motriz. El objetivo final es ser capaz de controlar un teléfono móvil por medio de control remoto, mediante el uso de un ordenador personal. Para ello se establece una conexión remota entre el terminal móvil y el ordenador personal, a través del protocolo de comunicación Bluetooth. De este modo, a través de la aplicación móvil se transmite la información de posición de las coordenadas, proporcionada por el acelerómetro del terminal, a un servicio instalado en el ordenador que se encarga de gestionar la información recibida, y así crear las interrupciones pertinentes en el sistema operativo para mover el puntero del ratón. Para controlar el teléfono móvil de forma remota se dispondrá de un emulador de telefonía móvil instalado en el ordenador que implemente las funciones básicas de control de llamadas. Por medio de comunicación Bluetooth, las acciones que realice el usuario en emulador serán invocadas en el propio terminal móvil. SUMMARY. The project presented develops a mobile application, which is intended to provide an added value to the already existing project Localiza, on-demand position system for people with severe disabilities. This project aims to facilitate the access to the personal computer and to the mobile telephony environment for disabled people. The main goal is to be able to control a mobile phone by remote control, using a personal computer. Thus, a remote connexion will to be established between the mobile device and the personal computer, through Bluetooth communication protocol. Thus, the mobile application will transmit the coordinate’s position, provided by the accelerometer of the mobile device, to a Bluetooth service running in the personal computer. That service will be in charge of managing the information received in order to create the interruptions on the operational system for moving the mouse pointer. The remote controlling of the mobile device is carried out using a mobile telephony emulator installed in the personal computer, which will implement the basic functionality of calling control. Using Bluetooth communication, the user actions done in the emulator interface will be invoked on the mobile device itself.

Propuesta de un núcleo estándar de metadatos para los recursos del Patrimonio Histórico Español

La amplitud con la que se define legalmente el Patrimonio Histórico Español según la Ley 16/1985 de 25 de junio es uno de los factores claves en la identificación y caracterización del mismo: “... inmuebles y objetos muebles de interés artístico, históricos, paleontológicos, arqueológicos, etnográficos, científico o técnico. Asimismo el patrimonio documental o bibliográfico, los yacimientos o zonas arqueológicas así como sitios naturales, jardines y parques que tengan valor artístico, histórico y antropológico”. Los recursos patrimoniales están catalogados y descritos mediante un conjunto de metadatos, que permiten su búsqueda en la red, así como la consulta de las propiedades básicas de los mismos. INSPIRE plantea la directiva donde quedan enmarcados los recursos digitales georreferenciados. En este sentido el patrimonio es un recurso con componente espacial y la propia directiva lo considera dentro de los datos de referencia (Sistemas de coordenadas de referencia, Sistema de cuadrículas geográficas, Nombres geográficos, Unidades administrativas, Direcciones, Parcelas catastrales, Redes de transporte, Hidrografía, Lugares protegidos) como se puede ver en el anexo I. Según INSPIRE: “Lugares protegidos: zonas designadas o gestionadas dentro de un marco legislativo internacional, comunitario o propio de los estados miembros, para la consecución de unos objetivos de conservación específicos”.

Estudio de la evolución del color de la piel y de la pulpa de aguacate cv. "Hass" y su aplicación para la evaluación de la madurez post-recolección.

Los frutos sufren tras la recolección diversos cambios fisicoquímicos que determinan su calidad intrínseca. Uno de los cambios más importantes y manifiestos que ocurren durante la maduración, para muchos frutos es el del color. El resaltado mis común de las modificaciones de color es la pérdida de la pigmentación verde como consecuencia de la degradación de la clorofila, que va asociada a la síntesis o al desenmascaramiento de otros pigmentos. Los cambios de color o desaparición del color verde de la piel constituye un buen indicativo no destructivo del grado de madurez para muchos frutos, incluyendo el aguacate "Hass". Actualmente se viene realizando una extensa investigación direccionada a la detección no destructiva del grado de madurez de los frutos, desarrollando sistemas y experimentando sistemas y técnicas basadas en propiedades físicas de los productos. Sin embargo, para los frutos que presentan el cambio de color de piel durante la maduración, el análisis de estos cambios presenta diversas ventajas sobre otras técnicas, como la manipulación mínima de los frutos y la no necesidad de acoplamiento de sensores en la superficie del fruto. Además del potencial de utilización en procesos reales con alta productividad. Así varios investigadores han estudiado el color de los frutos con el fin de evaluar el grado de madurez de los mismos. Así Rood (1957), citado por Delwiche (1987), concluye que los mejores índices de madurez para el melocotón son por el orden de importancia la firmeza de la pulpa, color de la piel, color de la pulpa y contenido de clorofila de la pulpa. Sims et al. (1963) concluyen que la firmeza de la pulpa y el color de fondo pueden ser utilizados como índices de madurez para melocotones y sugiere la utilización del colorímetro "tristimulus" para el desarrollo de una carta colorimétrica para la evaluación de madurez. De la Plaza (1973), define los conceptos de" umbral de color de fondo" y de "umbral de firmeza" como licites superior e inferior, respectivamente, para la madurez de consumo de pera "Dr Jules Guyot". Delwiche et al. (1983,1985,1987) estudian las correlaciones entre la firmeza de la pulpa y otros índices de madurez y el color de fondo proponiendo una carta de colores y madurez. Martínez-Javega y Otero (1989) concluyen que hay una correlación entre la firmeza de la pulpa y el índice cromático obtenido a partir de las coordenadas Hunter L,a,b para el aguacate cv. "HASS" con distintos tratamientos de conservación. Estudios semejantes para otros productos se han referido: Bittner y Noris (1986); Long et al. (1973); Kramer (1976); Robbins y Moore (1990); Pai y Sastry(1990)L El presente trabajo tiene como objetivo estudiar los índices cromáticos obtenidos a partir de las coordenadas Hunter L, a, b para la piel y la pulpa de aguacate cv." HASS ".

Análisis de las diferencias en el crecimiento vegetativo, en los componentes del rendimiento y en la composición de las uvas y vinos del cv. cabernet sauvignon (vitis vinifera L.), producidos por el déficit hídrico durante pre-envero y post-envero

La presente Tesis Doctoral se realizó con el fin de estimar conjuntamente la respuesta agronómica y fisiológica de la vid (Vitis vinifera L.), así como los efectos sobre la evolución de la maduración, composición y la calidad de la uva y del vino, bajo la aplicación de diferentes déficit hídricos en pre-envero y post-envero, dentro de un marco de referencia de cambio climático. La variación climática que prevén los estudios sobre el cambio climático, resulta un factor decisivo en la eficiencia del uso del agua en la vid. En zonas cálidas, las estrategias de cultivo del viñedo frente al cambio climático deben de ir dirigidas a atenuar sus efectos sobre el crecimiento y el desarrollo de la vid, haciéndose imprescindible el estudio pormenorizado del déficit hídrico como factor decisivo en la obtención de las uvas adecuadas, ya que son la clave indispensable para el éxito en la elaboración del vino, y de forma muy especial en los vinos enfocados a un sector de alta calidad. El ensayo se llevó a cabo en un viñedo comercial de Bodegas Licinia, en la Comunidad de Madrid, durante los años 2010 y 2011. La variedad estudiada fue Cabernet sauvignon / 41 B, plantada a un marco de plantación de 3 m x 1 m, con un guiado vertical de la vegetación. El dispositivo experimental fue totalmente al azar, y se establecieron 4 tratamientos experimentales con 4 grados de disponibilidad hídrica, déficit moderado continuo (T0,45-0,6), déficit severo continuo (T0-0,3), déficit severo después de envero (T0,45-0,3) y déficit severo antes de envero (T0-0,6). En cada tratamiento se distribuyeron 3 repeticiones. El año 2010 fue el más lluvioso de los años de ensayo, con 478 mm de precipitaciones anuales, lo que supuso 146 mm más que en el año 2011. Su distribución a lo largo del ciclo fue más homogénea en el año 2010, mientras que en 2011 las precipitaciones contabilizadas en el período de maduración de la uva fueron nulas. La temperatura media subió 0,9ºC en 2011, respecto a 2010 y en cuanto a la integral térmica eficaz, en 2011 se acumularon, desde el 1 de abril hasta el final de ciclo, 217 grados•día más que en 2010. El déficit hídrico en pre-envero, modificó notablemente el crecimiento vegetativo y la producción de cosecha de la parcela de ensayo, no así la fertilidad de las yemas. El tratamiento con mayor disponibilidad hídrica (T0,45-0,6) obtuvo el mayor peso de baya, y los tratamientos con menor déficit hídrico en pre-envero (T0,45-0,6 y T0,45-0,3) registraron los mayores rendimientos de cosecha, mientras que las menores tasas de cuajado correspondieron al tratamiento con un déficit severo continuo (T0-0,3). La parcela de ensayo se caracterizó por un exceso de vigor y un alto crecimiento vegetativo. El pH del mosto se vio afectado por el déficit hídrico, disminuyendo su valor en el tratamiento de déficit hídrico severo antes de envero (T0-0,6). Organolépticamente, no se percibieron diferencias significativas en los vinos elaborados en función del déficit hídrico, y respecto a su composición físico-química, solo existieron diferencias en la concentración de ácido L-Málico, con mayores concentraciones en los tratamientos sin déficit hídrico en pre-envero, T0,45-0,6 y T0,45-0,3. El déficit hídrico modificó notablemente el color del vino, aumentando los valores de las coordenadas CIELAB a* y b*, la luminosidad (L*), croma (C*) y tonalidad (H*), para los tratamientos con un déficit severo en pre-envero (T0-0,3 y T0-0,6) y disminuyendo estas en el tratamiento con mayor disponibilidad hídrica (T0,45-0,6). Del mismo modo, mediante el análisis de color por métodos tradicionales, IPT e IC de los vinos, aumentó en los tratamientos con mayor déficit hídrico en pre-envero (T0-0,3 y T0-0,6), respecto a los tratamientos de mayor disponibilidad (T0,45-0,6 y T0,45-0,3). La concentración de taninos de la baya en vendimia, no se vio afectada por el déficit hídrico, aunque sí estuvo relacionada positivamente con el tamaño de las bayas. Organolépticamente, los hollejos del año 2011 resultaron con menor frescura, acidez, afrutado, sensación herbácea e intensidad tánica, aunque con mayor astringencia respecto a 2010. Las pepitas fueron más astringentes y aromáticas pero menos crujientes, sin llegar a los niveles de madurez del año 2010. El catador relacionó los taninos con la calidad del vino, asociándolos con un mayor cuerpo, acidez, intensidad, equilibrio gustativo, amargor y menor astringencia en la fase gustativa. La concentración de taninos en los vinos se vio favorecida con el déficit hídrico en pre-envero y post-envero. Los tratamientos con mayor déficit hídrico en pre-envero, T0-0,6 y T0-0,3, obtuvieron las menores concentraciones de potasio en mostos y vinos. Las relaciones entre la concentración de potasio, ácido L-Málico y el porcentaje de color rojo puro (dA(%)) resultaron altamente significativas, de modo que las mayores tasas de potasio en el vino se asociaron a los valores más bajos de color rojo y a los mayores de ácido L-Málico. ABSTRACT The present Doctoral Thesis has been done in order to estimate the grapevine (Vitis vinifera L.) agronomic and physiologic performance or response as well as the impact in the grape and wine maturity, composition and quality evolution, with different water deficits. The variation in climate that the global warming studies for seen is a key factor for the grapevine water use efficiency. In warm areas the farming vineyards strategy to face the climatic change, should be focused on diminish the effects on the grapevine growth and development, so that the water deficit detailed analysis becomes decisive to obtain the appropriate grapes, that are the main subject for a successful wine production and especially for top quality wines. The trial was carried out in a commercial vineyard in Chinchón (Madrid), Licinia winery, during the 2010 and 2011 seasons. The grape variety studied was Cabernet Sauvignon grafted onto 41B with a vine spacing 3m x 1m trained as VSP. Experimental design consisted on 4 irrigation treatments with 3 replications totally randomized. Irrigation treatments were: moderate regulated deficit (T0,45-0,6), severe continuous deficit (T0-0,3), severe post-veraison deficit (T0,45-0,3) and severe pre veraison deficit (T0-0,6). The 2010 was rainier year than the 2011; Total annual rain in 2010 was 478 mm, which resulted in 146 mm more than in 2011. The distribution along the vine cycle was more homogeneous in the 2010, whereas precipitations in 2011 along the grape maturity period were nonexistent. The average temperature in 2011 was 0,9ºC higher than that of the 2010 and regarding to the thermal integral, in the 2011 from 1st April to the end of the growing cycle, was 217 degrees•day higher than that in 2010. Water deficit significantly modified the vegetative growth and yield but, it did not modified bud fertility. The treatment with the highest water availability (T0,45-0,6) got the highest berry size, the lowest berry set rates were found in the severe continuous deficit treatment (T0-0,3). The plot studied in this trial was characterized by both excessive vigour and vegetative growth. Water deficit modified the pH must by, reducing it in the severe water deficit during pre-veraison (T0-0,6). There were not differences in wine tasting between the water deficits treatments. Regarding to the physical-chemical composition, it only existed differences in the L-malic acid concentration, resulting higher concentrations in the water deficit pre-veraison treatments: T0,45-0,6 y T0,45-0,3. Water deficit significantly modified wine colour by, increasing the CIELAB coordinates a* and b*, the brightness (L*), croma (C*) and tonality (H*), in the lower water availability pre-veraison treatments (T0-0,3 y T0-0,6), and reducing them in the in the moderate continuous water deficit ones (T0,45-0,6). By means of traditional wine colour parameters analyses, red colour percentage, TPI, they became higher in the lower water availability pre-veraison treatments (T0-0,3 y T0-0,6), than in those with higher availability (T0,45-0,6 y T0,45-0,3). At harvest, berry tannins concentrations was not affected by the water deficit although it did in a positive way, in the berry size. Berry tasting in 2011, resulted in a lower freshness, acidity, fruity, herbaceous flavour and tannic intensity, but with higher astringency respect to the 2010 season. Seeds, in 2011, were more astringent and aromatic as in the 2010, but less crunchy, without getting to the point of maturity. The taster linked the tannins to wine quality, associating them with a higher bodiest wine, acidity, intensity, taste balance, bitterness and with a lower astringency in the tasting stage. Treatments with a higher water deficit up to veraison T0-0,6 y T0-0,3 got less musts and wines potassium concentration. The relation between L-malic acid and the full red color percentage (dA(%)), were highly related, resulting the higher potassium content the lower wine quality.

System supporting public transport based on geolocation

El propósito de este proyecto de fin de Grado es el estudio y desarrollo de una aplicación basada en Android que proporcionará soporte y atención a los servicios de transporte público existentes en Cracovia, Polonia. La principal funcionalidad del sistema será consultar la posición de un determinado autobús o tranvía y mostrar su ubicación con exactitud. Para lograr esto, necesitaremos tres fases de desarrollo. En primer lugar, deberemos implementar un sistema que obtenga las coordenadas geográficas de los vehículos de transporte público en cada instante. A continuación, tendremos que registrar todos estos datos y almacenarlos en una base de datos en un servidor web. Por último, desarrollaremos un sistema cliente que realice consultas a tiempo real sobre estos datos almacenados, obteniendo la posición para una línea determinada y mostrando su ubicación con un marcador en el mapa. Para hacer el seguimiento de los vehículos, sería necesario tener acceso a una API pública que nos proporcionase la posición registrada por los GPS que integran cada uno de ellos. Como esta API no existe actualmente para los servicios de autobús, y para los tranvías es de uso meramente privado, desarrollaremos una segunda aplicación en Android que hará las funciones del lado servidor. En ella podremos elegir mediante una simple interfaz el número de línea y un código específico que identificará a cada vehículo en particular (e.g. podemos tener 6 tranvías recorriendo la red al mismo tiempo para la línea 24). Esta aplicación obtendrá las coordenadas geográficas del teléfono móvil, lo cual incluye latitud, longitud y orientación a través del proveedor GPS. De este modo, podremos realizar una simulación de como el sistema funcionará a tiempo real utilizando la aplicación servidora desde dentro de un tranvía o autobús y, al mismo tiempo, utilizando la aplicación cliente haciendo peticiones para mostrar la información de dicho tranvía. El cliente, además, podrá consultar la ruta de cualquier línea sin necesidad de tener acceso a Internet. Almacenaremos las rutas y paradas de cada línea en la memoria del teléfono móvil utilizando ficheros XML debido al poco espacio que ocupan y a lo útil que resulta poder consultar un trayecto en cualquier momento, independientemente del acceso a la red. El usuario también podrá consultar las tablas de horarios oficiales para cada línea. Aunque en este caso si será necesaria una conexión a Internet debido a que se realizará a través de la web oficial de MPK. Para almacenar todas las coordenadas de cada vehículo en cada instante necesitaremos crear una base de datos en un servidor. Esto se resolverá mediante el uso de MYSQL y PHP. Se enviarán peticiones de tipo GET y POST a los servicios PHP que se encargarán de traducir y realizar la consulta correspondiente a la base de datos MYSQL. Por último, gracias a todos los datos recogidos relativos a la posición de los vehículos de transporte público, podremos realizar algunas tareas de análisis. Comparando la hora exacta a la que los vehículos pasaron por cada parada y la hora a la que deberían haber pasado según los horarios oficiales, podremos descubrir fallos en estos. Seremos capaces de determinar si es un error puntual debido a factores externos (atascos, averías,…) o si por el contrario, es algo que ocurre muy a menudo y se debería corregir el horario oficial. ABSTRACT The aim of this final Project (for University) is to develop an Android application thatwill provide support and feedback to the public transport services in Krakow. The main functionality of the system will be to track the position of a desired bus or tram line, and display its position on the map. To achieve this, we will need 3 stages: the first one will be to implement a system that sends the geographical position of the public transport vehicles, the second one will be to collect this data in a web server, and the last one will be to get the last location registered for the desired line and display it on the map. For tracking the vehicles, we would need to have access to a public API that should be connected with each bus/tram GPS. As this doesn’t exist in Krakow or at least is not available for public use, we will develop a second android application that will do the server side job. We will be able to choose in a simple interface the line number and a code letter to identify each vehicle (e.g. we can have 6 trams that belong to the line number 24 working at the same time). It will take the current mobile geolocation; this includes getting latitude, longitude and bearing from the GPS provider. Thus, we will be able to make a simulation of how the system works in real time by using the server app inside a tram and at the same time, using the client app and making requests to display the information of that tram. The client will also be able to check the path of the desired line without internet access. We will store the path and stops for each line locally in the phone memory using xml files due to the few requirements of available space it needs and the usefulness of checking a path when needed. This app will also offer the functionality of checking the timetable for the line, but in this case, it will link to the official Mpk website, so Internet access will be required. For storing all the coordinates for each vehicle at every moment we will need to create a database on a server. We have decided that the easiest way is to use Mysql and PHP for the deployment of the service. We will send GET and POST requests to the php files and those files will make the according queries to our database. Finally, based on all the collected data, we will be able to get some information about errors in the system of public transport timetables. We will check at what time a line was in each specific stop and compare it with the official timetable to find mistakes of time. We will determine if it is something that happens occasionally and related to external factors (e.g. traffic jams, breakdowns…) or if on the other hand, it is something that happens very often and the public transport timetables should be looked over and corrected.

Diseño de una interfaz de representación gráfica de la base de datos de HRTF

El actual proyecto consiste en la creación de una interfaz gráfica de usuario (GUI) en entorno de MATLAB que realice una representación gráfica de la base de datos de HRTF (Head-Related Transfer Function). La función de transferencia de la cabeza es una herramienta muy útil en el estudio de la capacidad del ser humano para percibir su entorno sonoro, además de la habilidad de éste en la localización de fuentes sonoras en el espacio que le rodea. La HRTF biaural (terminología para referirse al conjunto de HRTF del oído izquierdo y del oído derecho) en sí misma, posee información de especial interés ya que las diferencias entre las HRTF de cada oído, conceden la información que nuestro sistema de audición utiliza en la percepción del campo sonoro. Por ello, la funcionalidad de la interfaz gráfica creada presenta gran provecho dentro del estudio de este campo. Las diferencias interaurales se caracterizan en amplitud y en tiempo, variando en función de la frecuencia. Mediante la transformada inversa de Fourier de la señal HRTF, se obtiene la repuesta al impulso de la cabeza, es decir, la HRIR (Head-Related Impulse Response). La cual, además de tener una gran utilidad en la creación de software o dispositivos de generación de sonido envolvente, se utiliza para obtener las diferencias ITD (Interaural Time Difference) e ILD (Interaural Time Difference), comúnmente denominados “parámetros de localización espacial”. La base de datos de HRTF contiene la información biaural de diferentes puntos de ubicación de la fuente sonora, formando una red de coordenadas esféricas que envuelve la cabeza del sujeto. Dicha red, según las medidas realizadas en la cámara anecoica de la EUITT (Escuela Universitaria de Ingeniería Técnica de Telecomunicación), presenta una precisión en elevación de 10º y en azimut de 5º. Los receptores son dos micrófonos alojados en el maniquí acústico llamado HATS (Hats and Torso Simulator) modelo 4100D de Brüel&Kjaer. Éste posee las características físicas que influyen en la percepción del entorno como son las formas del pabellón auditivo (pinna), de la cabeza, del cuello y del torso humano. Será necesario realizar los cálculos de interpolación para todos aquellos puntos no contenidos en la base de datos HRTF, este proceso es sumamente importante no solo para potenciar la capacidad de la misma sino por su utilidad para la comparación entre otras bases de datos existentes en el estudio de este ámbito. La interfaz gráfica de usuario está concebida para un manejo sencillo, claro y predecible, a la vez que interactivo. Desde el primer boceto del programa se ha tenido clara su filosofía, impuesta por las necesidades de un usuario que busca una herramienta práctica y de manejo intuitivo. Su diseño de una sola ventana reúne tanto los componentes de obtención de datos como los que hacen posible la representación gráfica de las HRTF, las HRIR y los parámetros de localización espacial, ITD e ILD. El usuario podrá ir alternando las representaciones gráficas a la vez que introduce las coordenadas de los puntos que desea visualizar, definidas por phi (elevación) y theta (azimut). Esta faceta de la interfaz es la que le otorga una gran facilidad de acceso y lectura de la información representada en ella. Además, el usuario puede introducir valores incluidos en la base de datos o valores intermedios a estos, de esta manera, se indica a la interfaz la necesidad de realizar la interpolación de los mismos. El método de interpolación escogido es el de la ponderación de la distancia inversa entre puntos. Dependiendo de los valores introducidos por el usuario se realizará una interpolación de dos o cuatro puntos, siendo éstos limítrofes al valor introducido, ya sea de phi o theta. Para añadir versatilidad a la interfaz gráfica de usuario, se ha añadido la opción de generar archivos de salida en forma de imagen de las gráficas representadas, de tal forma que el usuario pueda extraer los datos que le interese para cualquier valor de phi y theta. Se completa el presente proyecto fin de carrera con un trabajo de investigación y estudio comparativo de la función y la aplicación de las bases de datos de HRTF dentro del marco científico y de investigación. Esto ha hecho posible concentrar información relacionada a través de revistas científicas de investigación como la JAES (Journal of the Audio Engineering Society) o la ASA (Acoustical Society of America), además, del IEEE ( Institute of Electrical and Electronics Engineers) o la “Web of knowledge” entre otras. Además de realizar la búsqueda en estas fuentes, se ha optado por vías de información más comunes como Google Académico o el portal de acceso “Ingenio” a los todos los recursos electrónicos contenidos en la base de datos de la universidad. El estudio genera una ampliación en el conocimiento de la labor práctica de las HRTF. La mayoría de los estudios enfocan sus esfuerzos en mejorar la percepción del evento sonoro mediante su simulación en la escucha estéreo o multicanal. A partir de las HRTF, esto es posible mediante el análisis y el cálculo de datos como pueden ser las regresiones, siendo éstas muy útiles en la predicción de una medida basándose en la información de la actual. Otro campo de especial interés es el de la generación de sonido 3D. Mediante la base de datos HRTF es posible la simulación de una señal biaural. Se han diseñado algoritmos que son implementados en dispositivos DSP, de tal manera que por medio de retardos interaurales y de diferencias espectrales es posible llegar a un resultado óptimo de sonido envolvente, sin olvidar la importancia de los efectos de reverberación para conseguir un efecto creíble de sonido envolvente. Debido a la complejidad computacional que esto requiere, gran parte de los estudios coinciden en desarrollar sistemas más eficientes, llegando a objetivos tales como la generación de sonido 3D en tiempo real. ABSTRACT. This project involves the creation of a Graphic User Interface (GUI) in the Matlab environment which creates a graphic representation of the HRTF (Head-Related Transfer Function) database. The head transfer function is a very useful tool in the study of the capacity of human beings to perceive their sound environment, as well as their ability to localise sound sources in the area surrounding them. The binaural HRTF (terminology which refers to the HRTF group of the left and right ear) in itself possesses information of special interest seeing that the differences between the HRTF of each ear admits the information that our system of hearing uses in the perception of each sound field. For this reason, the functionality of the graphic interface created presents great benefits within the study of this field. The interaural differences are characterised in space and in time, varying depending on the frequency. By means of Fourier's transformed inverse of the HRTF signal, the response to the head impulse is obtained, in other words, the HRIR (Head-Related Impulse Response). This, as well as having a great use in the creation of software or surround sound generating devices, is used to obtain ITD differences (Interaural Time Difference) and ILD (Interaural Time Difference), commonly named “spatial localisation parameters”. The HRTF database contains the binaural information of different points of sound source location, forming a network of spherical coordinates which surround the subject's head. This network, according to the measures carried out in the anechoic chamber at the EUITT (School of Telecommunications Engineering) gives a precision in elevation of 10º and in azimuth of 5º. The receivers are two microphones placed on the acoustic mannequin called HATS (Hats and Torso Simulator) Brüel&Kjaer model 4100D. This has the physical characteristics which affect the perception of the surroundings which are the forms of the auricle (pinna), the head, neck and human torso. It will be necessary to make interpolation calculations for all those points which are not contained the HRTF database. This process is extremely important not only to strengthen the database's capacity but also for its usefulness in making comparisons with other databases that exist in the study of this field. The graphic user interface is conceived for a simple, clear and predictable use which is also interactive. Since the first outline of the program, its philosophy has been clear, based on the needs of a user who requires a practical tool with an intuitive use. Its design with only one window unites not only the components which obtain data but also those which make the graphic representation of the HRTFs possible, the hrir and the ITD and ILD spatial location parameters. The user will be able to alternate the graphic representations at the same time as entering the point coordinates that they wish to display, defined by phi (elevation) and theta (azimuth). The facet of the interface is what provides the great ease of access and reading of the information displayed on it. In addition, the user can enter values included in the database or values which are intermediate to these. It is, likewise, indicated to the interface the need to carry out the interpolation of these values. The interpolation method is the deliberation of the inverse distance between points. Depending on the values entered by the user, an interpolation of two or four points will be carried out, with these being adjacent to the entered value, whether that is phi or theta. To add versatility to the graphic user interface, the option of generating output files in the form of an image of the graphics displayed has been added. This is so that the user may extract the information that interests them for any phi and theta value. This final project is completed with a research and comparative study essay on the function and application of HRTF databases within the scientific and research framework. It has been possible to collate related information by means of scientific research magazines such as the JAES (Journal of the Audio Engineering Society), the ASA (Acoustical Society of America) as well as the IEEE (Institute of Electrical and Electronics Engineers) and the “Web of knowledge” amongst others. In addition to carrying out research with these sources, I also opted to use more common sources of information such as Academic Google and the “Ingenio” point of entry to all the electronic resources contained on the university databases. The study generates an expansion in the knowledge of the practical work of the HRTF. The majority of studies focus their efforts on improving the perception of the sound event by means of its simulation in stereo or multichannel listening. With the HRTFs, this is possible by means of analysis and calculation of data as can be the regressions. These are very useful in the prediction of a measure being based on the current information. Another field of special interest is that of the generation of 3D sound. Through HRTF databases it is possible to simulate the binaural signal. Algorithms have been designed which are implemented in DSP devices, in such a way that by means of interaural delays and wavelength differences it is possible to achieve an excellent result of surround sound, without forgetting the importance of the effects of reverberation to achieve a believable effect of surround sound. Due to the computational complexity that this requires, a great many studies agree on the development of more efficient systems which achieve objectives such as the generation of 3D sound in real time.

Simulación en tiempo real de vehículos industriales con modelos multicuerpo de gran complejidad

El objetivo de esta Tesis ha sido la consecución de simulaciones en tiempo real de vehículos industriales modelizados como sistemas multicuerpo complejos formados por sólidos rígidos. Para el desarrollo de un programa de simulación deben considerarse cuatro aspectos fundamentales: la modelización del sistema multicuerpo (tipos de coordenadas, pares ideales o impuestos mediante fuerzas), la formulación a utilizar para plantear las ecuaciones diferenciales del movimiento (coordenadas dependientes o independientes, métodos globales o topológicos, forma de imponer las ecuaciones de restricción), el método de integración numérica para resolver estas ecuaciones en el tiempo (integradores explícitos o implícitos) y finalmente los detalles de la implementación realizada (lenguaje de programación, librerías matemáticas, técnicas de paralelización). Estas cuatro etapas están interrelacionadas entre sí y todas han formado parte de este trabajo. Desde la generación de modelos de una furgoneta y de camión con semirremolque, el uso de tres formulaciones dinámicas diferentes, la integración de las ecuaciones diferenciales del movimiento mediante métodos explícitos e implícitos, hasta el uso de funciones BLAS, de técnicas de matrices sparse y la introducción de paralelización para utilizar los distintos núcleos del procesador. El trabajo presentado en esta Tesis ha sido organizado en 8 capítulos, dedicándose el primero de ellos a la Introducción. En el Capítulo 2 se presentan dos formulaciones semirrecursivas diferentes, de las cuales la primera está basada en una doble transformación de velocidades, obteniéndose las ecuaciones diferenciales del movimiento en función de las aceleraciones relativas independientes. La integración numérica de estas ecuaciones se ha realizado con el método de Runge-Kutta explícito de cuarto orden. La segunda formulación está basada en coordenadas relativas dependientes, imponiendo las restricciones por medio de penalizadores en posición y corrigiendo las velocidades y aceleraciones mediante métodos de proyección. En este segundo caso la integración de las ecuaciones del movimiento se ha llevado a cabo mediante el integrador implícito HHT (Hilber, Hughes and Taylor), perteneciente a la familia de integradores estructurales de Newmark. En el Capítulo 3 se introduce la tercera formulación utilizada en esta Tesis. En este caso las uniones entre los sólidos del sistema se ha realizado mediante uniones flexibles, lo que obliga a imponer los pares por medio de fuerzas. Este tipo de uniones impide trabajar con coordenadas relativas, por lo que la posición del sistema y el planteamiento de las ecuaciones del movimiento se ha realizado utilizando coordenadas Cartesianas y parámetros de Euler. En esta formulación global se introducen las restricciones mediante fuerzas (con un planteamiento similar al de los penalizadores) y la estabilización del proceso de integración numérica se realiza también mediante proyecciones de velocidades y aceleraciones. En el Capítulo 4 se presenta una revisión de las principales herramientas y estrategias utilizadas para aumentar la eficiencia de las implementaciones de los distintos algoritmos. En primer lugar se incluye una serie de consideraciones básicas para aumentar la eficiencia numérica de las implementaciones. A continuación se mencionan las principales características de los analizadores de códigos utilizados y también las librerías matemáticas utilizadas para resolver los problemas de álgebra lineal tanto con matrices densas como sparse. Por último se desarrolla con un cierto detalle el tema de la paralelización en los actuales procesadores de varios núcleos, describiendo para ello el patrón empleado y las características más importantes de las dos herramientas propuestas, OpenMP y las TBB de Intel. Hay que señalar que las características de los sistemas multicuerpo problemas de pequeño tamaño, frecuente uso de la recursividad, y repetición intensiva en el tiempo de los cálculos con fuerte dependencia de los resultados anteriores dificultan extraordinariamente el uso de técnicas de paralelización frente a otras áreas de la mecánica computacional, tales como por ejemplo el cálculo por elementos finitos. Basándose en los conceptos mencionados en el Capítulo 4, el Capítulo 5 está dividido en tres secciones, una para cada formulación propuesta en esta Tesis. En cada una de estas secciones se describen los detalles de cómo se han realizado las distintas implementaciones propuestas para cada algoritmo y qué herramientas se han utilizado para ello. En la primera sección se muestra el uso de librerías numéricas para matrices densas y sparse en la formulación topológica semirrecursiva basada en la doble transformación de velocidades. En la segunda se describe la utilización de paralelización mediante OpenMP y TBB en la formulación semirrecursiva con penalizadores y proyecciones. Por último, se describe el uso de técnicas de matrices sparse y paralelización en la formulación global con uniones flexibles y parámetros de Euler. El Capítulo 6 describe los resultados alcanzados mediante las formulaciones e implementaciones descritas previamente. Este capítulo comienza con una descripción de la modelización y topología de los dos vehículos estudiados. El primer modelo es un vehículo de dos ejes del tipo chasis-cabina o furgoneta, perteneciente a la gama de vehículos de carga medianos. El segundo es un vehículo de cinco ejes que responde al modelo de un camión o cabina con semirremolque, perteneciente a la categoría de vehículos industriales pesados. En este capítulo además se realiza un estudio comparativo entre las simulaciones de estos vehículos con cada una de las formulaciones utilizadas y se presentan de modo cuantitativo los efectos de las mejoras alcanzadas con las distintas estrategias propuestas en esta Tesis. Con objeto de extraer conclusiones más fácilmente y para evaluar de un modo más objetivo las mejoras introducidas en la Tesis, todos los resultados de este capítulo se han obtenido con el mismo computador, que era el top de la gama Intel Xeon en 2007, pero que hoy día está ya algo obsoleto. Por último los Capítulos 7 y 8 están dedicados a las conclusiones finales y las futuras líneas de investigación que pueden derivar del trabajo realizado en esta Tesis. Los objetivos de realizar simulaciones en tiempo real de vehículos industriales de gran complejidad han sido alcanzados con varias de las formulaciones e implementaciones desarrolladas. ABSTRACT The objective of this Dissertation has been the achievement of real time simulations of industrial vehicles modeled as complex multibody systems made up by rigid bodies. For the development of a simulation program, four main aspects must be considered: the modeling of the multibody system (types of coordinates, ideal joints or imposed by means of forces), the formulation to be used to set the differential equations of motion (dependent or independent coordinates, global or topological methods, ways to impose constraints equations), the method of numerical integration to solve these equations in time (explicit or implicit integrators) and the details of the implementation carried out (programming language, mathematical libraries, parallelization techniques). These four stages are interrelated and all of them are part of this work. They involve the generation of models for a van and a semitrailer truck, the use of three different dynamic formulations, the integration of differential equations of motion through explicit and implicit methods, the use of BLAS functions and sparse matrix techniques, and the introduction of parallelization to use the different processor cores. The work presented in this Dissertation has been structured in eight chapters, the first of them being the Introduction. In Chapter 2, two different semi-recursive formulations are shown, of which the first one is based on a double velocity transformation, thus getting the differential equations of motion as a function of the independent relative accelerations. The numerical integration of these equations has been made with the Runge-Kutta explicit method of fourth order. The second formulation is based on dependent relative coordinates, imposing the constraints by means of position penalty coefficients and correcting the velocities and accelerations by projection methods. In this second case, the integration of the motion equations has been carried out by means of the HHT implicit integrator (Hilber, Hughes and Taylor), which belongs to the Newmark structural integrators family. In Chapter 3, the third formulation used in this Dissertation is presented. In this case, the joints between the bodies of the system have been considered as flexible joints, with forces used to impose the joint conditions. This kind of union hinders to work with relative coordinates, so the position of the system bodies and the setting of the equations of motion have been carried out using Cartesian coordinates and Euler parameters. In this global formulation, constraints are introduced through forces (with a similar approach to the penalty coefficients) are presented. The stabilization of the numerical integration is carried out also by velocity and accelerations projections. In Chapter 4, a revision of the main computer tools and strategies used to increase the efficiency of the implementations of the algorithms is presented. First of all, some basic considerations to increase the numerical efficiency of the implementations are included. Then the main characteristics of the code’ analyzers used and also the mathematical libraries used to solve linear algebra problems (both with dense and sparse matrices) are mentioned. Finally, the topic of parallelization in current multicore processors is developed thoroughly. For that, the pattern used and the most important characteristics of the tools proposed, OpenMP and Intel TBB, are described. It needs to be highlighted that the characteristics of multibody systems small size problems, frequent recursion use and intensive repetition along the time of the calculation with high dependencies of the previous results complicate extraordinarily the use of parallelization techniques against other computational mechanics areas, as the finite elements computation. Based on the concepts mentioned in Chapter 4, Chapter 5 is divided into three sections, one for each formulation proposed in this Dissertation. In each one of these sections, the details of how these different proposed implementations have been made for each algorithm and which tools have been used are described. In the first section, it is shown the use of numerical libraries for dense and sparse matrices in the semirecursive topological formulation based in the double velocity transformation. In the second one, the use of parallelization by means OpenMP and TBB is depicted in the semi-recursive formulation with penalization and projections. Lastly, the use of sparse matrices and parallelization techniques is described in the global formulation with flexible joints and Euler parameters. Chapter 6 depicts the achieved results through the formulations and implementations previously described. This chapter starts with a description of the modeling and topology of the two vehicles studied. The first model is a two-axle chassis-cabin or van like vehicle, which belongs to the range of medium charge vehicles. The second one is a five-axle vehicle belonging to the truck or cabin semi-trailer model, belonging to the heavy industrial vehicles category. In this chapter, a comparative study is done between the simulations of these vehicles with each one of the formulations used and the improvements achieved are presented in a quantitative way with the different strategies proposed in this Dissertation. With the aim of deducing the conclusions more easily and to evaluate in a more objective way the improvements introduced in the Dissertation, all the results of this chapter have been obtained with the same computer, which was the top one among the Intel Xeon range in 2007, but which is rather obsolete today. Finally, Chapters 7 and 8 are dedicated to the final conclusions and the future research projects that can be derived from the work presented in this Dissertation. The objectives of doing real time simulations in high complex industrial vehicles have been achieved with the formulations and implementations developed.

Optimización de la Exactitud Geométrica al Integrar Diversas Fuentes de Datos en Proyectos de Inventario Forestal Basados en Escaneo Laser Aerotransportado=Optimizing Geometric Accuracy in Data Source Integration for Airborne Laser Scanning-based Forest Inventory Projects

La mayoría de las aplicaciones forestales del escaneo laser aerotransportado (ALS, del inglés airborne laser scanning) requieren la integración y uso simultaneo de diversas fuentes de datos, con el propósito de conseguir diversos objetivos. Los proyectos basados en sensores remotos normalmente consisten en aumentar la escala de estudio progresivamente a lo largo de varias fases de fusión de datos: desde la información más detallada obtenida sobre un área limitada (la parcela de campo), hasta una respuesta general de la cubierta forestal detectada a distancia de forma más incierta pero cubriendo un área mucho más amplia (la extensión cubierta por el vuelo o el satélite). Todas las fuentes de datos necesitan en ultimo termino basarse en las tecnologías de sistemas de navegación global por satélite (GNSS, del inglés global navigation satellite systems), las cuales son especialmente erróneas al operar por debajo del dosel forestal. Otras etapas adicionales de procesamiento, como la ortorectificación, también pueden verse afectadas por la presencia de vegetación, deteriorando la exactitud de las coordenadas de referencia de las imágenes ópticas. Todos estos errores introducen ruido en los modelos, ya que los predictores se desplazan de la posición real donde se sitúa su variable respuesta. El grado por el que las estimaciones forestales se ven afectadas depende de la dispersión espacial de las variables involucradas, y también de la escala utilizada en cada caso. Esta tesis revisa las fuentes de error posicional que pueden afectar a los diversos datos de entrada involucrados en un proyecto de inventario forestal basado en teledetección ALS, y como las propiedades del dosel forestal en sí afecta a su magnitud, aconsejando en consecuencia métodos para su reducción. También se incluye una discusión sobre las formas más apropiadas de medir exactitud y precisión en cada caso, y como los errores de posicionamiento de hecho afectan a la calidad de las estimaciones, con vistas a una planificación eficiente de la adquisición de los datos. La optimización final en el posicionamiento GNSS y de la radiometría del sensor óptico permitió detectar la importancia de este ultimo en la predicción de la desidad relativa de un bosque monoespecífico de Pinus sylvestris L. ABSTRACT Most forestry applications of airborne laser scanning (ALS) require the integration and simultaneous use of various data sources, pursuing a variety of different objectives. Projects based on remotely-sensed data generally consist in upscaling data fusion stages: from the most detailed information obtained for a limited area (field plot) to a more uncertain forest response sensed over a larger extent (airborne and satellite swath). All data sources ultimately rely on global navigation satellite systems (GNSS), which are especially error-prone when operating under forest canopies. Other additional processing stages, such as orthorectification, may as well be affected by vegetation, hence deteriorating the accuracy of optical imagery’s reference coordinates. These errors introduce noise to the models, as predictors displace from their corresponding response. The degree to which forest estimations are affected depends on the spatial dispersion of the variables involved and the scale used. This thesis reviews the sources of positioning errors which may affect the different inputs involved in an ALS-assisted forest inventory project, and how the properties of the forest canopy itself affects their magnitude, advising on methods for diminishing them. It is also discussed how accuracy should be assessed, and how positioning errors actually affect forest estimation, toward a cost-efficient planning for data acquisition. The final optimization in positioning the GNSS and optical image allowed to detect the importance of the latter in predicting relative density in a monospecific Pinus sylvestris L. forest.

Levantamiento y red topográfica del parque "Juan Pablo II".

El objeto del presente proyecto es dotar al parque público “Juan Pablo II”, de una cartografía a escala 1/250, además de realizar el modelo digital del terreno de dicho parque. El trabajo va a consistir en un levantamiento topográfico del sector oeste del parque público “Juan Pablo II”. Dado que el parque en toda su dimensión se encuentra dividido en dos por la Avenida de Machupichu, siendo la oeste la que forma parte de este proyecto. El levantamiento se realizará mediante el método de observación G.N.S.S. denominado R.T.K. (Real Time Kinematic) que consiste en la obtención de coordenadas en tiempo real. Este levantamiento se realizará apoyándose en puntos de una red implantada, a la que se dotará de coordenadas por el método de observación G.N.S.S. denominado “Estático” y por metodología clásica. A la red se le dotará de coordenadas en el Sistema de Referencia Oficial. Se realizará la cartografía del parque a partir de los datos del levantamiento y con el uso de programas de CAD. Finalmente, se realizará un modelo digital del terreno. El parque público “Juan Pablo II” se encuentra situado en barrio de La Piovera, en el término municipal de Madrid, provincia de Madrid. El parque ocupa 6 hectáreas y es el elegido para nuestro proyecto. Este parque está limitado por el norte por la Avda. del Papa Negro, por el sur está limitado por la Avda. de las Piceas, por el este por la Avda. Machupichu, por el oeste por un pinar. 1.- Estudio de Redes para la obtención de la Red Básica: Por metodología G.N.S.S., se ha realizado cuatro simulaciones, teniendo como criterios: Longitud de la líneas base; Tiempo de observación; Precisiones obtenidas. Y comprobación de la Red Básica por metodología Clásica. 2.- Creación de aplicación para representar redes en Google Earth (Geodesicas-kml). Ante la falta de aplicaciones libres para esta opción de difusión de redes de manera fácil y gratuita, se creó esta aplicación. 3.- Diseño cartográfico, mapas a escala 1/250. Esta escala fue elegida porque se quería tener un buen detalle para la creación del MDT. 4.- Creación de MDT, y adaptado a visualización en Google Earth. Se eligió esta plataforma para que no se necesitase ninguna aplicación de pago y por su gran difusión, facilitando la divulgación de los MDT. 5.- Paseo virtual en el interior del parque. Película de un recorrido dentro del parque con todos los elementos más representativos.

Vial entre el municipio de Móstoles y el polígono industrial Valdearenal(Arroyomolinos).

el proyecto consiste en el diseño de un vial que une el municipio de Móstoles y el Polígono Industrial Valdearenal, perteneciente al término municipal de Arroyomolinos. Adicionalmente se pretende diseñar una glorieta en su enlace con las vías urbanas de Móstoles para enlazar el vial con las vías urbanas de Móstoles. El enlace con el polígono industrial se realiza en una glorieta ya presente en la zona. Para la realización de este proyecto se ha seguido un proceso que engloba varias fases. En primer lugar, se ha creado una red que enlaza con las redes geodésicas nacionales con el fin de crear un grupo de vértices iniciales (bases) con coordenadas suficientemente precisas para basar en ellas el resto del trabajo y basadas en un sistema de coordenadas oficial. Desde las bases obtenidas se ha realizado un levantamiento en el que se ha dado una gran cantidad de puntos con el fin de definir el terreno de la zona lo más ajustado posible a la realidad y definiendo los elemento relevantes que se encuentran en este tales como calles, caminos, vallas, etc. Esta nube de puntos se ha medido en una franja transversal a la carretera planificada inicialmente de unos 60m, 30 a cada lado. La longitud total del tramo levantado es de aproximadamente 1,6km. Finalmente en campo, se ha realizado una nivelación entre las bases para que el desnivel entre estas sea lo más preciso posible y las diferentes zonas de puntos del levantamiento se ajusten la mejor posible entre sí. Realizada la parte de campo, se ha pasado a la parte de gabinete,donde se ha comenzado ajustando todos los resultados obtenidos para que la nube de puntos se acerque lo máximo que se pueda a la superficie real del terreno. Con estos datos se ha realizado la cartografía, dibujando todos los elementos de la zona y las líneas de nivel para definir la zona sobre el plano. Una vez hecha la cartografía se ha diseñado el vial y la glorieta nueva sobre esta y los enlaces con las zonas urbanas, ajustando esto a la normativa legal existente y estableciendo la señalización. Sobre el vial ya diseñado se ha realizado el cálculo de expropiaciones y, por último, el cálculo de presupuestos. El sistema de referencia en el que se encuentra todo el proyecto es ETRS89, que utiliza como modelo de elipsoide el GRS80. El sistema de coordenadas es UTM, ubicándose la zona del proyecto en el huso 30 y en el hemisferio norte. Se ha trabajado con alturas ortométricas, que emplean como datum el nivel medio del mar en Alicante entre los años 1870 y 1882. La razón por la que todos estos parámetros se han usado es debido a que son los oficiales en España.