คำถามติดแท็ก raster

สั้น: ฉันต้องการเปลี่ยนความละเอียดของแรสเตอร์และปรับทางลาดสีเทาตามที่แสดงในภาพตะโกน การตั้งค่าคือการใช้ GDAL, PIL หรือ Numpy คำอธิบาย: ฉัน Kriging จุดลงไปในเชิงภาพที่มีความละเอียดการส่งออกของ 20 เมตรกับGeostatistical ห้องสมุดที่มีประสิทธิภาพสูง ฉันไม่ต้องการเปลี่ยนความละเอียดเอาต์พุตเนื่องจากเวลาการแก้ไขเพิ่มขึ้นเป็นทวีคูณ ด้วยความละเอียดนี้ภาพที่ส่งออกจะน่าเกลียด (pixelated และ aliased) ฉันไม่รู้ว่ามันถูกต้องในแนวความคิดหรือไม่ แต่ฉันต้องการภาพที่นุ่มนวลเหมือนในตัวอย่างร้อง มันเหมือนกับการ 'ปรับกำลังภาพ' ใหม่ให้เป็นภาพที่มีความละเอียดมากกว่า ฉันใช้ python ดังนั้นการตั้งค่าของฉันคือ GDAL, Python Imaging Library หรือ Numpy คำตอบอาจเป็นเชิงทฤษฎีเช่นชี้ให้เห็นชื่ออัลกอริทึมหรือแนวคิดของการดำเนินการชนิดนี้ ที่มา: ปลายทาง: แก้ไขผลลัพธ์ด้วย gdalwarp cubic spline:

15 python  raster  gdal 

ฉันมีการจำแนกประเภทที่ดินครอบคลุมซึ่งได้มาจากภาพถ่าย Landsat TM - ฉันสร้างการจำแนกประเภทนี้โดยใช้ ENVI ฉันต้องการแก้ไขรายบุคคลและกลุ่มของเซลล์ในบางพื้นที่ของการจัดหมวดหมู่เพื่อปรับปรุงความแม่นยำของผลลัพธ์โดยการเปรียบเทียบกับภาพถ่ายทางอากาศ มีวิธีใดบ้างที่ใช้ ArcGIS Desktop ในการแก้ไขกลุ่มเซลล์ขนาดใหญ่ด้วยตนเองโดยการเลือกและกำหนดให้คลาสใหม่หรือคลาสที่มีอยู่? ฉันทราบวิธีการรวมตามสถิติของชั้นเรียน แต่พื้นที่ที่ฉันจัดประเภทมีความสับสนทางสเปกตรัมในระดับสูงเนื่องจากมีเมฆปกคลุมบาง - ดูเหมือนว่าการแก้ไขด้วยตนเองเป็นวิธีเดียวที่จะไป

15 arcgis-desktop  raster  arcgis-10.0  editing 

