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

ฉันกำลังใช้ QGIS ฉันต้องการคลิปเลเยอร์แรสเตอร์แรสเตอร์โดยใช้เลเยอร์ขอบเขตผู้ดูแลระบบที่เป็นข้อมูลเวกเตอร์ อย่างไรก็ตามเครื่องมือการประมวลผลทางภูมิศาสตร์ดูเหมือนจะใช้งานได้กับข้อมูลเวกเตอร์เท่านั้น ฉันจะคลิปชั้นฝนนี้ได้อย่างไร

58 qgis  raster  vector  clip 

ข้อมูลแรสเตอร์และเวคเตอร์ในบริบท GIS คืออะไร ในแง่ทั่วไปการใช้งานกระบวนการหรือการวิเคราะห์แต่ละประเภทนั้นเหมาะสมกับอะไร? (และไม่เหมาะสำหรับ!) มีใครบ้างที่มีรูปภาพขนาดเล็กรัดกุมและมีประสิทธิภาพซึ่งสื่อความหมายและความเปรียบต่างของข้อมูลพื้นฐานทั้งสองนี้

47 raster  vector 

ฉันมี DEM ที่ฉันต้องการเรียบหรือพูดคุยทั่วไปเพื่อลบความสุดขั้วภูมิประเทศ (ตัดออกจากยอดเขาและหุบเขาเติม) เป็นการดีที่ฉันต้องการควบคุมรัศมีหรือระดับของ "ความพร่ามัว" ในท้ายที่สุดฉันจะต้องมีชุดแรสเตอร์ที่มีตั้งแต่เล็กน้อยเบลอไปจนถึงเบลอจริงๆ (ในทางทฤษฎีความเบลอจะเป็น raster คงที่ของค่าเฉลี่ยเลขคณิตของค่าทั้งหมด) มีเครื่องมือหรือวิธีการใดบ้างที่ฉันสามารถใช้ได้ (อิงจาก Esri, GDAL, GRASS) ฉันต้องกลับบ้านอบรูสัวแบบเกาส์ของตัวเองหรือไม่? ฉันสามารถใช้ฟิลเตอร์ low-pass (เช่นฟิลเตอร์ของ ArcGIS ) ได้หรือไม่และถ้าเป็นเช่นนั้นฉันต้องรันหลายครั้งเพื่อให้ได้ผลของรัศมีขนาดใหญ่หรือไม่?

46 raster  tools  generalization  smoothing 

GDAL รวมถึงวิธีการ resampling นอกเหนือจากการผสมปกติของเพื่อนบ้านที่ใกล้ที่สุด, bilinear, ลูกบาศก์และ splines: "การ resampling แบบสุ่มของ Lanczos แบบมีหน้าต่าง" ฉันเข้าใจว่ามันเป็นตัวกรองคอนโวลูน แต่แตกต่างจากภาพที่ผลลัพธ์มักจะเป็นอัตนัยการ resampling ที่ใช้สำหรับข้อมูลเชิงพื้นที่มีความหมายอื่น Lanczos คืออะไรและการใช้งานมีผลกับเอาต์พุตอย่างไร

46 raster  gdal  resampling 

มีกฎที่ยากและรวดเร็วเกี่ยวกับวิธีการแก้ไขที่เหมาะกับข้อมูลแรสเตอร์แต่ละประเภทหรือไม่?

40 raster  interpolation  gis-principle  resampling 

