topic Re: how to convert a feature layer to point in ArcGIS API for Python Questions

how to convert a feature layer to point

MVAC — Tue, 01 Mar 2022 15:20:25 GMT

Hello,

I'm using jupyter notebooks (localhost, not in pro or AGOL), and I want to convert a feature layer to point. I've been trying to do it using arcpy.FeatureToPoint_management() but I'm having some issues:

parcels_id = 'the id'

# get layer

parcels = gis.content.get(parcels_id)

# get the required layer
parcels_layer = parcels.layers[0]

# variables
inFeatures = parcels_layer
outFeatureClass =  i dont know how to address this part

# feature to point
arcpy.FeatureToPoint_management(inFeatures, outFeatureClass, "INSIDE")

For the outFeatureClass, how can I do it to save it in a folder i have in my AGOL? is the inFeatures ok?

Thanks

Re: how to convert a feature layer to point

ABishop — Tue, 01 Mar 2022 15:46:09 GMT

Hello MVAC,

Not sure about using the feature to point to save directly to AGOL. I have a feature that I create from a feature class in our versioned geodatabase with definition query applied and saved to a project in ArcGIS Pro. I then upload it using the script below. Hope this helps.

import arcpy, os, sys from arcgis.gis import GIS from datetime import datetime import ctypes ctypes.windll.kernel32.SetConsoleTitleW('AgUseParcelsFeatureServiceUpdate') ### Start setting variables # Set the path to the project prjPath = r"I:\users\abishop\PY\FeatureServiceUpdates\AgUseParcels\AgUseParcels.aprx" # Update the following variables to match: # Feature service/SD name in arcgis.com, user/password of the owner account sd_fs_name = "Agricultural_Use_Parcels" portal = "https://mcpafl.maps.arcgis.com/" # ArcGIS Online account user = "abishop_MCPAFL" password = "********" # Set sharing options shrOrg = True shrEveryone = False shrGroups = "" ### End setting variables # Local paths to create temporary content relPath = os.path.dirname(prjPath) sddraft = os.path.join(relPath, "WebUpdate.sddraft") sd = os.path.join(relPath, "WebUpdate.sd") # Create a new SDDraft and stage to SD print("Creating SD file") arcpy.env.overwriteOutput = True prj = arcpy.mp.ArcGISProject(prjPath) mp = prj.listMaps()[0] arcpy.mp.CreateWebLayerSDDraft(mp, sddraft, sd_fs_name, "MY_HOSTED_SERVICES", "FEATURE_ACCESS","", True, True, allow_exporting=True) arcpy.StageService_server(sddraft, sd) print("Connecting to {}".format(portal)) gis = GIS(portal, user, password) # Find the SD, update it, publish /w overwrite and set sharing and metadata print("Search for original SD on portal…") sdItem = gis.content.search("{} AND owner:{}".format(sd_fs_name, user), item_type="Service Definition")[0] print("Found SD: {}, ID: {} n Uploading and overwriting…".format(sdItem.title, sdItem.id)) sdItem.update(data=sd) print("Overwriting existing feature service…") fs = sdItem.publish(overwrite=True) if shrOrg or shrEveryone or shrGroups: print("Setting sharing options…") fs.share(org=shrOrg, everyone=shrEveryone, groups=shrGroups) print("Finished updating: {} – ID: {}".format(fs.title, fs.id)) #provide input to stop script - uncomment input below if you need to see the script final output in the command window input('Press ENTER to exit')

Re: how to convert a feature layer to point

DanPatterson — Tue, 01 Mar 2022 15:53:00 GMT

Moved this thread to the arcgis-api-for-python since it isn't python specific.

If you want to use the FeatureToPoint tool work local and upload the result. Especially for shapefiles, you can only upload a zip of its parts

Re: how to convert a feature layer to point

MVAC — Tue, 01 Mar 2022 15:59:26 GMT

Would it be possible to do it on the notebook? I just checked and I don't need to save the point layer. I need this point layer to use arcgis.features.analysis.plan_routes() because I want to automate some stuff on workforce.

Re: how to convert a feature layer to point

DanPatterson — Tue, 01 Mar 2022 16:06:36 GMT

You will have to wade through the help topic for it

arcgis.features.analysis module — plan_routes

Re: how to convert a feature layer to point

jcarlson — Tue, 01 Mar 2022 16:46:35 GMT

Here's a purely ArcGIS Python API way of doing this. Well, plus a little pandas.

Query the layer
1. Return as a dataframe
2. Return centroid instead of geometry
The centroid is just a JSON object, not a true shape, so we need to convert it.
1. Export the centroid column to a dict
2. Read the dict back into its own dataframe
3. Convert the x/y columns to a true shape
Export the spatial dataframe to a new layer in your portal

from arcgis import GIS, GeoAccessor import pandas as pd # connect to portal gis = GIS('your url', 'username') # get parcel layer parcels = gis.content.get('parcels itemid').layers[0] # set our spatial reference sr = 4326 # or whatever your spatial reference is # query features to dataframe df = parcels.query( as_df = True, return_centroid = True, return_geometry = False, out_sr = sr ) # convert centroid to dict # ordinarily I don't like long single lines, but there's no point in assigning all the intermediate outputs to their own variables sdf = GeoAccessor.from_xy(df.join(pd.DataFrame(df.loc[:, 'centroid'].to_dict()).T), 'x', 'y', sr=sr) # drop the centroid columns sdf.drop(columns=['x', 'y', 'centroid', 'objectid'], inplace=True) # export to new feature class sdf.spatial.to_featurelayer( title = 'My Awesome Parcel Centroids', gis = gis, tags = ['parcels', 'centroids', 'outputs'], folder = 'My Outputs', sanitize_columns = True, service_name = 'parcels_centroids' )

You can just as easily export to other formats if you like.

Alternatively, you can use the GeoSeriesAccessor to take a spatial polygon column and derive the true_centroid of each, then assign the dataframe's shape attribute to that, but I tend to find the text → dict → x/y method performs a lot faster.

Re: how to convert a feature layer to point

MVAC — Wed, 02 Mar 2022 11:26:57 GMT

Thank you very much!

Is there a way to overwrite the feature layer every time I run the code?