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

ฉันกำลังพยายามใช้ ArcGIS 10.2 เพื่อสร้างบัฟเฟอร์ของจุดตามพื้นที่ที่กำหนดไว้ล่วงหน้า (เช่น 400 ตารางกิโลเมตร) ยิ่งไปกว่านั้นบัฟเฟอร์ของบางจุดอยู่ใกล้กับแนวชายฝั่งที่ต้องการให้บัฟเฟอร์ถูกตัดที่แนวชายฝั่งและยังคงมีพื้นที่เดียวกับที่อยู่ในบก (400 ตารางกิโลเมตร) ไม่มีใครรู้ว่าสิ่งนี้สามารถทำได้ด้วยตัวสร้างแบบจำลองหรือ Arcpy? ฉันมีทักษะ จำกัด กับ Arcpy และ R แต่ยินดีที่จะใช้สคริปต์บางอย่างเพื่อแก้ปัญหานี้ โปรดดูภาพด้านล่างแสดงภาพกราฟิกของสิ่งที่ฉันพยายามที่จะบรรลุ

10 arcpy  arcgis-10.2  r  buffer 

ฉันได้นำเข้า rasters บางอย่างที่ควรจะเชื่อมต่ออย่างต่อเนื่อง แต่พวกเขาสร้างตารางที่มีขอบกั้นระหว่างพวกเขา ฉันอาจสันนิษฐานว่านี่เป็นเรื่องของการสอบเทียบระหว่างแรสเตอร์แต่ละอัน แต่ฉันไม่แน่ใจว่าจะรันมันอย่างไรในรูปแบบดาต้า มีใครเคยมีประสบการณ์เช่นนี้มาก่อนและอาจจะรู้วิธีแก้ปัญหาที่เป็นไปได้เพื่อให้ภาพรวมของ rasters ทั้งหมดจะราบรื่นขึ้น ทั้งหมดนี้ควรจะเชื่อมโยงกันในเชิงภูมิศาสตร์ ดังนั้นปัญหาของ 'ขอบเขต' ของสแควร์ ฉันใช้ QGIS บน windows

10 qgis  raster  windows  smoothing 

ฉันกำลังมองหาการสร้างเครื่องมือที่ดูเหมือนง่ายในWebAppBuilder Map ของฉัน พื้นฐานบางอย่าง: เนื้อหาถูกสร้างขึ้นใน ArcGIS 10.3.1 เนื้อหาถูกเผยแพร่ไปยังเซิร์ฟเวอร์ ArcGIS 10.3.1 เป็นบริการแคช แผนที่ได้รับการพัฒนาบน ArcGIS Online แผนที่ถูกส่งไปยัง WebAppBuilder วิดเจ็ตแบบใช้งานนอกกรอบใช้งานอยู่ในปัจจุบัน (มีWebAppBuilder เวอร์ชันนักพัฒนาซอฟต์แวร์ดังนั้นฉันสามารถทำงานกับโค้ดเพื่อสร้างแบบกำหนดเองได้หากต้องการ) ภารกิจ: เพื่อทำซ้ำฟังก์ชั่นของSelect Features แบบโต้ตอบใน WebAppBuilder ถึงวันที่ฉันได้ลองใช้ชุดรูปแบบของQuery Widgetซึ่งช่วยให้ผู้ใช้สามารถเลือกคุณลักษณะแบบโต้ตอบโดยใช้ "ใช้ตัวกรองเชิงพื้นที่เพื่อ จำกัด คุณสมบัติ" จากนั้นเลือก "เฉพาะคุณลักษณะที่สัมผัสพื้นที่ที่กำหนดโดยผู้ใช้" (ดูภาพด้านล่าง) อย่างไรก็ตามสิ่งนี้ไม่ได้พิสูจน์ว่าประสบความสำเร็จเนื่องจากเป็นสิ่งที่ฉันต้องการ ฉันแค่ต้องการให้ผู้ใช้สามารถเลือกรูปทรงเรขาคณิต (รูปหลายเหลี่ยมเส้นจุด) และให้พวกเขาเน้นบนแผนที่ เป็นโบนัส - ฉันต้องการให้ผู้ใช้สามารถควบคุมสีที่ไฮไลต์สิ่งที่เลือก - คล้ายกับตัวเลือกการเลือก (ดูภาพด้านล่าง)

