CN107945242A - It is a kind of towards IDL projection transform algorithms - Google Patents
It is a kind of towards IDL projection transform algorithms Download PDFInfo
- Publication number
- CN107945242A CN107945242A CN201711136796.3A CN201711136796A CN107945242A CN 107945242 A CN107945242 A CN 107945242A CN 201711136796 A CN201711136796 A CN 201711136796A CN 107945242 A CN107945242 A CN 107945242A
- Authority
- CN
- China
- Prior art keywords
- latitude
- longitude
- arctic
- global
- projection
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
The present invention relates to one kind towards IDL projection transform algorithms, and the present invention is into polar region normal projection by Arctic data conversion of the Global Temperature data of utm projection.Using the longitude and latitude (hereinafter referred to as global longitude and latitude) of the Global Temperature data of IDL language reading utm projection, the longitude and latitude (hereinafter referred to as arctic longitude and latitude) of Arctic of polar region normal projection, Global Temperature data；Defined variable, for storing transformed arctic temperature record；Calculate the distance of each grid of each grid away from global longitude and latitude of arctic longitude and latitude, assign arctic longitude and latitude index pointed variable grid the pointed temperature record of the corresponding global longitude and latitude index of minimum range, all grids of cycle calculations successively, finally obtain Arctic temperature record of polar region normal projection.The method of the present invention can cut polar region part directly from the data of global utm projection and be converted to polar region normal projection, the former data for solving the problems, such as certain global parameter change hardly possible to polar region towards IDL projection transform algorithms.
Description
Technical field
The present invention relates to one kind towards IDL projection transform algorithms.Belong to Remote Sensing Data Processing technical field.
Background technology
IDL programming languages are powerful, in aerospace, believe etc. that field is widely used to remote sensing.Many researchers exist
When studying polar region associated arguments such as temperature, Hai Wen, wind etc., a problem can be run into, temperature, sea this kind of data of gentle wind are big
All it is the global metadata of utm projection, as shown in figure 4, be Global Temperature distribution map, rather than the regional data of polar region part.
Researcher simply cannot cut arctic regions partial data again, because the projection of global metadata and specializing in polar region
The projection in area is different, this problem brings more difficulty to researcher.But there is presently no directly will be from global utm projection
Data in, the method that directly cuts and be converted to polar region normal projection.In view of the above problems, the present invention implements specific solution final accounts
Method, researcher can download to global longitude and latitude grid data and polar region grid data, based on these data, be using the present invention
It can be cut from the data of global utm projection and switch to the data of polar region normal projection.
The content of the invention
In view of the above-mentioned problems, the present invention proposes that one kind, can be directly from global utm projection towards IDL projection transform algorithms
Polar region part is cut in data and is converted to polar region normal projection.
The technical solution adopted by the present invention to solve the technical problems is：One kind towards IDL projection transform algorithms, including with
Lower step：
Read the global longitude and latitude, arctic longitude and latitude, Global Temperature data of utm projection；Global longitude and latitude is Global Temperature
Longitude and latitude corresponding to data；Arctic longitude and latitude is the longitude and latitude changed into corresponding to polar region normal projection data；
Defined variable, for storing transformed arctic temperature record；
Calculate the distance of each mesh point of each mesh point away from global longitude and latitude of arctic longitude and latitude；Minimum range is corresponded to
The pointed temperature record of global longitude and latitude index assign arctic longitude and latitude index pointed variable grid, finally obtain pole
Arctic temperature record of ground normal projection.
The defined variable type is floating type.
The temperature record that the corresponding global longitude and latitude index of minimum range is pointed assigns arctic longitude and latitude index
Arctic temperature record that pointed variable grid finally obtains polar region normal projection comprises the following steps：
The xth row of arctic temperature grid, the numerical value of the mesh point AT [x-1, y-1] of y rows are obtained with the method for comparing distance
Take：Arctic longitude and latitude data network lattice point (lon [x-1, y-1], lat [x-1, y-1]) is away from global all mesh points of longitude and latitude degrees of data
Distance in, the Global Temperature value at global longitude and latitude mesh point corresponding to minimum range is equal to polar region temperature mesh point AT
The temperature value at [x, y] place；
All data network lattice points of arctic longitude and latitude are traveled through, finally obtain Arctic temperature record of polar region normal projection
The distance is obtained by following formula：
Arctic longitude and latitude data grids xth row, mesh point (lon [x-1, y-1], lat [x-1, the y-1]) distance of y rows are complete
Ball warp latitude data X arranges, the distance L of the mesh point (LON [X-1, Y-1], LAT [X-1, Y-1]) of Y rows is expressed as
1st row, the 1st row are expressed as [0,0], and xth row, y rows are expressed as [x-1, y-1].
One kind is realized towards IDL projection transform algorithms with IDL programs.
The invention has the advantages that and advantage：
1. the method for the present invention can cut polar region portion towards IDL projection transform algorithms directly from the data of global utm projection
Divide and be converted to polar region normal projection, the former data for solving the problems, such as certain global parameter change hardly possible to polar region.
2. an algorithm is succinctly easy to operation.
3. the conversion between a variety of data for projection of an algorithm.
Brief description of the drawings
Fig. 1 is the flow chart of the present invention；
Fig. 2 is grid schematic diagram；
Fig. 3 is temperature grid data schematic diagram；
Fig. 4 is the air temperature distribution figure of global utm projection；
Fig. 5 is the arctic air temperature distribution figure of polar region normal projection；
Embodiment
The present invention will be further described below in conjunction with the accompanying drawings.
As shown in Figure 1, read the longitude and latitude of Global Temperature data of utm projection, polar region normal projection polar region warp
Latitude, Global Temperature data；Defined variable, for storing transformed polar region temperature record；By calculating polar region longitude and latitude
The distance of each each grid of the grid away from global longitude and latitude, by the pointed gas of the corresponding global longitude and latitude index of minimum range
Warm data assign arctic longitude and latitude index pointed variable grid, and all grids of cycle calculations, are finally obtaining polar region just successively
The arctic regions temperature record of axial projection.Example of the present invention is to turn Arctic data of the Global Temperature data of utm projection
Change polar region normal projection into.Comprise the following steps：
Read the global longitude and latitude, arctic longitude and latitude, Global Temperature data of utm projection；Global longitude and latitude is Global Temperature
Longitude and latitude corresponding to data；Arctic longitude and latitude is the longitude and latitude changed into corresponding to polar region normal projection data.
Defined variable, for storing transformed arctic temperature record.It need to define and arctic longitude and latitude sizing grid sum number
The consistent variable of mesh, accurately stores data, defined variable type of the present invention is floating type for the later stage.This example arctic longitude and latitude net
Lattice arrange for 304,448 rows.Therefore it is 304 row, 448 row floating type array to define arctic temperature variables AT, transformed for storing
Arctic temperature record.
Calculate the distance of each grid of each grid away from global longitude and latitude of arctic longitude and latitude.Assuming that global longitude and latitude the
With LON [0,0], (LON [0,0], LON represent global longitude grid to the grid values of a line first row, and [0,0] represents longitude grid
The 1st row the first row mesh point), dimension size is represented with LAT [0,0], i.e. (LON [0,0], LAT [0,0])；The arctic passes through
The grid values of latitude the first row first row represent with lon [0,0], dimension size with lat [0,0], i.e. (lon [0,0], lat
[0,0]).Then grid of the grid of arctic longitude and latitude the first row first row away from global longitude and latitude the first row first row distance L [0,
0] calculation formula is：
Mesh point (lon [0,0], lat [0,0]) and the distance of all mesh points of global longitude are obtained successively.The present invention
Example, global longitude and latitude grid arrange for 192,94 rows；Therefore (lon [0,0], lat [0,0]) and all mesh points of global longitude and latitude
Distance have 18048.
It is signified to assign the pointed temperature record of the corresponding global longitude and latitude index of minimum range to arctic longitude and latitude index
To variable grid.18048 distance values are compared, select the index corresponding to the numerical value of minimum.Because of minimum value
May be multiple, it is therefore desirable to judged.If the index corresponding to minimum value only has 1, by the temperature of this index direction
Data, assign the variable grid AT [0,0] pointed by arctic longitude and latitude index [0,0]；If the index corresponding to minimum value has more
Temperature record that is a, then being directed toward any of which index, the pointed variable grid AT of imparting arctic longitude and latitude index [0,
0], example of the present invention takes first index.
Using circulation, then calculate arctic longitude and latitude mesh point (lon [1,0], lat [1,0]) and all grids of global longitude
The distance of point, seeks calculation principle by AT [0,0], draws AT [1,0]；Through cycle calculations until obtaining AT [303,447], final
To Arctic temperature record AT of polar region normal projection.
One kind is realized towards IDL projection transform algorithms with IDL programs.
Remotely-sensed data is the form of remote sensing image a bit, and researcher can therefrom extract relevant information, some are grid datas,
The data of satellite collection are stored among regular grid in digital form.It is that Zhangbei County pole region positive axis is thrown such as Fig. 2
Shadow grid, black square are that can store data with mesh point, the inside.Temperature grid data example such as Fig. 3, since space limitation is cut
11 row, 31 row temperature record, unit K are taken.When being handled using IDL programs such grid data, be all using each grid as
What unit was calculated.It is the embodiment of algorithm below.
Documentation function, which is read, using IDL reads the longitude and latitude of Global Temperature data of utm projection, polar region normal projection respectively
The longitude and latitude of polar region, Global Temperature data；Longitude and latitude of the global longitude and latitude corresponding to Global Temperature data；Arctic longitude and latitude
It is the longitude and latitude changed into corresponding to polar region normal projection data.
The variable consistent with arctic longitude and latitude sizing grid and number is defined, accurately stores data for the later stage, the present invention is fixed
Adopted types of variables is floating type.This example arctic longitude and latitude grid arranges for 304,448 rows.Therefore defining arctic temperature variables AT is
304 row, 448 row floating type array, for storing transformed arctic temperature record.
Calculate the distance of each grid of each grid away from global longitude and latitude of arctic longitude and latitude.Assuming that global longitude and latitude the
With LON [0,0], (LON [0,0], LON represent global longitude grid to the grid values of a line first row, and [0,0] represents longitude grid
The 1st row the first row mesh point), dimension size is represented with LAT [0,0], i.e. (LON [0,0], LAT [0,0])；The arctic passes through
The grid values of latitude the first row first row represent with lon [0,0], dimension size with lat [0,0], i.e. (lon [0,0], lat
[0,0]).Then grid of the grid of arctic longitude and latitude the first row first row away from global longitude and latitude the first row first row distance L [0,
0] calculation formula is：
Mesh point (lon [0,0], lat [0,0]) and the distance of all mesh points of global longitude are obtained successively.The present invention
Example, global longitude and latitude grid arrange for 192,94 rows；Therefore (lon [0,0], lat [0,0]) and all mesh points of global longitude and latitude
Distance have 18048.
Minimum range numerical value is picked out from 18048 distance values using size discriminant function, minimum range is corresponded to
The pointed temperature record of global longitude and latitude index assign arctic longitude and latitude index pointed variable grid.Because minimum value can
Can be multiple, it is therefore desirable to judged.If the index corresponding to minimum value only has 1, by the temperature number of this index direction
According to the variable grid AT [0,0] pointed by imparting arctic longitude and latitude index [0,0]；If the index corresponding to minimum value have it is multiple,
The temperature record that then any of which index is directed toward, assigns arctic longitude and latitude index pointed variable grid AT [0,0], this
Invention example takes first index.
Circulated using for, first circulate row (448), arranged (304) in circulation, first mesh point and then progress has been calculated
Second mesh point, that is, calculate arctic longitude and latitude mesh point (lon [1,0], lat [1,0]) and all mesh points of global longitude
Distance, by the principle in claim 3 and 4, draws AT [1,0]；Through cycle calculations until obtaining AT [303,447], final
To Arctic temperature record AT of polar region normal projection.As shown in Figure 5.
Claims (5)
- It is 1. a kind of towards IDL projection transform algorithms, it is characterised in that to comprise the following steps：Read the global longitude and latitude, arctic longitude and latitude, Global Temperature data of utm projection；Global longitude and latitude is Global Temperature data Corresponding longitude and latitude；Arctic longitude and latitude is the longitude and latitude changed into corresponding to polar region normal projection data；Defined variable, for storing transformed arctic temperature record；Calculate the distance of each mesh point of each mesh point away from global longitude and latitude of arctic longitude and latitude；Minimum range is corresponding complete The pointed temperature record of ball warp latitude index assigns arctic longitude and latitude index pointed variable grid, is finally obtaining polar region just Arctic temperature record of axial projection.
- It is 2. according to claim 1 a kind of towards IDL projection transform algorithms, it is characterised in that the defined variable type is Floating type.
- It is 3. according to claim 1 a kind of towards IDL projection transform algorithms, it is characterised in that described to correspond to minimum range The pointed temperature record of global longitude and latitude index assign the pointed variable grid of arctic longitude and latitude index and finally obtain pole Arctic temperature record of ground normal projection comprises the following steps：The xth row of arctic temperature grid, the numerical value of the mesh point AT [x-1, y-1] of y rows are obtained with the method for comparing distance：North The distance of pole longitude and latitude data network lattice point (lon [x-1, y-1], lat [x-1, y-1]) away from global all mesh points of longitude and latitude degrees of data In, the Global Temperature value at global longitude and latitude mesh point corresponding to minimum range is equal to polar region temperature mesh point AT [x, y] The temperature value at place；All data network lattice points of arctic longitude and latitude are traveled through, finally obtain Arctic temperature record of polar region normal projection.
- It is 4. according to claim 3 a kind of towards IDL projection transform algorithms, it is characterised in that the distance is obtained by following formula Arrive：Arctic longitude and latitude data grids xth row, mesh point (lon [x-1, y-1], lat [x-1, y-1]) the distance whole world of y rows pass through Latitude data X arranges, the distance L of the mesh point (LON [X-1, Y-1], LAT [X-1, Y-1]) of Y rows is expressed as1st row, the 1st row are expressed as [0,0], and xth row, y rows are expressed as [x-1, y-1].
- It is 5. according to claim 1 a kind of towards IDL projection transform algorithms, it is characterised in that to be realized with IDL programs.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711136796.3A CN107945242B (en) | 2017-11-16 | 2017-11-16 | IDL-oriented projection conversion method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711136796.3A CN107945242B (en) | 2017-11-16 | 2017-11-16 | IDL-oriented projection conversion method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107945242A true CN107945242A (en) | 2018-04-20 |
| CN107945242B CN107945242B (en) | 2021-02-09 |
Family
ID=61931434
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711136796.3A Active CN107945242B (en) | 2017-11-16 | 2017-11-16 | IDL-oriented projection conversion method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107945242B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488392A (en) * | 2019-08-13 | 2019-11-22 | 中国科学院海洋研究所 | A kind of cyclone center's identification and radius evaluation method based on sea-level pressure data |
| CN114417579A (en) * | 2021-12-31 | 2022-04-29 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Different unit grid distance conversion method for troposphere electric wave environment numerical mode result |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080154493A1 (en) * | 2006-12-21 | 2008-06-26 | Thales | Distance estimation method for a moving object having a constrained vertical path profile |
| CN102831626A (en) * | 2012-06-18 | 2012-12-19 | 清华大学 | Visualization method for multivariable spatio-temporal data under polar region projection mode |
| CN104123343A (en) * | 2013-12-26 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Rendering-oriented spatial data real-time coordinate conversion/projection transformation method |
-
2017
- 2017-11-16 CN CN201711136796.3A patent/CN107945242B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080154493A1 (en) * | 2006-12-21 | 2008-06-26 | Thales | Distance estimation method for a moving object having a constrained vertical path profile |
| CN102831626A (en) * | 2012-06-18 | 2012-12-19 | 清华大学 | Visualization method for multivariable spatio-temporal data under polar region projection mode |
| CN104123343A (en) * | 2013-12-26 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Rendering-oriented spatial data real-time coordinate conversion/projection transformation method |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110488392A (en) * | 2019-08-13 | 2019-11-22 | 中国科学院海洋研究所 | A kind of cyclone center's identification and radius evaluation method based on sea-level pressure data |
| CN110488392B (en) * | 2019-08-13 | 2021-05-25 | 中国科学院海洋研究所 | Cyclone center identification and radius estimation method based on sea level air pressure data |
| CN114417579A (en) * | 2021-12-31 | 2022-04-29 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Different unit grid distance conversion method for troposphere electric wave environment numerical mode result |
| CN114417579B (en) * | 2021-12-31 | 2023-01-24 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Different unit grid distance conversion method for troposphere electric wave environment numerical mode result |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107945242B (en) | 2021-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11302060B2 (en) | Method and system for vector-raster overlay analysis of ground surface image area based on edge clipping | |
| CN103838829B (en) | Raster vectorization system based on hierarchical boundary-topology search model | |
| CN110322497A (en) | A kind of interactive point cloud object extraction method based on three-dimensional visualization | |
| CN107945242A (en) | It is a kind of towards IDL projection transform algorithms | |
| CN104715024A (en) | Multimedia hotspot analysis method | |
| CN111222539B (en) | Method for optimizing and expanding supervision classification samples based on multi-source multi-temporal remote sensing image | |
| CN103235810A (en) | Intelligent remote-sensing image control point data search method | |
| CN115965812B (en) | Evaluation method for classification of unmanned aerial vehicle images on wetland vegetation species and land features | |
| CN109829074A (en) | Target area orthography search method, storage medium, electronic equipment | |
| Kang et al. | Support vector machine classification of crop lands using sentinel-2 imagery | |
| Zhao et al. | An entropy and MRF model-based CNN for large-scale landsat image classification | |
| CN103970901A (en) | Geographic information graphic data integration method | |
| CN116109657A (en) | Geographic information data acquisition processing method, system, electronic equipment and storage medium | |
| CN114565833A (en) | Full convolution network building extraction method based on contour guidance and structure attention | |
| CN113762251A (en) | Target classification method and system based on attention mechanism | |
| JP5352435B2 (en) | Classification image creation device | |
| Sitokonstantinou et al. | A sentinel based agriculture monitoring scheme for the control of the cap and food security | |
| CN112668448B (en) | Ecological process change analysis method, device, medium and terminal equipment | |
| CN114485568B (en) | Mapping method and device, computer equipment and storage medium | |
| CN114066910A (en) | Automatic extraction method and system for model effective area of unmanned aerial vehicle photogrammetry | |
| Li et al. | Construction of a small sample dataset and identification of Pitaya trees (Selenicereus) based on UAV image on close-range acquisition | |
| Memduhoğlu | Identifying impervious surfaces for rainwater harvesting feasibility using unmanned aerial vehicle imagery and machine learning classification | |
| Yu et al. | Remote sensing image classification based on dynamic Co-evolutionary parameter optimization of SVM | |
| Pishnamaz Ahmadi et al. | Evaluating the pattern of forest cover changes using fuzzy object-oriented techniques (case study: Kaleybar county) | |
| Lowe et al. | Multi-source K-nearest neighbor, Mean Balanced forest inventory of Georgia |
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 |