Passing Python Variables into the Geospatial Modeling Environment?

Question

Hello,

I am trying to process a batch of rasters using the GME command isectpolyrst.  I am running ArcMap 10 and the newest standalone version of GME.  Does anyone know if it is possible to pass a list of .tif files using the command arcpy.ListRasters() from Python into the GME environment?  If so, how is that done?  The script that I started is below.  It runs without error but doesn't do anything.  I am a beginner when it comes to Python and all other programming so please take that into consideration.

Thank you!

Shawn O.
U.S. Fish and Wildlife Service (HAPET)
Great Falls, MT 59404

# import standard system modules
import arcpy, os, sys
from arcpy import env

# set workspace directory

Dir = "C:\Users\GISuser\Desktop\ONeil_2011\Scaup Folders\Habitat Attributes\2007"
env.workspace = Dir

# List all rasters in the workspace

rasters = arcpy.ListRasters()
grid = "C:\Users\GISuser\Desktop\ONeil_2011\Scaup Folders\all_grids.shp"

# Run each raster through the GME command "isectpolyrst"

for raster in rasters:
    os.system("C:\Program Files (x86)\SpatialEcology\GME\SEGME.exe isectpolyrst(in=\"grid\", raster = \"raster\", prefix = \"raster[:3])\", where \"Year=2007\";")

AlisonMontgomery · Answer

Hi Z,

Regarding that pop-up, I will also post in that forum...

As an update to my most recent post here, the script did run and it worked, sometimes it just takes a couple tries (annoyingly). Typically, with the scripts I wrote that integrate GME and Arc via Python, Arc and GME don't always play nice and I have a feeling the Arc software runs something that GME has a tough time handling. However, when I have had that pop-up appear, when I click ok, typically my scripts run, albeit on the 2nd or 3rd try.

I know the GME creator is switching everything to R so that he won't have to rely on the libraries. Already GME doesn't work with 10.2, or least I had an issue with it.

I created the scripts and use them as in a ArcToolbox, but they definitely run best when run directly via a Python IDE (I use PythonWin), and then are ok through ArcCatalog, and all the inputs have to be out of the map document if I run them in ArcMap.

The writing everything to a text file, then reading that into the run is definitely a lot cleaner and neater than trying to read in each variable. Pretty new to Python too, so thanks for the thumbs up!

Alison

ZacharyHart · Answer

Allison,

I know this thread is 2 months old now, but I'm scouring the forums in a desperate attempt to get to the bottom of an error that myself and (at least one other) user have experienced

Can you describe the exact behavior that occurs in GME when you execute the python script? Would it be possible for you to drop the '-c' switch to have a look?

Also, nice work. I had also thought about just taking all of my GME commands and stuffing them down into a text file to be run as a script for GME, but i'm not sure it will help me at this point given the error we're experiencing.

Thanks,
Z

ZacharyHart · Answer

Alison,

I know this thread is 2 months old now, but I'm scouring the forums in a desperate attempt to get to the bottom of an error that myself and (at least one other) user have experienced

Can you describe the exact behavior that occurs in GME when you execute the python script? Would it be possible for you to drop the '-c' switch to have a look?

Also, nice work. I had also thought about just taking all of my GME commands and stuffing them down into a text file to be run as a script for GME, but i'm not sure it will help me at this point given the error we're experiencing.

Thanks,
Z