Sueños proféticos en el Baladro de Burgos 1498

No parece necesario abundar sobre el peso específico del que gozó el sueño como vehículo de revelación divina dentro del sistema conceptual de las religiones grecorromana y judeocristiana; tampoco, sobre sus posibles productividad discursiva y significatividad ideológica para una tradición mucho más específica, la novela de caballerías castellana, entre cuyos resortes narrativos predilectos se cuentan los discursos de profetas, magos y adivinos. Son estos dos presupuestos los que explican la leve sorpresa de que El baladro del sabio Merlín con sus profecías del incunable burgalés de 1498 -que en su relativa brevedad como ejemplar del género caballeresco prodiga al lector ciento cincuenta textos proféticos- proponga sólo tres revelaciones oníricas, dos habidas por el rey Artur y una por el propio "sabio" Merlín; de hecho, la sección inicial del relato -que comprende hasta XVIII: 33 y que traduce una prosificación del hoy fragmentario Roman de Merlin en verso francés escrito por Robert de Boron- prescinde por completo de tal especie profética. El objetivo de esta comunicación es reflexionar brevemente sobre la funcionalidad y el sentido de esos tres casos que constituyen una raraavisdentro del enorme sistema proféticode nuestra novela.

Wais III. Indice de comprensión verbal : Descripción de una muestra para las edades 25-34 años en la ciudad de La Plata, según género

En el marco de la segunda etapa del Proyecto de Investigación para la Construcción de Normas para los subtests que conforman el índice de Comprensión Verbal (ICV) de la Escala de Inteligencia para Adultos de Wechsler-WAIS III, para la ciudad de La Plata, se analizan los resultados observados sobre una muestra probabilística de 381 sujetos, entre 25-34 años, de nivel terciario y universitario, según el género. En el subtest de Información el mejor desempeño se observó en los varones de ambos grupos de edad. En el subtest de Analogías, el mejor desempeño se observó en las mujeres del grupo de edad 25-29 años, seguido por el de los varones del grupo de edad 30-34 años.Las mujeres del grupo de edad 25-29 años mostraron el mejor desempeño en comprensión verbal. Sólo para este caso; el ICV calculado se ubicó levemente por encima del promedio estandarizado. Las mujeres del grupo de edad 30- 34 años, y los varones, en ambos grupos de edad, mostraron un ICV calculado levemente por debajo del promedio estandarizado. Las diferencias entre el puntaje promedio obtenido en cada uno de los subtests para ambos géneros en los dos grupos de edad, y el promedio esperado, no son estadísticamente significativos (p es menor que 0,05). La misma interpretación corresponde realizar para el ICV. El ICV calculado en las mujeres del grupo de edad 25-29 años, las ubica en el percentil 53 indicando que su desempeño específico observado en comprensión verbal, deja por debajo al 53de los sujetos de su grupo de referencia. Los rendimientos observados que dependen de los conocimientos adquiridos en la escuela y a través de la experiencia cultural, son de levemente superiores al promedio. Su capacidad para hacer uso del cuerpo de conocimiento general acumulado, para emitir juicios y resolver problemas tales como tareas de abstracción de las relaciones entre dos elementos distintos, o la de formar conceptos, son normales y alcanzan a superar el promedio del rendimiento de los sujetos tipificados para los mismos rangos de edad. El ICV calculado en las mujeres del grupo de edad 30-34 años, y en los varones de ambos grupos, ubica a los sujetos en el percentil 45 indicando que su desempeño específico observado en comprensión verbal, deja por debajo al 45de los sujetos de su grupo de referencia. Los rendimientos observados son normales aunque no llegan a superar el promedio del rendimiento de los sujetos tipificados para los mismos rangos de edad. El promedio de puntajes en cada uno de tres subtest que integran el ICV, así como el promedio en el ICV para ambos géneros, ubica a los sujetos dentro promedio esperado, según su grupos de edad (Paidós, 2002). Siguiendo los criterios establecidos en el Manual de la Prueba para la interpretación de las diferencias encontradas entre los puntajes e índices observados sobre una muestra de mujeres y varones instruidos, y los estandarizados, este Equipo considera que debería analizarse una muestra más amplia para determinar si dicha diferencia es inusual o resulta común en la población en general (Wechsler, 2002)

Principales efectos de la aplicación de una metodología destinada a inducir la muda en ponedoras comerciales sobre su producción cuantitativa y cualitativa durante el segundo ciclo de puesta

