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

ถามตอบสำหรับนักทำแผนที่นักภูมิศาสตร์และผู้เชี่ยวชาญด้าน GIS

1
การอ่านการแก้ไขและการเขียนพิกัดด้วย GDAL ในหลาม
ฉันพยายามเรียนรู้วิธีการประมวลผลภาพจากการรับรู้ระยะไกลโดยใช้การเชื่อมโยง Python GDAL และการถ่ายภาพจำนวนมาก เป็นครั้งแรกที่ฉันอ่านไฟล์ geotiff Landsat8 ทำการจัดการอย่างง่ายและเขียนผลลัพธ์ลงในไฟล์ใหม่ รหัสด้านล่างดูเหมือนว่าจะทำงานได้ดียกเว้นว่าแรสเตอร์ดั้งเดิมถูกทิ้งในไฟล์เอาต์พุตแทนที่จะเป็นแรสเตอร์ที่ได้รับการจัดการ ความคิดเห็นหรือข้อเสนอแนะใด ๆ ยินดีต้อนรับ แต่โดยเฉพาะอย่างยิ่งบันทึกเกี่ยวกับสาเหตุที่แรสเตอร์ที่จัดการไม่แสดงในผลลัพธ์ import os import gdal gdal.AllRegister() file = "c:\~\LC81980242015071LGN00.tiff" (fileRoot, fileExt) = os.path.splitext(file) outFileName = fileRoot + "_mod" + fileExt ds = gdal.Open(file) band = ds.GetRasterBand(1) arr = band.ReadAsArray() [cols, rows] = arr.shape arr_min = arr.Min() arr_max = …

2
การตั้งค่าหน่วยความจำในพื้นที่ทำงานใน ArcPy?
เป็นไปได้ไหมที่จะตั้งค่าenv.workspace = "in_memory"ใน ArcPy โดยใช้ทั้ง ArcGIS Pro และ ArcGIS 10.2.2 (หรือ 10.3) สิ่งที่ฉันพยายามทำคือรับเอาต์พุตของการดำเนินการsnap pour point ที่เขียนลงหน่วยความจำตรงข้ามกับดิสก์ ฉันรู้ว่าฉันสามารถเขียนลงดิสก์แล้วนำไปไว้ในหน่วยความจำได้ แต่ก็ไม่ได้ช่วยอะไร ฉันกำลังพยายามเพิ่มประสิทธิภาพชุดของงานลดการใช้การเขียนลงในสื่อทางกายภาพผลลัพธ์ของกระบวนการที่จำเป็นเพียงแค่เป็นอินพุตสำหรับกระบวนการถัดไป

1
ฉันจะเปลี่ยนชื่อโต๊ะใน SpatiaLite DB ได้อย่างไร
เลเยอร์พื้นที่ (ตาราง) ใน SpatiaLite ประกอบด้วยวัตถุและเมตาดาต้าที่รองรับต่าง ๆ รวมถึงทริกเกอร์ดัชนีและรายการใน (อย่างน้อยที่สุด) ในgeometry_columnsตาราง ฉันกำลังมองหา (ควร) GUI ที่จะจัดการกับการเปลี่ยนแปลงที่จำเป็นทั้งหมดในครั้งเดียวหรือ (ทางเลือก) เอกสารประกอบของการเปลี่ยนแปลงที่จำเป็นทั้งหมดเพื่อที่จะไม่ทำลายชั้นอวกาศ ตารางมี: geometry_columns.f_table_nameรายการใน ทริกเกอร์ชื่อห้าตัว[prefix]_[table_name]_geometryโดยที่คำนำหน้าระบุว่าเป็นการย้อนกลับธุรกรรมหรือการอัพเดตดัชนี ดัชนีอวกาศสามชื่อ idx_[table_name]_geometry[_suffix] ฉันได้ลองสิ่งนี้ในสองแอพพลิเคชั่นคือ QGIS DB Manager และ SpatiaLite-GUI QGIS DB Manager ปรากฏว่ามีผลกระทบต่อไปนี้ อัปเดตอย่างถูกต้องgeometry_columnsด้วยชื่อตารางใหม่ ห้ามเปลี่ยนชื่อทริกเกอร์ นิยามทริกเกอร์ที่มีการแก้ไขบางส่วนเพื่อให้BEFORE [INSERT|UPDATE|DELETE]หมายถึงชื่อตารางใหม่ geometry_columnsแต่ยังคงสภาพการค้นหาสำหรับชื่อเก่า ห้ามเปลี่ยนชื่อดัชนีเชิงพื้นที่ ไม่แน่ใจว่าสิ่งนี้สำคัญหรือไม่เพราะทริกเกอร์ยังคงอ้างถึงชื่อดัชนีเก่า เมื่อคุณเลือกการบำรุงรักษา→เปลี่ยนชื่อตารางใน SpatiaLite-GUI คุณเพิ่งได้สตับของALTER TABLEคำสั่งSQL นี่คือ SQL ที่ตรงไปตรงมาและทำได้น้อยกว่า QGIS DB Manager หากคุณกรอกชื่อตารางใหม่ตารางจะถูกเปลี่ยนชื่อ ผลกระทบอื่น ๆ …

