CN117422832A - Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite - Google Patents

Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite Download PDF

Info

Publication number
CN117422832A
CN117422832A CN202311463858.7A CN202311463858A CN117422832A CN 117422832 A CN117422832 A CN 117422832A CN 202311463858 A CN202311463858 A CN 202311463858A CN 117422832 A CN117422832 A CN 117422832A
Authority
CN
China
Prior art keywords
typhoon
source
data
satellite data
orbit satellite
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
CN202311463858.7A
Other languages
Chinese (zh)
Other versions
CN117422832B (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.)
Guangzhou Meteorological Satellite Ground Station Guangdong Meteorological Satellite Remote Sensing Center
Original Assignee
Guangzhou Meteorological Satellite Ground Station Guangdong Meteorological Satellite Remote Sensing Center
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 Guangzhou Meteorological Satellite Ground Station Guangdong Meteorological Satellite Remote Sensing Center filed Critical Guangzhou Meteorological Satellite Ground Station Guangdong Meteorological Satellite Remote Sensing Center
Priority to CN202311463858.7A priority Critical patent/CN117422832B/en
Publication of CN117422832A publication Critical patent/CN117422832A/en
Application granted granted Critical
Publication of CN117422832B publication Critical patent/CN117422832B/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
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Image Processing (AREA)

Abstract

The invention discloses an automatic matching and quantitative display method for a multi-source-orbit satellite typhoon warm center three-dimensional structure, which comprises the following steps: s1, searching typhoon path information to obtain the position of a typhoon real-time center point; s2, automatically matching the effective multi-source orbit satellite data of the corresponding time according to the position of the typhoon real-time center point; and S3, drawing a three-dimensional structure slice dynamic diagram and a vertical outline line diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, and quantitatively displaying the three-dimensional structure of the typhoon heating center. According to the invention, by researching the three-dimensional warm-core structure of typhoons, the change characteristics of the three-dimensional warm-core structure of typhoons in different development processes can be visually seen, and a remote sensing scientific monitoring means is provided for typhoons monitoring, forecasting and analyzing.

