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

มีคำที่สามารถใช้อธิบายกริดแรสเตอร์ตั้งแต่สองรายการขึ้นไปที่มีความละเอียดและขอบเขตเท่ากันหรือไม่? ฉันกำลังตกปลารอบ ๆ ด้วยคำว่า "ประจวบ" และ "coextensive" แต่ไม่มีอะไรที่รู้สึกค่อนข้างถูกต้อง และฉันอยากหาบางสิ่งที่มีขนาดกะทัดรัดกว่า "rasters ที่มีขนาดและความละเอียดเท่ากัน"

9 raster  terminology 

ฉันมี rasters หลายรูปแบบ TIFF ที่มีตารางสีให้พวกเขา: จากนั้นฉันสร้างชุดข้อมูล Mosaic พร้อมสเปคที่ตรงกับภาพและเพิ่ม tiffs ทั้งหมดของฉันลงไป (คุณสามารถเพิกเฉยต่อคำเตือนนี้เพื่อแสดงให้เห็นถึงสิ่งที่ฉันทำ): Colormap ไม่ปรากฏขึ้นและไม่ให้ฉันนำเข้า Colormap ของฉันจาก tiffs: มีวิธีการเพิ่ม / นำเข้า Colormaps ของฉันเพื่อให้เราสามารถดูข้อมูลตามค่าเหล่านั้นหรือไม่

9 arcgis-desktop  raster  arcmap  mosaic-dataset  colormap 

