Lithologic Control on the Scaling Properties of the First-Order Streams of Drainage Networks: A Monofractal Analysis

The movement of water through the landscape can be investigated at different scales. This study dealt with the interrelation between bedrock lithology and the geometry of the overlying drainage systems. Parameters of fractal analysis, such as fractal dimension and lacunarity, were used to measure and quantify this relationship. The interrelation between bedrock lithology and the geometry of the drainage systems has been widely studied in the last decades. The quantification of this linkage has not yet been clearly established. Several studies have selected river basins or regularly shaped areas as study units, assuming them to be lithologically homogeneous. This study considered irregular distributions of rock types, establishing areas of the soil map (1:25,000) with the same lithologic information as study units. The tectonic stability and the low climatic variability of the study region allowed effective investigation of the lithologic controls on the drainage networks developed on the plutonic rocks, the metamorphic rocks, and the sedimentary materials existing in the study area. To exclude the effect of multiple in- and outflows in the lithologically homogeneous units, we focused this study on the first-order streams of the drainage networks. The geometry of the hydrologic features was quantified through traditional metrics of fluvial geomorphology and scaling parameters of fractal analysis, such as the fractal dimension, the reference density, and the lacunarity. The results demonstrate the scale invariance of both the drainage networks and the set of first-order streams at the study scale and a relationship between scaling in the lithology and the drainage network.

Methodology of image analysis for study of the vertisols moisture content

The main problem to study vertical drainage from the moisture distribution, on a vertisol profile, is searching for suitable methods using these procedures. Our aim was to design a digital image processing methodology and its analysis to characterize the moisture content distribution of a vertisol profile. In this research, twelve soil pits were excavated on a ba re Mazic Pellic Vertisols ix of them in May 13/2011 and the rest in May 19 /2011 after a moderate rainfall event. Digital RGB images were taken from each vertisol pit using a Kodak? camera selecting a size of 1600x945 pixels. Each soil image was processed to homogenized brightness and then a spatial filter with several window sizes was applied to select the optimum one. The RGB image obtained were divided in each matrix color selecting the best thresholds for each one, maximum and minimum, to be applied and get a digital binary pattern. This one was analyzed by estimating two fractal scaling exponents box counting dimension D BC) and interface fractal dimension (D) In addition, three pre-fractal scaling coefficients were determinate at maximum resolution: total number of boxes intercepting the foreground pattern (A), fractal lacunarity (?1) and Shannon entropy S1). For all the images processed the spatial filter 9x9 was the optimum based on entropy, cluster and histogram criteria. Thresholds for each color were selected based on bimodal histograms.

Polymer interaction with macroscopic carbon nanotube fibres and fabrication of nanostructured composites

