CN104408750A - Method and system for calculating proportion of seawater to land - Google Patents

Method and system for calculating proportion of seawater to land Download PDF

Info

Publication number
CN104408750A
CN104408750A CN201410605763.9A CN201410605763A CN104408750A CN 104408750 A CN104408750 A CN 104408750A CN 201410605763 A CN201410605763 A CN 201410605763A CN 104408750 A CN104408750 A CN 104408750A
Authority
CN
China
Prior art keywords
land
point
seawater
calculate
grayscale image
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201410605763.9A
Other languages
Chinese (zh)
Other versions
CN104408750B (en
Inventor
李晴岚
汤婷婷
彭玉龙
陈城
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201410605763.9A priority Critical patent/CN104408750B/en
Publication of CN104408750A publication Critical patent/CN104408750A/en
Application granted granted Critical
Publication of CN104408750B publication Critical patent/CN104408750B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/62Analysis of geometric attributes of area, perimeter, diameter or volume
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10132Ultrasound image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30181Earth observation
    • G06T2207/30192Weather; Meteorology

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种计算海水陆地比例的方法,包括:根据陆地、海水的下垫面情况,绘制并剪裁得到欲计算海水陆地比例地区的二值灰度图像;检测得到上述二值灰度图像中各边界点的像素;在上述检测得到的像素范围内,计算得到任意点到位置点的距离;利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。本发明还涉及一种计算海水陆地比例的系统。本发明可以精确地计算出各个经纬度位置点其预设范围内的海水与陆地的面积比例,改善现阶段近海台风强度预报难度大的现状。

The invention relates to a method for calculating the ratio of seawater to land, which includes: drawing and cutting out a binary grayscale image of the area where the ratio of seawater to land is to be calculated according to the underlying surface of the land and seawater; Pixels of each boundary point; within the pixel range obtained by the above detection, calculate the distance from any point to the location point; use the distance from any point to the location point obtained above to calculate the area ratio of seawater to land within the Rad range. The invention also relates to a system for calculating the ratio of seawater to land. The present invention can accurately calculate the area ratio of seawater and land within the preset range of each longitude and latitude position point, and improves the current situation that the forecasting of offshore typhoon intensity is difficult at the present stage.

Description

一种计算海水陆地比例的方法及系统A method and system for calculating the ratio of seawater to land

技术领域 technical field

本发明涉及一种计算海水陆地比例的方法及系统。  The invention relates to a method and system for calculating the ratio of seawater to land. the

背景技术 Background technique

台风,亦称飓风,是形成于热带或副热带海面温度在26℃以上的广阔海面上的热带气旋。台风过境时常常带来狂风暴雨天气,引起海面巨浪,严重威胁航海安全。台风登陆后带来的狂风暴雨可摧毁庄稼、各种建筑设施等,造成人民生命、财产的巨大损失。因此,台风登录前及时准确地向公众预报台风路径、台风强度、降雨和大风的强度及分布情况,对保护生命安全、降低财产损失是至关重要的。  A typhoon, also known as a hurricane, is a tropical cyclone formed on a vast ocean with a tropical or subtropical sea surface temperature above 26°C. When a typhoon crosses the border, it often brings violent stormy weather, causing huge waves on the sea surface, and seriously threatening the safety of navigation. The violent wind and rain brought by the typhoon after landfall can destroy crops, various building facilities, etc., causing huge losses of people's lives and property. Therefore, timely and accurate forecasting of typhoon track, typhoon intensity, intensity and distribution of rainfall and strong winds to the public before typhoon lands is crucial to protecting life safety and reducing property losses. the

随着数值模式的发展,在台风的预报研究中,台风路径预报比台风强度和台风风雨分布的预报遥遥领先。在台风强度预报方面,过去20多年来,与世界各大台风预报中心一样,我国台风强度业务预报进展非常缓慢,目前业务预报中广泛应用的方法主要是一些气候持续性方法和统计动力模式,国外现在越来越多用统计方法去研究台风强度变化,如美国联合台风警报中心的台风强度统计预报(STIPS)、美国国家飓风中心的飓风强度统计预报模式(SHIFOR)和飓风强度统计预报方案(SHIPS)等,目前我国尚无类似的探讨研究。  With the development of numerical models, in typhoon forecast research, typhoon track forecast is far ahead of typhoon intensity and typhoon rain distribution forecast. In terms of typhoon intensity forecasting, in the past 20 years, like the major typhoon forecasting centers in the world, my country’s typhoon intensity operational forecasting has progressed very slowly. At present, the methods widely used in operational forecasting are mainly some climate persistence methods and statistical dynamic models. Now more and more statistical methods are used to study the change of typhoon intensity, such as the Statistical Typhoon Intensity Forecast (STIPS) of the United States Joint Typhoon Warning Center, the Statistical Hurricane Intensity Forecast Model (SHIFOR) and the Statistical Hurricane Intensity Forecasting Program (SHIPS) of the National Hurricane Center At present, there is no similar research in our country. the

近海台风强度演变的影响因素有:垂直风切变、海表温度、中高层环流特征、水汽输送等,其中,最重要的是下垫面的影响,因为海水是台风动力的来源。那么如何计算下垫面水陆比,目前没有相应的方法可供参考。  Factors affecting the evolution of offshore typhoon intensity include: vertical wind shear, sea surface temperature, mid-level and upper-level circulation characteristics, water vapor transport, etc. Among them, the most important is the influence of the underlying surface, because sea water is the source of typhoon power. So how to calculate the water-to-land ratio of the underlying surface, there is currently no corresponding method for reference. the

发明内容 Contents of the invention

有鉴于此,有必要提供一种计算海水陆地比例的方法及系统。  In view of this, it is necessary to provide a method and system for calculating the ratio of seawater to land. the

本发明提供一种计算海水陆地比例的方法,该方法包括如下步骤:a.根据陆地、海水的下垫面情况,绘制并剪裁得到欲计算海水陆地比例地区的二值灰度图像;b.检测得到上述二值灰度图像中各边界点的像素;c.在上述检测得到的像素范围内,计算得到任意点到位置点的距离;d.利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。  The present invention provides a method for calculating the ratio of seawater to land, the method comprising the following steps: a. drawing and tailoring the binary grayscale image of the area where the ratio of seawater to land is to be calculated according to the underlying surface of the land and seawater; b. detecting Obtain the pixels of each boundary point in the above-mentioned binary grayscale image; c. within the pixel range obtained by the above-mentioned detection, calculate the distance from any point to the position point; d. use the distance from any point to the position point obtained above to calculate Get the ratio of seawater to land area in the Rad range. the

其中,所述的步骤a包括:根据陆地、海水下垫面情况,绘制高分辨率的陆地、海水二维下垫面图像;从上述绘制的二维下垫面图像中,剪裁得到欲计算海水陆地比例地区的下垫面图像;将上述剪裁好的下垫面图像处理得到二值灰度图像;读取上述处理得到的二值灰度图像的二维矩阵形式。  Wherein, the step a includes: drawing a high-resolution two-dimensional image of the underlying surface of land and seawater according to the conditions of the underlying surface of the land and seawater; An image of the underlying surface of a land-scale area; process the clipped image of the underlying surface to obtain a binary grayscale image; read the two-dimensional matrix form of the binary grayscale image obtained through the above processing. the

所述的二值灰度图像通过PHOTOSHOP等绘图软件处理得到。  The binary grayscale image is obtained by processing drawing software such as PHOTOSHOP. the

所述的步骤c中计算A点到B点的距离为:  In the described step c, the distance from point A to point B is calculated as:

C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )

