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

วิธีที่ง่ายที่สุดในการ จำกัด หน่วยความจำที่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 

ฉันพยายามที่จะหลายเชิงภาพโมเสคเข้าไปในแรสเตอร์ขนาดใหญ่เดียวในอาร์การใช้สคริปต์ที่โพสต์ใน/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 

ฉันพยายามที่จะวางทับตัวแปรที่ระบุ (2 หมวดหมู่ - rainfed เทียบกับชลประทาน) ในแผนที่ choropleth เหนือตัวแปรตัวเลข เพื่อเน้นความแตกต่างของลวดลายในสองหมวดหมู่ฉันได้ระบุหมวดหมู่หนึ่งที่มีสีแตกต่างกัน (พื้นที่ชลประทานใน "สีแดง") ในขณะที่หมวดที่ไม่มีสีคือพื้นที่ที่มีน้ำฝน อย่างไรก็ตามแผนที่รู้สึกว่ารกและโครงร่างผสมกับเฉดสีเข้มกว่าของตัวแปรตัวเลข ฉันพยายามหลีกเลี่ยงการใช้ชุดรูปแบบการแรเงาสีที่แตกต่างกันสองแบบเพื่อเป็นตัวแทนของทั้งสองประเภทเนื่องจากฉันรู้สึกว่าจะซ่อนการกระจายโดยรวมของตัวแปรพื้นฐาน มีวิธีที่ดีกว่าในการเป็นตัวแทนของการจำแนกประเภท bivariate นี้โดยใช้แพ็คเกจ ggplots2 ใน R (ซึ่งตอนนี้ฉันใช้เพื่อสร้างแผนที่ที่แนบมา)? ในที่สุดแผนที่เหล่านี้จะถูกพิมพ์ลงบนกระดาษ

10 r  cartography  ggplot2  choropleth 

ฉันกำลังทดลองกับไฟล์โมเดลภูมิประเทศดิจิตอลในรูปแบบแรสเตอร์กริดของ ESRI ASCII แม้จะไม่มีประสบการณ์กับข้อมูลประเภทนี้ฉันพบว่ามันง่ายมากในการโหลดลงใน R ตามโค้ดด้านล่าง plot3Dฟังก์ชั่นจากrasterVisค่าเริ่มต้นแพคเกจดีมองพล็อต แต่รูปแบบค่อนข้างสะท้อนแสง เนื่องจากมันควรจะเป็นแนวนอนพื้นผิวมันวาวนี้จึงไม่เหมาะสม ฉันกำลังมองหาธรรมชาติมากขึ้น (หากคำดังกล่าวสามารถใช้เกี่ยวกับการเรนเดอร์ 3 มิติของชุดตัวเลข), พล็อตเรื่อง, ทอพอโลยี ฉันคิดว่าจะต้องมีวิธีแก้ไขปัญหานี้ แต่ฉันไม่มีประสบการณ์ในการใช้rasterและค้นหาตัวชี้ ฉันจะลดหรือลบการสะท้อนของพล็อตได้อย่างไร library(raster) library(rgdal) library(rasterVis) foo <- raster(readGDAL("my.dtm.asc")) plot3D(foo)

10 raster  r  dem  visualisation  color 

