ArcGIS 10 में "फील्ड मैपिंग"

मैंने एक पायथन स्क्रिप्ट लिखी है जो एक स्थानिक जुड़ाव और कुछ सरल गणनाएँ करती है। मेरी समस्या एक विशिष्ट क्षेत्र के लिए मर्ज नियम स्थापित करने और बाकी फ़ील्ड को छोड़ने के साथ है। उदाहरण के लिए, मेरे पास एक जनसंख्या क्षेत्र है जो स्थानिक स्थान से जुड़ने पर मर्ज नियम "प्रथम" का उपयोग करता है जो जनसंख्या गणना के पहले होने को पकड़ लेता है। मैं चाहता हूँ कि अन्य बहुभुज की स्थानिक सीमा में पाए जाने वाले सभी बहुभुजों के बीच जनसंख्या मूल्यों को जोड़ने के लिए "सुम" में मर्ज नियम स्थापित करने में सक्षम हो ।

मैंने फील्ड मैप्स और फील्ड मैपिंग ऑब्जेक्ट्स के साथ कुछ गहन छेड़छाड़ की है, लेकिन इसे ठीक से काम नहीं कर पा रहा हूं। विशेष रूप से मैंने विधि की कोशिश की है: mergeRule सेट करने के लिए popFieldMap.mergeRule = 'Sum' , लेकिन यह हमेशा "प्रथम" का सम्मान करता है।

कोई भी विचार कि मैं कैसे एक स्थानिक जुड़ाव में एक क्षेत्र के लिए मर्ज नियम को प्रोग्रामेटिक रूप से बदल सकता हूं?

धन्यवाद!

यहां मेरा कोड है (ध्यान रखें कि यह मेरे डेटा के लिए बहुत विशिष्ट है और इसमें स्क्रिप्ट के कुछ चरणों का परीक्षण करने के लिए लाइनें हैं):

import arcpy,sys,os

#Get the Files involved, set some variables.
SectorTable = sys.argv[1]
SectorShape = sys.argv[2]
MaxDev = sys.argv[3]
PopulationFC = sys.argv[4]
OutputFC = sys.argv[5]
DeviationField ="Angle_Deviation"
ID = "SectorID"
newID = "BP_ID"
mxd = arcpy.mapping.MapDocument('CURRENT')
df = arcpy.mapping.ListDataFrames(mxd)[0]

#Check to see if ID field types and name match
try:
    SectorShapeFields = arcpy.ListFields (SectorShape,ID)
    SectorTableFields = arcpy.ListFields (SectorTable,ID)
    arcpy.AddMessage("Finished Listing Fields")
    nameCheck = SectorShapeFields[0].name == SectorTableFields[0].name
    typeCheck = SectorShapeFields[0].type == SectorTableFields[0].type
except:
    arcpy.AddMessage("Failed to List Fields")

#If they do not match, add new fields to correct.
if not nameCheck:
    arcpy.AddMessage("Field Names do not match!  Adding new field to circumvent...")
if not typeCheck:
    arcpy.AddMessage("Field Types do not match!  Adding new field to circumvent...")
    if SectorShapeFields[0].type != SectorTableFields[0].type:
        try:
            arcpy.AddField_management(SectorShape, newID, SectorTableFields[0].type,10)
            arcpy.CalculateField_management(SectorShape, newID,"!"+ID+"!", "PYTHON")
            arcpy.RefreshTOC()
        except:
            arcpy.AddMessage("Error in Creating Field. Does it already exist?")


#Join the two together
arcpy.AddMessage("Performing Join")
arcpy.AddJoin_management( SectorShape, newID, SectorTable, ID)
arcpy.SelectLayerByAttribute_management (SectorShape,"NEW_SELECTION","Angle_Deviation>"+str(MaxDev))
df.zoomToSelectedFeatures()

#Field Mapping ...
myMap = arcpy.FieldMappings()
myMap.addTable(PopulationFC)
myMap.addTable(SectorShape)

#Verify the field merge rule for the pop10 field.
fIndex = myMap.findFieldMapIndex("POP10")
arcpy.AddMessage(str(fIndex))
popFieldMap = myMap.getFieldMap(fIndex)
arcpy.AddMessage(str(popFieldMap.mergeRule))
popFieldMap.mergeRule = 'Sum'
arcpy.AddMessage(str(popFieldMap.mergeRule))
popFieldMap2 = popFieldMap