,  ,

dist=R*arcos(C)  dist=R*arcos(C)

其中,A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R。  Among them, the longitude of point A is α1, and the latitude is β1; the longitude of point B is α2, and the latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R. the

所述的步骤d包括:计算得出所有距离小于半径Rad范围内的像素点的个数总和S;计算出Rad范围内的二维矩阵值累加之和T;计算T与S之比值,得到以Rad为半径范围内的海水与陆地的面积比值。  The step d includes: calculating the sum S of the number of all pixels whose distance is less than the radius Rad; calculating the cumulative sum T of two-dimensional matrix values within the Rad range; calculating the ratio of T to S to obtain Rad is the ratio of seawater to land area within the radius. the

本发明还提供一种计算海水陆地比例的系统,包括绘制模块、检测模块及计算模块,其中:所述绘制模块用于根据陆地、海水的下垫面情 况,绘制并剪裁得到欲计算海水陆地比例地区的二值灰度图像;所述检测模块用于检测得到上述二值灰度图像中各边界点的像素;所述计算模块用于在上述检测得到的像素范围内,计算得到任意点到位置点的距离;所述计算模块还用于利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。  The present invention also provides a system for calculating the proportion of seawater and land, including a drawing module, a detection module and a calculation module, wherein: the drawing module is used to draw and cut out the seawater and land to be calculated according to the conditions of the underlying surface of the land and seawater The binary grayscale image of the proportional area; the detection module is used to detect the pixels of each boundary point in the above binary grayscale image; the calculation module is used to calculate any point to The distance of the location point; the calculation module is also used to calculate the area ratio of seawater to land within the range of Rad by using the distance from any point to the location point obtained above. the

