nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

use python to create definition query

RafiqBasaria

Hi. Ive been trying to figure out how to accomplish this, but Ive run into some trouble. Im not sure if I am coding this correctly. My goal, overall, is to select features from a layer that intersect with features from another layer. The features in this second layer need to be defined by a definition query that I do not want to manually change. The selected features will then be exported to a new shapefile. I have been able to get the script to select and export, but I cannot get it to create the definition query. Any assistance will be greatly appreciated. Here is the code that I have, so far. Thanks!!

import arcpy

mxd = arcpy.mapping.MapDocument("C:/Rafiq_GIS/PYTHONTESTING/2011/2011_CountyComm/Data/DataPrep.mxd")

for lyr in arcpy.mapping.ListLayers(mxd):
    if lyr.name == "CountyCommission_TriCounty":
        layer.definitionQuery = "[DISTRICT] = '1' AND [COUNTY] = 'ORANGE'"

arcpy.SelectLayerByLocation_management("BusRoutes_Ridership","INTERSECT","CountyCommission_TriCounty","","NEW_SELECTION")

arcpy.CopyFeatures_management("BusRoutes_Ridership","C:/Rafiq_GIS/PYTHONTESTING/2011/2011_CountyComm/Data/Orange/Or_Com1.shp")

Encuentra más publicaciones etiquetadas con

Accepted answers

All comments

rzufelt


import arcpy
from arcpy import env
import os


mxd = arcpy.mapping.MapDocument("CURRENT")
df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]
lyr = arcpy.mapping.ListLayers(mxd, "RevisedGroundZones", df)[0]
lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''

So, similar to your code, could do it this way as well if you wanted to iterate through the list instead of just grabbing the first item (normally first of one):

for lyr in arcpy.mapping.ListLayers(mxd, "*"):
     if lyr.name == "RevisedGroundZones":
          lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''

The "List"Layers and "List"DataFrames create a list of matching objects, so, all Dataframes with name "Layers" would be appended to the df list, and all layers, within the specified dataframe, that match "RevisedGroundZones" would be appended to the lyr list. However, putting the [0] actually makes it grab the first item in the list. I am assuming here that you only have one dataframe named "Layers" and one layer within that named "RevisedGroundZones". That way, there would only be one item in each list, so putting the [0] on the end grabs it without having to iterate through a list of one.

