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

ฉันจะแยกรูทลูกบาศก์ของการคูณแรสเตอร์ 3 ตัวในเครื่องคิดเลขแรสเตอร์ได้อย่างไร เช่น cuberoot (ras1 * ras2 * ras3)

10 arcgis-desktop  raster  spatial-analyst  raster-calculator 

ฉันสับสนเล็กน้อยเกี่ยวกับการทำงานของเครื่องมือสถิติโซนใน QGIS (ทั้งเครื่องมือ 'สถิติสถิติ' ภายใต้เครื่องมือ Raster จากกล่องเครื่องมือทางภูมิศาสตร์ของ QGIS หรือเครื่องมือ 'สถิติ Raster จากรูปหลายเหลี่ยม' ภายใต้ SAGA Vector <-> เครื่องมือแรสเตอร์ ) โดยทั่วไปฉันมีเลเยอร์แผนที่แรสเตอร์และรูปหลายเหลี่ยมขนาดเล็กในชั้นแยกต่างหาก ฉันต้องการกำหนดค่าของรูปหลายเหลี่ยมด้านล่างแต่ละรูปหลายเหลี่ยม ในกรณีที่รูปหลายเหลี่ยมซ้อนกันหลาย ๆ เซลล์แรสเตอร์ฉันแค่ต้องการให้มีค่าสูงสุด แต่ผลลัพธ์ไม่สมเหตุสมผลและฉันไม่สามารถเข้าใจเหตุผลของผลลัพธ์เหล่านั้นได้ (ซึ่งเห็นได้ชัดว่าไม่ใช่จำนวนสูงสุดของเซลล์แรสเตอร์แบบขีดเส้นใต้): ความคิดใด ๆ ที่ฉันทำผิดหรือสิ่งที่ฉันยังไม่ได้พิจารณา? ผลกับ API หลามในqgis.analysis.QgsZonalStatistics()( ตามคำถามนี้ที่นี่ ) ที่ดีมาก แต่จะมีเพียงcount, meanและsum... แต่ด้วยเหตุผลบางอย่างไม่มีหรือmaxmin

10 qgis  raster  layers  saga  zonal-statistics 

ฉันต้องการใช้ไดอะล็อก QGIS Raster Style เพื่อแสดงค่าแรสเตอร์ที่ดัดแปลง โดยเฉพาะฉันต้องการเห็นการไล่ระดับสีโดยใช้สแควร์รูทของมูลค่าของแรสเตอร์ ฉันเข้าใจว่าเป็นไปได้ที่จะใช้เครื่องคำนวณแรสเตอร์เพื่อแปลงค่าแรสเตอร์ แต่ถ้าเป็นไปได้ฉันอยากข้ามขั้นตอนนั้นและใช้สูตรทางคณิตศาสตร์ในระหว่างการจัดแต่งทรงผมและหลีกเลี่ยงการสร้างสำเนาแรสเตอร์ดั้งเดิม การทำงานประเภทนี้สามารถทำได้โดยใช้ไฟล์เวกเตอร์โดยใช้ค่าที่กำหนดไว้สำหรับสี / ความกว้าง / ฯลฯ เป็นไปได้ไหมที่จะทำเช่นนี้กับ rasters?

10 qgis  raster  style 

ฉันได้นำเข้า rasters บางอย่างที่ควรจะเชื่อมต่ออย่างต่อเนื่อง แต่พวกเขาสร้างตารางที่มีขอบกั้นระหว่างพวกเขา ฉันอาจสันนิษฐานว่านี่เป็นเรื่องของการสอบเทียบระหว่างแรสเตอร์แต่ละอัน แต่ฉันไม่แน่ใจว่าจะรันมันอย่างไรในรูปแบบดาต้า มีใครเคยมีประสบการณ์เช่นนี้มาก่อนและอาจจะรู้วิธีแก้ปัญหาที่เป็นไปได้เพื่อให้ภาพรวมของ rasters ทั้งหมดจะราบรื่นขึ้น ทั้งหมดนี้ควรจะเชื่อมโยงกันในเชิงภูมิศาสตร์ ดังนั้นปัญหาของ 'ขอบเขต' ของสแควร์ ฉันใช้ QGIS บน windows