其中,所述绘制模块具体用于:根据陆地、海水下垫面情况,绘制高分辨率的陆地、海水二维下垫面图像;从上述绘制的二维下垫面图像中,剪裁得到欲计算海水陆地比例地区的下垫面图像;将上述剪裁好的下垫面图像处理得到二值灰度图像;读取上述处理得到的二值灰度图像的二维矩阵形式。  Wherein, the drawing module is specifically used for: drawing a high-resolution two-dimensional underlying surface image of land and sea water according to the conditions of the underlying surface of land and sea; The image of the underlying surface of the area with a proportion of seawater and land; process the clipped image of the underlying surface to obtain a binary grayscale image; read the two-dimensional matrix form of the binary grayscale image obtained through the above processing. the

所述的二值灰度图像通过PHOTOSHOP等绘图软件处理得到。  The binary grayscale image is obtained by processing drawing software such as PHOTOSHOP. the

所述计算模块计算A点到B点的距离为:  The calculation module calculates the distance from point A to point B as:

C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )

,  ,

dist=R*arcos(C)  dist=R*arcos(C)

其中,A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R。  Among them, the longitude of point A is α1, and the latitude is β1; the longitude of point B is α2, and the latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R. the

所述计算模块具体用于:计算得出所有距离小于半径Rad范围内的像素点的个数总和S;计算出Rad范围内的二维矩阵值累加之和T;计算T与S之比值,得到以Rad为半径范围内的海水与陆地的面积比值。  The calculation module is specifically used for: calculating the sum S of the number of pixels within the range of all distances less than the radius Rad; calculating the cumulative sum T of two-dimensional matrix values within the range of Rad; calculating the ratio of T to S to obtain Ratio of seawater to land area within a radius of Rad. the

本发明一种计算海水陆地比例的方法及系统,可以精确地计算出各个经纬度位置点其预设范围内的海水与陆地的面积比例,分析和预测台风近海登陆时海水面积以及海水温度对台风强度的影响,改善现阶段近海台风强度预报难度大的现状。  The present invention is a method and system for calculating the proportion of seawater and land, which can accurately calculate the area ratio of seawater and land within the preset range of each longitude and latitude position point, and analyze and predict the impact of seawater area and seawater temperature on typhoon intensity when a typhoon lands offshore. To improve the current situation that it is very difficult to forecast the intensity of offshore typhoons. the

附图说明 Description of drawings

图1为本发明一种计算海水陆地比例的方法的流程图;  Fig. 1 is a kind of flow chart of the method for calculating seawater land ratio of the present invention;

图2为本发明一种计算海水陆地比例的系统的硬件架构图。  FIG. 2 is a hardware architecture diagram of a system for calculating the ratio of seawater to land according to the present invention. the

具体实施方式 Detailed ways

下面结合附图及具体实施例对本发明作进一步详细的说明。  The present invention will be described in further detail below in conjunction with the accompanying drawings and specific embodiments. the

参阅图1所示,是本发明一种计算海水陆地比例的方法较佳实施例的作业流程图。  Referring to Fig. 1, it is a flow chart of the operation of a preferred embodiment of the method for calculating the proportion of seawater and land in the present invention. the

本实施例以华南地区下垫面图像为例进行说明:  In this embodiment, the image of the underlying surface in South China is taken as an example for illustration:

步骤S401,根据北半球陆地、海水的下垫面情况,绘制并剪裁得到华南地区的二值灰度图像。具体如下:  Step S401, according to the conditions of the underlying surface of land and sea in the northern hemisphere, draw and cut out a binary grayscale image of South China. details as follows:

第一步,根据北半球东南亚陆地、太平洋下垫面情况,绘制高分辨率的北半球东南亚陆地、太平洋二维下垫面图像。  The first step is to draw a high-resolution two-dimensional image of the land of Southeast Asia in the northern hemisphere and the underlying surface of the Pacific Ocean according to the conditions of the underlying surface of Southeast Asia in the northern hemisphere and the Pacific Ocean. the