##Test
fIndex = myMap.findFieldMapIndex("POP10")
arcpy.AddMessage(str(fIndex))
popFieldMap = myMap.getFieldMap(fIndex)
arcpy.AddMessage(str(popFieldMap.mergeRule))

#Perform Spatial Join
arcpy.AddMessage("Performing Spatial Join")
arcpy.SpatialJoin_analysis(SectorShape, PopulationFC, OutputFC,"JOIN_ONE_TO_ONE","",myMap,"INTERSECT")

#Add Result and Symbolize
arcpy.mapping.AddLayer(df,arcpy.mapping.Layer(OutputFC))
translayer = arcpy.mapping.ListLayers(mxd,"",df)[0]
translayer.transparency = 50

arcpy.RefreshActiveView()

EDIT - नीचे दिए गए समाधान के साथ कोड है!

import arcpy,sys,os

#Get the Files involved, set some variables.
SectorTable = sys.argv[1]
SectorShape = sys.argv[2]
MaxDev = sys.argv[3]
PopulationFC = sys.argv[4]
OutputFC = sys.argv[5]
DeviationField ="Angle_Deviation"
ID = "SectorID"
newID = "BP_ID"
mxd = arcpy.mapping.MapDocument('CURRENT')
df = arcpy.mapping.ListDataFrames(mxd)[0]

#Check to see if ID field types and name match
try:
    SectorShapeFields = arcpy.ListFields (SectorShape,ID)
    SectorTableFields = arcpy.ListFields (SectorTable,ID)
    arcpy.AddMessage("Finished Listing Fields")
    nameCheck = SectorShapeFields[0].name == SectorTableFields[0].name
    typeCheck = SectorShapeFields[0].type == SectorTableFields[0].type
except:
    arcpy.AddMessage("Failed to List Fields")

#If they do not match, add new fields to correct.
if not nameCheck:
    arcpy.AddMessage("Field Names do not match!  Adding new field to circumvent...")
if not typeCheck:
    arcpy.AddMessage("Field Types do not match!  Adding new field to circumvent...")
    if SectorShapeFields[0].type != SectorTableFields[0].type:
        try:
            arcpy.AddField_management(SectorShape, newID, SectorTableFields[0].type,10)
            arcpy.CalculateField_management(SectorShape, newID,"!"+ID+"!", "PYTHON")
            arcpy.RefreshTOC()
        except:
            arcpy.AddMessage("Error in Creating Field. Does it already exist?")


#Join the two together
arcpy.AddMessage("Performing Join")
arcpy.AddJoin_management( SectorShape, newID, SectorTable, ID)
arcpy.SelectLayerByAttribute_management (SectorShape,"NEW_SELECTION","Angle_Deviation>"+str(MaxDev))
df.zoomToSelectedFeatures()

#Field Mapping ...
myMap = arcpy.FieldMappings()
myMap.addTable(PopulationFC)
myMap.addTable(SectorShape)

#Verify the field merge rule for the pop10 field.
fIndex = myMap.findFieldMapIndex("POP10")
arcpy.AddMessage(str(fIndex))
popFieldMap = myMap.getFieldMap(fIndex)
arcpy.AddMessage(str(popFieldMap.mergeRule))
popFieldMap.mergeRule = 'Sum'
arcpy.AddMessage(str(popFieldMap.mergeRule))

myMap.replaceFieldMap(fIndex,popFieldMap)

##Test
fIndex = myMap.findFieldMapIndex("POP10")
arcpy.AddMessage(str(fIndex))
popFieldMap = myMap.getFieldMap(fIndex)
arcpy.AddMessage(str(popFieldMap.mergeRule))

#Perform Spatial Join
arcpy.AddMessage("Performing Spatial Join")
arcpy.SpatialJoin_analysis(SectorShape, PopulationFC, OutputFC,"JOIN_ONE_TO_ONE","",myMap,"INTERSECT")

#Add Result and Symbolize
arcpy.mapping.AddLayer(df,arcpy.mapping.Layer(OutputFC))
translayer = arcpy.mapping.ListLayers(mxd,"",df)[0]
translayer.transparency = 50

arcpy.RefreshActiveView()