El ensamblado de nanotubos de carbono (CNT) como una fibra macroscópica en la cual están orientados preferentemente paralelos entre sí y al eje de la fibra, ha dado como resultado un nuevo tipo de fibra de altas prestaciones derivadas de la explotación eficiente de las propiedades axiales de los CNTs, y que tiene un gran número de aplicaciones potenciales. Fibras continuas de CNTs se produjeron en el Instituto IMDEA Materiales mediante el proceso de hilado directo durante la reacción de síntesis por deposición química de vapores. Uno de los objetivos de esta tesis es el estudio de la estructura de estas fibras mediante técnicas del estado del arte de difracción de rayos X de sincrotrón y la elaboración de un modelo estructural de dicho material. Mediciones texturales de adsorción de gases, análisis de micrografías de electrones y dispersión de rayos X de ángulo alto y bajo (WAXS/SAXS) indican que el material tiene una estructura mesoporosa con una distribución de tamaño de poros ancha derivada del amplio rango de separaciones entre manojos de CNTs, así como una superficie específica de 170m2/g. Los valores de dimensión fractal obtenidos mediante SAXS y análisis Barrett-Joyner-Halenda (BJH) de mediciones texturales coinciden en 2.4 y 2.5, respectivamente, resaltando el carácter de red de la estructura de dichas fibras. La estructura mesoporosa y tipo hilo de las fibra de CNT es accesible a la infiltración de moléculas externas (líquidos o polímeros). En este trabajo se estudian los cambios en la estructura multiescala de las fibras de CNTs al interactuar con líquidos y polímeros. Los efectos de la densificación en la estructura de fibras secas de CNT son estudiados mediante WAXS/SAXS. El tratamiento de densificación junta los manojos de la fibra (los poros disminuyen de tamaño), resultando en un incremento de la densidad de la fibra. Sin embargo, los dominios estructurales correspondientes a la transferencia de esfuerzo mecánica y carga eléctrica en los nanotubos no son afectados durante este proceso de densificación; como consecuencia no se produce un efecto sustancial en las propiedades mecánicas y eléctricas. Mediciones de SAXS and fibra de CNT antes y después de infiltración de líquidos confirman la penetración de una gran cantidad de líquidos que llena los poros internos de la fibra pero no se intercalan entre capas de nanotubos adyacentes. La infiltración de cadenas poliméricas de bajo peso molecular tiende a expandir los manojos en la fibra e incrementar el ángulo de apertura de los poros. Los resultados de SAXS indican que la estructura interna de la fibra en términos de la organización de las capas de tubos y su orientación no es afectada cuando las muestras consisten en fibras infiltradas con polímeros de alto peso molecular. La cristalización de varios polímeros semicristalinos es acelerada por la presencia de fibras de CNTs alineados y produce el crecimiento de una capa transcristalina normal a la superficie de la fibra. Esto es observado directamente mediante microscopía óptica polarizada, y detectado mediante calorimetría DSC. Las lamelas en la capa transcristalina tienen orientación de la cadena polimérica paralela a la fibra y por lo tanto a los nanotubos, de acuerdo con los patrones de WAXS. Esta orientación preferencial se sugiere como parte de la fuerza impulsora en la nucleación. La nucleación del dominio cristalino polimérico en la superficie de los CNT no es epitaxial. Ocurre sin haber correspondencia entre las estructuras cristalinas del polímero y los nanotubos. Estas observaciones contribuyen a la compresión del fenómeno de nucleación en CNTs y otros nanocarbonos, y sientan las bases para el desarrollo de composites poliméricos de gran escala basados en fibra larga de CNTs alineados. ABSTRACT The assembly of carbon nanotubes into a macroscopic fibre material where they are preferentially aligned parallel to each other and to the fibre axis has resulted in a new class of high-performance fibres, which efficiently exploits the axial properties of the building blocks and has numerous applications. Long, continuous CNT fibres were produced in IMDEA Materials Institute by direct fibre spinning from a chemical vapour deposition reaction. These fibres have a complex hierarchical structure covering multiple length scales. One objective of this thesis is to reveal this structure by means of state-of-the-art techniques such as synchrotron X-ray diffraction, and to build a model to link the fibre structural elements. Texture and gas absorption measurements, using electron microscopy, wide angle and small angle X-ray scattering (WAXS/SAXS), and pore size distribution analysis by Barrett-Joyner-Halenda (BJH), indicate that the material has a mesoporous structure with a wide pore size distribution arising from the range of fibre bundle separation, and a high surface area _170m2/g. Fractal dimension values of 2.4_2.5 obtained from the SAXS and BJH measurements highlight the network structure of the fibre. Mesoporous and yarn-like structure of CNT fibres make them accessible to the infiltration of foreign molecules (liquid or polymer). This work studies multiscale structural changes when CNT fibres interact with liquids and polymers. The effects of densification on the structure of dry CNT fibres were measured by WAXS/SAXS. The densification treatment brings the fibre bundles closer (pores become smaller), leading to an increase in fibre density. However, structural domains made of the load and charge carrying nanotubes are not affected; consequently, it has no substantial effect on mechanical and electrical properties. SAXS measurements on the CNT fibres before and after liquid infiltration imply that most liquids are able to fill the internal pores but not to intercalate between nanotubes. Successful infiltration of low molecular weight polymer chains tends to expand the fibre bundles and increases the pore-opening angle. SAXS results indicate that the inner structure of the fibre, in terms of the nanotube layer arrangement and the fibre alignment, are not largely affected when infiltrated with polymers of relatively high molecular weight. The crystallisation of a variety of semicrystalline polymers is accelerated by the presence of aligned fibres of CNTs and results in the growth of a transcrystalline layer perpendicular to the fibre surface. This can be observed directly under polarised optical microscope, and detected by the exothermic peaks during differential scanning calorimetry. The discussion on the driving forces for the enhanced nucleation points out the preferential chain orientation of polymer lamella with the chain axis parallel to the fibre and thus to the nanotubes, which is confirmed by two-dimensional WAXS patterns. A non-epitaxial polymer crystal growth habit at the CNT-polymer interface is proposed, which is independent of lattice matching between the polymer and nanotubes. These findings contribute to the discussion on polymer nucleation on CNTs and other nanocarbons, and their implication for the development of large polymer composites based on long and aligned fibres of CNTs.