第二步,从上述绘制的二维下垫面图像中,剪裁得到华南地区的下垫面图像。因为华南地区的纬度范围为1~52度、经度范围为98~153度,故本实施例剪裁出纬度范围1~52度、经度范围98~153度的下垫面图像。  In the second step, the image of the underlying surface in South China is obtained by cropping the two-dimensional underlying surface image drawn above. Since the latitude range of South China is 1-52 degrees and the longitude range is 98-153 degrees, this embodiment cuts out the underlying surface image with the latitude range of 1-52 degrees and the longitude range of 98-153 degrees. the

第三步,将上述剪裁好的下垫面图像处理得到二值灰度图像。具体而言:本实施例将上述剪裁好的下垫面图像通过绘图软件处理成二值灰度图像。所述绘图软件包括:PHOTOSHOP及Windows系统的绘图软件。  The third step is to process the cropped image of the underlying surface to obtain a binary grayscale image. Specifically: In this embodiment, the cropped image of the underlying surface is processed into a binary grayscale image by drawing software. Described drawing software comprises: the drawing software of PHOTOSHOP and Windows system. the

第四步,读取上述处理得到的二值灰度图像。具体而言:本实施例读取上述二值灰度图像的结果为二维矩阵形式,其中“1”代表陆地,“0”代表海洋。  The fourth step is to read the binary grayscale image obtained by the above processing. Specifically: the result of reading the binary grayscale image in this embodiment is in the form of a two-dimensional matrix, where "1" represents land and "0" represents ocean. the

步骤S402,检测得到上述二值灰度图像中各边界点的像素,使其像素点范围与经纬度范围一一对应。具体而言:  Step S402, detecting and obtaining the pixels of each boundary point in the above-mentioned binary grayscale image, so that the pixel point range corresponds to the latitude and longitude range one by one. in particular:

本实施例使用Windows系统的绘图软件,检测得到上述二值灰度图像中各边界点的像素,通过程序代码实现遍历该像素点数值的范围,使 其与经纬度范围一一对应。  The present embodiment uses the drawing software of Windows system, detects and obtains the pixel of each boundary point in the above-mentioned binary gray scale image, realizes traversing the range of this pixel value by program code, makes it one-to-one correspondence with latitude and longitude range. the

步骤S403,在上述检测得到的像素范围内,计算得到任意点到位置点的距离。具体而言:  Step S403, within the pixel range obtained by the above detection, calculate the distance from any point to the position point. in particular:

本实施例利用下述方法得到球面上任意两点的距离:  This embodiment uses the following method to obtain the distance between any two points on the sphere:

假设A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R(Km)。同时约定:东经为正,西经为负;北纬为正,南纬为负。则A点到B点的距离为:  Assume that the longitude of point A is α1 and latitude is β1; the longitude of point B is α2 and latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R (Km). At the same time, it is agreed that east longitude is positive and west longitude is negative; north latitude is positive and south latitude is negative. Then the distance from point A to point B is:

C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )

公式(1)  Formula 1)

dist=R*arcos(C)  dist=R*arcos(C)

步骤S404,利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。具体如下:  Step S404, using the above obtained distance from any point to the location point, to calculate the area ratio of seawater to land within the range of Rad. details as follows:

本实施例通过比较可获得由公式(1)计算得出所有dist小于半径Rad范围内的像素点的个数总和S,再由步骤S401得出二维矩阵,计算出Rad范围内的二维矩阵值累加之和T,T与S之比值即可计算得到该位置点以Rad为半径范围内的海水与陆地的面积比值。  In this embodiment, the sum S of the number of pixels whose dist is less than the radius Rad can be calculated by formula (1) by comparison, and then the two-dimensional matrix is obtained by step S401, and the two-dimensional matrix in the Rad range is calculated. The sum T of value accumulation, the ratio of T and S can be calculated to obtain the area ratio of seawater and land within the radius of Rad at the position point. the

参阅图2所示,是本发明一种计算海水陆地比例的系统的硬件架构图。该系统包括绘制模块、检测模块及计算模块。  Referring to FIG. 2 , it is a hardware architecture diagram of a system for calculating the ratio of seawater to land according to the present invention. The system includes a drawing module, a detection module and a calculation module. the