ฉันกำลังใช้levelplotฟังก์ชั่นของrasterVisแพ็คเกจสุดเจ๋งเพื่อสร้างแผนที่ที่มีค่าต่างกันประมาณศูนย์ ฉันกำลังวางแผนไฟล์แรสเตอร์โดยใช้จานสีแดงเป็นสีน้ำเงิน แต่ฉันพยายามกำหนดค่าสีเทาให้เป็นศูนย์ในแผนที่ โดยเฉพาะฉันต้องการทำซ้ำสีของรูปนี้: ขอให้สังเกตว่ามีการใช้ระดับสีแดงถึงสีน้ำเงิน แต่ค่าศูนย์มีการใช้สีเทา ขณะนี้แผนที่ของฉันมีลักษณะดังนี้: และนี่คือรหัสที่จะทำซ้ำ (ไฟล์มีอยู่ที่https://www.dropbox.com/s/cypfdu1eaz2fuok/r.annual.tif?dl=0 ): # Load required packages library(rasterVis) # open file r.annual <- raster("Downloads/r.annual.tif") # Set color palette myTheme=rasterTheme(region=brewer.pal('RdBu', n=11)) # Plot levelplot(annual.mask, par.settings=myTheme, margin=F) ฉันจะกำหนดสีเทาให้กับค่าศูนย์ทั้งหมดในแผนที่ของฉันด้านบนได้อย่างไร

9 raster  r  legend  visualisation  color-ramp 

ฉันกำลังแปลงเวกเตอร์เป็นแรสเตอร์ใน 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 

การทำงานกับ QGIS 2.4.0 Chugiak ฉันมีคำถามเกี่ยวกับการตัดภาพแรสเตอร์ ฉันมีภาพถ่ายทางอากาศที่ฉันต้องการคลิปโดยใช้ extents ของ shapefile จนถึงตอนนี้ฉันใช้เครื่องมือแรสเตอร์ / การสกัด / clipper เพื่อทำสิ่งนี้ ปัญหาคือขนาดไฟล์ของเอาท์พุทแรสเตอร์ของฉันมีขนาดใหญ่กว่าไฟล์ต้นฉบับเกือบสามเท่า (เช่น 725.283 และ 249.693 kb ตามลำดับ) ในขณะที่มันครอบคลุมพื้นผิวที่เล็กลง ทำไมไฟล์ของฉันถึงใหญ่กว่า และมีวิธีการ 'ตั้งค่า' ขนาดไฟล์ที่ส่งออกหรือไม่? ฟังก์ชั่นความช่วยเหลือของปัตตาเลี่ยนอ้างถึงหน้าต่อไปนี้: http://www.gdal.org/gdal_translate.htmlแม้ว่าฉันจะไม่ได้ประโยชน์อะไรมากมาย ฉันค่อนข้างใหม่กับ QGIS

9 qgis  raster  clip  file-size 

ฉันพยายามฉายภาพแรสเตอร์ ใน R มีprojectRaster()ฟังก์ชั่นสำหรับสิ่งนี้ (ด้านล่างเป็นตัวอย่างที่ทำซ้ำได้อย่างสมบูรณ์): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') ซึ่งใช้งานได้ดี อย่างไรก็ตามมันค่อนข้างช้า เพื่อเพิ่มความเร็วฉันแม้ว่าจะใช้gdalwarpแทน (ด้วย SSD ค่าใช้จ่ายในการอ่านและเขียนจาก / …

9 raster  gdal  r  coordinate-system  gdalwarp 

ฉันต้องการโหลดไฟล์. adm ลงในอาร์ข้อมูลมาจากหน้านี้: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 ฉันลองใช้รหัสต่อไปนี้ที่ฉันพบหลังจากการวิจัยในอินเทอร์เน็ต ปัญหาคือว่าในคลาส RasterLayer ฉันได้รับค่า negativ ที่ไม่ควรมี ฉันไม่รู้ว่าทำไมสิ่งนี้จึงเกิดขึ้นดังนั้นหวังว่าใครบางคนจะสามารถช่วยฉันได้! รหัส: library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r <- raster(xx) ผลลัพธ์สำหรับ 'r' คือ: คลาส r: RasterLayer ขนาด: 2160, 4320, 9331200 (nrow, ncol, ncell) ความละเอียด: 0.08333333, 0.08333333 (x, y) ขอบเขต: -180, 180, -90, 90 (xmin, …

9 raster  data  r  esri-grid-format 

ฉันต้องการนำเข้าแรสเตอร์โมเสคหลาย ๆ อัน (แต่ละโมเสคประกอบด้วย 7 ชั้น) จากโฟลเดอร์ที่ระบุลงในอาร์ของฉันจากนั้นเข้าถึงพวกมันเป็นแรสเตอร์แบบหลายเลเยอร์เดียว ดังนั้นฉันจึง: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported rasters p.ex. raster n°8 from "s" raster stack s[[8]] และแรสเตอร์ของฉัน [[8]] มีเพียงชั้นเดียวเท่านั้นดังนั้นจึงไม่มีการนำเข้าโมเสคทั้งหมด! nlayers(s[[8]]) [[1]] ถ้าฉันอ่านแต่ละโมเสกแยกกันมันใช้งานได้: # read 1 mosaic (composed …

9 raster  r  remote-sensing  landsat 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันกำลังสร้างสคริปต์ใน python โดยใช้ OGR / GDAL ฉันมีชุดไฟล์รูปร่างและชุดไฟล์แรสเตอร์ GeoTiff ฉันต้องการให้สคริปต์ของฉันเพิกเฉยกับรูปร่างของไฟล์หากพวกเขาไม่ตัดกันกับพื้นที่แรสเตอร์ Shapefile ไม่ได้เป็นรูปสี่เหลี่ยมผืนผ้าดังนั้นฉันไม่สามารถเปรียบเทียบค่า xmin / xmax, ymin / ymax ที่ส่งกลับโดย layer เท่านั้น GetGetExtent () ฉันต้องการรูปหลายเหลี่ยมที่เกิดขึ้นจริงซึ่งเป็นรูปร่างโดยรวมแล้ววิธีพิจารณาว่ารูปหลายเหลี่ยมนั้นตัดกับสี่เหลี่ยมแรสเตอร์หรือไม่ ฉันคิดว่าฉันสามารถรวมรูปหลายเหลี่ยมทั้งหมดใน shapefile เข้ากับคุณลักษณะหนึ่งแล้วอ่านเรขาคณิตในคุณลักษณะนั้นแล้วเปรียบเทียบข้อมูลนั้นกับขอบเขตแรสเตอร์ อย่างไรก็ตามฉันไม่แน่ใจเกี่ยวกับวิธีการดำเนินการนี้โดยเฉพาะ วิธีดึงข้อมูลรูปหลายเหลี่ยมเส้นขอบจาก shapefile? วิธีการตรวจสอบว่ารูปหลายเหลี่ยมนั้นตัดกันพื้นที่สี่เหลี่ยมที่กำหนดหรือไม่

9 python  raster  shapefile  gdal  ogr 

ฉันจะหาค่าสูงสุดของเซลล์ได้จากแรสเตอร์แบบเรียงซ้อน Rmax <- maxValue(RAD1998.all[[1]]) ทำงานได้ดี แต่ Rmax <- maxValue(RAD1998.all[[2]]) ให้ NA แน่นอนไม่ได้อยู่ในแรสเตอร์ซ้อน นี่คือรหัสของฉัน: RAD1998 <- raster(paste(getwd(), "/1998bil/1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil", sep = "")) list.ras <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = T, pattern = ".asc.bil")) RAD1998.all <- stack(list.ras)

9 raster  r 