En la presente Tesis Doctoral se ha estudiado el efecto de una metodología para inducir la muda en dos estirpes de gallinas ponedoras comerciales, una ligera y otra semipesada, mediante el suministro de tres alimentos: salvado de trigo, cebada en grano y un pienso comercial de ponedoras aportado en cantidad restringida. Se realizaron dos pruebas en dos lotes diferentes de gallinas. En la primera de ellas se utilizaron 472 animales (236 de cada estirpe) alojados en jaulas con cuatro o seis gallinas por jaula, con una estructura factorial 2 x 3 x 2 (2 estirpes, 3 alimentos, 2 densidades) y una duración total de 32 semanas (4 de muda y 28 de producción posmuda). En la segunda prueba fueron 432 animales los utilizados (216 de cada estirpe), alojados en grupos de cuatro aves por jaula, con una estructura factorial 2 x 3 (2 estirpes, 3 alimentos) y una duración de 27 semanas (4 de muda y 23 de producción postmuda). En los dos experimentos realizados se estudió el efecto del uso de los alimentos citados para inducir de la muda sobre los resultados cuantitativos: pérdida de peso vivo durante la muda e intensidad de puesta, peso medio del huevo y masa de huevo diaria, durante y después de la muda, así como la distribución de la puesta en clases comerciales durante el segundo ciclo de puesta. Así como sobre los resultados cualitativos después de la muda: color de la cáscara de los huevos morenos, espesor de la cáscara, peso específico del huevo, altura del albumen, unidades Haugh y color de la yema. En la primera prueba se estudió, además, el efecto que la inducción de la muda mediante los tres alimentos considerados tuvo sobre la regresión del aparato reproductor de las gallinas durante este período de muda. En el primer experimento se observaron diferencias entre estirpes. Las gallinas ligeras tuvieron una más rápida regresión del aparato reproductor (ovario+oviducto) (P=0,003) que las semipesadas, aunque la regresión total no presentó diferencias significativas. Las gallinas semipesadas tuvieron mejores resultados después de la muda en intensidad de puesta (P<0,0001), en peso medio del huevo (P<0,0001) y en masa de huevo diaria (P=0,0002). También hubo diferencias significativas para las variables cualitativas espesor de la cáscara (mayor en huevos de gallinas semipesadas) mientras que los huevos procedentes de gallinas ligeras presentaron mejores valores de altura del albumen, de unidades Haugh y de color de yema; todas estas variables tuvieron un nivel de significación P<0,0001. El suministro restringido de pienso dio lugar a un mayor porcentaje de pérdida de peso vivo (P<0,0001) aunque la regresión de los órganos del aparato reproductor fue la más baja (P<0,003), no habiéndose encontrado diferencias entre los otros dos alimentos utilizados. Con este alimento también fue más lenta (P<0,0001) la disminución de la puesta durante la muda, aunque fue mayor la producción durante el segundo ciclo (P<0,0001). La única variable cualitativa afectada fue el espesor de cáscara (P<0,0001), con valores más altos en los huevos producidos por las gallinas mudadas con cebada. Los grupos de seis gallinas por jaula produjeron más huevos durante la muda (P<0,0001) aunque después de ésta la densidad de animales no tuvo efecto significativo, como tampoco lo hubo sobre los parámetros de calidad del huevo. En la segunda prueba las gallinas semipesadas experimentaron un menor porcentaje de pérdida de peso corporal (P<0,01) pero tuvieron mayores índices de puesta (P<0,001) y de huevos clasificables (P<0,001) durante la muda. En cambio, durante el segundo ciclo de producción las gallinas ligeras produjeron más huevos (P=0,0041), de menor peso (P<0,02) y con menor consumo de pienso (P<0,001). Los huevos puestos por las gallinas semipesadas tuvieron mayor espesor de cáscara y mayor color de yema, pero peor calidad de albumen (P<0,0001). La mayor pérdida de peso la experimentaron las gallinas mudadas con salvado (P<0,02). La producción durante la muda fue mayor (P<0,001) en las gallinas que consumieron pienso en cantidad restringida y también fueron las que tuvieron menor intensidad de puesta (P<0,006) y masa de huevo diaria (P<0,042) durante el segundo ciclo. El tratamiento de muda no tuvo efecto significativo sobre la calidad del huevo en esta segunda prueba. La principal conclusión que merece destacarse es que es posible inducir la muda a las gallinas ponedoras utilizando alimentos bajos en energía o en proteína, o altos en fibra, con un porcentaje de pérdida de peso vivo no tan alta como las recomendaciones tradicionales, y alcanzar buenos resultados productivos, tanto cuantitativos como cualitativos, durante el segundo ciclo de puesta. ABSTRACT Present Doctoral Thesis has studied the effect of a methodology to induce molting in two strains of commercial laying hens, one light and another semi-heavy one, through the provision of three feed: wheat bran, barley grain and a commercial laying hens feed provided in limited quantity. Two tests were performed in two different lots of layers. In the first 472 animals were used (236 of each strain) housed in cages with four or six hens per cage, with a structure 2 x 3 x 2 factorial (2 strains, 3 meals, 2 densities) and a total duration of 32 weeks (4 of molt and 28 of postmolting). In the second test 432 animals were used (216 each strain), housed in groups of four birds per cage, with a structure factorial 2 x 3 (2 strain, 3 meals) along 27 weeks (4 of molt and 23 of postmolting). In both experiments, we studied the effect of the use of above mentioned foods to induce molting on the quantitative results: body weight lost during the molt and laying index, average egg weight and egg mass daily during and after the molt, as well as on grading in commercial classes during the second laying cycle. As well as on qualitative outcomes after the molt: colour of Brown eggsshell, shell thickness, specific density, albumen height, Haugh units and yolk colour. In the first test was studied, in addition, the effect of induction of molting through the three feed considered on the regression of the reproductive tract of hens during molting period. In the first experiment, differences between strains were observed. Light hens had a faster regression of the reproductive tract (ovary+oviduct) (P=0,003) than semi-heavy hens, although the total regression did not present significant differences. Semi-heavy hens had better outcomes after the molt in laying index (P<0.0001), in average egg weight (P<0.0001) and daily egg mass (P=0, 0002). There were also significant differences for the qualitative variables (higher in semi-heavy hen eggs) as shell thickness while light chicken eggs showed better values of albumen height, Haugh units and yolk color; all these variables had the same level of significance (P<0.0001). Restricted supply of layer feed resulted in a greater percentage of live weight loss (P<0.0001) although the regression of the reproductive organs was the lowest (P<0.003), having not found differences between the other two feed used. With this food decreasing of laying during molting period was also slower (P<0.0001), although production was higher during the second cycle (P < 0.0001). The only qualitative variable affected was shell thickness (P<0.0001), with higher values in the eggs produced by hens molted with barley. Groups of six hens per cage produced more eggs during the moult (P<0.0001) but after this animal density had no significant effect, as neither had it on egg quality parameters. In the second trial hens semi-heavy experienced a lower percentage of body weight loss (P<0.01) but had higher rates of egg production (P<0.001) and of grading eggs (P<0.001) during the moult. On the other hand, during the second production cycle light hens produced more eggs (P=0,0041), of lower weight (P < 0.02) and with less feed intake (P<0.001). Eggs from semi-heavy hens had thicker shell and greater color yolk, but poorer quality of albumen (P<0.0001). The greater weight loss was experimented by the hens molted with wheat bran (P=0.02). Production during molting period was greater (P<0.001) in hens which consumed feed in restricted quantities and were also less the laying index (P=0.006) and daily egg mass (P=0.042) during the second cycle. The molting treatment had no significative effect on the quality of the egg in this second test. The main conclusion that deserves to stand out is that it is possible to induce molting hens using foods low in energy or in protein or high in fiber, with a percentage of live weight loss not as high as the traditional recommendations, and achieve good productive, both quantitative and qualitative results during the second implementation cycle.

Control y reducción de daños mecánicos en la manipulación de frutos cítricos

Durante las operaciones de recolección, transporte y manipulación de frutas y hortalizas se producen inevitablemente lesiones que progresivamente se acumulan sobre cada uno de los frutos, resultando disminuida su calidad, y por tanto su valor comercial. Los mercados actuales, tanto nacionales como internacionales, demandan frutas de alta calidad. Uno de los aspectos que constituyen esa calidad es la apariencia externa del producto, y dentro de ésta la ausencia de daños mecánicos, heridas, defectos de coloración etc. Es por tanto necesario revisar todos aquellos factores y procesos que afectan al producto, desde el momento de la recolección hasta la llegada del producto a manos del consumidor, con el fin de reducir las cuantiosas pérdidas (alrededor del 20% de la cosecha), debidas a la aplicación de diferentes tipos de cargas mecánicas (impacto, compresión, fricción, etc) a lo largo de los procesos, y mejorar el aspecto externo del fruto. Son bien conocidas las consecuencias del estado actual de los procesos de recolección y post-recolección a nivel de fruto. Cuando un fruto es sometido reiteradamente a cargas mecánicas exteriores, éste induce una respuesta fisiológica de autodefensa y, como consecuencia, puede desencadenar en los tejidos cambios tales como envejecimiento acelerado, degeneración y degradación. Todas estas reacciones redundan indiscutiblemente en las pérdidas de calidad arriba mencionadas. Una vez definido el problema, es necesario conocer en detalle las causas, los procesos en si mismos, ser capaces de detectar los puntos potencialmente peligrosos en los que los frutos están expuestos a impactos, compresiones y otras situaciones de estrés mecánico, cuestión que hasta hace poco ha sido mal resuelta. Hasta ahora la detección y caracterización de los puntos críticos se venía realizando de un modo intuitivo, basado en la simple observación de la ejecución de las operaciones por parte del responsable de cada proceso. Actualmente, con la ayuda de los frutos electrónicos simulados, es posible realizar la evaluación de cualquier procedimiento de forma rápida, precisa y objetiva. El desarrollo de los frutos electrónicos viene a cubrir una demanda del sector, tanto de los productores como de los agentes de comercialización, que necesitaban un método objetivo para la evaluación de los procedimientos y líneas de manipulación. Agencias de extensión agraria y algunas cooperativas han mostrado interés por estos dispositivos, empleándolos para la detección de puntos críticos en sus instalaciones. Incluso algunos agentes de grandes cadenas de supermercados han empezado a solicitar "certificados de calidad" de las líneas que manipulan los productos que ellos comercializan como garantía de su calidad.

Propiedades físicas de las variedades de tomate para la recolección mecánica

La utilización de máquinas para la recolección y manipulación de los productos agrícolas, aún los más frágiles (y perecederos) es un hecho. Sin embargo, se está lejos de llegar a una perfección en los resultados. Las relaciones entre los elementos que constituyen las máquinas, y los productos a manejar son desconocidas: Por un lado, falta el análisis de los elementos de los sistemas de recolección Por otro lado, se desconocen las causas intrínsecas o biológicas (histológicas, fisiológicas, bioquímicas) de la resistencia de los frutos a las distintas agresiones. Siendo el tomate uno de los frutos (de tipo "perecedero" o frágil) en los que más se ha avanzado en lo que respecta a recolección mecánica, sin haberse resuelto totalmente el problema, se cuenta con una experiencia, y un material que constituye una buena base de partida. La amplitud de la variabilidad de los materiales cultiva dos actualmente, es otro punto a favor del estudio del tomate. Aunque actualmente toda esta experiencia se basa en tomate para industria, la mayor parte de las conclusiones de estos estudios serán de inmediata aplicación en la mejora de variedades de tomate para fresca En éstos, ya existe actualmente la necesidad de técnicas de medida, así como de datos básicos, de las características físico-mecánicas de los frutos. En los estudios que se vienen realizando sobre cuajado de los frutos de tomate en invernadero son también necesarias estas técnicas y datos. El propósito de estos trabajos es conseguir técnicas para la diferenciación de variedades de tomate, respecto a las propiedades físicas que influyen en su resistencia a la manipulación mecánica. Dentro de ello, se intenta encontrar las relaciones causales entre las diversas características que puedan contribuir a un mayor - conocimiento de las propiedades mecánicas de estos frutos. La correlación de los ensayos de laboratorio con los sistemas de determinación de la calidad del trabajo - de recolección llevado a cabo por las máquinas contribuirá al desarrollo final de una óptima solución de mecanización. El estudio de los agentes o causas reales de los daños en los procesos de manipulación se inició en 1955, en que Halsey realizó un estudio preliminar sobre los - daños producidos a los frutos de tomate de mercado, in dicando ya que los daños eran producidos principalmente por impactos. A partir de los años 60 se han realizado numerosos estudios sobre daños producidos a los frutos de tomate en la recolección mecánica, dirigiéndose especialmente a los daños en el transporte a granel (aunque también a los producidos en la máquina cosechadora) (O'Brien - et al. 1965, id 1968, id. 1971, O'Brien 1974, Angle et al. 1974). Las características que confieren resistencia a los frutos en estos procesos suponen un complejo de factor e s , entre los que sobresalen la elasticidad del fruto - entero y la resistencia de la piel. En este proceso del transporte los estudios se basan en las características de vibración de los frutos, que dependen a su vez de los factores indicados (elasticidad, etc.).

