cookwest_georeferencepoint
Cook West Georeference Points
2019-09-24
2019-09-24
2019-09-25
David Huggins
USDA-ARS NSAR
Soil Scientist
<b>Methods</b><div><br /></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Produced a 30.48 by 30.48 meter grid running in a true north south orientation and superimposed over the <a href='https://usdaars.maps.arcgis.com/home/item.html?id=5f564666af3347cba4ac6b935340ce4c' target='_blank'>Cook West field boundary</a>. The field boundary area was used to clip the grid producing a grid with many cells on the boundary smaller than the 30.48 by 30.48 meters. One point was randomly placed within each full grid cell producing 292 potential sampling points. Points were then manually removed from many of the smaller cells formed along the boundary of the area resulting in 250 points. Points were also randomly moved within grid cells when the original placement overlapped roads, excavations conducted to map tile drains, known tile drain locations or any other infrastructure. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'> </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Sampling locations 250</p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Area = 23.4 hectares = 234000 square meters</p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>10.68 points per hectare, 1 point per 936 meters. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'> </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Randomized points for Cook West generated in ArcGIS Desktop version 10.4.05524:</p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'> </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Create a grid over the area of interest using the Fishnet tool in ArcGIS. The height and width of the grid should be predetermined to produce the desired density of sampling locations per area if one sampling point is randomly paced within each grid. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Convert the boundary for the area of interest into a polygon feature if it is a polyline feature. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Clip the grid using the boundary polygon feature of the area of interest, using the Clip (Analysis) tool. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>Convert the clipped polyline grid into a polygon feature using the Feature to Polygon tool. Use both the clipped poly line grid and the polyline boarder for the same area of interest. The result should be a polygon feature with multiple polygons representing the clipped grid. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'>There are a number of ways to generate random points. For this a custom ArcToolbox script was used "Generate Random Points" http://ianbroad.com/ , the built in Create Random Points tool can work as well. </p><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'> </p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>#----------------------------------------------------------------------</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'># Purpose: Generates a user defined number of random points inside</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'># each polygon of an input Polygon feature class.</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>#</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'># Author: Ian Broad</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'># Website: www.ianbroad.com</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>#</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'># Created: 02/27/2015</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>#----------------------------------------------------------------------</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>import arcpy</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>import random</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>import sys</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>arcpy.env.overwriteOutput = True</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>polygon = arcpy.GetParameterAsText(0)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>use_field = arcpy.GetParameterAsText(1)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>field = arcpy.GetParameterAsText(2)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>number_of_points = abs(int(arcpy.GetParameterAsText(3)))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>min_distance = abs(float(arcpy.GetParameterAsText(4)))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>attempts = abs(int(arcpy.GetParameterAsText(5)))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>keep_attempts = arcpy.GetParameterAsText(6)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>output = arcpy.GetParameterAsText(7)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>if "in_memory" in output:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> mem_name = output.split("\\")[-1]</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> mem_name = "mem_point"</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>mem_point = arcpy.CreateFeatureclass_management("in_memory", mem_name, "POINT", "", "DISABLED", "DISABLED", polygon)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>arcpy.AddField_management(mem_point, "PolygonOID", "TEXT")</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>result = arcpy.GetCount_management(polygon)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>features = int(result.getOutput(0))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>arcpy.SetProgressor("step", "Generating Random Points...", 0, features, 1)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>fields = ["SHAPE@", "OID@"]</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>if use_field == "YES":</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> fields.append(field)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>with arcpy.da.InsertCursor(mem_point, ("SHAPE@", "PolygonOID")) as insert:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> with arcpy.da.SearchCursor(polygon, (fields)) as search:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> for row in search:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> try:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> polygon_geom = row[0]</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> oid = row[1]</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count = 0</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> attempt = 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance = []</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> xmin, xmax= row[0].extent.XMin, row[0].extent.XMax</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> ymin, ymax = row[0].extent.YMin, row[0].extent.YMax</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if use_field == "YES":</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> number_of_points = int(row[2])</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> while count < number_of_points:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> random_x = random.uniform(xmin, xmax)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> random_y = random.uniform(ymin, ymax)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> point = arcpy.Point(random_x, random_y)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> point_geom = arcpy.PointGeometry(point, polygon_geom.spatialReference)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> contains_point = polygon_geom.contains(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if keep_attempts == "NO":</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if attempt < attempts:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if contains_point == True and count == 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> insert.insertRow((point_geom, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> elif contains_point == True and count > 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = True</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> for point in distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if point_geom.distanceTo(point) > min_distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> pass</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = False</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> attempt += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if distance_check == True:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> insert.insertRow((point_geom, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddError("Failed attempting to generate {0} random points for Polygon OID: {1}".format(attempts, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddError("Decrease number of random points or the minimum distance and try again.")</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> sys.exit(0)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if keep_attempts == "KEEP_ALL":</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if attempt < attempts:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if contains_point == True and count == 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> insert.insertRow((point_geom, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> elif contains_point == True and count > 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = True</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> for point in distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if point_geom.distanceTo(point) > min_distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> pass</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = False</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> attempt += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if distance_check == True:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> insert.insertRow((point_geom, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddWarning("Only generated {0} out of {1} random points for Polygon OID: {2}".format(count, number_of_points, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddWarning("Decrease number of random points or the minimum distance and try again.")</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> break</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if keep_attempts == "KEEP_ONLY_SUCCESSFUL":</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if attempt < attempts:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if contains_point == True and count == 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> elif contains_point == True and count > 0:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = True</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> for point in distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if point_geom.distanceTo(point) > min_distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> pass</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance_check = False</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> attempt += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if distance_check == True:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> distance.append(point_geom)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> count += 1</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> if int(len(distance)) == int(number_of_points):</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> for geom in distance:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> insert.insertRow((geom, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddWarning("Failed attempting to generate {0} random points for Polygon OID: {1}".format(number_of_points, oid))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddWarning("Decrease number of random points or the minimum distance and try again.")</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> break</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.SetProgressorPosition()</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> except Exception as e:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.AddMessage(str(e.message))</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>if "in_memory" in output:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.SetParameter(8, mem_point)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>else:</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.CopyFeatures_management(mem_point, output)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'> arcpy.Delete_management(mem_point)</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'><br /></font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>arcpy.ResetProgressor()</font></p></span></div><div><span><p style='margin-top: 0pt; margin-bottom: 0pt; line-height: 1.38;'><font size='3'>arcpy.GetMessages()</font></p></span></div><div><span><div><br /></div><div><span style='font-size: 12pt; font-family: "Times New Roman"; color: rgb(0, 0, 0); background-color: transparent; font-variant-numeric: normal; font-variant-east-asian: normal; vertical-align: baseline; white-space: pre-wrap;'><br /></span></div><div><span style='font-size: 12pt; font-family: "Times New Roman"; color: rgb(0, 0, 0); background-color: transparent; font-variant-numeric: normal; font-variant-east-asian: normal; vertical-align: baseline; white-space: pre-wrap;'><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><b>Processing</b></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><br /></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'>Methods for data processing from original files are found on Github: <a href='https://github.com/cafltar/CafGeospatialBoundariesAndGridPoints_CleanAndStandardize' style='color: rgb(0, 121, 193); text-decoration-line: none; font-family: inherit;'>https://github.com/cafltar/CafGeospatialBoundariesAndGridPoints_CleanAndStandardize</a></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><br /></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><b>Data Dictionary</b></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><br /></div><div style='font-family: "Avenir Next W01", "Avenir Next W00", "Avenir Next", Avenir, "Helvetica Neue", sans-serif; color: rgb(76, 76, 76); white-space: normal;'><table border='0' cellpadding='0' cellspacing='0' style='border-collapse: collapse; width: 192pt;' width='256'><tbody><tr height='19' style='height: 14.4pt;'><td height='19' style='height: 14.4pt; width: 48pt;' width='64'>FieldName</td><td style='width: 48pt;' width='64'>Units</td><td style='width: 48pt;' width='64'>Description</td><td style='width: 48pt;' width='64'>DataType</td></tr><tr height='19' style='height: 14.4pt;'><td height='19' style='height: 14.4pt;'>ID2</td><td>unitless</td><td>Numeric value used to identify georeferenced points for long-term sample collection. Values are unique among both Cook fields (CE and CW). Use 'ID2' instead of 'Id' for historic reasons.</td><td>Int</td></tr></tbody></table></div></span></div></span></div>
Georeference points for standard sampling locations for a long-term experiment. These points are used to collect data for the "business as usual" treatment of the LTAR Common Experiment.
LTARCAFLTAR
1
-117.1034
-117.053
46.7733
46.7994
2017-01-01T00:00:00
now
<p style='margin-top: 0px; margin-bottom: 16px; font-size: 16px; box-sizing: border-box; color: rgb(36, 41, 46); font-family: -apple-system, BlinkMacSystemFont, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji";'>As a work of the United States government, this project is in the public domain within the United States.</p><p style='margin-top: 0px; font-size: 16px; box-sizing: border-box; color: rgb(36, 41, 46); font-family: -apple-system, BlinkMacSystemFont, "Segoe UI", Helvetica, Arial, sans-serif, "Apple Color Emoji", "Segoe UI Emoji"; margin-bottom: 0px !important;'>Additionally, we waive copyright and related rights in the work worldwide through the CC0 1.0 Universal public domain dedication.</p>
North American Profile of ISO19115 2003
2019-09-24
11:49:20.89
2019-09-25
13:09:23.06
1.0
NAP
editor:esri.dijit.metadata.editor
2019-09-24
b0601d346e6a4cce83ca0c9544128fc0
Bryan Carlson
USDA-ARS NSAR
Ecoinformaticst