มีเครื่องมือหรือวิธีการใดใน ArcMap ที่สามารถวัดพื้นที่ของแต่ละระดับของความหนาแน่นเคอร์เนลเอาท์พุทได้หรือไม่? โปรดทราบว่านี่เป็นเพียงพล็อตความหนาแน่นเคอร์เนล (ไม่ใช่ไฟล์รูปร่างหรือรูปหลายเหลี่ยม) การวิเคราะห์ความหนาแน่นของเคอร์เนล (สีเขียวเข้ม (10%) ถึงสีแดง (90%) (1) ฉันพยายามทำบางสิ่ง แต่ต้องการความแม่นยำ

9 arcgis-desktop  raster  arcmap  kernel-density 

ฉันรวบรวมและติดตั้งต้นแบบ QGIS ปัจจุบันตามปกติและฉันนำเข้าข้อมูลแรสเตอร์ลงในฐานข้อมูล 9.3 ของฉันที่มีส่วนขยาย postgis 2.0.2 เมื่อใดก็ตามที่ฉันพยายามนำเข้าข้อมูลแรสเตอร์เป็น qgis โดยใช้ตัวจัดการฐานข้อมูลฉันได้รับข้อความนี้: GDAL provider: Cannot get GDAL raster band: Raster layer: Provider is not valid (provider: gdal, URI: PG: dbname=gis host=127.0.0.1 user=xxx password=xxx port=5432 mode=2 schema=public table=elevation ในบรรทัดคำสั่งฉันได้รับ ERROR 1: Computed PostGIS Raster dimension is invalid. You've probably specified unappropriate resolution. ฉันไม่รู้ว่าจะค้นหาข้อผิดพลาดได้ที่ไหน …

9 qgis  postgis  raster  gdal  raster2pgsql 

ฉันมีสองชั้นที่ซ้อนทับกัน. cw ในพื้นที่เดียวกันแม้ว่าจะมีเลเยอร์อยู่ตรงกลางโดยมีขอบสีขาว เมื่อฉันตั้งค่าพื้นหลังเป็น 255,255,255 และกำหนดเป็นกลวงมันจะกลวงเฉพาะค่า RGB นั้นเฉพาะเมื่อในความเป็นจริงมีค่า "สีขาว" อื่น ๆ ในเลเยอร์ของฉันเช่น: 254254254; 253,253,253 .. ไปจนถึง 249,249,249 มีวิธีในการกำหนดค่าพื้นหลังของจอแสดงผลเพื่อรวมชุดของค่า RGB ที่ไม่ซ้ำกันหนึ่งชุดหรือไม่ ฉันรู้ว่าฉันสามารถแก้ไขปัญหานี้ได้ (การตัด, เครื่องคำนวณแรสเตอร์และอื่น ๆ ) แต่ฉันต้องการเก็บไฟล์. cw ต้นฉบับเพียงไฟล์เดียวโดยไม่สร้างรายการซ้ำดังนั้นฉันจึงต้องการใช้ Arcmap symbology หรือเครื่องมือสร้างภาพ สิ่งที่ฉันต้องการก็เหมือนกับCustom transparency optionsใน QGIS แต่สำหรับ ArcMap

9 raster  arcmap  arcgis-10.2  ecw  display 

ฉันกำลังพยายามปรับปรุงกระบวนการเวกเตอร์ / หลามที่ยุ่งยากในขณะนี้สำหรับแบบจำลองความเป็นอันตรายตามธรรมชาติ ในขณะนี้เรามีสคริปต์ที่มีความยาวซึ่งสร้างเส้นระยะทาง / ทิศทางจากจุดที่กำหนดเพื่อพิจารณา: ประเภทของรูปหลายเหลี่ยมที่มันตัดกัน (เช่นป่าหญ้าป่าพรุ ฯลฯ ) ระยะทางถึงรูปหลายเหลี่ยมนั้น จำนวนเส้นเหล่านี้ตัดกันรูปหลายเหลี่ยมเพื่อกำหนดว่า 'ล้อมรอบ' เป็นเท่าใด มีส่วนร่วมมากขึ้น แต่นั่นคือส่วนสำคัญของมัน ฉันพยายามที่จะหาวิธีในการปรับปรุงและตอนนี้กำลังนิ่งงันในส่วนที่ 3 ความคิดคือการตรวจสอบว่าจุดที่ล้อมรอบด้วยรูปหลายเหลี่ยมอย่างสมบูรณ์ภายใน 200m ดังนั้นในภาพที่แนบมาของฉันฉันต้องการให้จุด A ถูกทำเครื่องหมายว่ามีความเสี่ยงสูงกว่าจุด B เนื่องจากล้อมรอบด้วยรูปหลายเหลี่ยมของฉัน สิ่งนี้ซ้ำไปซ้ำมาประมาณ 13 ล้านจุดดังนั้นจึงไม่ใช่งานเล็ก ๆ และฉันอยากมีพื้นผิวที่จะได้รับค่ามากกว่าการใช้งานสคริปต์ของเรา ฉันคิดว่าจะต้องมีการเปลี่ยนแปลงของเครื่องมืออุทกวิทยาหรือเส้นทางต้นทุนในการทำเช่นนี้ แต่ฉันไม่สามารถไปรอบ ๆ มันได้ ฉันจะไปเกี่ยวกับเรื่องนี้ได้อย่างไร

9 arcgis-10.0  arcgis-desktop  raster  arcpy 

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

9 qgis  raster  vector  elevation  transparency