Selection of Grids for Weed Mapping

This study aimed to assess the degree of similarity presented by thematic maps generated by different sampling grids of weed plants in a commercial agricultural area of 7.95 hectares. Monocotyledons and dicotyledons were counted on the 2012/2013 and 2013/2014 harvests, before soybean planting, in the fallow period after wheat harvest, in both years. A regular grid of 10 x 10 m was produced to sample the invasive plants, used as reference, and the counting was done in 1 m² of each sample point, totaling 795 samples in each year, compared to regular grids of 30 and 50 m, generated from the data exclusion of the standard grid. Twenty-two composite soil samples were taken at a depth of 0-20 cm to correlate soil properties with weeds occurrence. For the generation of the thematic maps, the Inverse Distance Weighting (IDW) for interpolation was used; when comparing the maps generated from each grid with the reference map, the kappa coefficient was used to assess the loss of quality of the maps as the number of sample points was reduced. It was observed that the map quality loss was lower in 2013 compared to 2012 when the sampling density of the points was reduced. The 30 x 30 m grids have satisfactorily described the infestation data of the dicotyledons and the 50 x 50 m grids have adequately described the monocotyledon weeds infestation, compared to the standard 10 x 10 m grids.

Self-adaptive grids for noise mapping refinement

Noise maps are usually represented as contour or isolines maps describing the sound levels in a region. Using this kind of representation the user can easily find the noise level assigned to every location in the map. But the acoustic calculations behind the map are not performed for every single location on it; they are only performed in a grid of receivers. The results in this calculation grid are interpolated to draw the isolines or contours. Therefore, the resolution of the calculation grid and the way it was created (rectangular, triangulated, random…) have an effect on the resulting map. In this paper we describe a smart iterative procedure to optimize the quality of the map at a really low additional computational cost, using self-adaptive grids for the acoustic calculations. These self-adaptive grids add new receivers to the sampling grid in those locations where they are expected to be more useful, so that the performance at the output of the interpolator is enhanced. Self-adaptive sampling grids can be used for minimizing the overall error of the map (improving its quality), or for reducing calculation times, and can be also applied selectively to target areas or contour lines. This can be done by the user customizing the maximum number of iterations, the number of new receivers for each iteration, the target isolines, the target quality…

Wheat yield and physical properties of a brown latosol under no-tillage in south-central Paraná

Soil management influences the chemical and physical properties of soil. Chemical conditions have been thoroughly studied, while the role of soil physical conditions regarding crop yield has been neglected. This study aimed to analyze the wheat yield and its relationship with physical properties of an Oxisol under no-tillage (NT). The study was carried out between 2010 and 2011, in Reserva do Iguaçu, State of Paraná, Brazil, on the Campo Bonito farm, after 25 years of NT management. Based on harvest maps of barley (2006), wheat (2007) and maize (2009) of a plot (150 ha), zones with higher and lower yield potential (Z1 and Z2, respectively) were identified. Sampling grids with 16 units (50 x 50 m) and three sampling points per unit were established. The wheat grain yield (GY) and water infiltration capacity (WIC) were evaluated in 2010. Soil samples with disturbed and undisturbed structure were collected from the 0.00-0.10 and 0.10-0.20 m layers. The former were used to determine soil organic carbon (Corg) levels and the latter to determine soil bulk density (BD), total porosity (TP), macroporosity (Mac), and microporosity (Mic). Soil penetration resistance (PR) and water content (SWC) were also evaluated. The wheat GY of the whole plot was close to the regional average and the yield between the zones differed significantly, i.e. 22 % higher in Z1 than in Z2. No significant variation in Mic was observed between zones, but Z1 had higher Corg levels, SWC, TP and Mac and lower BD than Z2 in both soil layers, as well as a lower PR than Z2 in the 0.00-0.10 m layer. Therefore, soil physical conditions were more restrictive in Z2, in agreement with wheat yield and zone yield potential defined a priori, based on the harvest maps. Soil WIC in Z1 was significantly higher (30 %) than in Z2, in agreement with the results of TP and Mac which were also higher in Z1 in both soil layers. The correlation analysis of data of the two layers showed a positive relationship between wheat GY and the soil properties TP, SWC and WIC.

