Lucy’s blog

プログラムで気になったことの覚書

Vimで改行コード^Mを\nに置換

Windowsから送られてきたASCII形式のファイルをVimで見ると、 改行コードが全て^Mとなっていたので、 Macの改行コード\nに置換します。Vimでファイルを開いて以下を普通に入力 :%s/^M/\n/gすると、エラーが返ってきました。 E486: Pattern not found: ^M^Mが…

dbfファイルをcsvに変換

国土数値情報からDLしたシェープファイル形式のデータの中にある、 dbfファイルのデータを使いたくて、 これを一気に(全国分だと170ファイルくらい)csvに変換したいときのシェルスクリプト(bash)。 少数ファイルだけの処理の場合はExcelでファイル形式を…

100mメッシュを緯度経度に変換する

今回利用したデータはこちら 土地利用細分メッシュデータ 参考にしたのはこちらの記事たちです。 標準メッシュの体系とコード white-bear.info 1/10(100m)メッシュコードを緯度経度に変換する関数をつくります。 環境:OSXSierra、Anaconda、python3.6 =====…

緯度経度を用いてmatplotlibで地図をplotする

ざっくりですが覚書です。 環境:OSX,Anaconda,Python3.5 #必要なものをインポート #matplotlibとBasemap (OSXではどちらもpip installで入手可能だった気がします。違ったらスミマセン。。) import pandas as pd import matplotlib.pyplot as pltfrom mat…

測地系のはなし

GoogleMapからある施設の緯度経度を求めて、Pythonでplotするという事案がありました。 そこでGoogleMapの測地系を調べてみたところ、WGS84(世界測地系)となっていました。 WGS84とは何か、JGD2000と違うのか、と思い調べるとこんな記事がありました。 4.…