Python– tag –
-
日の出・日の入時刻の計算
# 日の出・日の入時間の計算 # http://www.sci-museum.kita.osaka.jp/~egoshi/astronomy/python/python_moon.html import ephem import datetime import pytz # 定義 time_zone = 9 lon = '139.8090057' lat = '35.682998' date = '2023-01-01' # 現地日 horizon = '0:00' # 水平線角度 (D:M:S) # UTC 時刻へ変換 local = pytz.timezone("Asi... -
arcpy の動作を確認する方法
arcpy.Usage関数でジオプロセシングの内容を確認できる。Usage関数が実行できればArcPyが実行できていることになる。 import arcpy help(arcpy.Buffer_analysis) print(arcpy.Usage("Buffer_analysis")) -
ArcPy ユーザーのための Python 基本構文集
ArcPy サイト パッケージを使用するために必要な最低限の Python 基本構文です。ArcGIS Pro 2.1 の Pytyon 3.5、ArcMap の Python 2.7 共通で利用できる構文です。 変数 [python title="Python"] # 動的にデータ タイプが変化 a = 1 print(a) a = 10 print(type(a)) a = "10" print(type(a)) # 日本語や数字、予約語は変数名に使え... -
ArcGIS が認識できないラスター データセットをログに出力
別のロジックで大量に作成したラスター データセット（画像ファイル）が ArcGIS で正しく認識できるかどうかを検品するために作成しました。ファイル ジオデータベースのラスター データセットが認識できないことは無いはずなので対象外です。 [python] #ArcGIS が認識できないラスター データセットをログに出力 #パラメーター #input_dir:... -
Vricon 3D i3s をポータルにアップロードするためのバッチ リネーム処理
ArcGIS Online や Portal for ArGIS にコンピューター上のファイルをアイテムとして登録する際は、サイトにつきユニークなファイル名でなければいけません。 例えば、Vricon 3D という製品の i3s (*.spk) ファイルは、ポータルから登録するアップローダーのファイルサイズ上限がある関係で、購入範囲がタイル上に分割されています。分割され... -
ドーナツ ポリゴンの穴を埋めるジオプロセシング ツール
穴の空いたポリゴンをドーナツ ポリゴンと呼びますが、この穴を埋めるツールの紹介です。ラスターベクター変換を行うと不要な微小のドーナツ ポリゴンが生成されたので、それらを除外したいと思い作成しました。 はじめに ArcMap や ArcGIS Pro の標準 GUI では、1フィーチャずつ図形を編集することでドーナツの穴を埋めることはできますが、... -
ポリゴン フィーチャでクリップしたラスター レイヤーのセル最大値と最小値を取得
久しぶりに GIS らしい処理を ArcPy で実行したのでメモ。今回は、大きな範囲のラスター レイヤーの中にポリゴン フィーチャがいくつか存在し、各ポリゴン フィーチャの範囲内でセルの最大値・最小値を取得するというものです。 [python title="Python"] import arcpy #以下を ArcGIS Pro の Python ウィンドウに貼り付けて実行 arcpy.env.ad... -
2 点間の角度
[python title="Python"] import math radian = math.atan((Y2-Y1)/(X2-X1)) degrees = radian * 180 / math.pi [/python] -
Decimal 型の操作
[python title="Python"] import numpy import decimal x = 198576 / 1000.0 x = numpy.float64(x) y = decimal.Decimal(x).quantize(decimal.Decimal('.000'), rounding=decimal.ROUND_HALF_UP) z = decimal.Decimal(x) print type(y),y print type(z),z [/python] -
VARCHAR 型の日付を変換
日付と時刻の表記も国際規格があり、ISO 8601で決められている。最新版は 2004年の ISO 8601:2004。 https://ja.wikipedia.org/wiki/ISO_8601 これを Python のモジュールで変換するには以下のように記述する [python title="Python"] import dateutil.parser dt = dateutil.parser.parse("YYYY-MM-DDTHH:MM:SSZ") [/python] ただ...
