ระบบสารสนเทศทางภูมิศาสตร์

ภายในตารางคุณลักษณะฉันพยายามแทนที่สตริงหนึ่งด้วยอีกสตริงหนึ่งโดยใช้ฟังก์ชัน 'replace' ตัวอย่างเช่นมีคุณสมบัติหลายอย่างที่ฟิลด์ 'NAME' มีสตริง '(B)' และฉันต้องการแทนที่ด้วยคำว่า 'County' ฉันกำลังเลือกช่องทำเครื่องหมาย 'อัปเดตฟิลด์ที่มีอยู่' และใช้นิพจน์ต่อไปนี้ replace( '(B)','(B)','County') ผลลัพธ์ที่ได้คือฟิลด์ 'NAME' สำหรับทุกคุณสมบัติจะถูกแทนที่ด้วย 'County' โดยไม่คำนึงว่าฟิลด์ 'NAME' นั้นมีอยู่ในตอนแรกหรือไม่นั้นมีสตริง '(B)' ความช่วยเหลือใด ๆ ที่จะได้รับการชื่นชม

12 qgis  field-calculator 

เรามีโครงการ QGIS ที่มีสัญลักษณ์มากมายหลายบรรทัด ข้อมูลมาจากฐานข้อมูลทางภูมิศาสตร์ที่สามารถเข้าถึงได้ทั้ง QGIS และ ArcGIS ตอนนี้เราต้องการใช้สัญลักษณ์ของโครงการ QGIS เพื่อจัดรูปแบบคุณสมบัติในเอกสารแผนที่ ArcGIS วิธีแรกของฉันคือใช้ SLD รูปแบบ QGIS สามารถส่งออกไปยัง SLD ได้อย่างง่ายดาย แต่คุณจะจัดรูปแบบคุณสมบัติตาม SLD โดยใช้ ArcGIS Desktop ได้อย่างไร การตั้งค่า WMS แบบง่ายนั้นไม่เพียงพอทั้ง QGIS และ ArcGIS Desktop ต้องการความสามารถในการแก้ไขคุณสมบัติ ถ้าใครสามารถจัดเตรียม SLD ด้วย WFS (ในแบบ Esri-readable) ที่จะเป็นไปได้เช่นกัน ไม่พบสิ่งใดเลย

12 qgis  arcgis-desktop  symbology  sld  interoperability 

ฉันได้ลองใช้ IRC เพื่อดูว่าผู้คนสามารถช่วยฉันได้หรือไม่ แต่สื่อดูเหมือนจะเงอะงะเล็กน้อยสำหรับการตอบคำถามที่ชัดเจนดังนั้นฉันจะลองทำแทน ฉันต้องการตั้งโปรแกรมแอปพลิเคชันแบบสแตนด์อโลนโดยใช้ PyQGIS วิธีที่ฉันสามารถโปรแกรมแอปพลิเคชันแบบสแตนด์อโลนโดยใช้ ArcPy PyQGIS นั้นตรงไปตรงมา ฉันไม่มีปัญหาในการใช้งานสคริปต์ Python ภายในเอดิเตอร์คอนโซล QGIS อย่างไรก็ตามเพื่อที่จะใช้การผูก Python ในสคริปต์แบบสแตนด์อโลนฉันแนะนำให้ฉันต้องเปลี่ยนการตั้งค่าสภาพแวดล้อมโดยใช้ชุดสคริปต์ (ฉันใช้ Windows Vista ที่บ้านดังนั้นเราจึงกำลังพูดถึงการเขียนสคริปต์ชุดสำหรับ Windows) สิ่งนี้ส่งผลให้เกิดความยุ่งเหยิงในตอนท้ายของฉันเนื่องจากฉันไม่เข้าใจการเขียนสคริปต์แบบชุดและฉันไม่ชอบความคิดที่ว่าฉันกำลังรีเซ็ตสภาพแวดล้อมที่ต้องเป็นวิธีหนึ่งสำหรับแอปพลิเคชันอื่นที่ฉันใช้ ฉันมีปัญหาในการทำความเข้าใจว่าทำไมคำแนะนำในการสร้างการตั้งค่าเหล่านี้จึงมีความซับซ้อนเกินความจำเป็นและเหตุใดพวกเขาจึงไม่ถึงจำนวนหนึ่งบรรทัดของรหัสว่า "ไปยังที่ที่มีการเชื่อมโยงทำให้ไลบรารีที่ฉันต้องการ ) แน่นอนว่า ArcPy ไม่ต้องการให้ผู้คนยุ่งกับการตั้งค่าสภาพแวดล้อมของคอมพิวเตอร์ดังนั้นฉันจึงมีปัญหาในการเข้าใจว่าทำไม PyQGIS ทำเช่นนั้นดูเหมือนว่าฉันจะยอมแพ้แทนที่จะสอนภาษา Windows Batch Script ให้ตัวเอง ตอนนี้ฉันสังเกตเห็นว่าเมื่อฉันเรียกใช้สคริปต์ในตัวแก้ไขคอนโซลของ QGIS มันทำงานได้ ได้รับมันยังทำให้ส่วนที่เหลือของ QGIS ใช้งานไม่ได้หลังจากที่สคริปต์ทำงานดังนั้นฉันต้องรีสตาร์ท QGIS อีกครั้งหลังจากนั้น แต่สคริปต์ทำงานตรงตามที่ฉันได้ตั้งโปรแกรมให้ทำงาน แต่เมื่อฉันทำนอกเหนือจากแอปพลิเคชัน QGIS มันไม่ได้ ปัญหาดูเหมือนว่าฉันไม่ได้นำเข้าแพคเกจ Python บางอย่างที่นำเข้าโดยอัตโนมัติโดยแอปพลิเคชัน …