Recolección Mecánica de Pimiento de Pimentón: Ensayos de Campo.

Para mecanizar de forma integral el cultivo de pimientos de pimentón es necesario el desarrollo de la mecanización de la recolección. En este trabajo se analizan los resultados obtenidos al ensayar el sistema de desprendimiento por peinado del fruto utilizando dos cosechadoras de judía verde sobre el cultivo sembrado a 20 cm y a 40 cm entre líneas, con resultados positivos referidos tanto al porcentaje de frutos recogidos como a la capacidad de trabajo de la máquina.

Cuantificación de la susceptibilidad a magulladura en frutos de hueso en términos probabilísticos.

El objetivo de este estudio es definir una magnitud que caracterice la susceptibilidad a daños mecánicos en frutos de hueso (melocotón cv. Caterina y Sudanell y albaricoque cv Búlida y Pepito). La aparición de magulladura en frutos de hueso está más relacionada con la deformación máxima alcanzada durante una carga de compresión (indicativo de su respuesta mecánica) que con el nivel de carga (fuerza) en sí mismo. La aparición de magulladura va acompañada de elevadas superficies de daño lo que hace que los frutos al magullarse queden directamente fuera del rango de aceptación comunitario (0.5cm2). Este hecho indica la necesidad de caracterizar la susceptibilidad a magulladura como probabilidad de daño. Este estudio indica la posibilidad de establecer la probabilidad de daño mediante la combinación de la respuesta mecánica del fruto (deformación máxima) y el nivel de carga (fuerza). Asimismo la respuesta mecánica se puede establecer instrumentalmente mediante la "dureza al tacto" y la tipología de comportamiento reológico.

Modelo para la clasificación no destructiva de fruta en grados de firmeza a partir de medidas ópticas y mecánicas. Aplicación a melocotones amarillos de carne dura.

Los melocotones objeto del estudio son amarillo de carne dura tipo pavía. Durante la campaña 1997, se analizaron muestras representativas de un número importante de variedades de melocotón que llegaban a cooperativas en la Región de Murcia y a un hipermercado de Madrid. Las variedades fueron Caterina, BabyGold, Sudanell, Vesubio y Miraflores. El número de frutos, melocotones amarillos de carne dura, fue de 224. Los ensayos que se realizaron fueron: 1-º Ensayo destructivo de estimación de firmeza por penetrometría Magness Taylor, realizado mediante punzón metálico de 8mm de diámetro, a una velocidad de 20 mm por minuto. 2-. Ensayos no destructivos: Impacto, realizado mediante el impactador del Laboratorio de Propiedades Físicas. 2.- Ensayo de deformación mediante el empleo de durómetro tipo Durofel-10. Posee un cilindro metálico que emerge 3 mm de superficie metálica y plana. Dicho cilindro se aplica perpendicularmente a la superficie del fruto; está conectado a un resorte que registra la fuerza correspondiente a la deformación máxima. Medida de reflectancia en el espectro visible desde 400 a 700 nm, mediante el uso del espectrofotómetro Minolta CM-508Í. Se consideraron las reflectancia correspondientes a 450 nm y a 680 nm, por ser las que mejor se correlacionan con la presencia de carotenoides y clorofila respectivamente. Ambos pigmentos están relacionados con el proceso de maduración, en el cual también se reblandecen los frutos.

Lidar uncertainty in complex terrain development of a bias correction methodology

