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

ฉันจะคำนวณจำนวนเงินเฉลี่ย ฯลฯ ของคะแนนแรสเตอร์ (หลายแบนด์) ต่อรูปหลายเหลี่ยมของเวกเตอร์เลเยอร์ได้อย่างไร ฉันบอกว่านี่เรียกว่า "สถิติเชิงพื้นที่" ฉันลองกับ QGIS ก่อน มีวิธีที่จะทำ แต่มันช้าเกินไป (แปลงแรสเตอร์เป็นเวกเตอร์ตัดกับเวกเตอร์เลเยอร์ที่สองคำนวณเรขาคณิตส่งออกตัวเลขคำนวณสถิติด้วยสเปรดชีตหรือโปรแกรมอื่นนำเข้าผลลัพธ์ใช้เวลาตลอดไปสำหรับฉัน 350,000 คะแนนแรสเตอร์) ฉันยังได้รับคำใบ้ในการใช้ saga-gis ที่มี "สถิติโซน" แต่ขึ้นอยู่กับหมวดหมู่จากแรสเตอร์แบนด์ไม่ใช่รูปหลายเหลี่ยมจากเวกเตอร์เลเยอร์ ดังนั้นในการใช้สิ่งนี้ฉันจะต้องแปลงเวกเตอร์เลเยอร์ของฉันเป็นแรสเตอร์แล้วคำนวณสถิติ นี่เป็นวิธีที่ผิดในการแก้ไขปัญหานี้ จะไม่มีวิธีในการพิจารณาคะแนนแรสเตอร์ที่เป็นของรูปหลายเหลี่ยม 2 รูปขึ้นไปเพราะมันถูกตัดกันโดยเขตแดนรูปหลายเหลี่ยม ฉันคิดว่าสถิติที่มีรูปหลายเหลี่ยมควรสามารถจัดการกับสิ่งนี้ดังนั้นฉันจึงคิดว่าฉันยังไม่พบโมดูลที่ถูกต้อง Saga-gis มีโมดูลมากมาย โปรดแจ้งให้เราทราบว่าอันไหนเหมาะสมที่สุดสำหรับแอปพลิเคชันนี้

14 qgis  raster  statistics  saga 

ฉันมีแรสเตอร์ในเวลากลางวันของพื้นที่ศึกษาของฉัน ฉันมีเคาน์เตอร์มากกว่า 20 แห่งในพื้นที่ศึกษาที่นับ "สิ่ง" ด้วยเวลาประทับ ใช้PyEphemฉันกำหนดประเภทของแสงโดยรอบสำหรับแต่ละจุดข้อมูลที่เก็บรวบรวม สิ่งที่ฉันต้องการทำ: เพิ่มการซ้อนทับแบบโปร่งใสสีเทา / กึ่งเพื่อสะท้อนแสงโดยรอบ (พลบค่ำ, พระอาทิตย์ขึ้น, กลางวัน, พระอาทิตย์ตกและมืด) เป้าหมายสุดท้าย: สร้างการสร้างภาพข้อมูลตามเวลาของการรวบรวมข้อมูลที่เน้นเวลาของวัน หรือสามารถสร้างสแนปชอตในเวลาที่แน่นอนด้วยการซ้อนทับสีเทาที่เหมาะสม คิดว่านี่เป็นภาพ: แหล่ง

14 qgis  raster  arcgis-9.3  overlay 

ทุกคนสามารถแนะนำวิธีที่รวดเร็ว (และควรเป็นโอเพนซอร์ส) ในการสร้างภาพจาก Shapefile ตัวอย่างเช่นฉันมีรูปร่างไฟล์ที่มีขอบเขตของพัสดุฉันต้องการสร้างไฟล์ png เพื่อแสดงภาพเหล่านี้ สัญลักษณ์ที่จำเป็นสำหรับภาพนี้จะเป็นกระดูกเปลือยเปล่า (โครงร่างแข็งและเติม) ฉันรู้ว่า shp2img จาก MapServer จะทำเคล็ดลับ ฉันแค่อยากรู้ว่ามีอะไรอีกบ้าง การอ้างอิงถึงเครื่องมือที่คอมไพล์ล่วงหน้าหรือยินดีต้อนรับ api / sdks

