Latest Contributions by DavidPike

‎03-25-2023

This is exactly the problem, 0 is still a value. You need to change it to NODATA then re-run the Euclidean distance with that new layer. I'm sure it will work as intended then.

‎03-25-2023

Could be an environment setting of processing extent or output extent, best to check the model builder environment settings and the settings for each tool.

‎03-25-2023

what's the cell values of 1 on flow_rec1? can you change that symbology to a colour rather than blank?

‎03-25-2023

Can you share a screenshot of flow_recl1 ? When you click on a few pixels of the euclidean distance result and flow_rec1 what values do you get? To me it just looks like you've set the other raster cells in flow_rec1 to no symbology rather than actually reclassifying the data https://pro.arcgis.com/en/pro-app/latest/tool-reference/spatial-analyst/reclassify.htm. What are the cell values of 1 in the screenshot? This would explain the Euclidean distance result, as the tool doesn't care about how you symbolise the input data.

‎03-24-2023

ah glad it may work, I've no experience of doing it however. The documentation also seems a bit arcane. I'm assuming it's returning a set of points where they intersect? the doc here suggest a 'dimension' parameter https://developers.arcgis.com/python/guide/part3-spatial-operations-on-geometries/ so perhaps: yourGeometry.intersect(yourOtherGeom, dimension=2) a4. intersect() V.S. overlaps() The intersect() constructs a geometry that is the geometric intersection of the two input geometries. Different dimension values can be used to create different shape types. The intersection of two geometries of the same shape type is a geometry containing only the regions of overlap between the original geometries, and its arguments include: second_geometry: required arcgis.geometry.Geometry. A second geometry dimension: required Integer. The topological dimension (shape type) of the resulting geometry, which can be: 1: A zero-dimensional geometry (point or multipoint). 2: A one-dimensional geometry (polyline). 4: A two-dimensional geometry (polygon).

‎03-24-2023

What's your search radius? Do you have a population field? The Line Density tool probably doesn't work the way you think it does by counting the lines intersecting each search radius - it takes the length into account also https://pro.arcgis.com/en/pro-app/latest/tool-reference/spatial-analyst/how-line-density-works.htm What do you (exactly) need to achieve? You may be better using some aggregation methods using vector only.

‎03-24-2023

Hi all, I have published a hosted feature layer to ArcGIS Online AGOL using the data interoperability extension (in other words - FME). Note I am no expert on FME or Data Interoperability. My original data is in a United Kingdom format called 'British National Grid', 'BNG' for short. EPSG:27700 Airy (1880) Spheroid and OSGB1936 (also called OSGB36) Datum. The process works as intended and creates the hosted feature layer, however the service is EPSG 27700 (which is what I want) but the default transformation is probably sub-optimal. As the service could likely be used with a Web Mercator or WGS84 basemap, this may affect both performance and accuracy. This is the transformation which I believe is the default for the service: "datumTransformations" : [ { "transformDirection" : "forward", "transformJson" : "[{\"wkid\":1149,\"latestWkid\":1149,\"transformForward\":false,\"name\":\"ETRS_1989_To_WGS_1984\"},{\"wkid\":5338,\"latestWkid\":5338,\"transformForward\":false,\"name\":\"OSGB_1936_To_ETRS_1989_1\"}]", "fromSR" : { "wkid" : 4326, "latestWkid" : 4326 }, "toSR" : { "wkid" : 4277, "latestWkid" : 4277 } }, { "transformDirection" : "reverse", "transformJson" : "[{\"wkid\":5338,\"latestWkid\":5338,\"transformForward\":true,\"name\":\"OSGB_1936_To_ETRS_1989_1\"},{\"wkid\":1149,\"latestWkid\":1149,\"transformForward\":true,\"name\":\"ETRS_1989_To_WGS_1984\"}]", "fromSR" : { "wkid" : 4326, "latestWkid" : 4326 }, "toSR" : { "wkid" : 4277, "latestWkid" : 4277 } } ], If my interpretation is correct, for reprojection on-the-fly it is transforming OSGB36 datum into ETRS, then ETRS to WGS84. I think there must be better ways to set this. Questions: 1. For any UK folks or geodesy geeks - I think that transformation OSGB_1936_To_WGS_1984_7 "wkid": 5339 would be preferred. Is that a fair assumption? 2. Can this transformation be configured in Data Interoperability/FME (fully aware that this is more of an FME question, so will post in their forum also. Advice is appreciated however). 3. Can the default transformation be updated using the REST API?

‎03-24-2023

