return extrap_site #------------------------------------------------------------ #-------------------------------------------------------------------------- #-------------------------------------------------------------------------- #MAIN if __name__ == "__main__": try: #---- Create the Geoprocessing Object ---------------------- gp = arcgisscripting.create(9.3) #----- Create the Custom Error Object ---------------------- e = SvmpToolsError(gp) # Set some basic defaults for error handling e.debug = True e.full_tb = True #-------- unit conversion flag -------------------- unit_convert = "sf2m" def msg(msg): gp.AddMessage(msg) # ----------- PARAMETERS ---------------------------------- # Temporary stand-ins for parameters #dbPath = "c:/projects/dnr_svmp/output/eelgrass_monitoring/" #blPath = "c:/projects/dnr_svmp/output/eelgrass_monitoring/base_layers/"
'''Create dictionary of extrapolation type by site id using lookup from analysis/geo stratum to determine analysis and extrapolation type ''' site_extrap = {} for (siteid,data) in sites_strata.items(): site_extrap[siteid] = lookup_dict[tuple(data)] # invert this dictionary to get sites by extrapolation/analysis stratum extrap_site = invert_dict(site_extrap) return extrap_site #---- Create the Geoprocessing Object ---------------------- gp = arcgisscripting.create() #----- Create the Custom Error Object ---------------------- e = SvmpToolsError(gp) # Set some basic defaults for error handling e.debug = True e.full_tb = True #-------- unit conversion flag -------------------- unit_convert = "sf2m" #------------------ DATA QUERIES AND GROUPING -------------------------------- #------------------------------------------------------------------------------- #----------------- MATCHING SITES sampled in YEARS 1 & 2 ----------------------- #------------------------------------------------------------------------------- # Get list of sites sampled in both years and their strata sites_strata_2yr = sites_strata(gp,svmp_all_sites_fc,sample_group,year1,year2) sites2yr = sorted(sites_strata_2yr.keys()) # For use in query strings
def msg(msg): gp.AddMessage(msg) shapefile = "C:/projects/dnr_svmp/output/eelgrass_monitoring/site_folders/nps1446/video_transect_data/2008_nps1446_transect_data.shp" selStatement = '"tran_num" = 1' try: # Create the geoprocessing object gp = arcgisscripting.create() # Overwrite existing output data gp.OverWriteOutput = 1 # Create the custom error class # and associate it with the gp e = SvmpToolsError(gp) # Set some basic defaults for error handling e.debug = True e.full_tb = True #e.exit_clean = False #e.pdb = True msg(selStatement) rows = gp.UpdateCursor(shapefile,selStatement) row = rows.Next() while row: id = row.getValue("Id") date = row.getValue(utils.shpDateCol) msg(id) msg(date) row = rows.Next()
import os import sys import traceback import arcgisscripting from svmp_exceptions import SvmpToolsError DEBUG = True PDB_DEBUG = False # Create the geoprocessing object gp = arcgisscripting.create() # Create the custom error class # and associate it with the gp error = SvmpToolsError(gp) # Set some basic defaults for error handling error.debug = True error.pdb = True def main(): try: raise RuntimeError except Exception, E: error.call('An intentional RuntimeError occurred in the python script..') if __name__ == "__main__": try: main() except Exception, E: error.call(E)
#-------------------------------------------------------------------------- #-------------------------------------------------------------------------- #MAIN if __name__ == "__main__": try: # Create the geoprocessing object gp = arcgisscripting.create() # Overwrite existing output data (for testing only?) gp.OverWriteOutput = 1 # Create the custom error class # and associate it with the gp e = SvmpToolsError(gp) # Set some basic defaults for error handling e.debug = True e.full_tb = True #-- Get parameters ----- # Input Data Parent Directory # Transect ASCII files are located in subdirectories below here # Default: Environment - current workspace # outParentDir: Output Data Parent Directory. Default: Environment - scratch workspace # siteFile: Full Path of text file containing list of sites to process inParentDir,siteFile,outParentDir = get(0),get(2),get(3) # Input Data Coordinate System. Default: GCS_WGS_1984 # Ouput Data Coordinate System. Default: Default: NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet # Survey Year for data to be processed inCoordSys,outCoordSys, surveyYear = get(1),get(4),get(5)
#-------------------------------------------------------------------------- #-------------------------------------------------------------------------- #MAIN if __name__ == "__main__": try: # Create the geoprocessing object gp = arcgisscripting.create() # Overwrite existing output data gp.OverWriteOutput = 1 # Create the custom error class # and associate it with the gp e = SvmpToolsError(gp) # Set some basic defaults for error handling e.debug = True e.full_tb = True #e.exit_clean = False #e.pdb = True #Get parameters # Input Transect Point Shapefile Parent Directory # Transect ASCII files are located in subdirectories below here ptParentDir = gp.GetParameterAsText(0) # Input Sample Polygon Shapefile Parent Directory pyParentDir = gp.GetParameterAsText(1) # Control File Parent Directory ctlParentDir = gp.GetParameterAsText(2) # Full Path of text file containing list of sites to process