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

ฉันกำลังตรวจสอบวิธีที่ดีที่สุดในการตัดการเชื่อมต่อตามจุด สถานการณ์คือ: ถนนจำนวนมากต้องแบ่งส่วนตามจุดตัดกันแบบนี้: ฉันได้ ตารางการเคลื่อนไหว (เต็มจุดโดยไม่ได้คะแนน) ตารางคะแนน st_intersection ฉันต้องได้รับการตัดส่วนที่เป็นอิสระโดยตารางจุดตัด ฉันใช้ฟังก์ชั่น PostGIS และพบวิธีการต่าง ๆ แต่ทุกคนมีปัญหาบางอย่างให้ฉัน นี่คือสิ่งที่ฉันทดสอบแล้ว: 1 ตารางเส้น: 1 แถว, st_memunion ของ 1200 บรรทัดตารางจุด: 1,700 แถว (คะแนน) มีอะไรที่ไม่ดี: ใช้เวลาและหน่วยความจำอย่างมาก ไม่สามารถสร้างหลายตารางพร้อมกันทำให้หน่วยความจำไม่สามารถจัดการได้ และผลลัพธ์ก็สกปรกและยุ่งเหยิง แทนที่จะให้หมายเลขแถวที่ถูกต้องกับฉันและฉันจำเป็นต้องทำความสะอาดในภายหลัง (อธิบายอย่างดีที่นี่แยกบรรทัดที่จุดตัดกัน ) CREATE TABLE lines_with_messy_result AS ( SELECT ((ST_DUMP(ST_SPLIT(a.geom,b.ix))).geom) as geom FROM st_union_lines a INNER JOIN lots_of_points b ON …

10 postgis  postgresql  openstreetmap  splitting  streets 

ฉันสนใจที่จะคำนวณกราฟเส้นสำหรับกราฟอินพุตโดยเฉพาะอย่างยิ่งใน GRASS เนื่องจากฉันจะใช้เครื่องมือวิเคราะห์เครือข่าย GRASSอื่น ๆ ในภายหลัง ฉันไม่พบการอ้างอิงถึงกราฟเส้นในเอกสาร GRASS มีคำอื่นที่ฉันควรใช้ในการค้นหาหรือฟังก์ชันซ่อนเร้นที่จะสร้างกราฟเช่นนี้หรือไม่? สำหรับบริบท: กราฟเส้นของกราฟที่ไม่ได้บอกทิศทาง G คือกราฟ L (G) อีกอันที่แสดงถึงคำคุณศัพท์ระหว่างขอบของ G ที่มา: http://en.wikipedia.org/wiki/Line_graph

10 grass  network  graph 

ฉันกำลังทำงานในด้านของไอโซนีและอัลกอริธึมพื้นฐาน ตอนนี้สิ่งที่ทำให้เกิดปัญหาไม่ใช่การคำนวณถ้า isochrone ตัวเอง แต่การสร้างภาพของผลลัพธ์ ผลลัพธ์ของอัลกอริทึม isochrone ของฉันคือจุดและขอบ ในความเป็นจริงฉันมีวิธีแก้ปัญหาการทำงาน แต่สำหรับ 3873 edge และสำหรับ 1529 nodes สิ่งต่าง ๆ ดูเหมือนจะใช้ตลอดไป (ประมาณ 2.0 วินาทีบนแล็ปท็อป Lenovo T440 ของฉันซึ่งมี CPU Core i7 และ SSD ความเร็วสูง) แทนที่จะวินาทีฉันต้องการอะไรมากกว่าเช่น msec :-) บางทีใครบางคนสามารถช่วยฉันในการลดเวลาการคำนวณที่จำเป็นในการสร้างรูปหลายเหลี่ยมซึ่งมองเห็นพื้นที่ที่สามารถเข้าถึงได้ แต่รอ ... สิ่งแรกก่อน! นี่คือการมองเห็นของขอบที่ฉันเป็นผลการคำนวณของ isochrone ของฉัน: ขอบเหล่านี้จะถูกเก็บไว้ในตารางฐานข้อมูล PostGIS และเป็น linestrings ง่าย สิ่งที่ฉันต้องการแสดงต่อผู้ใช้มีลักษณะดังนี้: สังเกตบริเวณที่ไม่ได้เชื่อมต่อในทิศใต้และทิศตะวันออกของภาพ สิ่งเหล่านี้ควรวาดเป็นพื้นที่แยกต่างหาก (ดังนั้นจึงไม่อนุญาตการรวมที่นี่ :-)) …