Spatial estimation of foliar phosphorus in different species of the genus Coffea based on soil properties

Information underlying analyses of coffee fertilization systems should consider both the soil and the nutritional status of plants. This study investigated the spatial relationship between phosphorus (P) levels in coffee plant tissues and soil chemical and physical properties. The study was performed using two arabica and one canephora coffee variety. Sampling grids were established in the areas, and the points georeferenced. The assessed properties of the soil were levels of available phosphorus (P-Mehlich), remaining phosphorus (P-rem) and particle size, and of the plant tissue, phosphorus levels (foliar P). The data were subjected to descriptive statistical analysis, correlation analysis, cluster analysis, and probability tests. Geostatistical and trend analyses were only performed for pairs of variables with significant linear correlation. The spatial variability for foliar P content was high for the variety Catuai and medium for the other evaluated plants. Unlike P-Mehlich, the variability in P-rem of the soil indicated the nutritional status of this nutrient in the plant.

Soybean yield maps using regular and optimized sample with different configurations by simulated annealing

ABSTRACT This study aimed to compare thematic maps of soybean yield for different sampling grids, using geostatistical methods (semivariance function and kriging). The analysis was performed with soybean yield data in t ha-1 in a commercial area with regular grids with distances between points of 25x25 m, 50x50 m, 75x75 m, 100x100 m, with 549, 188, 66 and 44 sampling points respectively; and data obtained by yield monitors. Optimized sampling schemes were also generated with the algorithm called Simulated Annealing, using maximization of the overall accuracy measure as a criterion for optimization. The results showed that sample size and sample density influenced the description of the spatial distribution of soybean yield. When the sample size was increased, there was an increased efficiency of thematic maps used to describe the spatial variability of soybean yield (higher values of accuracy indices and lower values for the sum of squared estimation error). In addition, more accurate maps were obtained, especially considering the optimized sample configurations with 188 and 549 sample points.

Desenvolvimento de um perfilômetro laser para determinação da área e volume de contato entre o pneu e o solo: Cristiano Alberto de Paula. -

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Determinação da densidade amostral ótima para a geração de mapas temáticos na cafeicultura de precisão

Among the crops commercially exploited in Brazil, the coffee has a great economic importance, especially in the states of Minas Gerais and Espirito Santo. In the search for higher yield and lower environmental impact, farmers and researchers seek to develop new technologies that result in greater efficiency in various production processes of the coffee. For this, the adoption of precision agriculture in the management of operations in coffee crops, called precision coffee, has shown results that justify its use, by identifying the spatial variability of several variables, allowing its localized management and in the proper intensity. Unlike conventional management that is based on the average of observations in an area, precision agriculture uses a more detailed sampling, based on a sampling grid, which allows to represent in greater detail the reality of farming. Many previous studies have identified the spatial variability of the production of coffee system variables, but without worrying about the quality of information obtained due to the sampling grid used as precision and accuracy. Given the above, the objective of this study was to evaluate the quality of four different sampling grids for different variables and three times, in order to identify the most appropriate grid for use in precision coffee. Also aimed to compare the results between the precision coffee and conventional, according to reference values.

Variabilidade espacial de atributos de solos coesos do leste maranhense

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Variabilidade espacial de atributos dendrométricas do eucalipto e de um neossolo quartzarênico a partir da aplicação de calcário ou lama de cal + oxyfertil

Pós-graduação em Agronomia - FEIS

Do predators aggregate in response to pest density in agroecosystems? Assessing within-field spatial patterns

