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

ฉันต้องการใช้การฉายภาพใน QGIS ที่ฉันไม่เห็นการสนับสนุนโดยตรง ตัวอย่างจะเป็นอะไรที่คล้ายกับการฉายคิวบิกก้อนนี้ http://www.csiss.org/map-projections/Polyhedral_Globes/Cube_Gnomonic.pdf คำแนะนำในการค้นหารายละเอียดและใช้การคาดการณ์นี้หรือสิ่งที่คล้ายกัน?

9 qgis  cartography 

ฉันต้องการคลิปชุด polylines (เส้นสีดำในภาพด้านล่าง) ไปยังขอบเขตด้านนอกของรูปหลายเหลี่ยม ควรละเว้นช่องว่างใด ๆ ภายในรูปหลายเหลี่ยม ผลลัพธ์ที่ดีที่สุดของฉันคือเส้นสีเหลืองประ บรรทัดเริ่มต้นอาจจะหรืออาจจะไม่ตรง ภาพเป็นตัวอย่างที่เรียบง่าย แต่ในความเป็นจริงรูปหลายเหลี่ยมนั้นซับซ้อนกว่ามากและมีหลายร้อยบรรทัด ฉันไม่คิดว่าตัวเรือนูนจะทำงานได้ (แต่ฉันอาจผิด) ฉันเปิดให้แก้ปัญหาใน arcgis, qgis, arcpy, หุ่นดี ฯลฯ การเขียนโปรแกรมจะอยู่ในหลามโดยฉันเปิดให้ตัวเลือกอื่น ๆ หากจำเป็น Arcgis น่าจะดีกว่าที่จะทำให้เพื่อนร่วมงานของฉันแบ่งปันเครื่องมือได้ง่ายขึ้น แต่ไม่ใช่ข้อกำหนด สิ่งที่ดีที่สุดที่ฉันคิดได้ในตอนนี้คือการตัดเส้นแต่ละเส้นด้วยรูปหลายเหลี่ยมที่สร้างชุดของจุดที่จุดตัดของขอบเขตทั้งหมด เรียงคะแนนตามระยะทางถึงจุดเริ่มต้นของบรรทัด จุดที่อยู่ไกลและใกล้ที่สุด (FAC) จะเป็นขอบเขตด้านนอกของรูปหลายเหลี่ยม จากนั้นใช้จุด FAC เพื่อเลือกจุดยอดที่เหมาะสมจากบรรทัดเดิมและสร้างเส้นประสีเหลืองจากจุดที่เหมาะสม ควรใช้งานได้ แต่ดูเหมือนจะซับซ้อนเกินความจำเป็น ความคิดเพิ่มเติมบางประการ: เส้นดังกล่าวเป็นเส้นตรง "เพียงพอ" ที่การคำนวณระยะทางง่าย ๆ ระหว่างจุดควรใช้งานไม่จำเป็นต้องมีการอ้างอิงเชิงเส้น นี่จะเป็นเรื่องง่ายใน arcpy ถ้ามีเครื่องมือในการแยกเส้นตรงจุด แต่ฉันไม่สามารถหาได้ คิดว่าทุกคน?

9 qgis  arcgis-desktop  arcpy  shapely 

ฉันกำลังพยายามแมป CSV ที่มีจุด lat / lon กว่าล้านจุดในสหรัฐอเมริกาโดยไม่มีข้อมูลเมตา ฉันไม่เคยพยายามแมปข้อมูลจำนวนมากนี้และต้องการทราบขั้นตอนที่ฉันควรทำ ฉันรู้ว่ามันเป็นข้อมูลที่มากเกินไปสำหรับ QGIS ดังนั้นฉันจึงพยายามแยกมันออกเป็น 100,000 บรรทัด CSV แต่แล้วก็มีไฟล์จำนวนมากที่ฉันไม่รู้ว่าจะทำอย่างไร ฉันคิดว่ามันจะเป็นการดีที่สุดที่จะรวบรวมคะแนนจากการสำรวจสำมะโนประชากร แต่ไม่ทราบวิธีที่ดีที่สุดในการใช้ Javascript (ฉันไม่เคยใช้ Python) หรือเครื่องมือบรรทัดคำสั่ง? ข้อเสนอแนะใด ๆ

9 qgis  turf 