10 postgis  postgis-2.0  polygonize 

ฉันได้แปลงรูปร่างไฟล์ที่มีคุณสมบัติเป็น. kml สำหรับลูกค้า จากนั้น. kml และแอตทริบิวต์จะถูกแก้ไขภายใน google earth ฉันได้รับ. kml ที่อัปเดตแล้วและตอนนี้ฉันต้องกลับไปที่. shp ปัญหาคือ: คุณลักษณะทั้งหมดของฉันตอนนี้อยู่ในรูปแบบ html ซึ่งทำให้กล่องข้อมูลป๊อปอัปใน Google Earth เมื่อฉันพยายามแปลงกลับเป็น. shp คุณลักษณะของฉันจะถูกรวมเป็นฟิลด์ (สร้างโดย Google) "คำอธิบาย" คุณลักษณะและข้อมูลของฉันทั้งหมดถูกรวมเข้าในเซลล์เดียว มีวิธีดึงข้อมูลจากเขตข้อมูล "คำอธิบาย" แม้ว่าฉันจะทำเช่นนั้นมีวิธีที่ง่ายในการกำจัดรูปแบบ html และทำให้ใช้งานได้หรือไม่ นี่คือสิ่งที่ฉันใช้โดยใช้ "ตารางฟิวชั่น" ที่ฉันพบในบางฟอรัม คุณสามารถเห็น html ทางด้านซ้ายมันจะดำเนินต่อไปมากกว่าที่จะแสดง ฉันใช้ ArcGIS Desktop 10.2

10 arcgis-desktop  shapefile  kml  convert  google-earth 

บางครั้งฉันต้องจัดทำแผนที่เพื่อแสดงจุดที่น่าสนใจ ขั้นตอนแรกในการสร้างหน้าโดยใช้ mesh ปกติ: ฉันไม่ชอบวิธีแก้ปัญหาเพราะก) มีบางหน้าเว็บที่มีจุดเดียว (เช่นหน้า 25) อยู่ที่ขอบและ b) หน้ามากเกินไป ปัญหาแรกนั้นง่ายต่อการแก้ไขโดยใช้รหัส - ย้ายสี่เหลี่ยมผืนผ้าของขอบเขตหน้าไปที่กึ่งกลางของจุดที่เกี่ยวข้อง: ฉันยังคงไม่ชอบมันดูแออัดมากเพราะจำนวนหน้ายังคงเท่าเดิม จำไว้ว่าพวกเขาทั้งหมดกลายเป็นหน้ากระดาษ A3 จริงในรายงานหลายฉบับ! ดังนั้นฉันจึงจัดทำโค้ดเพื่อลดจำนวนหน้า ในตัวอย่างนี้จาก 45 ถึง 34 ฉันไม่แน่ใจว่านี่เป็นผลลัพธ์ที่ดีที่สุดที่สามารถทำได้หรือไม่ กลยุทธ์ที่ดีที่สุดคืออะไร (รหัสเทียม, สิ่งพิมพ์, ห้องสมุด Python) เพื่อสับเปลี่ยนจุดเพื่อลดจำนวนสี่เหลี่ยมขนาดที่กำหนดไว้เพื่อจับจุดทั้งหมด แน่นอนว่ามีคนค้นพบมันในทฤษฎีเกมศิลปะการทหารหรืออุตสาหกรรมการประมง นี่คือการปรับปรุงคำถามเดิม: สิ่งนี้แสดงขอบเขตและขนาดหน้าจริงที่ต้องการ: การย่อ / ขยายแสดง 10 จาก 164 หน้า: คลาสคุณลักษณะจุดตัวอย่าง ขนาดสี่เหลี่ยมผืนผ้าสามารถเปลี่ยนได้ทันทีที่มันอยู่ในขอบเขตที่ จำกัด เช่นขนาดเล็กก็ดี

