r-raster

    1熱度

    2回答

    我需要從grib2文件(我將其放入柵格堆棧)中有效地提取溫度數據。堆棧中的每個柵格圖層都表示一個時間點。 現在我需要爲每個觀察值(x,y,t)提取一個單值。下面的代碼可以完成這項工作,但這需要花費太多時間。任何提高效率的建議非常感謝。 files <- list.files(path="Weather/NCEP/temperature_3hour_forecast", full.names = T

    0熱度

    1回答

    我有一個具有非負值(每個像素都有一個高度值)的柵格矩陣,並且想要使用5x5移動窗口計算一些自定義高度指標,使用對稱填充矩陣的邊緣。我使用QGIS中的pktools(空間過濾器)做了一些基本的計算(平均值,標準差,最大值),這使我可以做到這一點(使用零填充,複製或對稱填充)。 現在我想使用自定義函數(即高度係數,偏度和峯度)來計算一些其他指標,我不能使用pktools。我一直試圖使用raster中的

    1熱度

    2回答

    我有從日常氣候數據創建的柵格堆棧。可以在這裏找到: #!/bin/bash wget -nc -c -nd http://northwestknowledge.net/metdata/data/tmmx_1982.nc 目標是從這些日常記錄中獲得每月95%的溫度值。每當我使用raster包中的calc時,它只返回一層而不是12(例如,12個月)我錯過了什麼? library(raster)

    0熱度

    2回答

    我正在使用R(非RStudio)3.3.3在Mac0S10.10.5中的R柵格軟件包中處理地理參考圖像。使用GUI將工作區保存到xxx.RData文件後,關閉會話重新啓動R(無需在啓動時讀取歷史記錄文件),然後使用GUI再次加載工作區。當我輸入ls()時,所有對象(包括我加載的原始完整全局柵格以及使用crop()創建的子集)都會出現。但是,當我繪製柵格時,只會顯示原始非裁剪圖像。裁剪的子集不繪製,

    -1熱度

    1回答

    我有一個帶有2個波段的柵格,我想計算3x3移動窗口中的協方差,並用協方差線代替焦點值。我正在研究柵格包的焦點函數,並且還提供了一些使用getValuesBlock for循環的提示。如果有人指出我解決這個問題,我會非常感激。

    0熱度

    2回答

    我有一個光柵,有50多個波段。我想要搜索柵格第一個波段內的所有像素== 0。以下我想在其他光柵波段中設置所有這些像素NA。因此,我不必再次搜索50次以上。這是我的例子: cl_input <- brick("sometif") for(i in 1:nlayers(cl_input)){ print(names(cl_input[[i]])) cl_input[[i]][cl

    6熱度

    2回答

    基本上,我用ASCII的形式計算了一個全局分佈概率模型,例如: gdpm。 gdpm的值全部爲0和1之間 然後我輸入從形狀文件的本地地圖: shape <- file.choose() map <- readOGR(shape, basename(file_path_sans_ext(shape))) 下一步,我光柵化gdpm,並裁剪使用本地地圖: ldpm <- mask(gdpm, ma

    0熱度

    1回答

    我已經收到一個.tif文件,其中包含有關每個值的顏色信息(例如,森林是綠色的,雲是黑色的).qml文件。 library(raster) library(rasterVis) sm <- raster("UT_classification_coded.tif") ########## Classify raster/build attribute table sm <- rat

    0熱度

    1回答

    好吧,我有女傭在科羅拉多州的地圖。我喜歡能夠遮蔽精選縣。我的數據給了我縣的名字,但我不知道如何將它應用到我的ggmap。 這是從我的多邊形 [email protected][["NAME_2"]] [1] "Adams" "Alamosa" "Arapahoe" "Archuleta" "Baca" [6] "Bent" "Boulder" "Broomfield" "Chaffe

    2熱度

    1回答

    library(raster) r <- raster('glc2000_v1_1') # http://forobs.jrc.ec.europa.eu/products/glc2000/products/glc2000_v1_1_Grid.zip extent(r) # class : Extent # xmin : -180.0045 # xmax : 179.9955 #