ฉันต้องการค้นหาอัลกอริทึมหรือวิธีการที่สามารถตรวจจับ latitude longitude จุดที่ผิดปกติในวิถีในระหว่างการโพสต์การประมวลผลซึ่งสามารถแก้ไขได้ (นำกลับไปสู่เส้นทางของวิถีตามเพื่อนบ้าน) เป็นตัวอย่างของคะแนนผิดปกติที่ฉันต้องการตรวจจับและแก้ไขฉันได้แนบรูปภาพที่แสดง: ฉันได้ลองใช้ตัวกรองคาลมานเพื่อรวบรวมข้อมูลให้ดีที่สุดเท่าที่จะเป็นไปได้ แต่ดูเหมือนว่ามันจะไม่ทำงานอย่างมีประสิทธิภาพเพียงพอสำหรับค่าผิดปกติที่รุนแรงมากขึ้น (ข้อมูลดิบเป็นสีน้ำเงินข้อมูลเรียบเป็นสีแดง): UKF ของฉันอาจไม่ได้รับการปรับเทียบอย่างถูกต้อง (แต่ฉันค่อนข้างแน่ใจว่าเป็น) วิถีการเคลื่อนที่นั้นเป็นของวอล์กเกอร์นักวิ่งนักปั่นจักรยาน - การเคลื่อนไหวที่ขับเคลื่อนด้วยมนุษย์ซึ่งสามารถเริ่มต้นและหยุดได้ แต่ไม่สามารถเปลี่ยนความเร็วหรือตำแหน่งอย่างรวดเร็วหรือฉับพลัน โซลูชันที่ไม่พึ่งพาข้อมูลเวลา (และข้อมูลตำแหน่งเท่านั้น) จะมีประโยชน์อย่างมาก (เนื่องจากข้อมูลที่กำลังประมวลผลอาจไม่ประกอบด้วยข้อมูลเวลา) อย่างไรก็ตามฉันทราบดีว่ามีวิธีการแก้ปัญหาแบบนี้ไม่น่าเป็นไปได้ดังนั้นฉันมีความสุขเท่ากันที่จะมีทางออกใด ๆ ! โดยหลักการแล้ววิธีการแก้ปัญหาจะตรวจจับค่าผิดปกติเพื่อให้สามารถแก้ไขได้ทำให้เกิดวิถีการแก้ไขที่ถูกต้อง: ทรัพยากรที่ฉันกรองผ่าน: Smooth GPS data- /programming/1134579/smooth-gps-data Common GPS and Geospatial Tracking Challenges and Solutions- http://www.toptal.com/gis/adventures-in-gps-track-analytics-a-geospatial-primer (ดูเหมือนว่าโซลูชันจะสูญเสียความแม่นยำของข้อมูล) ฉันควรใช้อัลกอริทึมใดเพื่อลบค่าผิดปกติในข้อมูลติดตาม

9 gps  latitude-longitude  noise  outlier 

เป็นไปได้ไหมที่จะกำหนดแป้นพิมพ์ลัดให้กับสคริปต์ (กำหนดเองหรือไม่ ... ) ของเครื่องมือประมวลผลใน QGIS? ฉันไม่พบข้อมูลใด ๆ เกี่ยวกับมัน

9 pyqgis  qgis-processing 

มีเอกสารอย่างเป็นทางการสำหรับโมดูล python GDAL หรือไม่ มีตำราหลายเล่ม แต่มีข้อมูลไม่เพียงพอในแต่ละฟังก์ชันหรืออาร์กิวเมนต์สำหรับฟังก์ชันพื้นฐานที่สุด

9 python  gdal  documentation 

