ArcMapでは、ラスターの一意の値を自動的に計算し、一意の値ごとに異なるスタイルを適用できます(最初の画像を参照)。
ただし、QGISでは、「Singleband pseudocolor」スタイルを使用するときにスタイル設定の値を手動で追加する必要があります(2番目の画像を参照)。ArcMapのように一意の値を自動入力する方法はありますか?
ArcMap:
QGIS(値リストに一意の値を自動入力する方法):
ArcMapでは、ラスターの一意の値を自動的に計算し、一意の値ごとに異なるスタイルを適用できます(最初の画像を参照)。
ただし、QGISでは、「Singleband pseudocolor」スタイルを使用するときにスタイル設定の値を手動で追加する必要があります(2番目の画像を参照)。ArcMapのように一意の値を自動入力する方法はありますか?
ArcMap:
QGIS(値リストに一意の値を自動入力する方法):
回答:
役に立つかもしれないアプローチを以下に示します。GRASSツールr.report
は、ラスターレイヤーの一意の値ごとに基本的な統計を計算できます。
出力は、次のようなプレーンなtxtファイルとして書き込むことができます(左:一意の値、右:各値でカバーされる領域):
+-----------------------------------------------------------------------------+
| Category Information | square|
| #|description | miles|
|-----------------------------------------------------------------------------|
|1111| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 1.249845|
|1113| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 0.062666|
[...]
|3412| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 0.013926|
|4111| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 4.713902|
|4211| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 0.083555|
|4212| . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | 0.135777|
|-----------------------------------------------------------------------------|
|TOTAL | 17.048727|
+-----------------------------------------------------------------------------+
次のジオプロセシングスクリプトは、このテキストファイルを入力として受け取り、ラスターレイヤーに適したレンダラー(便宜上、ランダムな色を使用)を設定します(LAWHEAD、J。に触発:QGIS Python Programming Cookbook、135ページ)。
##giswg=group
##thexml=file
##theraster=raster
from random import randint
from PyQt4.QtGui import QColor
from qgis.core import *
myraster = processing.getObject(theraster)
infile = open(thexml, 'r')
univalues = [l.split('|')[1] for l in [f for f in infile.readlines()][4:-4]]
s = QgsRasterShader()
c = QgsColorRampShader()
c.setColorRampType(QgsColorRampShader.EXACT)
i = []
for u in univalues:
i.append(QgsColorRampShader.ColorRampItem(float(u), QColor('#%06x' % randint(0, 2**24)), u))
c.setColorRampItemList(i)
s.setRasterShaderFunction(c)
ps = QgsSingleBandPseudoColorRenderer(myraster.dataProvider(), 1, s)
myraster.setRenderer(ps)
QGISの伝説的な処理機能のおかげで、これら2つをグラフィカルモデラーにまとめて、ロードされたラスターレイヤーを選択するためにラスターパラメーターを追加できます。
モデルを保存した後、処理ツールボックスでモデルをダブルクリックすると、通常のツールのように使用できます。
結果は次のようになります(ラスタライズされたnatura2000データセットの露頭):
または、レイヤーのプロパティで:
これにはユーザーエクスペリエンスの改善が必要だと思いますが、それは最初のステップです。