5
จะดาวน์โหลดภาพ MOD16 (Modis Evapotranspiration) ได้อย่างไร?
ฉันใช้สคริปต์ 'R' เพื่อดาวน์โหลดตัวดัดแปลง MODIS อย่างไรก็ตามมันไม่ได้ให้ผลิตภัณฑ์ MOD16 ฉันไม่พบผลิตภัณฑ์ MOD16 ในเว็บไซต์ earthexplorer แม้ว่าผลิตภัณฑ์ที่ต้องการจะมีให้ที่http://www.ntsg.umt.edu/project/mod16#data-productแต่ฉันไม่พบสคริปต์ที่จะดาวน์โหลด มีสคริปต์ใด ๆ สำหรับมันหรือไม่? หรือฉันต้องคลิกแต่ละภาพเพื่อดาวน์โหลด ตัวอย่างเช่นหากฉันต้องการดาวน์โหลดภาพ " h02v08 " จากแต่ละเดือนและทุกปีจะทำอย่างไร?
11 qgis  r  remote-sensing  modis 

1
แสดงรูปภาพถัดจากคุณสมบัติใน ArcMap หรือไม่
ฉันมีคุณลักษณะหลายประการของจุดที่แสดงถึงป้ายจราจรและแต่ละจุดมีรูปถ่ายของแต่ละเครื่องหมายที่เกี่ยวข้อง ฉันต้องการสร้างแผนที่ด้วยตำแหน่งของสัญญาณที่แสดงเป็นจุดและภาพถ่ายที่แสดงถัดจากจุดนั้น แต่ละหน้าใน mapbook จะมีที่ตั้งและรูปภาพหลายแห่ง มีวิธีที่ง่ายกว่าในการนำเข้ารูปภาพไปยัง MXD โดยอัตโนมัติในมุมมองข้อมูลถัดจากคุณสมบัติจุดหรือไม่? มีประมาณ 1,000 จึงต้องการระบบอัตโนมัติ หรืออย่างใด "ติดป้าย" คุณสมบัติที่มีรูปถ่าย?

1
การเปลี่ยนตัวคั่นนับพันจากเครื่องหมายจุลภาคเป็นช่องว่างในกริดที่วัดได้ของ ArcMap?
ใน ArcMap เป็นการง่ายที่จะรวมเครื่องหมายจุลภาคเป็นตัวคั่นหลักพันบนกริดที่วัดได้ดังที่เห็นในภาพด้านล่าง แต่ฉันไม่สามารถหาวิธีเปลี่ยนตัวคั่นที่ใช้จากการเป็นเครื่องหมายจุลภาคเป็นช่องว่างแทนได้ ดูเหมือนจะไม่เป็นการตั้งค่าขั้นสูง ArcMap ที่จะทำและฉันได้ค้นหาและไม่พบแนวคิด ArcGIS สำหรับสิ่งนี้ มีวิธีการที่ฉันมองเห็นไหม?

1
อนาคตของ CartoCSS คืออะไร (เช่นมี) [ปิด]
ปิด คำถามนี้เป็นคำถามความคิดเห็นตาม ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้สามารถตอบข้อเท็จจริงและการอ้างอิงได้โดยแก้ไขโพสต์นี้ ปิดให้บริการใน4 ปีที่แล้ว MapBox เพิ่งประกาศว่าพวกเขาจะย้ายออกจาก CartoCSS ไปยังส่วนต่อประสานที่มองเห็นได้อย่างสมบูรณ์ เห็นได้ชัดว่าสไตล์ชีตใหม่ที่หนุน UI นั้นจะเป็นพื้นฐานของ JSON ดังนั้นแม้จะมีความจริงที่ว่าอินเทอร์เฟซใหม่เป็นระบบจุดและคลิกหนึ่งในนักพัฒนาที่ทำงานในโครงการนี้เชื่อว่าในท้ายที่สุดมันจะดีกว่า : เครื่องมือแก้ไขใหม่https://t.co/i7uy2V8YGLนั้นดียิ่งขึ้นสำหรับ coders: รูปแบบ GL สามารถปรับเปลี่ยนเป็น JSON ในตัวแก้ไขและแบบไดนามิกใน JS- Tom MacWright (@tmcw) 29 กรกฎาคม 2558 อย่างไรก็ตามมันก็คุ้มค่าที่จะสังเกต (และหวังว่าฉันมีสิทธิ์นี้) ว่าสไตล์ JSON เหล่านี้ถูกสร้างขึ้นโดย WebGL ในขณะที่ CartoCSS ใช้ในการสร้าง XML ที่ Mapnik สร้างการแสดงผล ดังนั้นในขณะที่หนึ่งแทนที่สำหรับอย่างอื่นอย่างน้อยในความสามารถบางอย่างพวกเขาไม่ได้ทำสิ่งเดียวกัน ดังนั้นสิ่งนี้ทำให้ฉันสงสัยว่าจะทิ้ง CartoCSS ไว้ที่ไหน หลายรูปแบบแผนที่กระเบื้องเด่นใช้ภาษาในเวลานี้อาจจะเป็นที่สะดุดตาที่สุดรูปแบบเริ่มต้นของ OpenStreetMap …
11 carto  mapbox  cartocss 