10 arcgis-desktop  arcgis-server  arcgis-10.3  arcgis-online  arcgis-web-appbuilder 

มีวิธีการกำจัด "slivers" ขนาดเล็กระหว่างรูปหลายเหลี่ยมที่ใช้Rหรือไม่? วิธีการแก้ปัญหาที่ต้องการจะสร้างใหม่SpatialPolygonsDataFrameที่ขอบเขตร่วมกันระหว่างรูปหลายเหลี่ยมเป็นเรื่องบังเอิญ ฉันสนใจโซลูชันที่ใช้Rเป็นพิเศษมากกว่า ArcMap หรือ QGIS ฉันก็อยากจะได้ยินคำอธิบายว่าทำไมช่องว่างเหล่านี้อยู่ในตอนแรก นี่คือตัวอย่างของข้อมูลเชิงพื้นที่ที่ฉันทำงานด้วย: library(rgdal) library(sp) library(tigris) library(magrittr) library(leaflet) library(gplots) # This project will use WGS 84 projected coordinate system crs_proj <- CRS("+init=epsg:4326") # These are the FIPS codes of the specific block groups in my study area sel <- c("530330079005", "530330079001", "530330079004", "530330085002", …

10 polygon  leaflet  r  tiger  slivers 

ฉันกำลังวิเคราะห์ชุดข้อมูลจุด postGIS ใน QGIS เวอร์ชัน 2.8.2 ประมาณ 200k ฟีเจอร์โดยรวม ในการคำนวณค่าเฉลี่ยถ่วงน้ำหนักฉันใช้เทคนิคการแก้ไขค่าน้ำหนักระยะทางผกผันของ Inverse อย่างไรก็ตามฉันกังวลว่าค่า "สูง" ที่ฉันเห็นในผลลัพธ์นั้นเป็นผลมาจากค่าที่เกินขอบเขตจำนวนหนึ่งไม่ใช่พื้นฐานของรูปแบบเชิงพื้นที่ที่มีความหมายใด ๆ มีวิธีในปลั๊กอิน SAGA IDW - หรือเครื่องมือ IDW อื่น - เพื่อกำหนดจำนวนขั้นต่ำที่จำเป็นสำหรับพื้นผิวแรสเตอร์เพื่อคำนวณค่าเฉลี่ยถ่วงน้ำหนักวิธีเดียวกับที่จะพิจารณาจำนวนคะแนนสูงสุดหรือไม่

10 qgis  saga  inverse-distance-weighted 

เป็นไปได้หรือไม่ที่จะดึงขนาดของ shapefile โดยใช้ python และ arcpy? ถ้าเป็นเช่นนั้นได้อย่างไร

10 arcpy  shapefile  file-size 