14 open-source-gis  raster  shapefile 

นี่เป็นคำถามสองส่วนจริง: มีวิธีใดในการคำนวณความชันเฉลี่ยตามหน่วย (พูดกม. ²) สำหรับพื้นที่ที่กำหนด มีข้อกำหนดพิเศษเกี่ยวกับข้อมูลหรือไม่เช่นข้อมูลเพิ่มเติมนอกเหนือจากการวัดความสูง? GIS แบบโอเพ่นซอร์สใดที่ใช้วิธีการที่ใช้ในการคำนวณความชันเฉลี่ยบนกริด (เช่นไฟล์ TIF)

14 raster  open-source-gis  statistics  analysis  slope 

ฉันประสบความสำเร็จอย่าง จำกัด ในการสร้างรูปหลายเหลี่ยมเส้นเรียบจากชั้นแรสเตอร์ ฉันได้สร้างเส้นชั้นความสูงและเส้นเขตแดนและใช้ QGIS Union จากนั้นฉันก็ใช้ QGIS Polygonize จากนั้นใช้สถิติเชิงพื้นที่กับแรสเตอร์ดั้งเดิมเพื่อรับค่ารูปหลายเหลี่ยม วิธีที่สองที่ฉันพยายามคือการแปลงแรสเตอร์เป็นรูปหลายเหลี่ยม จากนั้นสร้างคอลัมน์คุณลักษณะใหม่ที่โค้งมนไปยังรูปร่างที่ฉันต้องการ จากนั้นฉันก็ยุบบรรทัดระหว่างค่าเดียวกัน จากนั้นใช้ v.generalize วิธีนี้เร็วกว่าเล็กน้อย แต่จะทำให้รูปหลายเหลี่ยมมีรอยหยัก / ไม่สลบ มีวิธีที่ง่ายกว่านี้หรือเปล่าที่ใช้งานได้จริงตลอดเวลา ดูเหมือนว่าควรมีคำสั่งง่าย ๆ ที่ทำกระบวนการทั้งหมด

14 qgis  raster  contour 

ฉันกำลังเขียนยูทิลิตี้ง่าย ๆ ในการครอบตัดไฟล์แรสเตอร์ geotiff หลายแบนด์ไปยังพื้นที่ (เล็กกว่า) เดียวกัน เมื่อใช้ gdalwarp ฉันสามารถครอบตัดไฟล์ได้อย่างง่ายดายโดยใช้รูปทรงการตัดรูปหลายเหลี่ยมเดียว: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif อย่างไรก็ตามพื้นที่จริงที่ฉันต้องการคลิปจะถูกกำหนดไว้ในขั้นต้นเสมอโดยไฟล์แรสเตอร์ geotiff อื่นไม่ใช่ไฟล์รูปร่าง มันจะดีถ้าฉันสามารถใช้ขอบเขตของแรสเตอร์นั้นเป็นไฟล์รูปวาด แต่ฉันไม่แน่ใจว่าจะทำอย่างไร น่าแปลกใจที่สิ่งต่อไปนี้ใช้งานไม่ได้ (มันไม่ได้ทำให้เกิดข้อผิดพลาด แต่ก็ไม่ได้สร้างอะไร): gdalwarp -cutline clipper.tif-crop_to_cutline input.tif output.tif ดังนั้นคำถามของฉันคือมีวิธีให้แรสเตอร์gdalwarp -cutline? อีกวิธีหนึ่งมีฟังก์ชัน gdal อื่นที่สามารถคลิปแรสเตอร์โดยใช้แรสเตอร์อื่นได้หรือไม่ หากไม่มีสิ่งเหล่านี้เป็นไปได้มีวิธีที่ง่ายมากในการสร้างรูปร่างไฟล์ที่มีรูปหลายเหลี่ยมเดียวที่กำหนดโดยขอบเขตของแรสเตอร์หรือไม่? รหัสนี้จะถูกห่อหุ้มด้วยสคริปต์ไพ ธ อนที่กว้างขวางกว่าดังนั้นฉันสามารถใช้ยูทิลิตี้บรรทัดคำสั่ง gdal หรือการผูกไพ ธ อนใด ๆ สำหรับ gdal ในฐานะที่เป็นบันทึกด้านข้างฉันรู้ว่าฉันสามารถสร้างรูปทรงการตัดได้อย่างง่ายดายซึ่งครอบคลุมขอบเขตของแรสเตอร์ของฉันใน QGIS ฉันอาจเลิกทำหากไม่พบวิธีแก้ปัญหาที่ตรงไปตรงมา แต่ท้ายที่สุดฉันจะเลิกใช้ยูทิลิตีนี้ในหลายสิบถ้าไม่ใช่พื้นที่นับร้อยซึ่งเป็นส่วนหนึ่งของการวิเคราะห์อัตโนมัติขนาดใหญ่ดังนั้นฉันจึงไม่อยากจะเบื่อ ขั้นตอนแบบแมนนวลแม้ว่าจะง่ายมาก