2
เหตุใดฉันจึงได้พื้นที่ที่ถูกต้องและพื้นที่ตัดกันเมื่อใช้การฉายภาพผิด
ฉันต้องคำนวณพื้นที่และพื้นที่สี่แยกสำหรับรูปหลายเหลี่ยม (วัตถุทางภูมิศาสตร์จริงบางอย่างเช่นทะเลสาบเมืองประเทศ ฯลฯ ) รูปหลายเหลี่ยมตั้งอยู่ในแคลิฟอร์เนียนิวซีแลนด์รัสเซียแอนนาดร์สวีเดน รูปหลายเหลี่ยมทั้งหมดอยู่ใน WGS84 สิ่งที่ฉันทำโดยใช้ GeoTool java API: ฉายรูปหลายเหลี่ยมทั้งหมดโดยใช้EPSG: 3488 , EPSG: NAD83 (NSRS2007) / California Albers และพื้นที่คำนวณและพื้นที่ทับซ้อน ทำเช่นเดียวกันโดยใช้World_MollweideและWorld_Eckert_IV เลือก "การฉายเฉพาะท้องถิ่น " สำหรับรูปหลายเหลี่ยมจากแคลิฟอร์เนียนิวซีแลนด์และอื่น ๆ ฉันคิดว่า # 3 เป็นผลลัพธ์ที่แม่นยำที่สุดเนื่องจากฉันเลือกการฉายที่ครอบคลุมพื้นที่รูปหลายเหลี่ยม ผลลัพธ์: '# 2 แสดงผลลัพธ์ที่แย่ที่สุดเมื่อเทียบกับ # 3 '# 1 และ # 3 พื้นที่และพื้นที่ทางแยกต่างกันน้อยกว่า 0.1% ทำไม? ฉันเลือกการคาดการณ์ที่ผิดอย่างผิดพลาด EPSG: 3488 (แคลิฟอร์เนีย) สำหรับรูปหลายเหลี่ยมจากสวีเดนและได้พื้นที่เดียวกันและพื้นที่สี่แยก …

3
วิธีการระบุชื่อไฟล์ gps?
ลองนึกภาพคุณมีการติดตาม GPS ที่คุณต้องการระบุตัวตนในแง่กฎหมาย คุณจะทำอย่างไร การจัดระยะห่างที่ใกล้ที่สุดและทำให้เวลาเพียงพอหรือไม่ มีมาตรฐานที่ตกลงกันในระดับสากลในเรื่องนี้หรือไม่? มีใครเขียนอัลกอริทึมเพื่อทำสิ่งนี้อยู่แล้ว? ถ้าไม่ใช่ฉันวางแผนที่จะทำฟังก์ชั่นในแพ็คเกจ stplanr ที่กำลังพัฒนาของฉันเพื่อทำสิ่งนี้ ตัวอย่างที่ทำซ้ำได้ (ใช้ฟังก์ชั่นการหมุนที่ยอดเยี่ยมจาก @geospacedman) จากข้อมูล 'ตัวตนของฉัน': library(rgdal) library(tmap) downloader::download("https://www.openstreetmap.org/trace/1619756/data", "test.gpx") r <-readOGR(dsn = "test.gpx", layer = "tracks") r <- spTransform(r, CRS("+init=epsg:27700")) rproj <- rotateProj(rs, 90) # rotate projection for plotting r <- spTransform(r, rproj) rs <- rgeos::gSimplify(r, 1000) # snap to …
11 r  gps  gpx  simplify 

1
วิธีเพิ่มความเร็วในการทำงานพร้อมกันกับชั้น PostGIS ใน QGIS?
เรากำลังทำงานบนเลเยอร์ใหญ่ในเวลาเดียวกัน (10-20 คน) เราบางคนใช้ QGIS 2.8.1 และ 2.6 อื่น ๆ ผู้ใช้ที่มี 2.8 มีข้อผิดพลาดจำนวนมาก (QGIS) ผู้ใช้ 2.6 คนไม่มีข้อขัดข้องนี้ แต่การเพิ่มคุณสมบัตินั้นช้าอย่างไม่น่าเชื่อ เซิร์ฟเวอร์ Postgres 9.4 ทำงานบนเครื่องเสมือน (Hyper V) ด้วย Win7 อาจเป็นเครือข่ายไม่ใช่ QGIS ใช่ไหม บางทีใครบางคนอาจมีความคิด จะดี :)