El autor ha trabajado como parte del equipo de investigación en mediciones de viento en el Centro Nacional de Energías Renovables (CENER), España, en cooperación con la Universidad Politécnica de Madrid y la Universidad Técnica de Dinamarca. El presente reporte recapitula el trabajo de investigación realizado durante los últimos 4.5 años en el estudio de las fuentes de error de los sistemas de medición remota de viento, basados en la tecnología lidar, enfocado al error causado por los efectos del terreno complejo. Este trabajo corresponde a una tarea del paquete de trabajo dedicado al estudio de sistemas remotos de medición de viento, perteneciente al proyecto de intestigación europeo del 7mo programa marco WAUDIT. Adicionalmente, los datos de viento reales han sido obtenidos durante las campañas de medición en terreno llano y terreno complejo, pertenecientes al también proyecto de intestigación europeo del 7mo programa marco SAFEWIND. El principal objetivo de este trabajo de investigación es determinar los efectos del terreno complejo en el error de medición de la velocidad del viento obtenida con los sistemas de medición remota lidar. Con este conocimiento, es posible proponer una metodología de corrección del error de las mediciones del lidar. Esta metodología está basada en la estimación de las variaciones del campo de viento no uniforme dentro del volumen de medición del lidar. Las variaciones promedio del campo de viento son predichas a partir de los resultados de las simulaciones computacionales de viento RANS, realizadas para el parque experimental de Alaiz. La metodología de corrección es verificada con los resultados de las simulaciones RANS y validadas con las mediciones reales adquiridas en la campaña de medición en terreno complejo. Al inicio de este reporte, el marco teórico describiendo el principio de medición de la tecnología lidar utilizada, es presentado con el fin de familiarizar al lector con los principales conceptos a utilizar a lo largo de este trabajo. Posteriormente, el estado del arte es presentado en donde se describe los avances realizados en el desarrollo de la la tecnología lidar aplicados al sector de la energía eólica. En la parte experimental de este trabajo de investigación se ha estudiado los datos adquiridos durante las dos campañas de medición realizadas. Estas campañas has sido realizadas en terreno llano y complejo, con el fin de complementar los conocimiento adquiridos en casa una de ellas y poder comparar los efectos del terreno en las mediciones de viento realizadas con sistemas remotos lidar. La primer campaña experimental se desarrollo en terreno llano, en el parque de ensayos de aerogeneradores H0vs0re, propiedad de DTU Wind Energy (anteriormente Ris0). La segunda campaña experimental se llevó a cabo en el parque de ensayos de aerogeneradores Alaiz, propiedad de CENER. Exactamente los mismos dos equipos lidar fueron utilizados en estas campañas, haciendo de estos experimentos altamente relevantes en el contexto de evaluación del recurso eólico. Un equipo lidar está basado en tecnología de onda continua, mientras que el otro está basado en tecnología de onda pulsada. La velocidad del viento fue medida, además de con los equipos lidar, con anemómetros de cazoletas, veletas y anemómetros verticales, instalados en mástiles meteorológicos. Los sensores del mástil meteorológico son considerados como las mediciones de referencia en el presente estudio. En primera instancia, se han analizado los promedios diez minútales de las medidas de viento. El objetivo es identificar las principales fuentes de error en las mediciones de los equipos lidar causadas por diferentes condiciones atmosféricas y por el flujo no uniforme de viento causado por el terreno complejo. El error del lidar ha sido estudiado como función de varias propiedades estadísticas del viento, como lo son el ángulo vertical de inclinación, la intensidad de turbulencia, la velocidad vertical, la estabilidad atmosférica y las características del terreno. El propósito es usar este conocimiento con el fin de definir criterios de filtrado de datos. Seguidamente, se propone una metodología para corregir el error del lidar causado por el campo de viento no uniforme, producido por la presencia de terreno complejo. Esta metodología está basada en el análisis matemático inicial sobre el proceso de cálculo de la velocidad de viento por los equipos lidar de onda continua. La metodología de corrección propuesta hace uso de las variaciones de viento calculadas a partir de las simulaciones RANS realizadas para el parque experimental de Alaiz. Una ventaja importante que presenta esta metodología es que las propiedades el campo de viento real, presentes en las mediciones instantáneas del lidar de onda continua, puede dar paso a análisis adicionales como parte del trabajo a futuro. Dentro del marco del proyecto, el trabajo diario se realizó en las instalaciones de CENER, con supervisión cercana de la UPM, incluyendo una estancia de 1.5 meses en la universidad. Durante esta estancia, se definió el análisis matemático de las mediciones de viento realizadas por el equipo lidar de onda continua. Adicionalmente, los efectos del campo de viento no uniforme sobre el error de medición del lidar fueron analíticamente definidos, después de asumir algunas simplificaciones. Adicionalmente, durante la etapa inicial de este proyecto se desarrollo una importante trabajo de cooperación con DTU Wind Energy. Gracias a esto, el autor realizó una estancia de 1.5 meses en Dinamarca. Durante esta estancia, el autor realizó una visita a la campaña de medición en terreno llano con el fin de aprender los aspectos básicos del diseño de campañas de medidas experimentales, el estudio del terreno y los alrededores y familiarizarse con la instrumentación del mástil meteorológico, el sistema de adquisición y almacenamiento de datos, así como de el estudio y reporte del análisis de mediciones. ABSTRACT The present report summarizes the research work performed during last 4.5 years of investigation on the sources of lidar bias due to complex terrain. This work corresponds to one task of the remote sensing work package, belonging to the FP7 WAUDIT project. Furthermore, the field data from the wind velocity measurement campaigns of the FP7 SafeWind project have been used in this report. The main objective of this research work is to determine the terrain effects on the lidar bias in the measured wind velocity. With this knowledge, it is possible to propose a lidar bias correction methodology. This methodology is based on an estimation of the wind field variations within the lidar scan volume. The wind field variations are calculated from RANS simulations performed from the Alaiz test site. The methodology is validated against real scale measurements recorded during an eight month measurement campaign at the Alaiz test site. Firstly, the mathematical framework of the lidar sensing principle is introduced and an overview of the state of the art is presented. The experimental part includes the study of two different, but complementary experiments. The first experiment was a measurement campaign performed in flat terrain, at DTU Wind Energy H0vs0re test site, while the second experiment was performed in complex terrain at CENER Alaiz test site. Exactly the same two lidar devices, based on continuous wave and pulsed wave systems, have been used in the two consecutive measurement campaigns, making this a relevant experiment in the context of wind resource assessment. The wind velocity was sensed by the lidars and standard cup anemometry and wind vanes (installed on a met mast). The met mast sensors are considered as the reference wind velocity measurements. The first analysis of the experimental data is dedicated to identify the main sources of lidar bias present in the 10 minute average values. The purpose is to identify the bias magnitude introduced by different atmospheric conditions and by the non-uniform wind flow resultant of the terrain irregularities. The lidar bias as function of several statistical properties of the wind flow like the tilt angle, turbulence intensity, vertical velocity, atmospheric stability and the terrain characteristics have been studied. The aim of this exercise is to use this knowledge in order to define useful lidar bias data filters. Then, a methodology to correct the lidar bias caused by non-uniform wind flow is proposed, based on the initial mathematical analysis of the lidar measurements. The proposed lidar bias correction methodology has been developed focusing on the the continuous wave lidar system. In a last step, the proposed lidar bias correction methodology is validated with the data of the complex terrain measurement campaign. The methodology makes use of the wind field variations obtained from the RANS analysis. The results are presented and discussed. The advantage of this methodology is that the wind field properties at the Alaiz test site can be studied with more detail, based on the instantaneous measurements of the CW lidar. Within the project framework, the daily basis work has been done at CENER, with close guidance and support from the UPM, including an exchange period of 1.5 months. During this exchange period, the mathematical analysis of the lidar sensing of the wind velocity was defined. Furthermore, the effects of non-uniform wind fields on the lidar bias were analytically defined, after making some assumptions for the sake of simplification. Moreover, there has been an important cooperation with DTU Wind Energy, where a secondment period of 1.5 months has been done as well. During the secondment period at DTU Wind Energy, an important introductory learning has taken place. The learned aspects include the design of an experimental measurement campaign in flat terrain, the site assessment study of obstacles and terrain conditions, the data acquisition and processing, as well as the study and reporting of the measurement analysis.

Desarrollo, implementación y análisis de un modelo cinemático local de velocidades para Costa Rica basado en la integración y homogenización de datos de estaciones GNSS continuas

