ABSTRACT This study aimed to identify whether there is an association between types of storage and categories of commercialization and use of grains, group the types of storage with the categories of commercialization, and group static and dynamic capacities of the units in the state of Paraná, Brazil. The data were obtained from the Brazilian Registry System of Storage Units for the 2014/2015 season. The association between variables under study was carried out with the chi-square test of independence and correspondence analysis. The cluster analysis consisted of the unweighted pair group method with arithmetic mean and considered a measure of mixed dissimilarity obtained for sets composed of qualitative and quantitative variables. A significant association was observed between the type of silo battery with the grain usage characteristic and cooperatives that commercialize the grain for use in the domestic market (CICOOPT); between the bulk warehouse and cooperatives that commercialize the grain in the foreign market (CIECOOPT); and between the silo and grain sellers who commercialize the grain in the domestic market (CI). Most types of storage units were grouped in Group 1, with a predominance of the CI characteristic and small to medium size static and dynamic capacity.

RESUMO No Estado do Paraná, Brasil, não há grandes mudanças nas áreas cultivadas com culturas anuais, principalmente devido a leis ambientais que não permitem expansões para novas áreas. Há grande contribuição das culturas anuais para a demanda doméstica de alimentos e econômica nas exportações. Assim, a área e distribuição das culturas anuais são informações de grande importância. Novas metodologias, como data mining, estão sendo testadas com o objetivo de analisar e melhorar seu potencial de uso para classificação do uso e cobertura da terra. Neste estudo, foram utilizados os classificadores decision tree e random forest com métricas temporais de Normalized Difference Vegetation Index (NDVI) em imagens do Operational Land Imager (OLI)/ Landsat-8. Os resultados foram comparados com os métodos tradicionais (imagens espectrais e classificador Maximum Likelihood Classifier - MLC). Inicialmente, foram mapeadas sete classes (corpos d’água, cana-de-açúcar, área urbana, culturas anuais, floresta, pastagem e áreas de reflorestamento) e posteriormente apenas duas classes foram consideradas (culturas anuais e outras classes). Ao classificar os sete alvos, ambos os métodos tiveram resultados correspondentes, mostrando exatidão global próxima a 84%. As métricas temporais de NDVI mostraram a acurácia do produtor e do usuário para a classe de cultura de 86 e 100%, respectivamente. No entanto, considerando-se apenas duas classes, as métricas temporais do NDVI alcançaram exatidão global próxima a 98% e a acurácia do produtor e do usuário acima de 94%.

ABSTRACT In the state of Paraná, Brazil, there are no major changes in areas cultivated with annual crops, mainly due to environmental laws that do not allow expansions to new areas. There is a great contribution of the annual crops to the domestic demand of food and economic demand in the exports. Thus, the area and distribution of annual crops are information of great importance. New methodologies, such as data mining, are being tested with the objective of analyzing and improving their potential use for classification of land use and land cover. This study used the classifiers decision tree and random forest with Normalized Difference Vegetation Index (NDVI) temporal metrics on images from Operational Land Imager (OLI)/Landsat-8. The results were compared with traditional methods spectral images and Maximum Likelihood Classifier (MLC). At first, seven classes were mapped (water bodies, sugarcane, urban area, annual crops, forest, pasture and reforestation areas); then, only two classes were considered (annual crops and other targets). When classifying the seven targets, both methods had corresponding results, showing global accuracy near 84%. NDVI temporal metrics showed producer’s and user’s accuracy for the annual crop class of 86 and 100%, respectively. However, if considering only two classes, the NDVI temporal metrics reached global accuracy of near 98% and producer’s and user’s accuracy above 94%.