12 qgis  pyqgis 

อะไรคือความแตกต่างหลักระหว่างความแตกต่างและเครื่องมือการประมวลผลทางภูมิศาสตร์คลิปใน QGIS?

12 qgis  clip  qgis-processing  differences 

ฉันต้องการแปลง Raster เป็น Point Vector ใน QGIS เป็นไปได้ไหม? QGIS มีตัวเลือกในการแปลงเป็นรูปหลายเหลี่ยม แต่ฉันไม่พบเครื่องมือใด ๆ ที่จะแปลงเป็นจุด ใครช่วยฉันได้บ้าง

12 qgis  raster 

ฉันได้เผยแพร่เลเยอร์ต่างๆโดยใช้ Geoserver 2.2.x พร้อมที่เก็บข้อมูล PostGIS 2.0.1 (PostgreSQL 9.1.x) ฟิลด์ที่มีค่าวันที่เท่านั้นและไม่ใช่ค่าที่เกี่ยวข้องกับเวลา / การประทับเวลาจะถูกเก็บไว้เป็นวันที่ หลังจากเผยแพร่เลเยอร์ของฉันเป็น wfs / wms และเปิดตารางในวันที่เบราว์เซอร์ vlues แสดงขึ้นพร้อมกับ 'Z' ต่อท้าย นี่เป็นผลิตภัณฑ์ของ api data store ของ geotools ที่ geoserver ใช้อยู่หรือไม่? ฉันคิดว่า api data store ของ geotools กำลังโยนผ้าห่มเหนือวันที่อ่านข้อมูลในรูปแบบ datetime, timestamp เป็นต้นและ 'z' จะปรากฏขึ้นตามเวลาของ zulu ฉันสามารถปิดโดยสิ้นเชิงกับสมมติฐานทางภูมิศาสตร์ของฉัน แต่ใน PostgreSQL 'z' หมายถึงซูลูหรือ utc หากค่าแหล่งข้อมูลถูกเก็บไว้เป็นวันที่และไม่มีส่วนประกอบเวลาทำไมจึงมี 'Z' …

12 postgis  geoserver  postgresql  postgis-2.0  geotools 

ฉันมีชุดแรสเตอร์ (.vrt) พร้อมข้อมูลความชื้นในดินทุกวัน ฉันต้องการที่จะรวมพิกเซลของแรสเตอร์ทั้งหมดเพื่อให้มีการวัดตามเดือน อย่างไรก็ตามโลกทั้งใบไม่ได้ครอบคลุมในแต่ละวันซึ่งส่งผลให้เกิดค่านิยมในสถานที่ที่มีการวัดสำหรับวันอื่น ๆ สิ่งที่ฉันต้องการทำคือผลรวมของค่าแรสเตอร์แต่ละอัน อย่างไรก็ตามดูเหมือนว่าในแต่ละครั้งที่พิกเซล nodata อยู่ในผลรวมพิกเซลที่ได้จะถูกจัดประเภทเป็นโนดาต้าโดยตรง ฉันอยากจะมีสิ่งที่ตรงกันข้าม: ไม่สนใจคุณค่าของสิ่งที่มีอยู่ทั้งหมดและรวมส่วนที่เหลือ ฉันคิดถึงวิธีแก้ปัญหา 2 วิธี: การรวม rasters โดยไม่สนใจค่าของโหนด การแปลงพิกเซลโนดาต้าให้เป็นค่า 0 จากนั้นรวมแรสเตอร์ทั้งหมด น่าเสียดายที่ฉันไม่พบเครื่องมือใด ๆ ในการทำเช่นนี้ มีใครช่วยฉันบ้าง

