Simple ModelBuilder

MattSeitz

Hi all,

I'm trying to create a process in modelbuilder that will allow me to load .shp into an existing sde database.

I'm sure I need to append the new .shp to the existing database, but I'm just not sure exactly how to go about it in modelbuilder. I also need to create a loop or make a condition where the .shp are in the correct projection. As they are sometimes received in the wrong or with no projection at all.

 

Thanks!

MicahBabinski

Hi Matt,

To convert from shapefile to SDE feature class I'd use the Feature Class To Feature Class tool.

The projection step introduces another layer of complexity. I'm not 100% sure how to check the spatial reference of a shapefile in ModelBuilder. This is one thing arcpy is great at. If you haven't tried Python/arcpy, I'd encourage you to give it a shot. Assuming your desired spatial reference is NAD 1983 UTM Zone 10N, your script might look something like this:

<SPAN class="comment token"># import required modules</SPAN>
<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> os

<SPAN class="comment token"># create variables pointing to the desired spatial reference, a folder containing your input shapefiles and your output SDE workspace</SPAN>
spRef <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SpatialReference<SPAN class="punctuation token">(</SPAN><SPAN class="number token">26910</SPAN><SPAN class="punctuation token">)</SPAN>
shapefile_folder <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\WorkArea\ShapefilesFromContractor"</SPAN>
sde_workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"C:\WorkArea\DatabaseConnections\MySDEWorkspace.sde"</SPAN>


<SPAN class="comment token"># list all of the shapefiles in the shapefile folder</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> shapefile_folder
shpList <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFeatureClasses<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>


<SPAN class="comment token"># loop through the shapefiles</SPAN>
<SPAN class="keyword token">for</SPAN> shp <SPAN class="keyword token">in</SPAN> shpList<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># test to see if the projection is undefined</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>spatialReference<SPAN class="punctuation token">.</SPAN>GSCName <SPAN class="operator token">==</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># define the projection</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>DefineProjection_management<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SpatialReference<SPAN class="punctuation token">(</SPAN><SPAN class="number token">4269</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># project to the desired spatial reference</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Project_management<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">,</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>sde_workspace<SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>baseName<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> spRef<SPAN class="punctuation token">)</SPAN>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># if it not undefined...</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># test to see if it is in the wrong projected coordinate system</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>spatialReference<SPAN class="punctuation token">.</SPAN>PCSCode <SPAN class="operator token">!=</SPAN> <SPAN class="number token">26910</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># project the shapefile into the correct spatial reference, outputting to sde</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>Project_management<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">,</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>sde_workspace<SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>baseName<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> spRef<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># otherwise, just move it to sde</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>FeatureClassToFeatureClass_conversion<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">,</SPAN> sde_workspace<SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>baseName<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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Of course, if the shapefiles come in with the spatial reference undefined, you may need to guess, sleuth out, or ask what geographic coordinate system was used to collect the data and possibly use a transformation. For instance, GPS data would likely have a geographic coordinate system of WGS 1984.

Here's the help for the geoprocessing tools that will get you there:

Feature Class To Feature Class

Project

Define Projection

Good luck! Hope this helps.

Micah