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

คำถามของฉันเกี่ยวกับการแยกค่าจาก rasters ที่ตำแหน่งของจุด ด้วยฟังก์ชั่นการแยกมันง่ายมากและฟังก์ชั่นให้ฉันดาต้าเฟรมที่มีค่าของตัวแปรทั้งหมดในจุด ฉันต้องการที่จะมีใน dataframe พิกัดของแต่ละจุด ฉันจะทำให้เกิดขึ้นได้อย่างไร เป็นไปได้ไหมที่จะบอกกับ R ว่าเมื่อทำการแยกค่าจากแรสเตอร์จะเพิ่มคอลัมน์ของตำแหน่งที่ตั้ง นี่คือขั้นตอนของฉัน: presencias=read.table("c:/SDM_R/presencias/P_lentiscus_pres.csv",header=TRUE,sep=";") lista_variables <-list.files(path="Variables_modelizacion/solo_ascii",pattern='*.asc',full.names=TRUE) variables <- stack(lista_variables) variables_presencia<-extract(variables, presencias) ผลลัพธ์เป็นดังนี้: > bio1 bio12 bio18 bio2 bio4 > 90 875 165 95 4886 > 115 1085 158 83 4075 > 135 1153 153 67 3402 > 85 1026 137 99 …

13 raster  r 

ฉันใช้ ArcGIS 10.1 และต้องการสร้างแรสเตอร์ใหม่โดยใช้แรสเตอร์สองรูปแบบมาก่อน RasterToNumPyArrayมีตัวอย่างที่ดีที่ฉันต้องการที่จะปรับตัว import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape = (myArray.shape[0],1) myArrayPerc = (myArray * 1.0)/ myArraySum newRaster = arcpy.NumPyArrayToRaster(myArrayPerc) newRaster.save("C:/output/fgdb.gdb/PercentRaster") ปัญหาคือมันแถบการอ้างอิงเชิงพื้นที่และขนาดของเซลล์ ฉันคิดว่ามันต้องทำ arcpy.env แต่ฉันจะตั้งพวกมันตามอินพุตแรสเตอร์ได้อย่างไร ฉันไม่สามารถคิดออก รับคำตอบของลุคนี่คือวิธีแก้ปัญหาเบื้องต้นของฉัน โซลูชันทั้งสองของลุคตั้งค่าการอ้างอิงขอบเขตและขนาดเซลล์อย่างถูกต้อง แต่วิธีแรกไม่ได้เก็บข้อมูลในอาเรย์อย่างถูกต้องและเอาท์พุทแรสเตอร์เต็มไปด้วยข้อมูลทั้งหมด วิธีที่สองของเขาใช้งานได้เป็นส่วนใหญ่ แต่ที่ฉันมีพื้นที่ขนาดใหญ่ของโนดาต้ามันเต็มไปด้วยศูนย์บล็อกและ 255s สิ่งนี้อาจเกี่ยวข้องกับวิธีที่ฉันจัดการเซลล์โนดาต้าและฉันก็ค่อนข้างแน่ใจว่าฉันกำลังทำมันอย่างไร (ควรเป็น Q อีกอันหนึ่ง) ฉันรวมภาพสิ่งที่ฉันพูดถึง #Setting the raster properties directly import …

13 arcpy  raster  arcgis-10.1  coordinate-system  numpy 

