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

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

10 qgis  vector  style  masking  qgis-custom-function 

ฉันมีถนน (เส้นกลางถนน) และสร้างตารางหลายเหลี่ยมในฐานข้อมูล PostgreSQL ของฉัน สถานการณ์ตัวอย่างเป็นดังนี้: ปัญหา: ฉันต้องคำนวณเส้นขนานตามถนนที่จุดตัดของบัฟเฟอร์ 50 เมตรรอบถนนและรูปหลายเหลี่ยมอาคารที่ใกล้ที่สุดทั้งสองด้าน สถานการณ์เอาท์พุทที่ต้องการคือ: สิ่งที่ฉันได้ลอง: แนวทางของฉันคือ: 1) Generate 50m buffer around street layer 2) get the intersection of buffer and polygons 3) Compute the distance 4) Draw offset curves (parallel lines) at both sides of street layer 5) Merge both curves to get …

10 postgis  postgresql 

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

10 qgis  print-composer  atlas 

ฉันกำลังพยายามกำหนด ortelius oval projection และ Nicolosi globular projection ใน QGIS โดยใช้ proj4 foramt ตามที่อธิบายไว้ในเอกสารประกอบ แต่ QGIS ไม่สามารถฉายภาพในการฉายภาพเหล่านี้แม้ว่าพวกเขาจะได้รับการสนับสนุนในห้องสมุด proj4 ... แผนที่ ไม่มีใครรู้ว่าสิ่งที่อาจเป็นปัญหาหรือไม่

10 qgis  proj 

เมื่อเร็ว ๆ นี้ฉันได้มองไปรอบ ๆpg_stat_user_tablesและรู้สึกประหลาดใจที่ได้เห็นการสแกนตามลำดับจำนวนมากในตารางอวกาศของฉัน แน่นอนตารางเหล่านี้ขาดดัชนีเชิงพื้นที่ ฉันจะค้นหาตารางทั้งหมดที่มีคอลัมน์รูปทรงเรขาคณิตที่ไม่ได้จัดทำดัชนีได้อย่างไร

10 postgis  postgresql 

คำนึงถึงพิกัด 31.96212, -103.004715 แปลง UTM ให้มันพิกัด UTM 13/R/FRเป็น ตัวแปลงตัวอย่างอยู่ที่นี่: http://www.rcn.montana.edu/resources/converter.aspx แต่มีหลายคนและพวกเขาให้คำตอบที่คล้ายกันสำหรับพิกัดเหล่านี้ พร้อมกันในชุดข้อมูล Sentinel-2 ที่นี่http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/ ฉันไม่พบFRไดเรกทอรีย่อย ใน google สถานที่นี้อยู่ที่นี่: และการค้นหาสถานที่เดียวกันในเบราว์เซอร์รูปภาพ Sentinel ที่ฉันเห็นกระเบื้องนั้นแตกต่างกัน ซึ่งหมายถึง13/S/FRเช่นเดียวกันUTMและสี่เหลี่ยม แต่วงดนตรีที่แตกต่างกัน เป็นไปได้อย่างไร? UPDATE KML ที่มีไทล์ Sentinel-2 จะรายงานSไทล์ในตำแหน่งที่กำหนด อัพเดท 2 ตามรูปนี้ นำมาจากที่นี่ที่FRตารางตั้งอยู่ครึ่งหนึ่งในSUTM โซนและครึ่งหนึ่งในRเขต เห็นได้ชัดว่าตัวแปลงอัตโนมัติส่วนใหญ่จะกำหนดตารางนี้ให้กับRโซนในขณะที่ Sentinel-2 จะพิจารณาSพื้นที่ ที่นี่มีความจริงหรือไม่? อัพเดท 3 รหัส Python ง่าย ๆ นำมาจากที่นี่https://gis.stackexchange.com/a/224994/32207 bandVals = "CDEFGHJKLMNPQRSTUVWXX" lon = …

10 coordinate-system  utm  sentinel-2  mgrs 

เมื่อฉันกำลังละลายรูปร่างไฟล์: ogr2ogr "b:\out.shp" "b:\Tele.shp" -dialect sqlite -sql "SELECT ST_Union(geometry), nev, jar FROM Tele GROUP BY jar" ฉันได้รับข้อความแสดงข้อผิดพลาดเช่น: คำเตือน 1: อักขระหนึ่งหรือหลายตัวไม่สามารถแปลงได้อย่างถูกต้องจาก UTF-8 ถึง ISO-8859-1 ฉันไม่เข้าใจเพราะ Tele.cpg มี "UTF-8" เมื่อลองทำสิ่งนี้: ogr2ogr "b:\out.shp" "b:\Tele.shp" -dialect sqlite -sql "SELECT ST_Union(geometry), nev, jar FROM Tele GROUP BY jar" ENCODING=UTF-8 ข้อความแสดงข้อผิดพลาดคือ: คำเตือน 1: ชื่อเลเยอร์จะถูกละเว้นพร้อมกับ -sql ข้อผิดพลาด …

10 ogr2ogr  utf-8  encode 

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

10 qgis  buffer 

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

10 qgis  symbology  style  visualisation  qgis-custom-function 

ใน QGIS right-clickคุณสามารถสลับเปิด / ปิดบางพื้นที่ในแผงควบคุมแถบเครื่องมือดังแสดงในภาพด้านล่าง (แผงข้อมูล GPS เป็นต้น) ฉันต้องทำสิ่งนี้โดยทางโปรแกรมดังนั้นจะสลับเปิด / ปิดพาเนลจาก PyQGIS ได้อย่างไร ฉันพบว่าQgsInterfaceชั้นเรียนมีบางอย่างaddDockWidgetและremoveDockWidgetวิธีการ แต่ฉันไม่สามารถคิดออกว่าเป็นวิธีที่เหมาะสมเพื่อให้บรรลุนี้ ...