10 arcgis-desktop  arcpy  algorithm  data-driven-pages  maximum-coverage-problem 

ฉันกำลังพยายามสร้างบรรทัดผู้นำแบบไดนามิกโดยใช้มุมมอง PostGIS เพิ่มเติมจากเครื่องมือ QGIS„ ย้ายป้ายกำกับ“ CREATE VIEW leader_line AS SELECT gid, ST_MakeLine(geom, ST_SetSRID(ST_MakePoint(xcord_label, ycord_label), SRID))::geometry(linestring, SRID) AS geom FROM point WHERE xcord_label IS NOT NULL; นี้ทำงานได้ดีสำหรับป้ายกำกับทั้งหมดแต่สร้างผิดมองเส้นผู้นำสำหรับป้ายชื่อWHERE ST_X(geom) < xcord_labelWHERE ST_X(geom) > xcord_label ไม่มีใครรู้วิธีการวางสายผู้นำสำหรับป้ายกำกับอย่างเหมาะสมWHERE ST_X(geom) > xcord_labelหรือไม่ มีวิธีใดบ้างในการอ้างถึงพิกัด xmax ของฉลาก?

10 qgis  postgis  labeling  coordinates  dynamic-layer 

ฉันมีแผนที่แผ่นพับที่ฉันวาง SVG "ด้านบนของ" โดยใช้บานหน้าต่างแผ่นพับซ้อนทับ ฉันต้องการเพิ่มแปรง D3 ลงในบานหน้าต่างวางซ้อนที่จะอัปเดตแผนภูมิแยกต่างหาก ฉันได้ดูตัวอย่างของตัวอย่าง D3.brush () และตัวอย่างแผ่นพับใด ๆ แต่ไม่พบสิ่งใดเลย สิ่งที่ใกล้เคียงที่สุดที่ฉันพบคือhttps://github.com/mbostock/d3/issues/1321และhttps://www.mapbox.com/mapbox.js/example/v1.0.0/timeline-scaled-markers/ (นี่เป็นการแปรงแผนภูมิไม่ใช่แผนที่) ฉันสงสัยว่า: ใครรู้ตัวอย่างหรือ ใครรู้ว่าสิ่งนี้เป็นไปได้ / ตรงไปตรงมาสมเหตุสมผลหรือ ถ้าเป็นไปได้ทุกคนมีคำแนะนำเกี่ยวกับปัญหาที่ฉันอาจพบ

10 leaflet  d3 

แนวคิดของคอมโพสิตฮิลล์เชดคือการรวมหลายฮิลส์เชดกับการวางแนวแสงแดดที่แตกต่างกันเพื่อหลีกเลี่ยงจุดบอด มันต้องรวมภาพฮิลส์ 3 ภาพด้วย315⁰ (แสงแดด NW, เลเยอร์เริ่มต้น), + 355⁰เลเยอร์, ​​+ 275⁰แสงตะวันตามลำดับ นี้รายละเอียดเพิ่มขึ้นและความสง่างามเต็มรูปแบบที่ไม่ใช่ GIS ดังแสดง below.A คำอธิบายอยู่ที่นี่ การใช้ไฟล์ srtm ใด ๆ ฮิลล์เชดเริ่มต้น 3 ตัวสามารถทำได้ผ่าน: gdaldem hillshade input.tif hillshades_A.tmp.tif -s 111120 -z 5 -az 315 -alt 60 -compute_edges gdaldem hillshade input.tif hillshades_B.tmp.tif -s 111120 -z 5 -az 355 -alt 60 -compute_edges gdaldem …

