R:使用连接和栅格包从gz文件中读取GeoTiff
R: 使用连接和栅格包从gz文件中读取GeoTiff
答案：
在R语言中，可以使用连接和栅格包（raster package）来从gz文件中读取GeoTiff格式的数据。连接包（rgdal package）提供了与GDAL（Geospatial Data Abstraction Library）的接口，可以处理各种地理空间数据格式，而栅格包（raster package）则提供了处理栅格数据的功能。
以下是一种完整的步骤，用于从gz文件中读取GeoTiff数据：
- 安装和加载所需的包：
install.packages("rgdal")
install.packages("raster")
library(rgdal)
library(raster)- 创建一个连接对象，用于打开gz文件：
gz_file <- "path/to/your/gz/file.gz"
gz_conn <- gzfile(gz_file, "rb")- 使用连接对象读取GeoTiff数据：
geotiff <- readGDAL(gz_conn)- 关闭连接对象：
close(gz_conn)现在，你可以使用raster包中的函数来处理和分析读取的GeoTiff数据了。例如，你可以使用plot函数来可视化数据：
plot(geotiff)这是一个简单的例子，展示了如何使用连接和栅格包从gz文件中读取GeoTiff数据。根据具体的应用场景和需求，你可以进一步探索和使用raster包中的其他函数和方法。
推荐的腾讯云相关产品：腾讯云地理空间计算服务（Tencent Cloud Spatial Computing Service），该服务提供了丰富的地理空间数据处理和分析能力，可以帮助用户高效处理和管理地理空间数据。
更多关于腾讯云地理空间计算服务的信息，请访问：腾讯云地理空间计算服务
相关·内容
2回答
R:使用连接和栅格包从gz文件中读取GeoTiff
r、r-raster、gzip、geotiff
我想在不下载的情况下从服务器上读取一个压缩后的GeoTIFF。我只是不想创建太多的临时文件，我必须在以后删除这些文件。
包含.csv.gz文件的。
使用download时，我使用以下方式：
library(raster)
link <- "ftp://ftp.glcf.umd.edu/glcf/SRTM/Degree_Tiles/n000/SRTM_ff03_n000e010
/SRTM_ff03_n000e010.tif.gz"
download.file(link, "test.tif.gz")
gunzip("test.tif.gz"
浏览 13提问于2017-02-07得票数 3
1回答
如何叠加和显示光栅？
r、raster、rgdal
我在一个文件夹中有12个光栅文件，我想在其中使用RasterVis的平面图。
kpacks <- c('tiff','rgdal','raster','sp','rasterVis')
new.packs <- kpacks[!(kpacks %in% installed.packages()[,"Package"])]
if(length(new.packs)) install.packages(new.packs)
lapply(kpacks, require, character.
浏览 4提问于2016-10-11得票数 0
回答已采纳
1回答
R:直接从web url读取geotiff数据(httr：：获取原始内容)
r、raster、httr、rgdal、geotiff
我想从服务器提供的RasterLayer数据创建一个GeoTIFF。我将使用httr::GET调用来查询服务器上的数据(数据是按需提供的，因此在应用程序中，不会有以.tif结尾的url，而是一个查询url)。
将此调用的结果作为GeoTIFF文件写入磁盘后，很容易从磁盘上生成的GeoTIFF文件创建RasterLayer：
library(httr)
library(raster)
url <- 'http://download.osgeo.org/geotiff/samples/gdal_eg/cea.tif'
geotiff_file <- tempfile
浏览 4提问于2018-09-05得票数 6
回答已采纳
1回答
通过R Shiny的downloadHandler下载多个GeoTIFF格式的栅格文件时出错
r、shiny、download、zip、raster
我正在尝试制作一个闪亮的基于web的应用程序，它可以接受多个二进制文件作为用户的输入，并为每个输入文件创建一个GeoTiff栅格文件供用户下载。但是，我无法使用R shiny的DownloadHandler将多个输出的GeoTiff文件作为压缩文件夹下载。我尝试过的代码如下所示。
#Import required packages
library(shiny)
library(raster)
library(rgdal)
library(zip)
library(stringr)
#User Interface
ui = fluidPage(
titlePanel(title = "
浏览 1提问于2021-06-24得票数 0
2回答
如何将ncdf (of Sentinel_5_L2)文件转换为geotiff？
r、raster、geotiff、ncdf4、sentinel-5p
我正在尝试使用netcdf文件创建geotiff文件。我有这个文件(Sentinel_5_L2)。但是我的代码没有运行。有人知道创建geotiff文件吗？ 文件https://wetransfer.com/downloads/7c5692bef082bc64999e10d8a7a91e8f20211116115821/7774f8 这是错误 Error in .local(x, ...) :
unused arguments (xmin = -179.996887207031, ymin = c(....) 这是代码。 library(ncdf4)
library(raster)
浏览 60提问于2021-11-17得票数 0
5回答
栅格软件包绘图函数的命名空间问题？
r、namespaces、raster
我遇到了一个关于'raster‘包的'plot’函数的奇怪问题。
library(raster)
ras <- raster(ncol=10, nrow=10)
编辑
values(ras) <- runif(ncell(ras))
结束编辑
plot(ras)
Erreur dans as.double(y) :
cannot coerce type 'S4' to vector of type 'double'
对于我在网上读到的内容，这个错误取决于用户，也可能取决于加载的包。在我的例子中，问题来自于r使用来自'gra
浏览 3提问于2013-01-21得票数 2
回答已采纳
1回答
将shapefile转换为栅格
r、shapefile、rasterize
我在栅格化shapefile以在0.5*0.5网格上生成点时遇到了问题。shapefile表示全球珊瑚礁的风险级别(低-0，中-100，高-1000，极高-1500)到综合威胁的分类。
我从另一个运行良好的示例中提取了代码，但是当我对我的数据尝试它时，我从plot函数中得不到任何东西。有关shapefile和我的代码的链接，请参见以下内容：
# Read shapefile into R
library(rgdal)
library(raster)
int.threat.2030 <- readOGR(dsn = "Global_Threats/Integrated_
浏览 7提问于2016-01-30得票数 3
回答已采纳
2回答
高效读取大型多波段栅格
r、raster、gdal
我正在使用rgdal包在R中编写一个图像分类脚本。所讨论的光栅是一个具有28个通道的PCIDSK文件:一个训练数据通道、一个验证数据通道和26个光谱数据通道。目标是填充一个数据帧，该数据帧包含训练数据通道中非零的每个像素的值，以及26个频带中的相关谱值。
在Python/Numpy中，我可以轻松地将整个图像的所有波段导入到多维数组中，然而，由于内存限制，R中唯一的选择似乎是逐块导入此数据块，这非常慢：
library(rgdal)
raster = "image.pix"
training_band = 2
validation_band = 1
BlockWidth = 5
浏览 2提问于2010-11-16得票数 2
回答已采纳
17回答
rgdal包安装
r、maps
这里的问题并不是如何通过R绘制地图，因为我已经发现了一个非常好的示例，而是如何使其工作。实际上，我无法加载库rgdal
library(rgdal)
Error in library(rgdal) : there is no package called ‘rgdal’
但是，当我尝试手动安装上面的包时，我得到了以下错误：
....
configure: error: proj_api.h not found in standard or given locations.
ERROR: configuration failed for package ‘rgdal’
* removing ‘/
浏览 4提问于2013-03-06得票数 113
回答已采纳
1回答
在R中打开栅格图像时，像素值以某种方式重新缩放时出现问题
r、r-raster
我使用光栅函数在R中打开了一张Geotiff图像。然而，图像似乎被重新缩放了。最小值和最大值应为13607和15461，但应为275%和305%。Geotiff图像在GIS中打开时是正确的，但在R中不正确。
这是我的代码：
脚本：
library(raster)
trial<-raster("MOD11A2.A2000049.h17v05.006.2015058135048.tif", datatype = "INT2U")
trial
plot(trial)
浏览 1提问于2020-03-23得票数 0
我有这样的数据集(pts)：
x <- seq(-124.25,length=115,by=0.5)
y <- seq(26.25,length=46,by=0.5)
z = 1:5290
longlat <- expand.grid(x = x, y = y) # Create an X,Y grid
pts=data.frame(longlat,z)
names(pts) <- c( "x","y","data")
我知道我可以通过执行以下操作将数据帧(Pts)映射到映射中：
library(sp)
li
浏览 1提问于2012-03-30得票数 1
回答已采纳
1回答
适用于readGDAL的plot3d？
r、plot、raster
我想在R中准备3D图形！使用栅格数据。我听说了光栅包中的命令plot3d。不幸的是，我有GeoTiff格式的数据，所以不可能直接通过光栅包读取它(而不是它，为了加载这个文件，我使用readGDAL命令)。
有没有机会使用这样的数据来准备3D图形？如何检索坐标？我知道如何使用as.matrix命令获取栅格值，但因为它是GeoTiff格式的，所以我无法使用xmin或xmax。
浏览 0提问于2011-09-12得票数 0
1回答
光栅:提取光栅图像的中心(geoTIFF)并绘制其边界
r、raster、r-raster
我想在geoTIFF中提取4个光栅子图像中心的坐标(并创建一个数据框)，并绘制最初在r2图像中拆分的2个子图像的边界，为此我尝试： library(raster)
library(rgeos)
#create GeoTIFF raster
r <- raster(ncol=100, nrow=100)
s <- stack(lapply(1:3, function(i) setValues(r, runif(ncell(r)))))
f1 <- file.path(tempdir(), "sl1.tif")
writeRaster(s,f1, form
浏览 23提问于2020-01-22得票数 0
回答已采纳
1回答
将路径/路由图写入R中的GeoTiff
r、gis、raster、tiff、geotiff
我正在用geom_path函数绘制R中的一些路线。我希望做的是将我正在绘制的数据转换为GeoTiff (它包括用于投影的GeoSpatial组件& Long角)，可以导入到中。
我所引用的工件可以在这里获得：
我举了一个非常简单的例子来说明我所拥有的，以及我想要做的事情：
library(rgdal)
library(ggplot2)
library(png)
library(raster)
library(tiff)
wrld <- readOGR("data" , "ne_110m_admin_0_countries")
浏览 1提问于2018-12-13得票数 4
回答已采纳
1回答
如何正确地将NCDF文件放入data.frame？
r、dataframe、spatial、cbind、ncdf4
几个星期以来，我一直在尝试将NetCDF文件放入数据帧中。尽管我成功地从ncdf文件中提取了变量/尺寸并绘制了切片，但在将其绑定到dataframe中然后绘制时，所有数据都会被隐藏。这些数据是来自哥白尼的天气数据，包括世界上每个经度和纬度点的数据。这里的最终目标是对数据帧进行栅格，以便能够按栅格对天气进行分类。 可以从https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=overview检索数据。 我的代码如下所示： library(ncdf4)
library(raste
浏览 32提问于2021-10-22得票数 0
1回答
在Rstudio (Mac )上安装"SSOAP包“-获取错误
xml、r、soap、rstudio、rcurl
我希望将我的R脚本连接到SOAP服务器，以便提取数据。我是一个R新手，但一个体面的程序员。
我安装了Rstudio
我安装了RCurl和XML2软件包。
现在我正在尝试安装
SSOAP包我从这里下载了.tar.gz文件：
下面是Rstudio，并使用了"Packages“下的"Install”按钮，请参见此处：
我得到这个错误后：
> library("RCurl", lib.loc="/Library/Frameworks/R.framework/Versions/3.2/Resources/library")
Lo
浏览 2提问于2015-08-30得票数 0
1回答
我试图绘制一个我在ArcGIS中制作的GeoTIFF作为一些地图的背景画布，但我在正确的颜色上遇到了困难。
ArcGIS输出：
但同样的GeoTIFF -使用readGDAL() {rgdal}读入R并使用image()绘制-会产生以下结果：
我尝试过颜色选择--例如，配色方案(rainbow，heat.colors，topo.colors，terrain.colors)是非常有限的，而grey()会产生一个统一的可定制色调的灰色盒子！
非常感谢任何关于如何绘制GeoTIFF的真实颜色的建议，或者其他不太花哨的东西。
(我使用的函数实际上就是我给出的形式--没有额外的参数或
浏览 0提问于2012-07-27得票数 0
回答已采纳
2回答
R光栅:裁剪光栅的上界总是水平的--为什么？
r、r-raster
我试着用一个小的多边形来裁剪一个大的多多边形形状文件。它使用st_intersection工作，但是这需要很长时间，所以我尝试将多多边形转换为栅格，然后用较小的多边形进行裁剪。
## packages - sorry if I've missed any!
library(raster)
library(rgdal)
library(fasterize)
library(sf)
## load files
shp1 <- st_read("pathtoshp", crs = 27700) # a large multipolygon shapefile to c
浏览 7提问于2022-10-26得票数 0
1回答
R：“光栅”不能读取“GeoTIFF”能够读取的readGDAL。
r、r-raster、sp、rgdal
问题
我和MODIS GeoTIFF一起工作。我想使用sf和raster库。这个文件很容易用rgdal::readGDAL()打开，这给了我一个Large SpatialGridDataFrame类型的对象。所以文件没有损坏。
但是，如果使用raster，则会出现以下错误：
raster::raster("../MCD64A1/TIFF/Win09/2014/MCD64monthly.A2014001.Win09.006.burndate.tif")
Error in .rasterObjectFromFile(x, band = band, objecttype = &#
浏览 3提问于2021-10-08得票数 0
1回答
使用R将大型csv文件转换为光栅层的内存问题
r、csv、matrix、raster、r-raster
我有一个大的CSV文件表示的地理光栅地图- 42011列18510行，约10 on在磁盘上。该文件是使用FFT函数作为卷积的输出创建的，其中涉及将原始光栅分解成六条，每条处理一条，并使用重叠添加方法将结果矩阵附加到csv文件中。
我现在试图使用R包GeoTIFF将结果写回光栅文件类型(例如ESRI、等)。CSV文件中的每个数字映射到光栅的等效单元格。
在R中这样做的明显方法是：
library(raster)
library(data.table)
data <- fread("data.csv") ##
浏览 4提问于2014-06-17得票数 1
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
活动推荐
腾讯云开发者
