SPI_Example/OGR_LatLon.tcl at main · kuohsienchang/SPI_Example · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
#!/bin/sh
# the next line restarts using tclsh \
exec $SPI_PATH/tclsh "$0" "$@"
#============================================================================
# Environnement Canada
# Centre Meteorologique Canadien
# 2100 Trans-Canadienne
# Dorval, Quebec
#
# Projet     : Exemple de scripts.
# Fichier    : OGR_LatLon.tcl
# Creation   : Novembre 2007 - J.P. Gauthier - CMC/CMOE
# Description: Generer un fichier de latlon a un interval specifique (Pour Lewis)
#
# Parametres :
#   <Delta>  : Delta en degre entre les latlon
#
# Retour:
#
# Remarques  :
#
#============================================================================

package require TclData
#package require TclGeoEER
package require Logger

Log::Start [info script] 0.1

set Delta [lindex $argv 0]

#----- Open the output shape file
catch { file delete -force $env(CI_DATA_OUT)/OGR_LatLon.shp }
ogrfile open FILE write $env(CI_DATA_OUT)/OGR_LatLon.shp "ESRI Shapefile"

#----- Create the ouput layer and its field
ogrlayer create FILE LATLON "LatLon"
ogrlayer define LATLON -field Coord Real

ogrgeometry create LINE "Line String"

set nb 0
for { set lon -180.0 } { $lon<180.0 } { set lon [expr $lon+$Delta] } {

   ogrlayer define LATLON -nb [incr nb]
   ogrlayer define LATLON -feature [expr $nb-1] Coord $lon

   ogrgeometry define LINE -points {}
   for { set lat [expr -90+$Delta] } { $lat<[expr 90.0-$Delta] } { set lat [expr $lat+$Delta] } {
      ogrgeometry define LINE -addpoint $lon $lat
   }
   ogrlayer define LATLON -geometry [expr $nb-1] False LINE
}

for { set lat [expr -90+$Delta] } { $lat<[expr 90.0-$Delta] } { set lat [expr $lat+$Delta] } {

   ogrlayer define LATLON -nb [incr nb]

   ogrgeometry define LINE -points {}
   for { set lon -180.0 } { $lon<180.0 } { set lon [expr $lon+$Delta] } {
      ogrgeometry define LINE -addpoint $lon $lat
   }

   ogrlayer define LATLON -feature [expr $nb-1] Coord $lat
   ogrlayer define LATLON -geometry [expr $nb-1] False LINE
}

ogrfile close FILE

Log::End