ABSTRACT The aim of this study is to identify homogeneous rainfall zones in the winter and summer 1st and 2nd crops, in the state of Paraná, Brazil. The zones were defined by clustering using the expectation-maximization (EM) algorithm to transform seasonal rainfall series. Monthly average rainfall data collected from 157 weather stations for 20 years (1996 to 2015) were employed. The results show that the number of homogeneous zones varied among growing seasons. The summer crop presented two clusters, with rainfall averages of 1489 and 1925 mm; the second crop presented four clusters, with averages of 1849, 1004, 1454, and 1182 mm; and the winter crop had three clusters, with averages of 969, 1498, and 1171 mm. Clustering was a useful instrument to identify geographical regions with similar rainfall regimes during different growing seasons in the state of Paraná. Rainfall distribution was more homogeneous in the summer crop. In all crops analyzed, the clusters with the lowest rainfall rate were present in the northwestern, northern center, and northern pioneer of the state of Paraná, whereas the clusters with the highest rainfall rate were found in the coastal regions.

ABSTRACT The application of nitrogen (N) fertilizer is complex and expensive, so its correct management has financial and environmental benefits. The use of optical proximity sensors is a promising technique. However, the movement of the agricultural machinery or of the person carrying the sensor will result in height differences and/or different tilt and twist angles with respect to the canopy. We considered whether these variations would affect the reflectance measurement. In this study, we took normalized difference vegetation index (NDVI) readings of a wheat canopy, to which 90 kg ha-1 of urea had been applied in stage 5, and observed the NDVI in stages 6, 8 and 10.5. We also tested soybeans, to which 90 kg ha-1 of urea had been applied in stage R1, and took NDVI readings in stages R2 and R5. Our goal was to study the effects of the position of an active reflectance sensor (GreenSeeker) on the NDVI index at different heights and at different angles to the canopy. We observed that the height of the sensor affected the NDVI depending on the stage of the plant and that angles up to 15° of the sensor did not directly affect the readings.

ABSTRACT The goal of this study was to use the spatial bootstrap method to model the spatial dependence structure of soybean yield and soil chemical attributes in an agricultural area. The study involved developing confidence intervals in probability plots to determine the probability distributions assumed by the data; determine the empirical distributions of the semivariances and model parameters, allowing to obtain statistics and confidence intervals; and to construct maps for the variables. The quantile-quantile plots indicated that the data follows a normal distribution. The confidence intervals for the semivariances helped to model the spatial dependence structure, and the descriptive statistics of the bootstrap replicates of the model parameters allowed to test the consistency of the estimates. The soil chemical attributes (calcium, potassium, and organic matter) were at levels suitable for soybean cultivation. However, the pH was below the ideal range in most of the study area, and water stress during cultivation decreased the mean yield. Therefore, according to the results, a recommendation to the farmer is to correct the soil pH to increase the yield.

ABSTRACT This work aimed to study the spatial autocorrelation of the total static capacity storage, the total number of warehouses in 2013/2014 (CONAB) and the average of the total grain production (soybean, corn 1st and 2nd crops and wheat) in the harvest years 2008/2009 to 2013/2014 (SEAB) in Paraná State, Brazil. The study was based on Moran's global autocorrelation index, Moran's local and Moran's bivariate correlation. It was possible to identify regions with low and high total grain production. There was positive spatial autocorrelation for the Total Static Storage Capacity (TSSC) and Total Quantity of Warehouses (TQW). For the total grain production, significant spatial autocorrelation were found. The total static storage capacity showed similarity between the studied regions. When evaluating the bivariate spatial correlation between Total Production of Harvested Grains (TPHG) in relation to the total static storage capacity and total quantity of warehouses, the presence of positive spatial correlation was observed. The results indicated that Moran's global autocorrelation and local indexes showed significant patterns of spatial autocorrelation, as well as the bivariate spatial correlation indexes in the studied variables.

