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

ฉันกำลังทำงานกับ DEM ซึ่งอยู่ในพื้นที่ที่มีความโล่งใจเล็กน้อย ฉันต้องการ 'เผา' เครือข่ายแม่น้ำลงใน DEM เพื่อให้ฉันสามารถคำนวณการสะสมการไหลและความยาวของการไหลได้อย่างถูกต้อง ฉันใช้ ArcGIS Desktop 10

15 arcgis-desktop  arcgis-10.0  dem  rasterization  hydrography 

เป็นไปได้ไหมที่จะแบ่ง shapefile ต่อคุณสมบัติใน python (ดีที่สุดจะเป็นวิธีการแก้ปัญหาที่ฉันสามารถบันทึกวัตถุเวกเตอร์ผลลัพธ์ลงในหน่วยความจำแทนลงดิสก์ได้ชั่วคราว) เหตุผล: ฉันต้องการใช้ฟังก์ชั่น gdal rasterizeLayer กับหลาย ๆ เซตย่อยของ shapefile ฟังก์ชั่นต้องการวัตถุ osgeo.ogr.Layer mkay ฉันลองไปรอบ ๆ แล้วมันสามารถทำงานได้ดังต่อไปนี้ คุณสามารถรับรูปทรงเรขาคณิตของวัตถุเลเยอร์ gdal ต่อคุณสมบัติดังต่อไปนี้ # Load shape into gdal shapefile=str(vectorPath) layer_source = ogr.Open(shapefile) lyr = layer_source.GetLayer(0) for i in range(0,lyr.GetFeatureCount()): feat = lyr.GetFeature(i) ge = feat.geometry() ตอนนี้ฉันแค่ต้องรู้วิธีการสร้างวัตถุ osgeo.ogr.layer ตามรูปทรงเรขาคณิตนี้ สำหรับการชี้แจง ฉันต้องการฟังก์ชั่นในรหัส ogr …

15 python  gdal  rasterization  osgeo  split-by-attribute 

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

14 raster  r  line  rasterization 

ฉันพยายามเบิร์นไฟล์ไปยังแรสเตอร์โดยใช้ RasterizeLayer ของ GDAL ฉันสร้างพื้นที่แรสเตอร์ที่น่าสนใจจากรูปร่างที่แตกต่างกันก่อนกำหนดขนาดพิกเซลที่เฉพาะเจาะจง AOI นี้ทำหน้าที่เป็นฐานสำหรับการแรสเตอร์ต่อไปนี้ทั้งหมด (จำนวน collumns และแถวเดียวกัน, การฉายภาพและ geotransform เดียวกัน) ปัญหาเกิดขึ้น แต่เมื่อฉันไปเผารูปร่างเพื่อแรสเตอร์ของตัวเองตามขนาดพิกเซลและประมาณการเดียวกัน ลิงค์ด้านล่าง (มีตัวแทนไม่เพียงพอที่จะโพสต์ภาพ) แสดงไฟล์รูปร่างดั้งเดิมในผิวสีแทนและสีชมพูเข้มที่ RasterizeLayer เขียนข้อมูล สีชมพูอ่อนเป็นค่า nodata สำหรับข้อมูลแรสเตอร์สีชมพูเข้ม สีเทาคือ AOI ตามการเบิร์นไฟล์รูปร่างเสร็จสมบูรณ์ เมื่อพิจารณาถึงขอบเขตของรูปหลายเหลี่ยม shapefile ฉันคาดว่าจะเห็นค่าการเบิร์นที่มุมด้านล่างทั้งสองรวมทั้งพิกเซลสองพิกเซลที่อยู่ใต้ข้อมูลที่แสดง เห็นได้ชัดว่านี่ไม่ใช่กรณี ดังต่อไปนี้เป็นรหัสที่ฉันใช้ในการสร้างสิ่งเหล่านี้ รูปร่างทั้งหมดถูกสร้างขึ้นโดยใช้ QGIS และถูกสร้างขึ้นทั้งหมดในการฉายภาพเดียวกัน (ควรสังเกตว่าการ gridding ในภาพที่แสดงเป็นเพียงเพื่อให้ความคิดเกี่ยวกับขนาดพิกเซลที่ฉันใช้) from osgeo import ogr from osgeo import gdal aoi_uri = 'AOI_Raster.tif' aoi_raster = …

12 python  raster  gdal  rasterization 

ฉันต้องการสร้างแรสเตอร์ที่มีความละเอียด 25 เมตร× 25 เมตรซึ่งแต่ละเซลล์มีระยะห่างจากแนวชายฝั่งที่ใกล้ที่สุดซึ่งคำนวณจากจุดศูนย์กลางของเซลล์ การทำเช่นนี้ทั้งหมดที่ฉันต้องเป็นshapefile ของแนวชายฝั่งของนิวซีแลนด์ ฉันได้ลองทำตามบทช่วยสอนของ Dominic Royeเพื่อใช้ใน R ซึ่งใช้งานได้ดี มันปรับลงไปที่ความละเอียดประมาณ 1 km × 1 km แต่ถ้าฉันพยายามที่จะใช้ RAM ที่สูงกว่านี้ก็ต้องใช้งานได้ดีกว่าบนพีซีของฉัน (ต้องใช้ RAM 70 ~ GB) หรืออื่น ๆ ที่ฉันมีการเข้าถึงด้วย ในการบอกว่าฉันคิดว่านี่เป็นข้อ จำกัด ของ R และฉันสงสัยว่า QGIS อาจมีวิธีการคำนวณที่มีประสิทธิภาพมากกว่าในการสร้างแรสเตอร์นี้ แต่ฉันใหม่สำหรับมันและฉันไม่สามารถเข้าใจได้ว่าจะทำอย่างไร ฉันได้ลองติดตามการสร้างแรสเตอร์ด้วยระยะทางถึงฟีเจอร์โดยใช้ QGIS แล้วหรือยัง เพื่อสร้างใน QGIS แต่มันกลับข้อผิดพลาดนี้: _core.QgsProcessingException: ไม่สามารถโหลดเลเยอร์แหล่งที่มาสำหรับ INPUT: C: /..../ Coastline / …

10 qgis  shapefile  r  rasterization 

ฉันกำลังแปลงเวกเตอร์เป็นแรสเตอร์ใน R อย่างไรก็ตามกระบวนการนี้ยาวเกินไป มีความเป็นไปได้ไหมที่จะนำสคริปต์ไปใช้ในการประมวลผลแบบมัลติเธรดหรือ GPU เพื่อให้เร็วขึ้น? สคริปต์ของฉันเป็นเวกเตอร์ที่แรสเตอร์ r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') r.raster คลาส: RasterLayer ขนาด: 9636, 11476, 110582736 (nrow, ncol, ncell) ความละเอียด: 10, 10 (x, y) ขอบเขต: 505755, 620515, 8555432, …

9 raster  r  rasterization  rgdal  parallel 

มี 100 SpatialPointsDataFrames ในพื้นที่ทำงานของฉัน ฉันโหลดมันด้วยวิธีนี้: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) พวกเขาถูกเรียกเช่นนี้: XYhectareTravelTimes_ez10.*shp XYhectareTravelTimes_ez11.*shp XYhectareTravelTimes_ez12.*shp เป็นต้น ฉันจะแปลงพวกมันเป็นแรสเตอร์วนซ้ำผ่านพื้นที่ทำงานได้อย่างไร ฉันเป็นผู้ใช้ R ใหม่มากและหวังว่าจะได้รับความช่วยเหลือ ขอบคุณมาก.

9 r  rasterization