Latest Contributions by TomGeo

When creating buffers, you have several options to dissolve the buffers, but what I am missing in Buffer, as well as in Multiple Ring Buffer, is the possibility to specify at least one attribute of the input features to be transferred onto the resulting buffers. Of course you can argue that I can do a spatial join, or intersect or something else to make the connection, but on one hand that means an extra geoprocessing tool run, and on the other hand, that only works as long as the object used to create the buffer is intersecting, or touching the buffer, or you search within a specific distance. Especially when creating multi ring buffers the intersecting and/or touching falls possibly short in case of all the ring buffers. The buffersize I use to create the multiple ring buffers depend on a feature value in the attribute table, which goes into a formula and the resulting distance for the buffers is calculated on the fly. Thereby, fishing for buffers in the vicinity is not much fun. Hence, I propose that the buffer tools functionality is extended by the possibility to transfer one, or multiple feature attributes from the initiating feature over to the resulting buffer(s).   ... View more

I need to project geometries in an spatial dataframe and try to do it with this code   import pandas as pd from arcgis import geometry from arcgis.features import GeoAccessor, GeoSeriesAccessor in_csv = r'C:\Path\to\points.csv' point_df = pd.read_csv(in_csv, sep=';') point_df = pd.DataFrame.spatial.from_xy(point_df, x_column='Longitude', y_column='Latitude', sr=4326) point_df['SHAPE'] = geometry.project(geometries=point_df['SHAPE'], in_sr=4326, out_sr=3006)   A check on the geometry column in the SeDF returns I am logged into our Enterprise portal when running the code, and I am getting: Exception: Invalid Token (Error Code: 498)   I am running ArcGIS Pro 3.0.2, our portal is in version 10.9 and arcgis is in version 2.0.1 PS: I am getting the same error when trying to find the transformation needed to do the projection   arcgis.geometry.find_transformation(in_sr=4326, our_sr=3006)     ... View more

I am just getting started with SeDF and I was wondering on how I can actually visualize objects with multiple geometries. Following the introduction, a buffer geometry is added to the dataframe. The map widget of my notebook is however, only showing the points. How can I visualize the buffer geometries? Is there a way to specify what geometry column should be used for the visualization?     # generates map widget visulizing the points, though info widgets places two objects there, when I click on a point point_df.spatial.plot() # AttributeError: 'Series' object has no attribute 'spatial' point_df.Buffer_2.spatial.plot() # Since plot generates a MapView, the object is not subscriptable point_df.spatial.plot()['Buffer_2'] # or has no attribute 'Buffer_2' point_df.spatial.plot().Buffer_2     ... View more

I was wondering about how to store several hundred LAS files in an ESRI environment. We have ArcGIS Enterprise running and I would like to make the point cloud data of our entire country available to everybody within my organization. What is the best way possible to do that, so that the data are easily and fast accessible for ArcGIS Pro users, and that the data at the same time can be used for analysis on a GeoAnalytics Server? ... View more

I am working on a little task, where I create some points, based on coordinates. Every point is then used to create a couple of polylines, where each point is supposed to be the start point for a set of polylines. The endpoints for the polylines are determined by angle and distance. So, the outcome should be a point, in the center of some polylines, looking like sunrays. I create the points and use the created points thereafter as point of origin for my polylines. Only problem, the created point is not where the polylines 'meet' and the polylines do not meet at the same point either. I wonder how that can be, if the point of origin does not change!   import arcpy import arcgis from arcgis.gis import GIS from arcgis.geometry import Point, Polyline, Geometry bearing_list = [(345,165),(315,135),(285,105),(255,75),(225,45),(195,15)] distance_2 = 2000 in_sr = 4326 out_sr = 3006 # create featureclass and populate with points from coordinates points_fc = arcpy.CreateFeatureclass_management(out_name='Sites', geometry_type='POINT', spatial_reference=out_sr) arcpy.AddField_management(points_fc, 'SiteName', 'TEXT') point_insert_cursor = arcpy.da.InsertCursor(points_fc, ('SiteName', 'SHAPE@XY')) pnt_list = [] for item in point_dict: y, x = point_dict[item] pnt = Point({ 'x': x, 'y': y, 'spatialReference':{ 'wkid': in_sr, 'latestWkid': in_sr } }) out_pnt = pnt.project_as(out_sr) pnt_list.append([item, out_pnt]) row_point = (item, out_pnt.centroid) point_insert_cursor.insertRow(row_point) del point_insert_cursor # iterate over points, create temporary line featureclass, temp_lines = arcpy.CreateFeatureclass_management(out_name='temp_lines', geometry_type='POLYLINE', spatial_reference=out_sr) arcpy.AddField_management(temp_lines, 'SiteName', 'TEXT') line_insert_cursor = arcpy.InsertCursor(temp_lines) line_length = distance_2 + 5 for pnt in pnt_list: for bearing_set in bearing_list: pnt_array = arcpy.Array() x = pnt[1].first_point.x y = pnt[1].first_point.y pg = arcpy.PointGeometry(arcpy.Point(x, y), arcpy.SpatialReference(out_sr)) start = pg.pointFromAngleAndDistance(bearing_set[0], line_length) start_pnt = arcpy.Point(start.firstPoint.X, start.firstPoint.Y) end = pg.pointFromAngleAndDistance(bearing_set[1], line_length) end_pnt = arcpy.Point(end.firstPoint.X, end.firstPoint.Y) pnt_array.append(start_pnt) pnt_array.append(end_pnt) line_feature = line_insert_cursor.newRow() line_feature.shape = arcpy.Polyline(pnt_array) line_feature.setValue("SiteName", pnt[0]) line_insert_cursor.insertRow(line_feature) del line_insert_cursor     The result for one of the points: the point is 10 centimeters away from where the lines 'meet'.   And it looks like this at the 'center' of the polylines   ... View more

I was creating a web map app, adding a web map that I styled just to realize that the app is interpreting the colors and apparently also the blending settings in a different way than the map is set to. To illustrate... this is my map with some bathymetry data styled and in blending mode 'Multiply', showing OpenStreetMap information in the waters. This is the app, showing its interpretation of the map above...  You can clearly tell, the blending option is disregarded, and I assume thereby the colors for the bathymetry are changed as well! I created the app with the Web AppBuilder on a 10.9.1 portal, where the map was created in the Map Viewer (not the classic). Did I miss a setting enabling blending, or is this a bug? On another note... one of the layers I show in the map has labels turned on and shows them nicely in the map. When I visualize the layer in the app, the labels are not shown by default, but instead the user has to manually enable them through some settings menu of the layer. Again, that is altering the Web Map when used in the App   ... View more

I recently experienced that while trying to create a layer from one of our portal services my AGP crashed. Turned out that the service I requested to create the layer from was stopped on the server. To avoid that, I would like to check on the PortalItem (Item) status, but I do not seem to have the possibility to do so (Item Class Members). Fingers crossed I am only looking at the wrong place. Can I check the PortalItem status? ... View more