10 qgis  raster  windows  smoothing 

ฉันกำลังพยายามจัดแนวแรสเตอร์สองกริดในอาร์เมื่อจัดตำแหน่งแล้วฉันต้องการที่จะเพิ่มพวกเขาเข้าด้วยกัน ฉันได้ลองตรวจสอบว่าการทำแบบstackนี้จะได้ผลหรือไม่: grid_snap <- stack(habi_sdw, Pop_sdw) และฉันได้รับข้อผิดพลาดต่อไปนี้: ข้อผิดพลาดใน comparRaster (x): ขอบเขตที่ต่างกัน ตะแกรงแรสเตอร์มีคุณสมบัติดังต่อไปนี้: show(habi_sdw) # class : RasterLayer # dimensions : 9187, 9717, 89270079 (nrow, ncol, ncell) # resolution : 0.00892857, 0.00892857 (x, y) # extent : -28.83706, 57.92186, -36.02464, 46.00214 (xmin, xmax, ymin, ymax) # coord. ref. : +proj=longlat +ellps=WGS84 …

10 raster  r  clip  snapping  resolution 

ฉันมี orthoimage RGB ที่มีพื้นที่สีดำแทน NoData ดังภาพที่นี่: เมื่อฉันใช้เครื่องมือ SetNull มันจะลบพื้นที่สีดำออกไป แต่มันจะเปลี่ยนเป็นสีเทาแบบออร์โธไมซ์ซึ่งไม่ใช่สิ่งที่ฉันต้องการ ฉันรู้ว่ามีวิธีแก้ไขปัญหานี้โดยใช้เครื่องมือ Clip (Data Management) เพื่อคลิปแรสเตอร์ให้กับตัวเองและในนั้นคุณสามารถตั้งค่า NoData เป็นค่าของเซลล์สีดำ (0) อย่างไรก็ตามฉันกำลังสร้างเวิร์กช็อปสำหรับนักเรียนและฉันต้องการใช้วิธีการที่ใช้งานง่ายกว่า แน่นอนว่ามีตัวเลือกในการเรียกใช้เครื่องมือ NoData ในแต่ละแบนด์แยกจากกันแล้วรวมสามแบนด์เข้าด้วยกัน แต่ฉันต้องมีเงื่อนไขที่ band1, band2 และ band3 เท่ากับ 0 ทั้งหมดเพื่อตั้งเป็น NoData ความสามารถในการทำเช่นนี้ใน ModelBuilder จะเหมาะ แต่การเขียนสคริปต์โดยใช้ ArcPy ก็จะดีเช่นกัน ฉันกำลังใช้เดสก์ท็อป ArcGIS 10.2 พร้อมใบอนุญาตขั้นสูง

10 arcgis-desktop  raster  arcgis-10.2  nodata  multi-band 

ฉันได้ให้ภาพ Geotiff และข้อมูล Lidar ที่สอดคล้องกัน (x, y, z) ในพิกัด UTM ฉันต้องรวมข้อมูล Lidar เข้ากับค่า RGB จากภาพ นั่นหมายความว่าในตอนท้ายฉันต้องพล็อต (3D) แต่ละจุดของสีคลาวด์ LiDAR ที่เข้ารหัสด้วยค่า RGB ที่สอดคล้องกันจากภาพ Geotiff ฉันแปลงข้อมูล Lidar เป็น shapefile โดยใช้ QGIS ฉันควรทำอย่างไรต่อไป ใน R ฉันลองใช้plot3Dฟังก์ชั่น แต่มันไม่ทำงาน ฉันกำลังแนบdoc text , shapefile และอิมเมจ tif แก้ไข: ฉันได้ทำโปรแกรมต่อไปนี้ตามที่แสดงด้านล่าง: require(raster) require(maptools) # to take shape files #require(car) …

