|
POST
|
Okay, I figured it out! The image server is using EPSG:3857 as default coordinate system! I created an imagery layer actively deleting the horizontal coordinate system in the property settings. Still, without any coordinate system provided, the server is using the default coordinate system and not the one provided by the input data from the data store! So, I specified EPSG:3035 in the properties and provided the transformation ETRS_1989_To WGS84... it worked. But please ESRI, can you change the default behavior to recognize the coordinate system provide by the input data, and not overriding them by the server defaults?
... View more
12-22-2022
12:40 AM
|
0
|
0
|
1266
|
|
POST
|
When creating an imagery layer in our portal I select 'image collection' as layer configuration and continue to define the imagery. All my images are tif files, results from satellite imagery classifications, and have spatial reference 3035. In the following I will illustrate the property settings Changing the default horizontal coordinate system (WGS 1984 Web Mercator (auxiliary sphere)) to ETRS 1989 LAEA has as consequence that the imagery layer is created, but will not visualize in one of the portal viewers, or in ArcGIS Pro. It will also not show any error when loading the layer, indicating something is wrong with it. If I keep the default horizontal coordinate system, the layer is created and shows right from the start with a thumbnail, and shows in ArcGIS Pro but will of course be in WGS 1984 Web Mercator (auxiliary sphere) and has to be reprojected on the fly. Can one of you explain to me why I cannot keep the imagery layer at the original coordinate system? The property settings for spatial reference do only point at the output layer, and since I do not want to have any other coordinate system applied on my data, I set the output coordinate system to the input coordinate system. Since, I do not change anything , I would a) not know why I should select a datum transformation, and b) not know what transformation to choose.
... View more
12-21-2022
09:36 AM
|
0
|
1
|
1309
|
|
POST
|
I was trying to use Streamlit with my arcgis pro python clone, and while I can install version 1.10.0 of Streamlit, I cannot import it. import streamlit as st returns the following File "C:\...\AppData\Local\ESRI\conda\envs\pro-py3-open30\lib\site-packages\pyarrow\__init__.py", line 69, in <module> from pyarrow.lib import (BuildInfo, RuntimeInfo, MonthDayNano, ImportError: cannot import name 'RuntimeInfo' from 'pyarrow.lib' (C:\...\AppData\Local\ESRI\conda\envs\pro-py3-open30\lib\site-packages\pyarrow\lib.cp39-win_amd64.pyd). Does anybody made better experiences using Streamlit? I also had to realize that Python 3.9.x is in 'security fix only' state and thereby the last version coming with a binary installer is 3.9.13... ArcGIS Pro 3.0.2 is using 3.9.11 and I was wondering with Python in stable release of 3.11.1, if ArcGIS Pro 3.1 will come with a more recent version of Python?
... View more
12-09-2022
05:26 AM
|
0
|
1
|
2648
|
|
POST
|
Thanks for all the valuable thoughts and inputs. In the end, since I had all the logic anyway in Python, I ended up using an update cursor. Time wise, it is almost the same if I run those two CalculateField runs, slightly slower... I tried the whole thing of course straight in ArcGIS Pro and there my code block works. However, I think @RogerDunnGIS has a fair point talking about the format of fields. The first CalculateField run is creating the buffer field, as 'TEXT'! But my second CalculateField run is using integers in the if-statement... I haven't had time to try it, but 1 should not be '1', and thereby the error source is once more found about 40 cm in front of the screen. 😅
... View more
12-02-2022
01:02 AM
|
0
|
0
|
1696
|
|
POST
|
I have one calculation where it works great getting a variable in: arcpy.management.CalculateField(sector_angle, 'Buffer', "getBuffer(!Area!)", "PYTHON3",
f"""def getBuffer(area):
if area < {median_area}:
return 1
if area > {median_area}:
return 2""",
"TEXT", "NO_ENFORCE_DOMAINS") And I have another calculation on the same feature class where it does not work. arcpy.management.CalculateField(sector_angle, 'Buffer_Radius', "setDistance(!Buffer!)", "PYTHON3",
f"""def setDistance(buffer):
if buffer == 1:
return {pnt[2]}
elif buffer == 2:
return {pnt[3]}""", "FLOAT") Some example on how the variables look: pnt[2] = 2043.8900146484375 pnt[3] = 4087.780029296875 The type of both variables is float, but the resulting 'Buffer_Radius' field contains only NaN. I tried all sorts of stuff on how to get the variables into the code block, but while it was so easy in the first function, I don't have luck with the second one. How do I get my variables into the 'Buffer_radius' field?
... View more
12-01-2022
06:02 AM
|
0
|
5
|
1790
|
|
POST
|
I need to analyze the area covered by some polygons. The analysis are based on mosaic and LAS datasets. Since the mosaic/LAS datasets contain several thousand raster and las files it seems to take ages to do the analysis. For example, to create my workflow I used a minimal subset of the LAS files to create a LAS dataset and the *StackProfile* analysis went within seconds. Now that I have more las files building the LAS dataset it takes 2 minutes to do the analysis. Hence, I was thinking to use the polygon to identify the files of the mosaic/LAS dataset intersecting with the polygon, extract their name (or better source path) and create smaller mosaic/LAS datasets on-the-fly. Can somebody help me out on how to query a mosaic/LAS dataset for the footprints and get the footprint name or source path returned?
... View more
11-29-2022
12:17 AM
|
0
|
0
|
971
|
|
POST
|
Yes, together with a colleague I tried today the set_geometry method and we experienced the same, also the issues with plotting the geometry into the map widget we can confirm. Thanks for looking into it and the test with geometry_type and testing against the geometries area.
... View more
11-28-2022
12:06 PM
|
0
|
0
|
2255
|
|
IDEA
|
Agreed @DanPatterson, it offers the possibility to link back, and join fields needed. Of course it does not work with the multi ring buffer tool, and it requires the join as an extra step in the users procedure, while adding fields and values while creating the buffers, seems much more straight forward, since the data are at that point at hand.
... View more
11-28-2022
12:25 AM
|
0
|
0
|
1058
|
|
IDEA
|
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
11-27-2022
01:38 PM
|
0
|
2
|
1099
|
|
POST
|
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
11-27-2022
02:15 AM
|
0
|
0
|
1168
|
|
POST
|
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
11-26-2022
06:27 AM
|
0
|
4
|
2325
|
|
POST
|
Certainly bookmarked your blog, and I am pretty sure I will steal some of your examples. 😉
... View more
11-20-2022
02:37 AM
|
0
|
0
|
2123
|
|
POST
|
Hi @DanPatterson, thanks for the advice about the coordinate system when creating the geometry. It didn’t change anything on the outcome though. 🙂 Neither did density, and the only thing that solved the issue with the point of origin was to create all lines from the point of origin, leaving me with 12 instead of 6 lines. I assume the culprit is actually the precision of the points, calculated by pointFromAngleAndDistance. Meaning calculating points on a circle with a distance of 2000 meters, and angles of 345 and 165 degrees will give you points, not marking a precise opposite of each other on a circle. Once you connect these two points by a line, the line has good chances not to cross your point of origin. For my purposes it is not important that the point where all polylines meet is spot on on top of the point used for the calculations (even though that’s what I would expect). But it is important that all polylines have a common cross point. Later in the process these polylines are used to segment a buffer around my calculation point into 12 equal circle segments. As long as the polylines do not have a common crossing point, I will create more than 12 segments, where an unknown number of them will simply be slivers, and in the end the 12 expected circle segment differ in their area size.
... View more
11-19-2022
11:27 PM
|
0
|
0
|
2133
|
|
POST
|
Hard to say what you are doing wrong. Since you are doing it in FME, this also might be the wrong community for the question. However, off the top of my head there would be three things. I would use the transformer to reproject the data. Setting the coordinate system does not change the features projection. When creating the writer, I would set the schema to ‘Automatic’. Default is to inherit from the reader In the Navigator sidebar you can set the reference system of the writer to EPSG:4326. The writer should reject thereafter features projected in other coordinate systems
... View more
11-19-2022
10:53 PM
|
0
|
0
|
1048
|
|
POST
|
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
11-19-2022
11:11 AM
|
0
|
0
|
911
|
| Title | Kudos | Posted |
|---|---|---|
| 1 | 10-24-2022 01:23 AM | |
| 1 | 09-15-2021 01:21 AM | |
| 1 | 08-30-2022 04:31 AM | |
| 2 | 04-24-2024 04:23 AM | |
| 3 | 12-12-2023 07:49 AM |
| Online Status |
Offline
|
| Date Last Visited |
11-08-2025
10:36 PM
|