ฉันพยายามที่จะตรวจสอบว่าการปรากฏตัวของกองกำลังทหารประจำการจำนวนมากในพื้นที่นั้นมีความสัมพันธ์เชิงพื้นที่กับอาชญากรรมที่รุนแรงในระดับที่สูงขึ้น / ต่ำลงหรือไม่ นั่นคือพื้นที่โดยรอบฐานทัพขนาดใหญ่มีความรุนแรงโดยเฉลี่ยมากกว่าพื้นที่ที่ไม่ได้อยู่ใกล้ฐานทัพหรือไม่ ฉันกำลังทำงานกับชุดข้อมูลสองชุดต่อไปนี้: (1)ชุดข้อมูลจุดของฐานทัพทหารในสหรัฐอเมริกาและระดับกองกำลังที่เกี่ยวข้อง: (2)ชุดข้อมูลทั่วประเทศเกี่ยวกับอัตราการเกิดอาชญากรรมรุนแรงตามเมือง / เมือง: ฉันรู้สึกว่าฉันกำลังมองหาโมเดลแรงโน้มถ่วงตามแบบที่ฟังก์ชั่น "มวล" ให้ระดับกองทหารในแต่ละฐาน ดังนั้นการปรากฏตัวของกองทหารขนาดใหญ่จะมีอิทธิพลเหนือพื้นที่ขนาดใหญ่และจะมีผลที่แข็งแกร่งกว่าใกล้กับจุดศูนย์กลางของมวล (เช่นที่ตั้งจุดในชั้น GIS) ฉันคิดว่าในความคิดมันจะมีลักษณะเช่นนี้: ในแผนภาพนี้X, Y, Zเป็นตัวแทนฐานทัพทหาร a, b, c, dแต่ละเมืองเป็นตัวแทน (แต่ละแห่งมีเขตข้อมูลอัตราความรุนแรงในตารางคุณลักษณะ) การไล่ระดับสีรอบฐานแสดงถึงเขตอิทธิพลซึ่งจะลดลงแบบเอกซ์โปเนนเชียลเมื่ออยู่ห่างจากศูนย์กลางเซนทรอยด์ การปรากฏตัวของกองทหารที่ใหญ่กว่านั้นเท่ากับรัศมีอิทธิพลที่ใหญ่กว่า (ด้วยระยะทางสูงสุดที่กำหนด) และยังมีอิทธิพลที่มากขึ้นใกล้กับศูนย์กลางที่สัมพันธ์กับพื้นที่ใกล้เคียงกับฐานที่เล็กกว่า แต่ละเมืองจะได้รับคะแนนตามการรวมขนาดของเวกเตอร์ "แรง" ทั้งหมดจากฐานโดยรอบทั้งหมดที่มีอิทธิพลต่อรัศมีที่พวกเขาอยู่เช่นในแผนภาพของฉันเมืองaจะมีคะแนน 0 เพราะอยู่ ด้านนอกของรัศมีฐานใด ๆ City bจะได้รับอิทธิพลจากฐานXเท่านั้น เมืองcจะได้รับอิทธิพลจากฐานZเท่านั้นและคะแนนจะต่ำกว่าbเนื่องจากXเป็นฐานที่ใหญ่กว่าZมาก ในที่สุด City dตั้งอยู่ภายในรัศมีของทั้งฐานXและYมันจะได้รับคะแนนจากการสรุปขนาดของอิทธิพลจากฐานทั้งสอง ฉันจะดูว่ามีความสัมพันธ์ระหว่างคะแนนที่สูงขึ้นสำหรับเมืองหนึ่งกับอัตราความรุนแรงที่สูงขึ้นหรือไม่ ฉันได้ดูโมเดลแรงโน้มถ่วงหลายแบบ ( Huff Modelsและอื่น ๆ ) แต่ไม่พบเท่าที่ QGIS …

10 qgis  spatial-statistics  statistics  correlation 

งาน ArcGIS Pro สนับสนุนเครื่องมือ Python Toolbox ซึ่งทำงานกับ Map frames หรือไม่? เหตุผลที่ฉันถามคือฉันติดเมื่อพยายามทำต่อไปนี้: เริ่ม ArcGIS Pro 1.1.1 สร้างโครงการใหม่ - ฉันเรียกว่า TestProject ของฉันและวางใน C: \ Temp ใช้บานหน้าต่างโครงการเพื่อเพิ่มการเชื่อมต่อโฟลเดอร์ไปยังตำแหน่งที่ฉันมีไฟล์รูปร่างของประเทศต่างๆในโลกจาก Natural Earth () ลากและวาง ne_10m_admin_0_countries.shp ลงในแผนที่เพื่อสร้างเลเยอร์ที่เรียกว่า ne_10m_admin_0_countries แทรกเค้าโครงใหม่ - ฉันใช้แนวนอน A3 แทรกกรอบแผนที่ใหม่ลงในเค้าโครง ในบานหน้าต่างโครงการสร้างกล่องเครื่องมือ Python ใหม่ในโฟลเดอร์ TestProject - ฉันเรียกว่า TestPYT ของฉัน คลิกขวาที่ TestPYT ในบานหน้าต่างโครงการเพื่อแก้ไข แทนที่รหัสด้วยด้านล่างเพื่อให้เครื่องมือ Python สองเครื่องมือที่เรียกว่าชิลีและสวิตเซอร์แลนด์ …

