ArcGISのarcpyで断面線を作成するには stack_Profile()を使用する。 このジオプロツールを使用すれば、断面の標高と始点からの距離をテーブルに出すこともできるし、断面図を作成できる（※1）。必要なデータはDEMと断面線だけ。 もちろん、複数の線にも対応…

ArcGISのフォーカル統計を使えば、ラスタの近傍統計解析ができる。 近傍統計解析ができれば、周囲の平均をとったり、周囲の得点を集計して真ん中の得点にする、ってことが楽にできる。 ArcGIS Help 特にこの関数（ツール）のよいところが、neighborhoodパラ…

実はnumpyのインデックス指定は、float型（浮動小数点）も使えました！ ってお話。 自動で型変換してくれていました。。。 こんな時はエラーを返してほしいよね…。 まぁ、変数管理のできていなかった俺が悪いんだけど。 以下、詳細 numpyはインデックス指定…

ラスタをSpatial analysysのツールを使って演算しようとすると、エラーになる。 エラー例） ・テーブルが見つかりません。と言われる。 ・（python上で）import文の箇所で、"ASCII codec can't encode ...."と言われる。 などなど… 調べてみても、 原因がよ…

python 四則演算する上での注意 pythonで同じ型同士の演算をすると、返り値の型は入力と同じになる。 3 * 4 int * int >> 12 int 1.0 + 2.0 float * float >> 3.0 float そのため、整数同士で割り算をすると... 29 / 2 int * int >> 14 int ！！！ 切り捨て…

先日、磐梯山に行ってきた。 素晴らしい山だった。 何が良かったのか？ 荒々しい磐梯山、赤く染まる磐梯山、浅葱色の絨毯、深い森が見れたからだと思う。 一つのルートでここまで景色が変わるところは、そうはない。 もし自分の百名山を決めるなら、ここは必…

背景>> 先日、こんなようなコードを書いた。 import, numpy,arcpy np = numpy.arange(9).reshape([3,3]) raster2 = arcpy.sa.Con(raster1 > 10, np[2,2], float(u"NaN")) するとRuntimeだかジオプロセシングエラー999999だかが出て、原因不明で落ちる。 原因…

システムドライブを入れているHDDを交換する方法 Cドライブの入ったHDDやSSDを交換する方法メモ システムイメージの作成方法 121ware.com > サービス＆サポート > Q&A > Q&A番号 015953 古いシステムイメージは削除できないか？ Windows8.1のシステムイメー…

○オープンソース or フリーウェア アイダホ州立大学のページ LiDAR Tools | Idaho LiDAR Consortium ポイントクラウドライブラリ http://derivecv.tumblr.com/post/10198738207 ○おまけ アジア航測 >> 何か専用ビューアー 中日本航空 >> Mierrun

TODO管理するツールを作ったとする。 テキストボックスなどにやる事や締め切りを記述して保存ボタンを押すとデータベースに登録されるようなツール。 保存するとき、入力内容のエラーチェックをすべきだと思うが、その結果はどのように返すべきか？ ①エラー…

GUIのツールを作る際、入力した値を取り込みたい場合がある。 そんな時どうすればいいか？ ○entry widgetの場合 self.<entry>.get() ・・・引数なし。 ○text widgetの場合 sel.<text>.get(<index_start>,<index_end>=None) ・・・どこからどこまで取得したいかを指定する。 ※indexとは ・ドキュメ</index_end></index_start></text></entry>…

mappingモジュールを用い、レイヤオブジェクトのデータソースを変更しようとした。 データソースをSDEファイルのデータに変更しようとしたらエラーが出たのでメモ。 レイヤのデータソースを変更するには、LayerオブジェクトのreplaceDataSourceを用いる。 ■…

arcpy.CheckProductを実行した際、 arcpy.CheckProduct("arcInfo") と旧名で指定するとエラーが返る。 arcpy.CheckProduct("arcAdvanced") と新名で指定しないと通らない。 ちなみに、standord（arcEditor）、basic（arcview）は旧名で指定しないと通らない。

○処理の一時中断 timeモジュールのsleep関数を用いる。 import time time.sleep(100) で処理が一時停止し、100秒後に再開される。 ○処理の入力待ち 組み込み関数のraw_input関数を用いる。 raw_input("何か入力してください >>") で、何か入力するまで処理が…

DeleteFields_management() の第二引数はフィールドオブジェクトではなく、フィールドの名前を指定する。 フィールドオブジェクトを指定すると、なんか変なRuntimeErrorが返ってくる。 ホームページには引数はフィールドオブジェクトと書いてあるし、なんじ…

■リストに整数をappendする。 >>> list=[] >>> list.append(1) >>> print(list) [1] リストが一つ増えた。 ■リストに整数をextendする。 >>> list.extend(2) Traceback (most recent call last): File "<interactive input>", line 1, in <module> TypeError: 'int' object is not itera</module></interactive>…

変数はどこまで適用されるのか ■関数の外部で定義した変数は、関数内で使用できるのか？ →使用可能。 def function(text): print(aaa) print(text) aaa="HELLO" function("WORLD") >>> HELLO WORLD ■内部でも同じ変数を定義していた場合はどうなるか？ →内部…

arcpyで投影変換を行う。出力にはあるGDB下のデータセットを指定する。 すると、データセット下にではなく、GDB下に投影変換が出力される。 なぜか？

どうすればプロパティを変更できるのか？ ## #空間参照のプロパティを変更する方法 ## #しかし、距離単位が変わらない・・・ ## sr.linearUnitName="Meter" ## sr.linearUnitCode=9001 ## ## print"*******変更結果********" ## print sr.linearUnitCode ## …

プログラムのテスト仕様を組み立てるうえで、何に注意すればいいのか？ 項目立ての注意点 ・仕様に沿っているか？ ・出力データはどんなデータであって欲しいか？ ・入力される可能性のあるデータはどんなデータか？ シナリオを作成する上での注意点 ・誰が…

距離単位が定義されていない時のXY座標精度 →投影座標が定義されていないshpをフィーチャクラスに変換した時、単位がメートルであるとしてXY座標精度が扱われる。 例） Point = arcpy.Point(12.34567,89.01234)を作成 arcpy.Describe(Point).spatialReferenc…

試しに書いてみる