El primer procesamiento estricto realizado con el software científico Bernese y contemplando las más estrictas normas de cálculo recomendadas internacionalmente, permitió obtener un campo puntual de alta exactitud, basado en la integración y estandarización de los datos de una red GPS ubicada en Costa Rica. Este procesamiento contempló un total de 119 semanas de datos diarios, es decir unos 2,3 años, desde enero del año 2009 hasta abril del año 2011, para un total de 30 estaciones GPS, de las cuales 22 están ubicadas en el territorio nacional de Costa Rica y 8 internaciones pertenecientes a la red del Sistema Geocéntrico para las Américas (SIRGAS). Las denominadas soluciones semilibres generaron, semana a semana, una red GPS con una alta exactitud interna definida por medio de los vectores entre las estaciones y las coordenadas finales de la constelación satelital. La evaluación semanal dada por la repetibilidad de las soluciones brindó en promedio errores de 1,7 mm, 1,4 mm y 5,1 mm en las componentes [n e u], confirmando una alta consistencia en estas soluciones. Aunque las soluciones semilibres poseen una alta exactitud interna, las mismas no son utilizables para fines de análisis cinemático, pues carecen de un marco de referencia. En Latinoamérica, la densificación del Marco Internacional Terrestre de Referencia (ITRF), está representado por la red de estaciones de operación continua GNSS de SIRGAS, denominada como SIRGAS-CON. Por medio de las denominadas coordenadas semanales finales de las 8 estaciones consideradas como vínculo, se refirió cada una de las 119 soluciones al marco SIRGAS. La introducción del marco de referencia SIRGAS a las soluciones semilibres produce deformaciones en estas soluciones. Las deformaciones de las soluciones semilibres son producto de las cinemática de cada una de las placas en las que se ubican las estaciones de vínculo. Luego de efectuado el amarre semanal a las coordenadas SIRGAS, se hizo una estimación de los vectores de velocidad de cada una de las estaciones, incluyendo las de amarre, cuyos valores de velocidad se conocen con una alta exactitud. Para la determinación de las velocidades de las estaciones costarricenses, se programó una rutina en ambiente MatLab, basada en una ajuste por mínimos cuadrados. Los valores obtenidos en el marco de este proyecto en comparación con los valores oficiales, brindaron diferencias promedio del orden de los 0,06 cm/a, -0,08 cm/a y -0,10 cm/a respectivamente para las coordenadas [X Y Z]. De esta manera se logró determinar las coordenadas geocéntricas [X Y Z]T y sus variaciones temporales [vX vY vZ]T para el conjunto de 22 estaciones GPS de Costa Rica, dentro del datum IGS05, época de referencia 2010,5. Aunque se logró una alta exactitud en los vectores de coordenadas geocéntricas de las 22 estaciones, para algunas de las estaciones el cálculo de las velocidades no fue representativo debido al relativo corto tiempo (menos de un año) de archivos de datos. Bajo esta premisa, se excluyeron las ocho estaciones ubicadas al sur de país. Esto implicó hacer una estimación del campo local de velocidades con solamente veinte estaciones nacionales más tres estaciones en Panamá y una en Nicaragua. El algoritmo usado fue el denominado Colocación por Mínimos Cuadrados, el cual permite la estimación o interpolación de datos a partir de datos efectivamente conocidos, el cual fue programado mediante una rutina en ambiente MatLab. El campo resultante se estimó con una resolución de 30' X 30' y es altamente constante, con una velocidad resultante promedio de 2,58 cm/a en una dirección de 40,8° en dirección noreste. Este campo fue validado con base en los datos del modelo VEMOS2009, recomendado por SIRGAS. Las diferencias de velocidad promedio para las estaciones usadas como insumo para el cálculo del campo fueron del orden los +0,63 cm/a y +0,22 cm/a para los valores de velocidad en latitud y longitud, lo que supone una buena determinación de los valores de velocidad y de la estimación de la función de covarianza empírica, necesaria para la aplicación del método de colocación. Además, la grilla usada como base para la interpolación brindó diferencias del orden de -0,62 cm/a y -0,12 cm/a para latitud y longitud. Adicionalmente los resultados de este trabajo fueron usados como insumo para hacer una aproximación en la definición del límite del llamado Bloque de Panamá dentro del territorio nacional de Costa Rica. El cálculo de las componentes del Polo de Euler por medio de una rutina programa en ambiente MatLab y aplicado a diferentes combinaciones de puntos no brindó mayores aportes a la definición física de este límite. La estrategia lo que confirmó fue simplemente la diferencia en la dirección de todos los vectores velocidad y no permitió reveló revelar con mayor detalle una ubicación de esta zona dentro del territorio nacional de Costa Rica. ABSTRACT The first strict processing performed with the Bernese scientific software and contemplating the highest standards internationally recommended calculation, yielded a precise field of high accuracy, based on the integration and standardization of data from a GPS network located in Costa Rica. This processing watched a total of 119 weeks of daily data, is about 2.3 years from January 2009 to April 2011, for a total of 30 GPS stations, of which 22 are located in the country of Costa Rica and 8 hospitalizations within the network of Geocentric System for the Americas (SIRGAS). The semi-free solutions generated, every week a GPS network with high internal accuracy defined by vectors between stations and the final coordinates of the satellite constellation. The weekly evaluation given by repeatability of the solutions provided in average errors of 1.7 mm 1.4 mm and 5.1 mm in the components [n e u], confirming a high consistency in these solutions. Although semi-free solutions have a high internal accuracy, they are not used for purposes of kinematic analysis, because they lack a reference frame. In Latin America, the densification of the International Terrestrial Reference Frame (ITRF), is represented by a network of continuously operating GNSS stations SIRGAS, known as SIRGAS-CON. Through weekly final coordinates of the 8 stations considered as a link, described each of the solutions to the frame 119 SIRGAS. The introduction of the frame SIRGAS to semi-free solutions generates deformations. The deformations of the semi-free solutions are products of the kinematics of each of the plates in which link stations are located. After SIRGAS weekly link to SIRGAS frame, an estimate of the velocity vectors of each of the stations was done. The velocity vectors for each SIRGAS stations are known with high accuracy. For this calculation routine in MatLab environment, based on a least squares fit was scheduled. The values obtained compared to the official values, gave average differences of the order of 0.06 cm/yr, -0.08 cm/yr and -0.10 cm/yr respectively for the coordinates [XYZ]. Thus was possible to determine the geocentric coordinates [XYZ]T and its temporal variations [vX vY vZ]T for the set of 22 GPS stations of Costa Rica, within IGS05 datum, reference epoch 2010.5. The high accuracy vector for geocentric coordinates was obtained, however for some stations the velocity vectors was not representative because of the relatively short time (less than one year) of data files. Under this premise, the eight stations located in the south of the country were excluded. This involved an estimate of the local velocity field with only twenty national stations plus three stations in Panama and Nicaragua. The algorithm used was Least Squares Collocation, which allows the estimation and interpolation of data from known data effectively. The algorithm was programmed with MatLab. The resulting field was estimated with a resolution of 30' X 30' and is highly consistent with a resulting average speed of 2.58 cm/y in a direction of 40.8° to the northeast. This field was validated based on the model data VEMOS2009 recommended by SIRGAS. The differences in average velocity for the stations used as input for the calculation of the field were of the order of +0.63 cm/yr, +0.22 cm/yr for the velocity values in latitude and longitude, which is a good determination velocity values and estimating the empirical covariance function necessary for implementing the method of application. Furthermore, the grid used as the basis for interpolation provided differences of about -0.62 cm/yr, -0.12 cm/yr to latitude and longitude. Additionally, the results of this investigation were used as input to an approach in defining the boundary of Panama called block within the country of Costa Rica. The calculation of the components of the Euler pole through a routine program in MatLab and applied to different combinations of points gave no further contributions to the physical definition of this limit. The strategy was simply confirming the difference in the direction of all the velocity vectors and not allowed to reveal more detail revealed a location of this area within the country of Costa Rica.

Can seed production and restricted dispersal limit recruitment in Pinus pinaster Aiton from the Spanish Northern Plateau?

Natural regeneration faces increasing difficulties in dry forests from the Mediterranean basin, including for normally well-regenerating species such as maritime pine (Pinus pinaster Aiton). In this paper, we studied female fertility, seed dispersal and spread rates in P. pinaster from the Spanish Northern Plateau, where natural regeneration failure is a main concern for forest managers. For this purpose we periodically collected data from seed traps and trees located at two core locations across several years. We found significant variation in interannual cone production, with the best seed trees being the same across years. In addition, we found highly skewed distributions of female reproductive effort and large fertility differences across stands located few kilometres away. Annual seed dispersal kernels fitted lognormal or 2Dt models depending on the stand analysed, with median dispersal distances between 14 and 25 m. Kernels fitted for maximum dispersal periods showed an outstanding intraseasonal variation of median dispersal distances, from 10 to 54 m, in association to variable patterns of rainfall and maximum wind speed. The amount of seed produced appeared to be enough to guarantee the natural regeneration of the stands during the typical 20-year regeneration period. Colonisation simulations concluded that Mediterranean maritime pine has a notable dispersion capacity, which is strongly influenced by levels of fecundity and, especially, by the number and frequency of long-distance dispersal events. The latter play a key role in tree dispersion processes through enlarging the occupied area and fostering the invasion of abandoned crop land.

Puente mixto de canto variable sobre la A-30 en Murcia

El puente mixto sobre la A-30 se engloba dentro de las obras para la construcción de la autovía A-33, que conecta 3 corredores de gran importancia entre el centro peninsular y levante, concretamente la A-31 (Madrid-Alicante), la A-30 (Albacete-Cartagena) y la A-35 (Albacete-Valencia); y permite la conexión entre las autovías A-30 y A-33. Para esta estructura se adoptó una tipología del tablero de viga de canto constante con cartelas de refuerzo sobre apoyos, más propia de puentes de hormigón que de puentes mixtos, para lograr una mayor ligereza y transparencia de la estructura. El canto necesario para salvar el gran vano central derivado del esviaje en el cruce de vías y el estricto gálibo previsto, implican la necesidad de disminuir el peso visual del tablero para equilibrarlo con la luz aparente del vano principal, que es la percibida por los usuarios de la vía, y la reducida altura de pilas.

Determinación de la amplificación sísmica local en el Campus Sur UPM mediante la técnica de sismología de microtremores ReMi.