10 arcpy  arcgis-pro  python-toolbox  tasks 

ฉันต้องการทราบว่าเครื่องมือเรขาคณิตของ QGIS -> polygon centroid คำนวณพิกัดจุดอย่างไร ยกตัวอย่างเช่นมันแบ่งนักเรียนประจำของรูปหลายเหลี่ยมออกเป็นจุดเล็ก ๆ หรือไม่ใช้สองคู่ของสองจุด (ความยาวและความกว้าง) ด้วยระยะทางสูงสุดและข้ามเส้นทั้งสองนี้เป็นเซนทรอยด์?

10 qgis  polygon  geometry  centroids 

ฉันทำ semivariogram ใน R โดยใช้ package gstat, variogram()function ฉันต้องการที่จะดูว่ามี autocorrelation เชิงพื้นที่ในแบบจำลองของฉัน (ความอุดมสมบูรณ์ของสปีชีส์เป็นหน้าที่ของที่อยู่อาศัยในพื้นที่ต่าง ๆ ห่างกันเพียงไม่กี่กม. ถึง 900 กม. โดยใช้ glmm) หน่วยของฉันอยู่ในกม. และดังนั้นการตีความของฉันก็คือช่วงนั้นอยู่ที่มากกว่า 100 กม. จนกว่าการเปลี่ยนแปลงเชิงพื้นที่สัมพันธ์จะไม่เป็นปัญหาอีกต่อไป ฉันสงสัยว่าบางคนสามารถอธิบายได้ว่าทำไมนักเก็ตดูจึงสูง นี่หมายความว่าแม้ในสถานที่ใกล้เคียงกัน แต่ก็มีความแตกต่างค่อนข้างสูงใช่ไหม? หรือรูปคลื่นที่แปรปรวนแบบนี้หมายความว่าฉันควรจะปรับจำนวนความล่าช้าและระยะทางที่ช้ากว่าของฉันจนกว่าฉันจะได้รูปร่างที่มากกว่าปกติหรือไม่? เพื่อที่จะตรวจสอบเพิ่มเติมอีกเล็กน้อยฉันยังใช้ฟังก์ชั่นvariog()ในแพคเกจ geoR และใช้breaks=seq(0,100,10)เพื่อลองดูที่ระยะทางที่ใกล้กว่า (ใช้จุดเดียวกันและส่วนที่เหลือแบบเดียวกัน) สิ่งนี้บ่งชี้ว่าจุดที่ใกล้ที่สุดนั้นแตกต่างกันมากขึ้นซึ่งก็ไม่สมเหตุสมผล บางทีนี่อาจบ่งบอกว่าไม่มีความสัมพันธ์เชิงพื้นที่อัตโนมัติและแบบจำลองของฉันได้ทำสิ่งนี้แล้ว ฉันพบแหล่งข้อมูลที่ยอดเยี่ยมนี้"Geostats ที่ไม่มีน้ำตา"และในหน้า 51 มีคำแนะนำที่ดีเกี่ยวกับรูปแปรที่เหมาะสม ตามคำแนะนำนี้คนแรกของฉันดูเหมือนจะมีช่วงที่ถูกต้อง ดังนั้นสิ่งนี้กลับไปที่คำถามแรก - ฉันจะตีความสิ่งนี้ได้อย่างไร

10 r  interpolation  spatial-statistics  autocorrelation  variogram 