ฉันต้องการความช่วยเหลือ ... ฉันมีปัญหากับ QGIS ฉันมีแรสเตอร์บางตัวที่ดูเหมือนจะมีค่า Null ถึงแม้ว่าฉันจะจัดประเภทใหม่เพื่อให้ได้ Raster แบบบูลีน :( ฉันจะเปลี่ยนค่า NULL เป็น 0 ได้อย่างไร

13 qgis  raster  reclassify  nodata 

ฉันมี DEM สร้างไฟล์ทั้งหมดที่เกี่ยวข้องกับลุ่มน้ำรวมถึงการสะสมทิศทางลำธาร ฯลฯ ... ฉันยังคงสร้างเขื่อน ฉันขอให้คุณโปรดแจ้งให้เราทราบวิธีการสร้าง 'เขื่อน' เทียมข้ามลำธาร / แม่น้ำ กรุณาแนะนำซอฟต์แวร์ / วิธีการใด ๆ

13 raster  dem  software-recommendations 

มีวิธีที่ตรงไปตรงมาในการสุ่มตัวอย่างแรสเตอร์เพื่อให้ผลลัพธ์ของกระบวนการเป็นแรสเตอร์หรือไม่? ฉันใช้ตัวอย่างที่ฉันพบในr-sig-geo รายการและฉันได้ลองใช้sampleRandomฟังก์ชันในrasterแพ็คเกจด้วย วิธีการทั้งสองนี้สร้างผลลัพธ์ที่ฉันไม่แน่ใจว่าจะเปลี่ยนเป็นแรสเตอร์ได้อย่างไร ฉันไม่สามารถค้นหาวิธีการหลังจากค้นหาชุดค่าผสมหลายชุดของ "SpatialPointsDataFrame raster" library(raster) # read in raster rasterSource <- 'landsat.TIF' r <- raster(rasterSource) # convert to spatial points data frame r.spgrd<-as(r,"SpatialPointsDataFrame") # elminate NA values r.spgrd = r.spgrd[!is.na(r.spgrd[[1]]),] # sample points selectedPoints = sample(1:length(r.spgrd[[1]]), 1000) r.sampled = r.spgrd[selectedPoints,] # try to make spgrd into a …

13 raster  r 

ฉันมีแรสต์ riff ในฐานข้อมูล PostGIS 2.0 ของฉัน ฉันต้องการรับค่าของแรสเตอร์ที่ตำแหน่ง lat / lon ที่ระบุ ฉันมี: PostgreSQL-9.1 postgis-2.0-svn พร้อมการสนับสนุนแบบราสเตอร์ raster ที่นำเข้าสู่ PostGIS SRID ของมันคือ 3035

13 postgis  raster  point 

การอ่านวรรณกรรม GIS บางครั้งก็รู้สึกว่าผู้เขียนใช้คำสองคำที่มีความหมายเฉพาะในใจ ฉันมักจะปฏิบัติต่อข้อกำหนดแทนกัน ฉันผิดหรือเปล่า? อ่านบรรทัดแรกในการเชื่อมโยง ดูว่าบริการรูปภาพให้การเข้าถึง raster (และรูปภาพ) ได้อย่างไรโดยบอกนัยถึงความหมายที่ต่างกันสำหรับคำสองคำนี้

13 raster  imagery  terminology 

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

13 raster  analysis  modelling 

ฉันต้องการสร้างรูปหลายเหลี่ยมสองรูป หนึ่งในนั้นคือรูปสี่เหลี่ยมผืนผ้าของแรสเตอร์ของฉัน ฉันรู้ว่าฉันสามารถได้รับขอบเขตของการใช้ raster r <- raster("band5.tif") e <- extent(r) plot(e)ให้ฉันนี้ ฉันจะสร้างหนึ่งรูปหลายเหลี่ยมของขอบเขตของแรสเตอร์ที่แสดงด้านล่างได้อย่างไร นี่คือลักษณะแรสเตอร์ของฉัน:

13 raster  r  rgdal 

ฉันต้องเปลี่ยนมันเป็น UTM เพื่อให้บัฟเฟอร์ทำงานได้ wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) ในระหว่างการสร้างบัฟเฟอร์ซึ่งเป็นบรรทัดสุดท้ายของรหัสมันให้คำเตือนนี้: Warning message: In couldBeLonLat(x) : raster has a longitude/latitude CRS, but coordinates do not match that นี่คือข้อมูล summary(wets.r) layer Min. 1 1st Qu. 1 Median 2 3rd Qu. 9 Max. 11 …

