I used to use SAGA GIS and automated simple tasks by dragging and dropping (or clicking something and copy and paste - don't recall sorry) to create a simple script that I could then run.  I've used this for example to mosaik a bunch of files and save as a new file.  I can then just update the list of files I want to mosaik and run the script.  I am working with a team who uses ArcGIS and was hoping that there was a simple way to automate simple tasks.  I know that you can use Python with ArcGIS, and I'm looking for an easy way to help them get started along the Python path.  I don't have access to ArcGIS itself yet, so can't just play and see.  I've checked out quite a few sites on this, but they focus on more advanced tasks.  I am somewhat guessing at the version being used - but it's a new system so I presume the most recent.
14130
Where can I start to learn ArcPy?
44247
I've used the integration of  with QGIS, and found it a very uesful tool. Having taken a look at the R toolbox  and the  it appears they both offer only point clustering to demonstrate the proof of concept.   Can other analysis methods available in R be implemented in ArcGIS Desktop through a python script (for example, implementing specific parts of the  package)?
45093
Has anyone used PostgreSQL with PostGIS successfully?   I am a semi-experienced ArcGIS user and want to build out a database and connect it to ArcCatalog.   I've scoured the internet for resources but nobody has really given a good walkthrough.
55497
I am very new to web-mapping, so this question may be a little misguided.  I want to have a very simple web interface using OpenLayers (or even Google Maps if that will be easier) which can display a number of point and polygon tables that are currently in PostGIS.  I want this to be as simple (and as low maintainence) as possible:   I don't want to have to create tiles (unless there is a very good reason), or have a massive, complex layer of middleware.   Performance doesn't need to be phenomenal, as this is for internal use and thus doesn't need to scale beyond a couple of concurrent users. Styling of the vector layers is of minimal importance. Open source tools are much preferred as I have almost no budget for this.   What is the recommended approach for this?
78606
I'm just starting to get into Python scripting for work.  I'm currently creating a script to automate a process.  Basically, it asks the user the client name, gets a projection if available, creates directory on C: drive for client, creates file geodatabase specific to client, creates required dataset, and creates feature classes specific to clients data. Eventually, it will also add the required fields to each feature class and probably some other things.  I started this not really knowing the proper etiquette of Python scripting for ArcMap. But what I've created so far will only run outside ArcMap I believe.   Is this acceptable?   Instead of getting user input through arcpy.getparamaterastext(), which I just found out about, I'm using raw_input().   Is this okay?   It does work, I'm just not sure if this is a proper way to do scripting.  Here's the code I have so far.  import sys import arcpy import os  #Records name of the client client = raw_input("Enter the name of the client: (letters and underscores only) \n")  #Records filepath of client to be created clientpath = "C:/" + client  #Inquires if projection file exists projection = raw_input("Is there a .prj or .shp available with correct projection? Y or N \n")  #Records the projection location if available if projection.upper() == "Y":     spatialr = raw_input("Drag the .prj or .shp here to record the filepath \n")     nspatialr = spatialr.replace('"', "") elif projection.upper() == "N":     alert = raw_input("You must add the spatial reference manually, hit enter to continue. \n") elif projection.upper() != "N" or "Y":     exit = raw_input("That is not a valid response. Try again. \n")     sys.exit()  #Checks if client folder exists; if not, creates one if not os.path.exists(clientpath):     os.makedirs(clientpath)  #Variable for file geodatabase location FGBpath = clientpath + "/" + client + ".gdb"  #Checks if client file geodatabase exists; if not, creates one if not arcpy.Exists(FGBpath):     arcpy.CreateFileGDB_management(clientpath, client)  #Variable for dataset location FDatasetpath = clientpath + "/" + client + ".gdb" + "/Network"  #Checks if dataset exists; if not, creates one if not arcpy.Exists(FDatasetpath):     if projection.upper() == "Y":         arcpy.CreateFeatureDataset_management(FGBpath, "Network", nspatialr)     elif projection.upper() == "N":         arcpy.CreateFeatureDataset_management(FGBpath, "Network")  #Variable for cable feature class location FCcablepath = clientpath + "/" + client + ".gdb" + "/Network" + "/cable"  #Checks if cable feature class exists; if not, creates one if not arcpy.Exists(FCcablepath):     if projection.upper() == "Y":         arcpy.CreateFeatureclass_management (FDatasetpath, "cable", "POLYLINE", "", "", "", nspatialr)     elif projection.upper() == "N":         arcpy.CreateFeatureclass_management (FDatasetpath, "cable", "POLYLINE")  #Variable for splice point feature class location FCsplicepath = clientpath + "/" + client + ".gdb" + "/Network" + "/splice_point"  #Checks if splice point feature class exists; if not, creates one if not arcpy.Exists(FCsplicepath):     if projection == 'Y' or projection == 'y':         arcpy.CreateFeatureclass_management (FDatasetpath, "splice_point", "POINT", "", "", "", nspatialr)     elif projection == 'N' or projection == 'n':         arcpy.CreateFeatureclass_management (FDatasetpath, "splice_point", "POINT")  exit = raw_input("\n\n File geodatabase, dataset, and the cable \n and splice point feature classes successfully created. \n\n Hit enter to exit.")   I still have some work to do, like adding the needed fields.
25432
I have been working with QGIS/PostGIS for a while now, without problems, but now I have to change QGIS for ArcGIS Desktop 10.1.  I am having trouble adding a layer from PostgreSQL and being able to edit that layer. I want to work directly with the database, and if I change something in the map it should be reflected on the database.  Can I do this without a Geodatabase?
106055
I am trying to make a script which selects every .png file in a folder beginning with the letters "LG". I then want the scipt create a shapefile, replacing the "LG" with "SH", and then i want the script to buffer that shapefile and rename the buffer with the first 2 letters being "SB".  I keep getting an error 99999 error message at line 37! ( gp.RasterToPolygon_conversion(INPUT_RASTER, Output_polygon_features, "SIMPLIFY", "VALUE") )  I have posted this on stackoverflow.com but probably here is a better place as scripting for arcmap is probably quite specialised!  Can anyone offer me any advice to why this script doesnt work? I am running ArcMap 9.3.1. and (attempting) to run the script from pyscripter!  # --------------------------------------------------------------------------- # RASTER2POLYGON.py # Created on: Wed Dec 01 2010 12:20:48 PM #   (generated by ArcGIS/ModelBuilder) # Usage: RASTER2POLYGON &lt;INPUT_RASTER&gt; &lt;Output_polygon_features&gt; &lt;Raster_Buffer_shp&gt; # ---------------------------------------------------------------------------  # Import system modules import sys, string, os, arcgisscripting  # Create the Geoprocessor object gp = arcgisscripting.create()  # Load required toolboxes... gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx") gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")  # Script arguments...  folder = "D:\\J04-0083\\IMAGEFILES" for root, dirs, filenames in os.walk(folder): # returms root, dirs, and files     for filename in filenames:         filename_split = os.path.splitext(filename) # filename and extensionname (extension in [1])         filename_zero = filename_split[0]     try:         first_2_letters = filename_zero[0] + filename_zero[1]     except:         first_2_letters = "XX"      if first_2_letters == "LG":         Output_polygon_features = "D:\\J04-0083\\ShapeFiles.gdb\\" + "SH_" + filename + ".shp"          # Process: Raster to Polygon...          INPUT_RASTER = os.path.join(root + "\\" + filename_zero + ".png")         gp.RasterToPolygon_conversion(INPUT_RASTER, Output_polygon_features, "SIMPLIFY", "VALUE")           Distance__value_or_field_ = "5 Meters"         Raster_Buffer_shp = "SB_" + filename + ".shp"          # Process: Buffer...         gp.Buffer_analysis(Output_polygon_features, Raster_Buffer_shp, Distance__value_or_field_, "FULL", "ROUND", "NONE", "")
110773
I have no skills in GIS, "map-making," or any other map program. Been a county planner for 25 yrs, where there was always a "drafting" staff.  I'm in another agency now and need to be able to produce basic and attractive site plans and area maps.  Is QGIS the right tool? I saw on a planners' website that it's free to download, along with GIMP and SketchUP.  I need some genuine advice on this, not something impossibly technical for my situation, nor a smug brush-off.
92026
I'm trying to write a script and as I want to automate it as much as possible.  Is it possible to set the workspace as a variable so that when it is run as a standalone script, or as a tool in ArcToolbox, that it will allow everything else to run?  Below, for example, I want the user to select a workspace and then for the script to pick up the shapefiles from this workspace and also create the new geodatabase there.   # Create file Geodatabase # Import ArcGIS modules import arcpy print arcpy.ProductInfo() # Check out the ArcGIS Spatial Analyst Extension arcpy.CheckOutExtension("spatial") # Need to be able to OverWrite Outputs arcpy.env.overwriteOutput = True # Set workspace arcpy.env.workspace = "C:/Prog_Data" # Set up variables out_folder_path = "C:/Prog_Data" out_name = "Prog.gdb" # Execute CreateFileGDB arcpy.CreateFileGDB_management(out_folder_path, out_name) # Set local variables before importing shapefiles to Geodatabase inFeatures =  ["Points.shp", "Extent.shp", "Rivers.shp", "Population.shp"] outLocation = "Prog.gdb" # Execute shapefile to Geodatabase arcpy.FeatureClassToGeodatabase_conversion(inFeatures, outLocation)
124210
Does ArcGIS Desktop or any of its native extensions allow me to create DEMs from stereo pairs (acquired by satellite, )?  If the answer is a clear no, brownie points will go to an answer that proposes a viable alternative using OpenSource.
77557
I'm working on a tool that converts raster layers to arrays for processing with NumPy, and ideally I would like to be able to work with rasters that come packaged in a .gdb without exporting them all (especially if this requires engaging ArcGIS or ArcPy).  Is this possible with the OpenFileGDB driver?   From what I can tell this driver seems to treat raster layers the same as vector layers, which gives you access to some data about the layer but doesn't give you the ReadAsArray functionality.
99847
I need to import Google Maps into ArcGIS for Desktop and apparently I need to activate a Google Maps extension. I'm working on ArcGIS 10.0 and there are no such extension. Do I have to install a new extension?
41476
I have written a script and am trying to run it in Python IDLE 2.7.6. I first try to import arcpy, but I get an error message saying:  ImportError: No module named arcpy  Why does it not work? I am using ArcGIS Desktop 10.2.
87188
One of my tasks for work is to divide parcels into groups. These groups will be used by agents to talk to property owners. The goal is to make the agent's job easy by grouping parcels that are near each other together, as well as divide the parcels into equal numbers so that the work is distributed evenly. The number of agents can fluctuate from a couple to 10+.  Currently I perform this task manually, but would like to automate the process if at all possible. I've explored various ArcGIS tools, but none seem to suit my need. I tried a script (in python) that makes use of near_analysis and selecting of polygons, but it's rather random and takes forever to accomplish a semi-correct result that then takes me longer to fix than if I just did everything manually from the start.  Is there a reliable method to automate this task?  Results example (hopefully without the division we see in yellow):
110167
I study computing on an advanced level and I need a program that will make it easier to create flowcharts and other diagrams efficiently. I don't need too many features; I just want something that works! I'm aware that there are many website that do exactly what I need but most of the times when I need to draw diagrams I don't have access to the internet
18936
I am trying to write a script to walk through all folders in my workspace and find and compact all geodatabases in my workspace.  I feel like I'm on the right track but I am stuck.  import acrpy, os from arcpy import env  arcpy.env.workspace "C:"  for (path, dirs, files) in os.walk(arcpy.env.workspace)     for file in files         if ".gdb" in file:         workspaces = arcpy.ListWorkspaces("*", "FileGDB")             for workspace in workspaces:             arcpy.Compact_management(workspace)   I am using ArcGIS 10.2.2 for Desktop.
38030
If I run this script:  import arcpy arcpy.AddMessage(arcpy.CheckExtension("Spatial")) from arcpy.sa import Con from arcpy import env fDir=r'd:\scratch\fdir' outFolder=r'd:\aerials\images' env.workspace = outFolder fDir=arcpy.Raster(fDir) Gter=Con(fDir&gt; 1,2) Gter.save("TEST")   from ArcGIS, the output is:    When I run it as stand-alone script, it is very different story:    How do I access Spatial Analyst extension from a stand alone script?