Monthly Archives: 4月 2017

QGIS #043 – レイヤの地物IDを取得するスクリプト

 

QGISで、レイヤの地物IDを取得するスクリプトは下記のように記述します。

 

Python
1
2
3
4
5
6
layer = iface.activeLayer()
features = layer.getFeatures()
print "features: " + str(features)
 
for f in features:
    print "featureID: " + str(f.id())

 

実行すると選択されたレイヤの地物IDがコンソールに表示されます。
QGIS_043_01

 

レイヤの全地物のオブジェクトを取得:

Python
1
2
features = layer.getFeatures()
print "features: " + str(features)

 

各地物のIDを取得:

Python
1
2
for f in features:
    print "featureID: " + str(f.id())

 

レイヤの地物IDを取得することができます。

 
 

  • 参考文献

QGIS

 

Leaflet #035 – 矩形を表示

 

Leafletで矩形を表示するには「L.rectangle」を利用します。

 

index.html

index.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
<!DOCTYPE html>
<html lang="ja">
 
<head>
    <meta charset="UTF-8">
    <title>Leaflet Sample</title>
 
    <script src="./library/leaflet-0.7.3/leaflet.js"></script>
    <link href="./library/leaflet-0.7.3/leaflet.css" rel="stylesheet" />
 
    <link href="./css/stylesheet.css" rel="stylesheet" />
 
</head>
<body>
 
    <div id="map"></div>
    <script src="./js/script.js"></script>
 
</body>
</html>

 

stylesheet.css

stylesheet.css
CSS
1
2
3
4
5
6
7
8
9
10
html, body {
    height: 100%;
    padding: 0;
    margin: 0;
}
 
#map {
    z-index: 0;
    height: 100%;
}

 

script.js

script.js
JavaScript
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
var m_mono = new L.tileLayer('https://tile.mierune.co.jp/mierune_mono/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var m_color = new L.tileLayer('https://tile.mierune.co.jp/mierune/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var t_pale = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/pale/{z}/{x}/{y}.png', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var t_ort = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/ort/{z}/{x}/{y}.jpg', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var o_std = new L.tileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
    attribution: '&copy; <a href="http://osm.org/copyright">OpenStreetMap</a> contributors'
});
 
var map = L.map('map', {
    center: [35.681,139.763],
    zoom: 14,
    zoomControl: true,
    layers: [m_mono]
});
 
var Map_BaseLayer = {
    "MIERUNE地図 color": m_color,
    "MIERUNE地図 mono": m_mono,
    "地理院地図 淡色": t_pale,
    "地理院地図 オルソ": t_ort,
    "OpenStreetMap 標準": o_std
};
 
var Rectangle = L.rectangle([
    [35.671,139.753],
    [35.691,139.773]
],{
    color: '#FF0000',
    weight: 3,
    opacity: 0.8,
    fillColor: '#FF0000',
    fillOpacity: 0.3
}).addTo(map);
 
L.control.scale({
    imperial: false,
    maxWidth: 300
}).addTo(map);
 
L.control.layers(Map_BaseLayer, null, {
    collapsed: true
}).addTo(map);

 

index.htmlを実行すると下記のようにブラウザで表示されます。
leaflet_035_01

example

 

スタイルを変更：

JavaScript
1
2
3
4
5
6
7
8
9
10
var Rectangle = L.rectangle([
    [35.671,139.753],
    [35.691,139.773]
],{
    color: '#FF0000',
    weight: 10,
    opacity: 0.8,
    fillColor: '#FF0000',
    fillOpacity: 0.1
}).addTo(map);

 

L.rectangleを利用すると地図上に矩形を表示することができます。

 
 

 

QGIS #042 – レイヤの地物を取得するスクリプト

 

QGISで、レイヤの地物を取得するスクリプトは下記のように記述します。

 

Python
1
2
3
4
5
6
layer = iface.activeLayer()
features = layer.getFeatures()
print "features: " + str(features)
 
for f in features:
    print "feature: " + str(f)

 

実行すると選択されたレイヤの地物オブジェクトがコンソールに表示されます。
QGIS_042_01

 

レイヤの全地物のオブジェクトを取得:

Python
1
2
features = layer.getFeatures()
print "features: " + str(features)

 

各地物のオブジェクトを取得:

Python
1
2
for f in features:
    print "feature: " + str(f)

 

レイヤ内にある地物を取得することができます。

 
 

  • 参考文献

QGIS

 

Python #002 – Pyminifierでファイル圧縮

 