ฉันมีชุดข้อมูลของค่าเหนือกริดกม. ในทวีปอเมริกาคอลัมน์คือ "ละติจูด", "ลองจิจูด" และ "การสังเกต" เช่น: "lat" "lon" "yield" 25.567 -120.347 3.6 25.832 -120.400 2.6 26.097 -120.454 3.4 26.363 -120.508 3.1 26.630 -120.562 4.4 หรือในฐานะกรอบข้อมูล R: mydata <- structure(list(lat = c(25.567, 25.832, 26.097, 26.363, 26.63), lon = c(-120.347, -120.4, -120.454, -120.508, -120.562), yield = c(3.6, 2.6, 3.4, 3.1, 4.4)), …

40 raster  convert  r 

ฉันมีอนุกรมเวลาของแรสเตอร์ที่เซลล์ของแรสเตอร์แต่ละตัวแทนค่าบางอย่างในเวลาหนึ่ง ฉันต้องการสร้างแผนที่ที่เน้นขนาดของแนวโน้มเชิงบวกหรือเชิงลบในค่านี้เมื่อเวลาผ่านไปสำหรับแต่ละเซลล์ วิธีการที่ค่อนข้างไร้เดียงสาของฉันคือการปรับให้พอดีกับการถดถอยเชิงเส้นอย่างง่าย (X = เวลาและ Y = ค่า) กับแต่ละเซลล์และส่งอาร์เรย์ของความลาดเอียงไปยังแรสเตอร์ (ตามภาพตัวอย่างด้านล่าง) สิ่งนี้สามารถกรองได้โดยการส่งออกผลลัพธ์ที่สำคัญเท่านั้น ฉันจะแสดงแนวโน้มในช่วงเวลาอื่นได้อย่างไร ฉันสนใจในเทคนิคทั่วไปไม่ใช่คำแนะนำเฉพาะซอฟต์แวร์

34 raster  time  gis-principle 

ฉันมีภาพถ่ายดาวเทียมประมาณ 1,000 ภาพในรูปแบบ TIFF และฉันต้องการสร้างไฟล์รูปร่างที่จะทำหน้าที่เป็นดัชนีของแรสเตอร์ นี่คือสิ่งที่คล้ายกับแคตตาล็อกแรสเตอร์ แต่ฉันไม่ต้องการที่จะสร้างแคตตาล็อกแรสเตอร์ อุปสรรคบางอย่างที่ฉันสามารถมองข้ามได้ก็คือภาพมีตำแหน่งทางภูมิศาสตร์ดังนั้นพวกเขาจึงไม่ได้เป็นรูปสี่เหลี่ยมผืนผ้า (ฉันกำลังพูดถึงพื้นที่ข้อมูล) เพื่อความกระจ่างแจ้งฉันต้องการรูปหลายเหลี่ยมเพื่อครอบคลุมเฉพาะพิกเซลที่ไม่ใช่ศูนย์ (หรือไม่มีข้อมูล) ของแรสเตอร์และไม่ใช่แรสเตอร์รูปสี่เหลี่ยมผืนผ้าทั้งหมด คำตอบส่วนใหญ่แล้วให้รูปหลายเหลี่ยมเป็นรูปสี่เหลี่ยมผืนผ้าซึ่งครอบคลุมข้อมูลเช่นเดียวกับพิกเซลที่ไม่ใช่ข้อมูล รูปภาพของฉัน ผลลัพธ์ที่ได้รับจากเครื่องมือที่ฉันตรวจสอบ (เช่น raster catlog, Arcscripts ต่างๆ, สคริปต์ Python ที่กำหนดเองที่กำหนดในคำตอบข้อใดข้อหนึ่ง): ผลลัพธ์ที่ฉันต้องการ:

33 qgis  arcgis-desktop  raster  arcgis-10.0  satellite 

ไม่มีใครมีคำแนะนำสำหรับกระบวนการทำงานสำหรับการถ่ายภาพทางอากาศที่มีความละเอียดสูงโดยใช้เครื่องมือโอเพนซอร์สหรือไม่? เปรียบเทียบความเป็นมิตรกับผู้ใช้และความทนทานระหว่าง Quantum GIS กับ GDAL?

33 raster  open-source-gis  georeferencing  imagery 

ฉันกำลังสร้างแผนที่สำหรับภาคตะวันออกเฉียงเหนือของสหรัฐอเมริกาพื้นหลังของแผนที่จำเป็นต้องเป็นแผนที่ระดับความสูงหรือแผนที่อุณหภูมิเฉลี่ยรายปี ฉันมีแรสเตอร์สองตัวจาก Worldclim.org ซึ่งให้ตัวแปรเหล่านี้แก่ฉัน แต่ฉันต้องคลิปพวกมันให้มากเท่าที่ฉันสนใจข้อเสนอแนะใด ๆ เกี่ยวกับวิธีการทำเช่นนี้ นี่คือสิ่งที่ฉันมี: #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj) #load data state<- data (stateMapEnv) elevation<-raster("alt.bil") meantemp<-raster ("bio_1.asc") #build the raw map nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut", "rhode island","new york","pennsylvania", "new jersey", "maryland", "delaware", "virginia", "west virginia") map(database="state", …

32 raster  r  clip  maps 

ฉันกำลังมองหาวิธีที่จะแทรกคำอธิบายสีอย่างต่อเนื่องสำหรับเลเยอร์แรสเตอร์ ฉันจัดการเพื่อสร้างสัญลักษณ์ที่ไม่ต่อเนื่องพร้อมค่าต่าง ๆ ตามแถบ แต่ฉันต้องการที่จะมีแถบต่อเนื่องที่มีค่าบางอย่างบนแถบคล้ายกับปลั๊กอินเก่า ๆ "One-Band Raster" ซึ่งมันไม่ทำงานกับรุ่น QGIS สูงกว่า 1.7 ฉันพยายามตั้งค่าในคุณสมบัติ -> สไตล์ซึ่งฉันเลือกประเภทเรนเดอร์ "Singleband pseudocolor" และดูกล่อง "สร้างทางลาดสีใหม่" ฉันลองทั้ง "ต่อเนื่อง" หรือ "ช่วงเวลาที่เท่าเทียมกัน" อย่างไรก็ตามใน Composer คำอธิบายจะแสดงด้วยสัญลักษณ์แบบไม่ต่อเนื่องเสมอ มีวิธีใดบ้างที่จะทำให้นักแต่งเพลง QGIS แสดงวงดนตรีต่อเนื่องที่ดีสำหรับเลเยอร์แรสเตอร์?

32 qgis  raster  legend  print-composer 

ฉันมีข้อมูล Raster สำหรับรัฐที่มีข้อมูลระดับความสูงเป็นพิกเซล ฉันยังมีเลเยอร์เวกเตอร์ซึ่งมีหน่วยงาน 300 รัฐ สิ่งที่ฉันต้องการคือฉันต้องการรับค่าระดับความสูงเฉลี่ยของแผนกเฉพาะจากข้อมูลแรสเตอร์โดยใช้ ArcGIS ฉันลองแยกโดยใช้หน้ากากใน ArcGIS เพื่อแยกโพลีกอนแต่ละอันเพื่อคลิปแรสเตอร์ลงใน 300 แรสเตอร์จากนั้นแต่ละแรสเตอร์ฉันได้รวบรวมระดับความสูงเฉลี่ย แต่ใช้เวลาหลายวันร่วมกันเพื่อทำงานให้เสร็จสมบูรณ์ มีวิธีอื่นง่าย ๆ ที่จะได้รับมัน

31 arcgis-desktop  raster  vector 

ฉันมีโมเดลการยกระดับใน R r <- raster("example.dem") plot(r, col = topo.colors(20)) ค่าบางค่าอยู่ต่ำกว่าระดับน้ำทะเล (0) และฉันต้องการให้ค่าเหล่านั้นแสดงเป็นสีแดง ฉันจะกำหนดช่วงเฉพาะให้กับสีเฉพาะในพล็อต () ได้อย่างไร

30 raster  r  legend  color 

ฉันใหม่กับ GIS ฉันมีรหัสบางอย่างที่แปลงภาพอินฟราเรดของดาวอังคารเป็นแผนที่ความเฉื่อยทางความร้อนซึ่งจะถูกเก็บไว้เป็นอาร์เรย์แบบสองมิติ ฉันบันทึกแผนที่เหล่านี้เป็นไฟล์ hdf5 แล้ว แต่ฉันต้องการบันทึกเป็นภาพแรสเตอร์จริงๆเพื่อให้สามารถประมวลผลได้ใน QGIS ฉันผ่านการค้นหาหลายครั้งเพื่อค้นหาวิธีการทำ แต่ไม่มีโชค ฉันได้ลองทำตามคำแนะนำในบทช่วยสอนที่http://www.gis.usu.edu/~chrisg/python/แต่ไฟล์ที่ฉันสร้างโดยใช้รหัสตัวอย่างของเขาเปิดเป็นกล่องสีเทาธรรมดาเมื่อฉันนำเข้า QGIS ฉันรู้สึกว่าถ้าใครบางคนสามารถแนะนำขั้นตอนที่ง่ายที่สุดเท่าที่เป็นไปได้สำหรับตัวอย่างที่เรียบง่ายของสิ่งที่ฉันต้องการจะทำแล้วฉันอาจจะสามารถก้าวหน้าได้บ้าง ฉันมี QGIS และ GDAL ฉันมีความสุขมากที่ได้ติดตั้งเฟรมเวิร์กอื่น ๆ ที่ทุกคนสามารถแนะนำได้ ฉันใช้ Mac OS 10.7 ดังนั้นถ้าเช่นฉันมีอาร์เรย์ของความเฉื่อยความร้อนที่ดูเหมือน: TI = ( (0.1, 0.2, 0.3, 0.4), (0.2, 0.3, 0.4, 0.5), (0.3, 0.4, 0.5, 0.6), (0.4, 0.5, 0.6, 0.7) ) และสำหรับแต่ละพิกเซลฉันมีละติจูดและลองจิจูด: lat = ( (10.0, …

30 python  raster  gdal  numpy 

ฉันกำลังแยกพื้นที่และร้อยละของการใช้ที่ดินประเภทต่าง ๆ จากแรสเตอร์ตามขอบเขตหลายเหลี่ยมหลายพัน ฉันพบว่าฟังก์ชั่นการดึงข้อมูลทำงานได้เร็วขึ้นมากถ้าฉันวนซ้ำรูปหลายเหลี่ยมและครอบตัดแต่ละอันจากนั้นปิดบังแรสเตอร์ให้มีขนาดเท่ากับรูปหลายเหลี่ยมที่เฉพาะเจาะจง อย่างไรก็ตามมันค่อนข้างช้าและฉันสงสัยว่าใครมีคำแนะนำในการปรับปรุงประสิทธิภาพและความเร็วของรหัสของฉันหรือไม่ สิ่งเดียวที่ฉันได้พบที่เกี่ยวข้องกับการนี้คือการตอบสนองนี้โดยโรเจอร์ Bivand ที่แนะนำให้ใช้GDAL.open()และGDAL.close()เช่นเดียวกับและgetRasterTable() getRasterData()ฉันดูสิ่งเหล่านี้ แต่เคยมีปัญหากับ gdal ในอดีตและไม่รู้จักดีพอที่จะรู้วิธีนำไปใช้ ตัวอย่างที่ทำซ้ำได้: library(maptools) ## For wrld_simpl library(raster) ## Example SpatialPolygonsDataFrame data(wrld_simpl) #polygon of world countries bound <- wrld_simpl[1:25,] #name it this to subset to 25 countries and because my loop is set up with that variable ## Example RasterLayer …

29 raster  r  performance  masking  clip