Latest Contributions by MatthewPoppleton

I completed an analysis using EBK with the "Power" semivariogram model and the "K-Bessel Detrended" semivariogram model.  The equation for power (γ(h)= Nugget + b|h|α) can be found in "What is empirical Bayesian kriging?", but the K-Bessel equation is not included. If understand correctly, the detrended model should would be the same equation as K-Bessel but with a first-order trend removal. Can anyone provide the equation for the K-Bessel semivariogram model, or suggest where I might be able to find this referenced? I'm not sure if K-Bessel would be an equation I could can reference elsewhere in literature or if it would be specific to this version of ArcGIS Pro EBK. I did find an answer to a similar question by @EricKrause from 2012 that references ArcGIS 9 Documentation, and includes a K-Bessel equation on page 257. Would this be the same model? ... View more

I am using Areal Interpolation to estimate population size using simulated and real data in the same study area. I've already completed the estimates with the simulated data without issues but have run into an error while using the real data. The tool "Areal Interpolation Layer to Polygons" interrupts with the error "ERROR 040065: Grid cell size is too small." I'm using the same input polygon features as I did with my simulated data, so I figure this must be a problem with my input layer. I've googled this error but have not found any related information - could someone clarify this error for me and what can be done to run the tool? ... View more

I am working on a simulation to compare different models' estimates of population counts in areas with known total observations. I am using a preset amount of points (e.g. 50,000) to represent observations within my study area, and run models using line transects that overlap with those points. I currently have a spatially balanced point feature class that acts as a uniform distribution, and another with point clustering based on land-cover type (stationary). I am trying to figure out a way to generate a third point feature class of 50,000 non-stationary randomly clustered points. What I've tried so far is using the "Generate Random Points" tool to generate 50,000 points within my boundary, and set a low minimum allowed distance (e.g. 1 meter). Then I used the "Density-based clustering" tool to detect areas of clustered points and copied them to a new feature class. This leaves me with the clustered areas and no noise points. However with the noise subtracted I am left with far less than the desired amount of points. The output looks good but I would like to set the total number of points to be clustered. See screenshots for reference. I could use a higher count of random points and use trial and error to end up with a clustered set around 50,000 but this wouldn't be ideal as I wouldn't have an exact count and it probably wouldn't be very replicable. I assume this has been done before but I haven't been able to find how. Is there a way to achieve this using other geoprocessing tools or perhaps arcpy?  ... View more

I am currently building a simulation using arcpy in ArcGIS Pro to analyze randomly generated points that intersect with preset line transects. The arcpy.CreateRandomPoints_management() function works as expected, but when I  run the arcpy.geoanalytics.OverlayLayers() function in python, the function fails and I get the error "Unable to get geoanalytics URL". However, when I run overlay layers using the geoprocessing toolbox using the exact same parameters I do not get any errors. I tried to google this error but have not found any related information. I have provided the code below: import arcpy randomName = "pythonRandomPoints_test" overlayName = "layer_test" number_of_points_or_field = 60000 arcpy.CreateRandomPoints_management("C:/Users/User/Documents/ArcGIS/Projects/Default.gdb",    randomName,    "shp_file",    "",    number_of_points_or_field,    "",   "POINT") try:    arcpy.geoanalytics.OverlayLayers('transects_t',       randomName,       overlayName,       "INTERSECT") except arcpy.ExecuteError:    print(arcpy.GetMessages()) #Start Time: Wednesday, July 29, 2020 14:28:09 #Running script Overlay Layers... #Unable to get geoanalytics URL. #Completed script Overlay Layers... #Failed to execute (OverlayLayers). #Failed at Wednesday, July 29, 2020 14:28:10 (Elapsed Time: 0.87 seconds) ... View more