CN106934793A - Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support - Google Patents
Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support Download PDFInfo
- Publication number
- CN106934793A CN106934793A CN201511032319.3A CN201511032319A CN106934793A CN 106934793 A CN106934793 A CN 106934793A CN 201511032319 A CN201511032319 A CN 201511032319A CN 106934793 A CN106934793 A CN 106934793A
- Authority
- CN
- China
- Prior art keywords
- spatial modeling
- water
- nuclear power
- waterfrom
- power plant
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10048—Infrared image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Radiation Pyrometers (AREA)
- Image Processing (AREA)
Abstract
The invention belongs to spatial modeling technical field, and in particular to the thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under a kind of spatial modeling technical support；During present invention aim to address monitoring thermal discharge waterfrom nuclear power plant situation with thermal infrared satellite remote-sensing image at present, complex operation, manual intervention is excessive, lacks efficiently this problem of implementation；The present invention is comprised the following steps：Step one：Data prepare；Step 2：Geometric accurate correction；Step 3：Cloud is recognized；Step 4：Land and water separates；Step 5：Spatial modeling is carried out with " mono window algorithm "；Step 6：Temperature anomaly is charted；Effect of the invention is：Artificial point position observation at present and both warm water discharge monitoring methods of aviation flight are contrasted, the method for general artificial point position observed temperature, certainty of measurement is high, can accurately reflect influence of the thermal discharge waterfrom nuclear power plant to water proximate.
Description
Technical field
The invention belongs to digital image processing techniques field, and in particular under a kind of spatial modeling technical support
Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method.
Background technology
The research of thermal discharge waterfrom nuclear power plant effect on environment is investigated for preventing thermal pollution, protection sea water quality and
Ecological environment is significant.Finite number website can only be obtained using traditional sea temperature investigation method
Temperature data, and have that area coverage is big, frequency is high using satellite data, can be synchronous in the world
The advantages such as observation, can draw the temperature conditions of water field of big area with inverting.Using the number of remote sensing sea temperature inverting
Mainly there is AVHRR and MODIS according to source, they have the Detection Using Thermal Infrared Channel of more than 2, main use is divided
Split-window algorithm inverting temperature.But the resolution ratio of this two satellite datas is relatively low (1000m), is not suitable for core
The electric less waters of warm water discharge influence area.
Therefore, the present invention is using sensor (TIRS, full name Thermal carried on Landsat-8 stars
Infrared Sensor, thermal infrared sensor), two Thermal infrared bands (10.6~11.2 μm,
11.5~12.5 μm), the warm water discharge to nuclear power station carries out remote sensing information extraction.Nuclear power station is extracted in satellite remote sensing
Warm water discharge work target be to realize day-to-day, rapid and facilitation, for this application target,
The present invention has designed and Implemented a kind of thermal discharge waterfrom nuclear power plant Satellite Remote Sensing side based on spatial modeling technology
The corresponding wave band of satellite thermal infrared data need to only be imported software by method, this method, computing be clicked on, with regard to energy
Realize corresponding warm water discharge monitoring.The function of " key operation " and " camera bellows treatment " is truly realized,
The difficulty of thermal discharge waterfrom nuclear power plant satellite monitoring work is greatly reduced, is effectively to manage thermal discharge waterfrom nuclear power plant
Work, for preventing the work such as thermal pollution, protection sea water quality and ecological environment from providing technical guarantee.
The content of the invention
The purpose of the present invention is, in view of the shortcomings of the prior art, provide a kind of solution being defended with thermal infrared at present
During star remote sensing image monitoring thermal discharge waterfrom nuclear power plant situation, complex operation, manual intervention is excessive, lacks quick
Implementation this problem spatial modeling technical support under thermal discharge waterfrom nuclear power plant Satellite Remote Sensing side
Method；
The technical scheme is that：
A kind of thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support；Including following
Step：
Step one：Data prepare；
Red wave band b4 and near infrared band b5 in nuclear power station geographic range is fabricated to single band image；
Step 2：Geometric accurate correction；
By choosing control point, in the range of above-mentioned wave band data fine correction to error permission；
Step 3：Cloud is recognized；
Step 4：Land and water separates；
Step 5：Spatial modeling is carried out with " mono window algorithm "；
Step 6：Temperature anomaly is charted；
According to the result of above-mentioned calculating, the temperature information in studied marine site is classified, is divided temperature rise
9 grades, form nuclear power station water temperature Satellite Remote Sensing Distribution of temperature rise figure.
The step 3 cloud identification includes：
3.1 calculate vegetation index " NDVI "；
Perform spatial modeling sentence：NDVI=(b5-b4)/(b5+b4).
3.2 calculate intermediate variable " Pv "；
Perform spatial modeling sentence：Pv=((b4-0.05)/(0.7-0.05)) * ((b4-0.05)/
(0.7-0.05))。
The step 4 land and water separates and comprises the following steps：
4.1 are separated land and water, obtain intermediate variable " e_water ".
Perform spatial modeling sentence：If Pv=21.0, e_water=0.995；Otherwise e_water=1.0；
Land is divided into building and the class of non-building two by 4.2, obtains non-building " e_surface "；
Perform spatial modeling sentence：If Pv=33.0, e_surface=0.9625+0.0614*e_water
-0.0461*e_water*e_water；Otherwise e_surface=e_water；
4.3 obtain building " e_build ", and " e " value is equal to " e_build " here；
Perform spatial modeling sentence：If e_water=23.0, e_build=0.9589+0.086
*e_surface-0.0671*e_surface*e_surface；Otherwise e_build=e_surface.
The step 5 is comprised the following steps：
1. intermediate variable " c " value is obtained；
Perform spatial modeling sentence：C=e_build* (0.1982007-0.109611*1.45)
2. intermediate variable " Tb " value is calculated；
Thermal infrared bands b10 in nuclear power station geographic range is fabricated to single band image, space is performed and is built
Mould sentence：Tb=1282.71/ (LOG10 ((666.09/ (0.0668235*b10))+1))
3. intermediate variable " D " value is calculated；
Perform spatial modeling sentence：D=(1-0.03926475) * (1+ (1-Tb) * 0.03926475)
4. calculate seawater thermal information and extract result " Ts " value；
Perform spatial modeling sentence：Ts=(- 67.355351* (1-c-D)+(0.458606* (1-c-D)+c+
D)*Tb-D*284.4822306)/Tb；
The beneficial effects of the invention are as follows：
1. contrast artificial point position observation at present and both warm water discharge monitoring methods of aviation flight, general artificial
The method of point position observed temperature, certainty of measurement is high, can accurately reflect thermal discharge waterfrom nuclear power plant to water proximate
Influence.But one-point measurement labor intensive material resources, acquisition range is restricted, and can only gather from
Scattered point data, it is impossible to the spatial variations of warm water discharge influence are showed from face；
2. aviation thermal infrared remote sensing not only has spatial resolution and temperature resolution higher, and can be with
Data acquisition time is adjusted as needed, with conveniently advantage；
3. requirement higher, and somewhat expensive are still needed in itself to airborne platform, easily by weather
The influence of condition；The present invention uses Thermal Infrared Remote Sensing data, real under spatial modeling technical support
The target of fast monitored thermal discharge waterfrom nuclear power plant work is showed, low with monitoring cost, manual intervention is few, energy
The space diffusion of the warm water discharge in reflection wide area and change in time and space situation；With very strong practicality,
The temperature of earth's surface within a short period of time, can be obtained, compared with conventional method, satellite remote sensing has been embodied
The significant advantage of observation.
Specific embodiment
A kind of thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support；Including following
Step：
Step one：Data prepare；
Red wave band b4 and near infrared band b5 in nuclear power station geographic range is fabricated to single band image；
Step 2：Geometric accurate correction；
By choosing control point, in the range of above-mentioned wave band data fine correction to error permission；
Step 3：Cloud is recognized；
Step 4：Land and water separates；
Step 5：Spatial modeling is carried out with " mono window algorithm "；
Step 6：Temperature anomaly is charted；
According to the result of above-mentioned calculating, the temperature information in studied marine site is classified, is divided temperature rise
9 grades, form nuclear power station water temperature Satellite Remote Sensing Distribution of temperature rise figure.
The step 3 cloud identification includes：
3.1 calculate vegetation index " NDVI "；
Perform spatial modeling sentence：NDVI=(b5-b4)/(b5+b4).
3.2 calculate intermediate variable " Pv "；
Perform spatial modeling sentence：Pv=((b4-0.05)/(0.7-0.05)) * ((b4-0.05)/
(0.7-0.05))。
The step 4 land and water separates and comprises the following steps：
4.1 are separated land and water, obtain intermediate variable " e_water ".
Perform spatial modeling sentence：If Pv=21.0, e_water=0.995；Otherwise e_water=1.0；
Land is divided into building and the class of non-building two by 4.2, obtains non-building " e_surface "；
Perform spatial modeling sentence：If Pv=33.0, e_surface=0.9625+0.0614*e_water
-0.0461*e_water*e_water；Otherwise e_surface=e_water；
4.3 obtain building " e_build ", and " e " value is equal to " e_build " here；
Perform spatial modeling sentence：If e_water=23.0, e_build=0.9589+0.086
*e_surface-0.0671*e_surface*e_surface；Otherwise e_build=e_surface.
The step 5 is comprised the following steps：
1. intermediate variable " c " value is obtained；
Perform spatial modeling sentence：C=e_build* (0.1982007-0.109611*1.45)
2. intermediate variable " Tb " value is calculated；
Thermal infrared bands b10 in nuclear power station geographic range is fabricated to single band image, space is performed and is built
Mould sentence：Tb=1282.71/ (LOG10 ((666.09/ (0.0668235*b10))+1))
3. intermediate variable " D " value is calculated；
Perform spatial modeling sentence：D=(1-0.03926475) * (1+ (1-Tb) * 0.03926475)
4. calculate seawater thermal information and extract result " Ts " value；
Perform spatial modeling sentence：Ts=(- 67.355351* (1-c-D)+(0.458606* (1-c-D)+c
+D)*Tb-D*284.4822306)/Tb；
In step one, the present invention uses Landsat-8 satellite datas, and data form is classical TIFF lattice
Formula.The transit time of geographical coordinate and required data according to nuclear power station region, obtains satellite data.
Including 11 wave band image files, 1 quality evaluation file and 1 metadata of TXT forms.From upper
The running environment information of extraction sensor in file is stated, mainly including shooting time, sun altitude, increasing
The information such as beneficial coefficient and longitude and latitude, forms the thermal discharge waterfrom nuclear power plant monitoring basic data of standard.
Geometric accurate correction parameter is set in step 2：Projection is WGS 84, projection zone 50,
Polynomial (polynomial correction), sampling interval 30m.
The method that the identification of step 3 medium cloud is used is that the reflection case according to different atural objects is different, in remote sensing
DN values on image are different, and the DN values of cloud layer are far longer than other atural objects, according to every width shadow
The mask of cloud layer is extracted as the minimum value of upper cloud layer DN values.
By extracting flood boundaries line in step 4, water-outlet body is individually extracted, be easy to that water-outlet body can be made
Temperature change echelon and temperature change distribution map so that the temperature change of water body becomes apparent from.The present invention
Using single band threshold method, the reflection of the other atural objects of luminance factor according to water body near infrared band
Rate is low, and its DN value is smaller, selects six wave bands, determines the maximum reflectivity value in waters, used as difference water
The threshold value of body and other atural objects, extracts the mask of water body.
The inverting of the bright temperature of water surface in step 5, from ripe " mono window algorithm ".Ginseng needed for the algorithm
Number is less, is easier to obtain, good stability.Algorithm had both considered the influence of seawater surface emissivity,
Consider the influence of atmospheric radiation；And the atmospheric parameter needed for refutation process is than traditional conduct radiation side
Journey method much less, it is only necessary near surface temperature and atmospheric water vapor content the two parameters, can be from meteorological portion
Door is easily obtained.
Claims (4)
1. the thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under a kind of spatial modeling technical support；It is characterized in that：Comprise the following steps：
Step one：Data prepare；
Red wave band b4 and near infrared band b5 in nuclear power station geographic range is fabricated to single band image；
Step 2：Geometric accurate correction；
By choosing control point, in the range of above-mentioned wave band data fine correction to error permission；
Step 3：Cloud is recognized；
Step 4：Land and water separates；
Step 5：Spatial modeling is carried out with " mono window algorithm "；
Step 6：Temperature anomaly is charted；
According to the result of above-mentioned calculating, the temperature information in studied marine site is classified, is divided 9 grades of temperature rise, formed nuclear power station water temperature Satellite Remote Sensing Distribution of temperature rise figure.
2. the thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under a kind of spatial modeling technical support as claimed in claim 1, it is characterised in that：The step 3 cloud identification includes：
3.1 calculate vegetation index " NDVI "；
Perform spatial modeling sentence：NDVI=(b5-b4)/(b5+b4).
3.2 calculate intermediate variable " Pv "；
Perform spatial modeling sentence：Pv=((b4-0.05)/(0.7-0.05)) * ((b4-0.05)/(0.7-0.05)).
3. the thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under a kind of spatial modeling technical support as claimed in claim 1, it is characterised in that：The step 4 land and water separates and comprises the following steps：
4.1 are separated land and water, obtain intermediate variable " e_water ".
Perform spatial modeling sentence：If Pv=21.0, e_water=0.995；Otherwise e_water=1.0；
Land is divided into building and the class of non-building two by 4.2, obtains non-building " e_surface "；
Perform spatial modeling sentence：If Pv=33.0, e_surface=0.9625+0.0614*e_water-0.0461*e_water*e_water；Otherwise e_surface=e_water；
4.3 obtain building " e_build ", and " e " value is equal to " e_build " here；
Perform spatial modeling sentence：If e_water=23.0, e_build=0.9589+0.086*e_surface-0.0671*e_surface*e_surfac e；Otherwise e_build=e_surface.
4. the thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under a kind of spatial modeling technical support as claimed in claim 1, it is characterised in that：The step 5 is comprised the following steps：
1. intermediate variable " c " value is obtained；
Perform spatial modeling sentence：C=e_build* (0.1982007-0.109611*1.45)
2. intermediate variable " Tb " value is calculated；
Thermal infrared bands b10 in nuclear power station geographic range is fabricated to single band image, spatial modeling sentence is performed：Tb=1282.71/ (LOG10 ((666.09/ (0.0668235*b10))+1))
3. intermediate variable " D " value is calculated；
Perform spatial modeling sentence：D=(1-0.03926475) * (1+ (1-Tb) * 0.03926475)
4. calculate seawater thermal information and extract result " Ts " value；
Perform spatial modeling sentence：Ts=(- 67.355351* (1-c-D)+(0.458606* (1-c-D)+c+D) * Tb-D*284.4822306)/Tb.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511032319.3A CN106934793A (en) | 2015-12-31 | 2015-12-31 | Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511032319.3A CN106934793A (en) | 2015-12-31 | 2015-12-31 | Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106934793A true CN106934793A (en) | 2017-07-07 |
Family
ID=59444778
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511032319.3A Pending CN106934793A (en) | 2015-12-31 | 2015-12-31 | Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106934793A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111611544A (en) * | 2020-05-12 | 2020-09-01 | 中国科学院上海技术物理研究所 | Thermal imager warm water drainage monitoring method for airborne large-view-field area array swinging |
| CN113610708A (en) * | 2021-07-28 | 2021-11-05 | 国家卫星气象中心(国家空间天气监测预警中心) | Mapping method and device for passive satellite remote sensing flood information |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103203A (en) * | 2011-01-19 | 2011-06-22 | 环境保护部卫星环境应用中心 | Environmental satellite 1-based surface temperature single-window inversion method |
| US20140016667A1 (en) * | 2012-07-13 | 2014-01-16 | Blue Water Satellite, Inc. | Method for monitoring water temperature |
| CN103761704A (en) * | 2014-01-15 | 2014-04-30 | 北京环境特性研究所 | Image generating method and system based on infrared remote sensing data |
-
2015
- 2015-12-31 CN CN201511032319.3A patent/CN106934793A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102103203A (en) * | 2011-01-19 | 2011-06-22 | 环境保护部卫星环境应用中心 | Environmental satellite 1-based surface temperature single-window inversion method |
| US20140016667A1 (en) * | 2012-07-13 | 2014-01-16 | Blue Water Satellite, Inc. | Method for monitoring water temperature |
| CN103761704A (en) * | 2014-01-15 | 2014-04-30 | 北京环境特性研究所 | Image generating method and system based on infrared remote sensing data |
Non-Patent Citations (6)
| Title |
|---|
| 丁凤 等: "TM 热波段图像的地表温度反演算法与实验分析", 《地球信息科学》 * |
| 于杰 等: "利用LandsatTM6数据反演大亚湾海水表层温度", 《国土资源遥感》 * |
| 宋挺 等: "Landsat 8 数据地表温度反演算法对比", 《遥感学报》 * |
| 林卉 王仁礼: "《数字摄影测量学》", 31 August 2015, 中国矿业大学出版社 * |
| 邓非 闫利: "《摄影测量实验教程》", 31 May 2012, 武汉大学出版社 * |
| 马文坡: "《航天光学遥感技术》", 31 January 2011, 中国科学技术出版社 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111611544A (en) * | 2020-05-12 | 2020-09-01 | 中国科学院上海技术物理研究所 | Thermal imager warm water drainage monitoring method for airborne large-view-field area array swinging |
| CN111611544B (en) * | 2020-05-12 | 2024-03-26 | 中国科学院上海技术物理研究所 | Thermal imaging system temperature drainage monitoring method for onboard large-field area array swing scanning |
| CN113610708A (en) * | 2021-07-28 | 2021-11-05 | 国家卫星气象中心(国家空间天气监测预警中心) | Mapping method and device for passive satellite remote sensing flood information |
| CN113610708B (en) * | 2021-07-28 | 2023-11-17 | 国家卫星气象中心(国家空间天气监测预警中心) | Imaging method and device for passive microwave remote sensing flood information |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wegmann et al. | Arctic moisture source for Eurasian snow cover variations in autumn | |
| Chen et al. | Comparisons of the circulation anomalies associated with extreme heat in different regions of eastern China | |
| Gómez-Gesteira et al. | The state of climate in NW Iberia | |
| Cha et al. | A dynamical initialization scheme for real-time forecasts of tropical cyclones using the WRF model | |
| Du et al. | Satellite and Argo observed surface salinity variations in the tropical Indian Ocean and their association with the Indian Ocean dipole mode | |
| Meredith et al. | The freshwater system west of the Antarctic Peninsula: spatial and temporal changes | |
| Xie et al. | Dynamical role of mode water ventilation in decadal variability in the central subtropical gyre of the North Pacific | |
| CN109509319B (en) | Power transmission line forest fire monitoring and early warning method based on static satellite monitoring data | |
| Zhang et al. | Evaluation of cloud effects on air temperature estimation using MODIS LST based on ground measurements over the Tibetan Plateau | |
| CN102253184A (en) | Remote sensing inversion method for land surface evapotranspiration of arid and semi-arid regions | |
| KR20200059085A (en) | A Method for Sea Surface Temperature Retrieval using Surface Drifter Temperature Data and Satellite Infrared Images | |
| Kaldybayev et al. | Glacier change in the Karatal river basin, Zhetysu (Dzhungar) Alatau, Kazakhstan | |
| Kobashi et al. | Interannual variability of the North Pacific Subtropical Countercurrent: Role of local ocean–atmosphere interaction | |
| Yuan et al. | Predecessor rain events over China’s low-latitude highlands associated with Bay of Bengal tropical cyclones | |
| CN106934793A (en) | Thermal discharge waterfrom nuclear power plant satellite remote-sensing monitoring method under spatial modeling technical support | |
| CN108594335B (en) | Global precipitation detection method | |
| Fang | Prediction and analysis of urban heat island effect in Dangshan by remote sensing | |
| Jun et al. | Fog detection over China’s adjacent sea area by using the MTSAT geostationary satellite data | |
| CN108897074B (en) | Global ocean rainfall inversion method | |
| Lin et al. | Recent methane surges reveal heightened emissions from tropical inundated areas | |
| Chu et al. | Rainfall nowcasting by blending of radar data and numerical weather prediction | |
| Yang et al. | FY-3A satellite microwave data assimilation experiments in tropical cyclone forecast | |
| Liang et al. | Climatological Averaging and Application of Optical Vegetation Index based on MODIS | |
| CN103743293B (en) | Reference diagram preparation method utilizing large-scale vegetation region forward-looking infrared guidance | |
| CN118052831A (en) | Sea ice dense area cloud top height satellite inversion method based on cloud shadow information |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170707 |
|
| RJ01 | Rejection of invention patent application after publication |