Description

Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite
Technical Field
The invention belongs to the typhoon observation field, and particularly relates to an automatic matching and quantitative display method for a multi-source orbit satellite typhoon heating core three-dimensional structure.
Background
Typhoons are extremely damaging disastrous weather systems, and strong wind and strong rainfall brought by typhoons have great influence on human life safety and socioeconomic development, and deviation of typhoons structure forecast can directly influence the forecast accuracy of the wind and rain range and strength caused by typhoons. In recent years, a plurality of scientists develop a discussion of the high-rise warm core structure of typhoons by utilizing satellite data, and find that differences and changes of the height, strength, range, shape and the like of the upper-rise warm region in the typhoons troposphere have indication significance for the development stage of typhoons. Therefore, the analysis, in particular the quantitative characterization, of the three-dimensional structure inside typhoons is of great significance for typhoons forecasting.
Disclosure of Invention
In order to solve the defects in the prior art, the invention provides an automatic matching and quantitative display method for a typhoon warm center three-dimensional structure of a multi-source-orbit satellite, which realizes the monitoring and display of the typhoon three-dimensional structure based on the polar orbit satellite remote sensing by analyzing multi-source-orbit satellite data such as FY3D, NPP, FY4A and the like, and can check and analyze important macroscopic physical characteristics such as the typhoon warm center structure, typhoon vertical gradient and the like through the display, thereby providing quantitative remote sensing combined diagnosis analysis of different scenes.
In order to achieve the above object, the method for automatically matching and quantitatively displaying a three-dimensional structure of a multi-source orbital satellite typhoon heater according to one embodiment of the present invention includes:
s1, searching typhoon path information to obtain the position of a typhoon real-time center point;
s2, automatically matching the effective multi-source orbit satellite data of the corresponding time according to the position of the typhoon real-time center point;
and S3, drawing a three-dimensional structure slice dynamic diagram and a vertical outline line diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, and quantitatively displaying the three-dimensional structure of the typhoon heating center.
Further, the step S1 of retrieving typhoon path information, obtaining a typhoon real-time center point position includes: and automatically searching typhoon path information, and acquiring the position of a typhoon real-time center point according to the searched typhoon positioning, and further acquiring corresponding time and longitude and latitude information.
Further, the step S2 of automatically matching valid multi-source satellite data corresponding to the time according to the typhoon real-time center point position includes:
s21, determining a validity check area according to the position of the typhoon real-time center point;
s22, acquiring multi-source orbit satellite data of corresponding time, and determining an effective polar orbit satellite data source according to a checking result of the validity checking area;
s23, satellite data of the effective polar orbit satellite data source are automatically matched around the typhoon real-time center point position.
Further, the validity checking area is a rectangular area with the typhoon real-time center point position as the center and the radius of 5 degrees.
Further, step S22 is to acquire multi-source orbital satellite data corresponding to the time, and determine an effective polar orbital satellite data source according to a checking result of the validity checking area, specifically: acquiring multi-source orbit satellite data of corresponding time, sequentially checking the number of effective points of satellite data of each polar orbit satellite data source in the validity checking area range, and determining the satellite data as the valid polar orbit satellite data source if the number of the effective points of the satellite data of each polar orbit satellite data source in the validity checking area range is more than a threshold value containing the total number of the validity checking area ranges.
Further, the polar orbit satellites include FY3D satellites, NPP satellites, and FY4A satellites.
Further, the step S3 of drawing a three-dimensional structure slice dynamic diagram and a vertical contour line diagram of the typhoon heating core according to the matched effective multi-source orbit satellite data, and quantitatively displaying the three-dimensional structure of the typhoon heating core includes:
s31, drawing a three-dimensional structure slice dynamic diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, wherein the dynamic diagram specifically comprises the following steps:
s311, acquiring multidimensional data of an effective polar orbit satellite data source in FY3D satellite data and NPP satellite data, and cutting the multidimensional data into rectangular 2D grid data with a plurality of specified ranges according to the space dimension information of the multidimensional data and the longitudinal direction, the latitudinal direction and the height direction;
s312, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in multiple longitudinal directions to obtain multiple longitudinal 3D static slice diagrams with 0.1 degree interval;
s313, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in the multiple latitudes to obtain multiple latitudinal 3D static slice diagrams with 0.1 degree interval;
s314, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in the multiple height directions, overlapping geographic information, and drawing multiple 3D static slice diagrams of the polar orbit satellite data height level;
s315, acquiring multidimensional data of an effective polar orbit satellite data source, and drawing a 3D graph overlapped with map information.
S316, after all slice images in a certain direction of longitude, latitude and altitude are obtained, sequentially synthesizing GIF dynamic images for all slice images in the same direction;
s317, acquiring multidimensional data of an FY-4A satellite data source corresponding to typhoon time, and drawing a comparison base map;
s318, acquiring all slice images in a certain direction according to three directions, and drawing a three-dimensional structure slice dynamic image of the typhoon heating core, wherein the three-dimensional structure slice dynamic image specifically comprises the following steps: after all slice images in a certain direction are obtained, each slice image is named sequentially, and all slices in the same direction are sequentially synthesized into the GIF dynamic image according to the named sequence through PIL based on Python.
S319, carrying out combined display on multiple pictures:
drawing a comparison base map of FY-4A satellite data and an upper map of FY3D satellite data or a 3D perspective map formed by combining the comparison base map of FY-4A satellite data and the upper map of NPP satellite data through an Axe 3D library, and displaying and combining two graphic surfaces in different colors through setting a coordinate axis range, adding a coordinate axis label and setting different color maps to form the 3D perspective map.
S32, drawing a typhoon warm center vertical profile line structure diagram according to the effective multi-source orbit satellite data.
Further, the threshold is 50%.
Further, step S315 is to acquire multidimensional data of the effective polar orbit satellite data source, and draw a 3D graph superimposed with map information, which specifically includes: acquiring multidimensional data of an effective polar orbit satellite data source, cutting the data into rectangular 2D grid data with longitude and latitude boundaries within an effectiveness check area according to the space dimension of the multidimensional data, drawing the processed grid data through a Python-based mplroolkits library, reading geographic information acquired by a map file in the process of drawing a 3D graph, and superposing the geographic information on the 3D graph through a copy expansion package in a visual display mode to generate the 3D graph with map information.
Further, step S32 is to draw a typhoon warm core vertical profile line structure according to the effective multi-source orbit satellite data, specifically: according to the effective multi-source orbit satellite data, polar orbit satellite data of each layer of height are obtained, regular grid data are obtained through interpolation of a griddata library, corresponding positive temperature range maximum value or temperature difference value is obtained according to longitude and latitude of a typhoon central point, and drawing is conducted to form a typhoon warm center vertical profile line composition comprising a positive temperature range flat profile and a temperature difference profile.
The beneficial effects of the invention are as follows:
1. the invention integrates a plurality of polar orbit satellite data such as FY3D, NPP, FY4A and the like to form complementation, thereby improving typhoon monitoring frequency;
2. according to the invention, by automatically matching the typhoon path with the polar orbit satellite data, the three-dimensional structure slice dynamic diagram of the typhoon heating core and the vertical profile line diagram of the typhoon heating core are drawn, so that the three-dimensional display of the quantitative typhoon development intensity is realized;
3. according to the invention, a threshold value of validity check is set through an automatic matching algorithm of a validity check area, if the remote sensing data of a certain polar orbit satellite is automatically matched with the threshold value (more than 50%), the data of the certain polar orbit satellite is automatically drawn, if the remote sensing data of the certain polar orbit satellite is lower than the threshold value, the drawing is not performed according to incomplete measurement or observation, and therefore the display quality of a three-dimensional structure slice dynamic diagram of a typhoon warm center and a vertical profile line diagram of the typhoon warm center is ensured;
4. according to the invention, the drawn graphic products are shared in an HTTP service mode, after typhoon path points are overlapped on different base maps of tropical cyclone or satellite cloud patterns through position information, whether animation images exist at the current path points or not is displayed in a bubble mode, if the animation images exist, on-line viewing of the 3D animation images can be realized by clicking live typhoons, so that quantitative remote sensing combined diagnosis analysis of different scenes is provided around cloud systems, thermal power, dynamic and water vapor conditions developed by typhoons.
Drawings
FIG. 1 is a flow chart of a method for automatically matching and quantitatively displaying a three-dimensional structure of a multi-source orbit satellite typhoon warm center;
FIG. 2 is a view of a plurality of 3D static slices spaced 0.1 degrees apart in the warp direction in accordance with an embodiment of the present invention;
FIG. 3 is a view of a plurality of 3D static slices spaced 0.1 degrees apart in the weft direction in accordance with an embodiment of the present invention;
FIG. 4 is a 3D static slice view of 4 height levels in height direction according to one embodiment of the invention;
FIG. 5 is a dynamic view of the GIF synthesized in three directions of warp, weft and height according to an embodiment of the invention;
FIG. 6 is a remote sensing temperature range plan of typhoon No. 03 "Siamese" 2022 in accordance with one embodiment of the present invention;
FIG. 7 is a diagram of typhoon No. 03 "remote kibble" light Wen Juping of 2022 according to an embodiment of the present invention;
FIG. 8 is a diagram of a typhoon warm core vertical profile line of one embodiment of the present invention;
FIG. 9 is a schematic diagram of a typhoon warm core structure characterizing typhoon strength in accordance with one embodiment of the present invention;
FIG. 10 is a comparison of remote sensing of NPP and FY-3D microwave thermometer warm core structures in accordance with an embodiment of the invention;
FIG. 11 is a typhoon three-dimensional stereo structure technical algorithm flow chart of one embodiment of the present invention;
FIG. 12 is a schematic representation of a three-dimensional structure of a multi-source satellite on a typhoon path according to one embodiment of the present invention.
In the figure: 1-warp 3D static slice; 2-weft 3D static slice diagram; 3-height direction 3D static slice; 4-base drawing; 5-a warm heart color code; 6-normal temperature distance flat profile; 7-temperature difference profile; 8-bubble point; 9-teli typhoon path; 10-Du Surui typhoon path.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the following description will be made with reference to fig. 1 to 12 and examples.
It is noted that relational terms such as first and second, and the like are used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Moreover, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
The warm core structure is the most obvious characteristic of typhoons different from temperate zone cyclone, and the establishment of the high-altitude warm core structure is one of main marks for typhoons. The relevant scholars think that the warm-core structure of typhoons can be extended to the whole troposphere in the vertical direction, and the warm-core strength is also maximum when the warm-core structure reaches the maximum near the middle-high layer 250 hPa. The central heating temperature distance before logging reaches the maximum value, and the central heating structure is smooth and deep in shape; after logging in, the temperature of the warm center gradually decreases from a flat value, the structural shape of the warm center is irregular, the thickness is reduced, and the height is reduced. There is a significant humidity and temperature gradient inside typhoons. Therefore, the three-dimensional warm-center vertical structure of typhoons is studied, the change characteristics of the three-dimensional warm-center structure of typhoons in different development processes can be seen, and a scientific monitoring means is provided for typhoons forecast analysis.
As an embodiment, as shown in fig. 1 and 11, the invention provides a method for automatically matching and quantitatively displaying a three-dimensional structure of a multi-source-orbit satellite typhoon warm center, which comprises the following steps:
s1, searching typhoon path information, and acquiring typhoon time, typhoon longitude and latitude information and typhoon real-time center point position;
the typhoons are retrieved by the urllib module automatically retrieving typhoon path information per hour. For each typhoon retrieved, the serial number, chinese name, longitude and latitude information, typhoon time and the like corresponding to the typhoon are recorded respectively, and the typhoon real-time center point position of the typhoon positioning report can be obtained according to the retrieval. Typhoon time is the time corresponding to the location of the typhoon real-time center point being retrieved.
Typhoon time varies with the movement of the location of the typhoon's real-time center point, typically calculated from the initial stage of typhoon formation and finally to typhoon dissipation time, which includes: typhoon formation time: typhoons are formed at a time point when typhoons begin to form, and generally refer to a process of developing and evolving a tropical low-voltage system into typhoons; typhoon active period: the active period of typhoons refers to the period in which typhoons persist and have strong wind power and precipitation activities. This period of time may last from days to weeks, depending on the characteristics of the typhoons and the environmental conditions; typhoon dissipation time: the dissipation time of typhoons refers to the point in time at which typhoons gradually fade and eventually disintegrate. The process of typhoon dissipation typically involves reduced wind forces, reduced precipitation, and collapse of the typhoon structure.
S2, matching the effective multi-source orbit satellite data corresponding to typhoon time according to the position of the typhoon real-time center point;
the polar orbit satellite can dynamically scan multiple areas, the function is realized through the orthogonal rotation of the earth and the motion of the polar orbit satellite, the image resolution is higher, and the polar orbit satellite is not overlapped with a typhoon path at any time, so that the automation is realized.
The multi-source orbital satellite data comprises data of FY3D satellites and NPP satellites, wherein:
the FY3D satellite refers to a China Fengyun No. three D polar orbit satellite, which is a polar orbit satellite for meteorological observation. The FY3D satellite carries a number of instruments including a multi-channel satellite scanning radiometer (MWTS), a microwave hygrometer (MWHS), a multi-channel satellite scanning dry balloon (MWRI), a Visible and Infrared Imaging Radiometer (VIIRS), and the like. The data of the FY3D satellite may provide approximately 43 levels of atmospheric and surface observations including, but not limited to, temperature, humidity, cloud cover, ocean surface temperature, vegetation index, and the like.
The NPP satellite is a polar orbit satellite developed by the united states national weather service (NOAA) and NASA, and is mainly used for global weather observation. The NPP satellite carries a plurality of sensors, including PTemp (Pathfinder Temperature) sensors. The PTemp sensor is a meteorological sensor on an NPP satellite and is mainly used for measuring the surface temperature. The method can acquire the surface temperature distribution condition in the global scope through a remote sensing technology. The PTemp sensor may provide high resolution surface temperature data of approximately 100 levels, including terrestrial and marine temperature information. The data has important application value in the fields of weather forecast, climate research, environmental monitoring and the like.
FY-4A (Fengyun No. A star) is a static weather polar orbit satellite independently developed in China, belongs to the fourth of the Fengyun series polar orbit satellites and is taken as an important component of a Chinese weather polar orbit satellite system. FY-4A satellites carry a plurality of sensors and instrumentation for acquiring various weather information of the earth's atmosphere and the earth's surface. The data has important application value in the aspects of weather forecast, climate monitoring, disaster early warning and the like. The data for the FY-4A satellite includes visible and infrared images: the images can be used for observing information such as cloud morphology, cloud quantity, cloud top temperature and the like, and are very important for analysis and prediction of a weather system.
FY3D, FY4A, NPP data names are (by way of example):
·FY3D_TSHSX_ORBT_L2_AVP_MLT_NUL_20220630_1701_033KM_MS.HDF
·NPR-MIRS-SND_v11r1_NPP_s202206301756080_e202206301808236_c202206301928570.nc
·FY4A-_AGRI--_N_DISK_1047E_L1-_FDI-_MULT_NOM_20220901040000_
20220901041459_4000M_V0001.HDF
further, the effective multi-source orbit satellite data corresponding to typhoon time is matched according to the typhoon real-time center point position, comprising the following steps:
s21, determining a validity check area according to the position of the typhoon real-time center point;
according to the position of the typhoon real-time center point, a validity check area is determined, and the validity check area is a rectangular area with the position of the typhoon center point as the center and the radius of 5 degrees. The radius of 5 degrees on the earth is about 550 km, and the range is that the structure of the typhoon is moderate to observe and display, and too large or too small can lead the typhoon to look smaller or not see the whole view of the typhoon.
Since the earth is a three-dimensional object that approximates an ellipsoid, the distance between latitude and longitude is not exactly the same, and an approximate calculation method is as follows:
determining a latitude range: the latitude of the earth ranges from 90 degrees in south latitude to 90 degrees in north latitude, and a rectangular range of 5 degrees radius means that the distance in the vertical direction should be 10 degrees (5 degrees on both north and south sides). Assuming that the latitude of the center point is lat_center, the latitude boundary of the range may be calculated as lat_min=lat_center-5 degrees, lat_max=lat_center+5 degrees.
Determining a longitude range: the longitude ranges from 180 degrees in the west to 180 degrees in the east, a rectangular range of 5 degrees radius means that the distance in the horizontal direction should be 10 degrees (5 degrees on each east and west side). Assuming that the longitude of the center point is lon_center, the longitude boundary of the range may be calculated as lon_min=lon_center-5 degrees, lon_max=lon_center+5 degrees.
S22, acquiring multi-source orbit satellite data of corresponding time, and determining an effective polar orbit satellite data source according to a checking result of the validity checking area;
and respectively acquiring FY3D satellite data, NPP satellite data and FY-4A satellite data according to typhoon corresponding time, and sequentially checking the number of effective points of satellite data of each polar orbit satellite data source in the validity checking area range. The method comprises the following steps: acquiring observation range data points of a satellite sensor: observation range data points for the satellite sensors are acquired, which typically represent geographic location information that the satellite sensors can cover. These data points may be latitude and longitude coordinates. Defining a rectangular range of a validity check area: and determining a rectangular range of the validity checking area, and searching for an intersection with the satellite sensor observation range in the range. The rectangular range may be defined using upper and lower bound longitude and latitude coordinates. Finding an intersection: the satellite sensor's observation range data points are compared to a rectangular range to find the intersection between the two. Geospatial computing methods, such as point-within-polygon determination algorithms, may be used to determine whether each data point is within a rectangular range. Extracting intersection data points: data points located at the intersection portion, which are data points of the observation range of the satellite sensor covered by the rectangular range, are extracted. If the total number of the validity check area ranges is more than 50%, the valid polar orbit satellite data source is determined. In the invention, if the remote sensing data is automatically matched with more than 50% of the threshold value as in the above case, the remote sensing data is automatically mapped, and is processed (not mapped) according to the lack of measurement or observation.
S23, satellite data of the effective polar orbit satellite data source are automatically matched around the typhoon real-time center point position.
S3, drawing a three-dimensional structure slice dynamic diagram and a vertical outline line diagram of the typhoon heater core according to the matched effective multi-source orbit satellite data, quantitatively displaying the three-dimensional structure of the typhoon heater core, and comprising the following steps:
s31, drawing a three-dimensional structure slice dynamic diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, wherein the dynamic diagram specifically comprises the following steps:
s311, acquiring multidimensional data of an effective polar orbit satellite data source in an FY3D satellite and an NPP satellite, and cutting the multidimensional data into rectangular 2D grid data with a plurality of specified ranges according to the space dimension information of the multidimensional data and the longitudinal direction, the latitudinal direction and the height direction;
acquiring multidimensional data of an effective polar orbit satellite data source in an FY3D satellite and an NPP satellite, cutting polar orbit satellite data with longitude and latitude boundaries within a range of an effectiveness check area and with heights within a polar orbit satellite data hierarchy into a plurality of flaky rectangular 2D grid data according to the polar orbit satellite data hierarchy at intervals of 0.1 degree in the longitude direction, 0.1 degree in the latitude direction and the height direction respectively, wherein the FY3D satellite has a height direction of 43 layers and the NPP satellite data has a height direction of 100 layers (similar to a cube with one longitude and latitude boundary within the range of the effectiveness check area and with heights within the polar orbit satellite data hierarchy, cutting the polar orbit satellite data into a plurality of rectangular 2D grid data at intervals of 0.1 degree in the longitude direction, cutting the polar orbit satellite data into a plurality of rectangular 2D grid data at intervals of 0.1 degree in the latitude direction and cutting the polar orbit satellite data into a plurality of rectangular 2D grid data according to the height direction).
S312, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in a plurality of longitudinal directions by using a Matplotlib library based on Python, and drawing a plurality of longitudinal 3D static slices with 0.1 degree intervals as shown in a figure 2, wherein a warm center color scale 5 for marking temperature colors is additionally attached;
s313, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in multiple latitudes based on a Matplotlib library of Python, and drawing a plurality of latitudinal 3D static slices with 0.1 degree intervals as shown in a figure 3, wherein a warm center color scale 5 for marking temperature colors is additionally attached;
s314, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in a plurality of height directions based on a Matplotlib library of Python, reading geographic information obtained by a map file in the process of drawing a graph, superposing the geographic information on a 3D graph in a visual display mode based on a support extension package of Python, and drawing 43 3D static slice diagrams in the FY3D satellite height direction and 100 3D static slice diagrams in the NPP satellite data height direction. As shown in fig. 4, the 3D static slice is shown in fig. 3 for 4 height directions, and a warm center color mark 5 for marking the color of the temperature is additionally attached.
S315, acquiring multi-dimensional data of an effective polar orbit satellite data source, and cutting the data into rectangular 2D grid data with an effective range (longitude and latitude boundaries are in a validity check area range and heights are in a polar orbit satellite data hierarchy) according to the space dimension of the multi-dimensional data; and finally, drawing the processed grid data through a Python-based mplroolkits library, reading geographic information obtained by a map file in the process of drawing the 3D graph, superposing the geographic information on the 3D graph through a Carcopy extension package in a visual display manner, and finally generating the 3D graph with the map information.
S316, acquiring multidimensional data of an FY-4A satellite data source corresponding to typhoon time, and drawing a comparison base map 4;
FY-4A (Fengyun No. A star) is a static weather polar orbit satellite independently developed in China, belongs to the fourth of the Fengyun series polar orbit satellites and is taken as an important component of a Chinese weather polar orbit satellite system. FY-4A satellites carry a plurality of sensors and instrumentation for acquiring various weather information of the earth's atmosphere and the earth's surface. The data has important application value in the aspects of weather forecast, climate monitoring, disaster early warning and the like. The data for the FY-4A satellite includes visible and infrared images: the images can be used for observing information such as cloud morphology, cloud quantity, cloud top temperature and the like, and are very important for analysis and prediction of a weather system.
The method comprises the steps of obtaining multidimensional data of an FY-4A satellite data source corresponding to typhoon time, obtaining visible light image data through an FY4A polar orbit satellite during daytime, cutting the data, processing the data to form three channel data of R, B and G, and fusing the three channels of R, G and B into a true color map through a cv2 library.
At night, infrared image data are acquired through FY4A, the data are cut, regular two-dimensional data are formed through processing, and the two-dimensional grid data are drawn into an infrared image.
S317, obtaining all slice images in a certain direction according to three directions, and drawing a three-dimensional structure slice dynamic image of the typhoon heating core, wherein the three-dimensional structure slice dynamic image specifically comprises the following steps:
after all slice images in a certain direction are obtained, each slice image is named sequentially, and all slices in the same direction are sequentially synthesized into the GIF dynamic image shown in figure 5 according to the named order through a PIL (Python Image Library) third party library based on Python.
S319, carrying out combined display on multiple pictures:
drawing a comparison base map 4 of FY-4A satellite data and an upper map of FY3D satellite data or a 3D perspective map formed by combining the comparison base map 4 of FY-4A satellite data and the upper map of NPP satellite data through an Axe 3D library, and displaying and combining two graphic surfaces in different colors through setting a coordinate axis range, adding a coordinate axis label and setting different color maps (cmap) to form the 3D perspective map.
The upper graph of FY3D satellite data and the upper graph of NPP satellite data can be slice static graphs or slice dynamic graphs.
According to the difference of polar orbit satellite data adopted in drawing, the three-dimensional structure diagram of the typhoon heating center can be a temperature diagram, a temperature range flat diagram, a bright temperature diagram, a bright Wen Juping diagram, a humidity diagram and the like.
And (3) performing range-level calculation on polar orbit satellite data: the average value of the effective points in the drawing range is calculated first, and then the average value is subtracted from the value of each grid point and recorded as a distance flat value.
The temperature range plan provides a visual understanding of typhoon structure by showing deviations from the long term average temperature. As shown in FIG. 6, which is a remote sensing temperature range plan for typhoons No. 03 of 2022, we can see the higher temperatures in the typhoons core area and the lower temperatures around the typhoon circulation system by observing the temperature range plan. The obvious temperature difference can clearly show the position of the surface typhoon heating center.
In the past, due to the limitation of the offshore observation means, the typhoon warm center structure is not researched by lacking enough observation data, and the temperature of the typhoon warm center structure can be calculated according to the intensity of electromagnetic radiation emitted by a celestial body along with the development of polar orbit satellite remote sensing technology, which is also called as bright temperature. In the remote sensing field, the bright temperature is also used for calculating the temperature distribution or other environmental parameters of the earth surface, which greatly compensates for the deficiency of the offshore observation data. As shown in figure 7, a typhoon No. 03 "yaoba" light Wen Juping in 2022 can be seen, and the characteristics of the change of the warm-center three-dimensional structure of typhoons in different development processes can be seen, so that a scientific monitoring means is provided for typhoons forecast analysis.
S32, drawing a typhoon warm center vertical profile line structure diagram according to effective multi-source orbit satellite data, wherein the structure diagram specifically comprises the following steps:
the typhoon warm core vertical profile line structure diagram can comprise a positive temperature distance flat profile 6 and a temperature difference profile 7. Wherein: the positive temperature distance flat profile 6 is a positive temperature distance maximum value vertical profile within a range of 1 degree of the typhoon center, and the temperature difference profile 7 is a difference vertical profile of a mean value within a range of 1 degree of the typhoon center radius and a mean value of the temperature within a range of 5 degrees of the radius. The profile is used for quantitatively checking the change of the warm heart intensity and the vertical distribution, and compared with the three-dimensional structure diagram, the three-dimensional structure diagram is qualitatively checked, and the profile diagram is quantitatively checked.
Positive temperature range refers to the positive value of the difference between the average air temperature in a region and the average air temperature over the years, over a period of time (typically a month or a quarter). The positive temperature range maximum value indicates what the positive temperature range maximum value of the region is during the period of time. The step of calculating the positive temperature distance maximum value is as follows: collecting the historical contemporaneous air temperature data of the region, and calculating the average air temperature of the historical contemporaneous air temperature; collecting air temperature data of the region in the period of time, and calculating the average air temperature in the period of time; calculating the positive temperature range in the period of time to be flat, namely subtracting the average temperature in the same period of history from the average temperature in the period of time, and taking the value if the result is positive, otherwise, taking the value as 0; repeating the step 3, calculating the positive temperature distance flat maximum value of each time period, and comparing to obtain the positive temperature distance flat maximum value.
According to the invention, after FY3D or NPP data files are read, polar orbit satellite data of each layer of height (FY 3D43 layer and NPP100 layer) are obtained according to effective multi-source polar orbit satellite data. Interpolation is carried out through the griddata library to obtain regular grid data, and corresponding positive temperature distance maximum value or temperature difference value is obtained according to longitude and latitude of a typhoon central point to draw a typhoon warm center vertical profile line graph shown in figure 8. Namely, according to the longitude and latitude of the typhoon central point, a corresponding positive temperature distance flat maximum value is obtained to draw a positive temperature distance flat profile 6, a corresponding temperature difference value is obtained to draw a temperature difference profile 7, and a typhoon warm center vertical profile line structure diagram shown in figure 8 is formed.
The three-dimensional structure automatic matching and quantitative display method of the multi-source-orbit satellite typhoon heating center can check and analyze important macroscopic physical characteristics such as the typhoon heating center structure, the typhoon vertical gradient and the like, and the typhoon heating center structure and the vertical profile can quantify the typhoon development intensity.
Typhoon warm heart structure can represent typhoon intensity: as shown in FIG. 9, the stronger warm-core structure is 7 months and 2 days (TY) of a certain year, the weaker warm-core structure is 6 months and 30 days (TS) of a certain year and 3 days (TS) of a certain year, and the warm-core structure height layer mainly occurs between 500hPa and 250 hPa.
Remote sensing monitoring and comparison of NPP and FY-3D microwave thermometer warm core structure: as shown in figure 10, the NPP has better remote sensing monitoring effect than the warm-core structure of the FY-3D microwave thermometer, and the data precision of the FY-3D microwave thermometer needs to be further improved.
As another aspect of the implementation of the invention, the invention also provides an online display function of the typhoon three-dimensional structure, namely, the drawn graphic products are shared in an HTTP service mode, after typhoon path points are overlapped on different display base graphs on a remote sensing satellite data display and analysis platform of an area, whether animation graphs exist on the current path points or not is displayed in a bubble mode, if so, the 3D animation graphs can be checked online by clicking the typhoon live bubble point 8. The show floor map may include a cloud top temperature map, a cloud top elevation map, a Yun Xiangtai map, a sea level temperature map, a cloud top barometric map, a cloud optical thickness map, a precipitation estimate map, a typhoon identification map, and a stationary cloud map and polar orbit cloud map of the show satellite cloud map class from different satellite (which may be FY4A/AGRI, himawari, FY B/AGR, FY2G/VISSR, and FY 2F/VISSR) sources. Different display base charts can provide different application scenes for typhoon analysis, for example, a cloud top temperature chart can analyze the convection development intensity of typhoons; the cloud top height map can analyze the development height of typhoon spiral cloud system; the sea surface temperature map can analyze the energy conditions of typhoons, and the invention can provide quantitative remote sensing combination diagnosis analysis of different scenes around the cloud systems, heat, power and water vapor conditions of typhoons.
As shown in fig. 12, taking the on-line display of remote sensing satellite data of typhoon "taili" and "Du Surui" in the south China area as an example, a base map showing tropical cyclone is provided on the analysis platform (which can be changed and adjusted according to the requirement), on the base map, according to the position of the real-time center point of the typhoon detected, according to the typhoon time, position and typhoon name, the position of the real-time center point of the typhoon is overlapped and drawn on the display base map, all the positions of the real-time center points of the typhoon are connected in series to form a typhoon path, after the taili typhoon paths 9 and Du Surui are overlapped with the typhoon path 10, whether a three-dimensional structure slice dynamic map of the typhoon warm center and a vertical profile line structure of the typhoon warm center exist at the current path point are displayed in a bubble mode, if the current path point has the map, the three-dimensional structure of the typhoon warm center can be checked on line by clicking the typhoon path bubble point 8, the cloud system, the thermal power, the conditions around the development of the typhoon can be visually seen, and the three-dimensional structure change characteristics of the typhoon in different development processes are intuitively monitored, and a scientific remote sensing monitoring means is provided for the analysis and prediction of the typhoon warm center.
Alternatively, in this embodiment, it will be understood by those skilled in the art that all or part of the steps in the methods of the above embodiments may be performed by a program for instructing a terminal device to execute the steps, where the program may be stored in a computer readable storage medium, and the storage medium may include: flash disk, read-Only Memory (ROM), random-access Memory (Random Access Memory, RAM), magnetic or optical disk, and the like.
Although embodiments of the present invention have been shown and described, it will be understood by those skilled in the art that various changes, modifications, substitutions and alterations can be made therein without departing from the principles and spirit of the invention, the scope of which is defined in the appended claims and their equivalents.
The invention and its embodiments have been described above with no limitation, and the actual construction is not limited to the embodiments of the invention as shown in the drawings. In summary, if one of ordinary skill in the art is informed by this disclosure, a structural manner and an embodiment similar to the technical solution should not be creatively devised without departing from the gist of the present invention.

