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

ในแพ็คเกจ HistData สำหรับ R ( https://r-forge.r-project.org/R/?group_id=574 ) ฉันมีชุดข้อมูลที่เกี่ยวข้องกับแผนที่ของ John Snow เกี่ยวกับการระบาดของอหิวาตกโรคในลอนดอนปี 1854 ฉันเชื่อว่าพวกเขา มีอำนาจมีการแปลงเป็นดิจิทัลอย่างระมัดระวังภายใต้การดูแลของวอลเตอร์ Tobler รายละเอียดบางอย่างเกี่ยวกับชุดข้อมูลเหล่านี้จะอธิบายโดยจอห์นแม็คเคนซี่ที่http://www1.udel.edu/johnmack/frec480/cholera/cholera2.html น่าเสียดายที่พิกัดแห่งความตายปั๊มและถนนใช้ระบบพิกัดโดยพลการไม่ใช่พิกัดแผนที่ที่เหมาะสมสำหรับแอปพลิเคชั่น GIS อื่น ๆ หรือซอฟต์แวร์แผนที่ใน R (แพ็คเกจอวกาศ ggmap ฯลฯ ) ในhttp://freakonometrics.hypotheses.org/19213อาร์เธอร์ Charpentier ใช้ ggmap กับรุ่นของข้อมูลจอห์นหิมะจาก http://www.rtwilson.com/downloads/SnowGIS_v2.zip Cholera_Deaths.shpไฟล์ แต่แสดงรายการเพียง 489 ตายไม่ได้ 578 HistData::Snow.deathsผมได้บันทึกไว้ใน แนวคิดหนึ่งคือการค้นหาความสัมพันธ์ระหว่างค่าเฉลี่ยและค่าเบี่ยงเบนมาตรฐานของพิกัด (x, y) และ rescale เป็นเส้นตรง แต่อาจจะมีวิธีที่ดีกว่านี้หรือ นี่คือสิ่งที่ฉันได้ลองไปแล้ว > data(Snow.deaths, package="HistData") > …

14 r  coordinate-system  scale 

ใน QGIS 2.14.5 LTR (และเมื่อเร็ว ๆ นี้ใน QGIS 3.0) ฉันได้ตั้งค่าการออกแบบแผนที่ที่สวยงามด้วยตัวแสดงภาพที่สำเร็จการศึกษาและตัวบ่งชี้เขตเวกเตอร์สำหรับแต่ละช่วงค่า (ความยาวลูกศรเป็นสัดส่วนโดยตรงกับความเร็วกระแสซึ่งอยู่ในช่วง 0.001 - ~ 5m / s) สิ่งนี้ใช้งานได้ดี แต่ฉันสงสัยว่าฉันพลาดการตั้งค่าเพื่อทำให้ลูกศรในตำนานสั้นลงหรือไม่ ด้วยหน้าต่างชั้นความกว้าง 'ปกติ' ค่าคลาสไม่สามารถมองเห็นได้: หน่วย: มิลลิเมตร เดาแรกของฉันคือว่ามันมีฏ จะทำอย่างไรกับการสลับความยาวลูกศรหน่วยความกว้างของเส้นจากหน่วยมิลลิเมตรไปยังหน่วยแผนที่ แต่สิ่งนี้จะไม่ส่งผลต่อตำนาน: หน่วย: หน่วยแผนที่ UPDATE ปัญหายังคงอยู่ใน QGIS 3.0 ฉันสร้าง screencast เพื่อให้เห็นภาพนี้: https://youtu.be/oxJ3tGXFB18 ในคุณสมบัติ / สัญลักษณ์สัญลักษณ์แสดงผลอย่างถูกต้อง: โปรดทราบว่ารายการเลเยอร์แสดงเครื่องหมายอย่างถูกต้องเมื่อทำซ้ำเลเยอร์หรือใช้การเปลี่ยนแปลงกับสไตล์ตราบใดที่ความกว้างของเลเยอร์แผงไม่เปลี่ยนแปลง

14 qgis  legend  markers  qgis-3  vector-field 

ฉันพยายามถ่ายสองบรรทัดโดยใช้ Shapely / Geopandas แต่ผลลัพธ์ของการจัดเรียงนั้นแปลกมาก ฉันเหนื่อย : import geopandas as gpd from shapely.geometry import * from shapely.ops import snap lines1 = gpd.GeoDataFrame.from_file('lines1.shp') lines1 = lines1.to_crs({'init': 'epsg:2227'}) lines2 = gpd.GeoDataFrame.from_file('lines2.shp') lines2 = lines2.to_crs({'init': 'epsg:2227'}) res = lines1 lines2_union = lines2.geometry.unary_union res.geometry = res.geometry.apply(lambda x: snap(x, lines2_union, 14)) res.to_file('result.shp', driver="ESRI Shapefile") และได้รับผลลัพธ์นี้: …

14 shapely  snapping  geopandas 

ฉันใหม่ในฟิลด์การรับรู้ระยะไกล บางครั้งฉันสับสนเมื่อฉันอ่านเกี่ยวกับการประมวลผลภาพล่วงหน้า บางคนช่วยอธิบายรายละเอียดเกี่ยวกับความแตกต่างและตัวอย่างวิธีการแก้ไขและการสอบเทียบด้วยคลื่นวิทยุได้ไหม? การแปลงค่า DN เป็นค่าการสะท้อนกลับสูงสุดของบรรยากาศ (TOA) จากภาพ 2 ชุด (เช่น Landsat 7) จากเวลาต่างกันนับเป็นการแก้ไขหรือการปรับเทียบหรือไม่? DOS (การลบวัตถุในความมืด) และส่วนการแก้ไขอาทิตย์ของการแก้ไขด้วยรังสีหรือไม่?

14 remote-sensing 

ในการรับค่าพิกเซลเพียงจุดเดียวใน raster โดยใช้ rasterio มีตัวอย่างอยู่ที่นี่: https://github.com/mapbox/rasterio/pull/275 อย่างไรก็ตามมี API โดยตรงภายใน rasterio (ไม่ใช่ cli) ซึ่งสามารถใช้เพื่อแยกค่าที่จุดเดียวใน raster หรือไม่? - แก้ไข with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') as src: x = (src.bounds.left + src.bounds.right) / 2.0 y = (src.bounds.bottom + src.bounds.top) / 2.0 vals = src.sample((x, y)) …

14 python  numpy  rasterio 

มีวิธีการรีเฟรชตารางคุณลักษณะของเลเยอร์หลังจากแก้ไขตาราง (เช่นเพิ่มคอลัมน์) ใน DB-Manager หรือไม่?

14 qgis  attribute-table  database 

ฉันมีแอปพลิเคชั่นที่ต้องทำการร้องขอ HTTP (ซึ่งจะจัดการการแปลงแผนที่) ก่อนดึงไฟล์ WMS จากเซิร์ฟเวอร์ มีวิธีหยุด TileWMS รีเฟรชอัตโนมัติหลังจากลากแผนที่หรือไม่ การตั้งค่าที่มองไม่เห็นหรือลบและเพิ่มเลเยอร์จะใช้งานได้ แต่ฉันไม่ต้องการให้เลเยอร์หายไปในขณะที่ดำเนินการตามคำขอ HTTP ฉันพยายามแล้วsetUrl('')แต่เลเยอร์ก็จะหายไปเช่นกัน ฉันพยายามอ่านซอร์สโค้ดของ OpenLayers เพื่อกำหนดคุณสมบัติที่จะตรวจสอบก่อนดำเนินการตามคำขอ แต่ฉันไม่พบรหัสที่ร้องขอ

14 wms  openlayers 

ดูเหมือนว่ามันจะช่วยได้มากในการจัดเก็บภาพให้มากที่สุดเท่าที่เป็นไปได้บน GPU สำหรับการแพนและซูม ฯลฯ เราจัดการกับภาพที่มีขนาดใหญ่มาก มีให้เลือกหลากหลายตั้งแต่ 8 GB ถึงหลายร้อย GB ที่บ้านฉันมีจอภาพขนาด 40 "4k และ GTX Titan ที่มีหน่วยความจำวิดีโอ 12 GB และมันทำงานได้ดีอย่างน่าอัศจรรย์สำหรับทุกสิ่งที่ฉันขว้างไปมันดูดีมากเมื่อดูภาพถ่าย งาน? หากว่าQGIS , ArcMapและPCI Geomaticaจะสามารถใช้ประโยชน์จากหน่วยความจำกราฟิกมากมายเมื่อจัดการกับภาพหลายกิกะไบต์และโมเสคภาพในช่วงหลายร้อยกิกะไบต์หรือไม่ แอปพลิเคชัน GIS ใดที่เห็นว่าได้รับประโยชน์มากที่สุดจากหน่วยความจำ GPU มากมายเมื่อจัดการกับภาพขนาดใหญ่

14 qgis  arcmap  performance  hardware 

ฉันมีแอพ django ที่มีแผ่นพับที่ส่วนหน้าและฉันจำเป็นต้องคำนวณ bbox ของวัตถุ geojson เพื่อให้ฉันสามารถส่งพิกัดไปยัง map.fitBounds() ฉันลองใช้ libs แล้วลองใช้getBoundsวิธีการของวัตถุบางอย่างใน Leaflet เช่นFeatureGroupแต่มันก็บ่นว่าขอบเขตไม่ได้ถูกกำหนดไว้ ใครช่วยชี้ให้ฉันเห็นวิธีแก้ปัญหาง่ายๆ

14 python  leaflet  geojson  django 

ฉันต้องรวม DEM rasters หลาย ๆ อันไว้ในแหล่งนี้: http://srtm.csi.cgiar.org/SELECTION/inputCoord.aspดังนั้นฉันจึงใช้เครื่องมือ Mosaic กับ New Raster ฉันตั้งค่าเครื่องมือ: ฉันป้อน rasters (ทั้งหมดจากแหล่งเดียวกันขนาดเดียวกันและไม่มีการฉายภาพ); ตั้งค่า Pixel Type (จุดลอยตัว 32 บิตเหมือนในภาพแรสเตอร์ดั้งเดิม); ตั้งค่าขนาดเซลล์เป็นแบบแรสเตอร์ดั้งเดิม จำนวนวง = 1 เป็น rasters ดั้งเดิม ผู้ประกอบการโมเสกฉันได้ทำมันด้วย BLEND และ MEAN (ฉันได้ผลลัพธ์เดียวกัน) ปัญหาที่ฉันมีคือโมเสคผลลัพธ์แสดงช่วงของค่าสูงสุดและต่ำสุดที่แตกต่างจากค่าสูงสุดและต่ำสุดของแต่ละแรสเตอร์เช่น raster 1 (-5123.8, 23.25), raster 2 (-5974.6, 40.09), แรสเตอร์ 3 (-57770.2, 38), แรสเตอร์ 4 (-2534.3, 23.55) …

14 arcgis-desktop  raster  arcgis-10.2  dem  mosaic 

ฉันต้องเพิ่มรายการเมนูใหม่ลงในเมนูระดับบนสุดของ QGIS Desktop (2.7.0); นั่นคือจะต้องอยู่ในระดับเดียวกับ "Project", "Edit", "View", "Layer" เป็นต้น เท่าที่ฉันรู้มีวิธีที่ดีในการซ่อนรายการเมนู (การตั้งค่า> การปรับแต่ง) และการใช้ QgisInterface ( http://qgis.org/api/classQgisInterface.html ) ฉันสามารถเพิ่มรายการใหม่ลงในเมนู QGIS ได้แล้ว ( addLayerMenu, addPluginToDatabaseMenu ฯลฯ ) จากรหัสปลั๊กอินของฉัน แต่ฉันต้องการรายการเมนูใหม่ที่เลเยอร์บนสุด (นี่คือคำขอจากลูกค้า)

14 qgis  pyqgis  qt  gui  context-menu 

ฉันเป็นผู้เริ่มต้นที่สมบูรณ์ของข้อมูลทางภูมิศาสตร์ดังนั้นโปรดยกโทษให้ฉันหากคำถามไม่เหมาะสม ฉันดาวน์โหลดข้อมูลจาก NCDC NARR และจัดการโหลดลงใน R โดยใช้rasterแพ็คเกจ ฉันต้องการรับรายการที่มีละติจูดลองจิจูดและค่า ฉันเข้าใจว่าrasterToPoints()ควรทำสิ่งที่ฉันต้องการอย่างแน่นอนอย่างไรก็ตามค่าละติจูดและลองจิจูดของฉันดูแปลก ๆ : r <- raster(myfile) data_matrix <- rasterToPoints(r) head(data_matrix) x y value [1,] -5405401 4347242 70 [2,] -5372938 4347242 88 [3,] -5340475 4347242 76 [4,] -5308012 4347242 85 [5,] -5275549 4347242 87 [6,] -5243086 4347242 88 ฉันคิดว่าฉันควรทำอะไรกับการฉายภาพซึ่งปัจจุบันเป็น Lambert Conformal Conic (LCC) …

14 raster  r 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน4 ปีที่แล้ว ฉันมีค่าใช้จ่ายจำนวนมากที่ขอให้รวมปลั๊กอินของ Google Earth เข้ากับเว็บไซต์ของเขา จากไซต์ปลั๊กอิน Google Earth ฉันอ่านสิ่งนี้: API ปลั๊กอินของ Google Earth เลิกใช้แล้วเมื่อวันที่ 12 ธันวาคม 2014 API จะยังคงทำงานกับเบราว์เซอร์ที่รองรับจนถึงวันที่ 12 ธันวาคม 2015 และจะปิดตัวลงในวันนั้น มีทางเลือกอะไรบ้าง? ฉันรู้ว่ามี Google Maps API, Google Maps Engine API, Google Street View API แต่ฉันสับสนเกี่ยวกับพวกเขาและที่หนึ่งสามารถ "ทดแทน" ปลั๊กอิน Google Earth คุณสมบัติที่ฉันต้องการคือสิ่งสำคัญตามลำดับ: โหลดไฟล์ KMZ …

14 google-earth  google-maps-api  api  software-recommendations  google-maps-engine 

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

14 geoprocessing  fme  clip 

ฉันมีตารางรูปหลายเหลี่ยมแสดงเวลาในการเดินทางในแต่ละวัน สำหรับจุดเริ่มต้นแต่ละจุดมีรูปทรงเรขาคณิต isochrone ห้า (เก็บไว้ในแถวที่แยกต่างหาก) สำหรับจุดเริ่มต้นแต่ละจุดฉันต้องการที่จะรวมห้า isochrones (ไบนารี NULL หรือ 1) แล้วรวมเข้ากับชั้นแรสเตอร์เดียว เลเยอร์แรสเตอร์นี้ต้องใช้พีชคณิตแผนที่อย่างง่าย ๆ : sum / 5 เพื่อให้ต้นกำเนิดแต่ละอันจะเชื่อมโยงกับเลเยอร์แรสเตอร์เดี่ยวที่มีค่าใน [NULL, 0.2, 0.4, 0.6, 0.8, 1] ขึ้นอยู่กับจำนวน ชั้นที่เป็นองค์ประกอบทับซ้อนกัน มันเป็นพื้นผิวที่น่าจะเป็น ข้อมูลของฉันถูกเก็บไว้ใน Postgres 9.3 (พร้อม PostGIS) ปัญหาของฉันคือในขณะที่ฉันต้องการเรียนรู้การใช้ PostGIS แรสเตอร์ดูเหมือนว่าจะมีช่วงการเรียนรู้ที่สูงชันและตัวอย่างทั้งหมดที่ฉันสามารถหาข้อตกลงกับชั้นแรสเตอร์เดียว ในตัวอย่างเลเยอร์นี้ใช้เป็นส่วนหนึ่งของรูปหลายเหลี่ยมซ้อนทับอาจเฉลี่ยค่าของแรสเตอร์สำหรับรูปหลายเหลี่ยมแต่ละรูป ฉันไม่พบตัวอย่างที่จำลองได้สำหรับการรวม: a) เวกเตอร์ -> raster b) พีชคณิตแผนที่ และ c) GROUP BY attribute ตามวรรคหนึ่งของฉัน …

14 postgis  raster  map-algebra  postgis-raster