10 pyqgis  panel 

ฉันมักจะสลับระหว่างตัวสร้างแบบกราฟิกและแอปพลิเคชัน QGIS หลักและสงสัยว่าทำไมฉันไม่สามารถใช้งาน QGIS ได้เมื่อเปิดหน้าต่างตัวสร้างแบบจำลอง การบันทึกแบบจำลองของฉันและการปิด / เปิดตัวสร้างแบบจำลองทำให้รำคาญฉันดังนั้น ... เหตุผลที่ลึกซึ้งกว่าสำหรับผู้สร้างแบบจำลองนั้นเป็นกิริยาช่วยหรือไม่? มีวิธีง่ายๆในการปิดพฤติกรรมนี้ (คอนโซลหลาม)?

10 qgis  qgis-modeler  modal 

ฉันมีปัญหาเล็กน้อยที่นี่กับสถานการณ์นี้ ฉันต้องการ "NumeroCamp" เท่านั้นที่จะแสดงตัวเลขใน "ข้อความ" มีฟังก์ชั่นใด ๆ ในเครื่องคำนวณภาคสนามที่สามารถช่วยฉันได้หรือไม่?

10 qgis  field-calculator  string 

ฉันมีชั้น MultiLineStringZM ในฐานข้อมูล sqlite และฉันพยายามที่จะเห็นภาพการวัดหรือค่า m ที่จุดสุดยอด ฉันได้ลองค้นหาข้อมูลเกี่ยวกับวิธีการทำสิ่งนี้ใน QGIS และสิ่งที่ฉันสามารถรวบรวมได้คือสิ่งนี้ไม่สามารถทำได้โดยตรงจากเลเยอร์ linestring และจุดนั้นต้องแยกออกเป็นเลเยอร์อื่น ฉันใช้Vector-> Geometry Tools-> Extract nodesเพื่อสร้างชั้นหลายจุดที่เป็นจุดสุดยอดของชั้น multlinestring ของฉัน แต่กระบวนการสูญเสียค่า m ของจุดยอด ฉันต้องการ m-values ​​ที่เก็บรักษาไว้โดยการบันทึก m-value เป็นคุณลักษณะของจุดหรืออย่างอื่น? ภายในเรามีเครื่องมือบรรทัดคำสั่งที่แปลง linestrings เป็น shapefile จุดด้วยค่า m ที่เก็บไว้เป็นแอตทริบิวต์ในแต่ละจุดและฉันใช้มันเพื่อตรวจสอบว่ามีค่า m- กำหนดให้กับจุดยอดและฉันสามารถใช้ ถ้าฉันต้องทำ แต่ถ้าเป็นไปได้มันจะดีถ้าทำสิ่งนี้ได้โดยตรงภายใน QGIS แก้ไข - ทำซ้ำสิ่งที่ฉันได้กล่าวข้างต้น แต่เน้นความจริงที่ว่าเรามีเครื่องมือบรรทัดคำสั่งที่สามารถบรรลุผลลัพธ์ที่ฉันกำลังมองหาที่ใช้ไลบรารี GDAL ดังนั้นวิธีแก้ปัญหาแสดงคำตอบเพียงบางส่วนใน PyQGIS ไม่ใช่คำตอบที่ฉันกำลังมองหา ฉันกำลังมองหาเครื่องมือในตัวปลั๊กอินที่พร้อมสำหรับ QGIS หรือสคริปต์ที่สมบูรณ์ที่สามารถแยก …

10 qgis  linestring  linear-referencing  extract  m-values 

ฉันจะสร้างบัฟเฟอร์ที่ละลายจากคุณสมบัติอินพุตหลายจุด ในตัวอย่างด้านล่างตารางอินพุตมี 4 ฟีเจอร์ คุณลักษณะ#2ประกอบด้วยรูปทรงเรขาคณิตสองจุด หลังจากสร้างบัฟเฟอร์ฉันได้รับรูปเรขาคณิต 4 เหลี่ยม: มีวิธีการจัดกลุ่มผลลัพธ์หรือไม่ บัฟเฟอร์ของจุด#1และ#2ถูกยุบและควรเป็นคุณลักษณะหลายรูปหลายเหลี่ยม ( a) สิ่งที่ฉันทำไปแล้ว: -- collect all buffers to a single multi-polygon feature -- dissolve overlapping polygon geometries CREATE TABLE public.pg_multibuffer AS SELECT row_number() over() AS gid, sub_qry.* FROM (SELECT ST_Union(ST_Buffer(geom, 1000, 8))::geometry(MultiPolygon, /*SRID*/) AS geom FROM public.multipoints) AS sub_qry; แก้ไข: …

10 postgis  buffer  multipart 

มือใหม่ที่นี่ต้องดิ้นรนกับ GIS ฉันพยายามที่จะแผนที่ออกวอร์ดสำหรับเมือง Milwuakee โดยใช้ shapefiles พบได้บนเว็บไซต์ของพวกเขาเขตเว็บไซต์เขต ฉันติดตามกระทู้ที่นี่ด้วยความสำเร็จ รหัสของฉันให้: from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire world, wgs 84...not obvious inProj = Proj(init='epsg:32054') outProj = Proj(init='epsg:4326') x1,y1 = 2560131.496875003, 406816.434375003 x2,y2 = transform(inProj,outProj,x1,y1) print(x2,y2) กับเอาท์พุท -65.70220967836329 43.08590211722421 ปัญหานี้เป็นสิ่งที่ผิด The lon / lat สำหรับมิลวอกีคือ …

10 python  shapefile  fiona  pyproj