本实施例以华南地区下垫面图像为例进行说明:  In this embodiment, the image of the underlying surface in South China is taken as an example for illustration:

所述绘制模块用于根据北半球陆地、海水的下垫面情况,绘制并剪裁得到华南地区的二值灰度图像。具体如下:  The drawing module is used to draw and cut out the binary grayscale image of South China according to the underlying surface of land and sea in the northern hemisphere. details as follows:

第一步,根据北半球东南亚陆地、太平洋下垫面情况,绘制高分辨率的北半球东南亚陆地、太平洋二维下垫面图像。  The first step is to draw a high-resolution two-dimensional image of the land of Southeast Asia in the northern hemisphere and the underlying surface of the Pacific Ocean according to the conditions of the underlying surface of Southeast Asia in the northern hemisphere and the Pacific Ocean. the

第二步,从上述绘制的二维下垫面图像中,剪裁得到华南地区的下垫面图像。因为华南地区的纬度范围为1~52度、经度范围为98~153度,故本实施例剪裁出纬度范围1~52度、经度范围98~153度的下垫面图像。  In the second step, the image of the underlying surface in South China is obtained by cropping the two-dimensional underlying surface image drawn above. Since the latitude range of South China is 1-52 degrees and the longitude range is 98-153 degrees, this embodiment cuts out the underlying surface image with the latitude range of 1-52 degrees and the longitude range of 98-153 degrees. the

第三步,将上述剪裁好的下垫面图像通过绘图软件处理得到二值灰度图像。具体而言:本实施例将上述剪裁好的下垫面图像通过绘图软件处理成二值灰度图像。所述绘图软件包括:PHOTOSHOP及Windows系统的绘图软件。  The third step is to process the clipped image of the underlying surface through drawing software to obtain a binary grayscale image. Specifically: In this embodiment, the cropped image of the underlying surface is processed into a binary grayscale image by drawing software. Described drawing software comprises: the drawing software of PHOTOSHOP and Windows system. the

第四步,读取上述处理得到的二值灰度图像。具体而言:本实施例读取上述二值灰度图像的结果为二维矩阵形式,其中“1”代表陆地,“0”代表海洋。  The fourth step is to read the binary grayscale image obtained by the above processing. Specifically: the result of reading the binary grayscale image in this embodiment is in the form of a two-dimensional matrix, where "1" represents land and "0" represents ocean. the

所述检测模块用于检测得到上述二值灰度图像中各边界点的像素,使其像素点范围与经纬度范围一一对应。具体而言:  The detection module is used to detect the pixels of each boundary point in the above-mentioned binary grayscale image, so that the pixel point range corresponds to the latitude and longitude range one by one. in particular:

本实施例使用Windows系统的绘图软件,检测得到上述二值灰度图像中各边界点的像素,通过程序代码实现遍历该像素点数值的范围,使其与经纬度范围一一对应。  This embodiment uses the drawing software of the Windows system to detect the pixels of each boundary point in the above-mentioned binary grayscale image, and realizes traversing the range of the pixel value through the program code, so that it corresponds to the latitude and longitude range one by one. the

所述计算模块用于在上述检测得到的像素范围内,计算得到任意点到位置点的距离。具体而言:  The calculation module is used to calculate the distance from any point to the position point within the detected pixel range. in particular:

本实施例利用下述方法得到球面上任意两点的距离:  This embodiment uses the following method to obtain the distance between any two points on the sphere:

假设A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R(Km)。同时约定:东经为正,西经为负;北纬为正,南纬为负。则A点到B点的距离为:  Assume that the longitude of point A is α1 and latitude is β1; the longitude of point B is α2 and latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R (Km). At the same time, it is agreed that east longitude is positive and west longitude is negative; north latitude is positive and south latitude is negative. Then the distance from point A to point B is:

C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )

公式(1)  Formula 1)

dist=R*arcos(C)  dist=R*arcos(C)

所述计算模块还用于利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。具体如下:  The calculation module is also used to calculate the area ratio of seawater to land within the range of Rad by using the distance from any point obtained above to the location point. details as follows:

本实施例通过比较可获得由公式(1)计算得出所有dist小于半径Rad范围内的像素点的个数总和S,再由步骤S401得出二维矩阵,计算出Rad范围内的二维矩阵值累加之和T,T与S之比值即可计算得到该 位置点以Rad为半径范围内的海水与陆地的面积比值。  In this embodiment, the sum S of the number of pixels whose dist is less than the radius Rad can be calculated by formula (1) by comparison, and then the two-dimensional matrix is obtained by step S401, and the two-dimensional matrix in the Rad range is calculated. The sum T of value accumulation, the ratio of T and S can be calculated to obtain the area ratio of seawater and land within the radius range of Rad at the position point. the