En el Campus Sur de la Universidad Politécnica de Madrid se ha llevado a cabo un proyecto para obtener una caracterización del subsuelo mediante ensayos ReMi, en colaboración con el departamento de Geofísica del Instituto Geográfico Nacional. La técnica ReMi (Refraction Microtremor) permite, mediante ensayos geofísicos realizados localmente sobre el terreno,obtener los parámetros físicos del mismo, que resultan de especial interés en el ámbito de la ingeniería civil. Esta técnica se caracteriza por englobarse dentro de la sísmica pasiva, muy empleada en prospección geofísica y basada en la obtención del modelo subyacente de distribución de velocidades de propagación de la onda S en función de la profundidad, con la ventaja de aprovechar el ruido sísmico ambiental como fuente de energía. Fue desarrollada en el Laboratorio Sismológico de Nevada (EEUU) por Louie (2001), con el objetivo de presentar una técnica innovadora en la obtención de las velocidades de propagación de manera experimental. Presenta ciertas ventajas, como la observación directa de la dispersión de ondas superficiales,que da un buen resultado de la velocidad de onda S, siendo un método no invasivo, de bajo coste y buena resolución, aplicable en entornos urbanos o sensibles en los que tanto otras técnicas sismológicas como otras variedades de prospección presentan dificultades. La velocidad de propagación de la onda S en los 30 primeros metros VS30, es ampliamente reconocida como un parámetro equivalente válido para caracterizar geotécnicamente el subsuelo y se halla matemáticamente relacionada con la velocidad de propagación de las ondas superficiales a observar mediante la técnica ReMi. Su observación permite el análisis espectral de los registros adquiridos, obteniéndose un modelo representado por la curva de dispersión de cada emplazamiento, de modo que mediante una inversión se obtiene el modelo de velocidad de propagación en función de la profundidad. A través de estos modelos, pueden obtenerse otros parámetros de interés sismológico. Estos resultados se representan sobre mapas isométricos para obtener una relación espacial de los mismos, particularmente conocido como zonación sísmica. De este análisis se extrae que la VS30 promedio del Campus no es baja en exceso, correspondiéndose a posteriori con los resultados de amplificación sísmica, período fundamental de resonancia del lugar y profundidad del sustrato rocoso. En última instancia se comprueba que los valores de amplificación sísmica máxima y el período al cual se produce posiblemente coincidan con los períodos fundamentales de resonancia de algunos edificios del Campus. ABSTRACT In South Campus at Polytechnic University of Madrid, a project has been carried out to obtain a proper subsoil description by applying ReMi tests, in collaboration with the Department of Geophysics of the National Geographic Institute. Through geophysical tests conducted locally, the ReMi (Refraction Microtremor) technique allows to establish the physical parameters of soil, which are of special interest in the field of civil engineering. This technique is part of passive seismic methods, often used in geophysical prospecting. It focuses in obtaining the underlying model of propagation velocity distribution of the shear wave according to depth and has the advantage of being able to use seismic ambient noise as a source of energy. It was developed in the Nevada Seismological Laboratory (USA) by Louie (2001) as an innovative technique for obtaining propagation velocities experimentally. It has several other advantages, including the direct observation of the dispersion of surface waves, which allows to reliably measure S wave velocity. This is a non-invasive, low cost and good resolution method, which can be applied in urban or sensitive environments where other prospection methods present difficulties. The propagation velocity of shear waves in the first 30 meters Vs30 is widely recognized as a valid equivalent parameter to geotechnically characterize the subsurface. It is mathematically related to surface wave's velocity of propagation, which are to observe using REMI technique. Spectral analysis of acquired data sets up a model represented by the dispersion curve at each site, so that, using an inversion process, propagation velocity model in relation to depth is obtained. Through this models, other seismologically interesting parameters can be obtained. These results are represented on isometric maps in order to obtain a spatial relationship between them, a process which is known as seismic zonation. This analysis infers that Vs30 at South Campus is not alarmingly low , corresponding with subsequent results of seismic amplification, fundamental period of resonance of soil and depth of bedrock. Ultimately, it's found that calculated values of soil's fundamental periods at which maximum seismic amplification occurs, may possibly match fundamental periods of some Campus buildings.

Estudio acústico de la sala 8 de los cines Kinépolis de Ciudad de la Imagen

Este proyecto consistira en la realization de un estudio aciistico sobre la sala 8 de los cines Kinepolis de Ciudad de la Imagen, que dispone de 408 butacas. Los cines Kinepolis es uno de los mayores complejos multisala de Europa. Cuenta con mas de 9.200 butacas en total distribuidas en 25 salas a las que se accede mediante dos grandes pasillos conectados por el hall. En 1998, ano de su apertura, el complejo recibio el Record Guinness a la sala cinematografica mas grande del mundo, que dispone de 996 butacas. El objetivo de este proyecto es conseguir caracterizar acusticamente una sala de cine a traves de la medicion de parametros acusticos de la sala y de un modelo virtual de la misma. Para llevar a cabo el proyecto, primero se van a realizar tanto una medicion geometrica como acustica de la sala mediante el sistema de medida DIRAC. Los resultados de estas mediciones nos serviran para construir y validar un modelo virtual de la sala real con el software de simulacion EASE. La medicion acustica se va a realizar con el sistema de medicion DIRAC. Este sistema nos dara information sobre una amplia variedad de parametros acusticos. En este proyecto no se va a trabajar con todos ellos, solo con los mas significativos. Estos se describen a continuacion en la introduccion teorica. La medicion geometrica nos va a servir para construir un modelo virtual que tenga las mismas dimensiones que la sala original. Esta medicion la realizaremos mediante un medidor laser y una cinta metrica. Una vez construido el modelo virtual, se procedera a su validacion. Este proceso se realiza ajustando el tiempo de reverberacion del modelo mediante la introduccion de distintos materiales acusticos en las superficies del mismo, de manera que, variando la absorcion de la sala, el tiempo de reverberacion promedio del modelo se asemeje lo mas posible al medido en la sala real. Este proceso tiene como objetivo comprobar que el modelo virtual tiene un comportamiento acustico similar al de la sala real. Es necesario validar adecuadamente el modelo para que las comparaciones y conclusiones sean fiables. Por ultimo, tras la simulacion acustica del modelo, se compararan los resultados simulados con los medidos en la sala. En este proceso se contrastaran algunos de los parametros que guardan relation con el tiempo de reverberacion. De esta manera se verificara si el tiempo de reverberacion es o no un parametro acustico fiable para la validacion de un modelo virtual de una sala de cine. Anteriormente se han realizado proyectos iguales de otras salas de diferente tamano de Kinepolis. El objetivo de realizar el mismo estudio en distintas salas, es comprobar si el tamano de la sala influye en la validacion de los modelos virtuales mediante el tiempo de reverberacion. ABSTRACT. This Project consists on the development of an acoustic research of the movie theater 8 of the Kinepolis complex in Ciudad de la Imagen, Madrid. This room has 408 spots. Kinepolis is one of the biggest multiplex complexes in Europe. It has 9,200 locations disposed in 25 rooms. There are two large corridors which give access to all of theaters. In the middle of the structure, there is the main hall that connects these corridors. In 1998, at the time when the complex was open, it was awarded with the Record Guinness for the biggest theater in the world, which has 996 locations. The target of this project is to successfully characterize the acoustics of a movie theater through reverberation time and a virtual model. In order to reach this goal, in the first place, we are going to perform both, an acoustic and a geometric measurement of the room using DIRAC measurement system. The results of these measures will allow us to build and validate a virtual model of the room, using the simulation software EASE. We are going to use the DIRAC system in order to accomplish the acoustic measure. This operation gives us a huge variety of acoustic parameters. Not all of these are going to be used for this research, only the most significant ones. These are described in the theoretical introduction. The geometric measure is essential to help us to build the virtual model, because the model has to be exactly equal as the real room. This measurement will be performed with an electronic distance meter and a measuring tape. Once the virtual model is finished, it will be proved. This validation process will be realized by adjusting the reverberation time in the model. We will change the walls materials, therefore, the overall absorption of the room will change. We want the model reverberation time resemble to the real one. This practice is going to ensure that the model acoustic performance is close to the real one. In addition, it has to be successfully validate of we want the future comparisons to be reliable. Finally, after the model virtual simulation, we will compare the simulated results with the measure in the room. In this process, we will compare not only the reverberation time, but others parameters that keep relation with the reverberation time. We will verify this way, if the reverberation time is or is not an appropriate acoustic parameter to validate a virtual model of a movie theater. There have been done others similar acoustic researches in different theaters with different sizes. The aim of performing similar researches in different rooms is to determine if the size of the room infers in the validation process.

Localización y mapeo visual monocular para robot móvil terrestre aplicado a la inspección ultrasónica aeronáutica

