#STACK is a raster stackras=STACK1 iii fun summary(lm(x ~ year))$coefficients} A=apply(STACK[iii], 1, fun) ras[iii]=A[1,] STACK[iii].intercept=ras ras[iii]=A[2,] STACK[iii].slope=ras ras[iii]=A[8,] STACK[iii].pval=ras

Raster １バンドのみで構成されるメッシュデータ r=raster("パス/◯◯.tif") Raster stack 複数バンドで構成されるメッシュデータ。brickと違い、データがすべてメモリ上に格納されるわけではないので、大規模データはstackを利用するほうがパンクしない。 s1=…

ここを参照 https://oku.edu.mie-u.ac.jp/~okumura/stat/colors.html #色を指定 #YnGnBlの連続色をかく。Hex codeをチェック brewer.pal(9, "YnGnBl")#色の定義 cols = colorRamp(c("#FFFFD9", "#EDF8B1", "#C7E9B4", "#7FCDBB", "#41B6C4", "#1D91C0", "#22…

今までやってこなかったので。 GRASSのi.groupで簡単に。i.group group=NAME in=image1,image2,image3,..image1,image2,image3,..のところはg.mlistで片付ける。 g.mlist pattern="AAA*" exclude="AAA1" separator=, 一枚のラスタ‐に大量のバンドが入ってる…

QGISでDissolveをするとテーブル内がわけわからないことになる。 （ランダムに行が選択されている？） GRASSでやるとテーブル内が消える。 ArcGISでやると和をとるか平均をとるかなど選択できる。が、選択が異様に面倒くさい。

LandsatのDLにて。 wget –nc –r –l 3 –A tar.gz http://XXX.XXX.XXX wget: Download data from http://XXX.XXX.XXX nc: Skip existing file which already download in the target directory. r: recursive. l 3: Directory level. Download only 3 level fr…

ルートのパーティションとホームのパーティションをddコマンドでimg化してバックアップし、その後リストアするメモ。１．DVDかUSBでUBUNTUを立ち上げる。Try UBUNTUを選択 ２．バックアップするデバイスをマウント 2-1. マウント先の作成 sudo mkdir /media/…

http://m884.jp/RzTutorial_Ja.htmlGUIの起動 library(Rz)

ここを参考に。 元shpデータ内テーブルと結合したいテーブル内にマッチングさせるIDが必要です。まず、もとshpをインポート A_shp=readShapePoly("A.shp", IDvar="OBJECTID")テーブルの確認 A_shp@data summary(A_shp)結合したいテーブルデータをインポート …

ここを参考にしています。 QGISでもPoint Sampling Toolがあるけども、ラスターの数が多いとしんどいですので、Rでやってみました。library(rgdal) library(maptools) library(raster) #ラスタースタックの作成 #実際は１００枚ほどのラスターをスタックrasS…

こちらを参考にしました。HOME=$HOME/【インストール先ディレクトリ】 /usr/bin/dropbox start -iこんだけでオッケっす。

簡単なことでも忘れてしまうのでメモ。tar.sh - #!/bin/shfor i in `ls | grep AA` dotar czvf ${i}.tar.gz ${i} rm ${i} done -