本发明已通过各种经纬度位置点的有效测试,经实际情况的对比发现,本发明可以精确地计算出各个位置的海水与陆地的面积比例。  The present invention has been effectively tested at various latitude and longitude locations, and found through comparison of actual conditions that the present invention can accurately calculate the area ratio of seawater and land at each location. the

虽然本发明参照当前的较佳实施方式进行了描述,但本领域的技术人员应能理解,上述较佳实施方式仅用来说明本发明,并非用来限定本发明的保护范围,任何在本发明的精神和原则范围之内,所做的任何修饰、等效替换、改进等,均应包含在本发明的权利保护范围之内。  Although the present invention has been described with reference to the current preferred embodiments, those skilled in the art should understand that the above-mentioned preferred embodiments are only used to illustrate the present invention, and are not used to limit the protection scope of the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and scope of principles shall be included in the protection scope of the present invention. the

Claims (10)

1.一种计算海水陆地比例的方法,其特征在于,该方法包括如下步骤:1. A method for calculating seawater and land ratios, characterized in that the method may further comprise the steps: a.根据陆地、海水的下垫面情况,绘制并剪裁得到欲计算海水陆地比例地区的二值灰度图像;a. According to the underlying surface of land and sea water, draw and cut out the binary grayscale image of the area where the ratio of sea water to land is to be calculated; b.检测得到上述二值灰度图像中各边界点的像素;b. Detecting the pixels of each boundary point in the binary grayscale image; c.在上述检测得到的像素范围内,计算得到任意点到位置点的距离;c. Within the range of pixels detected above, calculate the distance from any point to the position point; d.利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。d. Using the distance from any point obtained above to the location point, calculate the area ratio of seawater to land within the Rad range. 2.如权利要求1所述的方法,其特征在于,所述的步骤a包括:2. method as claimed in claim 1, is characterized in that, described step a comprises: 根据陆地、海水下垫面情况,绘制高分辨率的陆地、海水二维下垫面图像;Draw high-resolution two-dimensional images of land and seawater underlying surfaces according to the conditions of land and seawater underlying surfaces; 从上述绘制的二维下垫面图像中,剪裁得到欲计算海水陆地比例地区的下垫面图像;From the two-dimensional underlying surface image drawn above, crop the underlying surface image of the area where the proportion of seawater and land is to be calculated; 将上述剪裁好的下垫面图像处理得到二值灰度图像;Processing the image of the underlying surface cut out above to obtain a binary grayscale image; 读取上述处理得到的二值灰度图像的二维矩阵形式。Read the two-dimensional matrix form of the binary grayscale image obtained from the above processing. 3.如权利要求2所述的方法,其特征在于,所述的二值灰度图像通过PHOTOSHOP等绘图软件处理得到。3. The method according to claim 2, wherein the binary grayscale image is obtained by processing drawing software such as PHOTOSHOP. 4.如权利要求1所述的方法,其特征在于,所述的步骤c中计算A点到B点的距离为:4. the method for claim 1 is characterized in that, in the described step c, calculate the distance from point A to point B as: C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12)C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )                                                       ,                          dist=R*arcos(C)dist=R*arcos(C) 其中,A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R。Among them, the longitude of point A is α1, and the latitude is β1; the longitude of point B is α2, and the latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R. 5.如权利要求1所述的方法,其特征在于,所述的步骤d包括:5. method as claimed in claim 1, is characterized in that, described step d comprises: 计算得出所有距离小于半径Rad范围内的像素点的个数总和S;Calculate the sum S of the number of all pixels whose distance is less than the radius Rad; 计算出Rad范围内的二维矩阵值累加之和T;Calculate the cumulative sum T of the two-dimensional matrix values within the Rad range; 计算T与S之比值,得到以Rad为半径范围内的海水与陆地的面积比值。Calculate the ratio of T to S to obtain the ratio of the area of seawater to land within the radius of Rad. 6.一种计算海水陆地比例的系统,其特征在于,该系统包括绘制模块、检测模块及计算模块,其中:6. A system for calculating the proportion of seawater and land, characterized in that the system includes a drawing module, a detection module and a calculation module, wherein: 所述绘制模块用于根据陆地、海水的下垫面情况,绘制并剪裁得到欲计算海水陆地比例地区的二值灰度图像;The drawing module is used to draw and cut out the binary grayscale image of the area where the proportion of seawater and land is to be calculated according to the underlying surface of land and seawater; 所述检测模块用于检测得到上述二值灰度图像中各边界点的像素;The detection module is used to detect the pixels of each boundary point in the binary grayscale image; 所述计算模块用于在上述检测得到的像素范围内,计算得到任意点到位置点的距离;The calculation module is used to calculate the distance from any point to the position point within the range of pixels detected above; 所述计算模块还用于利用上述得到的任意点到位置点的距离,计算得到Rad范围内海水与陆地的面积比值。The calculation module is also used to calculate the area ratio of seawater to land within the range of Rad by using the distance from any point obtained above to the location point. 7.如权利要求6所述的系统,其特征在于,所述绘制模块具体用于:7. The system according to claim 6, wherein the drawing module is specifically used for: 根据陆地、海水下垫面情况,绘制高分辨率的陆地、海水二维下垫面图像;Draw high-resolution two-dimensional images of land and seawater underlying surfaces according to the conditions of land and seawater underlying surfaces; 从上述绘制的二维下垫面图像中,剪裁得到欲计算海水陆地比例地区的下垫面图像;From the two-dimensional underlying surface image drawn above, crop the underlying surface image of the area where the proportion of seawater and land is to be calculated; 将上述剪裁好的下垫面图像处理得到二值灰度图像;Processing the image of the underlying surface cut out above to obtain a binary grayscale image; 读取上述处理得到的二值灰度图像的二维矩阵形式。Read the two-dimensional matrix form of the binary grayscale image obtained from the above processing. 8.如权利要求7所述的系统,其特征在于,所述的二值灰度图像通过PHOTOSHOP等绘图软件处理得到。8. The system according to claim 7, wherein the binary grayscale image is obtained by processing drawing software such as PHOTOSHOP. 9.如权利要求6所述的系统,其特征在于,所述计算模块计算A点到B点的距离为:9. The system according to claim 6, wherein the calculation module calculates the distance from point A to point B as: C=sin(β1)*sin(β2)+cos(β1)*cos(β2)*cos(α12)C=sin(β 1 )*sin(β 2 )+cos(β 1 )*cos(β 2 )*cos(α 12 )                                              ,                      dist=R*arcos(C)dist=R*arcos(C) 其中,A点的经度是α1、纬度是β1;B点的经度是α2、纬度是β2,其中,α1,α2,β1,β2均为弧度,地球半径为R。Among them, the longitude of point A is α1, and the latitude is β1; the longitude of point B is α2, and the latitude is β2, where α1, α2, β1, and β2 are all radians, and the radius of the earth is R. 10.如权利要求6所述的系统,其特征在于,所述计算模块具体用于:10. The system according to claim 6, wherein the calculation module is specifically used for: 计算得出所有距离小于半径Rad范围内的像素点的个数总和S;Calculate the sum S of the number of all pixels whose distance is less than the radius Rad; 计算出Rad范围内的二维矩阵值累加之和T;Calculate the cumulative sum T of the two-dimensional matrix values within the Rad range; 计算T与S之比值,得到以Rad为半径范围内的海水与陆地的面积比值。Calculate the ratio of T to S to obtain the ratio of the area of seawater to land within the radius of Rad.
CN201410605763.9A 2014-10-30 2014-10-30 A kind of method and system for calculating seawater land ratio Active CN104408750B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410605763.9A CN104408750B (en) 2014-10-30 2014-10-30 A kind of method and system for calculating seawater land ratio

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410605763.9A CN104408750B (en) 2014-10-30 2014-10-30 A kind of method and system for calculating seawater land ratio