10 raster  r  3d  lidar  3d-model 

ฉันได้ค้นพบเพียงแค่การจัดหมวดหมู่ Hansen's น่าทึ่งเกี่ยวกับการสูญเสียป่าไม้และเรื่องการเจริญเติบโตได้ที่http://www.globalforestwatch.org/ตีพิมพ์ในวิทยาศาสตร์ 2013เมื่อ: แฮนเซน, MC, Potapov พีวีมัวร์ R ., Hancher, M. , Turubanova, SA, & Tyukavina, A. (2013) แผนที่ความละเอียดสูงระดับโลกของการเปลี่ยนแปลงการปกคลุมป่าในศตวรรษที่ 21 วิทยาศาสตร์, 342 (6160) (15 พฤศจิกายน), 850–854 ดอย: DOI: 10.1126 / science.1244693 อย่างไรก็ตามฉันไม่สามารถหาวิธีการที่แน่นอนในการทำซ้ำแผนที่ดังกล่าวในเว็บไซต์นี้ / บนเว็บไซต์เพื่อให้การจำแนกประเภทที่ Hansen ใช้? สิ่งเดียวที่ฉันสามารถหาได้คืออัลกอริทึมการเรียนรู้ภายใต้การดูแลใช้เพื่อระบุต้นไม้ปกคลุม แต่มันค่อนข้างกว้าง ถ้าเป็นไปได้ฉันต้องการใช้วิธีการเดียวกัน (แต่ใช้กับ 90 ปี) ดังนั้นก่อนการจัดประเภทของ Hansen ในพื้นที่ที่ฉันเลือก

10 raster  remote-sensing  classification  landsat  forest-ecology 

ฉันมีชุด rasters ในรูปแบบ DTED ซึ่งฉันโหลดลงในฐานข้อมูล PostGIS โดยใช้เครื่องมือบรรทัดคำสั่ง raster2pgsql แรสเตอร์แต่ละตัวจะถูกจัดเก็บไว้ในแถวและอธิบายด้วยค่าลบและค่าของรูปแบบแรสเตอร์ ตอนนี้ฉันต้องการสร้างฟังก์ชั่นฐานข้อมูลที่ใช้ลองจิจูดและละติจูดของจุดและส่งกลับค่าของพิกเซลที่สอดคล้องกับจุดนี้ ปัญหาที่ฉันมีคือมันใช้เวลาค่อนข้างมากในการทำเช่นนั้น (3-4 วินาที) เพราะฐานข้อมูลทำงานบนบอร์ด Odroid ฉันรู้ว่าชุดข้อมูลที่ฉันดำเนินการค่อนข้างใหญ่ (rasters ครอบคลุมพื้นที่ทั้งหมดในสหราชอาณาจักร) แต่เนื่องจากฉันไม่คุ้นเคยกับ PostgreSQL และ PostGIS มากฉันสงสัยว่ามันสามารถทำได้เร็วกว่า นี่คือสิ่งที่ฉันทำไปแล้ว: SELECT ST_Value(rast, ST_GeomFromText(CONCAT('POINT(', $1, ' ', $2, ')'), 4326)) FROM ( SELECT * FROM rasters WHERE rast && ST_GeomFromText(CONCAT('POINT(', $1, ' ', $2, ')'), 4326) ) …

10 postgis  raster  query  optimization 