Quantifying the relationship between drainage networks at hillslope scale and particle size distribution at pedon scale

Nowadays, translating information about hydrologic and soil properties and processes across scales has emerged as a major theme in soil science and hydrology, and suitable theories for upscaling or downscaling hydrologic and soil information are being looked forward. The recognition of low-order catchments as self-organized systems suggests the existence of a great amount of links at different scales between their elements. The objective of this work was to research in areas of homogeneous bedrock material, the relationship between the hierarchical structure of the drainage networks at hillslope scale and the heterogeneity of the particle-size distribution at pedon scale. One of the most innovative elements in this work is the choice of the parameters to quantify the organization level of the studied features. The fractal dimension has been selected to measure the hierarchical structure of the drainage networks, while the Balanced Entropy Index (BEI) has been the chosen parameter to quantify the heterogeneity of the particle-size distribution from textural data. These parameters have made it possible to establish quantifiable relationships between two features attached to different steps in the scale range. Results suggest that the bedrock lithology of the landscape constrains the architecture of the drainage networks developed on it and the particle soil distribution resulting in the fragmentation processes.

Quantifying the relationship between drainage networks at hillslope scale and particle size distribution at pedon scale

Nowadays, translating information about hydrologic and soil properties and processes across scales has emerged as a major theme in soil science and hydrology, and suitable theories for upscaling or downscaling hydrologic and soil information are being looked forward. The recognition of low-order catchments as self-organized systems suggests the existence of a great amount of links at different scales between their elements. The objective of this work was to research in areas of homogeneous bedrock material, the relationship between the hierarchical structure of the drainage networks at hillslope scale and the heterogeneity of the particle-size distribution at pedon scale. One of the most innovative elements in this work is the choice of the parameters to quantify the organization level of the studied features. The fractal dimension has been selected to measure the hierarchical structure of the drainage networks, while the Balanced Entropy Index (BEI) has been the chosen parameter to quantify the heterogeneity of the particle-size distribution from textural data. These parameters have made it possible to establish quantifiable relationships between two features attached to different steps in the scale range. Results suggest that the bedrock lithology of the landscape constrains the architecture of the drainage networks developed on it and the particle soil distribution resulting in the fragmentation processes.

Propriedades superficiais de filmes à base de gelatina

O objetivo do presente trabalho foi caracterizar as propriedades superficiais de filmes à base de gelatina. Para o qual foram elaborados filmes de: (i) Gelatina plastificado com glicerol (G) (gelatina: 5 g/100 g SFF; glicerol: 30 g/100 g de gelatina), (ii) Gelatina reforçado com montmorilhonita (G/MMT) (gelatina: 5 g/100 g SFF; glicerol: 30 g/100 g de gelatina; MMT: 5 g/100 g de gelatina) e Gelatina plastificado com citrato de acetiltributila (G/ATB) (gelatina: 5 g/100 g SFF; ATB: 50 g/100 de gelatina; lecitina de soja: 60 g/100 g de ATB; etanol: 20 g/100 g SFF). Os filmes foram produzidos mediante o uso de um aplicador automático de filmes \"Spreading\". Logo, os filmes foram submetidos a testes para determinação da espessura, umidade e propriedades óticas (brilho, cor e opacidade). Também foi caracterizada a microestrutura por microscopia eletrônica de varredura (MEV) e microscopia de força atômica (AFM); às imagens obtidas por MEV foi aplicado um analise de imagem mediante o programa Image J, para obter o valor da dimensão fractal (DF). Depois foram caracterizadas as propriedades superficiais de ângulo de contato (AC), molhabilidade ou coeficiente de espalhamento (Se), e energia livre superficial (ELS) mediante a medida do ângulo de contato pelo método da gota séssil (água: 5 µL e 1-Bromonaftaleno: 3 µL). Para o cálculo da ELS também foi aplicado o método de Owens-Wendt. Estas caracterizações foram feitas em ambos os lados do filme, lado ar e lado placa. A natureza do filme de G/ATB influenciou na umidade e as propriedades óticas, enquanto que os filmes de G e G/MMT apresentaram características similares. Em relação à microestrutura e rugosidade, o filme de G apresentou a superfície mais homogênea e lisa, contrario ao observado no filme de G/MMT, que apresentou a maior rugosidade seguida do filme de G/ATB. Foi observado que houve uma relação entre os valores de rugosidade e DF. De acordo com o valor do AC, os filmes apresentaram um caráter hidrofóbico, pois seus valores foram superiores a 65° (em ambos os lados dos filmes), na seguinte ordem: G/MMT > G > G/ATB; sendo que o Se seguiu a mesma tendência. Cabe mencionar também que não foi encontrada uma correlação significativa entre os valores de AC e rugosidade. Em função dos valores de AC, Se e ELS (especificamente a componente polar), o filme de G/ATB apresentou o menor caráter hidrofóbico, pois apresentou menores valores de AC e maiores valores de Se em comparação com os outros dois filmes. Os valores da componente polar da G/ATB foram os maiores, explicando de melhor maneira o caráter menos hidrofóbico deste filme. Pode-se concluir que os filmes à base de gelatina elaborados no presente trabalho têm propriedades hidrofóbicas (AC>65°), sendo a G/MMT o filme com melhores características hidrofóbicas.

