投稿

Python備忘録 Basic of ShapelyShapely

 Polygonで簡単なリボンを書いてみました。 #@title  デフォルトのタイトル テキスト #Import 'shapely.geometry.Point' class from  shapely.geometry  import  Point points3 =  [] for  i  in   range ( 0 , 2 ):      for  j  in   range ( 0 , 2 ):         points3.append (( i , j )) print ( points3 ) ribon = Polygon ( points3 ) ribon

GIS備忘録 緯度経度⇔ピクセル座標(タイル座標の内部座標)の変換式 2023年1月

イメージ
緯度経度の位置情報と国土地理院の標高タイルを利用して、任意の地点の標高を求めようとすると、緯度経度をピクセル座標に変換する必要に迫られました。 要は、ある緯度経度が、標高タイル(256×256のメッシュ)内のどのメッシュに該当するかを求める必要がございます。 最初は単純に当該タイルの北西端と南東端の緯度経度から、タイルの東西・南北の長さ(あるいは角度)を求めて、それを256で割ってやれば、1メッシュの大きさ(長さor角度の単位)が分かるので、当該タイル内の緯度経度が分かっていれば、 (当該タイルの北西端の北緯)+(当該タイルの北西端の北緯-任意の緯度)÷(1メッシュの緯度の角度増分)=(当該タイルの南北方向メッシュ番号) (当該タイルの北西端の東経)+(当該タイルの北西端の東経-任意の経度)÷(1メッシュの経度の角度増分)=(当該タイルの東西方向メッシュ番号) などのような式で単純に計算できるのかなぁと思っていました。 ただ、そのうち大きめのZoom(すなわち解像度が高い?)で複数(かなり大きな個数の)標高タイルを扱いたくなって色々とやっていると、前述の式ではめちゃめちゃ誤差が出ることが判明しました。 東西方向はあまり誤差がないのですが、南北方向ですごく大きな誤差が発生します。 で、ネット内を探しましたが、これがなかなか式が出てこない。 苦節数日、時間を見つけては検索を続けていたところ、ようやく超難しい換算式を発見しました。なんと、双曲線関数付きの難解な換算式でした。 URLはこちらです。 https://www.trail-note.net/tech/coordinate/ Zoomやタイル座標の説明からはじまり、ピクセル座標、そして、緯度経度の換算式を掲載いただいております。 GISを初心者は必見だと思いました。 (参考)PythonのGeoPandasを活用して、国土地理院の標高タイル(Zoom=12)を合成して、北海道全域の等高線図を書いてみました。

ランニング、グルメ、ドライブ 2023年1月〜

イメージ
2023年3月某日 松山道後温泉ラン風景。 伊予鉄道後温泉駅。坊っちゃん列車が飾られています。駅舎もレトロ。おしゃれなスタバが駅舎に入っています。 大街道のトワイライト。こちらにもスタバが… 道後温泉本館は現在耐震工事中。カラフルな仮設ウォールに囲まれています。手前は一分咲きくらいの桜です。 道後温泉本館の東側は壁が設置されておらず、昔ながらの佇まいを見ることができます。 一番星。マンション群の奥に松山城本丸が小さく見えています。 すこし露出を変えて撮ってみました。 2023年3月某日 道の駅源平の里むれにある食堂じゃこやさんでランチ!いつもはハマチの漬け丼をオーダーするのですが、今回は五剣山になぞらえた天丼、その名も「五剣山」をオーダー。 店員さんに「15分くらいかかりますよー」と念押しがありましたが、この日は時間的余裕があったのでオッケーサイン。 同伴者2名が漬け丼を食べ終わる頃に五剣山が完成。 メニューに載せている写真ほどのクオリティではないですが、高々と天に向かって聳え立つ海老天と茄子天が印象的。白身魚(キスかサヨリ)の天ぷらがとても美味しかったです。 1500円弱で少しお高めですが、ボリューム満点で、家族連れの方々などはシェアするみたいです。お皿はシェア用だそうです。 大満足です。 2023年3月某日 たもやさんで黒カレーうどんをオーダーしました。うどんはいわずもがなのハイクオリティですが、黒カレーが香り高くて、ちょうどよい辛さでめちゃめちゃ美味しかったです。今まではかけうどんとぶっかけうどんがメインでしたが、これは新発見! 2023年3月某日 久しぶりに朝ラン。池畔ラン。快晴です。 2023年3月某日 家庭菜園の脇にある河津桜が満開でした。 2023年3月某日 ジョブチューンを見て、久しぶりに王将に行ってみました。せっかくなので、全員合格のニンニク激増餃子と王将ラーメンを、パトナーは合格の揚げそばオーダー。ラーメンは万人受けする味だと思いました。ニンニク激増餃子はかなりニンニクのパンチが効いています。個人的には揚げそばが好みです。 2023年3月某日 ある春先の晴れた日曜日にランチを松山市駅周辺のランチを検索していたところ、評価の高いトマリギさんを発見。店内は狭いので、行列を覚悟しましたが、13時前で待つことなく、入店できました。ラッキー! 私はボロネー