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

ฉันจะหาค่าสูงสุดของเซลล์ได้จากแรสเตอร์แบบเรียงซ้อน 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 

ฉันได้คำนวณพื้นที่ผิวของการกระจายของสปีชีส์ (การรวมรูปหลายเหลี่ยมจากรูปร่างของไฟล์) แต่เนื่องจากพื้นที่นี้อาจประกอบด้วยรูปหลายเหลี่ยมที่ห่างไกลมากฉันจึงต้องการคำนวณการกระจายตัวของบางอย่าง สิ่งที่ฉันได้ทำไปแล้วคือการดึงเซนทรอยด์ของรูปหลายเหลี่ยมแต่ละอันคำนวณระยะห่างระหว่างพวกมันและใช้สิ่งเหล่านี้เพื่อคำนวณสัมประสิทธิ์การแปรผันตามตัวอย่างในตัวอย่างด้านล่าง require(sp) require(ggplot2) require(mapdata) require(gridExtra) require(scales) require(rgeos) require(spatstat) # Create the coordinates for 3 squares ls.coords <- list() ls.coords <- list() ls.coords[[1]] <- c(15.7, 42.3, # a list of coordinates 16.7, 42.3, 16.7, 41.6, 15.7, 41.6, 15.7, 42.3) ls.coords[[2]] <- ls.coords[[1]]+0.5 # use simple offset ls.coords[[3]] <- …

9 polygon  r  distance  centroids 

ทุกคนขอโทษที่รบกวน แต่ฉันค่อนข้างใหม่กับ r ต้องเผชิญกับ difficuty สำคัญ: ฉันต้องการสร้างแผนที่เคลื่อนไหวของ Russin กับการเปลี่ยนแปลงในการว่างงานกับปีที่แตกต่างเช่น ในภาพคุณสามารถดูข้อมูลเป็นเวลาหนึ่งปี require(sp) require(maptools) require(RColorBrewer) require(rgdal) rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData") print(load(rus)) unempl <- read.delim2(file="C:\\unempl1.txt", header = TRUE, sep = ";",quote = "", dec=",", stringsAsFactors=F) gadm_names <-gadm$NAME_1 total <- length(gadm_names) pb <- txtProgressBar(min = 0, max = total, style = 3) order <- vector() for (i …

9 cartography  r  animation 

ฉันพยายามลบจุดตัดออกจาก Poly2 ไปยัง Poly1 ด้านล่างนี้เป็น coords ของ 2 polygons > coords1 [,1] [,2] [1,] 992.0161 7462.531 [2,] 950.1962 7413.532 [3,] 902.3632 7421.207 [4,] 706.1985 7472.378 [5,] 654.0139 7711.760 [6,] 657.5960 7726.950 [7,] 786.2667 7844.380 [8,] 789.5935 7842.925 [9,] 981.7046 7498.659 [10,] 983.4246 7493.271 [11,] 990.6680 7469.888 [12,] 992.0161 7462.531 > …

9 r  polygon  intersection  symmetricdifference 

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

9 arcgis-desktop  python  arcpy  raster  r 

ภาพที่แนบมาแสดงให้เห็นช่องว่างของป่าที่มีต้นสนสีแดงแสดงเป็นวงกลมและต้นสนสีขาวที่แสดงเป็นรูปกากบาท ฉันสนใจที่จะพิจารณาว่ามีความสัมพันธ์เชิงบวกหรือเชิงลบระหว่างต้นสนสองสายพันธุ์หรือไม่ (เช่นมีการเติบโตในพื้นที่เดียวกันหรือไม่) ฉันตระหนักถึง Kcross และ Kmulti ในแพ็คเกจ R spatstat อย่างไรก็ตามเนื่องจากฉันมี 50 ช่องว่างในการวิเคราะห์และคุ้นเคยกับการเขียนโปรแกรมใน python มากกว่า R ฉันจึงต้องการหาวิธีการวนซ้ำโดยใช้ ArcGIS และ python ฉันยังเปิดให้โซลูชั่น R ฉันจะใช้ฟังก์ชั่น K ของริบลีย์ bivariate ได้อย่างไร?

9 arcgis-desktop  python  arcpy  r  spatial-statistics 