An introduction to alysidal algebra (V): phenomenological components

Purpose – This paper aims to refer to a subjective approach to a type of complex system: human ecosystems, referred to as deontical impure systems (DIS) to capture a set of properties fundamental to the distinction between human and natural ecosystems. There are four main phenomenological components: directionality, intensity, connection energy and volume. The paper establishes thermodynamics of deontical systems based on the Law of Zipf and the temperature of information. Design/methodology/approach – Mathematical and logical development of human society structure. Findings – A fundamental question in this approach to DIS is the intensity or forces of a relation. Concepts are introduced as the system volume and propose a system thermodynamic theory. It hints at the possibility of adapting the fractal theory by introducing the fractal dimension of the system. Originality/value – This paper is a continuation of other previous papers and developing the theory of DIS.

Particel size distribution and microstructures of sediments from ODP Leg 205 sites (Table 1)

Ocean Drilling Program Legs 170 and 205 offshore Costa Rica provide structural observations which support a new model for the geometry and deformation response to the seismic cycle of the frontal sedimentary prism and decollement. The model is based on drillcore, thin section, and electron microscope observations. The decollement damage zone is a few tens of meters in width, it develops mainly within the frontal prism. A clear cm-thick fault core is observed 1.6 km from the trench. The lower boundary of the fault core is coincident with the lithological boundary between the frontal prism and the hemipelagic and pelagic sediment of the Cocos plate. Breccia clast distributions in the upper portion of the decollement damage zone were studied through fractal analysis. This analysis shows that the fractal dimension changes with brecciated fragment size, implying that deformation was not accommodated by self-similar fracturing. A higher fractal dimensionality correlates with smaller particle size, which indicates that different or additional grain-size reduction processes operated during shearing. The co-existence of two distinct fracturing processes is also confirmed by microscopic analysis in which extension fracturing in the upper part of the damage zone farthest from the fault core is frequent, while both extension and shear fracturing occur approaching the fault core. The coexistence of extensional and shear fracturing seems to be best explained by fluid pressure variations in response to variations of the compressional regime during the seismic cycle. During the co-seismic event, sub-horizontal compression and fluid pressure increase, triggering shear fracturing and fluid expulsion. Fractures migrate upward with fluids, contributing to the asymmetric shape of the decollement, while slip propagates. In the inter-seismic interval the frontal prismrelaxes and fluid pressure drops. The frontal prismgoes into diffuse extension during the intervalwhen plate convergence is accommodated by creep along the ductile fault core. The fault core is typically a barrier to deformation, which is explained by its weak, but impermeable, nature. The localized development of a damage zone beneath the fault core is characterized by shear fracturing that appears as the result of local strengthening of the detachment.

Analyzation of marine snow and fecal pellets collected in sediment trap CBi-2 off Cape Blanc, Mauritania