Claims (10)

1. The automatic matching and quantitative display method for the three-dimensional structure of the typhoon warm center of the multi-source orbital satellite is characterized by comprising the following steps of:
s1, searching typhoon path information to obtain the position of a typhoon real-time center point;
s2, automatically matching the effective multi-source orbit satellite data of the corresponding time according to the position of the typhoon real-time center point;
and S3, drawing a three-dimensional structure slice dynamic diagram and a vertical outline line diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, and quantitatively displaying the three-dimensional structure of the typhoon heating center.
2. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warm center of the multi-source satellite according to claim 1, wherein the step S1 of retrieving typhoon path information and obtaining the position of the typhoon real-time center point comprises the following steps:
and automatically searching typhoon path information, and acquiring the position of a typhoon real-time center point according to the searched typhoon positioning, and further acquiring corresponding time and longitude and latitude information.
3. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warm center of the multi-source satellite according to claim 1, wherein the step S2 of automatically matching the valid multi-source satellite data of the corresponding time according to the position of the typhoon real-time center point comprises the following steps:
s21, determining a validity check area according to the position of the typhoon real-time center point;
s22, acquiring multi-source orbit satellite data of corresponding time, and determining an effective polar orbit satellite data source according to a checking result of the validity checking area;
s23, satellite data of the effective polar orbit satellite data source are automatically matched around the typhoon real-time center point position.
4. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warm center of the multi-source satellite according to claim 3, wherein the effectiveness checking area is a rectangular area with a radius of 5 degrees and with the position of the real-time center point of the typhoon as the center.
5. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warming center of the multi-source-orbit satellite according to claim 3, wherein the step S22 is to acquire multi-source-orbit satellite data corresponding to time, and determine an effective polar-orbit satellite data source according to a checking result of the validity checking area, specifically comprising the following steps:
acquiring multi-source orbit satellite data of corresponding time, sequentially checking the number of effective points of satellite data of each polar orbit satellite data source in the validity checking area range, and determining the satellite data as the valid polar orbit satellite data source if the number of the effective points of the satellite data of each polar orbit satellite data source in the validity checking area range is more than a threshold value containing the total number of the validity checking area ranges.
6. The method for automatically matching and quantitatively displaying a three-dimensional structure of a typhoon warm center of a multi-source orbiting satellite according to claim 1, wherein the polar orbiting satellites comprise FY3D satellites, NPP satellites and FY4A satellites.
7. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon-heating core of the multi-source-orbit satellite according to claim 1, wherein the step S3 of drawing a slice dynamic diagram and a vertical contour line diagram of the three-dimensional structure of the typhoon-heating core according to the matched effective multi-source-orbit satellite data for quantitatively displaying the three-dimensional structure of the typhoon-heating core comprises the following steps:
s31, drawing a three-dimensional structure slice dynamic diagram of the typhoon heating center according to the matched effective multi-source orbit satellite data, wherein the dynamic diagram specifically comprises the following steps:
s311, acquiring multidimensional data of an effective polar orbit satellite data source in FY3D satellite data and NPP satellite data, and cutting the multidimensional data into rectangular 2D grid data with a plurality of specified ranges according to the space dimension information of the multidimensional data and the longitudinal direction, the latitudinal direction and the height direction;
s312, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in multiple longitudinal directions to obtain multiple longitudinal 3D static slice diagrams with 0.1 degree interval;
s313, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in the multiple latitudes to obtain multiple latitudinal 3D static slice diagrams with 0.1 degree interval;
s314, sequentially carrying out slice drawing on the obtained rectangular 2D grid data in the multiple height directions, overlapping geographic information, and drawing multiple 3D static slice diagrams of the polar orbit satellite data height level;
s315, acquiring multidimensional data of an effective polar orbit satellite data source, and drawing a 3D graph overlapped with map information;
s316, after all slice images in a certain direction of longitude, latitude and altitude are obtained, sequentially synthesizing GIF dynamic images for all slice images in the same direction;
s317, acquiring multidimensional data of an FY-4A satellite data source corresponding to typhoon time, and drawing a comparison base map;
s318, acquiring all slice images in a certain direction according to three directions, and drawing a three-dimensional structure slice dynamic image of the typhoon heating core, wherein the three-dimensional structure slice dynamic image specifically comprises the following steps:
after all slice images in a certain direction are obtained, each slice image is named sequentially, and all slices in the same direction are sequentially synthesized into a GIF dynamic image according to the named sequence through PIL based on Python;
s319, carrying out combined display on multiple pictures:
drawing a comparison base map of FY-4A satellite data and an upper map of FY3D satellite data or a 3D perspective map formed by combining the comparison base map of FY-4A satellite data and the upper map of NPP satellite data through an Axe 3D library, displaying and combining two graphic surfaces in different colors through setting a coordinate axis range, adding a coordinate axis label and setting different color mapping to form a 3D perspective map;
s32, drawing a typhoon warm center vertical profile line structure diagram according to the effective multi-source orbit satellite data.
8. The method for automatically matching and quantitatively displaying a three-dimensional structure of a multi-source satellite typhoon warm center according to claim 5, wherein the threshold is 50%.
9. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warm center of the multi-source satellite according to claim 7, wherein the step S315 is to obtain multi-dimensional data of the effective polar orbit satellite data source, and draw a 3D graph superimposed with map information, specifically:
acquiring multidimensional data of an effective polar orbit satellite data source, cutting the data into rectangular 2D grid data with longitude and latitude boundaries within an effectiveness check area according to the space dimension of the multidimensional data, drawing the processed grid data through a Python-based mplroolkits library, reading geographic information acquired by a map file in the process of drawing a 3D graph, and superposing the geographic information on the 3D graph through a copy expansion package in a visual display mode to generate the 3D graph with map information.
10. The method for automatically matching and quantitatively displaying the three-dimensional structure of the typhoon warm center of the multi-source satellite according to claim 7, wherein the step S32 is to draw a vertical outline structure of the typhoon warm center according to the effective multi-source satellite data, specifically:
according to the effective multi-source orbit satellite data, polar orbit satellite data of each layer of height are obtained, regular grid data are obtained through interpolation of a griddata library, corresponding positive temperature range maximum value or temperature difference value is obtained according to longitude and latitude of a typhoon central point, and drawing is conducted to form a typhoon warm center vertical profile line composition comprising a positive temperature range flat profile and a temperature difference profile.
CN202311463858.7A 2023-11-06 2023-11-06 Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite Active CN117422832B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311463858.7A CN117422832B (en) 2023-11-06 2023-11-06 Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311463858.7A CN117422832B (en) 2023-11-06 2023-11-06 Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite

Publications (2)

Publication Number Publication Date
CN117422832A true CN117422832A (en) 2024-01-19
CN117422832B CN117422832B (en) 2024-04-12

Family

ID=89524662

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311463858.7A Active CN117422832B (en) 2023-11-06 2023-11-06 Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite

Country Status (1)

Country Link
CN (1) CN117422832B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006208195A (en) * 2005-01-28 2006-08-10 National Research Institute For Earth Science & Disaster Provention Device and method for estimating three-dimensional distribution of rainfall intensity and quantity
CN106443830A (en) * 2016-06-16 2017-02-22 杭州师范大学 Method for typhoon monitoring and evaluation of monitoring precision based on multi-source satellite data
CN112614230A (en) * 2021-01-06 2021-04-06 飞燕航空遥感技术有限公司 Typhoon dynamic path drawing method and drawing system
CN114896544A (en) * 2022-03-15 2022-08-12 国家卫星气象中心(国家空间天气监测预警中心) Three-dimensional cloud fusion analysis method for cooperation of wind cloud meteorological satellite and multiple data sources
CN116089546A (en) * 2021-11-05 2023-05-09 中国科学院深圳先进技术研究院 Typhoon cloud system identification method, typhoon cloud system identification system, typhoon cloud system identification terminal and storage medium
CN116931009A (en) * 2023-07-20 2023-10-24 国家卫星气象中心(国家空间天气监测预警中心) Method and equipment for drawing temperature logarithmic pressure diagram by satellite combined observation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006208195A (en) * 2005-01-28 2006-08-10 National Research Institute For Earth Science & Disaster Provention Device and method for estimating three-dimensional distribution of rainfall intensity and quantity
CN106443830A (en) * 2016-06-16 2017-02-22 杭州师范大学 Method for typhoon monitoring and evaluation of monitoring precision based on multi-source satellite data
CN112614230A (en) * 2021-01-06 2021-04-06 飞燕航空遥感技术有限公司 Typhoon dynamic path drawing method and drawing system
CN116089546A (en) * 2021-11-05 2023-05-09 中国科学院深圳先进技术研究院 Typhoon cloud system identification method, typhoon cloud system identification system, typhoon cloud system identification terminal and storage medium
CN114896544A (en) * 2022-03-15 2022-08-12 国家卫星气象中心(国家空间天气监测预警中心) Three-dimensional cloud fusion analysis method for cooperation of wind cloud meteorological satellite and multiple data sources
CN116931009A (en) * 2023-07-20 2023-10-24 国家卫星气象中心(国家空间天气监测预警中心) Method and equipment for drawing temperature logarithmic pressure diagram by satellite combined observation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王刚 等: "广东省气象卫星遥感数据资源池设计与应用", 气象科技, vol. 50, no. 4, 31 August 2022 (2022-08-31) *
王尚宏 等: "登陆过程中台风高层暖心结构演变特征分析", 气象, vol. 44, no. 5, 31 May 2018 (2018-05-31) *

