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

ฉันได้สร้างตารางเรขาคณิตเป็น postgis แล้ว ฉันต้องการโหลดรูปร่างหลายไฟล์ลงในตารางนั้นโดยใช้ ogr2ogr แต่ฉันไม่พบพารามิเตอร์ที่สามารถระบุตารางนั้น ogr2ogr -append -f "PostgreSQL" PG:"dbname=db and so on" shapefile.shp แค่นั้นแหละ. ฉันจะเพิ่มพารามิเตอร์เพื่อระบุตารางได้อย่างไร

13 postgis  ogr2ogr 

ฉันใช้gIntersectฟังก์ชันจากrgeosแพ็คเกจเพื่อกำหนดรูปหลายเหลี่ยมที่ตัดกันในสอง SpatialPolygonsDataFrames (spdfs) ผลลัพธ์เป็นเมทริกซ์เชิงตรรกะที่มี rownames ซึ่งประกอบด้วยรูปหลายเหลี่ยม ID จาก sfdf # 1 และ colnames ID รูปหลายเหลี่ยมจาก spdf # 2 ฉันต้องการที่จะใช้เมทริกซ์นี้ไปยังระบบย่อยทั้งสอง spdfs gIntersectionเดิมถนัดรูปหลายเหลี่ยมที่ไม่มีการทับซ้อนก่อนที่จะใช้ฟังก์ชั่น ฉันสามารถแปลงเมทริกซ์เชิงตรรกะได้อย่างง่ายดายเป็นสองเวกเตอร์ของรหัส (หนึ่งสำหรับแต่ละไฟล์ PDF) ที่แสดงรูปหลายเหลี่ยมเหล่านั้นอย่างน้อยทับซ้อนกัน คำถามของฉันคือ - ฉันสามารถใช้เวกเตอร์ของรูปหลายเหลี่ยม ID เพื่อเซตย่อย spdf ได้หรือไม่ ถ้าเป็นเช่นนั้นมันจะเซตค่าใน@dataสล็อตได้อย่างถูกต้องหรือไม่? ฉันรู้ว่าฉันสามารถย่อย spdfs โดยใช้ไวยากรณ์นี้ได้: spdf[spdf@data$myVariable[someVector], ] สิ่งนี้จะใช้ไม่ได้กับ ID เนื่องจากฝังอยู่ลึกเข้าไปในโครงสร้างรายการ spdf

13 r 