Publications (2)

Publication Number Publication Date
CN104408750A true CN104408750A (en) 2015-03-11
CN104408750B CN104408750B (en) 2017-12-05

Family

ID=52646378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410605763.9A Active CN104408750B (en) 2014-10-30 2014-10-30 A kind of method and system for calculating seawater land ratio

Country Status (1)

Country Link
CN (1) CN104408750B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388536A (en) * 2015-11-10 2016-03-09 中国科学院深圳先进技术研究院 Forecasting method and system of speed of instantaneous extremely strong wind caused by tropical cyclone at coastal region
CN107679760A (en) * 2017-10-16 2018-02-09 山东浪潮通软信息科技有限公司 A kind of monitoring and managing method and device
CN114942481A (en) * 2022-06-08 2022-08-26 中国气象局上海台风研究所(上海市气象科学研究所) Method and device for forecasting extreme value wind speed probability in tropical cyclone process and computer equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102542259A (en) * 2011-12-29 2012-07-04 宁波工程学院 Identification method for near-shore on-land water body
CN102708369A (en) * 2012-05-14 2012-10-03 大连理工大学 Sea ice parameter extraction method on basis of satellite image
CN104077777A (en) * 2014-07-04 2014-10-01 中国科学院大学 Sea surface vessel target detection method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102542259A (en) * 2011-12-29 2012-07-04 宁波工程学院 Identification method for near-shore on-land water body
CN102708369A (en) * 2012-05-14 2012-10-03 大连理工大学 Sea ice parameter extraction method on basis of satellite image
CN104077777A (en) * 2014-07-04 2014-10-01 中国科学院大学 Sea surface vessel target detection method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PATRICK J. FITZPATRICK: "Understanding and forecasting tropical cyclone intensity change with the typhoon intensity prediction scheme(TIPS)", 《WEATHERAND FORECASTING》 *
张杰: "卫星图片上估算点经、纬度及区域面积之方法", 《新疆气象》 *
韩忠民: "知经纬度计算两点精确距离", 《科学传播》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105388536A (en) * 2015-11-10 2016-03-09 中国科学院深圳先进技术研究院 Forecasting method and system of speed of instantaneous extremely strong wind caused by tropical cyclone at coastal region
CN107679760A (en) * 2017-10-16 2018-02-09 山东浪潮通软信息科技有限公司 A kind of monitoring and managing method and device
CN114942481A (en) * 2022-06-08 2022-08-26 中国气象局上海台风研究所(上海市气象科学研究所) Method and device for forecasting extreme value wind speed probability in tropical cyclone process and computer equipment