ฉันมีข้อมูลบางจุดซึ่งแสดงถึงตำแหน่ง lat-lon รายวันของสัตว์พร้อมการประทับเวลาที่เกี่ยวข้อง ฉันต้องการระบุทุกจุดที่ STATIONARY = TRUE จุดหนึ่งมีคุณสมบัติเป็นนิ่งถ้าบัฟเฟอร์ 100km รอบมันทับซ้อนกันเพิ่มเติม (พูด) 5 จุดที่อยู่ติดกันชั่วคราว ดังนั้นถ้าวันที่ 10เป็นจุดสนใจของฉันฉันต้องการถามว่า 5 วันที่อยู่ติดกันชั่วคราวนั้นอยู่ในระยะ 100km ของจุดนี้หรือไม่ หากวัน 5,6,7,8 & 9; หรือวัน 11,12,13,14 และ 15; หรือวัน 8,9,11,12,13 (ฯลฯ ) อยู่ในบัฟเฟอร์จากนั้น StationARY = TRUE อย่างไรก็ตามหากวันที่ 5,7,9,11 และ 13 อยู่ในบัฟเฟอร์ แต่ไม่ใช่วันที่เป็นทางเลือก (คู่) ในระหว่างนั้นให้ระบุว่า STATIONARY = FALSE ฉันคิดว่าการย้ายบัฟเฟอร์หน้าต่างบางชนิดจะช่วยแก้ปัญหานี้ แต่ฉันไม่รู้วิธีใช้สิ่งนี้ ฉันพยายามที่จะแก้ปัญหานี้ทั้งใน ArcGIS …

10 arcgis-desktop  r  point  buffer  proximity 

ฉันต้องการเรียกคืนชุดอาคารที่ฉันสามารถดูได้จากตำแหน่งที่ระบุ (จุด -lat, lng) และการวางแนว (ช่วงเวลาของมุม) + แน่นอนด้วยความลึกที่กำหนด (เช่น 100 เมตร) โดยหลักการแล้วฉันต้องการส่งออกบางอย่างเช่น: จากมุม 1 ถึงมุม 2 ผู้ใช้สามารถเห็นสิ่งปลูกสร้าง 1 จากมุม 2 ถึงมุม 3 เราจะเห็นอาคาร 2 Data =รูปหลายเหลี่ยมแสดงถึงสิ่งปลูกสร้างของเมืองใน shapefile อินพุต:ตำแหน่งและการวางแนว (ช่วงมุม) เพื่อแสดงการรับรู้ภาพของบุคคล เอาท์พุท : รุ่น 1 - รายการ {BuildingID} รุ่น 2 (ดีกว่า) - รายการ {BuildingID, ชุดของส่วนที่มองเห็นได้สำหรับแต่ละอาคาร (ซุ้มอาคาร)} รุ่น 3 (ผลลัพธ์ในอุดมคติ) - รายการ …

10 qgis  postgis  r  viewshed 

ฉันมีชุดข้อมูลสองจุดใน ArcGIS ซึ่งทั้งสองชุดจะอยู่ในพิกัด WGS84 lat / lon และจุดกระจายทั่วทั้งโลก ฉันต้องการค้นหาจุดที่ใกล้ที่สุดในชุดข้อมูล A ไปยังจุดแต่ละจุดในชุดข้อมูล B และหาระยะห่างระหว่างจุดเป็นกิโลเมตร ดูเหมือนว่าจะเป็นการใช้งานเครื่องมือใกล้เคียงอย่างสมบูรณ์แบบ แต่นั่นทำให้ฉันได้ผลลัพธ์ในระบบพิกัดของจุดอินพุต: นั่นคือองศาทศนิยม ฉันรู้ว่าฉันสามารถฉายข้อมูลอีกครั้ง แต่ฉันรวบรวม ( จากคำถามนี้ ) ว่ามันยาก (ถ้าไม่เป็นไปไม่ได้) เพื่อหาเส้นโครงที่จะให้ระยะทางที่แม่นยำทั่วโลก คำตอบสำหรับคำถามนั้นแนะนำให้ใช้สูตร Haversine เพื่อคำนวณระยะทางโดยใช้พิกัดละติจูด - ลองจิจูดโดยตรง มีวิธีการทำเช่นนี้และได้ผลลัพธ์เป็นกิโลเมตรโดยใช้ ArcGIS หรือไม่? ถ้าไม่เป็นวิธีที่ดีที่สุดในการเข้าถึงนี้คืออะไร?

10 arcgis-desktop  coordinate-system  r  distance  units 

