TWI465721B - A method for surface velocity estimation of river water - Google Patents
A method for surface velocity estimation of river water Download PDFInfo
- Publication number
- TWI465721B TWI465721B TW102101646A TW102101646A TWI465721B TW I465721 B TWI465721 B TW I465721B TW 102101646 A TW102101646 A TW 102101646A TW 102101646 A TW102101646 A TW 102101646A TW I465721 B TWI465721 B TW I465721B
- Authority
- TW
- Taiwan
- Prior art keywords
- river
- water
- actual length
- flow rate
- image data
- Prior art date
Links
Landscapes
- Measuring Volume Flow (AREA)
- Image Analysis (AREA)
Description
本發明係關於一種河水表面之流速估算方法,尤其是利用河水影像以進行流速估算之方法。The invention relates to a method for estimating the flow velocity of a river surface, in particular to a method for estimating a flow velocity by using a river water image.
河水表面之流速快慢資訊,與該河水之河床結構及水流量具有高度的相關性,亦即當該河水表面之流速較快時,可研判該河水之河床較淺或具有較高之傾斜度,且在相同之水流截面下,亦會具有較高的水流量。因此,在進行河床地勢判斷或河水流量估測時,河水表面的流速快慢資訊為相當重要的參考資料。The velocity information on the surface of the river has a high correlation with the structure of the river bed and the water flow, that is, when the flow velocity of the river surface is fast, the river bed of the river can be judged to be shallow or have a high inclination. And under the same water flow cross section, it will also have a higher water flow rate. Therefore, the information on the velocity of the river surface is a very important reference when conducting riverbed topographic judgment or river flow estimation.
對於河水表面之流速估算,習知係於一河道中投擲一漂流物,使該漂流物自該河道之一上游端漂流至一下游端,並統計該漂流物所行經之距離與所經歷之時間,再將該距離除以時間以得到該河水表面的流速。惟,該漂流物由上游端漂流至下游端的過程中,可能受到河水高低或河床礫石的影響,使該漂流物擱淺於河岸或是礫石群中,造成該次作業無效,無法準確估算該河水表面之流速。For the estimation of the flow velocity on the surface of the river, it is customary to throw a drifting object in a river channel, so that the drifting object drifts from the upstream end of the river channel to a downstream end, and the distance traveled by the drifting object and the elapsed time are counted. And divide the distance by the time to get the flow rate of the river surface. However, the drifting of the drift from the upstream end to the downstream end may be affected by the height of the river or the gravel of the river bed, causing the drifting object to hang on the river bank or the gravel group, causing the operation to be invalid and the surface of the river water cannot be accurately estimated. The flow rate.
另,河水表面之流速估算,習知亦可將一感應元件設置於靠近河水表面之流道中,以偵測該河水表面之流速。惟,該感應元件係感應位於該設置處之河道的河水表面流速,而該流速僅能代表該河道之一處的河水表面流速,無法代表該河道之一區段的河水表面流速,若欲得到該河道之一區段的河水表面流速,則必須於該河道中設置多個感應元件,此舉 將造成硬體成本上升。In addition, the flow velocity estimation on the surface of the river can also be performed by placing an inductive element in the flow channel near the surface of the river to detect the flow velocity of the surface of the river. However, the sensing element senses the flow velocity of the river surface of the river at the setting, and the flow velocity can only represent the surface velocity of the river water at one of the river channels, and cannot represent the surface velocity of the river water in one section of the river channel, if desired The surface velocity of the river in one section of the channel must have multiple sensing elements in the channel. Will cause hardware costs to rise.
本發明之主要目的係提供一種河水表面之流速估算方法,該方法可準確的估算一河道區域之河水表面流速。The main object of the present invention is to provide a method for estimating the flow rate of a river surface, which can accurately estimate the surface velocity of a river channel in a river channel region.
本發明之另一目的係提供一種河水表面之流速估算方法,該方法在進行河道區域之河水表面流速時,具有較低之硬體設置成本。Another object of the present invention is to provide a method for estimating the flow rate of a river surface which has a lower hardware installation cost when performing a river surface velocity in a river region.
為達到前述發明目的,本發明河水表面之流速估算方法係包含一影像比例計算步驟,係由一處理器讀取一影像資料,並依該影像資料計算得一像素所代表之實際長度;一傾斜角度計算步驟,係由該處理器依該單位像素所代表之實際長度,計算該影像資料中之一水面線、一基準線及一垂直高度的實際長度,並將該基準線及該垂直高度的實際長度代入一角度計算方程式,以計算一傾斜角度;一水流時間計算步驟,係由該處理器將該傾斜角度及該垂直高度的實際長度代入一時間計算方程式,以計算一水流時間;及一表面流速計算步驟,係由該處理器將該水面線的實際長度除以該水流時間,以得到一表面流速。In order to achieve the foregoing object, the method for estimating the flow rate of the surface of the river of the present invention comprises an image ratio calculation step of reading an image data by a processor and calculating an actual length represented by a pixel according to the image data; The angle calculation step is performed by the processor, according to the actual length represented by the unit pixel, an actual length of a water surface line, a reference line and a vertical height in the image data, and the reference line and the vertical height are calculated. The actual length is substituted into an angle calculation equation to calculate an inclination angle; a water flow time calculation step is performed by the processor by substituting the inclination angle and the actual length of the vertical height into a time calculation equation to calculate a water flow time; The surface flow rate calculation step is performed by the processor dividing the actual length of the water surface line by the water flow time to obtain a surface flow rate.
本發明河水表面之流速估算方法,其中上述影像資料之水面線的二端分別為一上游端及一下游端,該基準線係平行於地平線,該基準線與該水面線之間形成該傾斜角度,該基準線與該上游端形成該垂直高度。The method for estimating the flow rate of the surface of the river according to the present invention, wherein the two ends of the surface line of the image data are an upstream end and a downstream end, respectively, the reference line is parallel to the horizon, and the inclination angle is formed between the reference line and the water surface line. The reference line forms the vertical height with the upstream end.
本發明河水表面之流速估算方法,其中上述角度計算方程式如下所示:
其中θ 為該傾斜角度,tan-1 (‧)為一反正切函數,h為該垂直高度H的實際長度,b為該基準線B的實際長度。Where θ is the angle of inclination, tan -1 (‧) is an arctangent function, h is the actual length of the vertical height H, and b is the actual length of the reference line B.
本發明河水表面之流速估算方法,其中上述時間計算方程式
如下所示:
其中t為該水流時間,sin(‧)為一正弦函數,g為重力加速度值,θ 為該傾斜角度,h為該垂直高度H的實際長度。Where t is the water flow time, sin(‧) is a sinusoidal function, g is the gravitational acceleration value, θ is the inclination angle, and h is the actual length of the vertical height H.
本發明河水表面之流速估算方法,其中上述影像資料係由一影像擷取裝置拍攝一河道之影像,且該影像擷取裝置拍攝該河道之拍攝方向與該河道之水流方向呈垂直。The method for estimating the flow rate of the surface of the river according to the present invention, wherein the image data is captured by an image capturing device, and the image capturing device shoots the direction of the river perpendicular to the direction of the river.
1‧‧‧影像擷取裝置1‧‧‧Image capture device
2‧‧‧處理器2‧‧‧ Processor
A‧‧‧水面線A‧‧‧ surface line
A1‧‧‧上游端A1‧‧‧ upstream end
A2‧‧‧下游端A2‧‧‧ downstream end
B‧‧‧基準線B‧‧‧ baseline
H‧‧‧垂直高度H‧‧‧Vertical height
θ‧‧‧傾斜角度Θ‧‧‧ tilt angle
S1‧‧‧影像比例計算步驟S1‧‧‧image scale calculation step
S2‧‧‧傾斜角度計算步驟S2‧‧‧ Tilt angle calculation step
S3‧‧‧水流時間計算步驟S3‧‧‧Water flow time calculation steps
S4‧‧‧表面流速計算步驟S4‧‧‧ Surface flow rate calculation steps
第1圖:本發明之較佳實施裝置圖。Figure 1 is a view of a preferred embodiment of the present invention.
第2圖:本發明之影像資料示意圖。Figure 2: Schematic diagram of the image data of the present invention.
第3圖:本發明河水表面之流速估算方法流程圖。Figure 3 is a flow chart showing the method for estimating the flow rate of the surface of the river according to the present invention.
為讓本發明之上述及其他目的、特徵及優點能更明顯易懂,下文特舉本發明之較佳實施例,並配合所附圖式,作詳細說明如下:本發明所述之「上游端」及「下游端」,係根據一水流流經一河道的方向而定,亦即當該水流由該河道之第一端流至第二端時,該第一端即為「上游端」,該第二端即為「下游端」。The above and other objects, features and advantages of the present invention will become more <RTIgt; And the "downstream end" is determined by the direction in which a stream of water flows through a channel, that is, when the stream flows from the first end to the second end of the channel, the first end is the "upstream end". The second end is the "downstream end".
本發明所述之「像素」(pixels),係指組成一影像(image)資料之最小單位,並可用以表示該影像資料之解析度(resolution)。例如該 影像資料之解析度為1024×768,則代表該影像資料共有(1024×768=786432)個像素。The "pixels" as used in the present invention refer to the smallest unit constituting an image data and can be used to indicate the resolution of the image data. For example The resolution of the image data is 1024×768, which means that the image data has a total of (1024×768=786432) pixels.
請參照第1圖所示,其係本發明河水表面之流速估算方法之一較佳實施裝置,該裝置包含一影像擷取裝置1及一處理器2。Please refer to FIG. 1 , which is a preferred embodiment of a method for estimating a flow rate of a river surface according to the present invention. The device comprises an image capturing device 1 and a processor 2 .
該影像擷取裝置1係用以拍攝一河道之影像,以產生一影像資料。該影像擷取裝置1可為一般攝影機、網路攝影機或紅外線攝影機等裝置。The image capturing device 1 is configured to capture an image of a river channel to generate an image data. The image capturing device 1 can be a device such as a general camera, a network camera, or an infrared camera.
請參照第2圖所示,其係本發明之影像資料示意圖。在本實施例中,該影像擷取裝置1拍攝該河道之一拍攝方向,較佳係與該河道之水流方向呈垂直,以使該影像資料包含該河道之一水面線A,且該水面線A的二端分別為一上游端A1及一下游端A2。該影像資料另具有一基準線B,且該基準線B係平行於地平線,在本實施例中,該基準線B係為該影像資料之底緣。該基準線B與該水面線A之間形成一傾斜角度θ ,該基準線B與該上游端A1具有一垂直高度H。Please refer to FIG. 2, which is a schematic diagram of the image data of the present invention. In this embodiment, the image capturing device 1 captures a shooting direction of the river channel, preferably perpendicular to the water flow direction of the river channel, so that the image data includes a water surface line A of the river channel, and the water surface line The two ends of A are an upstream end A1 and a downstream end A2, respectively. The image data further has a reference line B, and the reference line B is parallel to the horizon. In the embodiment, the reference line B is the bottom edge of the image data. An inclination angle θ is formed between the reference line B and the water surface line A, and the reference line B has a vertical height H with the upstream end A1.
該處理器2可為一電腦或任何運算處理器,且電性連接該影像擷取裝置1,以接收該影像資料,並根據該影像資料進行運算統計等操作。The processor 2 can be a computer or any computing processor, and is electrically connected to the image capturing device 1 to receive the image data, and perform operations such as calculation and statistics according to the image data.
請參照第3圖所示,本發明所述之河水表面之流速估算方法係包含一影像比例計算步驟S1、一傾斜角度計算步驟S2、一水流時間計算步驟S3及一表面流速計算步驟S4。Referring to FIG. 3, the method for estimating the flow rate of the river surface according to the present invention comprises an image ratio calculation step S1, an inclination angle calculation step S2, a water flow time calculation step S3, and a surface flow rate calculation step S4.
該影像比例計算步驟S1,係由該處理器2讀取該影像資料,並依該影像資料計算得一單位像素所代表之實際長度。In the image ratio calculation step S1, the image data is read by the processor 2, and the actual length represented by one unit pixel is calculated according to the image data.
更詳言之,該處理器2讀取該影像資料後,可得到該影像資料之總像素值。接著,該處理器2便可計算得知該影像資料之該單位像素與實際景色之一影像比例,以供換算該影像資料之該單位像素所代表之實 際長度。例如,假設該影像擷取裝置1所拍攝之該影像資料中包含該河道,已知該影像資料中之河道實際長度為500公分,若該處理器2讀取該影像資料,並計算得該影像資料中之河道包含100個像素時,即可推算該單位像素代表5公分,亦即該影像比例為5(cm/pixel)。In more detail, after the processor 2 reads the image data, the total pixel value of the image data can be obtained. Then, the processor 2 can calculate the ratio of the unit pixel of the image data to the actual scene, so as to represent the unit pixel of the image data. Length. For example, it is assumed that the image data captured by the image capturing device 1 includes the river channel. The actual length of the channel in the image data is known to be 500 cm. If the processor 2 reads the image data, the image is calculated. When the channel in the data contains 100 pixels, it can be estimated that the unit pixel represents 5 cm, that is, the ratio of the image is 5 (cm/pixel).
該傾斜角度計算步驟S2,係由該處理器2依該單位像素所代表之實際長度,計算該影像資料中之該水面線A、該基準線B及該垂直高度H的實際長度,並將該基準線B及該垂直高度H的實際長度代入一角度計算方程式以計算該傾斜角度θ 。The tilt angle calculation step S2 is performed by the processor 2, according to the actual length represented by the unit pixel, the actual length of the water surface line A, the reference line B, and the vertical height H in the image data, and the The actual length of the reference line B and the vertical height H is substituted into an angle calculation equation to calculate the inclination angle θ .
更詳言之,由於在該影像比例計算步驟S1中,已得知該影像資料之該單位像素所代表之實際長度,因此,該處理器1便可針對該影像資料,計算該水面線A、該基準線B及該垂直高度H係由多少像素值組成,並進而計算該水面線A、該基準線B及該垂直高度H的實際長度,再接著計算該水面線A及基準線B間的傾斜角度θ
。該角度計算方程式表示如下:
其中,tan-1 (‧)為一反正切函數,h為該垂直高度H的實際長度,b為該基準線B的實際長度。Where tan -1 (‧) is an arctangent function, h is the actual length of the vertical height H, and b is the actual length of the reference line B.
藉此,便可計算出該傾斜角度θ 。Thereby, the inclination angle θ can be calculated.
該水流時間計算步驟S3,係由該處理器1將該傾斜角度θ 及該垂直高度H的實際長度代入一時間計算方程式,以計算一水流時間。The water flow time calculation step S3 is performed by the processor 1 by substituting the inclination angle θ and the actual length of the vertical height H into a time calculation equation to calculate a water flow time.
該時間計算方程式如下所示:
其中,t為該水流時間,sin(‧)為一正弦函數,g為重力加速度值,一般係為9.8(m/s2 )。Where t is the water flow time, sin(‧) is a sinusoidal function, and g is the gravitational acceleration value, which is generally 9.8 (m/s 2 ).
更詳言之,該水流時間係為該河道表面之單位水量行經該垂直高度H的實際長度所需之時間,亦可代表該河道表面之單位水量沿該水面線A之上游端流至該下游端所需的時間。More specifically, the water flow time is the time required for the unit water amount of the river surface to pass the actual length of the vertical height H, or the unit water amount of the river surface may flow along the upstream end of the water surface line A to the downstream. The time required at the end.
該表面流速計算步驟S4,係由該處理器2將該水面線A的實際長度除以該水流時間,以得到一表面流速。The surface flow rate calculation step S4 is performed by the processor 2 by dividing the actual length of the water surface line A by the water flow time to obtain a surface flow rate.
該表面流速的計算如下所示:
其中,a為該水面線A的實際長度。Where a is the actual length of the surface line A.
藉此,便可計算出該表面流速。Thereby, the surface flow rate can be calculated.
本發明河水表面之流速估算方法,係利用一影像擷取裝置1擷取該河道之影像,以產生該影像資料,當該影像資料傳至該處理器2後,並可藉由上述步驟計算得該河水表面之流速。The method for estimating the flow rate of the surface of the river of the present invention uses an image capturing device 1 to capture the image of the channel to generate the image data. When the image data is transmitted to the processor 2, the method can be calculated by the above steps. The flow rate of the river surface.
藉此,本發明係依據該河道之實際長度與傾斜角度,藉由該時間計算方程式計算該水流時間,並進而計算河水表面之流速,具有精確計算流速之功效。Thereby, the invention calculates the water flow time by calculating the equation according to the actual length and the inclination angle of the river channel, and further calculates the flow velocity of the river water surface, and has the effect of accurately calculating the flow velocity.
本發明係藉由該河道的影像資料進行流速計算,而該影像資料係包含該河道之一區域,故所計算出之該河水表面流速係為整體河道區域之平均流速,相較於將數個感應元件設置於河道中以計算整體河道區域 之平均流速,本發明具有降低硬體設置成本之功效。The invention calculates the flow velocity by using the image data of the river channel, and the image data system includes a region of the river channel, so the calculated surface velocity of the river water is the average flow velocity of the entire river channel region, compared to several Inductive components are placed in the river to calculate the overall channel area The average flow rate of the present invention has the effect of reducing the cost of setting the hardware.
雖然本發明已利用上述較佳實施例揭示,然其並非用以限定本發明,任何熟習此技藝者在不脫離本發明之精神和範圍之內,相對上述實施例進行各種更動與修改仍屬本發明所保護之技術範疇,因此本發明之保護範圍當視後附之申請專利範圍所界定者為準。While the invention has been described in connection with the preferred embodiments described above, it is not intended to limit the scope of the invention. The technical scope of the invention is protected, and therefore the scope of the invention is defined by the scope of the appended claims.
S1‧‧‧影像比例計算步驟S1‧‧‧image scale calculation step
S2‧‧‧傾斜角度計算步驟S2‧‧‧ Tilt angle calculation step
S3‧‧‧水流時間計算步驟S3‧‧‧Water flow time calculation steps
S4‧‧‧表面流速計算步驟S4‧‧‧ Surface flow rate calculation steps
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW102101646A TWI465721B (en) | 2013-01-16 | 2013-01-16 | A method for surface velocity estimation of river water |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW102101646A TWI465721B (en) | 2013-01-16 | 2013-01-16 | A method for surface velocity estimation of river water |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201430348A TW201430348A (en) | 2014-08-01 |
TWI465721B true TWI465721B (en) | 2014-12-21 |
Family
ID=51796846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW102101646A TWI465721B (en) | 2013-01-16 | 2013-01-16 | A method for surface velocity estimation of river water |
Country Status (1)
Country | Link |
---|---|
TW (1) | TWI465721B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4888694A (en) * | 1987-10-28 | 1989-12-19 | Quantum Medical Systems, Inc. | Ultrasound imaging system for relatively low-velocity blood flow at relatively high frame rates |
CN1281149A (en) * | 1998-06-10 | 2001-01-24 | 昌民技术株式会社 | Method of measuring overage flow rate of river water |
TWI230783B (en) * | 2001-05-23 | 2005-04-11 | Ge-Fei Liou | Method and system of average flowing velocity measurement in river |
CN202502107U (en) * | 2012-03-20 | 2012-10-24 | 河海大学 | River water surface flow rate monitoring device based on feature tracking |
US8340927B2 (en) * | 2010-04-30 | 2012-12-25 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Apparatus for measuring flow velocity and flow rate of groundwater leaking from earth surface and apparatus for monitoring the same |
-
2013
- 2013-01-16 TW TW102101646A patent/TWI465721B/en not_active IP Right Cessation
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4888694A (en) * | 1987-10-28 | 1989-12-19 | Quantum Medical Systems, Inc. | Ultrasound imaging system for relatively low-velocity blood flow at relatively high frame rates |
CN1281149A (en) * | 1998-06-10 | 2001-01-24 | 昌民技术株式会社 | Method of measuring overage flow rate of river water |
TWI230783B (en) * | 2001-05-23 | 2005-04-11 | Ge-Fei Liou | Method and system of average flowing velocity measurement in river |
US8340927B2 (en) * | 2010-04-30 | 2012-12-25 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Apparatus for measuring flow velocity and flow rate of groundwater leaking from earth surface and apparatus for monitoring the same |
CN202502107U (en) * | 2012-03-20 | 2012-10-24 | 河海大学 | River water surface flow rate monitoring device based on feature tracking |
Also Published As
Publication number | Publication date |
---|---|
TW201430348A (en) | 2014-08-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI448977B (en) | Method and apparatus for video analytics based object counting | |
JP6356266B2 (en) | Crowd monitoring system | |
WO2013132951A1 (en) | Distance calculation device and distance calculation method | |
JP6138861B2 (en) | Distance calculation device | |
JP6620881B2 (en) | Apparatus, system, method, and program for measuring the number of passengers, and apparatus, method, and program for calculating travel distance of vehicle | |
US11170229B2 (en) | Attitude estimation device | |
JP6217635B2 (en) | Fall detection device, fall detection method, fall detection camera, and computer program | |
JP2008116373A5 (en) | ||
JP2012103239A5 (en) | ||
CN107492123B (en) | Road monitoring camera self-calibration method using road surface information | |
KR101921610B1 (en) | Method and Apparatus for Monitoring Objects from Video | |
US11017552B2 (en) | Measurement method and apparatus | |
JP2011134012A5 (en) | ||
JP2014169939A5 (en) | ||
CN108509963B (en) | Target difference detection method and target difference detection equipment based on deep learning | |
JP2016217944A5 (en) | ||
JP2006318064A (en) | Image sensor | |
TWI465721B (en) | A method for surface velocity estimation of river water | |
JP2015534370A5 (en) | ||
JP7372616B2 (en) | Stone and gravel detection system, stone and gravel detection method and program | |
CN106796726A (en) | Method for testing motion and system | |
JPWO2021070415A5 (en) | ||
JP2015220623A (en) | Mobile body imaging system | |
JP5990042B2 (en) | Passing number estimation method and computer program | |
JP2011053378A5 (en) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Annulment or lapse of patent due to non-payment of fees |