0

curlパケットのgetBinaryURLを使用して、Webからラスターデータ（png、tiff、...）を受信し、データをディスクに保存します。このラスターデータは、後で空間分析に使用されます。現在、ラスター（およびrgdal）パケットを使用して、ラスターデータをロードし、ラスターオブジェクトを作成しています。

library(RCurl)
url<-"http://demo.mapserver.org/cgi-bin/wms?version=1.1.1&service=WMS&request=GetMap&VERSION=1.1.1&LAYERS=continents&WIDTH=500&HEIGHT=500&BBOX=0,45,15,60&SRS=EPSG:4326&FORMAT=image/png"
map_png = getBinaryURL(url)
#saving the png to disk
map_file <- file("C:/Users/.../map.png", "wb")
writeBin(map_png, map_file)
close(map_file)

#loading png as raster
library(raster)
map <- raster("C:/Users/.../map.png")
#projection and extend
projection(map)<-"+proj=longlat +datum=WGS84 +ellps=WGS84 towgs84=0,0,0"
map@extent@xmin<-0
map@extent@xmax<-15
map@extent@ymin<-45
map@extent@ymax<-60

それは機能しますが、生データの保存と読み込みの部分はあまり良くありません。だから私が好きなのは、「raw」（map_png）オブジェクトを「RasterLayer」（map）オブジェクトに直接渡すことです。このような：

map <- raster(map_png)

これをアーカイブする方法を知っている人はいますか？

-

pngをデコードしたり、packet pngを使用してデコードしたりできることは知っていますが、さまざまな入力形式があるため、これは望ましい方法ではありません。より具体的には、でフェッチした後にすでに利用可能なバイナリオブジェクトを直接操作する時間ゲインの例を示します。getBinaryURL()

system.time(
for(i in 1:200){
map_file <- file("C:/.../map.png", "wb")
writeBin(map_png, map_file)
close(map_file)
map <- raster("C:/.../map.png")
}
)

system.time(
for(i in 1:200){
xx <- readPNG(map_png)
map <- raster(xx[,,2])
}
)
4

1 に答える 1

0

ユーザー、あなたができることがいくつかあります。最良の方法は、呼び出しのmapserver formatパラメーターrgdalをインストールして、GeoTIFF（' ）などの地理的フォーマットに変更することです。パラメータオプションはここにあります。これにより、ジオリファレンスが保持されます。インストールすると、ファイルを非常に簡単に読み取ることができます。Gtiffrgdal

url <- "http://demo.mapserver.org/cgi-bin/wms?version=1.1.1&service=WMS&request=GetMap&VERSION=1.1.1&LAYERS=continents&WIDTH=500&HEIGHT=500&BBOX=0,45,15,60&SRS=EPSG:4326&FORMAT=Gtiff"
download.file(url, destfile="my.tiff")
map <- brick("my.tiff")
map2 <- raster(map, layer=2)

map
class       : RasterBrick 
dimensions  : 500, 500, 250000, 3  (nrow, ncol, ncell, nlayers)
resolution  : 0.03, 0.03  (x, y)
extent      : 0, 15, 45, 60  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : ~/my.tiff 
names       : my.1, my.2, my.3 
min values  :    0,    0,    0 
max values  :  255,  255,  255 

map2
class       : RasterLayer 
band        : 2 
dimensions  : 500, 500, 250000  (nrow, ncol, ncell)
resolution  : 0.03, 0.03  (x, y)
extent      : 0, 15, 45, 60  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : ~/my.tiff 
names       : my.2 
values      : 0, 255  (min, max)

ここに画像の説明を入力してください

png形式を使用する場合は、次のpngパッケージを使用すると簡単になります。

library(png)
library(RCurl)
url <-"http://demo.mapserver.org/cgi-bin/wms?version=1.1.1&service=WMS&request=GetMap&VERSION=1.1.1&LAYERS=continents&WIDTH=500&HEIGHT=500&BBOX=0,45,15,60&SRS=EPSG:4326&FORMAT=image/png"
map_png <- getBinaryURL(url)
xx <- readPNG(map_png)
map <- raster(xx[,,2]) # The png is read as an array, and your data is in the second layer.
于 2013-01-24T21:30:41.547 に答える