Plot 全球柵格文件中的 GAM 值

[英]Plot GAM values in a worlwide raster file

我正在嘗試制作一個 SDM（物種分布模型），為此我構建了一個 GAM model（假設物種 A： 在這種情況下，“存在”指的是二元因子結構，其中 1 表示存在，0 表示不存在（背景點）。 此二元因子結構與一組數字柵格預測變量（例如，高程、坡度和各種生物氣候變量）相結合。 注意：存在/不存在數據是全局 ...

R terra:: 累積和重置 function 跨堆棧層

[英]R terra:: cumulative sum reset function across layers of a stack

我有一個包含大量層的堆棧，我想在總和變為 <0 時通過重置跨層運行像素累積和 function。 我編譯了以下代碼： ，效果很好，但我有一大堆光柵（就 nlyrs 而言）。 我想知道是否有更快的方法來執行此任務。 ...

如何用.netCDF土地利用數據和Key Biodiversity Areas (.shp)數據進行空間數據分析

[英]How to conduct spatial data analysis with netCDF land use data and Key Biodiversity Areas (.shp) data

我有一個 .netCDF 文件（“SSP119.nc”），其中包含 2015 年到 2100 年的土地利用變化信息，包括在內（來源： https://luh.umd.edu/data.shtml ）。 我還有一個包含全球關鍵生物多樣性區域 (KBA) 多邊形的 shapefile。 我想知道在某些年 ...

從 shp 文件中獲取瑞士的救濟

[英]Get relief of Switzerland from a shp-file

我想要通過瑞士的 map 繪制浮雕的官方來源。 我下載了 https://cms.geo.admin.ch/ogd/topography/DHM25_BM_SHP.zip來使用 dhm25_p.shp-File 現在使用代碼 我會得到錯誤 ...

從 tiff 文件中解除瑞士

[英]relief of Switzerland from a tiff-file

我想要通過瑞士的 map 繪制浮雕的官方來源。 在 swisstopo 的其中一頁上，我找到了以下文件： https://data.geo.admin.ch/ch.swisstopo.pixelkarte-farbe-pk1000.noscale/swiss-map-raster1000_2022_ ...

循環中一個多邊形的柵格裁剪中的新柵格 R

[英]New rasters from crop of rasters by one polygon in loop R

從 19 個氣候柵格中，我嘗試僅裁剪多邊形中的柵格區域（在法國）。 為此，我嘗試生成一個循環，用多邊形逐一裁剪每個柵格，如下所示： 它有效，但現在我試圖找到一個循環來一次裁剪和遮罩 19 個柵格，並從我會命名的這個裁剪和遮罩中創建新柵格。 有什么建議么？ 我不知道是否清楚，如果您需要更多信息，請告 ...

使用 R 讀取 grib2 數據集而不丟失元數據

[英]Reading a grib2 dataset without loss of metata with R

我的問題遵循之前的問題，該問題已得到部分回答： Reading a grib2 dataset with 4 dimensions and 2 variables with R 我正在嘗試讀取帶有 R 的 GRIB2 文件。此文件是概率氣象預報：10 個變量、1 個提前期、17 個經度、23 個緯度 ...

將柵格范圍值從 0/100 更改為 1/100

[英]Change raster range values from 0/100 to 1/100

我有一個范圍從 0（最小值）到 100（最大值）的柵格圖層，我想將其值更改為 1 到 100 的范圍。重要的是數據不要改變，基本上我想要相同的柵格但不同范圍的值。 我可以在 R 和 QGIS 上工作。 謝謝大家的幫助：:) ...

哪個 package 應該用於大型數據集、terra 或柵格的算術函數？

[英]Which package should be used for arithmetic functions on large datasets, terra or raster?

我需要對很多大柵格（28000 個像元，181 層）進行計算。 我在一小部分（24 個單元，181 層）上嘗試了我的代碼。 我從這個論壇獲得了幫助，以盡可能多地進行優化。 現在我在光柵 package 中使用了磚塊，因為我讀到磚塊被加載到 memory 中並且處理起來更快。 但后來有人認為 ter ...

從 lapply 的結果中提取正確的列 R

[英]Extract the Correct Column from the Result of lapply in R

我有一個測深柵格和一個邊界網格化的 shapefile。 我無法隔離 lapply 生成的列表的“值”列。 目的是使用柵格平均值列作為 shapefile 中的新列 最后一行應該分配 Mean_Raster_Values 列表 object 中的整個第三列，但 [[3]] 只提供第三行如何從 Me ...

在 R 中的光柵堆棧/磚塊上應用復雜的數學函數，並創建兩個相互依賴的不同光柵堆棧