14 python  raster  gdal  clip  gdalwarp 

ฉันมีคำถามเกี่ยวกับการรวมตัวเชิงพื้นที่ใน R สิ่งที่ฉันพยายามทำคือการรวมชุดข้อมูลแบบจุดไปยังกริด ฉันไม่แน่ใจอย่างไรก็ตามวิธีการทำเช่นนี้เพราะฉันมีประสบการณ์น้อยกับสิ่งเหล่านี้ ฉันหวังว่าทุกคนในคุณอาจมีแนวทางที่เป็นประโยชน์ / แนวทางแก้ไขที่เป็นไปได้ จุดได้เปรียบของฉันคือชุดข้อมูลที่มีข้อมูลอ้างอิงทางภูมิศาสตร์สำหรับเหตุการณ์ความขัดแย้งในแอฟริกา (ดู www.acleddata.com) คะแนนจะถูกอ้างอิงทางภูมิศาสตร์ด้วยพิกัดละติจูด / ลองจิจูดและมีข้อมูลเกี่ยวกับประเภทและเวลาของเหตุการณ์ สิ่งที่ฉันต้องการทำคือรวมคะแนนเหล่านี้เข้ากับกริด 1x1 องศา ดังนั้นกริดเซลล์ควรมีข้อมูลของจุดข้อมูลถ้ามีเหตุการณ์เกิดขึ้นภายในกริดเซลล์นั้น ผลิตภัณฑ์สุดท้ายนี้ควรเป็น data frame หรือบางสิ่งที่ฉันสามารถส่งออกไปยัง csv-file ได้เนื่องจากข้อมูลนั้นมีวัตถุประสงค์เพื่อใช้ในชุด data panel สำหรับการวิเคราะห์ทางสถิติ จนถึงตอนนี้ฉันโหลดและลงจุดข้อมูลและ shapefile โดยใช้รหัสด้านล่าง ฉันเชื่อว่าฉันควรใช้ฟังก์ชันoverจากแพ็คเกจspเพื่อรวม แต่ฉันไม่ทราบวิธี หวังว่าหนึ่งในคุณสามารถช่วย รหัสฉันใช้เพื่อให้ห่างไกลสามารถพบได้ที่นี่กับผลภาพที่สอดคล้องกันมากกว่ามี ข้อเสนอแนะสำหรับการทำเช่นนี้ใน QGIS ก็ยินดีต้อนรับเช่นกัน

14 raster  r  point  aggregation 

ใน ArcGIS 9.3 ฉันมีไฟล์แรสเตอร์ที่ฉันต้องการ "ปกปิด" ด้วยรูปหลายเหลี่ยมให้ค่าเดียวกับเซลล์แรสเตอร์ทั้งหมดที่ซ้อนทับรูปหลายเหลี่ยม ฉันจะทำสิ่งนี้ได้อย่างไร นี่คือสกรีนช็อตที่รูปหลายเหลี่ยมสีดำเป็นเลเยอร์รูปหลายเหลี่ยมของฉันและส่วนที่เหลือเป็นแบบแรสเตอร์

14 arcgis-desktop  raster  polygon  arcgis-9.3  masking 