ฉันใช้PostGIS2.0ทำทางแยก raster / polygon ฉันมีปัญหาในการทำความเข้าใจการทำงานที่ฉันควรใช้และวิธีที่เร็วที่สุดในการแสดงคืออะไร ปัญหาของฉันเป็นดังนี้: ฉันมีรูปหลายเหลี่ยมและแรสเตอร์ ฉันต้องการค้นหาพิกเซลทั้งหมดที่อยู่ในรูปหลายเหลี่ยมและได้รับผลรวมของค่าพิกเซล และ (ปัญหาที่อัปเดต): ฉันได้รับค่าจำนวนมากสำหรับบางพิกเซลที่ไม่มีอยู่ในแรสเตอร์ดั้งเดิมเมื่อฉันทำการค้นหา ฉันมีความเข้าใจความยากลำบากไม่ว่าฉันควรใช้หรือST_Intersects() ST_Intersection()ฉันยังไม่รู้ด้วยว่าวิธีที่ดีที่สุดในการรวมพิกเซลของฉันคืออะไร นี่เป็นวิธีแรกที่ฉันได้ลอง (# 1): SELECT r.rast FROM raster as r, polygon as p WHERE ST_Intersects(r.rast, p.geom) ส่งคืนรายการของrastค่าซึ่งฉันไม่แน่ใจว่าจะทำอย่างไร ฉันพยายามคำนวณสถิติสรุปโดยใช้ST_SummaryStats()แต่ฉันไม่แน่ใจว่านี่คือผลรวมถ่วงน้ำหนักของพิกเซลทั้งหมดที่อยู่ในรูปหลายเหลี่ยม SELECT (result).count, (result).sum FROM ( SELECT ST_SummaryStats(r.rast) As result FROM raster As r, polygon As p WHERE ST_Intersects(r.rast, p.geom) ) …

15 raster  polygon  sql  spatial-database  postgis-2.0 

ฉันพยายามแปลง Shapefile เป็นแรสเตอร์ภายใน R วิธีการของฉันคือการอ่านในแรสเตอร์ดังต่อไปนี้: library(rgdal) # Loads SP package by default demo <- readOGR('F:/data/', 'shapefile') # Creates a SpatialPolygonsDataFrame class (sp) มันใช้งานได้ดีและฉันสามารถลงจุดได้ อย่างไรก็ตามมันเป็นไฟล์รูปร่างขนาดใหญ่และฉันต้องการแปลงให้เป็นแรสเตอร์ ฉันได้ลองทำสิ่งต่อไปนี้แล้ว: r <- raster(ncol=180, nrow=180) Demo_ras = rasterize(r, demo, 'pop') # pop is an integer here ข้อความผิดพลาด: Error in function (classes, fdef, mtable) : unable to …

15 raster  shapefile  r 

นี่คือรูปสวย ๆ ของแรสเตอร์ที่ฉันใช้ - โดยไม่มีเอฟเฟกต์ใด ๆ - เพียงแค่ปรับสีและฮิสโตแกรมเพื่อแสดงความโล่งใจ - แม้ว่าจะเป็นแบบ 2 มิติ! ฉันมีเลเยอร์แรสเตอร์ที่ฉันต้องการแสดงพร้อม Hillshade เป็นค่าเริ่มต้น ฉันสามารถทำเช่นนี้ในรหัสและผลก็เหมือนกับการเลือกตัวเลือก Hillshade ใน ArcMap (ดูด้านล่าง): เมื่อคุณทำเช่นนี้ - มันดูน่าผิดหวัง: โดยเฉพาะอย่างยิ่งถ้าคุณซูมเข้าไปในแรสเตอร์คุณจะเห็นสิ่งประดิษฐ์ที่ gridded ทุกชนิด - ซึ่งดูไม่ดีเลย! อย่างไรก็ตามถ้าฉันสร้างเลเยอร์ Hillshade ที่เหมาะสมโดยใช้แรสเตอร์ด้วย Spatial Analyst และวางทับด้วยความใส 85% มันดูดีและแตกต่างกันมาก!: ฉันจะได้รับผลกระทบนี้ใน ArcMap และแทนที่การทำงานเริ่มต้นของตัวเลือก Hillshade ได้อย่างไร (ฉันคิดว่ามันเป็นไปไม่ได้หรือมันเป็นการตั้งค่าที่ฉันพลาดไปซึ่งทำให้ raster แสดงได้ไม่ดี) ถ้ามันเป็นไปได้ที่จะแทนที่หรือเลือกการตั้งค่าที่ถูกต้อง - เป็นไปได้ที่จะทำใน ArcObjects หากตัวเลือกที่ 1 …

15 arcobjects  arcmap  raster  hillshade 

เวอร์ชั่น ArcGIS ล่าสุด (10) มีเครื่องมือ arctoolbox สองชุดสำหรับตรรกะคลุมเครือ ( Fuzzy MembershipและFuzzy Overlay ) ไม่มีใครรู้ว่าฉันสามารถหาตัวอย่างที่แสดงวิธีการที่พวกเขาสามารถนำไปใช้ที่ไหน? (จากวิธีใช้เดสก์ท็อป 10 )

15 arcgis-10.0  raster  spatial-analyst 

ฉันนิ่งงันและการค้นหาจำนวนมากไม่ได้มีประโยชน์อะไรเลย ฉันอยู่ในขั้นตอนแรกของสิ่งที่ในที่สุดจะเป็นแผนที่ที่มีสามชั้นแรสเตอร์และสามชั้นของเครื่องหมาย (พร้อมการควบคุมเพื่อสลับระหว่าง rasters และข้อมูลที่เกี่ยวข้องกับป๊อปอัป / เคล็ดลับเครื่องมือสำหรับข้อมูล) แต่ฉันไม่สามารถ รับไทล์มิลล์เพื่อแสดงอย่างถูกต้อง (ฉันสร้าง rasters ใน QGIS วิ่งผ่าน GDAL เพื่อทำสีและฉายมันอีกครั้งแล้วนำไปใส่ใน TileMill) หากต้องการปัญญา: นั่นคือใน Chrome นี่คือ Safari: และนี่คือรหัสซึ่งดูไม่น่ากลัวสำหรับฉัน (ฉันได้แสดงความคิดเห็นทุกอย่างยกเว้นกระเบื้อง ณ จุดนี้ดังนั้นนั่นคือทั้งหมดที่ฉันได้รวมไว้): <html> <head> <title> - Leaflet demo</title> <script type='text/javascript' src="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.js"></script> <style type='text/css'> html, body{ width:100%; height: 100%, margin:0; padding:0; } #map{ width:700px; height:600px } img …

15 raster  javascript  leaflet  tilemill  browser 

ฉันมีรหัส Python ที่ทำให้ Raster Attribute Tables สำหรับ GeoTiffs โดยใช้ GDAL ใน ArcGIS ฉันสามารถเห็นตารางคุณลักษณะเหล่านี้ได้ดีและพวกเขาประพฤติตัวตามที่ฉันชอบ (ในการจำแนกประเภท) แต่ใน QGIS แนวคิดทั้งหมดของตารางแอตทริบิวต์ดูเหมือนจะหายไป มีวิธีการเข้าถึง RAT ใน QGIS หรือไม่?

15 qgis  raster  attribute-table 

ฉันต้องรวม 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 

ฉันเป็นผู้เริ่มต้นที่สมบูรณ์ของข้อมูลทางภูมิศาสตร์ดังนั้นโปรดยกโทษให้ฉันหากคำถามไม่เหมาะสม ฉันดาวน์โหลดข้อมูลจาก 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 

ฉันมีตารางรูปหลายเหลี่ยมแสดงเวลาในการเดินทางในแต่ละวัน สำหรับจุดเริ่มต้นแต่ละจุดมีรูปทรงเรขาคณิต 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 

ฉันมีไฟล์รูปร่างเส้นเป็นตัวแทนเครือข่ายถนน ฉันต้องการ rasterize ข้อมูลนี้ด้วยค่าผลลัพธ์ในแรสเตอร์ที่แสดงความยาวรวมของบรรทัดที่อยู่ในเซลล์แรสเตอร์ ข้อมูลอยู่ในการประมาณการกริดแห่งชาติของอังกฤษดังนั้นหน่วยจะเป็นเมตร เป็นการดีที่ฉันต้องการใช้การดำเนินการนี้Rและฉันคาดเดาว่าrasterizeฟังก์ชั่นจากrasterแพคเกจจะมีส่วนร่วมในการบรรลุเป้าหมายนี้ฉันแค่ไม่สามารถทำงานออกได้ว่าควรจะใช้ฟังก์ชั่นใด

14 raster  r  line  rasterization 

ฉันมีภาพแรสเตอร์ (GeoTIFF) ที่มีค่าจำนวนเต็ม ฉันรู้วิธีการมองเห็นค่าแรสเตอร์ใน QGIS โดยใช้ colormaps หรือคล้ายกัน แต่ฉันต้องการกำหนดสัญลักษณ์ให้กับค่าสีแต่ละค่าแทน - เช่นเดียวกับ colormap แต่มีสัญลักษณ์ ตัวอย่าง: คำถาม: เป็นไปได้หรือไม่และในกรณีนี้เป็นอย่างไร?

14 qgis  raster  symbology  visualisation 

ฉันมี DEM ขนาด 1x1 เมตรและ LiDAR point cloud ต้นฉบับเป็น. Atlas นั่นคือ DEM ฉันต้องการแยกแม่น้ำโขง (จุดสูงสุดของเขื่อน) ไปยังคุณลักษณะของเวกเตอร์ (จุด, เส้นหลายเหลี่ยม) มีความคิดเกี่ยวกับอัลกอริทึมหรือเครื่องมือที่มีอยู่หรือไม่? เขื่อนกั้นน้ำมีสีเทาอ่อนในภาพที่ 1 และด้านล่างเป็นภาพเมฆจุดสำหรับพื้นที่ตัวอย่างที่มีเขื่อนกั้นน้ำ แม่น้ำตรงกลางเป็นสีน้ำเงิน เครื่องมือใดสำหรับการสกัดคุณสมบัติเชิงเส้นจาก DEM? ฉันเชื่อว่าฉันต้องการบางอย่างเช่นโมดูลการดึงข้อมูลเชิงพื้นที่ของ ENVI (หน้า 7) แต่ไม่มีส่วนของ ENVI เนื่องจากฉันไม่สามารถจ่ายได้ :)

