utm2bl-mtx.xl (Converting UTM Images to Ellipse Latitude/Longitude Coordinates)
Author: Hirohisa Mori / joshua@globalbase.org +
[Transrate]Reiko Inoue Bendtsen c/o Suzaku Translations /
(http://www.suzaku-translations.com/) +
PROTOTYPE
xl utm2bl-mtx.xl / [dest-mtx] [longitude-min] [longitude-max] [latitude-min] [latitude-max] [longitude-resolution] [latitude-resolution] [src-directory] [elipsoid-a] [elipsoid-b]
ARGUMENTS
[dest-mtx] Conversion source matrix file name
[longitude-min] Range of bl created, minimum value of longitude -180degrees to 180degrees
[longitude-max] Range of bl created, maximum value of longitude -180degrees to 540degrees
[latitude-min] Range of bl created, minimum value of latitude -90degrees to 90degrees
[latitude-min] Range of bl created, maximum value of latitude -90degrees to 90degrees
[longitude-resolution] Resolution in the longitude direction, dot/degree
[latitude-resolution] Resolution in the latitude direction, dot/degree
[src-directory] Conversion source directory
[elipsoid-a] Long diameter of ellipse (m)
[elipsoid-b] Short diameter of ellipse (m)
option
-
leveloing=on,off The default is on .
-
part-step=[integer] The default is 1.
-
part-offset=[integer] The default is 1/
ENVIRONMENT
Agent xl [UNDEF REF (xl)]
EXPLANATION
This command transforms matrix data represented in the UTM reference system into data represented in the ellipsoidal latitude/longitude coordinate system.
[longitude-min] and [longitude-max] indicate the range of the image after conversion in the longitude direction. Since it is also possible to convert data to a range including the date line, [longitude-max] can be up to However, [longitude-max]-[longitude-min] cannot exceed Similarly, [latitude-min] and [latitude-max] indicate the range in the latitude direction.
[longitude-resolution] and [latitude-resolution] determines the resolutions in the longitude and latitude directions, respectively. Note that it is not desirable to process high latitude areas and low latitude areas with the same resolution in the longitude direction. The resolution in the longitude direction must be lowered for high latitude areas. For this reason, it is normally desirable to convert a matrix file after dividing it into a southern high latitude area, an area around the equator and a northern high latitude area, and adjusting [latitude-resolution] for each area.
[src-directory] indicates the directory in which conversion source UTM data is saved. Each UTM zone data must be saved in a subdirectory with the zone name. Each subdirectory may contain img.crd, img.map, img.lst and img.mtx files and conform to the UTM coordinate format. Note that the name of the subdirectories must be composed as [zone number] or [zone number]S,[elipsoid-a]s, [elipsoid-a]N, or [elipsoid-a]n. Names without alphabetic characters at the end of a number or those with N or n indicate zones in the northern hemisphere while the names with S or s indicates zones in the southern hemisphere.
[elipsoid-a] and [elipsoid-b] determines the long and short diameters of a speroid. Set these values after checking the spheroid on which UTM is based.
If the option leveling=off is specified, leveling at the last stage is not performed. It is possible to use this option when executing utm2bl-mtx.xl several times and setting layers in a batch at the end.
The options part-step and part-offset are used to acquire UTM files in intervals rather than acquiring all the UTM files. If you specify part-step=2, for instance, every other UTM file is converted. If you specify part-offset=3, conversion is started from UTM3, etc. In general, since distortion of the UTM coordinate system is small at high latitudes, this option can be used for northern latitudes of 60degrees or higher, where even number UTM data is tied to odd number UTM.
Conversion method
The method of conversion when UTM data is given as data saved in multiple world files is explained briefly here.
-
Put together world file data in each UTM zone.
-
Convert world files of each zone to a single matrix data file set using the worldfile-mtx.xl script. At this point, the conversion destination directory must be set to the subdirectory name above.
It is not necessary that the generated UTM data can be viewed. If you use the leveling=off option at this conversion, the last task of data leveling is not performed and the processing time can be shortened.
-
Then, use this script and convert the data to a single data list.
[geosage]$ ls
52 54 53n bl
[geosage]$ xl utm2bl-mtx.xl - - / bl/img0.mtx -180 180 -60 60 3600 3600 . 6378137 6356752
....
[geosage]$ xl utm2bl-mtx.xl - - / bl/img1.mtx -180 180 60 90 1800 3600 . 6378137 6356752 part-step=2 part-offset=1
....
[geosage]$ xl utm2bl-mtx.xl - - / bl/img2.mtx -180 180 -90 -60 1800 3600 . 6378137 6356752 part-step=2 part-offset=1
....
[geosage]$ ls bl
img0.mtx
[geosage]$
-
Finally, create a single crd file, and create a map file and lst file according to img0.mtx to img2.mtx and map all of them to the aforementioned crd file.
[UP]
Go To Page Top
RETURN VALUE
ERRORS
REFERENCE
BUGS
[UP]
Go To Page Top