ฉันมีแรสเตอร์สองตัว: แบบจำลองระดับความสูงแบบดิจิทัล (DEM) และวงดนตรี Landsat ความละเอียดกริดของทั้งคู่อยู่ที่ 30 ม. แต่กริดไม่เหมือนกัน (ดูภาพ: แรสเตอร์สีเข้มด้านบนคือ DEM; แรสเตอร์แสงที่อยู่ด้านล่างคือแถบแลนดัต) ฉันต้องการทำการแปลงรูปแบบของวง Landsat เพื่อให้กริดของมันสอดคล้องกับกริด DEM อย่างสมบูรณ์ ฉันได้ลองใหม่อีกครั้งกับวง Landsat ถึง 30 เมตรใน ArcGIS โดยใช้ DEM เป็น Snap Raster ในการตั้งค่าสภาพแวดล้อม แต่สิ่งนี้ไม่มีผลใด ๆ ฉันได้ลองตัดวง Landsat ด้วยขนาดของรูปหลายเหลี่ยมเวกเตอร์ (ซึ่ง DEM นั้นถูกตัดด้วย) ในขณะที่ตั้งค่า DEM เป็น Snap Raster แต่จะไม่มีประโยชน์ ฉันสงสัยว่าสิ่งนี้เกิดขึ้นเพราะแรสเตอร์ทั้งสองมีความละเอียดเท่ากันอยู่แล้ว (เฉพาะกับกริดที่แตกต่างกัน) แต่ฉันต้องการให้กริดของพวกเขาอยู่ในแนวเดียวกันเพื่อให้สามารถทำการวิเคราะห์เพิ่มเติมได้ ฉันสามารถทำได้ใน ArcGIS Desktop …

14 arcgis-desktop  raster  alignment 

ฉันมีแรสเตอร์จำนวนเต็มที่ฉันต้องการสร้างตารางแอตทริบิวต์แรสเตอร์โดยใช้ Python และ GDAL ฉันสามารถสร้างตารางแอตทริบิวต์ GDAL raster ใน Python ดังนี้ >>> rat = gdal.RasterAttributeTable() ใช้งานได้ดีอย่างที่เราเห็น: >>> rat <osgeo.gdal.RasterAttributeTable; proxy of <Swig Object of type 'GDALRasterAttributeTableShadow *' at 0x0000000002A53D50> > ตารางที่สร้างขึ้นจึงไม่มีแถวหรือคอลัมน์: >>> rat.GetRowCount() 0 >>> rat.GetColumnCount() 0 ฉันสร้างคอลัมน์ชื่อ "ค่า" เพื่อเก็บค่าที่ไม่ซ้ำใน raster: >>> rat.CreateColumn("Value", gdalconst.GFT_Integer, gdalconst.GFU_MinMax) 0 นี่เป็นเรื่องปกติและมีการอัปเดตจำนวนคอลัมน์: >>> rat.GetColumnCount() 1 ตอนนี้ฉันต้องเพิ่มค่า …