สิ่งที่ฉันพยายามทำ: วนลูปผ่านไฟล์รูปร่างจุดและเลือกแต่ละจุดที่อยู่ในรูปหลายเหลี่ยม โค้ดต่อไปนี้ได้รับแรงบันดาลใจจากตัวอย่างการสืบค้นเชิงพื้นที่ที่ฉันพบในหนังสือ: mitte_path = r"D:\PythonTesting\SelectByLocation\mitte.shp" punkte_path = r"D:\PythonTesting\SelectByLocation\punkte.shp" polygon = QgsVectorLayer(mitte_path, 'Mitte', 'ogr') points = QgsVectorLayer(punkte_path, 'Berlin Punkte', 'ogr') QgsMapLayerRegistry.instance().addMapLayer(polygon) QgsMapLayerRegistry.instance().addMapLayer(points) polyFeatures = polygon.getFeatures() pointsCount = 0 for poly_feat in polyFeatures: polyGeom = poly_feat.geometry() pointFeatures = points.getFeatures(QgsFeatureRequest().setFilterRect(polyGeom.boundingBox())) for point_feat in pointFeatures: points.select(point_feat.id()) pointsCount += 1 print 'Total:',pointsCount มันใช้งานได้และมันเลือกชุดข้อมูล แต่ปัญหาก็คือมันเลือกตามขอบกล่องดังนั้นจึงทำให้เห็นได้ชัดว่าฉันไม่สนใจจุดกลับมา: …

9 python  pyqgis  select-by-location  spatial-query  qgis-2.10 

หลังจากใช้ SRC ที่ไม่ปกติเพื่อวัตถุประสงค์ในการสาธิตฉันต้องการล้างรายการ "ระบบพิกัดที่ใช้ล่าสุดของฉัน" เพื่อประโยชน์ของ praticity เป็นไปได้ไหม? ฉันใช้ QGIS 2.10.1

9 qgis  coordinate-system  qgis-2.10 

ฉันต้องการสร้างชุดข้อมูลเวกเตอร์แบบจุด 10,000 จุด (หรือใหญ่กว่า) ภายในไบนารีแรสเตอร์ซึ่งจุดควรถูก จำกัด พื้นที่ที่ค่าแรสเตอร์เป็น 1 ฉันลองทำตามขั้นตอนต่อไปนี้ รูปหลายเหลี่ยมแรสเตอร์ QGIS: Vector -> เครื่องมือวิจัย -> คะแนนสุ่ม วิธีนี้ใช้งานได้ดีมากถึง 2,000 คะแนน แต่สิ่งที่เหนือกว่าเพียงแค่ทำให้ QGIS ผิดพลาด มีวิธีในการสร้างชุดข้อมูลเวกเตอร์ที่มีคุณสมบัติจุดจำนวนมากที่ จำกัด โดยไบนารีแรสเตอร์ (หรือรุ่นที่มีรูปหลายเหลี่ยมของมัน) หรือไม่? เครื่องมือต่อไปนี้อยู่ในการกำจัดของฉันจัดอันดับจากที่นิยมมากไปน้อย: QGIS, Python, R, ArcGIS นี่คือสิ่งที่ฉันจะทำเพียง 10 เท่าที่มีคุณสมบัติจุด

9 qgis  arcgis-desktop  python  r  point-in-polygon 

การทำงานกับ QGIS 2.4.0 Chugiak ฉันมีคำถามเกี่ยวกับการตัดภาพแรสเตอร์ ฉันมีภาพถ่ายทางอากาศที่ฉันต้องการคลิปโดยใช้ extents ของ shapefile จนถึงตอนนี้ฉันใช้เครื่องมือแรสเตอร์ / การสกัด / clipper เพื่อทำสิ่งนี้ ปัญหาคือขนาดไฟล์ของเอาท์พุทแรสเตอร์ของฉันมีขนาดใหญ่กว่าไฟล์ต้นฉบับเกือบสามเท่า (เช่น 725.283 และ 249.693 kb ตามลำดับ) ในขณะที่มันครอบคลุมพื้นผิวที่เล็กลง ทำไมไฟล์ของฉันถึงใหญ่กว่า และมีวิธีการ 'ตั้งค่า' ขนาดไฟล์ที่ส่งออกหรือไม่? ฟังก์ชั่นความช่วยเหลือของปัตตาเลี่ยนอ้างถึงหน้าต่อไปนี้: http://www.gdal.org/gdal_translate.htmlแม้ว่าฉันจะไม่ได้ประโยชน์อะไรมากมาย ฉันค่อนข้างใหม่กับ QGIS

9 qgis  raster  clip  file-size 

ฉันได้สร้าง data.frame ใน R ฉันต้องการรับ SpatialLinesDataFrame ดังนั้นฉันเข้าใจว่าฉันต้องแปลง data.frame เป็นบรรทัดบรรทัดลงใน SpatialLines และ SpatialLines เป็น SpatialLinesDataFrame ฉันทำอย่างนี้ filedata.frame=data.frame(matrix(file),ncol=14, byrow=T)) file=Lines(filedata.frame) ฉันได้รับ data.frame ที่ฉันต้องการ แต่ไม่ใช่บรรทัด ฉันแน่ใจว่าพื้นฐาน แต่ทำไมฉันไม่สามารถรับสายได้

