I have two input files: - a feature class with a number of polygons - a raster I want to do an Extract by Mask for each of the polygons in my feature class. Due to my data structure, I sometimes do not have an overlap between one of the polygons and the raster. Polygon 3 in the picture below does not have an overlap with the raster. I already have a procedure that gives me the desired results, but I still have to do ExtractByMask for all of the polygons. All I do is not saving the temporary raster if the result is empty. Is there a way or a function that can check if a polygon and a raster intersect before I do the extract by mask procedure? Below is my - probably way too complicated - code. It does the job, but it is very slow. (The way I check if the raster is empty is by checking if the mean value and the standard Deviation is equal to 0) import arcpy arcpy . env . workspace = "my_workspace" raster = "my_raster" fc_polygons = "my_featureClass" # loop over polygons in the feature class fc_polygons flds = [ "OBJECTID" , "SHAPE@" ] with arcpy . da . SearchCursor ( fc_polygons , flds ) as curs : for row in curs : objectId = row [ 0 ] arcpy . AddMessage ( "Extracting by mask for polygon: {0}" . format ( objectId ) ) polygon = row [ flds . index ( "SHAPE@" ) ] # extracting for the polygon new_raster = arcpy . sa . ExtractByMask ( raster , polygon ) # if new_raster is empty discard the temporary raster, otherwise save it if new_raster . mean == 0 and new_raster . standardDeviation == 0 : arcpy . Delete_management ( new_raster ) arcpy . AddMessage ( "resulting raster is empty and is, therefore, discarded." ) else : outname = "ras_buf_{0}" . format ( objectId ) new_raster . save ( outname ) arcpy . AddMessage ( " - saved raster as {0}" . format ( outname ) ) Edit: Curtis suggested using the extent of the raster, but my raster is not continuous. Below is another picture of how I imagine an intermediary step to look like. Only the polygons that overlap with the raster are selected:

I would run the aggregate tool on your raster with maximum to make a coarse (ie small size) raster, convert this to a polygon. Select By Location against with polygon to select the polygons that overlap with your raster data cells. Note this workflow does require a Spatial Analyst license. (UPDATE fixed a typo in the code the guaranteed writable temp folder path is env.scratchFolder) from arcpy import env from arcpy . sa import * lyrPoly = arcpy . MakeFeatureLayer_management ( "my_featureClass" , "lyrPoly" ) # create an generalized extent polygon of your raster tmpAgg = os . path . join ( env . scratchFolder , "xagg.tif" ) tmpPoly = os . path . join ( env . scratchFolder , "xaggpoly.shp" ) agg = SetNull ( IsNull ( Aggregate ( "my_raster" , 5 , "MAXIMUM" , "EXPAND" ) ) , 1 ) agg . save ( tmpAgg ) arcpy . RasterToPolygon_conversion ( tmpAgg , tmpPoly ) # Select polygons that intersect the extent of the raster data areas arcpy . SelectLayerByLocation ( lyrPoly , "INTERSECT" , tmpPoly ) arcpy . Delete_management ( tmpAgg ) arcpy . Delete_management ( tmpPoly )

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; } } }

How to check for overlap: Raster and polygon

BennoPerchermeier

I have two input files:

- a feature class with a number of polygons

- a raster

I want to do an Extract by Mask for each of the polygons in my feature class. Due to my data structure, I sometimes do not have an overlap between one of the polygons and the raster. Polygon 3 in the picture below does not have an overlap with the raster.

I already have a procedure that gives me the desired results, but I still have to do ExtractByMask for all of the polygons. All I do is not saving the temporary raster if the result is empty. Is there a way or a function that can check if a polygon and a raster intersect before I do the extract by mask procedure?

Below is my - probably way too complicated - code. It does the job, but it is very slow.

(The way I check if the raster is empty is by checking if the mean value and the standard Deviation is equal to 0)

<SPAN class="keyword token">import</SPAN> arcpy

arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"my_workspace"</SPAN>
raster <SPAN class="operator token">=</SPAN> <SPAN class="string token">"my_raster"</SPAN>
fc_polygons <SPAN class="operator token">=</SPAN> <SPAN class="string token">"my_featureClass"</SPAN>

<SPAN class="comment token"># loop over polygons in the feature class fc_polygons</SPAN>
flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">"OBJECTID"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">]</SPAN>
<SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_polygons<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> curs<SPAN class="punctuation token">:</SPAN>
    <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> curs<SPAN class="punctuation token">:</SPAN>
        objectId <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
        arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Extracting by mask for polygon: {0}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>objectId<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
        polygon <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN>flds<SPAN class="punctuation token">.</SPAN>index<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
        <SPAN class="comment token"># extracting for the polygon</SPAN>
        new_raster <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa<SPAN class="punctuation token">.</SPAN>ExtractByMask<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> polygon<SPAN class="punctuation token">)</SPAN>
        <SPAN class="comment token"># if new_raster is empty discard the temporary raster, otherwise save it</SPAN>
        <SPAN class="keyword token">if</SPAN> new_raster<SPAN class="punctuation token">.</SPAN>mean <SPAN class="operator token">==</SPAN> <SPAN class="number token">0</SPAN> <SPAN class="operator token">and</SPAN> new_raster<SPAN class="punctuation token">.</SPAN>standardDeviation <SPAN class="operator token">==</SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN>
            arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>new_raster<SPAN class="punctuation token">)</SPAN>
            arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"resulting raster is empty and is, therefore, discarded."</SPAN><SPAN class="punctuation token">)</SPAN>
        <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
            outname <SPAN class="operator token">=</SPAN> <SPAN class="string token">"ras_buf_{0}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>objectId<SPAN class="punctuation token">)</SPAN>
            new_raster<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">)</SPAN>
            arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN><SPAN class="string token">" - saved raster as {0}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>outname<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Edit: Curtis suggested using the extent of the raster, but my raster is not continuous. Below is another picture of how I imagine an intermediary step to look like. Only the polygons that overlap with the raster are selected:

Vind meer berichten getagd met

Accepted answers

curtvprice

I would run the aggregate tool on your raster with maximum to make a coarse (ie small size) raster, convert this to a polygon. Select By Location against with polygon to select the polygons that overlap with your raster data cells. Note this workflow does require a Spatial Analyst license.

(UPDATE fixed a typo in the code the guaranteed writable temp folder path is env.scratchFolder)

<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>

lyrPoly <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_featureClass"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lyrPoly"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># create an generalized extent polygon of your raster</SPAN>
tmpAgg <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>scratchFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xagg.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
tmpPoly <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>scratchFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xaggpoly.shp"</SPAN><SPAN class="punctuation token">)</SPAN>
agg <SPAN class="operator token">=</SPAN> SetNull<SPAN class="punctuation token">(</SPAN>IsNull<SPAN class="punctuation token">(</SPAN>Aggregate<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_raster"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"EXPAND"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
agg<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToPolygon_conversion<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># Select polygons that intersect the extent of the raster data areas</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation<SPAN class="punctuation token">(</SPAN>lyrPoly<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"INTERSECT"</SPAN><SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpPoly<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

BennoPerchermeier

Awesome thanks! This idea works. The aggregation is relatively slow but still faster than trying to do an Extract by mask on all the empty polygons.

I had to amend your script to make it work, so here it is:

<SPAN class="keyword token">import</SPAN> os
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

raster <SPAN class="operator token">=</SPAN> <SPAN class="string token">"my_raster"</SPAN>
fc_polygons <SPAN class="operator token">=</SPAN> <SPAN class="string token">"my_featureClass"</SPAN>
tempFolder <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Path\to\temporary\folder'</SPAN>

<SPAN class="comment token"># create temporary layer</SPAN>
lyrPoly <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN>fc_polygons<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lyrPoly"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># create an generalized extent polygon of your raster</SPAN>
tmpAgg <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>tempFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xagg2.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
tmpPoly <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>tempFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xaggpoly.shp"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># create aggregated raster with larger cell sizes</SPAN>
agg <SPAN class="operator token">=</SPAN> SetNull<SPAN class="punctuation token">(</SPAN>IsNull<SPAN class="punctuation token">(</SPAN>Aggregate<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"EXPAND"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
agg<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># convert to polygon feature class </SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToPolygon_conversion<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># Select polygons that intersect the extent of the raster data areas</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation_management<SPAN class="punctuation token">(</SPAN>lyrPoly<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"INTERSECT"</SPAN><SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpPoly<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I had to specify a temporary folder where I could save the temporary raster and the temporary polygon. For one, because I could not use the command "env.scratchRaster". For the other, because I could not save the tiff file to the workspace, because can't be saved to GDBs as far as I know.

Anyways, thank you very much for your answer, Curtis!!

curtvprice

(UPDATE fixed a typo in the code the guaranteed writable temp folder path is env.scratchFolder)

<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN> <SPAN class="operator token">*</SPAN>

lyrPoly <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_featureClass"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lyrPoly"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="comment token"># create an generalized extent polygon of your raster</SPAN>
tmpAgg <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>scratchFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xagg.tif"</SPAN><SPAN class="punctuation token">)</SPAN>
tmpPoly <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>scratchFolder<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"xaggpoly.shp"</SPAN><SPAN class="punctuation token">)</SPAN>
agg <SPAN class="operator token">=</SPAN> SetNull<SPAN class="punctuation token">(</SPAN>IsNull<SPAN class="punctuation token">(</SPAN>Aggregate<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_raster"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"MAXIMUM"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"EXPAND"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>
agg<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>RasterToPolygon_conversion<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># Select polygons that intersect the extent of the raster data areas</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation<SPAN class="punctuation token">(</SPAN>lyrPoly<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"INTERSECT"</SPAN><SPAN class="punctuation token">,</SPAN> tmpPoly<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpAgg<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN>tmpPoly<SPAN class="punctuation token">)</SPAN>
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

BennoPerchermeier

Thank you for your answer Curtis! I like your idea because it seems much faster than what I am doing. Though I did not explain the issue properly. My raster is actually not continuous and there are empty polygons in the middle.

Running the script above selects all of the polygons:

I will update my question accordingly. Do you have another idea, how to do this?

curtvprice

I suggest using the extent of the raster as a polygon to select those that overlap before you start processing the polys. Then use the layer (with just the polygons that overlap) for your cursor loop.

lyrPoly <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_featureClass"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"lyrPoly"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># make a polygon from extent</SPAN>
ext <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"my_raster"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>extent
array <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
array<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMin<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMin<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
array<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMin<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMax<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
array<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMax<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMax<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
array<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMas<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMin<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
array<SPAN class="punctuation token">.</SPAN>add<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>ext<SPAN class="punctuation token">.</SPAN>XMin<SPAN class="punctuation token">,</SPAN> ext<SPAN class="punctuation token">.</SPAN>YMin<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
extPoly <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polygon<SPAN class="punctuation token">(</SPAN>array<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token"># Select polygons that intersect the extent of the raster</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation<SPAN class="punctuation token">(</SPAN>lyrPoly<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"INTERSECT"</SPAN><SPAN class="punctuation token">,</SPAN> extPoly<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

http://desktop.arcgis.com/en/arcmap/latest/analyze/python/using-geometry-objects-with-geoprocessing-tools.htm

Doe mee aan het gesprek

Aangemelde leden kunnen berichten plaatsen, updates volgen en meer. Nieuw hier? Registreer een gratis account.

Snelle Links

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