Abstract The uniformity of water application is an important factor in the evaluation of sprinkler irrigation systems. This uniformity depends on the type of sprinkler and its operating conditions, such as the arrangement and spacing between the sprinklers in the area; velocity and wind direction during the period of water application and the pressure variation of the irrigation system. The objective of this study was to model, analyze and compare the structure of spatial dependence, as well as the spatial variability of the water depths applied by a sprinkler irrigation system with compensating and non-compensating sprinklers, using geostatistical methods and measurements of accuracy or similarity between the applied water depth maps. The experiment was carried out in an agricultural area, in the city of Cascavel-Paraná-Brazil. A total area of 10 × 10 m was used, with 04 compensating and 04 non-compensating sprinklers installed at a height of 1.5 m. For each type of sprinkler, water levels were measured in 100 collectors spaced 1 × 1 m in the study area in 32 trials. On each test sprinkling was carried out for one hour. The conditions of wind, temperature and air humidity were evaluated at the beginning of each test and at 10-minute intervals with a climatological station. As the geostatistical analysis showed the existence of directional trends, the coordinates were incorporated as covariates to the linear spatial model in the study of the spatial dependence of the average depth of the irrigation water for the two types of sprinklers. The spatial dependence structure that best fits the data when using the compensating sprinklers was the Gaussian model and when the non-compensating sprinklers were used, it was the exponential model. The spatial variability maps of average irrigation water depth (mm) of the trials, obtained by universal kriging, revealed that for both sprinklers there was an increase in the mean level average values in the northwest-southeast direction (135° in the azimuth system) in the area under study, influenced by wind direction and velocity during the execution of the experiment.

ABSTRACT This study aims to quantify the uncertainties associated to the parameters of a Gaussian spatial linear model (GSLM) and the assumption of normality residuals in the modeling of the spatial dependence of the soybean yield as a function of soil chemical attributes. The spatial bootstrap methods were used to determine the point and interval estimators associated with the model parameters. Hypothesis tests were carried out on the significance of the model parameters and the quantile-quantile probability plot was elaborated to verify the data normality. The uncertainties associated to the parameters of the spatial dependence structure were quantified and the potassium content, phosphorus content and soil pH covariates were significant to explain the soybean yield mean. These covariates were used in the elaboration of a new model, which provided the elaboration of a contour map of soybean yield. Analysis of the quantile-quantile plot indicated that soybean yield data follow a normal probability distribution.

RESUMO Aplicações de sensoriamento remoto na agricultura se apresentam como realidade muito promissora, contudo ainda há necessidade de pesquisas para a correta utilização dos dados espectrais. Objetivou-se, neste estudo, avaliar os padrões espectro-temporais de onze cultivares de trigo (Triticum aestivum L.). O experimento foi conduzido, em Cascavel, PR, no ano de 2014. Com o auxílio dos sensores terrestres GreenSeeker e FieldSpec 4 foram determinadas assinaturas espectrais e criados os perfis temporais do Índice de Vegetação por Diferença Normalizada (NDVI). Verificou-se diferenças estatísticas entre as cultivares de trigo por meio da análise de variância (ANOVA) e do teste de médias de Scott-Knott. Correlacionaram-se as produtividades de grãos obtidas com os fatores INSEY (IN-Season Estimate of Yield). O NDVI normalizado por graus-dias acumulados dos estádios de desenvolvimento Feekes 2 e 8, mostrou-se mais consistente na estimativa da produtividade de grãos, explicando aproximadamente 70% da variação. No estádio Feekes 10.1, as cultivares de trigo apresentaram padrões espectrais diferentes nas faixas do infravermelho próximo e médio. Isto sugere que essas bandas espectrais podem ser utilizadas para diferenciar as cultivares de trigo.

ABSTRACT Remote sensing applications in agriculture are presented as a very promising reality, but research is still needed for the correct use of spectral data. The objective of this study was to evaluate the spectral-temporal patterns of eleven wheat cultivars (Triticum aestivum L.). The experiment was conducted in Cascavel, PR, in the year 2014. With the help of the GreenSeeker and FieldSpec 4 terrestrial sensors, spectral signatures were determined and the temporal profiles of the Normalized Difference Vegetation Index (NDVI) were created. Statistical differences between wheat cultivars were analysed using analysis of variance (ANOVA) and Scott-Knott test. Grain yields obtained with INSEY (In-Season Estimate of Yield) factors were correlated. NDVI normalized by degree-days accumulated from the Feekes growth stages 2 and 8 showed to be more consistent in the estimation of grain yield, explaining approximately 70% of the variation. At the Feekes stage 10.1, wheat cultivars presented different spectral patterns in the near and medium infrared bands. This suggests that these spectral bands can be used to differentiate wheat cultivars.