9 r  line 

การเข้าถึงข้อมูลการค้นหาเคอร์เซอร์ อะไรคือความแตกต่างระหว่างโทเค็น"SHAPE@XY"และ"SHAPE@TRUECENTROID"อะไร พวกเขาคืนค่า coordnates เดียวกันบนรูปหลายเหลี่ยมทดสอบ: >>> cursor = arcpy.da.SearchCursor("test", "SHAPE@XY") >>> for row in cursor: ... print row[0] ... (559389.3838043335, 4239093.201390337) >>> cursor = arcpy.da.SearchCursor("test", "SHAPE@TRUECENTROID") >>> for row in cursor: ... print row[0] ... (559389.3838043335, 4239093.201390337) มันทำงานในลักษณะเดียวกันสำหรับคุณลักษณะรูปหลายเหลี่ยมหลายส่วน มันเคยให้ผลลัพธ์ที่แตกต่างกันหรือไม่?

9 arcpy  centroids 

ฉันพยายามค้นหาพิกัดของจุดสนใจของฉัน (จุด X, เครื่องหมายสีน้ำเงิน) ฉันใช้อุปกรณ์ GPS ในรถของฉันเพื่อรวบรวมพิกัดตามตำแหน่งที่ฉันจอดยานพาหนะของฉันทุกครั้งที่ฉันไปที่จุด x ดังนั้นหลังจากทำแบบฝึกหัดนี้เป็นเวลา 16 วันฉันก็สามารถหาพิกัดได้ 16 ชุดกระจายไปทั่วจุดที่ฉันสนใจ หลังจากพล็อตพิกัดเหล่านี้บนแผนที่ฉันสังเกตสิ่งต่อไปนี้สองครั้งหรือสามครั้งจากสิบครั้งอุปกรณ์ GPS ของฉันให้พิกัดที่ผิดซึ่งกลายเป็นเงียบไกลจากจุด X นอกจากนี้เนื่องจากการจราจรบางครั้งฉัน ไม่สามารถจอดใกล้จุด x และด้วยเหตุนี้ในกรณีนี้เช่นกันพิกัดที่ได้อยู่ไกลจากจุด X ปัญหา: จากพิกัด 16 ชุดที่ได้รับฉันจะใช้กระบวนการใดในการ จำกัด พิกัดหนึ่งชุดซึ่งอยู่ใกล้กับจุดสนใจของฉัน (จุด X)

9 gps  spatial-statistics 

สคริปต์ของฉันตัดกันกับเส้นหลายเหลี่ยม เป็นกระบวนการที่ยาวนานเนื่องจากมีมากกว่า 3000 บรรทัดและมากกว่า 500,000 โพลิกอน ฉันดำเนินการจาก PyScripter: # Import import arcpy import time # Set envvironment arcpy.env.workspace = r"E:\DensityMaps\DensityMapsTest1.gdb" arcpy.env.overwriteOutput = True # Set timer from datetime import datetime startTime = datetime.now() # Set local variables inFeatures = [r"E:\DensityMaps\DensityMapsTest.gdb\Grid1km_Clip", "JanuaryLines2"] outFeatures = "JanuaryLinesIntersect" outType = "LINE" # Make lines …

9 arcpy  geoprocessing  performance 

มีวิธีใดที่จะละลาย shapefile แต่เก็บฟิลด์คุณสมบัติไว้ในเวลาเดียวกัน? นี่คือวิธีที่ฉันต้องการเก็บรักษาเขตข้อมูล: มีสองเขตข้อมูลที่ฉันต้องการอนุรักษ์ใน shapefile คอลัมน์แรกคือ Net_ID ที่มีชนิดยาวและคอลัมน์ที่สองคือธรณีวิทยาที่มีประเภทสตริง ฉันต้องการยุบ shapefile โดยใช้ฟิลด์ NET_ID นอกจากนี้ฉันต้องการเก็บฟิลด์ธรณีวิทยาไว้ดังรูปด้านล่าง

9 arcgis-desktop  shapefile  arcgis-10.2  dissolve