arcgis-10.0 arcpy field-mapping

— पिक्सेल
स्रोत

क्या आप अपनी स्क्रिप्ट पोस्ट कर सकते हैं?

— blah238

ज़रूर, मुझे इसे पोस्ट करने दो। ध्यान रखें कि मेरे कोड के कुछ हिस्से हैं जिन्हें मैंने परीक्षण चरणों के रूप में डाला है।

— पिक्सेल

क्या आपने populationफ़ील्ड के डेटा प्रकार की जाँच की है ?

— एलेक्स मार्कोव

मेरा पहला झुकाव यह है कि संभवतः पॉप फ़ील्ड एक पाठ प्रकार है।

— ब्रैड नेसोम

फ़ील्ड प्रकार "लॉन्ग" है, इसलिए यह गणितीय योग करने के लिए योग्य होना चाहिए, सही?

— पिक्सेल

मुझे लगता है कि आपको FieldMappings.replaceFieldMapइसे बनाए रखने के लिए वास्तव में उपयोग करने की आवश्यकता है । उदाहरण के लिए मदद विषय से आउटपुट फ़ील्ड पर इनपुट फ़ील्ड मैपिंग :

# First get the TRACT2000 fieldmap. Then add the TRACTCODE field
#   from Blocks2 as an input field. Then replace the fieldmap within
#   the fieldmappings object.
#
fieldmap = fieldmappings.getFieldMap(fieldmappings.findFieldMapIndex("TRACT2000"))
fieldmap.addInputField(fc2, "TRACTCODE")
fieldmappings.replaceFieldMap(fieldmappings.findFieldMapIndex("TRACT2000"), fieldmap)

एक अन्य विचार यह है कि फील्ड मैप्स का परीक्षण करते समय, मैं एक स्क्रिप्ट टूल और कस्टम स्क्रिप्ट टूल व्यवहार का उपयोग करना पसंद करता हूं ताकि मैं उनके द्वारा बनाए गए लंबे लंबे स्ट्रिंग्स को डिकोड करने की कोशिश करने के बजाय फील्ड मैप का निरीक्षण कर सकूं।

उदाहरण ToolValidatorवर्ग मैं यह निष्कर्ष निकालता था कि हां, replaceFieldMapइसके लिए आवश्यक है कि इसे पैरामीटर मान में बरकरार रखा जाए:

class ToolValidator:
  """Class for validating a tool's parameter values and controlling
  the behavior of the tool's dialog."""

  def __init__(self):
    """Setup arcpy and the list of tool parameters."""
    import arcpy
    self.params = arcpy.GetParameterInfo()

  def initializeParameters(self):
    """Refine the properties of a tool's parameters.  This method is
    called when the tool is opened."""
    return

  def updateParameters(self):
    """Modify the values and properties of parameters before internal
    validation is performed.  This method is called whenever a parmater
    has been changed."""
    if (not self.params[0].hasBeenValidated
    or not self.params[1].hasBeenValidated):
        targetFeatures = self.params[0].value
        joinFeatures = self.params[1].value
        fieldMappings = arcpy.FieldMappings()
        if targetFeatures:
            fieldMappings.addTable(targetFeatures)
            if joinFeatures:
                fieldMappings.addTable(joinFeatures)
                idx = fieldMappings.findFieldMapIndex("PRIORITY")
                fieldMap = fieldMappings.getFieldMap(idx)
                fieldMap.mergeRule = 'Sum'
                fieldMappings.replaceFieldMap(idx, fieldMap) # if this line is commented out, the merge rule reverts to 'First'
        self.params[3].value = fieldMappings.exportToString()
    return

  def updateMessages(self):
    """Modify the messages created by internal validation for each tool
    parameter.  This method is called after internal validation."""
    return

— blah238
स्रोत

अच्छा लगा, मुझे इसे देखने दें और देखें कि मैं इसे अपनी स्क्रिप्ट में कहां लागू करूंगा। मैं अपने निष्कर्षों के साथ वापस आऊंगा

— पिक्सेल

यह एक महान जवाब था जिसने समस्या को हल किया। इसके अलावा, समस्या के कारण के रूप में महान विस्तार प्रदान करने के लिए धन्यवाद! :)

— पिक्सेल

ArcGIS 10 में "फील्ड मैपिंग" - ArcPy