ฉันกำลังสร้างแผนที่ด้วยArcGIS Desktop v10.2.2โดยแสดงพิพิธภัณฑ์เป็นจุด เลเยอร์พิพิธภัณฑ์ไม่มีสัญลักษณ์และป้ายกำกับทำจากวงกลมที่เต็มไปด้วยสีที่มี ID ภายในซึ่งนำมาจาก[ID]สนาม ฉันต้องการเพิ่มรายชื่อพิพิธภัณฑ์ทั้งหมดในตำนานแผนที่เลย์เอาต์มีรายการสัญลักษณ์ของฉันพร้อมรหัสประจำตัวและคำอธิบาย (ชื่อของพวกเขา) ที่นำมาจาก[LABEL]สนามโดยไม่ต้องสร้างมันขึ้นมาเอง ฉันรู้ว่าโซลูชันแบบแมนนวลที่เป็นไปได้คือการตั้งสัญลักษณ์ของเลเยอร์เป็น "ค่าที่ไม่ซ้ำ" จากนั้นเพิ่มค่าทั้งหมด (มากกว่า 100!) และตั้งค่าหนึ่งครั้งสำหรับแต่ละสัญลักษณ์สัญลักษณ์ป้ายกำกับที่เหมือนกัน แต่ฉันหวังว่าจะมีทางออกที่ดีกว่า! ลองดูแผนที่ของฉันด้านล่าง เลเยอร์คือ "Musei" (พิพิธภัณฑ์) และอย่างที่คุณเห็นในตำนานทางด้านซ้ายไม่มีสัญลักษณ์ ด้านขวาคุณจะเห็นสัญลักษณ์ตำนานอธิบายวงกลมสีเทาเข้มพร้อมข้อความสีแดงด้านในถ่ายจาก[ID]สนาม ฉันต้องการเพิ่มคอลัมน์ที่มีวงกลมทั้งหมดที่มี ID และป้ายกำกับที่เกี่ยวข้องในเค้าโครงของฉัน - อัปเดต - นี่คือการตั้งค่า Field Calculator ที่ฉันใช้หลังจากทำตามวิธีแก้ปัญหาที่เสนอโดย @FelixIP แต่คะแนนไม่ได้อยู่ในแนวตั้งตามที่ฉันคาดหวัง แต่พวกมันขยับไปมาตรงนี้ ฉันลองใช้ค่าขั้นตอนที่ต่างกัน: 10, 100, 1,000, 10,000

10 arcgis-desktop  arcmap  arcgis-10.2  symbology  legend 

ฉันต้องตรวจสอบการสังเกตการณ์ของนกในระยะเวลานานขึ้นเพื่อดูรายการที่ซ้ำกัน / ทับซ้อนกัน ผู้สังเกตการณ์จากจุดต่าง ๆ (A, B, C) ทำการสังเกตการณ์และทำเครื่องหมายไว้บนแผนที่กระดาษ เส้นเหล่านั้นซึ่งนำมาสู่คุณสมบัติเส้นที่มีข้อมูลเพิ่มเติมสำหรับสปีชีส์จุดสังเกตและช่วงเวลาที่เห็น โดยปกติผู้สังเกตการณ์จะสื่อสารกันทางโทรศัพท์ในขณะที่สังเกต แต่บางครั้งพวกเขาก็ลืมดังนั้นฉันจึงได้เส้นที่ซ้ำกัน ฉันลดข้อมูลลงไปเป็นเส้นเหล่านั้นซึ่งแตะวงกลมดังนั้นฉันไม่จำเป็นต้องทำการวิเคราะห์เชิงพื้นที่ แต่เพียงเปรียบเทียบช่วงเวลาสำหรับแต่ละสปีชีส์และมั่นใจได้ว่าเป็นบุคคลเดียวกันที่พบโดยการเปรียบเทียบ . ตอนนี้ฉันกำลังมองหาวิธีใน R เพื่อระบุรายการเหล่านั้นซึ่ง: จะทำในวันเดียวกันกับช่วงเวลาที่ทับซ้อนกัน และมันอยู่ที่ไหนสายพันธุ์เดียวกัน และซึ่งทำจากจุดสังเกตที่แตกต่างกัน (A หรือ B หรือ C หรือ ... )) ในตัวอย่างนี้ฉันพบว่ามีรายการที่ซ้ำกันของบุคคลเดียวกัน จุดสังเกตแตกต่างกัน (A <-> B) สปีชีส์เหมือนกัน (Sst) และช่วงเวลาของการเริ่มต้นและสิ้นสุดของเวลาทับซ้อนกัน ตอนนี้ฉันจะสร้างเขตข้อมูลใหม่ "ทำซ้ำ" ใน data.frame ของฉันโดยให้ทั้งสองแถวมีรหัสทั่วไปเพื่อที่จะสามารถส่งออกและตัดสินใจในภายหลังว่าจะทำอย่างไร ฉันค้นหาวิธีแก้ปัญหาที่มีอยู่มากมาย แต่ไม่พบว่าเกี่ยวข้องกับความจริงที่ว่าฉันต้องเซตย่อยกระบวนการสำหรับสปีชีส์ (โดยเฉพาะอย่างยิ่งที่ไม่มีลูป) และต้องเปรียบเทียบแถวสำหรับจุดสังเกต 2 + x ข้อมูลบางส่วนที่จะเล่นด้วย: …