[英]Applying complex mathematical functions on Raster Stacks/Bricks in R and creating two different Raster stacks which are dependent on each other

我有一個光柵磚（ncell=28536 和 nlayers=181）。 我需要在原始積木上運行數學函數並創建另外兩塊相同大小的積木。 兩個 output 塊相互依賴的地方。 inputBrick 有 181 層，每層有 28536 個單元格。 outputBrick1 將通過分析 outputB ...

在 ggplot 中調整重疊柵格的 alpha

[英]Adjusting alpha on overlaid rasters in ggplot

我正在嘗試使用 ggplot 突出顯示 DEM 的一部分——本質上，只是一個彩色的、半透明的區域子集。 只要我不向疊加子集添加“alpha”參數，一切看起來都很好——我可以將我的 DEM plot 作為 geom_raster，並在頂部添加我的第二層，作為另一個 geom_raster 或 geo ...

錯誤：列中缺少數據：使用護林員運行隨機森林回歸時彈出 package

[英]Error: Missing data in columns: pop when running random forest regression using the ranger package

我正在嘗試使用 R 中的ranger R實現隨機森林 (RF) 回歸，但我收到此錯誤： Error: Missing data in columns: pop （ pop是我的自變量）在運行predict function 時。 作為參考，當使用randomForest package 時，我可以 ...

當經度為 c(-180, 180) 但 griddap 數據集經度為 0-360° 時，如何使用 library(griddap) 中的 rerddap function 獲取數據？

[英]How to get data using the rerddap function in library(griddap) when longitude is c(-180, 180) but the griddap dataset longitude is 0-360°?

R 中的庫 (rerddap) 中的 griddap function 提供了一個優秀的 function，用於以簡單直接的方式從 ERDDAP 服務器訪問衛星數據。 下面是一個可重現的示例，它使用ncdcOisst21Agg_LonPM180數據集提取大西洋的最新海面溫度 (NOAA OISS ...

如何更改柵格網格單元值的比例

[英]How to change the proportion of raster gridcell values

在示例柵格r ，如何將像素值84的比例從9%更改為25% ，如您在per_84中看到的那樣？ 即我想在以下條件下增加像素值為84的像素數並減少其他值的像素數。 84的比例從9%增加到25%如下：柵格中已經有9% （共84個）。 在剩余的16%中， 6%來自像素值90 ， 5%來自像素值85 ，剩余 ...

Sentinel-1 & 2：在不同圖塊的圖像中創建重合點

[英]Sentinel-1 & 2: create coincident points in images of different tiles

我有一個 Sentinel-1 圖像和 Sentinel-2 圖像，我想讓 Sentinel-1 中的像素與 Sentinel-2 重合。 更簡單的解決方案是使用S2 <- resample(S2, S1)將兩個圖像放在相同的幾何形狀（空間范圍和分辨率）中，但我擔心光譜信息的完整性。 更好地重 ...

打開帶有光柵的 a.nc 文件時出錯 package

[英]Error when opening a .nc file with raster package

我是 R 中光柵 package 的新手，我試圖用 package 光柵打開一個 .nc 文件，但彈出了一些錯誤。 如果您想嘗試，我使用的是來自哥白尼的 2018 年和 2019 年每月海水鹽度數據集（網格是魁北克聖勞倫斯河口和加斯佩西海岸）。 以前打開過類似的數據文件，但從來沒有出現過這個錯誤， ...

堆疊柵格並計算每個像素的最大值，然后在 R 中的 rest 個圖層中保留該值

[英]Stack rasters and calculate each pixel's max value then keep that value for rest of the layers in R

我每天堆疊 181 個 ndvi 光柵。 我想為每個像素找到峰值，然后將峰值分配給 rest。原始數據曲線： 這就是我要的： 此結果適用於單個像素。 我已經在 excel 中完成了，但為了方便起見，我想在 R 中創建一個代碼。 現在每個像素在不同的日期達到峰值（因此堆棧中的不同層）。 這是我試過 ...

將刻度顏色更改為 R 中的特定數字

[英]Changing scale colour to specific number in R

我已將我的數據上傳到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 設置我的數據： 以下是結果圖： 正如您所看到的，盡管它們具有不同的值，但比例都是相同的。 我怎樣才能改變它，使一張圖表上的 4.0 與 ...

如何交換 x 軸和 y 軸上的緯度和經度？

[英]How to swap latitude and longitude on x and y axes?

我已將我的數據上傳到 github： https://github.com/dschmidt42/Benguela_Chla_Danielle_Schmidt 。 它涉及非洲東海岸的 chl a 數據，時間跨度為 26 年。 我試圖將數據分成每個月，然后比較幾個月內的差異。 這是我第一次使用 .n ...