3
ลูกเรือนำทางอย่างไรก่อนการฉายของ Mercator (1536) [ปิด]
ปิด. คำถามนี้เป็นคำถามปิดหัวข้อ ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้เป็นไปตามหัวข้อสำหรับระบบแลกเปลี่ยนข้อมูลทางภูมิศาสตร์ของ Stack ปิดให้บริการใน2 ปีที่ผ่านมา การฉายภาพ Mercator (1536) จะต้องยอดเยี่ยมสำหรับการนำทางทะเล ช่วยให้ผู้ใช้สามารถค้นหาตลับลูกปืนและทำให้ค่าคงที่นี้ควรถึงเป้าหมาย (หากไม่อยู่ในเส้นที่สั้นที่สุด) แต่คนนำทางก่อนหน้านี้ได้อย่างไร ในฐานะที่เป็นตัวอย่างมาเจลลันโคจรรอบโลกแคลิฟอร์เนีย 1520

3
ไฟล์รูปร่างเมืองฟรี
ฉันกำลังมองหาไฟล์รูปร่างจุดที่เรียบง่ายซึ่งระบุว่าเมืองของโลกอยู่ที่ไหน ฉันดูเหมือนจะไม่พบที่ไม่มีค่าใช้จ่ายและตรงกับความต้องการของฉัน ฉันมีสิ่งนี้: http://techcenter.jefferson.kctcs.edu/data/แต่ฉันต้องการเมืองเล็ก ๆ ด้วย ... ความคิดใดที่ฉันสามารถรับได้


2
จุด GPS ที่จับคู่แผนที่ไปยังเครือข่ายถนน!
ฉันค่อนข้างใหม่สำหรับ Postgres และ PostGIS ฉันกำลังพยายามทำการจับคู่แผนที่อย่างง่าย ๆ (หรือการอนุมานเส้นทางถ้าคุณต้องการ) ของบางจุด GPS ฉันใช้postgres (PostGIS) ฉันมีตารางในฐานข้อมูลที่มีคะแนน GPS ของฉัน: gpspoints (Lat, Long, time, speed, azimuth, geometry, ... ) ฉันได้นำเข้าไฟล์รูปร่างของเครือข่ายถนนลงในฐานข้อมูลของฉันด้วย: เส้นทาง (gid, idrte, รุ่น, nomrte, norte, clsrte, geom, ... ) - ฉันไม่มี azimuth ของลิงก์และฉันเป็น ไม่แน่ใจว่าฉันจะคำนวณได้อย่างไรโดยใช้ฟังก์ชั่น postgis ST_azimuth ทั้งหมดที่ฉันมีคือ shapefile ของเครือข่ายถนนที่มีคอลัมน์ที่ระบุข้างต้น ฉันต้องการเชื่อมโยงจุด GPS แต่ละจุดกับลิงก์ที่ใกล้ที่สุด (ในบัฟเฟอร์ประมาณ 20 เมตรรอบจุด) …

4
บันทึกมุมมองนักแต่งเพลงพิมพ์ / แผนที่ QGIS เป็น PNG / PDF โดยใช้ Python (โดยไม่ต้องเปลี่ยนอะไรในเค้าโครงที่มองเห็นได้)?
ฉันเปิดมุมมองผู้แต่ง / พิมพ์ QGIS พร้อมองค์ประกอบทั้งหมดที่ปรับตามที่ฉันต้องการ ตอนนี้ฉันต้องพิมพ์ / บันทึก / ส่งออกเป็นไฟล์ PNG หรือ PDF จาก Python console / Python script ฉันไม่ต้องการเปลี่ยนแปลงอะไรในเค้าโครงปัจจุบัน ตัวอย่างส่วนใหญ่ที่ฉันพบ (ตัวอย่างเช่น: สิ่งนี้ ) เปลี่ยนตำแหน่งบนแผนที่หรือขนาดในเอาต์พุต PDF เปรียบเทียบกับสิ่งที่ฉันเห็นในมุมมองนักแต่งเพลงปัจจุบันของฉัน ฉันต้องการที่จะได้รับผลตรงเช่นเดียวกับฉันจะได้รับเมื่อฉันคลิกที่พิมพ์ -> ส่งออกเป็นภาพ ฉันจะทำสิ่งนั้นได้อย่างไร Atlas ไม่ใช่ทางออกสำหรับฉัน

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.