ฉันพยายามทำให้ร่างแสดงข้อมูล azimuthal ด้วยช่วงความไม่แน่นอนที่แตกต่างกันในแต่ละจุด รูปเก่าแก่ของโรงเรียนจากกระดาษ 1991 จับแนวคิด "bowtie plot" ที่ฉันตั้งเป้าไว้: ข้อเสนอแนะใด ๆ เกี่ยวกับวิธีที่ฉันจะทำให้ร่างที่คล้ายกัน? ฉันเป็นมือใหม่ที่เกี่ยวข้องกับ GIS แต่ฉันสามารถเข้าถึง ArcGIS ผ่านมหาวิทยาลัยของฉันได้ ประสบการณ์ Arc ของฉันถูก จำกัด ในการสร้างแผนที่ทางธรณีวิทยาดังนั้นฉันจึงไม่ต้องทำอะไรแปลกใหม่เกินไป ฉันแหย่ไปในตัวเลือกสัญลักษณ์ใน Arc และ QGIS แต่ไม่เห็นการตั้งค่าใด ๆ ที่ฉันคิดว่าจะทำงาน โปรดทราบว่านี่ไม่ได้เป็นเพียงเรื่องของการหมุนสัญลักษณ์รูปโบว์โดย azimuth; ช่วงเชิงมุมของ "bowtie" แต่ละอันจะต้องแตกต่างกัน ฉันให้คะแนนทักษะ Python ของฉันเป็น 'strong ตัวกลาง' และทักษะ R ของฉันเป็น 'low middle' ดังนั้นฉันไม่รังเกียจที่จะแฮ็คข้อมูลบางอย่างร่วมกับmatplotlibและmpl_toolkits.basemapหรือไลบรารีที่คล้ายคลึงกันหากจำเป็น แต่ฉันคิดว่าฉันขอคำแนะนำที่นี่ก่อนที่จะลงไปที่ถนนในกรณีที่มีวิธีแก้ปัญหาที่ง่ายกว่าจาก GIS-land ที่ฉันไม่เคยได้ยินมาก่อน

10 python  algorithm  r  visualisation  spherical-geometry 

ฉันกำลังพยายามดึงค่า bathymetry ของพื้นที่ที่ฉันสนใจจากเลเยอร์แรสเตอร์บทกวีโลกโดยใช้ฟังก์ชัน 'rasterize' ในแพ็คเกจ {sp} * การแก้ไข: ฉันพบฟังก์ชัน 'แยกข้อมูล' ซึ่งดูเหมือนจะเป็นสิ่งที่ฉันกำลังมองหา นี่คือสิ่งที่ฉันทำไปแล้ว: > class(subarea0) #This is my area of interest (Eastern Canadian Arctic Sea) [1] "SpatialPolygons" attr(,"package") [1] "sp" > extent(subarea0) class : Extent xmin : -82.21997 xmax : -57.21667 ymin : 60.2 ymax : 78.16666 library(marelac) data("Bathymetry")#World bathymetric data …

10 raster  r  overlay 

ฉันต้องการทำงานกับข้อมูลแรสเตอร์ใน R แต่ฉันมีปัญหาในการส่งออกชุดข้อมูลแรสเตอร์ของฉันจาก ArcMap โดย "export data" / GRID ฉันได้รับไฟล์ที่มีนามสกุล * .aux สิ่งที่ฉันต้องการคือ * .grd ใครช่วยฉันออกได้ไหม

10 arcmap  raster  r 

ฉันพยายามอ่านไฟล์แรสเตอร์ในรูปแบบ. DEM บน windows โดยใช้แพ็คเกจ 'raster' ใน R ฉันพบปัญหากับค่า NA เมื่อโหลดข้อมูลลงใน R ใน Windows 7 แต่ฉันไม่มีปัญหาใน Mac ที่มี OSX Lion บน windows ค่า NA ดูเหมือนจะอ่านไม่ถูกต้อง คำถามคือทำไมสิ่งนี้เกิดขึ้น? ไฟล์แรสเตอร์ที่ใช้ถูกดาวน์โหลดจาก USGS ด้วยรหัส R ต่อไปนี้: download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz') untar('e020n90.tar.gz') จากนั้นฉันก็อ่านแรสเตอร์เป็น R โดยใช้แพ็คเกจ 'raster' ใน OSX Lion และ R64 เวอร์ชัน 2.13.1 ค่า NA จะได้รับการยอมรับ: > onMac …