14 python  raster  gdal  attribute-table 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา พื้นหลังนี่เป็นคำถามที่สองของฉันที่เกี่ยวข้องกับการอ้างอิงแผนที่เชิงภูมิศาสตร์เพื่อแสดงภาพบนระบบพิกัดที่แตกต่างกันและร่วมกับชั้นข้อมูลอื่น ๆ คำถามก่อนหน้านี้อยู่ที่การแปลงรูปแผนที่โดยไม่ใช้ข้อมูล meta-data โดยพลการเป็นโครงการ QGIS ปัญหาเป้าหมายของฉันคือการอ้างอิงทางภูมิศาสตร์แผนที่นี้: ดูเหมือนว่าจะไม่เป็น Plate-Carrée ดังนั้นใน QGIS ฉันสร้างจุดควบคุมที่สมเหตุสมผลหลายจุดซึ่งเพื่อความสมบูรณ์ฉันได้แนบไว้ที่ด้านล่าง [อ้างอิง: 1] ฉันให้ QGIS Georeferencer SRS เป้าหมายเดียวกันกับไฟล์โครงการของฉัน EPSG: 4326 ฉันได้รับผลลัพธ์ที่แย่เป็นพิเศษกับ Helmert และการแปลงพหุนาม แต่ได้ภาพที่สมเหตุสมผลพร้อมแผ่นโค้งบาง ๆ (ซึ่งทำให้ผลลัพธ์ทางภูมิศาสตร์ผ่านจุดควบคุมของฉัน) อย่างไรก็ตามแม้ผลลัพธ์นี้จะแย่เช่นที่ละติจูดสูงขึ้น (ดูชายฝั่งรัสเซียทางตอนเหนือของญี่ปุ่น) นี่เป็นสกรีนช็อตของหน้าจอ QGIS ของฉันโดยใช้พื้นหลังของ Natural Earth ทางเลือกอื่นฉันลองทำแบบฝึกหัดที่คล้ายกันด้วยเครื่องมือที่ใช้งานง่ายกว่ามากที่ MapWarper: ดูผลลัพธ์และจุดควบคุมที่http://mapwarper.net/maps/758#Preview_Map_tabซึ่งฉันได้ผลลัพธ์ที่แย่ลง ที่ฉันเพิ่มจุดควบคุมน้อยลง) คำถามสั้น ๆ มีกลอุบายใด ๆ ที่ฉันพลาดในการรับตำแหน่งทางภูมิศาสตร์ที่ดีหรือไม่? …

13 qgis  raster  coordinate-system  georeferencing 

ฉันกำลังทำงานกับไฟล์ TRMM 3B31 เกี่ยวกับการตกตะกอนของ NASA เมื่อฉันโหลดมันลงบน QGIS เลือกการฉาย WGS84 พร้อมกับไฟล์รูปร่างอื่น ๆ แผนที่จะแสดงในแนวตั้งเช่น QGIS ตีความเป็นละติจูดว่าอะไรคือเส้นแวงและในทางกลับกัน ฉันได้ลองใช้ตัวเลือก gdal_translate -a_ullr แล้วเลือกพิกัดที่ดี แต่มันไม่ได้แก้ปัญหา: ภาพถูกบีบอัดและยืดออกเพื่อตอบสนองความต้องการของฉัน แต่นี่ไม่ใช่สิ่งที่ฉันต้องการ เมื่อคุณโหลดไฟล์. csv คุณสามารถเลือกว่าจะต้องเป็นละติจูดและลองจิจูดอะไรมีวิธีการทำเช่นเดียวกันสำหรับ rasters หรือไม่ เนื่องจากฉันมีแรสเตอร์ 193 มันจะดีมากถ้าคุณสามารถแนะนำคำสั่งที่สามารถทำซ้ำได้ ภาพแสดงแผนที่ของประเทศอินเดียและแถบแรกของไฟล์ 3B31 เช่นการตกตะกอนของพื้นผิว เส้นสีดำสองเส้นคือเส้นขนาน 40 และ -40 ฉันต้องการให้พวกมันเป็นแนวนอนและด้วยเหตุนี้ข้อมูลฝนจึงคาบเกี่ยวกับอินเดีย

13 qgis  raster  latitude-longitude  trmm 