We analyzed size-specific dry mass, sinking velocity, and apparent diffusivity in field-sampled marine snow, laboratory-made aggregates formed by diatoms or coccolithophorids, and small and large zooplankton fecal pellets with naturally varying content of ballast materials. Apparent diffusivity was measured directly inside aggregates and large (millimeter-long) fecal pellets using microsensors. Large fecal pellets, collected in the coastal upwelling off Cape Blanc, Mauritania, showed the highest volume-specific dry mass and sinking velocities because of a high content of opal, carbonate, and lithogenic material (mostly Saharan dust), which together comprised ~80% of the dry mass. The average solid matter density within these large fecal pellets was 1.7 g cm**-3, whereas their excess density was 0.25 ± 0.07 g cm**-3. Volume-specific dry mass of all sources of aggregates and fecal pellets ranged from 3.8 to 960 µg mm**-3, and average sinking velocities varied between 51 and 732 m d**-1. Porosity was >0.43 and >0.96 within fecal pellets and phytoplankton-derived aggregates, respectively. Averaged values of apparent diffusivity of gases within large fecal pellets and aggregates were 0.74 and 0.95 times that of the free diffusion coefficient in sea water, respectively. Ballast increases sinking velocity and, thus, also potential O2 fluxes to sedimenting aggregates and fecal pellets. Hence, ballast minerals limit the residence time of aggregates in the water column by increasing sinking velocity, but apparent diffusivity and potential oxygen supply within aggregates are high, whereby a large fraction of labile organic carbon can be respired during sedimentation.

The comparative characterization of structural heterogeneity of mesoporous activated carbon fibers (ACFs)

The comparative analysis of the most widely used methods of mesoporosity characterization of two activated carbon fibers is presented. Not only the older methods are used, i.e. Barrett, Joyner and Halenda (BJH), Dubinin (the so-called first variant-D-1ST and the so-called second variant-D-2ND), Dollimore and Heal (DH), and Pierce (P) but the recently developed ones, i.e. the method of Nguyen and Do (ND) and that developed by Do (Do) are also applied. Additionally, the method of the characterization of fractality is put to use (fractal analog of FHH isotherm). The results are compared and discussed. (C) 2002 Elsevier Science B.V. All fights reserved.

Impacts of morphological, physical and chemical properties of sludge flocs on dewaterability of activated sludge

This study examined how the floc characteristics affect dewaterability of activated sludge. The floc properties were characterized by morphological parameters (floc size distribution, fractal dimension and filament index), physical properties (flocculating ability, surface charge, relative hydrophobicity and viscosity), and chemical constituents in sludge and extracted extracellular polymeric substances (EPS), including the polymeric compounds protein, humic substances, carbohydrates and the ions Ca2+, Mg2+, Fe3+ and Al3+. The dewaterability was defined in terms of the bound water content and capillary suction time (CST). The bound water and CST corresponded to a similar indication with respect to dewaterability of activated sludge. The floc physical parameters were the most important factors which effect significantly on the water binding ability of the sludge flocs. The morphological characteristics had relatively weak impact on the dewaterability. The polymeric components protein and carbohydrate had a significant contribution to enhance the water binding ability of the sludge flocs. The effect of humic substances in the sludge on the dewaterability was, however, insignificant. The CST had good statistical correlations with the polymeric constituents measured in both sludge and the extracted EPS, and the bound water was only correlated well with the individual polymers measured in the sludge. High concentration of Ca2+, Mg2+, Fe3+ and Al3+ had significant improvement for dewaterability. (C) 2003 Elsevier B.V. All rights reserved.

A fuzzy approach to texture segmentation

The texture segmentation techniques are diversified by the existence of several approaches. In this paper, we propose fuzzy features for the segmentation of texture image. For this purpose, a membership function is constructed to represent the effect of the neighboring pixels on the current pixel in a window. Using these membership function values, we find a feature by weighted average method for the current pixel. This is repeated for all pixels in the window treating each time one pixel as the current pixel. Using these fuzzy based features, we derive three descriptors such as maximum, entropy, and energy for each window. To segment the texture image, the modified mountain clustering that is unsupervised and fuzzy c-means clustering have been used. The performance of the proposed features is compared with that of fractal features.

Applications of fractals to image data compression

Digital image processing is exploited in many diverse applications but the size of digital images places excessive demands on current storage and transmission technology. Image data compression is required to permit further use of digital image processing. Conventional image compression techniques based on statistical analysis have reached a saturation level so it is necessary to explore more radical methods. This thesis is concerned with novel methods, based on the use of fractals, for achieving significant compression of image data within reasonable processing time without introducing excessive distortion. Images are modelled as fractal data and this model is exploited directly by compression schemes. The validity of this is demonstrated by showing that the fractal complexity measure of fractal dimension is an excellent predictor of image compressibility. A method of fractal waveform coding is developed which has low computational demands and performs better than conventional waveform coding methods such as PCM and DPCM. Fractal techniques based on the use of space-filling curves are developed as a mechanism for hierarchical application of conventional techniques. Two particular applications are highlighted: the re-ordering of data during image scanning and the mapping of multi-dimensional data to one dimension. It is shown that there are many possible space-filling curves which may be used to scan images and that selection of an optimum curve leads to significantly improved data compression. The multi-dimensional mapping property of space-filling curves is used to speed up substantially the lookup process in vector quantisation. Iterated function systems are compared with vector quantisers and the computational complexity or iterated function system encoding is also reduced by using the efficient matching algcnithms identified for vector quantisers.