I have successfully ran this and it applies the definition query (of course, this is for text field, if your RevisedGroundZones attribute is a number field, the quotes need to dissapear from the definition query.

So, if you only have one dataframe name "Layers" (or "layers" as list function here is not case sensitive), and within that only one layer named RevisedGroundZones, it should be applying the definition query. If you have more than one with this name, it would only apply to the first on in the list. Would then need to iterate through if you wanted it applied to all with that name.

Of course, you can always throw a mxd.save() in there after the definition query, then open the mxd with arc map and see what the actual def query being applied is (and if there are any drawing errors related).

Once you get the def query working, then let it move on to the clipping.

R_

JonathanCusick

Thanks for the help. I tried running the script with the updates you made, and while the script runs all the way (doesn't receive any error messages), it doesn't actually apply the definition query. When I look at the newly clipped layer in ArcMap, it is clipped to the entire Ground Zones layer as opposed to just ground zone 1 or 2 for example.

Also, what is the function of the [0] at the end of the 6th and 7th lines?

rzufelt

Ah I didn't realize I had to do so. That's all of the code referring to the query. The rest of it begins to run the clip function along with others down the line. So I have to define the mxd that this is happening in in order to call upon the attributes of the Ground Zones feature class? Is there a way to run this just in Catalog?

I tried defining the mxd as shown in the post at the top of the thread, but am still having problems. Whenever I try to set mxd to the full pathname of the map, I get an error saying Invalid MXD Filename. When just using the keyword "CURRENT" to define the mxd, it at least ran, but then I got that same 'str' object has no attribute 'name'.
import arcpy
from arcpy import env
import os

mxd = arcpy.mapping.MapDocument("CURRENT")
df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]
lyr = arcpy.mapping.ListLayers(mxd, "RevisedGroundZones", df)[0]
lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''

Think you just need to define the layer. If it is just one layer, no need to iterate.

R_

JonathanCusick

Ah I didn't realize I had to do so. That's all of the code referring to the query. The rest of it begins to run the clip function along with others down the line. So I have to define the mxd that this is happening in in order to call upon the attributes of the Ground Zones feature class? Is there a way to run this just in Catalog?

I tried defining the mxd as shown in the post at the top of the thread, but am still having problems. Whenever I try to set mxd to the full pathname of the map, I get an error saying Invalid MXD Filename. When just using the keyword "CURRENT" to define the mxd, it at least ran, but then I got that same 'str' object has no attribute 'name'.

import arcpy
from arcpy import env
import os

mxd = arcpy.mapping.MapDocument("CURRENT")
inzones = (.../RevisedGroundZones)
for lyr in inzones:
     if lyr.name == "RevisedGroundZones":
          lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''

MathewCoyle

I've been trying to use a definition query just like rbasaria, and have used the syntax that clm42 suggested, but keep running into trouble. For now, I just want to be able to query out various zones from a larger map, and then run further geoprocessing functions on them.
import arcpy
from arcpy import env
import os

inzones = (.../RevisedGroundZones)
for lyr in inzones:
     if lyr.name == "RevisedGroundZones":
          lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''
When running it from the python window in ArcMap, I get an error message that says AttributeError: 'str' object has no attribute 'name' and when I run it as a script tool in ArcCatalog, I get an error saying that there is no definition Query function for 'lyr'

Any help would be very appreciated, thanks!

Can you post the rest of your code? How you create the mxd object etc. You need to make a layer object to get the name property from it if you are running this out of process.

JonathanCusick

I've been trying to use a definition query just like rbasaria, and have used the syntax that clm42 suggested, but keep running into trouble. For now, I just want to be able to query out various zones from a larger map, and then run further geoprocessing functions on them.

import arcpy
from arcpy import env
import os

inzones = (.../RevisedGroundZones)
for lyr in inzones:
     if lyr.name == "RevisedGroundZones":
          lyr.definitionQuery = ' "Grounds_GSC_Zones" = \'2\''

When running it from the python window in ArcMap, I get an error message that says AttributeError: 'str' object has no attribute 'name' and when I run it as a script tool in ArcCatalog, I get an error saying that there is no definition Query function for 'lyr'

Any help would be very appreciated, thanks!

justinperez

sql = "OWNER='" + "NuStar GP Holdings LLC" + "'"

That's how I do it where Owner is the field name and NuStar GP Holdings LLC is the value...this is on SQL Server 2008..on Oracle I think it may differ a bit.

If I want to vary the value...notice how I dropped the double quotes around the value

sql = "OWNER='" + value_variable + "'"

Note the single quote after Owner= and also inside the last double quotes.

justinperez

sql = "OWNER=" "'" + "NuStar GP Holdings LLC" + "'"

That's how I do it where Owner is the field name and NuStar GP Holdings LLC is the value...this is on SQL Server 2008..on Oracle I think it may differ a bit. If I was going to use a variable for the value it would be something like..

sql = "OWNER=" "'" + value_variable + "'"

DaveVerbyla

rbasaria;87366 wrote:
Hi. Ive been trying to figure out how to accomplish this, but Ive run into some trouble. Im not sure if I am coding this correctly. My goal, overall, is to select features from a layer that intersect with features from another layer. The features in this second layer need to be defined by a definition query that I do not want to manually change. The selected features will then be exported to a new shapefile. I have been able to get the script to select and export, but I cannot get it to create the definition query. Any assistance will be greatly appreciated. Here is the

Query = "[DISTRICT] = '1' AND [COUNTY] = 'ORANGE'"
arcpy.MakeFeatureLayer(inFC,layername,Query)

ChrisMathers

Your welcome. All where clauses and definition queries need to be strings so sometimes escaping is needed to make it work.

RafiqBasaria

That worked perfectly. Thanks!! I was unaware that this could be done. Thank you.

ChrisMathers

Just escape some of hte quotes and wrap it in more quotes:

'"DISTRICT" = \'1\' AND "COUNTY" = \'ORANGE\''

RafiqBasaria

Been working on this more. I chaged a couple of parts, but am now running into the issue of how the layer is seeing the definition query. Since the layer reads it using single and double quotes, python cannot correctly enter the query. Here is the new code.

import arcpy

mxd = arcpy.mapping.MapDocument("C:/Rafiq_GIS/PYTHONTESTING/2011/2011_CountyComm/Data/DataPrep.mxd")

for lyr in arcpy.mapping.ListLayers(mxd):
    if lyr.name == "CountyCommission_TriCounty":
        lyr.definitionQuery = "DISTRICT" = '1' AND "COUNTY" = 'ORANGE'

arcpy.SelectLayerByLocation_management("BusRoutes_Ridership","INTERSECT","CountyCommission_TriCounty","","NEW_SELECTION")

arcpy.CopyFeatures_management("BusRoutes_Ridership","C:/Rafiq_GIS/PYTHONTESTING/2011/2011_CountyComm/Data/Orange/Or_Com1.shp")

As you can see, the lyr.definitionQuery = line is not correctly syntaxed. I cannot figure out how to get python to read that entire query definition as a string to insert into the definition query of the layer.

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.