ฉันต้องการที่จะ georeference แรสเตอร์ที่ใช้และpython GDALแนวทางปัจจุบันของฉันคือการโทรgdal_translateและgdalwarpใช้os.systemและรายการจุดควบคุมภาคพื้นดินที่น่าเกลียด ฉันต้องการวิธีการทำสิ่งนี้ภายในpythonจริงๆ นี่เป็นกระบวนการปัจจุบันที่ฉันใช้: import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"') os.system('gdalwarp -r bilinear -tps -co COMPRESS=NONE "outraster2.tif" "outraster3.tif"') มีก่อนหน้านี้คือคำถามและคำตอบจาก 2012 ซึ่งระบุว่าสามารถเข้าถึงได้หลังจากการนำเข้าgdal_translate gdalฉันไม่แน่ใจว่าล้าสมัยหรือว่าผิด แต่เมื่อฉันรันfrom osgeo import gdalฉันไม่เห็นgdal.gdal_translateตัวเลือก ฉันไม่รู้ว่ามันมีอยู่หรือไม่ แต่ฉันจะรักถ้าฉันสามารถแปลและปฏิเสธ rasters ในแบบ pythonic ตัวอย่างเช่น: # translate gcp_points = [(1251.92, …

10 python  raster  gdal  coordinate-system  gdalwarp 

มีวิธีตรวจสอบเพื่อดูว่าเลเยอร์แรสเตอร์ที่กำหนด 2 รายการมีเนื้อหาเหมือนกันหรือไม่? เรามีปัญหาเกี่ยวกับปริมาณพื้นที่เก็บข้อมูลที่ใช้ร่วมกันในองค์กรของเรา: ขณะนี้มีขนาดใหญ่มากซึ่งใช้เวลามากกว่า 3 วันในการสำรองข้อมูลเต็มรูปแบบ การตรวจสอบเบื้องต้นพบว่าหนึ่งในผู้ร้ายที่กินพื้นที่มากที่สุดคือเปิด / ปิดแรสเตอร์ที่ควรเก็บเป็นเลเยอร์ 1 บิตด้วยการบีบอัด CCITT อิมเมจตัวอย่างนี้เป็น 2 บิตในปัจจุบัน (ดังนั้นค่าที่เป็นไปได้ 3 ค่า) และบันทึกเป็นไฟล์บีบอัด LZW ขนาด 11 MB ในระบบไฟล์ หลังจากแปลงเป็น 1 บิต (เพื่อให้มี 2 ค่าที่เป็นไปได้) และใช้การบีบอัด CCITT Group 4 เราได้รับการลดลงเหลือ 1.3 MB เกือบเต็มตามลำดับความสำคัญของการออม (อันที่จริงแล้วเป็นพลเมืองที่ประพฤติตัวดีมีคนอื่นเก็บไว้เป็น 32 บิตลอย!) นี่เป็นข่าวที่น่าอัศจรรย์! อย่างไรก็ตามมีเกือบ 7,000 ภาพที่จะใช้เช่นกัน มันจะตรงไปตรงมาในการเขียนสคริปต์เพื่อบีบอัดพวกเขา: for old_img in …

10 python  raster  compression  differences 

วิธีที่ง่ายที่สุดในการ จำกัด หน่วยความจำที่rasterแพ็คเกจใช้ในRRAM 5GB คืออะไร เป็นการดีที่ฉันต้องการแรสเตอร์ใด ๆ ที่มีขนาดใหญ่กว่านี้เพื่อเขียนลงในฮาร์ดไดรฟ์และวิเคราะห์จากไฟล์ชั่วคราว ฉันใช้ Windows 7 (64Bit) และรุ่น R 3.0.3 นี่คือรหัสบางส่วนที่ฉันเล่นด้วยบนเครื่องที่มี RAM ขนาด 16GB แต่มันก็ยังขยายออกไปเรื่อย ๆ เมื่อฉันจัดการกับ rasters ขนาดใหญ่: require(raster) memory.limit(size = 5000) rasterOptions(format="GTiff",overwrite=TRUE,datatype="INT1S", tmpdir="C:/Research/BIN",tmptime=1.1,progress="text",chunksize=1000, maxmemory=1000) rasterTmpFile("delete_me_") r <- raster("myraster.tif") r[r==0] <- NA มันอยู่ในขั้นตอนนี้ในสคริปต์ที่การใช้ RAM เพิ่มขึ้นเพื่อให้ 16GB ทั้งหมดแล้วหยุดทำงานและยกเลิกฟังก์ชั่น ฉันต้องเปลี่ยนอะไรที่นี่