Also Published As

Publication number Publication date
CN117422832B (en) 2024-04-12

Similar Documents

Publication Publication Date Title
Barrett The estimation of monthly rainfall from satellite data
Unger Connection between urban heat island and sky view factor approximated by a software tool on a 3D urban database
Du et al. Simultaneous investigation of surface and canopy urban heat islands over global cities
Matzarakis et al. Sky view factor as a parameter in applied climatology-rapid estimation by the SkyHelios model
CN111553245A (en) Vegetation classification method based on machine learning algorithm and multi-source remote sensing data fusion
An et al. Three-dimensional point cloud based sky view factor analysis in complex urban settings.
Sun-Mack et al. Regional apparent boundary layer lapse rates determined from CALIPSO and MODIS data for cloud-height determination
CN113343777B (en) Satellite rapid planning method and system for natural disaster emergency response
CN111241698A (en) Convective stratum top wind field inversion method and device, storage medium and computer equipment
Chen et al. Modeling of mean radiant temperature based on comparison of airborne remote sensing data with surface measured data
Kidd et al. Derivation of sky-view factors from lidar data
Mustaquim Utilizing remote sensing data and ArcGIS for advanced computational analysis in land surface temperature modeling and land use property characterization
Slättberg et al. Spatial and temporal patterns of planetary boundary layer height during 1979–2018 over the Tibetan Plateau using ERA5
Chen et al. Quantifying the carbon storage in urban trees using multispectral ALS data
Liu et al. Microclimatic measurements in tropical cities: Systematic review and proposed guidelines
Menon Remote sensing applications in agriculture and forestry
CN117422832B (en) Automatic matching and quantitative display method for typhoon warm center three-dimensional structure of multi-source orbit satellite
Voogt Image representations of complete urban surface temperatures
Erickson Comparison of developing vs. non-developing tropical disturbances
CN112629666B (en) Method suitable for constructing surface temperature daily transformation model of thermal infrared unmanned aerial vehicle
Caspari Tracking the Cold: Remote Sensing for Glacial Archaeology
Sefercik et al. Point-based and model-based geolocation analysis of airborne laser scanning data
Santosh Estimation of daytime planetary boundary layer height (PBLH) over the tropics and subtropics using COSMIC-2/FORMOSAT-7 GNSS–RO measurements
Ni et al. A method for the registration of multiview range images acquired in forest areas using a terrestrial laser scanner
Wang et al. Ground sampling methods for surface soil moisture in heterogeneous pixels

Legal Events

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