Latest Contributions by JohannesLindner

Sorry, forgot to mention... I took your sample data for the expression. In your real expression, you should of course load the featureset from your portal. So instead of lines 2&3 in my code above, you would do something like this: var fs = FeaturesetByPortalItem( Portal("https://my-server.de/portal"), // your Portal's url "bb11a82308a5439291c17a477755d27a", // the service id 0, // the sub-layer id ["*"], // the fields you want to load false // do you want to load geometries )   You can get all the neccessary arguments from your table's url:   ... View more

Your first script fails because you forgot the index in the failing line. So you're trying to use the non-existing union method of an array (row) instead of the geometry (row[0]). But the script would fail there anyway, because arcpy.management.CopyFeatures() does not what you think it does. You seem to want to extract the geometries, but this tool copies a feature class from one location to another.   You got the gist of what you need to do: Extract the geometries Iterate over the rows with an UpdateCursor union the current geometry with all extracted geometries write it into the feature class optionally: you probably should delete every row after the first one, so you don't get overlapping features.   zones = "TestPolygons" # layer name or fc path zone_name_field = "TextField" # name of the field that determines the "sameness" only_keep_first = True # only keep the first row of each zone, delete the following rows zone_names = ["A"] # Names you want to union # You can also do it automatically for every zone name: #zone_names = {r[0] for r in arcpy.da.SearchCursor(zones, [zone_name_field])} for zone_name in zone_names: sql = f"{zone_name_field} = '{zone_name}'" # extract the geometries into a list zone_geometries = [r[0] for r in arcpy.da.SearchCursor(zones, ["SHAPE@"], sql)] first_row = True # iterate over the fc with arcpy.da.UpdateCursor(zones, ["SHAPE@"], sql) as cursor: for r in cursor: if first_row: # union the current geometry with every extracted geometry geo = r[0] for other_geo in zone_geometries: geo = geo.union(other_geo) cursor.updateRow([geo]) # delete the following rows if only_keep_first == True first_row = not only_keep_first else: cursor.deleteRow()   Before:   With only_keep_first == True:   With only_keep_first == False:   ... View more

Then you should be able to use this simple script in ArcGIS Pro's Python Window: # edit these variables layer_name = "Layer" output_folder = "N:/.../temp" from pathlib import Path attachment_table = arcpy.Describe(layer_name).catalogPath + "__ATTACH" with arcpy.da.SearchCursor(attachment_table, ["DATA", "ATT_NAME"]) as cursor: for data, name in cursor: print(f"currently exporting {name}") out_file = Path(output_folder, name) out_file.write_bytes(data)   ... View more

Hello I am looking for tutorials for layout from start to finish My thanks in advance The official documentation: Mises en page dans ArcGIS Pro—ArcGIS Pro | Documentation  A tutorial on Youtube: Getting Started with ArcGIS Pro: Building a Layout - YouTube ... View more

There are ways, yes. The easiest way would be to find the author of the map (someone in the U.S. Department of Agriculture) and ask them for the source data. These people don't bite, and the worst they can do is say no 🙂   Extracting the points yourself leads you into quite advanced topics, which are probably way out of scope for a school project. But, for the sake of an example, I'll show one way to do it:   Luckily, the map you found is not in fact a png, but a pdf with vector data. You can recognize that by zooming in. If the content of the map gets pixely, it is a raster image (eg a png, jpg, gif, whatever). If it redraws at the new scale, it's vector data. We can extract vector data from the pdf to load it into ArcGIS and further process it there. You can do that with Incscape, or by converting the pdf to a dxf file. There are many online converters available, they all produce different outputs. In this case, I got usable results with onlineconvertfree.   So, convert the pdf to dxf and load that file into ArcGIS Pro:   The dxf file actually contains multiple layers with the different geometry types. Sadly, the points of your map aren't actually stored as points, but as polygons. So we will have to convert them. Also, the data doesn't have a coordinate system, so we have to do something about that.   Let's take a look at the pdf and see if there is anything in there about the coordinate system. We'll do this in Python.   import PyPDF2 from pprint import pprint # read the first (and only) page of the pdf pdf = "N:/bla/cattle_map.pdf" reader = PyPDF2.PdfFileReader(pdf) page = reader.pages[0] # let's take a look pprint(page) #{'/Contents': IndirectObject(1, 0), # '/MediaBox': [0, 0, 792, 612], # '/Parent': IndirectObject(37, 0), # '/Resources': IndirectObject(28, 0), # '/Type': '/Page', # '/VP': [{'/BBox': [58.20162, 469.49288, 692.41922, 43.1716], # '/Measure': IndirectObject(31, 0), # '/Name': 'US Data Frame', # '/Type': '/Viewport'}, # {'/BBox': [54.0015, 575.9982, 212.4059, 475.19317], # '/Measure': IndirectObject(33, 0), # '/Name': 'AK Data Frame', # '/Type': '/Viewport'}, # {'/BBox': [54.0015, 97.17429, 144.004, 35.97124], # '/Measure': IndirectObject(35, 0), # '/Name': 'HI Data Frame', # '/Type': '/Viewport'}]} # The /VP key looks interesting. It includes the viewports for the continental US; Alaska, and Hawaii. Let's look at the US. The Measure parameter is still hidden. pprint(page["/VP"][0]["/Measure"]) #{'/Bounds': [0, 1, 0, 0, 1, 0, 1, 1], # '/GCS': IndirectObject(30, 0), # '/GPTS': [21.6176, # -118.83703, # 48.90833, # -128.98979, # 49.20632, # -64.25761, # 21.82903, # -74.06214], # '/LPTS': [0, 1, 0, 0, 1, 0, 1, 1], # '/Subtype': '/GEO', # '/Type': '/Measure'} # GCS is the abbreviation of Geographic Coordinate System. That should be what we're after. pprint(page["/VP"][0]["/Measure"]["/GCS"]) #{'/Type': '/PROJCS', # '/WKT': 'PROJCS["USA_Custom_Albers_Equal_Area_Conic_US",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",23.0],UNIT["Statute_Mile",1609.344]]'} # Yep, there's the definition of the coordinate system. Let's store that in a variable. wkt = page["/VP"][0]["/Measure"]["/GCS"]["/WKT"]     Now that we know the coordinate system, lets apply it to the map:   active_map = arcpy.mp.ArcGISProject("current").activeMap active_map.spatialReference = arcpy.SpatialReference(text=wkt)     The coordinates in the pdf are wrong, so we need to georeference the dxw layer.   Now we can extract the polygons. There are probably good reasons for the custom projection, but it might be better to change the projection to something more widely used. I'm using WGS 84 here, which is the standard system used in GPS. We don't need to project every layer, only the "Polygon Group / C1E-64-B3" layer, which contains the points (as polygons).   arcpy.management.Project("cattle-map-Polygon Group/C1E-64-B3", "CattlePolygons", 4326)     Now we should probably do a little bit of cleanup: Every polygon appears twice in the featureclass (some remnant from converting from pdf to dxf), and the feature class includes the polygons from Hawaii and Alaska, which are now completely false. So we delete those polygons and then delete identical features:   arcpy.management.DeleteIdentical("CattlePolygons", ["Shape"])   That tool is actually not licensed for me and it would take too long to change my license right now, so I'll do it in the next step...   So now, everything that's left to do is to create points from the polygons:   # create the point feature class arcpy.management.CreateFeatureclass("", "CattlePoints", "POINT", spatial_reference=arcpy.SpatialReference(4326)) # for each polygon, insert its centroid into the point fc, if it isn't already there existing_coordinates = [] with arcpy.da.InsertCursor("CattlePoints", ["SHAPE@"]) as i_cursor: with arcpy.da.SearchCursor("CattlePolygons", ["SHAPE@"]) as s_cursor: for polygon, in s_cursor: point = polygon.centroid coordinates = [point.X, point.Y] if coordinates in existing_coordinates: continue existing_coordinates.append(coordinates) i_cursor.insertRow([point])       I'll attach the resulting shapefile to this post. This only includes the points for the continental USA. If you need Hawaii and Alaska, you'll have to ask for the data or do it yourself. Just start with georeferenceing the dxf to the respective state. ... View more