1. The spatial heterogeneity of predator populations is an important component of ecological theories pertaining to predator-prey dynamics. Most studies within agricultural fields show spatial correlation (positive or negative) between mean predator numbers and prey abundance across a whole field over time but generally ignore the within-field spatial dimension. We used explicit spatial mapping to determine if generalist predators aggregated within a soybean field, the size of these aggregations and if predator aggregation was associated with pest aggregation, plant damage and predation rate. 2. The study was conducted at Gatton in the Lockyer Valley, 90 km west of Brisbane, Australia. Intensive sampling grids were used to investigate within-field spatial patterns. The first row of each grid was located in a lucerne field (10 m from interface) and the remaining rows were in an adjacent soybean field. At each point on the grid the abundance of foliage-dwelling and ground-dwelling pests and predators was measured, predation rates [using sentinel Helicoverpa armigera (Hubner) egg cards] and plant damage were estimated. Eight grids were sampled across two summer cropping seasons (2000/01, 2001/02). 3. Predators exhibited strong spatial patterning with regions of high and low abundance and activity within what are considered to be uniform soybean fields. Ground-dwelling and foliage-dwelling predators were often aggregated in patches approximately 40 m across. 4. Lycosidae (wolf spiders) displayed aggregation and were consistently more abundant within the lucerne, with a decreasing trap catch with distance from the lucrene/soybean interface. This trend was consistent between subsequent grids in a single field and between fields. 5. The large amount of spatial variability in within-field arthropod abundance (pests and predators) and activity (egg predation and plant damage) indicates that whole field averages were misleading. This result has serious implications for sampling of arthropod abundance and pest management decision-making based on scouting data. 6. There was a great deal of temporal change in the significant spatial patterns observed within a field at each sampling time point during a single season. Predator and pest aggregations observed in these fields were generally not stable for the entire season. 7. Predator aggregation did not correlate consistently with pest aggregation, plant damage or predation rate. Spatial patterns in predator abundance were not associated consistently with any single parameter measured. The most consistent positive association was between foliage-dwelling predators and pests (significant in four of seven grids). Inferring associations between predators and prey based on an intensive one-off sampling grid is difficult, due to the temporal variability in the abundance of each group. 8. Synthesis and applications. This study demonstrated that generalist predator populations are rarely distributed randomly and field edges and adjacent crops can have an influence on within-field predator abundance. This must be considered when estimating arthropod (pest and predator) abundance from a set of samples taken at random locations within a field.

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.

Influence of seasons and sampling strategy on assessment of bioaerosols in sewage treatment plants in Switzerland

An assessment of sewage workers' exposure to airborne cultivable bacteria, fungi and inhaled endotoxins was performed at 11 sewage treatment plants. We sampled the enclosed and unenclosed treatment areas in each plant and evaluated the influence of seasons (summer and winter) on bioaerosol levels. We also measured personal exposure to endotoxins of workers during special operation where a higher risk of bioaerosol inhalation was assumed. Results show that only fungi are present in significantly higher concentrations in summer than in winter (2331 +/- 858 versus 329 +/- 95 CFU m(-3)). We also found that there are significantly more bacteria in the enclosed area, near the particle grids for incoming water, than in the unenclosed area near the aeration basins (9455 +/- 2661 versus 2435 +/- 985 CFU m(-3) in summer and 11 081 +/- 2299 versus 2002 +/- 839 CFU m(-3) in winter). All bioaerosols were frequently above the recommended values of occupational exposure. Workers carrying out special tasks such as cleaning tanks were exposed to very high levels of endotoxins (up to 500 EU m(-3)) compared to routine work. The species composition and concentration of airborne Gram-negative bacteria were also studied. A broad spectrum of different species within the Pseudomonadaceae and the Enterobacteriaceae families were predominant in nearly all plants investigated. [Authors]

Sampling strategies for tropical forest nutrient cycling studies: a case study in São Paulo, Brazil

