Latest Contributions by CCWeedcontrol

So i have two layers and a script that uses arcpy.SpatialJoin_analysis to do a spatial join. Basically it spacial joins a zoning layer to parcels but upon further review it seems as i get weird results after the spatial join. In the "FZONE_CODE" field i get some that have 3 different zones. Which in some case can be the case not always and i find it through out spatial joined layer. The 1st example(image1) i have is the parcel is joined with the zoning and it should just be Ind(blue) but  if you look at it Yellow(Res) or Red(Com) doesn't go into the that parcel but the "FZONE_CODE" field states it's ComIndRes. Second example (image2) after the spatial join has a parcel with ResResRes but Res(yellow) doesn't go into that parcel either but some how it has ResResRes in the "FZONE_CODE" field. So i thought maybe i didn't snap to the end point or vertex correctly so i zoomed in as close as i can and i notice that vertex of the Zoning layer is not snapped to the vertex or end of the parcel  so i try to move the zoning  feature to the end/vertex of the parcel but it doesn't snap to it, it snaps off to the side(image3). So i am thinking this is why i get those weird results i am talking about. Ever time i try to move it just snaps the the side of the end or vertex, so my question is this normal? ... View more

I am not sure what's is going on here but if use arcpy.SpatialJoin_analysis and with INTERSECT as the match_option with two layers i get something that doesn't look right. The layer Fzon3 has zoning COM, IND or RES i do the spatial join to Taxparcels but the Spatial Join with INTERSECT doesn't look right, see picture. After the script is ran and i add the layers to arcmap and i query the FZon2 layer for "COM"(COM  is the red) you can see that the spatial join didn't spatial join correctly because there should be more parcels with "COM". Any help would be great thanks! #WORKING SCRIPT # Import arcpy module import sys, arcpy, traceback, os from arcpy import env import win32com.client from os import path as p from datetime import datetime as d startTime = d.now() # Local variables: arcpy.env.overwriteOutput = True arcpy.env.workspace = "C:\TEMP\Test.gdb" arcpy.env.qualifiedFieldNames = False if arcpy.Exists("C:\TEMP\Test.gdb\FZon2"):     arcpy.Delete_management("C:\TEMP\Test.gdb\FZon2") Tax = r"C:\TEMP\Test|Taxparcels" Tax1 = "C:\TEMP\Test.gdb" arcpy.FeatureClassToGeodatabase_conversion(Tax, Tax1) Padm1 = "C:\TEMP\Test.gdb" Taxparcels = "C:\TEMP\Test.gdb\Taxparcels" # Process: Make Feature Layer arcpy.MakeFeatureLayer_management(Taxparcels, "In_memory\TaxPar") TaxPar1 = "In_memory\TaxPar" arcpy.MakeFeatureLayer_management("In_memory\TaxPar", TaxPar1) FZon3 = "C:\TEMP\Test.gdb\FUTURE_LAND_USE_ZONING" def Layers0(TaxPar1, FZon3):     FieldMapString = "" \                 + """DXF_TEXT "DXF_TEXT" true true false 18 Text 0 0 ,First,#,""" + TaxPar1 + """,DXF_TEXT,-1,-1;"""\                 + """ACCOUNT "ACCOUNT" true true false 11 Text 0 0 ,First,#,""" + TaxPar1 + """,ACCOUNT,-1,-1;"""\                 + """ACRES "ACRES" true true false 4 Double 10 3 ,First,#,""" + TaxPar1 + """,ACRES,-1,-1;"""\                 + """Instrument "Instrument" true true false 10 Text 0 0 ,First,#,""" + TaxPar1 + """,Instrument,-1,-1;"""\                 + """SiteAddress "SiteAddress" true true false 106 Text 0 0 ,First,#,""" + TaxPar1 + """,SiteAddress,-1,-1;"""\                 + """SiteCity "SiteCity" true true false 32 Text 0 0 ,First,#,""" + TaxPar1 + """,SiteCity,-1,-1;"""\                 + """SiteZip "SiteZip" true true false 10 Text 0 0 ,First,#,""" + TaxPar1 + """,SiteZip,-1,-1;"""\                 + """SubName "SubName" true true false 20 Text 0 0 ,First,#,""" + TaxPar1 + """,SubName,-1,-1;"""\                 + """OwnerName "OwnerName" true true false 100 Text 0 0 ,First,#,""" + TaxPar1 + """,OwnerName,-1,-1;"""\                 + """OwnerSec "OwnerSec" true true false 64 Text 0 0 ,First,#,""" + TaxPar1 + """,OwnerSec,-1,-1;"""\                 + """Address "Address" true true false 64 Text 0 0 ,First,#,""" + TaxPar1 + """,Address,-1,-1;"""\                 + """City "City" true true false 32 Text 0 0 ,First,#,""" + TaxPar1 + """,City,-1,-1;"""\                 + """State "State" true true false 4 Text 0 0 ,First,#,""" + TaxPar1 + """,State,-1,-1;"""\                 + """ZipCode "ZipCode" true true false 10 Text 0 0 ,First,#,""" + TaxPar1 + """,ZipCode,-1,-1;"""\                 + """Legal "Legal" true true false 254 Text 0 0 ,First,#,""" + TaxPar1 + """,Legal,-1,-1;"""\                 + """Sec "Sec" true true false 2 Text 0 0 ,First,#,""" + TaxPar1 + """,Sec,-1,-1;"""\                 + """Twp "Twp" true true false 2 Text 0 0 ,First,#,""" + TaxPar1 + """,Twp,-1,-1;"""\                 + """Rng "Rng" true true false 2 Text 0 0 ,First,#,""" + TaxPar1 + """,Rng,-1,-1;"""\                 + """Quarter "Quarter" true true false 2 Text 0 0 ,First,#,""" + TaxPar1 + """,Quarter,-1,-1;"""\                 + """TaxCode "TaxCode" true true false 7 Text 0 0 ,First,#,""" + TaxPar1 + """,TaxCode,-1,-1;"""\                 + """FireDist "FireDist" true false false 64 Text 0 0 ,First,#,""" + TaxPar1 + """,FireDist,-1,-1;"""\                 + """HighwayDist "HighwayDist" true false false 74 Text 0 0 ,First,#,""" + TaxPar1 + """,HighwayDist,-1,-1;"""\                 + """SchoolDist "SchoolDist" true false false 64 Text 0 0 ,First,#,""" + TaxPar1 + """,SchoolDist,-1,-1;"""\                 + """Acres1 "Acres1" true false false 4 Double 10 3 ,First,#,"""+ TaxPar1 + """,Acres1,-1,-1;"""\                 + """SITE_ADDR "SITE_ADDR" true false false 50 Text 0 0 ,First,#,""" + TaxPar1 + """,SITE_ADDR,-1,-1;"""\                 + """FZONE_CODE "FZONE_CODE" true false false 10 Text 0 0 ,First,#,"""+ FZon3 + """,ZONE_CODE,-1,-1;"""\                                     fieldmappings = arcpy.FieldMappings()     fieldmappings.loadFromString(FieldMapString)     return fieldmappings arcpy.SpatialJoin_analysis(TaxPar1, FZon3, "C:\TEMP\Test.gdb\FZon2", "JOIN_ONE_TO_ONE", "KEEP_ALL",Layers0(TaxPar1, FZon3), "INTERSECT") try:     print '(Elapsed time: ' + str(d.now() - startTime)[:-3] + ')' except Exception, e:     # If an error occurred, print line number and error message     import traceback, sys     tb = sys.exc_info()[2]     print "Line %i" % tb.tb_lineno     print e.message ... View more