14 raster  dem  lidar  algorithm  feature-extraction 

ฉันใหม่มากที่ใช้ข้อมูล GIS และมีประสบการณ์เพียงเล็กน้อยกับ R. ฉันได้อ่านเกี่ยวกับวิธีการวิเคราะห์ข้อมูลเชิงพื้นที่โดยใช้หนังสือ PDF spatial-analyst.net ดังนั้นฉันจึงไม่หลงทาง แต่ฉันคิดว่าฉันสามารถอธิบายได้ ปัญหาของฉันและผู้คนอาจแนะนำแนวคิด ฉันมีชุดข้อมูลที่มีการวัดประมาณ 2,000 ชุดที่พิกัดละติจูด / ลองจิจูดที่แตกต่างกันถึงแม้ว่าฉันอาจจะแบ่งชุดข้อมูลนี้ออกเนื่องจากมีการรวบรวมข้อมูลมากกว่า 3 ปีและมีการเปลี่ยนแปลงเงื่อนไขเมื่อเวลาผ่านไป ลองเรียกตัวแปรที่วัดว่า "IP" ฉันต้องการสร้างแผนที่ IP ในพื้นที่เต็มรูปแบบภายใต้คำถามโดยใช้ Kriging หรือวิธีการแก้ไขอื่น ๆ ในข้อมูลตัวอย่าง จากนั้นฉันต้องการสร้างฮิสโตแกรมที่วัดปริมาณที่ดินในถังเก็บไอพีแบบต่างๆ ฉันจะต้องสร้างฮิสโตแกรมซึ่งแสดงจำนวนตัวอย่างในถังแต่ละอัน (หมายเหตุตัวอย่างอาจมี IP จริงที่สูงกว่าหรือต่ำกว่าสิ่งที่ kriging คาดการณ์ไว้สำหรับที่ดินของตน) ฉันติดตามวิธีการโหลดข้อมูลลงใน SpatialPointsDataFrame และเรียกใช้การวิเคราะห์ที่น่าสนใจซึ่งฉันมีปัญหาคือวิธีการแปลงข้อมูลนั้นเป็นดาต้าเฟรมที่มี gridded เพื่อให้ฉันสามารถวิเคราะห์ฮิสโตแกรมได้ มีข้อเสนอแนะสำหรับการแปลงคะแนนเป็นกริดหรือไม่?

14 raster  r  kriging