ABSTRACT Outliers and inliers often arise during sample data acquisition. While outliers represent anomalous observations, inliers are erroneous data points within the main body of the dataset. It was aimed to conduct a systematic literature study (SLS) to survey methods and software employed for outlier and inlier removal, particularly within exploratory data analysis. The study was conducted in three phases: (i) systematic literature mapping (SLM), (ii) snowballing (SB), and (iii) SLR. Initially, 772 scientific studies were identified, subsequently narrowed down to 86 after applying selection criteria. Backward (BSB) and forward (FSB) snowballing further yielded 16 studies, resulting in a final pool of 102 studies for analysis. It was identified three outlier removal techniques (Chebyshev’s inequality, boxplot, and principal component analysis), one inlier removal technique (local Moran’s index), and thirteen commonly used software. acquisition observations dataset SLS (SLS analysis phases i (i SLM, SLM , (SLM) ii (ii SB, SB (SB) iii (iii SLR Initially 77 8 criteria BSB (BSB FSB (FSB 1 10 Chebyshevs Chebyshev s inequality boxplot analysis, analysis) local Morans Moran index, index index) (SLM (SB 7

RESUMO Outliers e inliers aparecem frequentemente na aquisição de dados amostrais. Outliers são observações anômalas e inliers são dados errôneos no interior do conjunto de dados. Esta pesquisa teve como objetivo realizar um estudo sistemático da literatura (SLS) para levantar os métodos utilizados para remoção de outliers e inliers e os softwares utilizados na análise exploratória dos dados. Este estudo sistemático da literatura foi realizado em três etapas: (i) mapeamento sistemático da literatura (SLM), (ii) bola de neve (SB) e (iii) revisão sistemática da literatura (SLR). Setecentos e setenta e dois estudos científicos foram obtidos e reduzidos para 86 após seleção. Foram acrescentados mais dezessete estudos selecionados por bola de neve (bola de neve para trás (BSB) e bola de neve para frente (FSB)), o que resultou em 102 estudos utilizados nesta pesquisa. Foram observadas três técnicas de remoção de outliers (desigualdade de Chebyshev, boxplot e análise de componentes principais), uma única técnica de remoção de inliers (índice de Moran local) e treze softwares. amostrais SLS (SLS etapas i (i SLM, SLM , (SLM) ii (ii SB (SB iii (iii SLR. SLR . (SLR) 8 seleção BSB (BSB FSB, FSB (FSB)) 10 desigualdade Chebyshev principais, principais principais) índice local (SLM (SLR (FSB) 1 (FSB

ABSTRACT Irrigated agriculture is considered one of the most important techniques for agricultural production in the world, despite being the sector that uses water resources the most. Thus, technological advances have been constantly developed and applied in the sector as a way to improve water management and reduce its consumption, considering that some cultures are affected by economic costs and feasibilities of production, because the water cost has increased significantly in recent years. Remote sensing is a technology that has provided good results for estimating evapotranspiration (ETr) in producing areas, enabling more efficient irrigation management for low costs. This study aimed to answer the question about how the generation of irrigation management zones (IMZs) can be used to optimize water use and improve irrigation systems and which ETr models are better for it. Three ETr estimation models (Surface Energy Balance Algorithm for Land [SEBAL], Mapping Evapotranspiration at High Resolution and with Internalized Calibration [METRIC] and Simple Algorithm for Evapotranspiration Retrieving [SAFER]) were used to design IMZs in center pivot cultivated with cotton. The design of IMZs with ETr data proved to be a viable alternative, making it possible to improve water management in irrigated systems, reducing costs with viability irrigation for cotton and others cultures. The METRIC method displayed the greatest ease in obtaining the data used in the ETr estimation for generating the IMZs.

RESUMO Esse estudo teve como objetivo avaliar e identificar a variabilidade espacial da erosão laminar em áreas que utilizam agricultura de precisão. Para tanto, foram utilizados dados de solo de três propriedades da região oeste do Paraná, uma no município de Céu Azul (área A) e duas em Serranópolis do Iguaçu (áreas B e C). Para identificar os dados discrepantes (outliers), foi realizada uma análise da dispersão dos quartis por meio de um gráfico box-plot. A verificação da normalidade dos dados foi realizada com base na aplicação do teste de Kolmogorov-Smirnov. Em seguida, foi realizada uma análise espacial, utilizando o software AgDataBox-Map. Após a obtenção dos mapas temáticos, foram estimados os parâmetros da Equação Universal de Perdas de Solo. Após essa etapa, foram calculados os fatores de estimativa das perdas de solo e gerados mapas temáticos, os quais foram integrados por álgebra cartográfica, resultando no modelo de identificação de áreas suscetíveis à erosão. A área A apresentou maior suscetibilidade (estimativa de perda de solo = 0-200 t ha-1 por ano) quando comparada às áreas B e C (estimativa de perda de solo = 0-150 t ha-1 por ano), porém, todas as áreas apresentam baixa suscetibilidade à erosão.

ABSTRACT The study aimed to identify and evaluate the spatial variability in laminar erosion in areas using precision agriculture tools. Soil data from three properties in the western region of Paraná state, Brazil, were used: one in the municipality of Céu Azul (area A) and two in Serranópolis do Iguaçu (areas B and C). To identify discrepant data (outliers), analysis of the dispersion of quartiles was performed using a box-plot graph. Data normality was verified using the Kolmogorov-Smirnov test. A spatial analysis was performed using AgDataBox-Map software. The parameters of the universal soil loss equation were estimated and used with map algebra to produce a model to identify areas susceptible to erosion. Area A (soil loss estimate = 0-200 t ha-1 per year) presented greater susceptibility to erosion than areas B and C (soil loss estimate = 0-150 t ha-1 per year); however, all areas had a low susceptibility to erosion.

ABSTRACT Spatial variability management of soil chemical attributes is one of the approaches to be employed in the face of the constant challenge of increasing agricultural yield to meet world demand. In this sense, precision agriculture has as one of its tools the application of inputs at varying rates, which seeks to determine the ideal amount of fertilizer at each point of the crop, contrary to the conventional recommendation approach based on average values. In this context, this work studied the fertilizer recommendation methods used in site-specific nutrient management and the calculation methodologies for N, P, and K recommendations. For this purpose, a systematic literature study (SLS), consisting of systematic literature mapping, snowballing, and systematic literature review was performed. The analyzed studies were grouped into five domains (precision agriculture, soil fertility, site-specific nutrient application, fertilizer recommendation methods, and recommendation software for site-specific nutrient application). As a result, the SLS identified 12 methods for recommending N, nine for recommending P, and six for recommending K, in addition to five computer programs for precision agriculture that perform fertilizer recommendations at varying rates.

ABSTRACT: Optimization of N management is one of the great challenges to be overcome in grain production, as it is directly related to productivity and can also cause environmental damage. Precision agriculture aims to solve this problem by applying nitrogen fertilizer at varying rates. Reflectance sensors are instruments capable of estimating N needs in various crops, including grain crops. However, it is not clear how these sensors perform under varying solar radiation and cloud cover, due to a lack of research on their temporal variability. Thus, this study examined the temporal variability of the NDVI (normalized difference vegetation index), as measured by an active reflectance sensor, in both soybean and wheat crops. The NDVI data were collected using a GreenSeeker sensor every 15 minutes over 12 or 14 consecutive hours. Incident solar radiation was recorded using an Instrutherm MES-100 pyranometer. In all experiments in soybean and wheat, NDVI was negatively influenced by irradiation, showing higher values at the beginning and end of the day. Changes in cloud cover also affected NDVI values during the experiments.

Este trabalho objetivou quantificar a dependencia espacial de atributos quimicos do solo, indice de clorofila em estadios V3, V6 e R1 e produtividade da cultura do milho, e analisar a correlacao espacial entre estes. Foram avaliados atributos do solo nos anos 2010 e 2011 e indice de clorofila (com clorofilometro) e produtividade nos anos de 2011 e 2012.Os atributos estudados apresentaram estrutura de dependencia espacial, permitindo o manejo localizado. Tambem se observou a existencia de correlacao espacial entre os atributos, demonstrando a importancia de estudos nesta area, identificando-se quais atributos exercem influencia no desenvolvimento da cultura.

This study aimed to quantify the spatial dependence of soil chemical properties, chlorophyll index in stages V3, V6 and R1 and corn yield, and analyze the spatial correlation between them. Soil properties were evaluated in 2010 and 2011 and the rate of chlorophyll (with clorofilometro) and productivity in the years 2011 and 2012. The studied attributes presented spatial dependence structure, allowing the site-specific management. It also observed the existence of spatial correlation between attributes, demonstrating the importance of studies in this area, identifying which attributes influence the development of culture.

Este estudio tuvo como objetivo cuantificar la dependencia espacial de las propiedades quimicas del suelo, indice de clorofila en las etapas V3, V6 y R1 y productividad de maiz, y analizar la correlacion espacial entre ellas. Las propiedades del suelo fueron evaluadas en 2010 y 2011 y el indice de clorofila (con clorofilometro) y la productividad en los anos 2011 y 2012. Los atributos presentaron dependencia espacial, lo que permite la aplicacion localizada. Tambien se observo la existencia de correlacion espacial entre atributos, demostrando la importancia de estudios en esta area, buscando identificar atributos que influyen en el desarrollo de la cultura.

RESUMO A agricultura de precisão proporciona aos agricultores identificar e tratar de forma adequada as variações encontradas na área agrícola. As zonas de manejo (ZMs) permitem a implantação da agricultura de precisão de forma viável e relativamente mais econômica. A forma mais importante para definir ZMs é usando o algoritmo fuzzy C-means. Um problema consiste em como selecionar a variável a ser usada como layer de entrada no processo fuzzy. Assim, o objetivo deste trabalho, foi testar a hipótese de que variáveis redundates podem prejudicar o processo de delineamento de ZMs. Este trabalho foi desenvolvido em uma área de 19,6 ha e 15 agrupamentos de ZMs foram gerados por meio do o algoritmo fuzzy C-means, dividindo-se em duas a cinco classes. Cada agrupamento usou uma composição diferente de variáveis, que são os atributos cobre, silte, argila, e altitude. Foi encontrado que algumas combinações dessas variáveis produziu melhores ZMs. Nenhuma combinação de variáveis produziu desempenho estatisticamente melhor que a ZM gerada apenas com as variáveis não redundantes. Assim, as variáveis redundantes podem ser descartadas. O agrupamento com todas as variáveis não forneceu maior separação e organização dos dados entre as classes de ZM, não sendo recomendado.

ABSTRACT Precision agriculture (PA) allows farmers to identify and address variations in an agriculture field. Management zones (MZs) make PA more feasible and economical. The most important method for defining MZs is a fuzzy C-means algorithm, but selecting the variable for use as the input layer in the fuzzy process is problematic. BAZZI et al. (2013) used Moran’s bivariate spatial autocorrelation statistic to identify variables that are spatially correlated with yield while employing spatial autocorrelation. BAZZI et al. (2013) proposed that all redundant variables be eliminated and that the remaining variables would be considered appropriate on the MZ generation process. Thus, the objective of this work, a study case, was to test the hypothesis that redundant variables can harm the MZ delineation process. BAZZI This work was conducted in a 19.6-ha commercial field, and 15 MZ designs were generated by a fuzzy C-means algorithm and divided into two to five classes. Each design used a different composition of variables, including copper, silt, clay, and altitude. Some combinations of these variables produced superior MZs. None of the variable combinations produced statistically better performance that the MZ generated with no redundant variables. Thus, the other redundant variables can be discredited. The design with all variables did not provide a greater separation and organization of data among MZ classes and was not recommended.

During the last few years, yield maps have become economically feasible for farmers due to technological advances in precision agriculture. However, evidence of yield profitability is still uncertain, and variability in yield has seldom been correlated to variability in profits. Differently from yield maps, profit maps can supply additional information about the economic return for each particular area of a field. The objective of the present work was to study how management decisions can be facilitated by transforming yield-map datasets into profit maps and the importance of the selection of interpolator type. Yield and profit maps were generated for each data set (three soybean fields and one corn field) using the inverse of the distance (ID), the inverse of the square of the distance (IDS) and kriging (KRG) as interpolation methods. It can be concluded that profit maps are an important tool for the diagnosis of the spatial variability of economic return because they can assist farmers with management decision-making. The impact of the interpolator type was less than 200 kg ha-1 for the yield and US$ 30 ha-1 for the profit, indicating that the choice of interpolator type is of secondary importance. In addition, the profit maps showed large variability that would not be easily found if only yield maps were analyzed.

Durante los últimos años, los mapas de rendimiento se han convertido económicamente viables para los agricultores, debido a los avances tecnológicos en la agricultura de precisión. Sin embargo, la evidencia de la rentabilidad del rendimiento, aún es desconocida, y la variabilidad del rendimiento rara vez ha sido correlacionada con la variabilidad de la rentabilidad. A diferencia de los mapas de rendimiento, los mapas de rentabilidad pueden suministrar información adicional relacionada con el retorno económico para cada área particular del campo. El objetivo del presente trabajo ha sido estudiar como las decisiones de gestión pueden ser facilitadas por la transformación de los conjuntos de datos de los mapas del rendimiento en mapas de beneficio, bien como la importancia de la selección de los mapas de interpolación. Para cada conjunto de datos (tres zonas de soya e un de maíz), mapas de rendimiento e de beneficio, fueron generados utilizando los métodos de interpolación de distancia inversa (DI), de distancia cuadrada inversa (DCI) y de Krigagem (KRG). Se puede concluir que los mapas de beneficio y rentabilidad son herramientas importantes para el diagnóstico de la variabilidad espacial del rendimiento económico, ya que ayudan a los agricultores en la gestión de la toma de decisiones. El impacto del tipo de interpolación fue de menos de 200 kg por hectárea en el rendimiento, US$ 30 en el beneficio, significando que la elección puede ser colocada en el segundo plano. En adicción, los mapas de beneficio mostraron una gran variabilidad que no haría fácil encontrar, solamente por el análisis de los mapas de productividad.

As plantas de cobertura são importantes na melhoria da qualidade do solo, cujos atributos frequentemente apresentam dependência espacial. Assim, o objetivo deste trabalho foi avaliar o efeito de plantas de cobertura em manejo de inverno nos atributos físicos do solo e a produtividade da soja por meio de mapas temáticos. Os tratamentos foram cinco manejos de inverno: aveia-preta; consórcio 1 (nabo forrageiro e aveia-preta); consórcio 2 (nabo forrageiro, aveia-preta e ervilhaca comum); trigo e testemunha (pousio). Avaliaram-se a macroporosidade, a microporosidade, a porosidade total, a densidade e o teor de água do solo, na profundidade de 0 - 0,1 m após o manejo das coberturas de inverno. A soja foi semeada em toda a área no verão, após o manejo das coberturas de inverno e determinada sua produtividade para cada tratamento. Elaboraram-se mapas temáticos para os tratamentos, os quais foram comparados com a testemunha, utilizando-se do coeficiente de desvio relativo (CDR). As plantas de cobertura melhoraram a macroporosidade e a porosidade total do solo em algumas regiões dentro da área em estudo. A aveia-preta foi mais eficiente em manter mais elevado o teor de água do solo e demonstrou maior potencial para reduzir a densidade do solo.

Cover crops are important for improving soil quality. However, soil properties usually have some spatial dependence. Thus, this study aimed to evaluate the effect of winter cover crops on physical properties of soil and soybean yields using thematic maps. Five winter treatments were used: black oats; intercropping 1 (forage turnips and black oats); intercropping 2 (forage turnips, black oats and common vetch); wheat; and control. Macroporosity, microporosity, total porosity, bulk density and water content of the soil from 0 - 0.1 m depths were evaluated after the winter cover crop management. Soybeans were sown over the entire area in the summer after the winter cover crop management, and the soybean yield was determined for each treatment. Maps for each treatment were created and compared to the control treatment using the relative deviation coefficient (RDC). The cover crops improved the total macroporosity of the soil in some regions of the study area. The black oats were more efficient at maintaining higher water content of the soil, and it can be used to decrease the bulk density.

Diversos equipamentos e metodologias vêm sendo desenvolvidos para tornar a agricultura de precisão disponível, especialmente considerando o alto custo de sua implantação e de amostragem. Uma possibilidade interessante é definir a área em unidades menores de produção que podem ser tratadas de maneira diferente, servindo como fonte de recomendação e análise. Assim, o presente estudo utilizou propriedades físicas e químicas do solo e de produtividade visando à geração de unidades de manejo, a fim de identificar se estas podem ser usadas como recomendação e análise. Unidades de manejo foram geradas pelo algoritmo Fuzzy C-Means, e sua avaliação foi realizada por meio da redução da variância e realização de testes de comparação de médias. A divisão da área em duas unidades de manejo foi considerada adequada, e as médias apresentaram-se distintas da maioria das propriedades do solo e produtividade; a metodologia utilizada permitiu a geração de unidades de manejo que podem servir como fonte de recomendação e análise do solo. Apesar de a eficiência relativa demonstrar que houve redução da variância para todos os atributos na divisão em três sub-regiões, a ANOVA não apresentou diferenças significativas entre as unidades de manejo.

Several equipments and methodologies have been developed to make available precision agriculture, especially considering the high cost of its implantation and sampling. An interesting possibility is to define management zones aim at dividing producing areas in smaller management zones that could be treated differently, serving as a source of recommendation and analysis. Thus, this trial used physical and chemical properties of soil and yield aiming at the generation of management zones in order to identify whether they can be used as recommendation and analysis. Management zones were generated by the Fuzzy C-Means algorithm and their evaluation was performed by calculating the reduction of variance and performing means tests. The division of the area into two management zones was considered appropriate for the present distinct averages of most soil properties and yield. The used methodology allowed the generation of management zones that can serve as source of recommendation and soil analysis; despite the relative efficiency has shown a reduced variance for all attributes in divisions in the three sub-regions, the ANOVA did not show significative differences among the management zones.

Os atributos químicos e físicos do solo apresentam, frequentemente, dependência espacial. Portanto, é imprescindível estudar sua variabilidade espacial visando ao manejo mais preciso da lavoura, a redução do impacto ambiental, a otimização do uso de máquinas e o aumento da eficiência das equipes de trabalho. Assim, o objetivo deste trabalho foi avaliar a variabilidade espacial de atributos físicos de um Latossolo Vermelho distroférrico típico após o cultivo de soja. Os dados foram coletados em uma área de 2 ha cultivada sob sistema de plantio direto utilizando-se uma grade amostral de 30 x 33 m, em que foram avaliados os parâmetros: macroporosidade, microporosidade, porosidade total, densidade e umidade do solo, na camada de 0 - 0,1 m; após as análises exploratória e geoestatística, os modelos teóricos foram ajustados aos semivariogramas para cada atributo; já a interpolação dos dados foi realizada utilizando-se a krigagem ordinária, verificando-se que todas as variáveis estudadas apresentaram dependência espacial, que variou entre moderada (macroporosidade, porosidade total e umidade do solo) e forte (microporosidade e densidade do solo), em função do coeficiente de efeito pepita.

The chemical and physical properties of a soil often show a spatial dependence. Thus, it is important to study their spatial variability, aiming to improve the crop management, to reduce the environmental impact, to optimize the machinery use, and to increase the efficiency of staff. So, this study aimed to evaluate the spatial variability of the physical properties of a distroferric Red Latosol (Typic Haplorthox) after soybean crop. The data were collected in an area of 2 ha managed under no tillage system, using a sampling grid of 30 x 33 m, where it was evaluated the macro, micro and the total porosity, the bulk density and the soil water content in 0 - 0.1 m layer. After the exploratory and geostatistical analyses, the theoretical models were set to the semivariograms for each attribute. The data interpolation was performed using the ordinary kriging. All the variables showed spatial dependence, which ranged from moderate (macroporosity, total porosity and soil water content) to strong (microporosity and bulk density), according to the nugget effect coefficient.

O mapa de produtividade representa a variabilidade espacial das características de uma área cultivada e permite intervir na produção dos anos posteriores, na aplicação diferenciada de insumos. Este trabalho teve por objetivo verificar a influência da densidade amostral e do tipo de interpolação na exatidão dos mapas de produtividade, gerados a partir da amostragem manual de grãos, solução que pode ser adotada quando um monitor não pode ser utilizado. Um mapa de produtividade foi obtido com monitor de colheita comercial em lavoura de milho, a partir do qual foram estabelecidas 84 grades e, por meio de três interpoladores, o inverso da distância ao quadrado, inverso da distância e krigagem geraram-se, assim, 252 mapas de produtividade, que foram então comparados com o original, utilizando-se o coeficiente de desvio relativo (CDR) e o índice kappa. A perda de informação do mapa de produtividade aumentou à medida em que se diminuiu a densidade amostral e foi dependente do método de interpolação utilizado. Um modelo de regressão múltipla foi ajustado à variável CDR, em função das variáveis: índice de variabilidade espacial e densidade amostral. Este modelo teve a finalidade de auxiliar o agricultor na definição da densidade amostral e permitir obter-se, manualmente, o mapa de produtividade, caso ocorram problemas eventuais no monitor de colheita.

Yield mapping represents the spatial variability concerning the features of a productive area and allows intervening on the next year production, for example, on a site-specific input application. The trial aimed at verifying the influence of a sampling density and the type of interpolator on yield mapping precision to be produced by a manual sampling of grains. This solution is usually adopted when a combine with yield monitor can not be used. An yield map was developed using data obtained from a combine equipped with yield monitor during corn harvesting. From this map, 84 sample grids were established and through three interpolators: inverse of square distance, inverse of distance and ordinary kriging, 252 yield maps were created. Then they were compared with the original one using the coefficient of relative deviation (CRD) and the kappa index. The loss regarding yield mapping information increased as the sampling density decreased. Besides, it was also dependent on the interpolation method used. A multiple regression model was adjusted to the variable CRD, according to the following variables: spatial variability index and sampling density. This model aimed at aiding the farmer to define the sampling density, thus, allowing to obtain the manual yield mapping, during eventual problems in the yield monitor.