found this for posterity https://support.esri.com/en/technical-article/000018479

‎03-24-2023

I believe they are stored in your local user profile i.e. %appdata% in a 'user.config' file.

‎03-24-2023

It's in German as its more efficient. I think you just want to clear the workspace environment like the sample, try it exactly as below: arcpy.ClearEnvironment("workspace") # prints None print(arcpy.env.workspace)

‎03-24-2023

try intersect()

‎03-23-2023

From the docs I'd say that LengthGeodetic() input is the feature, not a length field from the feature (which wouldn't make any real sense). https://developers.arcgis.com/arcade/function-reference/geometry_functions/ var test = Round(LengthGeodetic($datapoint, 'feet'),2) LengthGeodetic LengthGeodetic(inputGeometry, unit?) -> Number Since version 1.7 Function bundle: Geometry Returns the geodetic length of the input geometry or Feature in the given units. This is more reliable measurement of length than Length() because it takes into account the Earth's curvature. Support is limited to geometries with a Web Mercator (wkid 3857) or a WGS 84 (wkid 4326) spatial reference. Feature geometries in the visualization and labeling profiles are generalized according to the view's scale resolution to improve drawing performance. Therefore, using a feature's geometry (i.e. $feature) as input to any geometry function in these contexts will return different results at each scale level. Other profiles, such as popup, provide the full resolution geometry. Parameters inputGeometry: Geometry | Feature | Array<Point> - The geometry for which to calculate the geodetic length. unit (Optional😞 Text | Number - Measurement unit of the return value. For the visualization, labeling, and popup profiles, the default unit is the map's spatial reference. In other profiles, like field calculation, the default is based on the spatial reference of the data. Possible values: feet | kilometers | miles | nautical-miles | meters | yards Return value: Number Example Returns the geodetic length of the feature in kilometers LengthGeodetic($feature, 'kilometers')

‎03-23-2023

I won't pretend to know much about the heuristic algorithm but I do think the hierarchy is causing the strange results in even a relatively simple solve. I'd recommend analysing without a hierarchy to test this - also having a hierarchy seems counter-intuitive if your only question is which is the shortest/fastest route. from: https://desktop.arcgis.com/en/arcmap/latest/extensions/network-analyst/network-analysis-with-hierarchy.htm Using hierarchy in network analysis It makes more sense to use a time-based attribute as impedance when performing a hierarchical network analysis. If the impedance is not based on time, using a hierarchy may not produce realistic results. For analyses that don't minimize time, consider solving without using the hierarchy

‎03-23-2023

Can you share the output of your prints? I think it may be something to do with the raw strings then how you're concatenating them, it misreads the string and puts it into the workspace, I much prefer to use os.path.join perhaps try https://pro.arcgis.com/en/pro-app/latest/arcpy/functions/clearenvironment.htm and see if it fails with a helpful message when trying to write to that location.

‎03-22-2023

Is the script from a previous answer? Unless I'm misinterpreting, you just need the syntax for an if statement? import arcpy docPath = r"*PATH*" docName = r"Test.aprx" layoutName = 'name of your layout' p = arcpy.mp.ArcGISProject(docPath+"\\"+docName) # ge a list of all layouts in your map dox lLayout = p.listLayouts() for l in lLayout: print(l.name) #if statement to print ONLY a specific layout if l.name == layoutName: l.exportToJPEG(docPath+"\\"+l.name) print("Process Completed")

Online Status	Offline
Date Last Visited	‎07-06-2026 06:42 AM

My Ideas

Latest Contributions by DavidPike

Re: Eucledian Distance not working

Re: tool created for compactness is having issues

Re: Eucledian Distance not working

Re: Eucledian Distance not working

Re: Clipping Polyline from Polygon - Geometry

Re: Line Density by Line Count

ArcGIS Online Hosted feature service default transformation parameters

Re: Where is the file stored which contains ArcGIS Pro General Options choices?

Re: Where is the file stored which contains ArcGIS Pro General Options choices?

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

Re: Clipping Polyline from Polygon - Geometry

Re: Difficulty getting indicator to utilize length geodetic

Re: Network analysis find shortest route question

Re: ERROR 000875: I cannot get my raster into the dedicated output workspace while using ZIP- loop

Re: How do I call out a specific layout window to export?

Re: Calculate precise population using 2020 Census...

Re: Define Projection crack-down

Re: Generate polygons that represent historic air ...

Re: Enterprise geodatabase connection slow - machi...

Re: Discrepancy between measure tool and expected ...

GIS Jobs