COLOR(red){SIZE(30){rgdalで地球環境データを見る}} #contents *rgdal とは [#d7f706e6] rgdal はオープンソースの汎用ラスターデータ変換ライブラリ [[GDAL:http://www.remotesensing.org/gdal/]] および同じくオープンソースのベクトルデータ変換ライブラリ [[OGR:http://www.remotesensing.org/gdal/ogr/]] を使用したラスター/ベクトル型GISデータ用パッケージ。GISデータ (USGS の DEM(Digital Elevetaion Model, デジタル標高モデル) ) ばかりでなく、GIF, PNG,JPEG などの一般的な画像データにも対応している。sp パッケージとのリンケージがとれている。また、rmap のインポート・エクスポート技術を取り込んでいる。 [[rgdal(GDAL(Geospatial Data Abstraction Library, 地理空間データ抽象ライブラリ)のバインディング)パッケージ中のオブジェクト一覧]]。 *rgdalの入手先 [#gfcaa7e5] rgdal は CRANで入手できる。 *地球環境データの入手 [#s294488c] -[[GRID 全球データセット:http://www.grid.unep.ch/data/index.php]] [[日本語リスト:http://www-cger.nies.go.jp/grid-j/griddl-j.html]]。 -衛星画像 -DEM(Digital Elevation Model:数値標高モデル) --SRTM(Shuttle Radar Topography Mission)~ 毛利衛さんで有名。~ 分解能約30m~ SRTMのオリジナルフォーマットをサポートしていないが、 以下より GEOTIFF形式および Arc Grid 形式を入手できる。~ http://glcf.umiacs.umd.edu/data/srtm/index.shtml --[[ASTER全球3次元地形データ(G-DEM):http://202.32.29.90/J/index.html]] 日本の経済産業省がNASAと協力して GEOTIFF 形式および GTOPO30 形式で提供。~ 分解能約30m。~ 全地球をカバーの予定。現在[[西日本地域:]]をサンプル提供中。~ [[他のデータとの比較:http://202.32.29.90/J/2.html]]。 --GTOPO30~ USGS(アメリカ地質所)が提供。~ 分解能1km。 -土地被覆 *表示 [#l7a68b4a] - パッケージのロード library(rgdal) - 利用可能なラスターフォーマットの表示 > getGDALDriverNames() name long_name create copy 1 AAIGrid Arc/Info ASCII Grid FALSE TRUE 2 ADRG ARC Digitized Raster Graphics TRUE FALSE 3 AIG Arc/Info Binary Grid FALSE FALSE 4 AirSAR AirSAR Polarimetric Image FALSE FALSE 5 BMP MS Windows Device Independent Bitmap TRUE FALSE 6 BSB Maptech BSB Nautical Charts FALSE FALSE 7 BT VTP .bt (Binary Terrain) 1.3 Format TRUE FALSE 8 CEOS CEOS Image FALSE FALSE 9 COASP DRDC COASP SAR Processor Raster FALSE FALSE 10 COSAR COSAR Annotated Binary Matrix (TerraSAR-X) FALSE FALSE 11 CPG Convair PolGASP FALSE FALSE 12 DIMAP SPOT DIMAP FALSE FALSE 13 DIPEx DIPEx FALSE FALSE 14 DOQ1 USGS DOQ (Old Style) FALSE FALSE 15 DOQ2 USGS DOQ (New Style) FALSE FALSE 16 DTED DTED Elevation Raster FALSE TRUE 17 EHdr ESRI .hdr Labelled TRUE TRUE 18 ELAS ELAS TRUE FALSE 19 ENVI ENVI .hdr Labelled TRUE FALSE 20 ERS ERMapper .ers Labelled TRUE FALSE 21 ESAT Envisat Image Format FALSE FALSE 22 FAST EOSAT FAST Format FALSE FALSE 23 FIT FIT Image FALSE TRUE 24 FujiBAS Fuji BAS Scanner Image FALSE FALSE 25 GenBin Generic Binary (.hdr Labelled) FALSE FALSE 26 GFF Ground-based SAR Applications Testbed File Format (.gff) FALSE FALSE 27 GIF Graphics Interchange Format (.gif) FALSE TRUE 28 GS7BG Golden Software 7 Binary Grid (.grd) FALSE FALSE 29 GSAG Golden Software ASCII Grid (.grd) FALSE TRUE 30 GSBG Golden Software Binary Grid (.grd) TRUE TRUE 31 GSC GSC Geogrid FALSE FALSE 32 GTiff GeoTIFF TRUE TRUE 33 GXF GeoSoft Grid Exchange Format FALSE FALSE 34 HFA Erdas Imagine Images (.img) TRUE TRUE 35 IDA Image Data and Analysis TRUE FALSE 36 ILWIS ILWIS Raster Map TRUE TRUE 37 INGR Intergraph Raster TRUE TRUE 38 ISIS2 USGS Astrogeology ISIS cube (Version 2) FALSE FALSE 39 ISIS3 USGS Astrogeology ISIS cube (Version 3) FALSE FALSE 40 JAXAPALSAR JAXA PALSAR Product Reader (Level 1.1/1.5) FALSE FALSE 41 JDEM Japanese DEM (.mem) FALSE FALSE 42 JPEG JPEG JFIF FALSE TRUE 43 L1B NOAA Polar Orbiter Level 1b Data Set FALSE FALSE 44 LAN Erdas .LAN/.GIS FALSE FALSE 45 Leveller Leveller heightfield TRUE FALSE 46 MEM In Memory Raster TRUE FALSE 47 MFF Vexcel MFF Raster TRUE TRUE 48 MFF2 Vexcel MFF2 (HKV) Raster TRUE TRUE 49 NDF NLAPS Data Format FALSE FALSE 50 NITF National Imagery Transmission Format TRUE TRUE 51 PAux PCI .aux Labelled TRUE FALSE 52 PCIDSK PCIDSK Database File TRUE TRUE 53 PCRaster PCRaster Raster File FALSE TRUE 54 PDS NASA Planetary Data System FALSE FALSE 55 PNG Portable Network Graphics FALSE TRUE 56 PNM Portable Pixmap Format (netpbm) TRUE FALSE 57 RIK Swedish Grid RIK (.rik) FALSE FALSE 58 RMF Raster Matrix Format TRUE FALSE 59 RPFTOC Raster Product Format TOC format FALSE FALSE 60 RS2 RadarSat 2 XML Product FALSE FALSE 61 RST Idrisi Raster A.1 TRUE TRUE 62 SAR_CEOS CEOS SAR Image FALSE FALSE 63 SDTS SDTS Raster FALSE FALSE 64 SGI SGI Image File Format 1.0 FALSE FALSE 65 SRTMHGT SRTMHGT File Format FALSE TRUE 66 Terragen Terragen heightfield TRUE FALSE 67 TSX TerraSAR-X Product FALSE FALSE 68 USGSDEM USGS Optional ASCII DEM (and CDED) FALSE TRUE 69 VRT Virtual Raster TRUE TRUE 70 XPM X11 PixMap Format FALSE TRUE - 利用可能なベクトルフォーマットの表示 > ogrDrivers() name write 1 AVCBin FALSE 2 BNA TRUE 3 CSV TRUE 4 DGN TRUE 5 ESRI Shapefile TRUE 6 GeoJSON TRUE 7 GML TRUE 8 GMT TRUE 9 GPX TRUE 10 KML TRUE 11 MapInfo File TRUE 12 Memory TRUE 13 REC FALSE 14 S57 TRUE 15 SDTS FALSE 16 TIGER TRUE 17 UK .NTF FALSE 18 VRT FALSE - ラスターファイルの読み込み x <- new("GDALReadOnlyDataset", ファイル名) -画像データの要約 summary(画像オブジェクト) - ラスター・データの表示 plot(getPixmapGDAL(画像オブジェクト)) plot(getPixmapGDAL(x)) - 値の参照 - KML 出力 属性値を標高データにできない - PostGIS [[PostGISとの接続例:http://wiki.intamap.org/index.php/PostGIS]] - Cartographic Modelling or Map Algebra -地図投影変換 - ベクトルデータとの重ね合わせ