Services for Satellite Data Processing

Data processing services for Meteosat geostationary satellite are presented. Implemented services correspond to the different levels of remote-sensing data processing, including noise reduction at preprocessing level, cloud mask extraction at low-level and fractal dimension estimation at high-level. Cloud mask obtained as a result of Markovian segmentation of infrared data. To overcome high computation complexity of Markovian segmentation parallel algorithm is developed. Fractal dimension of Meteosat data estimated using fractional Brownian motion models.

Radiotherapy Treatment Assessment using DCE-MRI

Abstract

The goal of modern radiotherapy is to precisely deliver a prescribed radiation dose to delineated target volumes that contain a significant amount of tumor cells while sparing the surrounding healthy tissues/organs. Precise delineation of treatment and avoidance volumes is the key for the precision radiation therapy. In recent years, considerable clinical and research efforts have been devoted to integrate MRI into radiotherapy workflow motivated by the superior soft tissue contrast and functional imaging possibility. Dynamic contrast-enhanced MRI (DCE-MRI) is a noninvasive technique that measures properties of tissue microvasculature. Its sensitivity to radiation-induced vascular pharmacokinetic (PK) changes has been preliminary demonstrated. In spite of its great potential, two major challenges have limited DCE-MRI’s clinical application in radiotherapy assessment: the technical limitations of accurate DCE-MRI imaging implementation and the need of novel DCE-MRI data analysis methods for richer functional heterogeneity information.

This study aims at improving current DCE-MRI techniques and developing new DCE-MRI analysis methods for particular radiotherapy assessment. Thus, the study is naturally divided into two parts. The first part focuses on DCE-MRI temporal resolution as one of the key DCE-MRI technical factors, and some improvements regarding DCE-MRI temporal resolution are proposed; the second part explores the potential value of image heterogeneity analysis and multiple PK model combination for therapeutic response assessment, and several novel DCE-MRI data analysis methods are developed.

I. Improvement of DCE-MRI temporal resolution. First, the feasibility of improving DCE-MRI temporal resolution via image undersampling was studied. Specifically, a novel MR image iterative reconstruction algorithm was studied for DCE-MRI reconstruction. This algorithm was built on the recently developed compress sensing (CS) theory. By utilizing a limited k-space acquisition with shorter imaging time, images can be reconstructed in an iterative fashion under the regularization of a newly proposed total generalized variation (TGV) penalty term. In the retrospective study of brain radiosurgery patient DCE-MRI scans under IRB-approval, the clinically obtained image data was selected as reference data, and the simulated accelerated k-space acquisition was generated via undersampling the reference image full k-space with designed sampling grids. Two undersampling strategies were proposed: 1) a radial multi-ray grid with a special angular distribution was adopted to sample each slice of the full k-space; 2) a Cartesian random sampling grid series with spatiotemporal constraints from adjacent frames was adopted to sample the dynamic k-space series at a slice location. Two sets of PK parameters’ maps were generated from the undersampled data and from the fully-sampled data, respectively. Multiple quantitative measurements and statistical studies were performed to evaluate the accuracy of PK maps generated from the undersampled data in reference to the PK maps generated from the fully-sampled data. Results showed that at a simulated acceleration factor of four, PK maps could be faithfully calculated from the DCE images that were reconstructed using undersampled data, and no statistically significant differences were found between the regional PK mean values from undersampled and fully-sampled data sets. DCE-MRI acceleration using the investigated image reconstruction method has been suggested as feasible and promising.

