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 clip map sheets to a Map Gird

JosephPlessis

I have multiple geo referenced scanned maps from the 1950's, Like 60 of them.

I need to clip some white space off the edges so I can apply a fishnet.

So the main question is I have a 50k nts grid where the maps sit perfectly within each grid tile, When i try to clip the map to the grid it does nothing because it is reading the outside border of the nts grid. how would i get the clip to read each individual grid tile for a clip?

thanks

Encuentra más publicaciones etiquetadas con

Accepted answers

All comments

JosephPlessis

THIS WORKS!

<SPAN class="keyword token">import</SPAN> arcpy
<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>
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"3D"</SPAN><SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

<SPAN class="comment token"># Set the current workspace</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> <SPAN class="string token">"Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\soil"</SPAN>
Dir <SPAN class="operator token">=</SPAN> <SPAN class="string token">"Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Mapsheets\\ClippedMaps"</SPAN>


GRID<SPAN class="operator token">=</SPAN><SPAN class="string token">"Q:\\josephplessis\\GIS 307 Resource Managment\\Project\\Geodatabases\\Grid_points.gdb\\GRID"</SPAN>

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


<SPAN class="comment token"># Get and print a list of TIFs from the workspace</SPAN>

rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> raster <SPAN class="keyword token">in</SPAN> rasters<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> raster
&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> Dir <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN><SPAN class="string token">"clp_"</SPAN><SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> name
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>RasterDomain_3d<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/temp"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"POLYGON"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: Central Feature</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>MeanCenter_stats<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory/temp"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/temp1"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Process: Select Layer By Location</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation_management<SPAN class="punctuation token">(</SPAN>GRIDlyr<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"intersect"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"in_memory/temp1"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"NEW_SELECTION"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Clip Raster</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Clip_management<SPAN class="punctuation token">(</SPAN>raster<SPAN class="punctuation token">,</SPAN>"<SPAN class="comment token">#", name,GRIDlyr, "#", "ClippingGeometry")</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># Delete in_memory</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Delete_management<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"in_memory"</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"processing "</SPAN> <SPAN class="operator token">+</SPAN> raster <SPAN class="operator token">+</SPAN> <SPAN class="string token">" complete..."</SPAN>


&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"All processing is now finished"</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>‍</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>

JosephPlessis

is there any good resources to script this

DarrenWiens2

Yes, you could certainly do this with a Python script, looping through rasters, calculating raster centroid from extent, selecting the grid cell covering the centroid, doing the clip.

JosephPlessis

I did selct a grid portion and created a layer form selcted data then proceeded to do the clip.

Your way sounds easier.

Is there a way to script this for I have a lot of maps to do.

DarrenWiens2

Have you tried select (the grid polygon), then clip?

Ú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.