00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076 #include "tranmerc.h"
00077 #include "utm.h"
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089 #define PI 3.14159265358979323e0
00090 #define MIN_LAT ( (-80.5 * PI) / 180.0 )
00091 #define MAX_LAT ( (84.5 * PI) / 180.0 )
00092 #define MIN_EASTING 100000
00093 #define MAX_EASTING 900000
00094 #define MIN_NORTHING 0
00095 #define MAX_NORTHING 10000000
00096
00097
00098
00099
00100
00101
00102 static double UTM_a = 6378137.0;
00103 static double UTM_f = 1 / 298.257223563;
00104 static long UTM_Override = 0;
00105
00106
00107
00108
00109
00110
00111
00112
00113 long Set_UTM_Parameters(double a,
00114 double f,
00115 long override)
00116 {
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128 double inv_f = 1 / f;
00129 long Error_Code = UTM_NO_ERROR;
00130
00131 if (a <= 0.0)
00132 {
00133 Error_Code |= UTM_A_ERROR;
00134 }
00135 if ((inv_f < 250) || (inv_f > 350))
00136 {
00137 Error_Code |= UTM_INV_F_ERROR;
00138 }
00139 if ((override < 0) || (override > 60))
00140 {
00141 Error_Code |= UTM_ZONE_OVERRIDE_ERROR;
00142 }
00143 if (!Error_Code)
00144 {
00145 UTM_a = a;
00146 UTM_f = f;
00147 UTM_Override = override;
00148 }
00149 return (Error_Code);
00150 }
00151
00152
00153 void Get_UTM_Parameters(double *a,
00154 double *f,
00155 long *override)
00156 {
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166 *a = UTM_a;
00167 *f = UTM_f;
00168 *override = UTM_Override;
00169 }
00170
00171
00172 long Convert_Geodetic_To_UTM (double Latitude,
00173 double Longitude,
00174 long *Zone,
00175 char *Hemisphere,
00176 double *Easting,
00177 double *Northing)
00178 {
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194 long Lat_Degrees;
00195 long Long_Degrees;
00196 long temp_zone;
00197 long Error_Code = UTM_NO_ERROR;
00198 double Origin_Latitude = 0;
00199 double Central_Meridian = 0;
00200 double False_Easting = 500000;
00201 double False_Northing = 0;
00202 double Scale = 0.9996;
00203
00204 if ((Latitude < MIN_LAT) || (Latitude > MAX_LAT))
00205 {
00206 Error_Code |= UTM_LAT_ERROR;
00207 }
00208 if ((Longitude < -PI) || (Longitude > (2*PI)))
00209 {
00210 Error_Code |= UTM_LON_ERROR;
00211 }
00212 if (!Error_Code)
00213 {
00214 if (Longitude < 0)
00215 Longitude += (2*PI);
00216 Lat_Degrees = (long)(Latitude * 180.0 / PI);
00217 Long_Degrees = (long)(Longitude * 180.0 / PI);
00218
00219 if (Longitude < PI)
00220 temp_zone = (long)(31 + ((Longitude * 180.0 / PI) / 6.0));
00221 else
00222 temp_zone = (long)(((Longitude * 180.0 / PI) / 6.0) - 29);
00223 if (temp_zone > 60)
00224 temp_zone = 1;
00225
00226 if ((Lat_Degrees > 55) && (Lat_Degrees < 64) && (Long_Degrees > -1)
00227 && (Long_Degrees < 3))
00228 temp_zone = 31;
00229 if ((Lat_Degrees > 55) && (Lat_Degrees < 64) && (Long_Degrees > 2)
00230 && (Long_Degrees < 12))
00231 temp_zone = 32;
00232 if ((Lat_Degrees > 71) && (Long_Degrees > -1) && (Long_Degrees < 9))
00233 temp_zone = 31;
00234 if ((Lat_Degrees > 71) && (Long_Degrees > 8) && (Long_Degrees < 21))
00235 temp_zone = 33;
00236 if ((Lat_Degrees > 71) && (Long_Degrees > 20) && (Long_Degrees < 33))
00237 temp_zone = 35;
00238 if ((Lat_Degrees > 71) && (Long_Degrees > 32) && (Long_Degrees < 42))
00239 temp_zone = 37;
00240
00241 if (UTM_Override)
00242 {
00243 if ((temp_zone == 1) && (UTM_Override == 60))
00244 temp_zone = UTM_Override;
00245 else if ((temp_zone == 60) && (UTM_Override == 1))
00246 temp_zone = UTM_Override;
00247 else if (((temp_zone-1) <= UTM_Override) && (UTM_Override <= (temp_zone+1)))
00248 temp_zone = UTM_Override;
00249 else
00250 Error_Code = UTM_ZONE_OVERRIDE_ERROR;
00251 }
00252 if (!Error_Code)
00253 {
00254 if (temp_zone >= 31)
00255 Central_Meridian = (6 * temp_zone - 183) * PI / 180.0;
00256 else
00257 Central_Meridian = (6 * temp_zone + 177) * PI / 180.0;
00258 *Zone = temp_zone;
00259 if (Latitude < 0)
00260 {
00261 False_Northing = 10000000;
00262 *Hemisphere = 'S';
00263 }
00264 else
00265 *Hemisphere = 'N';
00266 Set_Transverse_Mercator_Parameters(UTM_a, UTM_f, Origin_Latitude,
00267 Central_Meridian, False_Easting, False_Northing, Scale);
00268 Convert_Geodetic_To_Transverse_Mercator(Latitude, Longitude, Easting,
00269 Northing);
00270 if ((*Easting < MIN_EASTING) || (*Easting > MAX_EASTING))
00271 Error_Code = UTM_EASTING_ERROR;
00272 if ((*Northing < MIN_NORTHING) || (*Northing > MAX_NORTHING))
00273 Error_Code |= UTM_NORTHING_ERROR;
00274 }
00275 }
00276 return (Error_Code);
00277 }
00278
00279
00280 long Convert_UTM_To_Geodetic(long Zone,
00281 char Hemisphere,
00282 double Easting,
00283 double Northing,
00284 double *Latitude,
00285 double *Longitude)
00286 {
00287
00288
00289
00290
00291
00292
00293
00294
00295
00296
00297
00298
00299
00300
00301 long Error_Code = UTM_NO_ERROR;
00302 double Origin_Latitude = 0;
00303 double Central_Meridian = 0;
00304 double False_Easting = 500000;
00305 double False_Northing = 0;
00306 double Scale = 0.9996;
00307
00308 if ((Zone < 1) || (Zone > 60))
00309 Error_Code |= UTM_ZONE_ERROR;
00310 if ((Hemisphere != 'S') && (Hemisphere != 'N'))
00311 Error_Code |= UTM_HEMISPHERE_ERROR;
00312 if ((Easting < MIN_EASTING) || (Easting > MAX_EASTING))
00313 Error_Code |= UTM_EASTING_ERROR;
00314 if ((Northing < MIN_NORTHING) || (Northing > MAX_NORTHING))
00315 Error_Code |= UTM_NORTHING_ERROR;
00316 if (!Error_Code)
00317 {
00318 if (Zone >= 31)
00319 Central_Meridian = ((6 * Zone - 183) * PI / 180.0 + 0.00000005);
00320 else
00321 Central_Meridian = ((6 * Zone + 177) * PI / 180.0 + 0.00000005);
00322 if (Hemisphere == 'S')
00323 False_Northing = 10000000;
00324 Set_Transverse_Mercator_Parameters(UTM_a, UTM_f, Origin_Latitude,
00325 Central_Meridian, False_Easting, False_Northing, Scale);
00326 if (Convert_Transverse_Mercator_To_Geodetic(Easting,
00327 Northing,
00328 Latitude,
00329 Longitude))
00330 Error_Code |= UTM_NORTHING_ERROR;
00331 if ((*Latitude < MIN_LAT) || (*Latitude > MAX_LAT))
00332 {
00333 Error_Code |= UTM_NORTHING_ERROR;
00334 }
00335 }
00336 return (Error_Code);
00337 }