Latest Contributions by Jay_Gregory

I'm trying to implement a fairly straightforward multiprocessing script that takes every polygon in a shapefile, tessellates into a grid, takes the center points of each grid square, and then clips the points by the original feature (since the tessellation generates a minimum bounding rectangle).  The tessellation feature class is never created so next step (FeatureToPoint) fails.  I don't understand - can someone help? Even if I comment out all the arcpy commands except GenerateTessellation, the output is never actually created.     import multiprocessing import os import time import arcpy startTime = time.time() basedir = r"C:\Users\Jay\Documents\ArcGIS\Projects\DataPrep" scratchGDB = 'C:\\Jay\\data\data.gdb' output = os.path.join(basedir, "grids.gdb") features = os.path.join(basedir,"features.shp") def multifunction(ft): import arcpy name = ft[0] extent = ft[1].extent print('Working on {}'.format(name)) fl = arcpy.management.MakeFeatureLayer(features, "{}fl".format(name), where_clause=f"Name='{name}'") gdb = os.path.join(basedir, "{}.gdb".format(name)) arcpy.management.CreateFileGDB(basedir, "{}.gdb".format(name)) output=os.path.join(gdb, "{}grid".format(name)) grid = arcpy.management.GenerateTessellation(output, Extent=extent, Shape_Type="SQUARE", Size="900 SquareMeters") print('tellessation') pointgrid = arcpy.management.FeatureToPoint(grid, os.path.join(scratchGDB, "{}pointgrid".format(apt))) arcpy.analysis.Clip(pointgrid, fl, os.path.join(output, f"{apt}Points")) arcpy.management.Delete(fl) arcpy.management.Delete(grid) arcpy.management.Delete(pointgrid) def main(): processList = [feature for feature in arcpy.da.SearchCursor(features, ['Name', 'SHAPE@'])] pool = multiprocessing.Pool(1) pool.map(multifunction, processList[0:10]) pool.close() pool.join() if __name__ == "__main__": print("Running") main() executionTime = (time.time() - startTime) print('Execution time in seconds: ' + str(executionTime))   ... View more

I have an S3 bucket with over 2000 geotiffs containing elevation data (I believe it is 1arcsec SRTM), and would like to publish an elevation service to my ArcGIS Enterprise, but am unclear on the process.   Do I need to mosaic all the geotiffs into one dataset and then just publish an image service, or can I just load all 2000 images into Pro and publish as is?  Is there any other data prep I would need to do?   I know I would need to register the S3 bucket as a cloud store to reference the images in place, but unsure on other steps and best practices.  Thanks! ... View more

I would like advice on how to visualize my project in 3D.   I have 12 datasets which are the result of a model, each dataset represents the model output at an elevation "slice" (i.e. dataset 1 is for 0 AGL, dataset 2 is for 10 m AGL, dataset 3 is for 50 m AGL, etc. etc.). Each dataset was generated by running a model against the same 123k features (the only difference between each model run is the elevation).  Each feature represents a 30m x 30m area (so yes it is a lot of data).  So my model output could be 12 point feature classes (by taking the center of each square), 12 polygon feature classes, or 12 raster datasets.  But what I have is 12 distinct products that I can only visualize one at a time.   How can I combine these into a 3D visualization - first to visualize in ArcGIS Pro, and then if possible publish to ArcGIS Enterprise as a Scene Service or some other mechanism.  Any advice would be helpful.     ... View more