En el ámbito de la robótica de servicio, actualmente no existe una solución automatizada para la inspección ultrasónica de las partes de material compuesto de una aeronave durante las operaciones de mantenimiento que realiza la aerolínea. El desarrollo de las nuevas técnicas de acoplamiento acústico en seco en el método de inspección no destructiva por ultrasonidos, está conduciendo a posibilitar su uso con soluciones de menor coste respecto a las técnicas tradicionales, sin perder eficacia para detectar las deficiencias en las estructuras de material compuesto. Aunque existen aplicaciones de esta técnica con soluciones manuales, utilizadas en las fases de desarrollo y fabricación del material compuesto, o con soluciones por control remoto en sectores diferentes al aeronáutico para componentes metálicos, sin embargo, no existen con soluciones automatizadas para la inspección no destructiva por ultrasonidos de las zonas del avión fabricadas en material compuesto una vez la aeronave ha sido entregada a la aerolínea. El objetivo de este trabajo fin de master es evaluar el sistema de localización, basado en visión por ordenador, de una solución robotizada aplicada la inspección ultrasónica estructural de aeronaves en servicio por parte de las propias aerolíneas, utilizando las nuevas técnicas de acoplamiento acústico en seco, buscando la ventaja de reducir los tiempos y los costes en las operaciones de mantenimiento. Se propone como solución un robot móvil autónomo de pequeño tamaño, con control de posición global basado en técnicas de SLAM Visual Monocular, utilizando marcadores visuales externos para delimitar el área de inspección. Se ha supuesto la inspección de elementos de la aeronave cuya superficie se pueda considerar plana y horizontal, como son las superficies del estabilizador horizontal o del ala. Este supuesto es completamente aceptable en zonas acotadas de estos componentes, y de cara al objetivo del proyecto, no le resta generalidad. El robot móvil propuesto es un vehículo terrestre triciclo, de dos grados de libertad, con un sistema de visión monocular completo embarcado, incluyendo el hardware de procesamiento de visión y control de trayectoria. Las dos ruedas delanteras son motrices y la tercera rueda, loca, sirve únicamente de apoyo. La dirección, de tipo diferencial, permite al robot girar sin necesidad de desplazamiento, al conseguirse por diferencia de velocidad entre la rueda motriz derecha e izquierda. El sistema de inspección ultrasónica embarcado está compuesto por el hardware de procesamiento y registro de señal, y una rueda-sensor situada coaxialmente al eje de las ruedas motrices, y centrada entre estas, de modo que la medida de inspección se realiza en el centro de rotación del robot. El control visual propuesto se realiza mediante una estrategia “ver y mover” basada en posición, ejecutándose de forma secuencial la extracción de características visuales de la imagen, el cálculo de la localización global del robot mediante SLAM visual y el movimiento de éste mediante un algoritmo de control de posición-orientación respecto a referencias de paso de la trayectoria. La trayectoria se planifica a partir del mapa de marcas visuales que delimitan el área de inspección, proporcionado también por SLAM visual. Para validar la solución propuesta se ha optado por desarrollar un prototipo físico tanto del robot como de los marcadores visuales externos, a los que se someterán a una prueba de validación como alternativa a utilizar un entorno simulado por software, consistente en el reconocimiento del área de trabajo, planeamiento de la trayectoria y recorrido de la misma, de forma autónoma, registrando el posicionamiento real del robot móvil junto con el posicionamiento proporcionado por el sistema de localización SLAM. El motivo de optar por un prototipo es validar la solución ante efectos físicos que son muy complicados de modelar en un entorno de simulación, derivados de las limitaciones constructivas de los sistemas de visión, como distorsiones ópticas o saturación de los sensores, y de las limitaciones constructivas de la mecánica del robot móvil que afectan al modelo cinemático, como son el deslizamiento de las ruedas o la fluctuación de potencia de los motores eléctricos. El prototipo de marcador visual externo utilizado para la prueba de validación, ha sido un símbolo plano vertical, en blanco y negro, que consta de un borde negro rectangular dentro del cual se incluye una serie de marcas cuadradas de color negro, cuya disposición es diferente para cada marcador, lo que permite su identificación. El prototipo de robot móvil utilizado para la prueba de validación, ha sido denominado VINDUSTOR: “VIsual controlled Non-Destructive UltraSonic inspecTOR”. Su estructura mecánica ha sido desarrollada a partir de la plataforma comercial de robótica educacional LEGO© MINDSTORMS NXT 2.0, que incluye los dos servomotores utilizados para accionar las dos ruedas motrices, su controlador, las ruedas delanteras y la rueda loca trasera. La estructura mecánica ha sido especialmente diseñada con piezas LEGO© para embarcar un ordenador PC portátil de tamaño pequeño, utilizado para el procesamiento visual y el control de movimiento, y el sistema de captación visual compuesto por dos cámaras web de bajo coste, colocadas una en posición delantera y otra en posición trasera, con el fin de aumentar el ángulo de visión. El peso total del prototipo no alcanza los 2 Kg, siendo sus dimensiones máximas 20 cm de largo, 25 cm de ancho y 26 cm de alto. El prototipo de robot móvil dispone de un control de tipo visual. La estrategia de control es de tipo “ver y mover” dinámico, en la que se realiza un bucle externo, de forma secuencial, la extracción de características en la imagen, la estimación de la localización del robot y el cálculo del control, y en un bucle interno, el control de los servomotores. La estrategia de adquisición de imágenes está basada en un sistema monocular de cámaras embarcadas. La estrategia de interpretación de imágenes está basada en posición tridimensional, en la que los objetivos de control se definen en el espacio de trabajo y no en la imagen. La ley de control está basada en postura, relacionando la velocidad del robot con el error en la posición respecto a las referencias de paso de una trayectoria. La trayectoria es generada a partir del mapa de marcadores visuales externo. En todo momento, la localización del robot respecto a un sistema de referencia externo y el mapa de marcadores, es realizado mediante técnicas de SLAM visual. La auto-localización de un robot móvil dentro de un entorno desconocido a priori constituye uno de los desafíos más importantes en la robótica, habiéndose conseguido su solución en las últimas décadas, con una formulación como un problema numérico y con implementaciones en casos que van desde robots aéreos a robots en entornos cerrados, existiendo numerosos estudios y publicaciones al respecto. La primera técnica de localización y mapeo simultáneo SLAM fue desarrollada en 1989, más como un concepto que como un algoritmo único, ya que su objetivo es gestionar un mapa del entorno constituido por posiciones de puntos de interés, obtenidos únicamente a partir de los datos de localización recogidos por los sensores, y obtener la pose del robot respecto al entorno, en un proceso limitado por el ruido de los sensores, tanto en la detección del entorno como en la odometría del robot, empleándose técnicas probabilísticas aumentar la precisión en la estimación. Atendiendo al algoritmo probabilístico utilizado, las técnicas SLAM pueden clasificarse en las basadas en Filtros de Kalman, en Filtros de Partículas y en su combinación. Los Filtros de Kalman consideran distribuciones de probabilidad gaussiana tanto en las medidas de los sensores como en las medidas indirectas obtenidas a partir de ellos, de modo que utilizan un conjunto de ecuaciones para estimar el estado de un proceso, minimizando la media del error cuadrático, incluso cuando el modelo del sistema no se conoce con precisión, siendo el más utilizado el Filtro de Kalman Extendido a modelos nolineales. Los Filtros de Partículas consideran distribuciones de probabilidad en las medidas de los sensores sin modelo, representándose mediante un conjunto de muestras aleatorias o partículas, de modo que utilizan el método Montecarlo secuencial para estimar la pose del robot y el mapa a partir de ellas de forma iterativa, siendo el más utilizado el Rao-Backwell, que permite obtener un estimador optimizado mediante el criterio del error cuadrático medio. Entre las técnicas que combinan ambos tipos de filtros probabilísticos destaca el FastSLAM, un algoritmo que estima la localización del robot con un Filtro de Partículas y la posición de los puntos de interés mediante el Filtro de Kalman Extendido. Las técnicas SLAM puede utilizar cualquier tipo de sensor que proporcionen información de localización, como Laser, Sonar, Ultrasonidos o Visión. Los sensores basados en visión pueden obtener las medidas de distancia mediante técnicas de visión estereoscópica o mediante técnica de visión monocular. La utilización de sensores basados en visión tiene como ventajas, proporcionar información global a través de las imágenes, no sólo medida de distancia, sino también información adicional como texturas o patrones, y la asequibilidad del hardware frente a otros sensores. Sin embargo, su principal inconveniente es el alto coste computacional necesario para los complejos algoritmos de detección, descripción, correspondencia y reconstrucción tridimensional, requeridos para la obtención de la medida de distancia a los múltiples puntos de interés procesados. Los principales inconvenientes del SLAM son el alto coste computacional, cuando se utiliza un número elevado de características visuales, y su consistencia ante errores, derivados del ruido en los sensores, del modelado y del tratamiento de las distribuciones de probabilidad, que pueden producir el fallo del filtro. Dado que el SLAM basado en el Filtro de Kalman Extendido es una las técnicas más utilizadas, se ha seleccionado en primer lugar cómo solución para el sistema de localización del robot, realizando una implementación en la que las medidas de los sensores y el movimiento del robot son simulados por software, antes de materializarla en el prototipo. La simulación se ha realizado considerando una disposición de ocho marcadores visuales que en todo momento proporcionan ocho medidas de distancia con ruido aleatorio equivalente al error del sensor visual real, y un modelo cinemático del robot que considera deslizamiento de las ruedas mediante ruido aleatorio. Durante la simulación, los resultados han mostrado que la localización estimada por el algoritmo SLAM-EKF presenta tendencia a corregir la localización obtenida mediante la odometría, pero no en suficiente cuantía para dar un resultado aceptable, sin conseguir una convergencia a una solución suficientemente cercana a la localización simulada del robot y los marcadores. La conclusión obtenida tras la simulación ha sido que el algoritmo SLAMEKF proporciona inadecuada convergencia de precisión, debido a la alta incertidumbre en la odometría y a la alta incertidumbre en las medidas de posición de los marcadores proporcionadas por el sensor visual. Tras estos resultados, se ha buscado una solución alternativa. Partiendo de la idea subyacente en los Filtros de Partículas, se ha planteado sustituir las distribuciones de probabilidad gaussianas consideradas por el Filtro de Kalman Extendido, por distribuciones equi-probables que derivan en funciones binarias que representan intervalos de probabilidad no-nula. La aplicación de Filtro supone la superposición de todas las funciones de probabilidad no-nula disponibles, de modo que el resultado es el intervalo donde existe alguna probabilidad de la medida. Cómo la efectividad de este filtro aumenta con el número disponible de medidas, se ha propuesto obtener una medida de la localización del robot a partir de cada pareja de medidas disponibles de posición de los marcadores, haciendo uso de la Trilateración. SLAM mediante Trilateración Estadística (SLAM-ST) es como se ha denominado a esta solución propuesta en este trabajo fin de master. Al igual que con el algoritmo SLAM-EKF, ha sido realizada una implementación del algoritmo SLAM-ST en la que las medidas de los sensores y el movimiento del robot son simulados, antes de materializarla en el prototipo. La simulación se ha realizado en las mismas condiciones y con las mismas consideraciones, para comparar con los resultados obtenidos con el algoritmo SLAM-EKF. Durante la simulación, los resultados han mostrado que la localización estimada por el algoritmo SLAM-ST presenta mayor tendencia que el algoritmo SLAM-EKF a corregir la localización obtenida mediante la odometría, de modo que se alcanza una convergencia a una solución suficientemente cercana a la localización simulada del robot y los marcadores. Las conclusiones obtenidas tras la simulación han sido que, en condiciones de alta incertidumbre en la odometría y en la medida de posición de los marcadores respecto al robot, el algoritmo SLAM-ST proporciona mejores resultado que el algoritmo SLAM-EKF, y que la precisión conseguida sugiere la viabilidad de la implementación en el prototipo. La implementación del algoritmo SLAM-ST en el prototipo ha sido realizada en conjunción con la implementación del Sensor Visual Monocular, el Modelo de Odometría y el Control de Trayectoria. El Sensor Visual Monocular es el elemento del sistema SLAM encargado de proporcionar la posición con respecto al robot de los marcadores visuales externos, a partir de las imágenes obtenidas por las cámaras, mediante técnicas de procesamiento de imagen que permiten detectar e identificar los marcadores visuales que se hallen presentes en la imagen capturada, así como obtener las características visuales a partir de las cuales inferir la posición del marcador visual respecto a la cámara, mediante reconstrucción tridimensional monocular, basada en el conocimiento a-priori del tamaño real del mismo. Para tal fin, se ha utilizado el modelo matemático de cámara pin-hole, y se ha considerado las distorsiones de la cámara real mediante la calibración del sensor, en vez de utilizar la calibración de la imagen, tras comprobar el alto coste computacional que requiere la corrección de la imagen capturada, de modo que la corrección se realiza sobre las características visuales extraídas y no sobre la imagen completa. El Modelo de Odometría es el elemento del sistema SLAM encargado de proporcionar la estimación de movimiento incremental del robot en base a la información proporcionada por los sensores de odometría, típicamente los encoders de las ruedas. Por la tipología del robot utilizado en el prototipo, se ha utilizado un modelo cinemático de un robot tipo uniciclo y un modelo de odometría de un robot móvil de dos ruedas tipo diferencial, en el que la traslación y la rotación se determinan por la diferencia de velocidad de las ruedas motrices, considerando que no existe deslizamiento entre la rueda y el suelo. Sin embargo, el deslizamiento en las ruedas aparece como consecuencia de causas externas que se producen de manera inconstante durante el movimiento del robot que provocan insuficiente contacto de la rueda con el suelo por efectos dinámicos. Para mantener la validez del modelo de odometría en todas estas situaciones que producen deslizamiento, se ha considerado un modelo de incertidumbre basado en un ensayo representativo de las situaciones más habituales de deslizamiento. El Control de Trayectoria es el elemento encargado de proporcionar las órdenes de movimiento al robot móvil. El control implementado en el prototipo está basado en postura, utilizando como entrada la desviación en la posición y orientación respecto a una referencia de paso de la trayectoria. La localización del robot utilizada es siempre de la estimación proporcionada por el sistema SLAM y la trayectoria es planeada a partir del conocimiento del mapa de marcas visuales que limitan el espacio de trabajo, mapa proporcionado por el sistema SLAM. Las limitaciones del sensor visual embarcado en la velocidad de estabilización de la imagen capturada han conducido a que el control se haya implementado con la estrategia “mirar parado”, en la que la captación de imágenes se realiza en posición estática. Para evaluar el sistema de localización basado en visión del prototipo, se ha diseñado una prueba de validación que obtenga una medida cuantitativa de su comportamiento. La prueba consiste en la realización de forma completamente autónoma de la detección del espacio de trabajo, la planificación de una trayectoria de inspección que lo transite completamente, y la ejecución del recorrido de la misma, registrando simultáneamente la localización real del robot móvil junto con la localización proporcionada por el sistema SLAM Visual Monocular. Se han realizado varias ejecuciones de prueba de validación, siempre en las mismas condiciones iniciales de posición de marcadores visuales y localización del robot móvil, comprobando la repetitividad del ensayo. Los resultados presentados corresponden a la consideración de las medidas más pesimistas obtenidas tras el procesamiento del conjunto de medidas de todos los ensayos. Los resultados revelan que, considerando todo el espacio de trabajo, el error de posición, diferencia entre los valores de proporcionados por el sistema SLAM y los valores medidos de posición real, se encuentra en el entorno de la veintena de centímetros. Además, los valores de incertidumbre proporcionados por el sistema SLAM son, en todos los casos, superiores a este error. Estos resultados conducen a concluir que el sistema de localización basado en SLAM Visual, mediante un algoritmo de Trilateración Estadística, usando un sensor visual monocular y marcadores visuales externos, funciona, proporcionando la localización del robot móvil con respecto al sistema de referencia global inicial y un mapa de su situación de los marcadores visuales, con precisión limitada, pero con incertidumbre conservativa, al estar en todo momento el error real de localización por debajo del error estimado. Sin embargo, los resultados de precisión del sistema de localización no son suficientemente altos para cumplir con los requerimientos como solución robotizada aplicada a la inspección ultrasónica estructural de aeronaves en servicio. En este sentido, los resultados sugieren que la posible continuación de este trabajo en el futuro debe centrarse en la mejora de la precisión de localización del robot móvil, con líneas de trabajo encaminadas a mejorar el comportamiento dinámico del prototipo, en mejorar la precisión de las medidas de posición proporcionadas por el sensor visual y en optimizar el resultado del algoritmo SLAM. Algunas de estas líneas futuras podrían ser la utilización de plataformas robóticas de desarrollo alternativas, la exploración de técnicas de visión por computador complementarias, como la odometría visual, la visión omnidireccional, la visión estereoscópica o las técnicas de reconstrucción tridimensional densa a partir de captura monocular, y el análisis de algoritmos SLAM alternativos condicionado a disponer de una sustancial mejora de precisión en el modelo de odometría y en las medidas de posición de los marcadores.