10 gdal  dem  hillshade  gdal-calc 

เนื่องจาก QGIS ใช้ GDAL แผนที่ใดที่ส่งออกเป็น PDF กลายเป็น Geospatial PDF หรือไม่ ฉันหวังว่าจะนำเข้าแผนที่จาก QGIS ไปยังแอพ Avenza PDF Map ในโทรศัพท์ของฉัน โพสต์ก่อนหน้ากล่าวถึงการส่งออกเป็น Geospatial PDF และ GeoTiff แต่ไม่มีคำตอบหรือกระบวนการที่แน่นอน

10 qgis  geotiff-tiff  pdf  geospatial-pdf  avenza 

ในฐานะที่เป็นความคิดเห็นในโพสต์บล็อกของ Esri ที่มีชื่อว่าการประกาศโปรแกรม ArcGIS Pro เบต้าเมื่อเดือนพฤษภาคมที่ผ่านมามีคนชื่อ mahenderraj ถามว่า: ฉันคิดว่ามันเป็นไปไม่ได้ที่จะแชร์แผนที่เป็น mapservice ไปยังเซิร์ฟเวอร์ ArcGIS จาก ArcGIS Pro เพื่อที่เราจะต้องใช้ผลิตภัณฑ์ ArcMap Desktop ปกติ ฉันดูเหมือนจะไม่ได้เห็นคำตอบนี้และในอินเทอร์เฟซ ArcGIS Pro 1.0 มันดูเหมือนง่ายต่อการเผยแพร่แผนที่ไปยัง ArcGIS Online หรือพอร์ทัลสำหรับ ArcGIS แต่ ArcGIS Server ไม่ได้กล่าวถึง เป็นไปได้ไหมที่จะเผยแพร่แผนที่จาก ArcGIS Pro ไปยัง ArcGIS Server เป็น Map Service

10 arcgis-server  map-service  arcgis-pro 

เมื่อฉันเปิด QGIS เพิ่มเลเยอร์และคำนวณพื้นที่ของ shapefile ผ่านเครื่องคิดเลขภาคสนามฉันได้พื้นที่ที่แตกต่างจากเมื่อฉันเปิด QGIS และตรวจสอบ "เปิดใช้งานการแปลง CRS ทันที" และคำนวณพื้นที่ นี่คือการทำให้แน่ใจว่าโครงการและเลเยอร์มีระบบพิกัดเดียวกัน (หมายเลข EPSG เดียวกัน) ผมทำอะไรผิดหรือเปล่า? ฉันมีรูปร่างไฟล์ที่มีการคำนวณพื้นที่ที่ทำกับ ArcGIS (ไม่ใช่ฉันข้อมูลถูกส่งมาให้ฉันและฉันไม่มีเงื่อนงำที่ CRS พื้นที่คำนวณด้วย ArcGIS) เลเยอร์ไฟล์ CRS คือ EPSG: 21781 (สวิตเซอร์แลนด์) ใน QGIS ถ้าฉันไม่เปลี่ยนการตั้งค่า OTF และปล่อยให้ CRS โครงการเป็น EPSG: 4326 (WGS84) ฉันได้รับค่าเดียวกันกับค่าพื้นที่ ArcGIS อย่างไรก็ตามถ้าฉันเปลี่ยน OTF ก่อนเพิ่มเลเยอร์เป็น EPSG: 21781 ฉันจะได้รับค่าพื้นที่ที่แตกต่างกัน ดังที่ฉันเข้าใจสิ่งนี้แนะนำว่า ArcGIS Area คำนวณด้วย …

10 qgis  coordinate-system  area 

