大型 RasterStack 上的錯誤計算分位數
[英]Error computing quantiles on large RasterStack
我想在 R 上的一組堆疊的 RasterLayers 上跨層獲取每個單元格的分位數。 我的 RasterStack 看起來像這樣尺寸：721、1440、1038240、18628（nrow、ncol、ncell、nlayers） 我運行以下行以獲得每個單元格的值為第 25 個和第 75 個百分位數的 ...
如何讀取 GeoTIFF 文件並以 R 語言返回特定數據類型？
[英]How to read a GeoTIFF file and return in a specific datatype in R language?
我正在嘗試通過raster package 讀取 R 中無符號字節數據格式的大型 GeoTIFF 文件（18000×18000×40）。我希望將其讀取為 R 中的raw數據類型（無符號字節），然后由Rcpp處理。 但是， raster package 的函數似乎只返回numeric數據類型（8 字 ...
Function 用於對焦方法
[英]Function for the focal method
我正在嘗試使用 focal 運行一個簡單的 function，它在每個 window 中取 5 個隨機值，然后根據這些值計算方差系數 (sd/mean)。 下面的代碼當前返回錯誤：評估錯誤：第一個參數無效。 提前致謝！ ...
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 年每月海水鹽度數據集（網格是魁北克聖勞倫斯河口和加斯佩西海岸）。 以前打開過類似的數據文件，但從來沒有出現過這個錯誤， ...
使用邊界定義的焦點填充 NA 柵格單元
[英]fill NA raster cells using focal defined by boundary
賞金將在 2 天后到期。 此問題的答案有資格獲得+50聲望賞金。 89_Simple想讓更多人關注這個問題。 我有一個柵格和一個 shapefile。 柵格包含 NA，我正在使用focal function 填充 NA 但是，還有一些 NA 仍然存在。 所以我這樣做： 填滿所有 NA 后，我再次 ...
與多個人創建空間線 Object
[英]Creating a Spatial Line Object with multiple individuals
我正在嘗試為我的數據集中的多個人創建一條空間線，然后跨線提取柵格值。 我已經嘗試了一些事情，但不幸的是我正在使用的代碼將所有線路分配給一個人。 WD是其中有37個人的數據集。 x 和 y 是坐標。 如果有人知道如何執行此操作，我將不勝感激。 謝謝！ 我希望能夠沿線提取柵格值。 ...
為 R 中的每個 GPS 位置從 RasterBrick 的多層中提取值
[英]Extract values from multiple layers of a RasterBrick for each GPS location in R
我正在嘗試計算 GPS 位置的 USDA CropScape 作物價值，該位置的日歷年、上一個日歷年和下一個日歷年（例如，如果 GPS 位置的日期是 10/1/2016 ，我想提取該位置 2015、2016 和 2017 年的作物價值）。 這是我用來准備 GPS 和 CropScape 數據的代碼： ...
如何使用單層 shapefile 屏蔽 multi-layer.netCDF 或柵格？
[英]How do I mask a multi-layer netCDF or raster with a single-layer shapefile?
我目前正在處理 .netCDF 格式的每日降水數據。 數據分辨率為 4 公里，覆蓋美國。 但是，我想為特定地理區域（大約一個縣的大小）使用更高分辨率的 shapefile 來屏蔽/剪輯數據。 最終，我希望 output 成為更小區域的每日降水數據，無論是高分辨率還是原始 4km 分辨率。 我嘗試了 ...
levelplot 未在圖例上顯示分類類
[英]levelplot not showing categorical classes on legend
嘗試在 R 中使用 levelplot map 分類柵格，但圖例始終顯示數字 ID，而不是分類類。 請參見下面的示例。 任何人遇到這個問題並且知道可能是什么問題？ 光柵 map 使用 levelplot ...
嘗試使用 raster::projectRaster 重新分配 crs 時出錯：rbind 中的錯誤：矩陣的列數必須匹配
[英]Error when attempting to reassign crs with raster::projectRaster: error in rbind:number of columns of matrices must match
我有以下名為 pred_th 的 SpatRaster 對象。 我正在嘗試將 crs 重新分配給 ESPG 2970。當我嘗試使用raster函數projectRaster執行以下操作時： 我收到此錯誤消息： 我很困惑，如果有人可以提供幫助，我將不勝感激。 ...
使用 R 將柵格 DEM 分辨率從 30m 更改為 500m 時，最佳使用方法是什么？
[英]What is the best way to use when changing raster DEM resolution from 30m to 500m using R?
我正在嘗試使用 R 以 500 米的分辨率使用 SRTM-DEM 計算高程、坡度和坡向。SRTM-DEM 坐標系在 UTM 中，分辨率為 30m。 有人可以告訴我如何將分辨率更改為 500m 並以相同的范圍計算高程、坡度和坡向。 我想知道我是否可以使用聚合方法獲得准確的 500m 分辨率。 ...
R中如何將正弦遙感數據導入柵格？
[英]How to import sinusoidal remotesensed data into raster in R?
我正在嘗試將來自氣候變化倡議 ( CCI ) 的 NetCDF 文件讀入 R 和terra package。鑒於數據不在規則網格上，我試圖找到將此類數據投影到規則網格上的正確方法.library(terra) #> terra 1.6.47 # Read the data r <- ...
R ggplot 繪圖 map 圓形柵格 - 如何刪除投影區域外的數據？
[英]R ggplot plotting map raster with rounded shape - How to remove data outside projected area?
我正在嘗試 plot 在協調系統中投影的柵格，該系統像大多數非 WGS84 投影一樣遵循地球曲率。 問題是地球環繞數據的地方不應該繪制在地球之外。 我意識到 ggplot 不能做一個圓形/橢圓形 plot 但我如何屏蔽或自動刪除地球以外的數據？ 我必須 plot 超過 100 張地圖，我不能手動執 ...
R terra：跨層滾動求和
[英]R terra: rolling sum across layers
擁有大量時間序列柵格，我想創建另一個 rasterstack，它代表圖層中每個像素的 5 天滾動總和。 鑒於光柵堆棧的大小， raster::calc會為此任務無休止地運行，因此我想知道是否可以使用terra::lapp更快地實現此操作。 然而： 返回錯誤： 如果有人可以建議錯誤出在哪里，或者建議另 ...
是否有 R function 將 rasterstack 應用於數據幀列表
[英]Is there an R function to apply rasterstack to a list of dataframes
我是全新的 TR(studio) 並嘗試構建光柵堆棧。 我的輸入大約有 600jp2 個文件，它們都在同一個目錄中，看起來像這個目錄中文件的圖像 ( https://i.stack.imgur.com/Mcvpn.png ) 文件已排序，因此波段編號 2、3、4 和 8 始終排在 4 行中。據此，我 ...
將 long/lat 點添加到 map 的 landsat 數據
[英]Add long/lat points to map of landsat data
我有一個用於陸地衛星圖像的 RasterBrick： 我想將經度/緯度點添加到 landsat_stack_brick 的landsat_stack_brick 這些點應該都沿着島嶼的海岸線，但這是我嘗試將它們添加到 map 時的樣子： 到目前為止使用的代碼是：[編輯：應用 John Polo 提 ...
使用 R 中的 sp、terra、raster 包將阿爾伯斯等面積轉換為緯度/經度
[英]Using sp, terra, raster packages in R to convert Albers Equal Area to Lat/Long
我有一個關於地球上某個區域的地面溫度的.tif 文件，我想找出該區域的坐標。 這是我正在處理的文件的鏈接： 2017001D.tif 在 R 中使用raster R ，我能夠將 .tif 文件加載為RasterLayer class。 我可以將坐標信息提取為 我還使用terra package 將其 ...
如何計算 terra package 中的聚合事實參數？
[英]How to calculate aggregate fact parameter in terra package?
terra package 具有聚合 function，允許創建具有較低分辨率（較大單元格）但需要事實參數的新 SpatRaster。 轉換很多柵格時，每次都需要計算fact，有沒有辦法根據另一個柵格的目標分辨率傳遞fact參數？ 其他函數將現有柵格作為輸入，例如 function(r1,r2) ...
空間點到多邊形查詢
[英]Spatial Points to Polygons query
我有一個關於將空間數據轉換為 R 並將其從 R 導入 QGIS 的問題。 我有南極海冰濃度的 GeoTiff，可從以下鏈接下載： https://seaice.uni-bremen.de/databrowser/#day=13&month=10&year=2022&img= ...