12 raster  gdal  qgis-2 

ฉันกำลังอ่าน shapefile โดยใช้ฟิโอน่า: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] สิ่งนี้ให้ฉันตามที่มีรายการของพิกัดและพิมพ์ 'MultiPolygon' อย่างไรก็ตามฉันไม่แน่ใจว่าจะใช้ข้อมูลพิกัดเพื่อสร้าง Shapely MultiPolygon ได้อย่างไรเนื่องจากต้องการ a sequence of exterior ring and hole list tuples: [((a1, ..., aM), [(b1, ..., bN), ...]), ...]. มีวิธีที่ 'ถูกต้อง' ในการทำเช่นนี้โดยใช้ Fiona และ Shapely หรือไม่?

12 python  shapely  fiona 

ฉันต้องการใช้ ArcGIS Desktop เพื่อรวมรูปหลายเหลี่ยมที่อยู่ติดกัน ไม่มีเขตข้อมูลทั่วไประหว่างรูปหลายเหลี่ยมและดังนั้นเครื่องมือการละลายไม่ได้เป็นทางออกที่ดีที่สุด รูปหลายเหลี่ยมที่ฉันต้องการผสานเข้าด้วยกันนั้นอาจจะอยู่ติดกันหรืออยู่ในระยะที่ห่างกัน ฉันต้องการหลีกเลี่ยงการใช้ ArcObjects อย่างไรก็ตามยินดีต้อนรับสคริปต์ Python เครื่องมือการละลายผสานรูปหลายเหลี่ยมที่อยู่ติดกันเมื่อแอตทริบิวต์ที่จะละลายถูกปล่อยออกมา (และสร้างหลายส่วนถูกปิดใช้งาน) อย่างไรก็ตามมันไม่ได้มีการตั้งค่าสำหรับการยอมรับเชิงพื้นที่ นอกจากนี้การละลายไม่สามารถใช้ได้กับรูปหลายเหลี่ยมที่มีจุดยอดเหมือนกัน ฉันได้ลองวิธีการด้วยเครื่องมือการรวมตามที่กล่าวไว้ในหนึ่งในคำตอบ ปัญหาคือว่ารูปร่างของคุณสมบัติสุดท้ายที่จะถูกสร้างขึ้นโดยการรวมคุณสมบัติหลายอย่างไม่เหมือนกับรูปร่างที่จะเกิดขึ้นหากขอบเขตระหว่างคุณสมบัติถูกละลาย

12 arcgis-desktop  arcgis-10.0  merge  dissolve  adjacency 

ฉันมีตารางจุดที่คอลัมน์รูปทรงเรขาคณิตมีค่าทั้ง 3 มิติและ 2 มิติฉันต้องการแปลงรูปทรงเรขาคณิต 3 มิตินี้เป็น 2 มิติมีฟังก์ชั่นใดบ้างใน PostGIS โปรดช่วยฉันแก้ปัญหานี้ด้วย ขอแสดงความนับถือ Sreesha.TA

12 postgis  dimensions 

ฉันใช้ Python พยายามเขียนไฟล์ vector layer ไปยังไฟล์ csv ด้วยรหัสนี้: QgsVectorFileWriter.writeAsVectorFormat(mylayer, r'c:\temp\xyz.csv', "utf-8", None, "CSV") มันสามารถส่งออกไปยัง xyz.csv แต่แสดงเฉพาะคุณลักษณะใน csv ไม่ใช่คอลัมน์รูปทรงเรขาคณิต ฉันจะส่งออกทั้งแอตทริบิวต์และข้อมูลเชิงพื้นที่ไปยังไฟล์ csv ได้อย่างไร