ในบางครั้งเราใช้ข้อมูล LiDAR เพื่อสร้างแบบจำลองใน ArcGIS ก่อนหน้านี้เราจะให้ข้อมูล. xyz หรือ .ascii ของข้อมูล แต่ไม่นานมานี้เรายังได้รับไฟล์. Atlas เพิ่มเติมจาก. xyz และ ascii ผมได้อ่านข้อความนี้บทความที่ระบุ"การมี LIDAR ในรูปแบบ LAS อาจจะเป็นที่ชัดเจนที่จะริเริ่ม แต่ไม่ให้ผู้ที่ใหม่กับการใช้ข้อมูล LIDAR"และ"มันไบนารีที่มีประสิทธิภาพ, การสนับสนุนอย่างกว้างขวางและ ArcGIS รูปแบบการทำงานที่ดีที่สุดด้วย" ฉันเป็นหนึ่งในคนเหล่านั้นที่อยู่ในหมวดหมู่ใหม่ในการใช้รูปแบบ. Atlas ตามปกติด้วย. xyz และ ascii ฉันจะใช้ Global Mapper เพื่อแปลงเป็นตารางแรสเตอร์ (สำหรับใช้ในการสร้างแบบจำลอง) และดูเหมือนว่าฉันจะต้องแปลงไฟล์. Atlas ด้วย (เป็นจุดของคลาวด์) ฉันเดาคำถามที่ฉันถามคือ: มีประโยชน์ของการใช้. Atlas เปรียบเทียบกับ. xyz หรือ ascii หรือไม่หากฉันยังต้องการแปลงเป็นรูปแบบแรสเตอร์อื่น (GRID หรือ. …

13 arcgis-desktop  raster  arcgis-10.1  lidar  las 

ฉันพยายามสร้างภาพแรสเตอร์จากฐานข้อมูลจุดที่เว้นระยะไม่สม่ำเสมอ ข้อมูลดูเหมือนว่า - > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6 255094.8 2633781 88.90244 ฉันต้องการค่า 'z' เหล่านี้ภายในตาข่ายที่ฉันสร้างขึ้น # Create a fine mesh grid my_mesh=expand.grid(seq(min(s100_ras$Y),max(s100_ras$Y),l=100), seq(min(s100_ras$X),max(s100_ras$X),l=100)) ฉันต้องการให้ค่า z ถูกกำหนดเป็น 'NA' สำหรับจุดตาข่ายที่อยู่นอกจุดข้อมูล คะแนนเหนือตาข่ายมีลักษณะดังนี้: https://drive.google.com/file/d/0B6GUNg-8d30vYzlwTkhvaHBFTnc/edit?usp=sharingเมื่อฉันลงจุด …

13 raster  r  vector-grid 

ดูลิงค์นี้สำหรับรายละเอียดเพิ่มเติม ปัญหา: ฉันต้องการวนรอบraster อย่างต่อเนื่อง (ที่ไม่มีตารางคุณลักษณะ), เซลล์โดยเซลล์และรับค่าของเซลล์ ฉันต้องการใช้ค่าเหล่านั้นและเรียกใช้ conditionals กับพวกเขาเลียนแบบขั้นตอนพีชคณิตแผนที่รายละเอียดด้านล่างโดยไม่ต้องใช้เครื่องคิดเลขแรสเตอร์ ตามคำขอของความคิดเห็นด้านล่างฉันได้เพิ่มรายละเอียดที่ให้ข้อมูลพื้นฐานเกี่ยวกับปัญหาและให้เหตุผลถึงความจำเป็นในการใช้วิธีการดังกล่าวในส่วนด้านล่างที่เรียกว่า "การวิเคราะห์ที่จำเป็น:" การวิเคราะห์ที่เสนอด้านล่างในขณะที่เกี่ยวข้องกับปัญหาของฉันโดยการให้ข้อมูลพื้นฐานไม่จำเป็นต้องมีคำตอบ ขอบเขตของคำถามเกี่ยวข้องกับการวนซ้ำผ่านแรสเตอร์อย่างต่อเนื่องเพื่อรับ / ตั้งค่าของเซลล์ การวิเคราะห์ที่จำเป็น: หากเงื่อนไขใด ๆ ต่อไปนี้เป็นจริงให้กำหนดค่าเซลล์เอาต์พุตเป็น 1 ให้ค่าเซลล์เอาต์พุตเป็น 0 เท่านั้นหากไม่มีเงื่อนไขใดที่จะสำเร็จ เงื่อนไข 1:หากค่าของเซลล์มากกว่าเซลล์บนและล่างให้ค่าเป็น 1: Con("raster" > FocalStatistics("raster", NbrIrregular("C:\filepath\kernel_file.txt"), "MAXIMUM"), 1, 0) ไฟล์เคอร์เนลมีลักษณะดังนี้: 3 3 0 1 0 0 0 0 0 1 0 เงื่อนไข 2:หากค่าของเซลล์มากกว่าเซลล์ซ้ายและขวาให้ค่าเป็น 1: Con("raster" > …

13 python  raster  arcobjects  gdal  c#