The precise sampling of soil, biological or micro climatic attributes in tropical forests, which are characterized by a high diversity of species and complex spatial variability, is a difficult task. We found few basic studies to guide sampling procedures. The objective of this study was to define a sampling strategy and data analysis for some parameters frequently used in nutrient cycling studies, i. e., litter amount, total nutrient amounts in litter and its composition (Ca, Mg, Κ, Ν and P), and soil attributes at three depths (organic matter, Ρ content, cation exchange capacity and base saturation). A natural remnant forest in the West of São Paulo State (Brazil) was selected as study area and samples were collected in July, 1989. The total amount of litter and its total nutrient amounts had a high spatial independent variance. Conversely, the variance of litter composition was lower and the spatial dependency was peculiar to each nutrient. The sampling strategy for the estimation of litter amounts and the amount of nutrient in litter should be different than the sampling strategy for nutrient composition. For the estimation of litter amounts and the amount of nutrients in litter (related to quantity) a large number of randomly distributed determinations are needed. Otherwise, for the estimation of litter nutrient composition (related to quality) a smaller amount of spatially located samples should be analyzed. The determination of sampling for soil attributes differed according to the depth. Overall, surface samples (0-5 cm) showed high short distance spatial dependent variance, whereas, subsurface samples exhibited spatial dependency in longer distances. Short transects with sampling interval of 5-10 m are recommended for surface sampling. Subsurface samples must also be spatially located, but with transects or grids with longer distances between sampling points over the entire area. Composite soil samples would not provide a complete understanding of the relation between soil properties and surface dynamic processes or landscape aspects. Precise distribution of Ρ was difficult to estimate.

The use of scaled semivariograms to plan soil sampling in sugarcane fields

Spatial sampling designs used to characterize the spatial variability of soil attributes are crucial for science studies. Sample planning for the interpolation of a regionalized variable may use several criteria, which could be best selected from an estimated semivariogram from a previously established grid. The objective of this study was to optimize the procedure for scaled semivariogram use to plan soil sampling in sugarcane fields in the Alfisol and Oxisol regions of Jaboticabal Town in So Paulo State, Brazil. A scaled semivariogram for several soil chemical attributes was estimated from the data obtained from two grids positioned on a sugarcane field area, sampled at a depth of 0.0-0.5 m. The research showed that regular grids with uniform intervals did not express the real spatial variability of the soil attributes of Oxisols and Alfisols in the study area. The calculated final sampling density based on the scaled parameters of the semivariogram was one sample for each 2 ha in Area 1 (convex landscape) and one sample for each 1 ha in Area 2 (linear landscape), as indicated by SANOS 0.1 software. The combined use of the simulation programs and scaled semivariograms can be used to define sampling points. These results may help in soil fertility mapping and thereby improve nutrient management in sugarcane crops.

Design of nest access grids and perches in front of the nests: Influence on the behavior of laying hens.

In aviary systems for laying hens, it is important to provide suitable nest access platforms in front of the nests, allowing hens to reach and explore each of the nests easily. This access platform is needed to achieve good nest acceptance by the hens and thereby prevent mislaid eggs. In the present experiment, the behavior of hens using 2 different nest access platforms, a plastic grid and 2 wooden perches, was examined. Furthermore, the nests were placed on both sides of the aviary rack (corridor side and outdoor side), either integrated into the aviary rack itself (integrated nest; IN) or placed on the walls of the pens (wall nest; WN), resulting in a 2 × 2 factorial design Four thousand five hundred white laying hens were housed in 20 test pens. The eggs in the nests and mislaid eggs were collected daily, and the behavior of hens on the nest accesses was filmed during wk 25 and 26, using focal observation and scan sampling methods. More balancing, body contact, and agonistic interactions were expected for nests with perches, whereas more walking and nest inspections were expected for nests with grids. There were more mislaid eggs and balancing found in pens equipped with nests with wooden perches. More agonistic interactions and balancing, less standing, and a longer duration of nest inspection were found with the WN compared with the IN. Interactions between platform design and position of the nests were found for duration of nest visits, body contact, and walking, with the highest amount for WN equipped with plastic grids. Nests on the corridor side were favored by the hens. Nest-related behaviors, such as nest inspection, standing, and walking, decreased over time as did the number of hens on the nest accesses, whereas sitting increased. These results indicate that the hens had more difficulties in gripping the perches as designed. The lower number of hens on the nest access platforms in front of IN may be due to a better distribution around nests and tier changes within the aviary rack. Based on these results, grids rather than perches provide for improved nesting behavior.