12 pyqgis  fields-attributes  geometry  csv 

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

12 arcgis-desktop  intersection  clip 

ฉันพยายามสร้างแผนที่โดยใช้ Javascript API ฉันต้องการใช้แผนที่ฐานของตัวเอง แต่ฉันไม่สามารถหารหัสที่ง่ายที่สุดในการทำเช่นนี้ ฉันมักจะเห็นสิ่งนี้: function init() { map = new esri.Map("mapDiv", { basemap: "satellite", center: [-97.395, 37.537], zoom: 11 }); แต่ฉันไม่ต้องการใช้ ArcGIS basemaps ออนไลน์ รหัสที่จะทำให้บริการแผนที่ฉันเป็นเจ้าภาพ basemap ของฉันคืออะไร

12 arcgis-server  arcgis-javascript-api  arcgis-online-basemaps 

ฉันมีสคริปต์ arcpy ง่าย ๆ ในการอัปเดตเขตข้อมูลในรูปแบบไฟล์จุดพร้อมข้อมูลจากคุณสมบัติรูปหลายเหลี่ยมที่อยู่ภายใน ใช้เวลา 9 นาทีในการทำ 100 คะแนนใน arcpy แต่การเข้าร่วมเชิงพื้นที่ใน arcmap นั้นเกิดขึ้นทันที ฉันแน่ใจว่ามีวิธีการแก้ไขปัญหานี้อย่างรวดเร็ว ใครบางคนชี้ให้ฉันในทิศทางที่ถูกต้อง? import took 0:00:07.085000 extent took 0:00:05.991000 one pt loop took 0:00:03.780000 one pt loop took 0:00:03.850000 one pt loop took 0:00:03.791000 import datetime t1 = datetime.datetime.now() import arcpy t2 = datetime.datetime.now() print "import took …

12 arcpy  performance 

ทุกเดือนฉันต้อง QA / QC ข้อมูลถนนสำหรับเขตทั้งหมด (มากกว่า 47,000 บันทึก) ขั้นตอนแรกคือการเปรียบเทียบสองไฟล์ ปัจจุบันฉันใช้เครื่องมือที่สร้างขึ้นเองผ่าน VBA เมื่อประมาณ 5 ปีที่แล้ว มันทำงานเปรียบเทียบสองไฟล์และรายงานความแตกต่างกลับมาตามปัจจัยสองสามประการ (เพิ่มบันทึก, ลบเรคคอร์ด, เปลี่ยนที่อยู่, เปลี่ยนเมือง, เปลี่ยนความยาวเซ็กเมนต์ / ยอดจุดเปลี่ยน) จากนั้นฉันจะส่งออกสิ่งเหล่านี้ไปยังไฟล์รูปร่างส่วนบุคคล นอกจากนี้เครื่องมือทำงานได้ArcGIS 9.3แต่จะไม่ทำงานอีกต่อไปArcGIS 10.1 เครื่องมือใช้งานได้ดี แต่ใช้เวลาประมาณหนึ่งชั่วโมงจึงจะเสร็จสมบูรณ์ ฉันเริ่มทำงานเพื่อหวังว่าจะได้พัฒนาเครื่องมือนี้เพื่อเพิ่มประสิทธิภาพ ผมมองเป็นความหลากหลายของวิธีการรวมทั้งและCompare Feature ToolSelect by Location เครื่องมือทั้งสองนี้ทำงานได้เร็วขึ้นอย่างมาก แต่ดูเหมือนว่าพวกเขาจะไม่จับทุกสิ่งที่ฉันต้องการ ปัญหาใหญ่ที่ฉันมีCompare Feature Toolคือฉันไม่สามารถเลือกฟิลด์ที่จะเปรียบเทียบ มันเริ่มต้นที่ ObjectID เมื่อฉันต้องการเปรียบเทียบโดย Segment_ID ไม่มีตัวเลือกให้เลือก ใครบ้างมีทางออกหรือแนวคิดที่ดีเกี่ยวกับวิธีการเปรียบเทียบแบบเต็มเช่นที่ฉันอธิบาย ด้วยเหตุผลบางอย่าง ArcGIS ไม่รู้จักการเข้าร่วมที่ฉันพยายามทำตามคำตอบของ RyanDalton แต่ …

12 arcgis-desktop  arcgis-10.1  modelbuilder  geoprocessing  data-quality