10 raster  r  memory 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน2 ปีที่ผ่านมา ฉันพยายามรวม 14 geotiff แบบนี้: geotiff แต่ละอันมีค่าประมาณ 50Mb ฉันต้องการ geotiff ที่เอาต์พุต เวิร์กโฟลว์ของฉัน: gdalbuildvrt -input_file_list list.txt test.vrt (ที่รายการของฉันมีชื่อของ tifs) จากนั้น: gdal_translate -of Gtiff test.vrt test.tif Input file size is 79841, 59955 มันใช้งานได้ แต่ผลที่ได้คือ geotiff ที่ 13,3 Gb! สำหรับไฟล์ 14 ไฟล์แต่ละ 50 Mb ฉันพยายาม Geotiff ที่ …

10 raster  gdal  geotiff-tiff  gdal-translate  gdalbuildvrt 

ฉันพยายามที่จะหลายเชิงภาพโมเสคเข้าไปในแรสเตอร์ขนาดใหญ่เดียวในอาร์การใช้สคริปต์ที่โพสต์ใน/programming/15287807/how-can-i-create-raster-mosaic-using-list-of-rasters แต่ฉันได้รับข้อความเตือนและข้อความแสดงข้อผิดพลาด rasters1 <- list.files("F:\\MOD15A2_LAI_1km\\MOD15A2_LAI_2009", pattern = "mod15a2.a2009001.*.005.*.img$", full.names = TRUE, recursive = TRUE) mos1 <-mosaic(rasters1, fun=mean) มันเป็นรายงานข้อผิดพลาดดังต่อไปนี้ Error in (function (classes, fdef, mtable) : unable to find an inherited method for function ‘extent’ for signature ‘"character"’ จากนั้นฉันลองรุ่นอื่น rasters1.mosaicargs <- rasters1 rasters1.mosaicargs$fun <- mean แต่นี่คือข้อความเตือนบางส่วนดังต่อไปนี้ Warning message: In rasters1.mosaicargs$fun <- …

10 raster  r  mosaic 

เมื่อประมวลผลข้อมูล DEM ขนาด 30 ม. และ 10 ม. ที่ดาวน์โหลดผ่าน National Map Viewer จากชุดข้อมูลระดับชาติเราสังเกตเห็นการสตริปในแนวนอนและแนวตั้งไม่เพียง แต่ในผลการผลิต แต่ในการวิเคราะห์เฉดสี DEM ดิบเท่านั้น ไม่มีใครรู้แหล่งที่มา? หากไม่ใช่แหล่งที่มาอาจเป็นวิธีการลบสิ่งประดิษฐ์เหล่านี้หรือไม่ สิ่งประดิษฐ์เหล่านี้เด่นชัดมากเมื่อใช้ DEMs ในการคำนวณดัชนีภูมิประเทศ สิ่งประดิษฐ์เหล่านี้ยังคงอยู่แม้จะเกิดภาวะซึมเศร้าเกิดขึ้น ภาพด้านล่างแสดงภาพสตริปในข้อมูล 30 ม. และ 10 ม. จากแหล่งต้นน้ำทั้งในเพนซิลเวเนียและโคโลราโดและการคำนวณดัชนีภูมิประเทศเสร็จแล้วแสดงสิ่งประดิษฐ์สำหรับสันปันน้ำในซีราคิวส์นิวยอร์ก โคโลราโด - HUC8 - 10190004 - 10m โคโลราโด - HUC8 - 10190004 - 30m เพนซิลเวเนีย - HUC8 - 02040103 - …

10 arcgis-desktop  raster  dem  saga  united-states