ในเอกสารของspฉันเห็นว่าการcoordsโต้แย้งของSpatialPolygonsฟังก์ชั่นยอมรับเพียง 2 คอลัมน์ ไม่สามารถสร้างรูปหลายเหลี่ยม 3 มิติได้หรือไม่ ฉันทำ kriging spatio กาลใช้แพคเกจโดยการเพิ่มเวลาเป็นมิติเชิงพื้นที่ในสามของgstat SpatialPointsDataFrameฉันต้องการที่จะขยายที่บล็อก kriging กับรูปหลายเหลี่ยมเป็นบล็อก มีวิธีใดที่จะทำเช่นนั้น?

9 r  interpolation  spatial-statistics  kriging  gstat 

ฉันต้องการ Spatial Grid เป็นกริดหลักสำหรับแผนที่เฉพาะเรื่องที่หลากหลาย ฉันจะสร้าง Spatial Grid จากแรสเตอร์เพื่อยกเลิกพิกเซล NA ทั้งหมดได้อย่างไร

9 r 

มี 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 

อ้างถึงคำถามก่อนหน้าของฉันฉันจะ "geocode" / แปลง - ใช้ R - ภาพแรสเตอร์ที่สร้างขึ้นในรูปแบบที่สามารถแสดงเป็นเลเยอร์บนแผนที่ได้อย่างไร ในระยะสั้นฉันตั้งใจจะสร้างผลลัพธ์ของซอฟต์แวร์ SpatialKey ( ภาพหน้าจอ ) อีกครั้ง ฉันใช้อินสแตนซ์ GeoServer เพื่อโฮสต์ข้อมูลสำหรับแผนที่ แต่ฉันจะรวมการฉายภาพและการประสานงานข้อมูลที่ต้องการลงในภาพได้อย่างไรเพื่อที่จะจัดตำแหน่งให้ถูกต้องกับตำแหน่งทางภูมิศาสตร์ที่แท้จริง Heatmap ความหนาแน่นถูกสร้างขึ้นโดยใช้ข้อมูลเชิงพื้นที่

9 raster  r  georeferencing 

อาจเป็นคำถามง่ายๆ ฉันมีรายการที่ตั้งทางภูมิศาสตร์ที่มีศักยภาพนับพัน (lat lon) และจากนั้นฉันต้องเลือก 200 แห่งที่ "แพร่กระจายมากที่สุด" ฉันคิดว่านั่นน่าจะเป็น 200 คะแนนที่มีระยะทางเฉลี่ยโดยรวมสูงสุด คิดถึงร้านค้าในเมือง มีวิธีการที่กำหนดไว้หรือไม่? บางทีในแพ็คเกจ R ขอบคุณทุกคนที่ทำให้ที่นี่เป็นสถานที่ที่ยอดเยี่ยมในการเรียนรู้! / คริส

9 r  spatial-statistics  business 

นี่คือcrosspost จาก StackOverflow ฉันรู้น้อยมากเกี่ยวกับซอฟต์แวร์ GIS และกำลังทำแผนที่ทั้งหมดของฉันใน R. Apologies ล่วงหน้าหากนี่เป็นคำถามพื้นฐานเกินไป สมมติว่าฉันมีสองรูปร่างไฟล์จากแหล่งที่แตกต่างกัน แต่มีคุณสมบัติที่แตกต่างกัน สมมติว่าข้อแรกมีไว้สำหรับขอบเขตการจัดการของเท็กซัส ( boundaries.shp) และอีกอันสำหรับแม่น้ำเท็กซัส ( rivers.shp) ฉันยังมีไฟล์ที่สามที่towns.csvแสดงที่ตั้งของเมืองในรัฐ หลังจากอ่านในไฟล์ฉันสามารถระบุตำแหน่งของเมืองมากเกินไปบนขอบเขตของการสั่งการในmaptoolsแพ็คเกจ: plot(boundaries); points(towns$lon, towns$lat) แต่ฉันจะซ้อนทับทั้งสามได้อย่างไร แน่นอนมีวิธีง่ายๆในการทำเช่นนี้?

9 r  overlay 

ฉันมีชุดข้อมูลในรูปแบบ csv เช่น: latitude, longitude, value -45, 45, 10 .... ฉันต้องการคำนวณค่าระดับเขตจากข้อมูลเหล่านี้โดยใช้หญ้า ArcGIS หรือ R + ฉันมีเลเยอร์ ArcGIS USA Counties แต่ฉันไม่สามารถทำการเข้าร่วมเชิงพื้นที่ให้สำเร็จได้ ฉันจะทำสิ่งนี้ได้อย่างไร

9 arcgis-desktop  data  geoprocessing  r