Second, for high temporal resolution DCE-MRI, a new PK model fitting method was developed to solve PK parameters for better calculation accuracy and efficiency. This method is based on a derivative-based deformation of the commonly used Tofts PK model, which is presented as an integrative expression. This method also includes an advanced Kolmogorov-Zurbenko (KZ) filter to remove the potential noise effect in data and solve the PK parameter as a linear problem in matrix format. In the computer simulation study, PK parameters representing typical intracranial values were selected as references to simulated DCE-MRI data for different temporal resolution and different data noise level. Results showed that at both high temporal resolutions (<1s) and clinically feasible temporal resolution (~5s), this new method was able to calculate PK parameters more accurate than the current calculation methods at clinically relevant noise levels; at high temporal resolutions, the calculation efficiency of this new method was superior to current methods in an order of 102. In a retrospective of clinical brain DCE-MRI scans, the PK maps derived from the proposed method were comparable with the results from current methods. Based on these results, it can be concluded that this new method can be used for accurate and efficient PK model fitting for high temporal resolution DCE-MRI.

II. Development of DCE-MRI analysis methods for therapeutic response assessment. This part aims at methodology developments in two approaches. The first one is to develop model-free analysis method for DCE-MRI functional heterogeneity evaluation. This approach is inspired by the rationale that radiotherapy-induced functional change could be heterogeneous across the treatment area. The first effort was spent on a translational investigation of classic fractal dimension theory for DCE-MRI therapeutic response assessment. In a small-animal anti-angiogenesis drug therapy experiment, the randomly assigned treatment/control groups received multiple fraction treatments with one pre-treatment and multiple post-treatment high spatiotemporal DCE-MRI scans. In the post-treatment scan two weeks after the start, the investigated Rényi dimensions of the classic PK rate constant map demonstrated significant differences between the treatment and the control groups; when Rényi dimensions were adopted for treatment/control group classification, the achieved accuracy was higher than the accuracy from using conventional PK parameter statistics. Following this pilot work, two novel texture analysis methods were proposed. First, a new technique called Gray Level Local Power Matrix (GLLPM) was developed. It intends to solve the lack of temporal information and poor calculation efficiency of the commonly used Gray Level Co-Occurrence Matrix (GLCOM) techniques. In the same small animal experiment, the dynamic curves of Haralick texture features derived from the GLLPM had an overall better performance than the corresponding curves derived from current GLCOM techniques in treatment/control separation and classification. The second developed method is dynamic Fractal Signature Dissimilarity (FSD) analysis. Inspired by the classic fractal dimension theory, this method measures the dynamics of tumor heterogeneity during the contrast agent uptake in a quantitative fashion on DCE images. In the small animal experiment mentioned before, the selected parameters from dynamic FSD analysis showed significant differences between treatment/control groups as early as after 1 treatment fraction; in contrast, metrics from conventional PK analysis showed significant differences only after 3 treatment fractions. When using dynamic FSD parameters, the treatment/control group classification after 1st treatment fraction was improved than using conventional PK statistics. These results suggest the promising application of this novel method for capturing early therapeutic response.

The second approach of developing novel DCE-MRI methods is to combine PK information from multiple PK models. Currently, the classic Tofts model or its alternative version has been widely adopted for DCE-MRI analysis as a gold-standard approach for therapeutic response assessment. Previously, a shutter-speed (SS) model was proposed to incorporate transcytolemmal water exchange effect into contrast agent concentration quantification. In spite of richer biological assumption, its application in therapeutic response assessment is limited. It might be intriguing to combine the information from the SS model and from the classic Tofts model to explore potential new biological information for treatment assessment. The feasibility of this idea was investigated in the same small animal experiment. The SS model was compared against the Tofts model for therapeutic response assessment using PK parameter regional mean value comparison. Based on the modeled transcytolemmal water exchange rate, a biological subvolume was proposed and was automatically identified using histogram analysis. Within the biological subvolume, the PK rate constant derived from the SS model were proved to be superior to the one from Tofts model in treatment/control separation and classification. Furthermore, novel biomarkers were designed to integrate PK rate constants from these two models. When being evaluated in the biological subvolume, this biomarker was able to reflect significant treatment/control difference in both post-treatment evaluation. These results confirm the potential value of SS model as well as its combination with Tofts model for therapeutic response assessment.

In summary, this study addressed two problems of DCE-MRI application in radiotherapy assessment. In the first part, a method of accelerating DCE-MRI acquisition for better temporal resolution was investigated, and a novel PK model fitting algorithm was proposed for high temporal resolution DCE-MRI. In the second part, two model-free texture analysis methods and a multiple-model analysis method were developed for DCE-MRI therapeutic response assessment. The presented works could benefit the future DCE-MRI routine clinical application in radiotherapy assessment.