10 raster  r  windows  osx 

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

10 spatial-analyst  spatial-statistics  r  clustering 

ฉันทำงานกับนักพัฒนาเว็บที่ทำงานกับ Google Maps API พยายามกำหนดขั้นตอนการประมวลผลล่วงหน้าที่ฉันควรทำ (เช่นใน R และ GRASS) และสิ่งที่สามารถทำได้ภายใน google maps API นี่ทำให้ฉันมีคำถามต่อไปนี้: มีวิธีการแปลงแรสเตอร์เป็นเวกเตอร์ข้อมูลใน Google Maps API หรือไม่? มีเครื่องมือทางสถิติใด ๆ ที่เกี่ยวข้องกับ API หรือไม่

10 google  r 

ฉันมีไฟล์รูปร่าง ESRI ที่มีขนาดใหญ่มาก (> 8000 polygons, .shp ไฟล์> 32MB) ซึ่งฉันต้องการอ่านเป็น R บน Windows 7 ฉันคุ้นเคยกับ rgdal และ readOGR () อย่างไรก็ตามด้วยขนาดของ shapefile การอ่าน shapefile ทั้งหมดในหน่วยความจำจึงเป็นกระบวนการที่ช้ามาก นอกจากนี้ฉันสนใจเพียงไม่กี่โพลิกอนคลัสเตอร์ (<100) จาก 8,000 ที่มีอยู่ในรูปร่างไฟล์ กระบวนการที่ฉันต้องทำคือทำซ้ำอย่างมากนั่นคือต้องโหลดรูปร่างไฟล์หลายครั้งดังนั้นจึงช้ามาก ดังนั้นฉันสงสัยว่ามันเป็นไปได้ที่จะอ่านเฉพาะในส่วนของ shapefile เพื่อเร่งกระบวนการ ดูเหมือนว่าไม่น่าจะโหลดข้อมูลจำนวนมหาศาลเช่นนั้นทิ้งส่วนใหญ่

10 shapefile  r  rgdal 

เมื่อฉันต้องการที่จะดำเนินการ kriging บางครั้งก็ใช้งานได้ขึ้นอยู่กับค่าที่ฉันใช้ใน DataTable ของฉัน เป็นผลมาจากฟังก์ชั่น krige ฉันได้รับvar1.pred: NA NA NA ...และvar1.var: NA NA NA ...(แต่เฉพาะเมื่อฉันใช้ค่า "ผิด" ใน DataTable ของฉัน) ตัวอย่างเช่น: มันทำงานได้เสมอ (จนถึงตอนนี้) เมื่อฉันใช้เพียง 10 ค่า มันใช้งานได้เมื่อฉันใช้ 50 ค่า แต่ใช้ได้กับบางค่าเท่านั้น มันใช้งานไม่ได้เมื่อฉันใช้ 50 ค่าและค่า "ผิด" ใช้งานได้เมื่อฉันใช้ 25 ค่าและค่าที่กล่าวถึงก่อนหน้านี้ "ผิด" ฉันไม่เข้าใจว่าทำไมบางครั้งมันถึงทำงานได้และบางครั้งก็ไม่ สิ่งที่แปลกคือเมื่อฉันเพิ่มลงZwiesel;49.02999878;13.22999954;2.2ใน DataTable มันจะทำงานเมื่อฉันใช้น้อยกว่า ~ 20 ค่า แต่มันไม่ทำงานเมื่อฉันใช้มากกว่า 50 ค่า ... ความผิดพลาดของฉันอยู่ที่ไหน …

10 r  interpolation  kriging  gstat