10 r  time  comparison  duplication 

ฉันกำลังพยายามจัดแนวแรสเตอร์สองกริดในอาร์เมื่อจัดตำแหน่งแล้วฉันต้องการที่จะเพิ่มพวกเขาเข้าด้วยกัน ฉันได้ลองตรวจสอบว่าการทำแบบstackนี้จะได้ผลหรือไม่: grid_snap <- stack(habi_sdw, Pop_sdw) และฉันได้รับข้อผิดพลาดต่อไปนี้: ข้อผิดพลาดใน comparRaster (x): ขอบเขตที่ต่างกัน ตะแกรงแรสเตอร์มีคุณสมบัติดังต่อไปนี้: show(habi_sdw) # class : RasterLayer # dimensions : 9187, 9717, 89270079 (nrow, ncol, ncell) # resolution : 0.00892857, 0.00892857 (x, y) # extent : -28.83706, 57.92186, -36.02464, 46.00214 (xmin, xmax, ymin, ymax) # coord. ref. : +proj=longlat +ellps=WGS84 …

10 raster  r  clip  snapping  resolution 

ฉันต้องการเรียกใช้สคริปต์ python ผ่านqgis --code myscript.pyแล้วออกจากทันที ฉันใช้iface.actionExit().trigger()ซึ่งฆ่า QGIS เมื่อฉันเรียกใช้จากคอนโซลหลาม --codeแต่ไม่เมื่อใส่ในสคริปต์ที่ผมผ่านไป วิธีที่ถูกต้องในการออกจากทันทีคืออะไร ฉันใช้ QGIS 2.0.1 อัปเดต : ฉันได้ลองsys.exit()แล้ว QGIS จับมันและดึงหน้าต่างขึ้นมาด้วยสิ่งต่อไปนี้: An error occured during execution of following code: execfile('myscript.py') Traceback (most recent call last): File "", line 1, in File "myscript.py", line 14, in sys.exit() SystemExit อัปเดต : ใช้os.kill(os.getpid(), 9)งานได้ แต่มันแฮ็คสกปรกและฉันกำลังมองหาสิ่งที่ดีกว่า

10 qgis  python  pyqgis 

ผมนำเข้าหลายชุดข้อมูลภูมิศาสตร์เข้า PostGIS และพวกเขามีความแตกต่างกันSRIDของ (บางคนมีEPSG:3857บางEPSG:4326อย่างบางอย่างอื่น) ฉันต้องการสร้างเพิ่มเติมgeometry columnเช่น the_geom_mercatorด้วยSRID EPSG:3857และยังคงgeomคอลัมน์เดิมในสิ่งที่SRIDมันเข้ามา ฉันจะใช้ฟังก์ชัน PostGIS ได้อย่างไร

10 postgis  coordinate-system  srid  mercator