มีวิธีเลือกตัวคั่นฟิลด์เมื่อส่งออก CSV ที่มีรูปทรงเรขาคณิตเป็น wkt ใน QGIS หรือไม่ ฉันมีไฟล์รูปร่างที่มีรูโดนัทและฉันต้องการส่งออกเป็น CSV สำหรับคู่ค้าที่ฐานข้อมูลไม่รองรับข้อมูลเชิงพื้นที่ ฉันทำตามวิธีการของ Greg Krakowเพื่อเพิ่มเขตข้อมูล WKT ไปยังไฟล์ที่ส่งออกของฉัน (เมื่อส่งออกไปยัง csv เพิ่มGEOMETRY=AS_WKTในฟิลด์ "Layer" แต่ตัวคั่นเขตข้อมูลที่ส่งออกเป็นเครื่องหมายจุลภาค: ปัญหาคือว่าแต่ละพิกัดของรูปหลายเหลี่ยมของฉัน (และฉันคิดว่ามันจะเหมือนกันสำหรับบรรทัด) จะถูกคั่นด้วยเครื่องหมายจุลภาคด้วยดังนั้นฉันจึงไม่สามารถนำเข้าข้อมูลกลับเข้าไปใน QGIS ด้วย shapefile จุดมันทำงานเหมือนมีเสน่ห์ ไม่มีใครรู้ว่าฉันสามารถเลือกเขตข้อมูลของฉันเมื่อส่งออกเป็น CSV หรือไม่ ยกตัวอย่างเช่นเลือกเครื่องหมายอัฒภาค? ฉันค้นหาบนเว็บและเอกสารประกอบของ QGIS แต่ไม่พบอะไรเลย

13 qgis  export  csv  wkt 

วิธีเดียวที่ดีในการสร้างแวดวงใน ArcMap คือการใช้ส่วนโค้งจริง (เส้นโค้งผสม) แต่ฉันจำเป็นต้องแปลงวงกลมเหล่านั้นให้เป็นจุดสุดยอด (ทำให้หนาแน่น) เนื่องจากปัญหาความเข้ากันได้ ฉันต้องการข้อมูลสำหรับ geodatabase และ SDO_geometry (arcsde) และใช่มีวิธี ฉันใช้: UPDATE layer1 a SET arctype = 'compound' WHERE has_compound_curves(a.shape) <> 0; จากนั้นใช้ ET Geotools เพื่อเพิ่มความหนาของชั้นโดยชั้น แต่ควรมีวิธีที่ง่ายกว่าในการทำตรงไปตรงมาใน Oracle หรือ ArcMap และฉันแค่ต้องการเพิ่มความหนาแน่นของฟีเจอร์ที่มีส่วนโค้งที่แท้จริง ฉันไม่ต้องการจุดยอดที่ไม่จำเป็นบนเส้นตรงเป็นต้น ฉันจะระบุส่วนโค้งที่แท้จริงใน ArcMap ได้อย่างไรและฉันจะทำให้พวกมันหนาแน่นได้เร็วและง่ายที่สุดได้อย่างไร เช่นเดียวกับ Oracle SQL

13 arcpy  arcmap  enterprise-geodatabase  oracle  density 

การใช้ QGIS ถ้าฉันโหลดเลเยอร์รูปหลายเหลี่ยมและสร้างคอลัมน์ใหม่ฉันจะเติมคอลัมน์นั้นด้วยวิธี: พื้นที่ของรูปหลายเหลี่ยมแต่ละอันและอีกคอลัมน์หนึ่งด้วย ความยาวเส้นรอบวงของรูปหลายเหลี่ยมแต่ละอัน?

13 qgis  polygon  area 

ฉันกำลังมองหาที่จะรวมจำนวนของแต่ละตารางลงในตารางใหม่ใน PostGIS นี่เป็นงานง่ายเมื่อทำงานกับข้อมูล Shapefile แต่ฉันไม่แน่ใจว่าจะทำอย่างไรใน PostGIS ความช่วยเหลือใด ๆ ที่จะได้รับการชื่นชมอย่างมาก. ฉันคิดว่าฉันใช้ผนวก แต่ฉันไม่แน่ใจว่าจะเริ่มต้นที่ไหน

13 postgis  merge 

ฉันพยายามประมวลผลข้อมูลแรสเตอร์โดยใช้ ogr / gdal และฉันไม่สามารถใช้ประโยชน์จากคอร์ทั้งหมดในเครื่องได้ เมื่อฉันเรียกใช้กระบวนการบนแกนเดียวฉันได้รับประโยชน์ 100% ของแกนนั้น เมื่อฉันพยายามที่จะแบ่งออกเป็นมัลติคอร์ (ในตัวอย่างด้านล่างโดย chunking x ชดเชยและวางไว้ในคิว) ฉันได้รับการใช้ประโยชน์ที่น่าสงสารในแต่ละ 8 แกนของฉัน ดูเหมือนว่าจะเพิ่มการใช้ประโยชน์ได้สูงสุด 100% ในแต่ละคอร์ (เช่น 12.5% ​​สำหรับแต่ละคอร์) ฉันกังวลว่าการใช้แหล่งข้อมูลเดียวกันคือคอขวด แต่ฉันก็ทำซ้ำไฟล์แรสเตอร์ต้นแบบสำหรับแต่ละคอร์ ... และการใช้งานคอร์ยังคงอึอยู่ สิ่งนี้ทำให้ฉันเชื่อว่า ogr หรือ gdal ทำตัวเหมือนทรัพยากรที่ใช้ร่วมกันในคอขวด แต่ฉันไม่พบสิ่งใดทางออนไลน์เกี่ยวกับเรื่องนั้น ความช่วยเหลือใด ๆ ที่จะได้รับการชื่นชมมาก! นี่คือฟังก์ชัน "ผู้ช่วยเหลือ" ที่ทำงานภายในแต่ละเธรดผู้ทำงาน: def find_pixels_intersect_helper(datasource, bounds_wkt, x_min, x_max): bounds = ogr.CreateGeometryFromWkt(bounds_wkt) rows_to_write = [] for …

13 python  gdal  ogr  parallel-processing 

ฉันกำลังพยายามแสดงคุณสมบัติบางอย่างในเลเยอร์ที่อิงกับ Shapefile โดยใช้ QGIS ฉันได้พบฟังก์ชั่นชุดย่อยในแผงคุณสมบัติซึ่งทำสิ่งที่ฉันต้องการ แต่น่าเบื่อเกินไปเพราะฉันต้องร่อนด้วยตนเองผ่านคุณสมบัติจำนวนมาก มีวิธีที่เร็วกว่าในการทำเช่นนี้? ฉันใช้ทั้งรูปร่างและ PostGIS ที่ดีที่สุด Dietmar

13 qgis  postgis  query  rendering  filter 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา มีแหล่งที่มาของแผนที่อินเทอร์เน็ต (ซูมได้ปูกระเบื้องเช่น "google like") ที่มีแม่น้ำหลายสายหรือไม่? Ie ติดป้ายอย่างดีแม้ในเครื่องชั่งขนาดใหญ่? ดังนั้นที่ผมยกตัวอย่างเช่นสามารถที่จะติดตามแม่น้ำ Donau (มีป้าย) โยเช่นนี้ แผนที่ google ขั้นพื้นฐานนั้นแย่มากแม่น้ำไม่แสดงผลดีมักสูญหายในเครื่องชั่งขนาดใหญ่และในกรณี 90% ที่ไม่มีป้ายกำกับ OpenStreetMap - คล้ายกันยิ่งกว่านั้นพวกเขายังระบุว่าถนนเป็นสีฟ้า (!!!) ซึ่งสร้างความสับสนมากยิ่งขึ้น! (การปฏิบัติที่เลวร้ายจริงๆ) OpenCycleMap - ยิ่งแย่กว่า ...

13 data  web-mapping  layer-resource 

ฉันกำลังมองหาฐานข้อมูลที่มีสนามบินโลกทั้งหมด (มากที่สุด ... มากที่สุดเท่าที่จะทำได้) ฉันรู้ว่ามีชุดดังกล่าวอยู่หนึ่งชุด - http://www.world-airport-database.com/database.php - แต่ข้อตกลงระบุว่ามีไว้สำหรับการใช้ภายในเท่านั้น ... ชนิดของการเอาชนะวัตถุประสงค์ ข้อมูลจะต้องมีรหัสสนามบิน, ชื่อ, สถานที่ (ประเทศ, เมือง, รัฐ / สถานที่, ไปรษณีย์ / ไปรษณีย์, lat / lon) คุณรู้ชุดดังกล่าวหรือไม่? แน่นอนว่ายินดีที่จะซื้อชุดดังกล่าวฟรีก็ยอดเยี่ยมเช่นกัน

13 data  transportation 

ฉันกำลังเขียนเว็บ API ไปยังข้อมูล GIS ของเราซึ่งปัจจุบันกำลังปรับปรุงบริการการเข้ารหัสภูมิศาสตร์ของเรา ฉันกำลังสร้างตัวเลือกสำหรับผู้ใช้เพื่อระบุประเภทผลลัพธ์เป็น GeoJSON ดังนั้นฉันจึงพยายามทำตามข้อมูลจำเพาะ ข้อมูลของเราถูกเก็บไว้ใน nad83 datum UTM zone 12 north projection GeoJSON สเปคกล่าวว่าการเพิ่ม CRS wkid คือ 26912 แต่ OGC CRS URN คืออะไร? ฉันสามารถเข้าใจได้ว่า OGC CRS URN ถูกแบ่งออกเป็นสองส่วน โกศ: OGC: def: CRS: OGC: 1.3: CRS84 urn เป็นตัวระบุ ogc เป็นองค์กร def เป็นอีกข้อตกลงคงที่ crs เป็นประเภท (ระบบอ้างอิงพิกัด), OGC เป็นหน่วยงานที่มีอำนาจ …

13 coordinate-system  geojson  ogc  specification 

ฉันรู้ว่าการที่จะหาที่ 5 หลักเขตแดนรหัสไปรษณีย์ มีขอบเขตรหัสไปรษณีย์ 9 หลักอยู่หรือไม่ ถ้าเป็นเช่นนั้นฉันจะหาได้ที่ไหน

13 data  reverse-geocoding  zip-codes 

เมื่อฉันประมวลผลแรสเตอร์ใน QGIS มันจะสร้างแถบที่สอง ปรากฏว่าแถบนั้นไม่มีข้อมูลที่เป็นประโยชน์ ฉันชอบใช้ Raster Colors Tool และจะไม่ยอมรับแรสเตอร์ที่มีอะไรมากกว่าในวง ฉันจะลบวงเสริมได้อย่างไร

13 raster  qgis 

ฉันต้องการเลือกตามเดือนที่เฉพาะเจาะจงจากเขตข้อมูลวันที่ / เวลาของฉันใน ArcGIS และดูเหมือนจะไม่พบคำสั่งที่ใช้งานได้ ฉันมีหลายปีและหลายวัน แต่ต้องการกรองออกและเลือกบันทึกจากเดือนพฤษภาคมเท่านั้น

13 arcgis-desktop  select-by-attribute  date 

เมื่อใช้ OpenLayers ฉันเพิ่มเลเยอร์ WFS (บน GeoServer) ด้วยตัวกรองที่ส่งคืนฟีเจอร์ทั้งหมด (สีดำ) ที่ตัดกันรูปหลายเหลี่ยมของฉัน (สีเหลือง) ที่วางอยู่เหนือประเทศละตินอเมริกาบางแห่งภายในวันที่กำหนด อย่างไรก็ตามคุณสมบัติที่ข้ามในแนวนอนบนแผนที่ไม่ได้ตัดกับรูปหลายเหลี่ยมของฉัน คุณลักษณะนี้อยู่ที่ไหนสักแห่งในมหาสมุทรแปซิฟิกระหว่างฮาวายและฟิจิและไม่ใช่ในละตินอเมริกา ปัญหาคือแทนที่จะข้ามเส้นวันที่ระหว่างประเทศกำลังแสดงผลบนแผนที่โดยล้อมรอบทั้งโลก คุณสมบัติ problamatic ถูกกำหนด: POLYGON ((- 179.700417 14.202717, -178.687422 13.992875,179.024138 8.24716, -179.98241 8.035567, -179.700417 14.202717) ฉันมีคุณสมบัติของเส้นวันที่ที่เป็นปัญหาเช่นนี้ แต่มีการ จำกัด ให้แคบลงเป็นตัวอย่างนี้ ฉันไม่สามารถเพิกเฉยได้ในแอปพลิเคชันของฉันเพราะฉันมีหลายอย่าง ฉันได้ลองใช้ "wrapDateLine: true" ในชั้นฐานและชั้น WFS ด้วยผลลัพธ์เดียวกัน ไม่แน่ใจว่านี่จะเป็นปัญหา GeoServer หรือปัญหา OpenLayers ไม่มีใครรู้วิธีการแก้ปัญหาบรรทัดวันที่ระหว่างประเทศของฉัน?

13 geoserver  openlayers-2  wfs  antimeridian