Also Published As

Publication number Publication date
CN104408750B (en) 2017-12-05

Similar Documents

Publication Publication Date Title
Schenkman et al. Impact of CASA radar and Oklahoma Mesonet data assimilation on the analysis and prediction of tornadic mesovortices in an MCS
US10725205B2 (en) Forensic weather system
Plomaritis et al. Variability in storm climate along the Gulf of Cadiz: the role of large scale atmospheric forcing and implications to coastal hazards
Viale et al. Summer precipitation events over the western slope of the subtropical Andes
Gesch Consideration of vertical uncertainty in elevation-based sea-level rise assessments: Mobile Bay, Alabama case study
Needham et al. Correlating storm surge heights with tropical cyclone winds at and before landfall
CN111340668B (en) Typhoon disaster assessment system
Andrade et al. Coastal flooding hazard related to swell events in Cartagena de Indias, Colombia
Li et al. Storm surge prediction: present status and future challenges
CN105653826A (en) Improved maritime search and rescue region predicting method and system
CN106199606A (en) A kind of multi thresholds squall line recognition methods based on radar return 3 d mosaics
CN105388536B (en) Tropical cyclone triggers coastal area instantaneous pole strong wind wind speed forecasting method and system
CN104408750B (en) A kind of method and system for calculating seawater land ratio
US20180156940A1 (en) Forensic weather system
Nurfaida et al. Intensifying swells and their impacts on the south coast of Java, Indonesia
Lee et al. Dual-Doppler radar analysis of a near-shore line-shaped convective system on 27 July 2011, Korea: a case study
Lee et al. A study on a predictive index for dangerous rip current generation at Haeundae Beach
Kim et al. Localizing local vulnerabilities in urban areas using crowdsourced visual data from participatory sensing
US10345485B2 (en) Forensic weather system
Mohapatra Tropical cyclone forecast verification by India Meteorological Department for north Indian Ocean: A review
CN113031012B (en) Storm surge monitoring method for shore-based GNSS reflected signals
CN106447126A (en) Determining and forecasting methods and determining and forecasting systems for visibility change, caused by tropical cyclone, in coastal area
Choi et al. A rip current warning system based on real-time observations for Haeundae Beach, Korea
Shamsad et al. Sea surface temperature anomaly in the Bay of Bengal in 2010
Gomes et al. Nearshore wave analysis in the Brazilian Northeast based on observations and numerical models

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant