คำถามติดแท็ก python

ฉันเพิ่งติดตั้งระบบปฏิบัติการใหม่ทั้งหมดและติดตั้ง ActivePython 2.7.2 ฉันต้องการใช้การติดตั้ง Python นี้เพื่อใช้กับ ArcGIS 10.1 ฉันเพิ่งติดตั้ง ArcGIS 10.1 และติดตั้ง Python ใหม่ ฉันต้องทำอย่างไรเพื่อให้ Python ให้ ArcGIS ติดตั้งใน Python ปัจจุบันของฉัน

9 arcgis-desktop  python  arcgis-10.1  installation 

ฉันสนใจที่จะสร้างพื้นผิวแรสเตอร์ที่แสดงความเข้มแสงอาทิตย์เฉลี่ยต่อหน่วยพื้นที่สำหรับพล็อตสวนขนาดเล็ก ฉันนึกภาพอินพุตในโมเดลจะเป็นดังต่อไปนี้: พื้นที่แปลง - เพื่อเห็นแก่ความเรียบง่ายสมมติว่าสวนเป็นแบบวงกลม ความสูงของต้นไม้โดยเฉลี่ยของขอบป่าโดยรอบ ความเข้มของแสงอาทิตย์ในมุมต่าง ๆ ของดวงอาทิตย์ มุมดวงอาทิตย์เฉลี่ยที่สัมพันธ์กับพล็อตสวนในช่วงเวลาหนึ่ง ในที่สุดฉันต้องการแสดงตำแหน่งที่ดีที่สุดสำหรับการเจริญเติบโตของพืชในแปลงสวนที่ล้อมรอบด้วยต้นไม้ ใครช่วยชี้ทางแก้ปัญหาให้ฉันได้บ้าง? ฉันแน่ใจว่ามีคนทำสิ่งนี้แล้วแม้ว่าฉันจะไม่สามารถหาข้อมูลอ้างอิงได้ แก้ไข: ตัวอย่างเช่นรูปภาพที่แนบมาแสดงให้เห็นว่ามีพื้นที่สวนขนาด 1 เอเคอร์ที่ล้อมรอบด้วยกิ่งสนที่มีอายุมาก เพื่อให้การวิเคราะห์นี้ง่ายขึ้นสมมติว่าเราสนใจเฉพาะสุดยอดแสงอาทิตย์และความเข้มในช่วงฤดูทำสวนตั้งแต่มิถุนายนถึงตุลาคมเท่านั้น การตัดโค่นต้นไม้ไม่ใช่ทางเลือก!

9 arcgis-desktop  python  arcpy  raster  r 

ภาพที่แนบมาแสดงให้เห็นช่องว่างของป่าที่มีต้นสนสีแดงแสดงเป็นวงกลมและต้นสนสีขาวที่แสดงเป็นรูปกากบาท ฉันสนใจที่จะพิจารณาว่ามีความสัมพันธ์เชิงบวกหรือเชิงลบระหว่างต้นสนสองสายพันธุ์หรือไม่ (เช่นมีการเติบโตในพื้นที่เดียวกันหรือไม่) ฉันตระหนักถึง Kcross และ Kmulti ในแพ็คเกจ R spatstat อย่างไรก็ตามเนื่องจากฉันมี 50 ช่องว่างในการวิเคราะห์และคุ้นเคยกับการเขียนโปรแกรมใน python มากกว่า R ฉันจึงต้องการหาวิธีการวนซ้ำโดยใช้ ArcGIS และ python ฉันยังเปิดให้โซลูชั่น R ฉันจะใช้ฟังก์ชั่น K ของริบลีย์ bivariate ได้อย่างไร?

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

มีคำถามมากมายและคำตอบที่เกี่ยวข้องสำหรับวิธีการสร้างคุณลักษณะหลายส่วนเลเยอร์เป็นคุณสมบัติส่วนเดียวโดยใช้ QGIS: จะแบ่งรูปหลายเหลี่ยมเดียวออกเป็นหลายโพลีได้อย่างไร จะแยกคุณสมบัติที่ไม่ต่อเนื่องใน QGIS ได้อย่างไร? แต่โซลูชั่นทั้งหมดที่ฉันเห็นสร้างเลเยอร์ใหม่ทั้งหมดและจะแปลงรูปเรขาคณิตหลายส่วนทั้งหมดเป็นแบบเดี่ยว ดังนั้นคำถามคือจะมีวิธีการอย่างไรในช่วงเซสชั่นการแก้ไขเลือกคุณสมบัติหลายส่วนให้เป็นคุณสมบัติส่วนเดียวโดยไม่จำเป็นต้องสร้างเลเยอร์ใหม่ มันจะตรงกันข้ามกับเครื่องมือ "ผสานคุณสมบัติที่เลือก" และคล้ายกับการระเบิดใน ArcGIS มีปลั๊กอินทำเช่นนี้? วิธีนี้สามารถทำได้โดยใช้คอนโซลหลาม?

9 qgis  python  qgis-plugins  pyqgis 

ฉันมีรูปหลายเหลี่ยมสองอัน POLYGON ((- 5 -5, -5 5,5 5,5 -5, -5 -5)) POLYGON ((1 1, 1 7, 7, 7 7, 7 1, 1 1)) และต้องการแยกพวกเขาอย่างเท่าเทียมกัน ผลลัพธ์ควรแยกตามเส้นที่มีเครื่องหมายสีแดง ฉันใช้ qgis / python shapely และต้องการใช้ python

9 qgis  python  postgis-1.5  shapely 

ฉันกำลังเขียนเว็บแอปพลิเคชั่นที่เน้นข้อมูลที่ส่งผ่าน apache คำถามของฉันเกี่ยวกับวิธีจัดการการประมวลผลที่ดีที่สุดเนื่องจากมีหลายตัวเลือก ฉันมี OpenLayers / JQuery / Javascript, PostGIS / Postgresql (พร้อม pgsql), python / psycopg2, php ฐานข้อมูลมีแถวประมาณ 3 ล้านแถวและต้นแบบปัจจุบันทำงานดังนี้: ผู้ใช้คลิกที่จุดบนหน้าต่าง OpenLayers พิกัดถูกส่งเป็นการร้องขอ AJAX ไปยังฟังก์ชัน python บนเซิร์ฟเวอร์ ขณะนี้ใบสมัครของฉันไม่มีสัญชาติ psycopg2 ของ Python ใช้เพื่อเรียกโพรซีเดอร์ที่เก็บไว้ pgsql และชุดค่า WKT ที่ยอดเยี่ยม (และฟิลด์ข้อมูล) จะถูกส่งกลับไปยังโมดูลหลาม ฟิลด์ข้อมูลใช้เพื่อจัดหมวดหมู่เรคคอร์ด WKT ในไพ ธ อนดังนี้: ค่า WKT ทั้งหมดถูกจัดหมวดหมู่เป็นหนึ่งใน 5 กลุ่ม ประมาณ …

9 python  openlayers-2  postgis  postgresql  design 

มีฉากที่น่าจะเป็น 3 ฉากที่ฉันพยายามจับภาพในระยะใกล้ สถานีรถไฟใต้ดินแลกเปลี่ยนซึ่งมีสถานีใกล้เคียง 2 แห่งขึ้นไป นั่นคือสถานีที่มีปัญหาเชื่อมต่อ 2 เส้นทางหลักหรือมากกว่าและมีสถานีใกล้เคียง 2 แห่งขึ้นไป สถานีรถไฟใต้ดินเทอร์มินัลซึ่งมีสถานีใกล้เคียงเพียง 1 แห่ง นี่คือสถานีที่ส่วนท้ายของบรรทัด สถานีรถไฟใต้ดินแบบอินไลน์ซึ่งมีสถานีใกล้เคียง 2 แห่งซึ่งเป็นหนึ่งในแนวทางใดวิธีหนึ่ง ฉันพยายามคำนวณค่าที่อาจเรียกว่า "ระยะทางเฉลี่ยระหว่างสถานีใกล้เคียง" arcpy.GenerateNearTable_analysis()สามารถจัดการกับสองตัวเลือก: ระยะทางคุณลักษณะที่อยู่ใกล้และระยะห่างระหว่างคุณลักษณะทั้งหมด ใครบ้างมีวิธีที่ชาญฉลาดสำหรับการแก้ไขสถานการณ์เหล่านี้? โปรดทราบว่าแต่ละสถานีได้รับการกำหนดให้เป็น "Interchange", "Terminal" หรือ "Inline" ในตารางคุณสมบัติใต้ช่อง "StationType" ที่เพิ่ม: นี่คือรหัส psuedo บางส่วนตามคำแนะนำของ @ whuber ในความคิดเห็น ฉันยังไม่มีเวลาคิดออกดังนั้นถ้าใครอยากแทงมันคุณจะได้รับเครื่องหมายถูก! ;) ฉันดูห้องสมุดNetworkXและดูเหมือนว่าจะทำงานได้ตามที่ต้องการ รับกราฟ: A —― B ―― C ―― D | E …

9 python  arcpy  network 

ฉันต้องการต่อไปนี้เพื่อเลือก (ไฮไลต์) เลเยอร์ "โครงการ" ที่ตัดกันตารางแล้วดึงระเบียนจาก "กริด" ที่เลือก ฉันไม่ต้องการบันทึกเลเยอร์ หมายเหตุ: การเรียกใช้สคริปต์จากภายใน ArcMap ในหน้าต่าง Python Console import arcpy # Make a layer from the feature class arcpy.MakeFeatureLayer_management("project","project_lyr") # Then add a selection to the layer based on location to features in another feature class arcpy.SelectLayerByLocation_management("project_lyr","INTERSECT","Grid") รหัสนี้ทำงานและไม่ได้รับข้อผิดพลาด แต่ก็ไม่ได้เลือก (เน้น) คุณลักษณะใด ๆ ความคิดใด ๆ …

9 python  arcmap  select-by-attribute 

ฉันกำลังพยายามสร้างแรสเตอร์โดยใช้ไลบรารี python gdal และฉันได้ไปยังจุดที่ข้อมูลกำลังถูกส่งออก แต่ข้อมูลเอาต์พุตถูกพลิกบนแกน x ของจุดกำเนิด ฉันรู้ว่าฉันต้องมองข้ามบางสิ่งบางอย่าง แต่ฉันไม่สามารถเข้าใจได้ว่าฉันกำลังทำอะไรอยู่ ความคิดใด ๆ เมื่อสร้างแรสเตอร์ฉันจะตั้งค่า x / y บนซ้ายและอาเรย์ดูเหมือนจะได้รับการจัดทำดัชนีจากบนซ้ายและดำเนินการต่อไปจนถึงล่างขวา ในรหัสด้านล่างฉันเติมอาร์เรย์ด้วยค่าของแถว เมื่อพิมพ์อาเรย์ออกมามันจะเป็นดังนี้: [[ 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. 1. …

9 python  raster  gdal 

ฉันมีข้อมูลพัสดุตัดข้อมูลเส้น ในข้อมูลพัสดุมีพัสดุบางส่วนที่ไม่ตัดกัน ฉันจะหาทางโปรแกรมได้อย่างไรว่าพัสดุที่ไม่ตัดกันนั้นอยู่ทางด้านขวาหรือด้านซ้ายของเส้น ขอบคุณ

9 arcgis-10.0  arcobjects  python  arcpy  geoprocessing 

ฉันต้องการทราบว่าเป็นไปได้หรือไม่ที่จะใช้ Python web framework ที่มีอยู่เช่น Django / GeoDjango เป็นส่วนหน้าเพื่อทำงานกับ Geoservices REST API ของ Esri ฉันมองเห็นการพัฒนาเอ็นจิ้นฐานข้อมูลที่เข้าถึงจุดปลายทาง REST ใน JSON ในฐานะฐานข้อมูล "ตาราง" ผู้ดูแลระบบ Django และแบบฟอร์ม Django ได้รับการออกแบบมาเป็นอย่างดีและฉันคิดว่าพวกเขาสามารถทำหน้าที่เป็นอินเทอร์เฟซที่ยอดเยี่ยมสำหรับการเพิ่มการเปลี่ยนแปลงและการลบตารางเชิงพื้นที่และที่ไม่ใช่เชิงพื้นที่โดยเฉพาะอย่างยิ่งผู้ที่มีความสัมพันธ์ แก้ไข: ฉันคิดเพิ่มเติมเกี่ยวกับเรื่องนี้อาจไม่เกี่ยวกับการเปลี่ยนซอร์สโค้ดของ Django แต่แทนที่จะออกแบบเว็บแอปพลิเคชันเพื่อรองรับความต้องการ โมเดลจัดเก็บลิงก์บริการแผนที่และ Views จัดเตรียมฟังก์ชัน CRUD ผ่านการโทรไปยังปลายทาง REST

9 python  arcgis-server  rest  geodjango  geoservices-rest-api 

เพื่อนร่วมงานและฉันเพิ่งเปิดตัวโครงการโอเพ่นซอร์สที่เรียกว่าMapotecaซึ่งเป็นชุดของ Python Scripts ที่รวมเข้ากับการติดตั้ง GDAL / FWTools ซึ่งจะช่วยให้ไดเรกทอรีต่างๆสามารถค้นหาข้อมูลเชิงพื้นที่และแยกขอบเขตเรขาคณิตจำนวนองค์ประกอบ และผู้ใช้เขียนลงในไฟล์ csv เรากำลังทำงานในส่วนที่อธิบายไฟล์ DBF ที่เกี่ยวข้องและเราต้องการสร้างภาพขนาดย่อของแต่ละไฟล์รูปร่างคุณรู้วิธีสร้างภาพขนาดย่อเหล่านี้ใน jpeg หรือ png โดยใช้ Python หรือไม่?

9 python  gdal  ogr  jpg 

หมายเหตุ: ในขณะที่คำถามนี้มีคำตอบเคล็ดลับเพิ่มเติมใด ๆ สำหรับการปรับกระบวนการเคอร์เซอร์ให้เหมาะสมจะได้รับการชื่นชมอย่างมาก ฉันจะตรวจสอบการอัปเดตใด ๆ ปัจจุบันเจ้านายของฉัน (ที่ทำงานในอเวนิว) และฉัน (ทำงานใน Python) กำลังพยายามแก้ปัญหาเดียวกัน แต่เราทั้งคู่แก้ไขมันได้ แต่ความเร็วที่โซลูชันของเราดำเนินการคือ ... ไม่เชื่อมโยงกันเพื่อพูดอย่างน้อยที่สุด สิ่งที่กระบวนการสคริปต์ของเขาใน 2 ชั่วโมงสามารถขุดได้ถึง 6 ความแตกต่างที่แท้จริงของไวยากรณ์และการใช้งานในตรรกะนั้นมาจาก Bitmaps ของ 3.x และ Cursors 10.x เราทั้งคู่: 1) จัดเก็บค่าจากตารางที่ 1 2) ใช้ค่าเหล่านั้นเพื่อค้นหาแถวในตารางที่ 2 3) จัดเก็บค่าจากตารางที่ 2 สำหรับการแทรกลงในตารางที่ 3 เป็นแถวใหม่ ในสคริปต์ทั้งสองกระบวนการเหล่านี้จะเสร็จสมบูรณ์ในสองลูปซ้อนกัน ก่อนที่ฉันจะเริ่มขุดสู่โลกมหัศจรรย์ของการเพิ่มประสิทธิภาพโค้ดนี่เป็นเหตุการณ์ที่คาดหวังเมื่อเปรียบเทียบประสิทธิภาพสคริปต์ของ Avenue กับ Python หรือไม่ นี่ไม่ใช่ครั้งแรกที่สคริปต์ของเขามีประสิทธิภาพสูงกว่าฉันอย่างมากในแง่ของเวลาในการทำงานดังนั้นฉันอยากรู้ว่ามีบางสิ่งที่ฉันควรทราบก่อนที่ฉันจะตรึงตัวเองเพื่อสคริปต์ที่น่ากลัว นี่คือสคริปต์ของฉันมีบิตที่ไม่เกี่ยวข้อง: import arcpy …

9 arcgis-10.0  python  arcview-3 

ฉันต้องการตั้งค่าไฟล์ python เดียวที่มีตัวแปรที่มีที่ตั้งของแหล่งข้อมูลทั้งหมดของฉัน จากนั้นสคริปต์อื่นของฉันทั้งหมดนี้จะถูกใช้งานจากนั้นเมื่อแหล่งข้อมูลเปลี่ยนแปลงฉันต้องแก้ไขไฟล์เดียว Data_sources.py ของฉันจะมีลักษณะเช่นนี้: BC_BEC = "Database Connections\\BC.sde\\FOREST_BC.BC_Eco\\BC.BC_BEC" BC_TFL = "Database Connections\\Forest_BC.sde\\BC.BC_Admin\\BC.TFL_Boundary" และอื่น ๆ ... ฉันจะทำให้สคริปต์อื่น ๆ ทั้งหมดโหลดหรือเข้าถึงตัวแปรเหล่านี้ใน Data_sources.py ได้อย่างไร

9 python 

ฉันโพสต์คำถามเมื่อสัปดาห์ที่แล้วเกี่ยวกับการปรับแต่งคลาส ToolValidator และได้รับคำตอบที่ดีมาก ในการทำงานกับโซลูชันที่เสนอฉันได้สร้างโมดูลแบบกำหนดเองที่ดำเนินการเคียวรีบน db และจะถูกเรียกโดยทั้งชั้น ToolValidator (เพื่อให้ค่าสำหรับรายการดรอปดาวน์) และภายหลังในสคริปต์การประมวลผลทางภูมิศาสตร์ (เพื่อรับอื่น ๆ พารามิเตอร์ตามรายการที่เลือกในรายการแบบหล่นลง) อย่างไรก็ตามฉันไม่สามารถเรียกโมดูลที่กำหนดเองได้ในคลาส ToolValidator ฉันพยายามที่จะผนวกกับเส้นทางโดยไม่มีโชค เมื่อฉันพยายามใช้การเปลี่ยนแปลงเหล่านี้กับสคริปต์ฉันได้รับข้อผิดพลาดรันไทม์: [Errno 9] ตัวอธิบายไฟล์ไม่ถูกต้อง ถ้าฉันใส่ความคิดเห็นออกสายนำเข้าจะไม่มีข้อผิดพลาด sys.path.append('my_custom_module_directory') import my_custom_module พวกคุณหลายคนอาจถามว่าทำไมฉันไม่ใช้เครื่องมือที่กำหนดเองกับ ArcObjects เหตุผลก็คือผู้ใช้ของฉันไม่มีสิทธิ์ที่จำเป็นในการลงทะเบียน DLLs ใด ๆ บนคอมพิวเตอร์ของพวกเขา UPDATE: สิ่งนี้เกิดขึ้นกับฉันใน ArcGIS 10 น่าสนใจพอฉันเริ่มผนวกเส้นทางภายในฟังก์ชัน initialiazeParameters ของคลาส ToolValidator ถ้าฉันผนวกส่วนนอก (เช่นด้านบน) ของคลาส ToolValidator ทุกอย่างจะทำงานอย่างที่คาดไว้ sys.path.append('C:/Working/SomeFolder') import somescript -------->THIS WORKS class ToolValidator: …

9 arcgis-10.0  python  geoprocessing