ฉันมี orthoimage RGB ที่มีพื้นที่สีดำแทน NoData ดังภาพที่นี่: เมื่อฉันใช้เครื่องมือ SetNull มันจะลบพื้นที่สีดำออกไป แต่มันจะเปลี่ยนเป็นสีเทาแบบออร์โธไมซ์ซึ่งไม่ใช่สิ่งที่ฉันต้องการ ฉันรู้ว่ามีวิธีแก้ไขปัญหานี้โดยใช้เครื่องมือ Clip (Data Management) เพื่อคลิปแรสเตอร์ให้กับตัวเองและในนั้นคุณสามารถตั้งค่า NoData เป็นค่าของเซลล์สีดำ (0) อย่างไรก็ตามฉันกำลังสร้างเวิร์กช็อปสำหรับนักเรียนและฉันต้องการใช้วิธีการที่ใช้งานง่ายกว่า แน่นอนว่ามีตัวเลือกในการเรียกใช้เครื่องมือ NoData ในแต่ละแบนด์แยกจากกันแล้วรวมสามแบนด์เข้าด้วยกัน แต่ฉันต้องมีเงื่อนไขที่ band1, band2 และ band3 เท่ากับ 0 ทั้งหมดเพื่อตั้งเป็น NoData ความสามารถในการทำเช่นนี้ใน ModelBuilder จะเหมาะ แต่การเขียนสคริปต์โดยใช้ ArcPy ก็จะดีเช่นกัน ฉันกำลังใช้เดสก์ท็อป ArcGIS 10.2 พร้อมใบอนุญาตขั้นสูง

10 arcgis-desktop  raster  arcgis-10.2  nodata  multi-band 

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

10 qgis  python  intersection  line  crossing 

ฉันได้ให้ภาพ Geotiff และข้อมูล Lidar ที่สอดคล้องกัน (x, y, z) ในพิกัด UTM ฉันต้องรวมข้อมูล Lidar เข้ากับค่า RGB จากภาพ นั่นหมายความว่าในตอนท้ายฉันต้องพล็อต (3D) แต่ละจุดของสีคลาวด์ LiDAR ที่เข้ารหัสด้วยค่า RGB ที่สอดคล้องกันจากภาพ Geotiff ฉันแปลงข้อมูล Lidar เป็น shapefile โดยใช้ QGIS ฉันควรทำอย่างไรต่อไป ใน R ฉันลองใช้plot3Dฟังก์ชั่น แต่มันไม่ทำงาน ฉันกำลังแนบdoc text , shapefile และอิมเมจ tif แก้ไข: ฉันได้ทำโปรแกรมต่อไปนี้ตามที่แสดงด้านล่าง: require(raster) require(maptools) # to take shape files #require(car) …

10 raster  r  3d  lidar  3d-model 

ฉันมีเลเยอร์บรรทัดที่คาดการณ์ (Shapefile) ใน QGIS สำหรับทุกส่วนของเส้นฉันต้องการสร้างรูปหลายเหลี่ยมโดยไม่มีการปิดท้ายบรรทัด Polygone ที่เป็นผลลัพธ์ไม่ควรถูกยุบไปเป็นอันเดียว ฉันลองใช้เครื่องมือบัฟเฟอร์ fTools สิ่งนี้สามารถสร้างบัฟเฟอร์ได้โดยไม่ต้องละลาย แต่เป็นไปไม่ได้ที่จะกำหนดแคปสิ้นสุดบรรทัด ฉันลอง v.buffer ของ GRASS ด้วย สิ่งนี้สามารถสร้างบัฟเฟอร์ด้วยสิ้นสุดบรรทัดที่กำหนด แต่ฉันไม่สามารถปิดการละลาย มีความเป็นไปได้ที่จะมีทั้งสองอย่าง: ไม่ละลายและไม่มีตัวพิมพ์ใหญ่ที่ปลายบรรทัด? ผลลัพธ์ควรเป็น Shapefile ที่มีรูปหลายเหลี่ยม (บัฟเฟอร์) สำหรับทุกส่วนของบรรทัด

10 qgis  buffer  dissolve