Pythonでファイル圧縮・コメント削除・難読化が手軽にできないか探してみました。
Webのフロントエンド(HTML・CSS・JS)とかだとgulp等を利用してファイル圧縮が手軽にできますがPythonだとあまり見当たりません。。。

Pythonでファイル圧縮したい時はPyminifierを利用するとよさそうです。

 

Pyminifierはpip(パッケージ管理システム)でインストールすることができます。

Shell
1
sudo pip install pyminifier

Python_002_01

 

インストールできているか確認してみます。

Shell
1
pyminifier --version

Python_002_02

 

今回は、以前作成したJapanElevationPluginというプラグインを材料にファイル圧縮してみます。

下記コマンドを実行すると、指定フォルダに圧縮された全てのファイルがコピーされます。

Shell
1
pyminifier --use-tabs --destdir=../dist/ ./JapanElevationPlugin/*.py

Python_002_03

Python_002_04

 

下記コマンドを実行すると、指定フォルダに圧縮された全てのファイル(スペースインデント)がコピーされます。

Shell
1
pyminifier --destdir=../dist2/ ./JapanElevationPlugin/*.py

Python_002_05

Python_002_06

 

下記コマンドを実行すると、指定フォルダに難読化された全てのファイルがコピーされます。

Shell
1
pyminifier --obfuscate --use-tabs --destdir=../dist3/ ./JapanElevationPlugin/*.py

Python_002_07

Python_002_08

 

Pyminifier利用するとPythonでファイル圧縮やコメント削除や難読化が手軽にできます。

 
 

  • 参考文献

Pyminifier

 

Leaflet #034 – 固定値のサークルを表示

 

Leafletで固定値のサークルを表示するには「L.circleMarker」を利用します。

 

index.html

index.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
<!DOCTYPE html>
<html lang="ja">
 
<head>
    <meta charset="UTF-8">
    <title>Leaflet Sample</title>
 
    <script src="./library/leaflet-0.7.3/leaflet.js"></script>
    <link href="./library/leaflet-0.7.3/leaflet.css" rel="stylesheet" />
 
    <link href="./css/stylesheet.css" rel="stylesheet" />
 
</head>
<body>
 
    <div id="map"></div>
    <script src="./js/script.js"></script>
 
</body>
</html>

 

stylesheet.css

stylesheet.css
CSS
1
2
3
4
5
6
7
8
9
10
html, body {
    height: 100%;
    padding: 0;
    margin: 0;
}
 
#map {
    z-index: 0;
    height: 100%;
}

 

script.js

script.js
JavaScript
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
var m_mono = new L.tileLayer('https://tile.mierune.co.jp/mierune_mono/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var m_color = new L.tileLayer('https://tile.mierune.co.jp/mierune/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var t_pale = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/pale/{z}/{x}/{y}.png', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var t_ort = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/ort/{z}/{x}/{y}.jpg', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var o_std = new L.tileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
    attribution: '&copy; <a href="http://osm.org/copyright">OpenStreetMap</a> contributors'
});
 
var map = L.map('map', {
    center: [35.681,139.763],
    zoom: 14,
    zoomControl: true,
    layers: [m_mono]
});
 
var Map_BaseLayer = {
    "MIERUNE地図 color": m_color,
    "MIERUNE地図 mono": m_mono,
    "地理院地図 淡色": t_pale,
    "地理院地図 オルソ": t_ort,
    "OpenStreetMap 標準": o_std
};
 
var CircleMarker = L.circleMarker([35.681,139.763], {
    color: '#FF0000',
    weight: 3,
    opacity: 0.8,
    fillColor: '#FF0000',
    fillOpacity: 0.3,
    radius: 50
}).addTo(map);
 
L.control.scale({
    imperial: false,
    maxWidth: 300
}).addTo(map);
 
L.control.layers(Map_BaseLayer, null, {
    collapsed: true
}).addTo(map);

 

index.htmlを実行すると下記のようにブラウザで表示されます。
leaflet_034_01

example

 

スタイルを変更：

JavaScript
1
2
3
4
5
6
7
var CircleMarker = L.circleMarker([35.681,139.763], {
    color: '#001AFF',
    weight: 6,
    opacity: 0.8,
    fillColor: '#001AFF',
    fillOpacity: 0.7
}).addTo(map);

 

L.circleMarkerを利用すると地図上に固定値のサークルを表示することができます。

 
 

 

QGIS #041 – レイヤの地物を全て選択解除するスクリプト

 

QGISで、レイヤの地物を全て選択解除するスクリプトは下記のように記述します。

 

Python
1
2
layer = iface.activeLayer()
layer.removeSelection()

 

実行すると選択したレイヤの地物を全て選択解除されます。
QGIS_041_01

 

レイヤの地物を全て選択解除:

Python
1
layer.removeSelection()

 

レイヤの地物を全て選択解除することができます。

 
 

  • 参考文献

QGIS

 

Leaflet #033 – サークルを表示

 

Leafletでサークルを表示するには「L.circle」を利用します。

 

index.html

index.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
<!DOCTYPE html>
<html lang="ja">
 
<head>
    <meta charset="UTF-8">
    <title>Leaflet Sample</title>
 
    <script src="./library/leaflet-0.7.3/leaflet.js"></script>
    <link href="./library/leaflet-0.7.3/leaflet.css" rel="stylesheet" />
 
    <link href="./css/stylesheet.css" rel="stylesheet" />
 
</head>
<body>
 
    <div id="map"></div>
    <script src="./js/script.js"></script>
 
</body>
</html>

 

stylesheet.css

stylesheet.css
CSS
1
2
3
4
5
6
7
8
9
10
html, body {
    height: 100%;
    padding: 0;
    margin: 0;
}
 
#map {
    z-index: 0;
    height: 100%;
}

 

script.js

script.js
JavaScript
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
var m_mono = new L.tileLayer('https://tile.mierune.co.jp/mierune_mono/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var m_color = new L.tileLayer('https://tile.mierune.co.jp/mierune/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var t_pale = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/pale/{z}/{x}/{y}.png', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var t_ort = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/ort/{z}/{x}/{y}.jpg', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var o_std = new L.tileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
    attribution: '&copy; <a href="http://osm.org/copyright">OpenStreetMap</a> contributors'
});
 
var map = L.map('map', {
    center: [35.681,139.763],
    zoom: 14,
    zoomControl: true,
    layers: [m_mono]
});
 
var Map_BaseLayer = {
    "MIERUNE地図 color": m_color,
    "MIERUNE地図 mono": m_mono,
    "地理院地図 淡色": t_pale,
    "地理院地図 オルソ": t_ort,
    "OpenStreetMap 標準": o_std
};
 
var Circle = L.circle([35.681,139.763], 500,{
    color: '#FF0000',
    weight: 3,
    opacity: 0.8,
    fillColor: '#FF0000',
    fillOpacity: 0.3
}).addTo(map);
 
L.control.scale({
    imperial: false,
    maxWidth: 300
}).addTo(map);
 
L.control.layers(Map_BaseLayer, null, {
    collapsed: true
}).addTo(map);

 

index.htmlを実行すると下記のようにブラウザで表示されます。
leaflet_033_01

example

 

スタイルを変更：

JavaScript
1
2
3
4
5
6
7
var Circle = L.circle([35.681,139.763], 500,{
    color: '#001AFF',
    weight: 6,
    opacity: 0.8,
    fillColor: '#001AFF',
    fillOpacity: 0.7
}).addTo(map);

 

Leafletのver1以上の記述方法：

JavaScript
1
2
3
4
5
6
7
8
var Circle = L.circle([35.681,139.763], {
    color: '#FF0000',
    weight: 3,
    opacity: 0.8,
    fillColor: '#FF0000',
    fillOpacity: 0.3,
    radius: 500
}).addTo(map);

 

L.circleを利用すると地図上にサークルを表示することができます。

 
 

 

QGIS #040 – レイヤのフィールド名を取得するスクリプト

 

QGISで、レイヤのフィールド名を取得するスクリプトは下記のように記述します。

 

Python
1
2
3
4
5
6
layer = iface.activeLayer()
fields = layer.pendingFields()
print "fields: " + str(fields)
 
for f in fields:
    print "fieldname: " + str(f.name())

 

実行すると選択されたレイヤのフィールド名がコンソールに表示されます。
QGIS_040_01

 

全てのレイヤのフィールドをオブジェクトで取得:

Python
1
2
fields = layer.pendingFields()
print "fields: " + str(fields)

 

.name()でフィールド名を取得:

Python
1
2
for f in fields:
    print "fieldname: " + str(f.name())

 

レイヤのフィールド名を調べる事ができます。

 
 

  • 参考文献

QGIS

 

Leaflet #032 – ポリゴンを表示

 

Leafletでポリゴンを表示するには「L.polygon」を利用します。

 

index.html

index.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
<!DOCTYPE html>
<html lang="ja">
 
<head>
    <meta charset="UTF-8">
    <title>Leaflet Sample</title>
 
    <script src="./library/leaflet-0.7.3/leaflet.js"></script>
    <link href="./library/leaflet-0.7.3/leaflet.css" rel="stylesheet" />
 
    <link href="./css/stylesheet.css" rel="stylesheet" />
 
</head>
<body>
 
    <div id="map"></div>
    <script src="./js/script.js"></script>
 
</body>
</html>

 

stylesheet.css

stylesheet.css
CSS
1
2
3
4
5
6
7
8
9
10
html, body {
    height: 100%;
    padding: 0;
    margin: 0;
}
 
#map {
    z-index: 0;
    height: 100%;
}

 

script.js

script.js
JavaScript
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
var m_mono = new L.tileLayer('https://tile.mierune.co.jp/mierune_mono/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var m_color = new L.tileLayer('https://tile.mierune.co.jp/mierune/{z}/{x}/{y}.png', {
    attribution: "Maptiles by <a href='http://mierune.co.jp/' target='_blank'>MIERUNE</a>, under CC BY. Data by <a href='http://osm.org/copyright' target='_blank'>OpenStreetMap</a> contributors, under ODbL."
});
 
var t_pale = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/pale/{z}/{x}/{y}.png', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var t_ort = new L.tileLayer('http://cyberjapandata.gsi.go.jp/xyz/ort/{z}/{x}/{y}.jpg', {
    attribution: "<a href='http://www.gsi.go.jp/kikakuchousei/kikakuchousei40182.html' target='_blank'>国土地理院</a>"
});
 
var o_std = new L.tileLayer('http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
    attribution: '&copy; <a href="http://osm.org/copyright">OpenStreetMap</a> contributors'
});
 
var map = L.map('map', {
    center: [35.681,139.763],
    zoom: 14,
    zoomControl: true,
    layers: [m_mono]
});
 
var Map_BaseLayer = {
    "MIERUNE地図 color": m_color,
    "MIERUNE地図 mono": m_mono,
    "地理院地図 淡色": t_pale,
    "地理院地図 オルソ": t_ort,
    "OpenStreetMap 標準": o_std
};
 
var Polygon = L.polygon([
    [35.675949251235025,139.76617813110352],
    [35.67410157813001,139.77188587188718],
    [35.67455478492641,139.77227210998535],
    [35.683757812281115,139.77862358093262],
    [35.68431553740134,139.77343082427979],
    [35.68469897115985,139.77094173431396],
    [35.679923346539084,139.76871013641357],
    [35.675949251235025,139.76617813110352]
],{
    "color": "#E92D63",
    "weight": 3,
    "opacity": 0.8,
    "fillColor": "#562DE9",
    "fillOpacity": 0.4
}).addTo(map);
 
L.control.scale({
    imperial: false,
    maxWidth: 300
}).addTo(map);
 
L.control.layers(Map_BaseLayer, null, {
    collapsed: true
}).addTo(map);

 

index.htmlを実行すると下記のようにブラウザで表示されます。
leaflet_032_01

example

 

スタイルを変更：

JavaScript
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
var Polygon = L.polygon([
    [35.675949251235025,139.76617813110352],
    [35.67410157813001,139.77188587188718],
    [35.67455478492641,139.77227210998535],
    [35.683757812281115,139.77862358093262],
    [35.68431553740134,139.77343082427979],
    [35.68469897115985,139.77094173431396],
    [35.679923346539084,139.76871013641357],
    [35.675949251235025,139.76617813110352]
],{
    "color": "#E200FF",
    "weight": 1,
    "opacity": 0.2,
    "fillColor": "#E200FF",
    "fillOpacity": 0.2
}).addTo(map);

 

L.polygonを利用すると地図上にポリゴンを表示することができます。

 
 

 

QGIS #039 – レイヤの地物の種類を取得するスクリプト

 

QGISで、レイヤの地物の種類を取得するスクリプトは下記のように記述します。

 

Python
1
2
3
4
layer = iface.activeLayer()
g = layer.geometryType()
 
print "geometryType: " + str(g)

 

実行すると選択されたレイヤの地物の種類がコンソールに表示されます。
QGIS_039_02

 

レイヤの地物の種類を取得する場合：

Python
1
layer.geometryType()

0 = Point
1 = LineString
2 = Polygon

 

レイヤの地物の種類を調べる事ができます。

 
 

  • 参考文献

QGIS