ABSTRACT: The traditional per-pixel classification methods consider only spectral information, and may be limited. Object-based classifiers, however, also consider shape and texture, firstly segmenting the image, and then classifying individual objects. Thus, a Geographic Object-Based Image Analysis (GEOBIA) was compared in conjunction with data mining techniques and a traditional per-pixel method. A cut of Landsat-8, bands 2 to 7, orbit/point 223/77, located between the municipalities of Cascavel, Corbélia, Cafelândia and Tupãssi, in the west part of the state of Paraná, from 12/18/2013 was used. In the GEOBIA approach was realized image segmentation, spatial and spectral attribute extraction, and classification using the decision tree supervised algorithm, J48. For the per-pixel method, we used the supervised Maximum Likelihood Classifier. Both approaches presented equivalent results, with Kappa Index of 0.75 and Global Accuracy (GA) of 78.97% for the approach by GEOBIA and Kappa Index of 0.72 and GA of 77.44% for the perpixel classification. The classification by GEOBIA showed better accuracy for the soil, forest and soybean classes, and did not show the splash aspect, which visually improves the classification result.

ABSTRACT: Due to the difficulty in discriminating soybean and corn in mappings obtained by the time series of satellite images, this study aimed to apply the data mining techniques to separate soybean and corn. Pure pixels selection from Landsat-8 were extracted and used to build a standard spectro-temporal EVI profile for both crops. These profiles were obtained with the Timesat software and, further incorporated in the Weka software. Five out of eleven variables of the standard spectro-temporal EVI profile for each crop were found through the decision tree, a data mining technique. These five variables were sufficient to achieve the separation of soybean and corn crops with an accuracy of 96.3% and a kappa index of 0.92.

RESUMO O conhecimento de estimativas confiáveis de áreas cultivadas de cana-de-açúcar é imprescindível para o agronegócio brasileiro por auxiliar no desenvolvimento de políticas públicas, na determinação dos preços aos produtores pelas usinas e permitir estabelecer a logística de escoamento da produção. O objetivo deste trabalho foi realizar o mapeamento de área da cultura de cana-de-açúcar para o estado do Paraná a partir de imagens dos satélites Landsat/TM/OLI e IRS/LISS-3, para as safras de 2010/2011 a 2013/2014. Os mapeamentos foram realizados por meio da classificação supervisionada de Máxima verossimilhança (Maxver) obtendo-se, em média, uma exatidão global de 94,13% e índice kappa de 0,82. As correlações com os dados oficiais do IBGE variaram de moderada a forte (0,64 ≤ rs ≤ 0,80) com concordância (dr) média de 0,81. Houve aumento de 2,73% (18.630 ha) de área com cana-de-açúcar no Paraná entre 2010/2011 e 2013/2014.

ABSTRACT The knowledge on reliable estimates of areas under sugarcane cultivation is essential for the Brazilian agribusiness, since it helps in the development of public policies, in determining prices by sugar mills to producers and allows establishing the logistics of production disposal. The objective of this work was to develop a methodology for mapping the sugarcane crop area in the state of Paraná, Brazil, using images from the Landsat/TM/OLI and IRS/LISS-3 satellites, for the crop years from 2010/2011 to 2013/2014. The mappings were conducted through the supervised Maximum likelihood classification (Maxver) achieving, on average, an overall accuracy of 94.13% and kappa index of 0.82. The correlation with the official data of the IBGE ranged from moderate to strong (0.64 ≤ rs ≤ 0.80) with average agreement (dr) of 0.81. There was an increase of 2.73% (18,630 ha) in the area with sugarcane in Paraná between 2010/2011 and 2013/2014.