13 raster  coordinate-system  shapefile  r 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันมีรูปร่างถนนขนาดใหญ่ (~ 70MB) และต้องการแปลงเป็นแรสเตอร์ด้วยความหนาแน่นของถนนในแต่ละเซลล์ เป็นการดีที่ฉันต้องการทำสิ่งนี้ใน R พร้อมกับเครื่องมือบรรทัดคำสั่งของ GDAL หากจำเป็น วิธีการเริ่มต้นของฉันคือการคำนวณความยาวของส่วนของเส้นตรงในแต่ละเซลล์ตามเธรดนี้โดยตรง สิ่งนี้ให้ผลลัพธ์ที่ต้องการ แต่ค่อนข้างช้าแม้สำหรับรูปร่างที่เล็กกว่าของฉันมาก นี่คือตัวอย่างที่ง่ายมากซึ่งค่าเซลล์ที่ถูกต้องชัดเจน: require(sp) require(raster) require(rgeos) require(RColorBrewer) # Create some sample lines l1 <- Lines(Line(cbind(c(0,1),c(.25,0.25))), ID="a") l2 <- Lines(Line(cbind(c(0.25,0.25),c(0,1))), ID="b") sl <- SpatialLines(list(l1,l2)) # Function to calculate lengths of lines in given raster …

13 raster  r  road  kernel-density  spatstat 

ฉันใหม่สำหรับ R และใช้แพ็คเกจแรสเตอร์ ฉันมีปัญหาในการแยกรูปหลายเหลี่ยมจากไฟล์แรสเตอร์ที่มีอยู่ ถ้าฉันใช้ extract(raster, poly_shape) ฟังก์ชั่นบนแรสเตอร์มันมักจะสร้างรายการที่มีข้อมูล สิ่งที่ฉันต้องการคือการแยกไฟล์แรสเตอร์อื่นที่ฉันสามารถโหลดด้วย ArcGIS อีกครั้ง หลังจากอ่านอีกเล็กน้อยฉันคิดว่าฟังก์ชั่นการครอบตัดคือสิ่งที่ฉันต้องการจริงๆ แต่เมื่อลองใช้ฟังก์ชั่นนี้ crop(raster, poly_shape) ฉันได้รับข้อผิดพลาดนี้: Error in .local(x, y, ...) : extents do not overlap In addition: Warning message: In intersect(extent(x), extent(y)) : Objects do not overlap ไฟล์ raster และ poly_shape เหมือนกันทั้งสองฟังก์ชั่น คุณช่วยบอกฉันทีว่าอะไรที่นี่ผิด มันถูกต้องไหมที่ฟังก์ชั่นการครอบตัดสร้างอีกแบบแรสเตอร์และไม่ใช่รายการ? แก้ไข : ฟังก์ชั่นขอบเขต () ไม่ทำงานสำหรับฉัน …

13 raster  shapefile  r  line 

เมื่อฉันประมวลผลแรสเตอร์ใน QGIS มันจะสร้างแถบที่สอง ปรากฏว่าแถบนั้นไม่มีข้อมูลที่เป็นประโยชน์ ฉันชอบใช้ Raster Colors Tool และจะไม่ยอมรับแรสเตอร์ที่มีอะไรมากกว่าในวง ฉันจะลบวงเสริมได้อย่างไร

13 raster  qgis 

มีวิธีตั้งค่าพิกเซล 0 ค่าทั้งหมดในแรสเตอร์เป็นโนดาต้าโดยใช้ GDAL หรือไม่ ฉันต้องการใช้ fillnodata เพื่อเติมหลุมบางส่วน แต่หลุมถูกตั้งค่าเป็น 0 ดังนั้นฉันต้องลบออกก่อน

13 raster  gdal 

ฉันพยายามใช้ QGIS v2.0.1 เพื่อทำงานเดียวกันกับที่ฉันสามารถทำได้โดยใช้ ArcMap v10.1 และฉันพบปัญหาในการค้นหาเครื่องมือที่เทียบเท่า ใน ArcMap ฉันใช้เครื่องมือ Focal Statistics (หมายถึง) นักวิเคราะห์เชิงพื้นที่และฉันไม่ทราบว่า GRASS หรือ SAGA มีความเทียบเท่าใน QGIS ไม่มีใครรู้ว่าเครื่องมือนี้จะถูกเรียกถ้ามันมีอยู่?

13 qgis  raster  grass  spatial-analyst