{"id":82071,"date":"2018-01-22T17:00:27","date_gmt":"2018-01-22T17:00:27","guid":{"rendered":"http:\/\/www.esri.com\/arcgis-blog\/products\/product\/uncategorized\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\/"},"modified":"2021-08-03T00:24:47","modified_gmt":"2021-08-03T07:24:47","slug":"new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","status":"publish","type":"blog","link":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","title":{"rendered":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons"},"author":7471,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"open","ping_status":"closed","template":"","format":"standard","meta":{"_acf_changed":false,"_searchwp_excluded":""},"categories":[23341],"tags":[35581,35591,35601,23391,25581],"industry":[],"product":[36561],"class_list":["post-82071","blog","type-blog","status-publish","format-standard","hentry","category-analytics","tag-ebk","tag-empirical-bayesian-kriging","tag-geostatistical-analyst","tag-spatial-analytics","tag-spatial-statistics","product-arcgis-pro"],"acf":{"short_description":"A new\u00a0subsetting algorithm has been developed in Geostatistical Analyst for ArcGIS Pro 2.1, Generate Subset Polygons, as a geoprocessin...","flexible_content":[{"acf_fc_layout":"content","content":"<h1><\/h1>\n<p>A new\u00a0subsetting algorithm has been developed in Geostatistical Analyst for ArcGIS Pro 2.1, <a href=\"http:\/\/pro.arcgis.com\/en\/pro-app\/tool-reference\/geostatistical-analyst\/generate-subset-polygons.htm\">Generate Subset Polygons<\/a>, as a geoprocessing tool.\u00a0The purpose of this tool is to break down the <a href=\"https:\/\/www.esri.com\/en-us\/arcgis\/products\/spatial-analytics-data-science\/overview\">spatial data<\/a> into small, nonoverlapping subsets.<!--more--><\/p>\n<p>The new tool is intended to be used to create <em>Subset Polygons <\/em>in <a href=\"http:\/\/pro.arcgis.com\/en\/pro-app\/tool-reference\/geostatistical-analyst\/ebk-regression-prediction.htm\">EBK Regression Prediction<\/a> and any future tools that allow you to define subsets using polygons.<\/p>\n<p><strong>Why do we need a new subsetting algorithm?<\/strong><\/p>\n<p>The current subsetting algorithm implemented in <a href=\"http:\/\/pro.arcgis.com\/en\/pro-app\/tool-reference\/geostatistical-analyst\/empirical-bayesian-kriging.htm\">Empirical Bayesian Kriging<\/a> (EBK) and EBK Regression Prediction often encounters problems with clustered data. For example, in the figure below, the current subsetting algorithm often combines data far away from each other into the same subset. This is not desirable because the subsets should be as compact as possible.<\/p>\n<p><figure id=\"attachment_97937\" aria-describedby=\"caption-attachment-97937\" style=\"width: 300px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig21.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-97937\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig21-300x253.png\" alt=\"Figure 1: Overlapping subsetting polygon using EBK Regression Prediction, where the blue dots are the point location of rainfall stations and the selected polygon (cyan) shows the overlapping nature. The data for this analysis have been taken from [1].\" width=\"300\" height=\"253\" \/><\/a><figcaption id=\"caption-attachment-97937\" class=\"wp-caption-text\"><em>Figure 1:<\/em> Overlapping subsetting polygon using EBK Regression Prediction, where the blue dots are the point location of rainfall stations and the selected polygon (cyan) shows the overlapping nature. The data for this analysis have been taken from [1].<\/figcaption><\/figure><strong>The new algorithm<\/strong><\/p>\n<p>In the new algorithm, for each subset\u00a0<em>S<\/em><sub><span style=\"font-size: small;\">i<\/span><\/sub>, the number of points <em>n<\/em><sub>i<\/sub>\u00a0satisfies the constraint <em>min\u00a0\u2264 <em>n<\/em><sub>i<\/sub>\u00a0\u2264max<\/em>\u00a0and minimizes the sum of the pairwise squared deviations within the subsets \u2211<em><sub>i<\/sub><\/em>\u00a0\u2211<sub><em>x<\/em>\u2208<em>S<\/em><sub>i<\/sub><\/sub>\u00a0\u2211<sub><em>y<\/em>\u2208<em>S<\/em><sub>i<\/sub><\/sub>\u00a0||<em>x-y<\/em>||<sup>2<\/sup>. This could also be reorganized as the sum of weighted variances within the subsets 2\u2211<em><sub>i\u00a0<\/sub><\/em><em>n<\/em><sub>i\u00a0<\/sub>\u2211<sub><em>x<\/em>\u2208<em>S<\/em><sub>i<\/sub><\/sub> ||<em>x<\/em>&#8211; <em>c<\/em><sub>i<\/sub>||<sup>2<\/sup>, where <em>c<\/em><sub>i<\/sub>\u00a0is the center of the subset <em>S<\/em><sub>i<\/sub>. Note, this algorithm has a harsher penalty on the number of points in each subset as compared to the K-mean clustering.<\/p>\n<p>The new algorithm performs the following three steps:<\/p>\n<p>Step 1: Connects all points and form a closed curve.<\/p>\n<p>Step 2: Cuts the curve into subsets.<\/p>\n<p>Step 3: Finds all overlaps and resolves them.<\/p>\n<p><strong><em>Step 1<\/em><\/strong><\/p>\n<p>The first step connects all points and forms a closed curve. To achieve this, we form a grid in the work space and put all points that are regarded as clusters into\u00a0corresponding\u00a0grid cells.<\/p>\n<p><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Zeren.gif\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-full wp-image-97954\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Zeren.gif\" alt=\"\" width=\"600\" height=\"338\" \/><\/a><\/p>\n<p><em>Figure 2:<\/em> Example of how the points are connected to form a closed curve.<\/p>\n<p>If two clusters are in the same cell, we merge them into one. A hash table is used to locate the clusters in the grid for efficient searching. We merge nearby clusters to form the curve. Initially, the cell size is extremely fine to ensure that points are connected locally. As the size of the cluster increases, the cell size also increases. The time complexity of this step is <em>O(n \u00b7 log n)<\/em>, where <em>n<\/em>\u00a0is the total number of points.<\/p>\n<p><strong><em>Step 2<\/em><\/strong><\/p>\n<p>The second step cuts the curve into subsets. While cutting the curve, we satisfy the constraint on the number of points in each subset and minimize the sum of the pairwise squared deviations within the subsets. A dynamic programming algorithm is applied to solve this task. The complexity of the algorithm is <em>O((max-min) \u00b7 n)<\/em>.<\/p>\n<figure id=\"attachment_97712\" aria-describedby=\"caption-attachment-97712\" style=\"width: 300px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/part_2_result.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-97712\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/part_2_result-300x225.png\" alt=\"\" width=\"300\" height=\"225\" \/><\/a><figcaption id=\"caption-attachment-97712\" class=\"wp-caption-text\"><em>Figure 3:<\/em> The curve is cut into subsets, where the number of points in each subset satisfies the minimum and maximum requirement, and the sum of the pairwise squared deviations within the subsets are minimized.<\/figcaption><\/figure>\n<p><strong><em>Step 3<\/em><\/strong><\/p>\n<p>The third step finds all overlaps and resolves them while further minimizing the sum of the pairwise squared deviations within the subsets. At present, we are using a brute force search in the third step, which accounts for most of the total execution time (with complexity\u00a0<em>O(n<sup>2<\/sup>)<\/em>), to find overlapping subset pairs. The overlapping is resolved by finding a partitioning between the two subsets, minimizing the sum of the pairwise squared deviations within each subset, while maintaining the required minimum and maximum number of points in each subset. This is performed by projecting points to a set of directions in n-dimension covering the space evenly. For each projection, the optimal division can be found by a dynamic programming approach. Among all projections, the one with the minimum penalty is chosen. The algorithm iterates until no overlaps are detected or the number of iterations reaches the upper bound. The figure below shows a single overlap being removed.<\/p>\n<figure id=\"attachment_97714\" aria-describedby=\"caption-attachment-97714\" style=\"width: 640px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig5.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-large wp-image-97714\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig5-1024x660.png\" alt=\"Figure 5: Finding and resolving overlaps within the subsets.\" width=\"640\" height=\"412\" \/><\/a><figcaption id=\"caption-attachment-97714\" class=\"wp-caption-text\"><em>Figure 4:<\/em> Finding and resolving overlaps within the subsets.<\/figcaption><\/figure>\n<p>The new subsetting algorithm classifies each point into a subset (Figure 5) and creates polygons that wrap around each individual subset (Figure 6).\u00a0 Thus, for each polygon, all points inside the polygon belong to the same subset.<\/p>\n<figure id=\"attachment_97935\" aria-describedby=\"caption-attachment-97935\" style=\"width: 300px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig6.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-97935\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig6-300x240.png\" alt=\"Figure 6: The new subsetting algorithm produces non-overlapping subset.\" width=\"300\" height=\"240\" \/><\/a><figcaption id=\"caption-attachment-97935\" class=\"wp-caption-text\"><em>Figure 5:<\/em> The new subsetting algorithm produces non-overlapping subset.<\/figcaption><\/figure>\n<figure id=\"attachment_97936\" aria-describedby=\"caption-attachment-97936\" style=\"width: 300px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig7.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-97936\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Fig7-300x253.png\" alt=\"Figure 7: The non-overlapping subsets produces the final tool output which are represented by non-overlapping polygons.\" width=\"300\" height=\"253\" \/><\/a><figcaption id=\"caption-attachment-97936\" class=\"wp-caption-text\"><em>Figure 6:<\/em> The non-overlapping subsets produces the final tool output which are represented by non-overlapping polygons.<\/figcaption><\/figure>\n<p><strong>Performance analysis of the new algorithm<\/strong><\/p>\n<p>Overall, the time complexity of the whole algorithm is <em>O(n<sup>2<\/sup>)<\/em>, given the high computational complexity of searching for overlapped subsets.<\/p>\n<p>The figure below shows the time complexity of the three steps as well as each function plotted on the same graph.\u00a0 In the current implementation, the third step takes the vast majority of the computation time for large numbers of points.<\/p>\n<figure id=\"attachment_97721\" aria-describedby=\"caption-attachment-97721\" style=\"width: 640px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Time.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-large wp-image-97721\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/Time-1024x775.png\" alt=\"Figure 8: Time complexity of the algorithm in the three stages.\" width=\"640\" height=\"484\" \/><\/a><figcaption id=\"caption-attachment-97721\" class=\"wp-caption-text\"><em>Figure 7:<\/em> Time complexity of the algorithm in the three stages.<\/figcaption><\/figure>\n<p>The space complexity of the algorithm is <em>O(n)<\/em>, which means that the required memory is proportional to the number of points.\u00a0The following image shows the memory allocation is linear with the number of input points.<\/p>\n<figure id=\"attachment_97722\" aria-describedby=\"caption-attachment-97722\" style=\"width: 300px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/memory.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-medium wp-image-97722\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/memory-300x225.png\" alt=\"Figure 9: Space complexity of the algorithm.\" width=\"300\" height=\"225\" \/><\/a><figcaption id=\"caption-attachment-97722\" class=\"wp-caption-text\"><em>Figure 8:<\/em> Space complexity of the algorithm.<\/figcaption><\/figure>\n<p>The following graphs show the analysis of the algorithm in the third step (resolution of overlapping subsets). When two subsets overlap, resolving their overlap reduces the total sum of the pairwise squared deviations. The graph on the left shows the decrease in the total sum after each resolution. The graphs on the right shows the number of resolutions in each iteration. For this data, the algorithm converges after about 14 iterations.<\/p>\n<figure id=\"attachment_97723\" aria-describedby=\"caption-attachment-97723\" style=\"width: 640px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/sum_of_err.png\"><img loading=\"lazy\" decoding=\"async\" class=\"size-large wp-image-97723\" src=\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/12\/sum_of_err-1024x576.png\" alt=\"Figure 10: Analysis of the algorithm in the third step where overlapping subsets are resolved.\" width=\"640\" height=\"360\" \/><\/a><figcaption id=\"caption-attachment-97723\" class=\"wp-caption-text\"><em>Figure 9:<\/em> Analysis of the algorithm in the third step where overlapping subsets are resolved.<\/figcaption><\/figure>\n<p><strong>Application<\/strong><\/p>\n<p>This new subsetting algorithm can work efficiently with clustered datasets with more than a billion points. The quality of the constructed subsets in the new subsetting algorithm is often better than the algorithm currently used in EBK. Though it does not yet support 3D points, the methodology can be easily extended to multiple dimensions.<\/p>\n<p><em>Part of this new subsetting algorithm and blog content was created by Zeren Shui, advised by Alexander Gribov, during his internship with the Geostatistical Analyst Team in summer 2017. Zeren is currently pursuing his Master\u2019s degree in Data Science at College of Science and Engineering, University of Minnesota \u2013 Twin Cities. His research interests are Bayesian Statistics, Machine Learning, and Data Mining. For additional questions, you can comment here or contact him at <\/em><a href=\"mailto:shuix007@umn.edu\"><em>shuix007@umn.edu<\/em><\/a><em>.<\/em><\/p>\n<p><strong>References<\/strong><\/p>\n<p>[1] S. D. Lynch, <a href=\"http:\/\/www.wrc.org.za\/Knowledge%20Hub%20Documents\/Research%20Reports\/1156-1-04.pdf\">Development of a raster database of annual, monthly and daily rainfall for southern Africa<\/a>, WRC Report (1156\/1\/04) (2004) 78.<\/p>\n"}],"authors":[{"ID":7471,"user_firstname":"Alexander","user_lastname":"Gribov","nickname":"Alexander Gribov","user_nicename":"alexander-gribov","display_name":"Alexander Gribov","user_email":"agribov@esri.com","user_url":"","user_registered":"2018-03-21 18:21:20","user_description":"","user_avatar":"<img alt='' src='https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=96&#038;d=blank&#038;r=g' srcset='https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=192&#038;d=blank&#038;r=g 2x' class='avatar avatar-96 photo' height='96' width='96' loading='lazy' decoding='async'\/>"}],"related_articles":"","card_image":false,"wide_image":false},"yoast_head":"<!-- This site is optimized with the Yoast SEO Premium plugin v25.9 (Yoast SEO v25.9) - https:\/\/yoast.com\/wordpress\/plugins\/seo\/ -->\n<title>New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons<\/title>\n<meta name=\"robots\" content=\"index, follow, max-snippet:-1, max-image-preview:large, max-video-preview:-1\" \/>\n<link rel=\"canonical\" href=\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\" \/>\n<meta property=\"og:locale\" content=\"en_US\" \/>\n<meta property=\"og:type\" content=\"article\" \/>\n<meta property=\"og:title\" content=\"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons\" \/>\n<meta property=\"og:url\" content=\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\" \/>\n<meta property=\"og:site_name\" content=\"ArcGIS Blog\" \/>\n<meta property=\"article:publisher\" content=\"https:\/\/www.facebook.com\/esrigis\/\" \/>\n<meta property=\"article:modified_time\" content=\"2021-08-03T07:24:47+00:00\" \/>\n<meta name=\"twitter:card\" content=\"summary_large_image\" \/>\n<meta name=\"twitter:site\" content=\"@ESRI\" \/>\n<script type=\"application\/ld+json\" class=\"yoast-schema-graph\">{\"@context\":\"https:\/\/schema.org\",\"@graph\":[{\"@type\":[\"Article\",\"BlogPosting\"],\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#article\",\"isPartOf\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\"},\"author\":{\"name\":\"Alexander Gribov\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/d2b862be484b092ad545bd5b9fb44b81\"},\"headline\":\"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons\",\"datePublished\":\"2018-01-22T17:00:27+00:00\",\"dateModified\":\"2021-08-03T07:24:47+00:00\",\"mainEntityOfPage\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\"},\"wordCount\":9,\"commentCount\":0,\"publisher\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#organization\"},\"keywords\":[\"EBK\",\"Empirical Bayesian Kriging\",\"Geostatistical Analyst\",\"spatial analytics\",\"spatial statistics\"],\"articleSection\":[\"Analytics\"],\"inLanguage\":\"en-US\",\"potentialAction\":[{\"@type\":\"CommentAction\",\"name\":\"Comment\",\"target\":[\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#respond\"]}]},{\"@type\":\"WebPage\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\",\"url\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\",\"name\":\"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons\",\"isPartOf\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#website\"},\"datePublished\":\"2018-01-22T17:00:27+00:00\",\"dateModified\":\"2021-08-03T07:24:47+00:00\",\"breadcrumb\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#breadcrumb\"},\"inLanguage\":\"en-US\",\"potentialAction\":[{\"@type\":\"ReadAction\",\"target\":[\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons\"]}]},{\"@type\":\"BreadcrumbList\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#breadcrumb\",\"itemListElement\":[{\"@type\":\"ListItem\",\"position\":1,\"name\":\"Home\",\"item\":\"https:\/\/www.esri.com\/arcgis-blog\/\"},{\"@type\":\"ListItem\",\"position\":2,\"name\":\"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons\"}]},{\"@type\":\"WebSite\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#website\",\"url\":\"https:\/\/www.esri.com\/arcgis-blog\/\",\"name\":\"ArcGIS Blog\",\"description\":\"Get insider info from Esri product teams\",\"publisher\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#organization\"},\"potentialAction\":[{\"@type\":\"SearchAction\",\"target\":{\"@type\":\"EntryPoint\",\"urlTemplate\":\"https:\/\/www.esri.com\/arcgis-blog\/?s={search_term_string}\"},\"query-input\":{\"@type\":\"PropertyValueSpecification\",\"valueRequired\":true,\"valueName\":\"search_term_string\"}}],\"inLanguage\":\"en-US\"},{\"@type\":\"Organization\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#organization\",\"name\":\"Esri\",\"url\":\"https:\/\/www.esri.com\/arcgis-blog\/\",\"logo\":{\"@type\":\"ImageObject\",\"inLanguage\":\"en-US\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/logo\/image\/\",\"url\":\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/04\/Esri.png\",\"contentUrl\":\"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/04\/Esri.png\",\"width\":400,\"height\":400,\"caption\":\"Esri\"},\"image\":{\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/logo\/image\/\"},\"sameAs\":[\"https:\/\/www.facebook.com\/esrigis\/\",\"https:\/\/x.com\/ESRI\",\"https:\/\/www.linkedin.com\/company\/5311\/\"]},{\"@type\":\"Person\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/d2b862be484b092ad545bd5b9fb44b81\",\"name\":\"Alexander Gribov\",\"image\":{\"@type\":\"ImageObject\",\"inLanguage\":\"en-US\",\"@id\":\"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/image\/\",\"url\":\"https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=96&d=blank&r=g\",\"contentUrl\":\"https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=96&d=blank&r=g\",\"caption\":\"Alexander Gribov\"},\"url\":\"\"}]}<\/script>\n<!-- \/ Yoast SEO Premium plugin. -->","yoast_head_json":{"title":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons","robots":{"index":"index","follow":"follow","max-snippet":"max-snippet:-1","max-image-preview":"max-image-preview:large","max-video-preview":"max-video-preview:-1"},"canonical":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","og_locale":"en_US","og_type":"article","og_title":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons","og_url":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","og_site_name":"ArcGIS Blog","article_publisher":"https:\/\/www.facebook.com\/esrigis\/","article_modified_time":"2021-08-03T07:24:47+00:00","twitter_card":"summary_large_image","twitter_site":"@ESRI","schema":{"@context":"https:\/\/schema.org","@graph":[{"@type":["Article","BlogPosting"],"@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#article","isPartOf":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons"},"author":{"name":"Alexander Gribov","@id":"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/d2b862be484b092ad545bd5b9fb44b81"},"headline":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons","datePublished":"2018-01-22T17:00:27+00:00","dateModified":"2021-08-03T07:24:47+00:00","mainEntityOfPage":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons"},"wordCount":9,"commentCount":0,"publisher":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/#organization"},"keywords":["EBK","Empirical Bayesian Kriging","Geostatistical Analyst","spatial analytics","spatial statistics"],"articleSection":["Analytics"],"inLanguage":"en-US","potentialAction":[{"@type":"CommentAction","name":"Comment","target":["https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#respond"]}]},{"@type":"WebPage","@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","url":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons","name":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons","isPartOf":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/#website"},"datePublished":"2018-01-22T17:00:27+00:00","dateModified":"2021-08-03T07:24:47+00:00","breadcrumb":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#breadcrumb"},"inLanguage":"en-US","potentialAction":[{"@type":"ReadAction","target":["https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons"]}]},{"@type":"BreadcrumbList","@id":"https:\/\/www.esri.com\/arcgis-blog\/products\/analytics\/analytics\/new-subsetting-tool-in-geostatistical-analyst-generate-subset-polygons#breadcrumb","itemListElement":[{"@type":"ListItem","position":1,"name":"Home","item":"https:\/\/www.esri.com\/arcgis-blog\/"},{"@type":"ListItem","position":2,"name":"New Subsetting Tool in Geostatistical Analyst: Generate Subset Polygons"}]},{"@type":"WebSite","@id":"https:\/\/www.esri.com\/arcgis-blog\/#website","url":"https:\/\/www.esri.com\/arcgis-blog\/","name":"ArcGIS Blog","description":"Get insider info from Esri product teams","publisher":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/#organization"},"potentialAction":[{"@type":"SearchAction","target":{"@type":"EntryPoint","urlTemplate":"https:\/\/www.esri.com\/arcgis-blog\/?s={search_term_string}"},"query-input":{"@type":"PropertyValueSpecification","valueRequired":true,"valueName":"search_term_string"}}],"inLanguage":"en-US"},{"@type":"Organization","@id":"https:\/\/www.esri.com\/arcgis-blog\/#organization","name":"Esri","url":"https:\/\/www.esri.com\/arcgis-blog\/","logo":{"@type":"ImageObject","inLanguage":"en-US","@id":"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/logo\/image\/","url":"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/04\/Esri.png","contentUrl":"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2018\/04\/Esri.png","width":400,"height":400,"caption":"Esri"},"image":{"@id":"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/logo\/image\/"},"sameAs":["https:\/\/www.facebook.com\/esrigis\/","https:\/\/x.com\/ESRI","https:\/\/www.linkedin.com\/company\/5311\/"]},{"@type":"Person","@id":"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/d2b862be484b092ad545bd5b9fb44b81","name":"Alexander Gribov","image":{"@type":"ImageObject","inLanguage":"en-US","@id":"https:\/\/www.esri.com\/arcgis-blog\/#\/schema\/person\/image\/","url":"https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=96&d=blank&r=g","contentUrl":"https:\/\/secure.gravatar.com\/avatar\/b803c48dfbaf6aded2605f2a4ca72e086357ff7dfc21934f8008c8a972652aa5?s=96&d=blank&r=g","caption":"Alexander Gribov"},"url":""}]}},"text_date":"January 22, 2018","author_name":"Alexander Gribov","author_page":false,"custom_image":"https:\/\/www.esri.com\/arcgis-blog\/app\/uploads\/2025\/08\/Newsroom-Keyart-Wide-1920-x-1080.jpg","primary_product":"ArcGIS Pro","tag_data":[{"term_id":35581,"name":"EBK","slug":"ebk","term_group":0,"term_taxonomy_id":35581,"taxonomy":"post_tag","description":"","parent":0,"count":1,"filter":"raw"},{"term_id":35591,"name":"Empirical Bayesian Kriging","slug":"empirical-bayesian-kriging","term_group":0,"term_taxonomy_id":35591,"taxonomy":"post_tag","description":"","parent":0,"count":3,"filter":"raw"},{"term_id":35601,"name":"Geostatistical Analyst","slug":"geostatistical-analyst","term_group":0,"term_taxonomy_id":35601,"taxonomy":"post_tag","description":"","parent":0,"count":7,"filter":"raw"},{"term_id":23391,"name":"spatial analytics","slug":"spatial-analytics","term_group":0,"term_taxonomy_id":23391,"taxonomy":"post_tag","description":"","parent":0,"count":348,"filter":"raw"},{"term_id":25581,"name":"spatial statistics","slug":"spatial-statistics","term_group":0,"term_taxonomy_id":25581,"taxonomy":"post_tag","description":"","parent":0,"count":129,"filter":"raw"}],"category_data":[{"term_id":23341,"name":"Analytics","slug":"analytics","term_group":0,"term_taxonomy_id":23341,"taxonomy":"category","description":"","parent":0,"count":1363,"filter":"raw"}],"product_data":[{"term_id":36561,"name":"ArcGIS Pro","slug":"arcgis-pro","term_group":0,"term_taxonomy_id":36561,"taxonomy":"product","description":"","parent":0,"count":2122,"filter":"raw"}],"primary_product_link":"https:\/\/www.esri.com\/arcgis-blog\/?s=#&products=arcgis-pro","_links":{"self":[{"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/blog\/82071","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/blog"}],"about":[{"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/types\/blog"}],"author":[{"embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/users\/7471"}],"replies":[{"embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/comments?post=82071"}],"version-history":[{"count":0,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/blog\/82071\/revisions"}],"wp:attachment":[{"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/media?parent=82071"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/categories?post=82071"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/tags?post=82071"},{"taxonomy":"industry","embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/industry?post=82071"},{"taxonomy":"product","embeddable":true,"href":"https:\/\/www.esri.com\/arcgis-blog\/wp-json\/wp\/v2\/product?post=82071"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}