ABSTRACT The survey information from growing regions, the interaction with the vegetation index and climatic variables is of great importance in the search for soybean productivity increase. Paraná is the second largest soybean producer in Brazil and presents great spatial variability, both in periods of the crop cycle as in soil and climate. The objective of this study was to analyze the spatial correlation of soybean productivity, the enhanced vegetation index (EVI) and agrometeorological variables (water balance, global radiation and average temperature) in the state of Paraná, on a decendial scale, using the Moran global autocorrelation index between the 2010/2011 and 2012/2013 crop years. Similarity was found in the average productivities in 2010/2011 and 2012/2013. In 2011/2012 the state average was 2.38 t ha−1 lower in 10.19% compared to the national average, caused by the water deficit in flowering and grain filling phases. As a consequence, spatial autocorrelation indicated a higher similarity in productivity among municipalities with a Moran index of 0.735. The use of vegetation indices and agrometeorological variables allowed the identification of different sowing periods between regions and great climatic variability, influencing the soybean productivity.

ABSTRACT Precision Agriculture (PA) uses technologies with the aim of increasing productivity and reducing the environmental impact by means of site-specific application of agricultural inputs. In order to make it economically feasible, it is essential to improve the current methodologies as well as proposing new ones, in which data regarding productivity, soil, and compound indicators are used to determine Management Areas (MAs). These units are heterogeneous areas within the same region. With these methodologies, data mining (DM) techniques and algorithms may be used. In order to integrate DM techniques to PA, the aim of this study was to associate MAs created for soy productivity using the Fuzzy C-means algorithm by SDUM software over a 9.9-ha plot as the reference method. It was in opposition to the grouping of 2, 3, and 4 clusters obtained by the K-means classification algorithms, with and without the Principal Component Analysis (PCA), and the EM algorithm using chemical and physical data of the soil samples collected in the same area during the same period. The EM algorithm with PCA modeling had a superior performance than K-means based on hit rates. It is noteworthy that the greater the number of analyzed MAs, the lower the percentage of hits, in agreement with the result shown by SDUM, which shows that two MAs compose the best configuration for this studied area.

RESUMO O Estado do Paraná caracteriza-se por uma grande variabilidade de épocas de semeadura (DS) e, consequentemente, pelo desenvolvimento máximo vegetativo (DMDV), colheita (DC) e ciclo (CI) para a cultura da soja. O objetivo deste trabalho foi estimar essas datas para o período de primavera-verão do ano-safra de 2011/2012, por meio de séries temporais de imagens do Índice de Vegetação Realçado (do inglês Enhanced Vegetation Index - EVI) do sensor Modis (Moderate Resolution Imaging Spectroradiometer). Gerou-se um perfil espectrotemporal médio de EVI, considerando todos os pixels mapeados como soja dentro de cada município. Estes dados serviram de entrada no software Timesat para estimar os decêndios do ciclo da cultura (DS, DMDV, DC e CI) por municípios. Os resultados mostraram que existe grande variabilidade de datas de plantio em diferentes mesorregiões do Estado. Verificaram-se também divergências entre os resultados encontrados e os dados oficiais de DS e DC. A maior parte da semeadura (65,16%) esteve entre o terceiro decêndio de outubro e o primeiro decêndio de novembro. A maior parte da área de soja do Estado do Paraná (65,46%) teve seu DMDV em janeiro e colheita em março (53,92%).

ABSTRACT Soybean sowing can be performed at a various dates (DS) in Paraná state, with subsequent maximum vegetative development (DMDV), harvest (DC) and cycle (CI). The purpose of this study was to estimate these dates within spring-summer season from the 2011/ 2012 crop harvest, by means of temporal series of the Enhanced Vegetation Index (EVI) from Moderate Resolution Imaging Spectroradiometer (Modis) sensor. We generated an average EVI spectro-temporal profile considering all pixels mapped as soybean within each municipality. These data were used as input to the Timesat software for estimates of 10-day crop cycle periods (DS, DMDV, DC and CI). The results demonstrated great variability of dates at different areas of the Paraná state. We could also observe differences between our results and official estimates for DS and DC. Sowing was majorly (65.16%) performed in the third 10-day period of October and in the first 10-day period of November. Most of the Southern Paraná soybean area (65.46%) had their maximum vegetative development (MDV) in January and harvest in March (53.92%).