damyarou

python, GMT などのプログラム

2019-01-01から1年間の記事一覧

Python jpg画像を縮小保存しhtmlで表示

タイトル通り、jpg画像を縮小保存しhtmlで表示する。 import glob import os from PIL import Image import os def select_pic(): path = '*.jpg' # 現フォルダ内にあって拡張子がJPGのファイル名を取得 file_list = glob.glob(path) print(file_list) retur…

雑記 また機械学習の本を買ってしまった

記事の最後に行く 昨日(2019.08.11)、また機械学習の本を買ってしまった。 オライリー・ジャパン PythonによりAIプログラミング入門 ディープラーニングを始める前に身に着けておくべき15の基礎技術 あと、久しぶりにpipで一括アップデートしたら、matplot…

雑記 機械学習の勉強を始めてみる

記事の最後に行く いよいよ7月。 昨年(2018年10月10日)、マレーシアでの業務が完了し、帰国した。 それ以降、以下のような感じで過ごしてきた。 2018年10月 ひま 2018年11月 ひま 2018年12月 忙しい 2019年01月 忙しい 2019年02月 忙しい(ベトナム出張2…

設計 Pythonによる水文量頻度解析用自作関数とその利用

記事の最後に行く 自作関数 以下の確率分布に対するパラメータ推定と予測値計算の部分を自作関数としてプログラム本体から分離している。 含まれる関数(確率分布) 関数名説明 ln33変数対数正規分布(LN3) lp3対数ピアソンIII型分布(LP3) gev一般化極値…

FEM 格子桁構造解析プログラム

FEM

記事の最後に行く はじめに 格子桁構造解析プログラムを紹介します。 このプログラムは、平面骨組構造解析プログラムの軸力・軸方向変位をねじりモーメント・ねじり角に書き換えたものです。 構造は平面形状とし、x - y 平面上で定義します。 格子桁解析は、…

設計 RC長方形断面の曲げ耐荷力計算

記事の最後に行く はじめに マレーシア時代に使っていた、RC長方形断面の曲げ耐荷力計算プログラムを掲載しています。 部分安全係数、コンクリートおよび鉄筋の物性は BS 8110-1 に準拠したものとしています。 日本のコンクリート標準示方書での計算も、部分…

設計 確率雨量の算定(4)追補

記事の最後に行く 平方根指数型最大値分布による確率雨量計算 最近ある方から、昔作った文書http://civilyarou.web.fc2.com/WANtaroHP_html5_win/f90_ENGI/dir_HFA/suimon.pdfの平方根指数型最大値分布に関連して、メールを頂いた。 そのころは、C あるいは …

設計 仮排水路設計における Flood Routine の活用(3)

記事の最後に行く ここでは、「設計 仮排水路設計における Flood Routine の活用(2)」で述べた事項に関連する計算・作図プログラムのうち、以下の4つを掲載する。 トンネル標準断面作図プログラム トンネル通水量計算・作図プログラム 貯水池水位容量曲線…

設計 仮排水路設計における Flood Routine の活用(2)

記事の最後に行く ここでは、「設計 仮排水路設計における Flood Routine の活用(2)」で述べた事項に関連する計算・作図プログラムのうち、以下の2つを掲載する。 Flood Routine 解析プログラム Flood Routine 作図プログラム Flood Routine 解析プログ…

設計 仮排水路設計における Flood Routine の活用(1)

記事の最後に行く 概要 仮排水路トンネル断面寸法と仮締切堤高さを決定するための設計検討を行う。 具体的には、トンネル径を変化させ、Flood Routine により、洪水流入波形、流出流量、貯水池水位の関係を求める。 Conditions for Flood Routine Analysis F…

matplotlib RC円形圧力トンネルモデル図

記事の最後に行く はじめに 以下に示す、「設計 RC円形圧力トンネルの配筋設計(1)」で示した図を作成するプログラムである。 ポイント アノーテーション用矢印とテキスト描画 説明用のボックス内テキストと矢印を描画する。 matplotlib では、annotate に…

設計 RC円形圧力トンネルの配筋設計(2)

記事の最後に行く 内水圧を受けるRC圧力トンネル 出力 ra覆工内半径 (mm) rb覆工外半径 (mm) r0岩盤外縁半径位 (mm) ta内側鉄筋等価板厚 (mm) tb外側鉄筋等価板厚 (mm) da内側鉄筋かぶり (mm) db外側鉄筋かぶり (mm) Pa内水圧 (MPa) T温度変化量(マイナスは…

設計 RC円形圧力トンネルの配筋設計(1)

記事の最後に行く 概要 均等な内水圧および外水圧を受けるRC円形圧力トンネルを、トンネル軸方向に均一な厚肉円筒として、平面ひずみ状態でモデル化する。 モデル化においては、以下の考え方を採用した。 圧力水路の構造は鉄筋コンクリート構造とする。 内水…

設計 圧縮材のせん断と引張に対する安全率

記事の最後に行く Safety Factor for Shear Strength and Tensile Strength The definitions of symbols are shown below in this discussion considering compression members such as concrete or rock material. cohesion of material internal friction …

設計 圧縮材のせん断強度と圧縮強度・引張強度の関係

記事の最後に行く Relationship between Shear Strength and Uniaxial Strength The definitions of symbols are shown below in this discussion considering compression members such as concrete or rock material. cohesion of material internal frict…

設計 取水口における渦防止のための被り水深

記事の最後に行く この記事は、取水口における渦防止のための被り水深を提案している以下の論文の提案式(ft系)を、m系に直したものを記載している。 Vortices at intakes by J. L. Gordon (Water Power April 1970) Gordonは、論文の中で、取水口における…

Python フォルダ内のjpgファイルの余白削除

記事の最後に行く フォルダ内の全 jpg ファイルの余白を削除し、同じ名前で保存します。 トリム関数は以下に掲載のものを使わせていただきました。 トリム:https://gist.github.com/mattjmorrison/932345 from PIL import Image, ImageChops import glob, o…

Python フォルダ内のJPGファイルを縮小しhtmlで表示

記事の最後に行く フォルダ内の全JPGファイルを、幅600pixelに縮小し、600_***.jpg という名前で保存します。 また、0_small.html というファイルを作成し、縮小画像を一覧形式で閲覧できるようにします。 import glob import os.path from PIL import Image…

Python 経度・緯度で与えられた2点間距離計算

記事の最後に行く はじめに この記事は、Qiita に投稿した以下の記事の再アレンジ版です。 https://qiita.com/damyarou/items/9cb633e844c78307134a 経度・緯度で与えられた2点間距離を計算するプログラムです。 Lambert-Andoyerの方法を用いています。 デ…

GMT ETOPO1で地形図作成

GMT

記事の最後に行く はじめに この記事は、Qiita に投稿した以下の記事の再アレンジ版です。 https://qiita.com/damyarou/items/b50b6230b103cbab02ba ETOPO1 は,National Geophysical Data Center (Colorado) により作成された,全球の1分刻みの地形モデル…

設計 制水口式調圧水槽のサージング解析

記事の最後に行く はじめに ここで紹介しているのは、制水口式調圧水槽のサージングの基礎微分方程式(連立微分方程式)を Runge-Kutta 法により解き、水面振動の時刻歴を求めるプログラムです。以下のような出力が得られます。 計算理論 Fundamental Differ…

設計 制水口式調圧水槽の基本設計

記事の最後に行く はじめに 制水口式調圧水槽の設計において用いられるVogt-Forchheimer式は、指定した制水口径および立坑径に対し貯水池水位からの最大上昇水位を求めるものであり、制水口径および立坑径をパラメータとして最大上昇水位を計算し、条件を満…

設計 ダムの安定計算

記事の最後に行く はじめに 図に示すような単純形状のダムの安定計算プログラムです。 このような計算は擁壁などの場合も含めよく行うのですが、その都度プログラムを作っていました。 でも簡単な事例をどこかに載せておけば何かと便利なのでここに事例とし…

GMT Windows10にインストールして使う

GMT

記事の最後に行く はじめに GMT (Generic mapping Tools)をWindows10環境にインストール必要が生じたので、その時のメモを残す。 GMTのインストール インストール GMT本家サイト http://gmt.soest.hawaii.edu/projects/gmt/wiki/Download から、gmt-5.4.5-wi…

matplotlib 構造物上にプロット

記事の最後に行く はじめに 描画出力に示すような、構造物の形状を描き、その中にデータをプロットします。 入力データファイルはエクセルで作成しています。 エクセルファイルは、シートを持つことができるので、このようなプロットの入力ファイルとして便…

matplotlib 時系列グラフの軸設定

記事の最後に行く はじめに 時系列グラフの事例を紹介する。 このグラフ作成のポイントは、時間軸の設定方法である。 この記事はQiitaに投稿した以下の記事を再アレンジしたものである。 https://qiita.com/damyarou/items/19f19658b618fd05b3b6 出力画像 プ…

設計 確率雨量の推定(3)プログラミング

記事の最後に行く はじめに 便術記事「Python 確率雨量の推定(2)」で示したデータと、「Python 確率雨量の推定(1)」で示した計算式に基づき、確率雨量の推定を行った。 ここでは、2つの方法で確率雨量の推定を行っている。 各年の日最大雨量を小さい…

設計 確率雨量の推定(2)データ概観

記事の最後に行く はじめに 確率雨量の推定を試験実施するため、札幌・前橋・宮崎の3地点の気象データを、気象庁のホームページより取り出し整理した。 データファイルは空白区切りのテキストファイルであり、列の構成は以下の通りである。 Year : 年 RFy :…

設計 確率雨量の推定(1)計算式

記事の最後に行く はじめに 確率雨量の推定方法を紹介します。 この記述の元ネタは、以下の論文によります。 星清:水文統計解析,開発土木研究所月報,No.540,1998年5月 (https://thesis.ceri.go.jp/db/files/0005005050.pdf) 星清:現場のための水文統…

matplotlib 河川流況表示(2)

記事の最後に行く 概要 河川の流況を示す作図プログラムを紹介します。 対象年数が多くなった場合の事例です。 流況図 縦軸に日平均流量、横軸に超過確率(厳密には超過率)をとった年間の流況を示す図です。 ここでは、1年毎に17年分の線を引き、赤い太…