การเลือกรูปหลายเหลี่ยม ArcSDE ทีละจุดใน ArcGIS Desktop โดยใช้ ArcPy?

ฉันคิดอยู่เสมอว่าจะต้องมีบางอย่างหายไป แต่ดูเหมือนว่าจะไม่มีเครื่องมือใน ArcGIS 10 ที่จะเลือกคุณสมบัติ (โดยเฉพาะรูปหลายเหลี่ยม) จากเลเยอร์ที่ตำแหน่ง (X, Y) ผ่าน ArcPy พารามิเตอร์สำหรับเครื่องมือเช่นนี้จะเป็นชื่อเลเยอร์และตำแหน่ง XY เท่านั้น

ในขณะนี้ฉันแก้ไขปัญหานี้โดยการสร้างคุณลักษณะระดับจุดที่มีจุดและดำเนินการ SelectLayerByLocation บนมัน อย่างไรก็ตามเมื่อคลาสคุณลักษณะรูปหลายเหลี่ยมอยู่ใน Oracle (เข้าถึงผ่าน ArcSDE 9.x) และมี 3.5 ล้านรูปหลายเหลี่ยมเวลาที่ใช้ในการเลือกอาจมากกว่า 5 นาทีเมื่อฉันคิดว่าสองหรือสอง (ด้วยรหัสน้อยกว่า) จะเป็น เหมาะสมกว่า คลาสคุณลักษณะมีดัชนีเชิงพื้นที่และฉันได้ลองใช้ arcpy.env.extent (ซึ่ง SelectLayerByLocation ดูเหมือนจะไม่สนใจ) เพื่อ จำกัด การเข้าถึงพื้นที่ทางภูมิศาสตร์ แต่ประสิทธิภาพยังคงแย่มาก

มีวิธีที่รวดเร็วกว่าในการทำเช่นนี้โดยใช้ ArcGIS Desktop 10 และ ArcPy?

อีกวิธีในการดำเนินการนี้คือการใช้เครื่องมือ Spatial Join ใช้จุดเป็นเลเยอร์คุณลักษณะการป้อนข้อมูลของคุณตามด้านบนและเลเยอร์รูปหลายเหลี่ยมเป็นคุณสมบัติตัวตนของคุณ
ซึ่งแตกต่างจาก SelectLayerByLocation, SpatialJoin ไม่เคารพสภาพแวดล้อมขอบเขต

targetlayer = layername
joinlayer=arcpy.PointGeometry(arcpy.Point(x, y))
fieldmappings = arcpy.FieldMappings()
fieldmappings.addTable(targetlayer)
arcpy.SpatialJoin_analysis(targetlayer, joinlayer, outputlayer, "JOIN_ONE_TO_MANY", "KEEP_COMMON", fieldmappings)

JOIN_ONE_TO_MANY อาจดูเหมือนเคาน์เตอร์ง่าย แต่เนื่องจากคุณมีเพียงหนึ่งคุณสมบัติการเข้าร่วมฟังก์ชั่นหลักของตัวเลือกนี้คือการปิดการรวมและกฎการรวม KEEP_COMMON จะตรวจสอบให้แน่ใจว่าผลลัพธ์ของคุณถูก จำกัด เฉพาะรูปหลายเหลี่ยมที่ตัดกันจุดของคุณ การแมปฟิลด์จะ จำกัด แอตทริบิวต์เอาท์พุทให้กับรูปร่างและคุณลักษณะของชั้นรูปหลายเหลี่ยมเท่านั้น ค่าเริ่มต้นจะรวมถึงคุณลักษณะของเลเยอร์ของจุดด้วย

ส่วนที่เหลือของค่าเริ่มต้นจะทำงานได้ดีดังนั้นคุณสามารถละทิ้งอาร์กิวเมนต์ที่เหลือได้

ฉันเพิ่งรู้อะไรบางอย่าง ...

หากคุณกำลังใช้สิ่งนี้เพื่อดำเนินการบริการประมวลผลทางภูมิศาสตร์คุณสามารถใช้งานบริการคุณลักษณะด้วยรูปหลายเหลี่ยมเลเยอร์ของคุณแทนและใช้การดำเนินการค้นหาบนบริการคุณลักษณะ

คุณสามารถใช้ตำแหน่ง XY แบบง่าย ๆ ได้พร้อมกับการทำงานของ Service Service Query รวมถึงควบคุมแอตทริบิวต์เอาต์พุตรวมถึงรูปร่าง

คุณอาจมีข้อ จำกัด เกี่ยวกับจำนวนที่คุณสามารถเปิดเผยคลาสรูปหลายเหลี่ยมได้ แต่หากคุณสามารถเข้าถึงบริการการประมวลผลทางภูมิศาสตร์ได้แล้วคุณควรสร้างบริการคุณลักษณะได้เช่นกัน

คำตอบนี้มาจากเก่าฟอรั่ม ArcGIS อภิปราย

ขอบคุณ Jason Scheirer สำหรับโค้ดที่กระชับกว่านี้:

SelectLayerByLocation(in_layer=arcpy.PointGeometry(arcpy.Point(x, y)), select_features="mylayer") 

และโดยเฉพาะ Chris Snyder สำหรับเคล็ดลับการแสดง:

การทำงานที่เร็วขึ้นอาจเป็นการบัฟเฟอร์จุดของคุณเล็กน้อยจากนั้นใช้ขอบเขตบัฟเฟอร์เป็นขอบเขตการวิเคราะห์เพื่อทำสำเนา in_memory (เครื่องมือ CopyFeatures) ของข้อมูล SDE ของคุณจากนั้นเลือก SelectByLocation บนชุดข้อมูล in_memory ที่มีขนาดเล็กและท้องถิ่น ด้วยวิธีการที่คุณจัดทำเครื่องมือ SelectByLocation ตามสภาพแวดล้อมขอบเขตการวิเคราะห์ซึ่งปกติแล้วจะไม่ทำ BTW: คุณสมบัติใด ๆ ที่ทับซ้อนกับขอบเขตการวิเคราะห์จะถูกคัดลอกด้วยเครื่องมือ CopyFeatures ฉันต้องการให้เครื่องมือ SelectByLocation และเมธอดเคอร์เซอร์ให้เกียรติขอบเขตการวิเคราะห์ ...