由于 某些原因 ,国内地图总会有些偏移。虽说这是个相对偏移,但总会对需要坐标定位的应用造成几百米的漂移。若定位坐标在国外,则不会存在此问题。因此我们需要判别当前坐标是否在国内,然后决定是否进行坐标修正。

判定坐标是否在国内,有下列几种方法可选:

  1. 通过网络地图接口查询
  2. 根据国家行政边界判定
  3. 根据国家行政边界近似判定

方法1,通过网络地图接口查询

网络接口:http://maps.google.com/maps/api/geocode/xml?latlng=xxx,xxx&sensor=false,其中latlng参数填写的是坐标的经纬度参数。

例如进行如下调用:http://maps.google.com/maps/api/geocode/xml?latlng=39.829632,117.059326&sensor=false

返回结果如下:

<?xml version="1.0" encoding="UTF-8"?>
<GeocodeResponse><status>OK</status><result><type>route</type><formatted_address>274 Provincial Rd, Sanhe, Langfang, Hebei, China</formatted_address><address_component><long_name>274 Provincial Rd</long_name><short_name>S274</short_name><type>route</type></address_component><address_component><long_name>Sanhe</long_name><short_name>Sanhe</short_name><type>sublocality</type><type>political</type></address_component><address_component><long_name>Langfang</long_name><short_name>Langfang</short_name><type>locality</type><type>political</type></address_component><address_component><long_name>Hebei</long_name><short_name>Hebei</short_name><type>administrative_area_level_1</type><type>political</type></address_component><address_component><long_name>China</long_name><short_name>CN</short_name><type>country</type><type>political</type></address_component><geometry><location><lat>39.8483120</lat><lng>117.0741900</lng></location><location_type>APPROXIMATE</location_type><viewport><southwest><lat>39.8292258</lat><lng>117.0728091</lng></southwest><northeast><lat>39.8674069</lat><lng>117.0757736</lng></northeast></viewport><bounds><southwest><lat>39.8292258</lat><lng>117.0728091</lng></southwest><northeast><lat>39.8674069</lat><lng>117.0757736</lng></northeast></bounds></geometry></result><result><type>sublocality</type><type>political</type><formatted_address>Sanhe, Langfang, Hebei, China</formatted_address><address_component><long_name>Sanhe</long_name><short_name>Sanhe</short_name><type>sublocality</type><type>political</type></address_component><address_component><long_name>Langfang</long_name><short_name>Langfang</short_name><type>locality</type><type>political</type></address_component><address_component><long_name>Hebei</long_name><short_name>Hebei</short_name><type>administrative_area_level_1</type><type>political</type></address_component><address_component><long_name>China</long_name><short_name>CN</short_name><type>country</type><type>political</type></address_component><geometry><location><lat>39.9827180</lat><lng>117.0782950</lng></location><location_type>APPROXIMATE</location_type><viewport><southwest><lat>39.8081545</lat><lng>116.7572022</lng></southwest><northeast><lat>40.0850111</lat><lng>117.2599867</lng></northeast></viewport><bounds><southwest><lat>39.8081545</lat><lng>116.7572022</lng></southwest><northeast><lat>40.0850111</lat><lng>117.2599867</lng></northeast></bounds></geometry></result><result><type>locality</type><type>political</type><formatted_address>Langfang, Hebei, China</formatted_address><address_component><long_name>Langfang</long_name><short_name>Langfang</short_name><type>locality</type><type>political</type></address_component><address_component><long_name>Hebei</long_name><short_name>Hebei</short_name><type>administrative_area_level_1</type><type>political</type></address_component><address_component><long_name>China</long_name><short_name>CN</short_name><type>country</type><type>political</type></address_component><geometry><location><lat>39.5380470</lat><lng>116.6837520</lng></location><location_type>APPROXIMATE</location_type><viewport><southwest><lat>39.4355307</lat><lng>116.5841675</lng></southwest><northeast><lat>39.6029123</lat><lng>116.8444061</lng></northeast></viewport><bounds><southwest><lat>38.4703170</lat><lng>116.1166769</lng></southwest><northeast><lat>40.0850111</lat><lng>117.2599867</lng></northeast></bounds></geometry></result><result><type>administrative_area_level_1</type><type>political</type><formatted_address>Hebei, China</formatted_address><address_component><long_name>Hebei</long_name><short_name>Hebei</short_name><type>administrative_area_level_1</type><type>political</type></address_component><address_component><long_name>China</long_name><short_name>CN</short_name><type>country</type><type>political</type></address_component><geometry><location><lat>38.0370570</lat><lng>114.4686650</lng></location><location_type>APPROXIMATE</location_type><viewport><southwest><lat>36.0482067</lat><lng>113.4653687</lng></southwest><northeast><lat>42.6197178</lat><lng>119.8542560</lng></northeast></viewport><bounds><southwest><lat>36.0482067</lat><lng>113.4653687</lng></southwest><northeast><lat>42.6197178</lat><lng>119.8542560</lng></northeast></bounds></geometry></result><result><type>country</type><type>political</type><formatted_address>China</formatted_address><address_component><long_name>China</long_name><short_name>CN</short_name><type>country</type><type>political</type></address_component><geometry><location><lat>35.8616600</lat><lng>104.1953970</lng></location><location_type>APPROXIMATE</location_type><viewport><southwest><lat>18.1535216</lat><lng>73.4994137</lng></southwest><northeast><lat>53.5609740</lat><lng>134.7728100</lng></northeast></viewport><bounds><southwest><lat>18.1535216</lat><lng>73.4994137</lng></southwest><northeast><lat>53.5609740</lat><lng>134.7728100</lng></northeast></bounds></geometry></result>
</GeocodeResponse>

查询结果各字段的含义可以参考Google Map API的Reverse Geocoding (Address Lookup)章节。简化后可以直接判断返回的字符串中是否存在“<short_name>CN</short_name>”作为判定所查询坐标是否在中国的依据。

方法2,根据国家行政边界判定

由于国家行政边界是由多个简单多边形组成的。对多个多边形进行命中测试 (Hit Testing),只要其中一个多边形测试成功,则表示整个命中测试成功。但多边形的命中测试有些麻烦,更简单的方法是判别点是否在多边形内(Point in polygon(PIP))。这些算法中简单高效的莫过于是光线投射算法 (Ray casting algorithm)。

光线投射算法的要义是:在二维空间的任意一点射出一条射线,射线穿过多边形边界奇数次则这个点在多边形内,否则在多边形外。此方法适用于复杂多边形。实现此算法往往会限制射线的方向是向右的,这样的好处是简化了许多计算。此算法的一个实现参考1987年Paul Bourke教授所写“Determining if a point lies on the interior of a polygon”的Solution 1。

在 Windows Phone 7 中的实现,首先需要获取到国家行政边界的WGS84坐标数据。Alexey Strakh(链接有恶意软件警告)所开发的gMaps应用中包含了此边界数据。因此,自己可以把数据抠出来。抠出的数据由四部分组成我国的行政边界:大陆、台湾省、海南省、崇明岛。此数据相对来说比较准确,但是边界附近误差较大。每次判别坐标是否在此区域内的计算次数较多,大约需要几百次的几何计算。

此方法的实现代码如下:

#region polygonList<Point[]> polygon = new List<Point[]>();private void InitPolygon()
{polygon.Add(new Point[]  //Mainland{new Point(27.32083, 88.91693),new Point(27.54243, 88.76464),new Point(28.00805, 88.83575),new Point(28.1168, 88.62435),new Point(27.86605, 88.14279),new Point(27.82305, 87.19275),new Point(28.11166, 86.69527),new Point(27.90888, 86.45137),new Point(28.15805, 86.19769),new Point(27.88625, 86.0054),new Point(28.27916, 85.72137),new Point(28.30666, 85.11095),new Point(28.59104, 85.19518),new Point(28.54444, 84.84665),new Point(28.73402, 84.48623),new Point(29.26097, 84.11651),new Point(29.18902, 83.5479),new Point(29.63166, 83.19109),new Point(30.06923, 82.17525),new Point(30.33444, 82.11123),new Point(30.385, 81.42623),new Point(30.01194, 81.23221),new Point(30.20435, 81.02536),new Point(30.57552, 80.207),new Point(30.73374, 80.25423),new Point(30.96583, 79.86304),new Point(30.95708, 79.55429),new Point(31.43729, 79.08082),new Point(31.30895, 78.76825),new Point(31.96847, 78.77075),new Point(32.24304, 78.47594),new Point(32.5561, 78.40595),new Point(32.63902, 78.74623),new Point(32.35083, 78.9711),new Point(32.75666, 79.52874),new Point(33.09944, 79.37511),new Point(33.42863, 78.93623),new Point(33.52041, 78.81387),new Point(34.06833, 78.73581),new Point(34.35001, 78.98535),new Point(34.6118, 78.33707),new Point(35.28069, 78.02305),new Point(35.49902, 78.0718),new Point(35.50133, 77.82393),new Point(35.6125, 76.89526),new Point(35.90665, 76.55304),new Point(35.81458, 76.18061),new Point(36.07082, 75.92887),new Point(36.23751, 76.04166),new Point(36.66343, 75.85984),new Point(36.73169, 75.45179),new Point(36.91156, 75.39902),new Point(36.99719, 75.14787),new Point(37.02782, 74.56543),new Point(37.17, 74.39089),new Point(37.23733, 74.91574),new Point(37.40659, 75.18748),new Point(37.65243, 74.9036),new Point(38.47256, 74.85442),new Point(38.67438, 74.35471),new Point(38.61271, 73.81401),new Point(38.88653, 73.70818),new Point(38.97256, 73.85235),new Point(39.23569, 73.62005),new Point(39.45483, 73.65569),new Point(39.59965, 73.95471),new Point(39.76896, 73.8429),new Point(40.04202, 73.99096),new Point(40.32792, 74.88089),new Point(40.51723, 74.8588),new Point(40.45042, 75.23394),new Point(40.64452, 75.58284),new Point(40.298, 75.70374),new Point(40.35324, 76.3344),new Point(41.01258, 76.87067),new Point(41.04079, 78.08083),new Point(41.39286, 78.39554),new Point(42.03954, 80.24513),new Point(42.19622, 80.23402),new Point(42.63245, 80.15804),new Point(42.81565, 80.25796),new Point(42.88545, 80.57226),new Point(43.02906, 80.38405),new Point(43.1683, 80.81526),new Point(44.11378, 80.36887),new Point(44.6358, 80.38499),new Point(44.73408, 80.51589),new Point(44.90282, 79.87106),new Point(45.3497, 81.67928),new Point(45.15748, 81.94803),new Point(45.13303, 82.56638),new Point(45.43581, 82.64624),new Point(45.5831, 82.32179),new Point(47.20061, 83.03443),new Point(46.97332, 83.93026),new Point(46.99361, 84.67804),new Point(46.8277, 84.80318),new Point(47.0591, 85.52257),new Point(47.26221, 85.70139),new Point(47.93721, 85.53707),new Point(48.39333, 85.76596),new Point(48.54277, 86.59791),new Point(49.1102, 86.87602),new Point(49.09262, 87.34821),new Point(49.17295, 87.8407),new Point(48.98304, 87.89291),new Point(48.88103, 87.7611),new Point(48.73499, 88.05942),new Point(48.56541, 87.99194),new Point(48.40582, 88.51679),new Point(48.21193, 88.61179),new Point(47.99374, 89.08514),new Point(47.88791, 90.07096),new Point(46.95221, 90.9136),new Point(46.57735, 91.07027),new Point(46.29694, 90.92151),new Point(46.01735, 91.02651),new Point(45.57972, 90.68193),new Point(45.25305, 90.89694),new Point(45.07729, 91.56088),new Point(44.95721, 93.5547),new Point(44.35499, 94.71735),new Point(44.29416, 95.41061),new Point(44.01937, 95.34109),new Point(43.99311, 95.53339),new Point(43.28388, 95.87901),new Point(42.73499, 96.38206),new Point(42.79583, 97.1654),new Point(42.57194, 99.51012),new Point(42.67707, 100.8425),new Point(42.50972, 101.8147),new Point(42.23333, 102.0772),new Point(41.88721, 103.4164),new Point(41.87721, 104.5267),new Point(41.67068, 104.5237),new Point(41.58666, 105.0065),new Point(42.46624, 107.4758),new Point(42.42999, 109.3107),new Point(42.64576, 110.1064),new Point(43.31694, 110.9897),new Point(43.69221, 111.9583),new Point(44.37527, 111.4214),new Point(45.04944, 111.873),new Point(45.08055, 112.4272),new Point(44.8461, 112.853),new Point(44.74527, 113.638),new Point(45.38943, 114.5453),new Point(45.4586, 115.7019),new Point(45.72193, 116.2104),new Point(46.29583, 116.5855),new Point(46.41888, 117.3755),new Point(46.57069, 117.425),new Point(46.53645, 117.8455),new Point(46.73638, 118.3147),new Point(46.59895, 119.7068),new Point(46.71513, 119.9315),new Point(46.90221, 119.9225),new Point(47.66499, 119.125),new Point(47.99475, 118.5393),new Point(48.01125, 117.8046),new Point(47.65741, 117.3827),new Point(47.88805, 116.8747),new Point(47.87819, 116.2624),new Point(47.69186, 115.9231),new Point(47.91749, 115.5944),new Point(48.14353, 115.5491),new Point(48.25249, 115.8358),new Point(48.52055, 115.8111),new Point(49.83047, 116.7114),new Point(49.52058, 117.8747),new Point(49.92263, 118.5746),new Point(50.09631, 119.321),new Point(50.33028, 119.36),new Point(50.39027, 119.1386),new Point(51.62083, 120.0641),new Point(52.115, 120.7767),new Point(52.34423, 120.6259),new Point(52.54267, 120.7122),new Point(52.58805, 120.0819),new Point(52.76819, 120.0314),new Point(53.26374, 120.8307),new Point(53.54361, 123.6147),new Point(53.18832, 124.4933),new Point(53.05027, 125.62),new Point(52.8752, 125.6573),new Point(52.75722, 126.0968),new Point(52.5761, 125.9943),new Point(52.12694, 126.555),new Point(51.99437, 126.4412),new Point(51.38138, 126.9139),new Point(51.26555, 126.8176),new Point(51.31923, 126.9689),new Point(51.05825, 126.9331),new Point(50.74138, 127.2919),new Point(50.31472, 127.334),new Point(50.20856, 127.5861),new Point(49.80588, 127.515),new Point(49.58665, 127.838),new Point(49.58443, 128.7119),new Point(49.34676, 129.1118),new Point(49.4158, 129.4902),new Point(48.86464, 130.2246),new Point(48.86041, 130.674),new Point(48.60576, 130.5236),new Point(48.3268, 130.824),new Point(48.10839, 130.6598),new Point(47.68721, 130.9922),new Point(47.71027, 132.5211),new Point(48.09888, 133.0827),new Point(48.06888, 133.4843),new Point(48.39112, 134.4153),new Point(48.26713, 134.7408),new Point(47.99207, 134.5576),new Point(47.70027, 134.7608),new Point(47.32333, 134.1825),new Point(46.64017, 133.9977),new Point(46.47888, 133.8472),new Point(46.25363, 133.9016),new Point(45.82347, 133.4761),new Point(45.62458, 133.4702),new Point(45.45083, 133.1491),new Point(45.05694, 133.0253),new Point(45.34582, 131.8684),new Point(44.97388, 131.4691),new Point(44.83649, 130.953),new Point(44.05193, 131.298),new Point(43.53624, 131.1912),new Point(43.38958, 131.3104),new Point(42.91645, 131.1285),new Point(42.74485, 130.4327),new Point(42.42186, 130.6044),new Point(42.71416, 130.2468),new Point(42.88794, 130.2514),new Point(43.00457, 129.9046),new Point(42.43582, 129.6955),new Point(42.44624, 129.3493),new Point(42.02736, 128.9269),new Point(42.00124, 128.0566),new Point(41.58284, 128.3002),new Point(41.38124, 128.1529),new Point(41.47249, 127.2708),new Point(41.79222, 126.9047),new Point(41.61176, 126.5661),new Point(40.89694, 126.0118),new Point(40.47037, 124.8851),new Point(40.09362, 124.3736),new Point(39.82777, 124.128),new Point(39.8143, 123.2422),new Point(39.67388, 123.2167),new Point(38.99638, 121.648),new Point(38.8611, 121.6982),new Point(38.71909, 121.1873),new Point(38.91221, 121.0887),new Point(39.09013, 121.6794),new Point(39.2186, 121.5994),new Point(39.35166, 121.7511),new Point(39.52847, 121.2283),new Point(39.62322, 121.533),new Point(39.81138, 121.4683),new Point(40.00305, 121.881),new Point(40.50562, 122.2987),new Point(40.73874, 122.0521),new Point(40.92194, 121.1775),new Point(40.1961, 120.4468),new Point(39.87242, 119.5264),new Point(39.15693, 118.9715),new Point(39.04083, 118.3273),new Point(39.19846, 117.889),new Point(38.67555, 117.5364),new Point(38.38666, 117.6722),new Point(38.16721, 118.0281),new Point(38.1529, 118.8378),new Point(37.87832, 119.0355),new Point(37.30054, 118.9566),new Point(37.14361, 119.2328),new Point(37.15138, 119.7672),new Point(37.35228, 119.8529),new Point(37.83499, 120.7371),new Point(37.42458, 121.58),new Point(37.55256, 122.1282),new Point(37.41833, 122.1814),new Point(37.39624, 122.5586),new Point(37.20999, 122.5972),new Point(37.02583, 122.4005),new Point(37.01978, 122.5392),new Point(36.89361, 122.5047),new Point(36.84298, 122.1923),new Point(37.00027, 121.9566),new Point(36.75889, 121.5944),new Point(36.61666, 120.7764),new Point(36.52638, 120.96),new Point(36.37582, 120.8753),new Point(36.42277, 120.7062),new Point(36.14075, 120.6956),new Point(36.0419, 120.3436),new Point(36.26345, 120.3078),new Point(36.19998, 120.0889),new Point(35.95943, 120.2378),new Point(35.57893, 119.6475),new Point(34.88499, 119.1761),new Point(34.31145, 120.2487),new Point(32.97499, 120.8858),new Point(32.63889, 120.8375),new Point(32.42958, 121.3348),new Point(32.11333, 121.4412),new Point(32.02166, 121.7066),new Point(31.67833, 121.8275),new Point(31.86639, 120.9444),new Point(32.09361, 120.6019),new Point(31.94555, 120.099),new Point(32.30638, 119.8267),new Point(32.26277, 119.6317),new Point(31.90388, 120.1364),new Point(31.98833, 120.7026),new Point(31.81944, 120.7196),new Point(31.30889, 121.6681),new Point(30.97986, 121.8828),new Point(30.85305, 121.8469),new Point(30.56889, 120.9915),new Point(30.33555, 120.8144),new Point(30.39298, 120.4586),new Point(30.19694, 120.15),new Point(30.31027, 120.5082),new Point(30.06465, 120.7916),new Point(30.30458, 121.2808),new Point(29.96305, 121.6778),new Point(29.88211, 122.1196),new Point(29.51167, 121.4483),new Point(29.58916, 121.9744),new Point(29.19527, 121.9336),new Point(29.18388, 121.8119),new Point(29.37236, 121.7969),new Point(29.19729, 121.7444),new Point(29.29111, 121.5611),new Point(29.1634, 121.4135),new Point(29.02194, 121.6914),new Point(28.9359, 121.4908),new Point(28.72798, 121.6113),new Point(28.84215, 121.1464),new Point(28.66993, 121.4844),new Point(28.34722, 121.6417),new Point(28.13889, 121.3419),new Point(28.38277, 121.1651),new Point(27.98222, 120.9353),new Point(28.07944, 120.5908),new Point(27.87229, 120.84),new Point(27.59319, 120.5812),new Point(27.45083, 120.6655),new Point(27.20777, 120.5075),new Point(27.28278, 120.1896),new Point(27.14764, 120.4211),new Point(26.89805, 120.0332),new Point(26.64465, 120.128),new Point(26.51778, 119.8603),new Point(26.78823, 120.0733),new Point(26.64888, 119.8668),new Point(26.79611, 119.7879),new Point(26.75625, 119.5503),new Point(26.44222, 119.8204),new Point(26.47388, 119.5775),new Point(26.33861, 119.658),new Point(26.36777, 119.9489),new Point(25.99694, 119.4253),new Point(26.14041, 119.0975),new Point(25.93788, 119.354),new Point(25.99069, 119.7058),new Point(25.67996, 119.5807),new Point(25.68222, 119.4522),new Point(25.35333, 119.6454),new Point(25.60649, 119.3149),new Point(25.42097, 119.1053),new Point(25.25319, 119.3526),new Point(25.17208, 119.2726),new Point(25.2426, 118.8749),new Point(24.97194, 118.9866),new Point(24.88291, 118.5729),new Point(24.75673, 118.7631),new Point(24.52861, 118.5953),new Point(24.53638, 118.2397),new Point(24.68194, 118.1688),new Point(24.44024, 118.0199),new Point(24.46019, 117.7947),new Point(24.25875, 118.1237),new Point(23.62437, 117.1957),new Point(23.65919, 116.9179),new Point(23.355, 116.7603),new Point(23.42024, 116.5322),new Point(23.23666, 116.7871),new Point(23.21083, 116.5139),new Point(22.93902, 116.4817),new Point(22.73916, 115.7978),new Point(22.88416, 115.6403),new Point(22.65889, 115.5367),new Point(22.80833, 115.1614),new Point(22.70277, 114.8889),new Point(22.53305, 114.8722),new Point(22.64027, 114.718),new Point(22.81402, 114.7782),new Point(22.69972, 114.5208),new Point(22.50423, 114.6136),new Point(22.55004, 114.2223),new Point(22.42993, 114.3885),new Point(22.26056, 114.2961),new Point(22.36736, 113.9056),new Point(22.50874, 114.0337),new Point(22.47444, 113.8608),new Point(22.83458, 113.606),new Point(23.05027, 113.5253),new Point(23.11724, 113.8219),new Point(23.05083, 113.4793),new Point(22.87986, 113.3629),new Point(22.54944, 113.5648),new Point(22.18701, 113.5527),new Point(22.56701, 113.1687),new Point(22.17965, 113.3868),new Point(22.04069, 113.2226),new Point(22.20485, 113.0848),new Point(21.8693, 112.94),new Point(21.96472, 112.824),new Point(21.70139, 112.2819),new Point(21.91611, 111.8921),new Point(21.75139, 111.9669),new Point(21.77819, 111.6762),new Point(21.61264, 111.7832),new Point(21.5268, 111.644),new Point(21.52528, 111.0285),new Point(21.21138, 110.5328),new Point(21.37322, 110.3944),new Point(20.84381, 110.1594),new Point(20.84083, 110.3755),new Point(20.64, 110.3239),new Point(20.48618, 110.5274),new Point(20.24611, 110.2789),new Point(20.2336, 109.9244),new Point(20.4318, 110.0069),new Point(20.92416, 109.6629),new Point(21.44694, 109.9411),new Point(21.50569, 109.6605),new Point(21.72333, 109.5733),new Point(21.49499, 109.5344),new Point(21.39666, 109.1428),new Point(21.58305, 109.1375),new Point(21.61611, 108.911),new Point(21.79889, 108.8702),new Point(21.59888, 108.7403),new Point(21.93562, 108.4692),new Point(21.59014, 108.5125),new Point(21.68999, 108.3336),new Point(21.51444, 108.2447),new Point(21.54241, 107.99),new Point(21.66694, 107.7831),new Point(21.60526, 107.3627),new Point(22.03083, 106.6933),new Point(22.45682, 106.5517),new Point(22.76389, 106.7875),new Point(22.86694, 106.7029),new Point(22.91253, 105.8771),new Point(23.32416, 105.3587),new Point(23.18027, 104.9075),new Point(22.81805, 104.7319),new Point(22.6875, 104.3747),new Point(22.79812, 104.1113),new Point(22.50387, 103.9687),new Point(22.78287, 103.6538),new Point(22.58436, 103.5224),new Point(22.79451, 103.3337),new Point(22.43652, 103.0304),new Point(22.77187, 102.4744),new Point(22.39629, 102.1407),new Point(22.49777, 101.7415),new Point(22.20916, 101.5744),new Point(21.83444, 101.7653),new Point(21.14451, 101.786),new Point(21.17687, 101.2919),new Point(21.57264, 101.1482),new Point(21.76903, 101.099),new Point(21.47694, 100.6397),new Point(21.43546, 100.2057),new Point(21.72555, 99.97763),new Point(22.05018, 99.95741),new Point(22.15592, 99.16785),new Point(22.93659, 99.56484),new Point(23.08204, 99.5113),new Point(23.18916, 98.92747),new Point(23.97076, 98.67991),new Point(24.16007, 98.89073),new Point(23.92999, 97.54762),new Point(24.26055, 97.7593),new Point(24.47666, 97.54305),new Point(24.73992, 97.55255),new Point(25.61527, 98.19109),new Point(25.56944, 98.36137),new Point(25.85597, 98.7104),new Point(26.12527, 98.56944),new Point(26.18472, 98.73109),new Point(26.79166, 98.77777),new Point(27.52972, 98.69699),new Point(27.6725, 98.45888),new Point(27.54014, 98.31992),new Point(28.14889, 98.14499),new Point(28.54652, 97.55887),new Point(28.22277, 97.34888),new Point(28.46749, 96.65387),new Point(28.35111, 96.40193),new Point(28.525, 96.34027),new Point(28.79569, 96.61373),new Point(29.05666, 96.47083),new Point(28.90138, 96.17532),new Point(29.05972, 96.14888),new Point(29.25757, 96.39172),new Point(29.46444, 96.08315),new Point(29.03527, 95.38777),new Point(29.33346, 94.64751),new Point(29.07348, 94.23456),new Point(28.6692, 93.96172),new Point(28.61876, 93.35194),new Point(28.3193, 93.22205),new Point(28.1419, 92.71044),new Point(27.86194, 92.54498),new Point(27.76472, 91.65776),new Point(27.945, 91.66277),new Point(28.08111, 91.30138),new Point(27.96999, 91.08693),new Point(28.07958, 90.3765),new Point(28.24257, 90.38898),new Point(28.32369, 89.99819),new Point(28.05777, 89.48749),new Point(27.32083, 88.91693)});polygon.Add(new Point[]  //Taiwan{new Point(25.13474, 121.4441),new Point(25.28361, 121.5632),new Point(25.00722, 122.0004),new Point(24.85028, 121.8182),new Point(24.47638, 121.8397),new Point(23.0875, 121.3556),new Point(21.92791, 120.7196),new Point(22.31277, 120.6103),new Point(22.54044, 120.3071),new Point(23.04437, 120.0539),new Point(23.61708, 120.1112),new Point(25.00166, 121.0017),new Point(25.13474, 121.4441)});polygon.Add(new Point[]  //Hainan{new Point(19.52888, 110.855),new Point(19.16761, 110.4832),new Point(18.80083, 110.5255),new Point(18.3852, 110.0503),new Point(18.39152, 109.7594),new Point(18.19777, 109.7036),new Point(18.50562, 108.6871),new Point(19.28028, 108.6283),new Point(19.76, 109.2939),new Point(19.7236, 109.1653),new Point(19.89972, 109.2572),new Point(19.82861, 109.4658),new Point(19.99389, 109.6108),new Point(20.13361, 110.6655),new Point(19.97861, 110.9425),new Point(19.63829, 111.0215),new Point(19.52888, 110.855)});polygon.Add(new Point[]  //Chongming{new Point(31.80054, 121.2039),new Point(31.49972, 121.8736),new Point(31.53111, 121.5464),new Point(31.80054, 121.2039)});
}
private bool IsInsideChina(GeoCoordinate positionToCheck)
{if (positionToCheck == null || positionToCheck.IsUnknown)return false;return polygon.Any(s => IsInsidePolygon(s, positionToCheck.Latitude, positionToCheck.Longitude));
}
private bool IsInsidePolygon(Point[] poly, double x, double y)
{var index = 0;var inside = false;var prePoint = poly[0];for (index = 1; index < poly.Length; index++){var nextPoint = poly[index];if (y > Math.Min(prePoint.Y, nextPoint.Y)&& y <= Math.Max(prePoint.Y, nextPoint.Y)&& x <= Math.Max(prePoint.X, nextPoint.X)&& prePoint.Y != nextPoint.Y){double xinters = (y - prePoint.Y) * (nextPoint.X - prePoint.X) / (nextPoint.Y - prePoint.Y) + prePoint.X;if (prePoint.X == nextPoint.X || x <= xinters)inside ^= true;}prePoint = nextPoint;}return inside;
}private void MainPage_Loaded(object sender, RoutedEventArgs e)
{InitPolygon();foreach (var points in polygon){var poly = new MapPolyline();poly.Locations = new LocationCollection();for (int i = 0; i < points.Length; i++){poly.Locations.Add(new GeoCoordinate(points[i].X, points[i].Y));}poly.Stroke = new SolidColorBrush(Colors.Red);myMap.Children.Add(poly);}myMap.SetView(new GeoCoordinate(39.4916649164344, 108.745576123357), 3.7090892722116307);
}#endregion

方法3,根据国家行政边界近似判定

有时候需要一种快速的方法判定,这种快速判定的方法是从Nokia Maps中挖掘出来的。基本思路是:把整个行政区域划分为几个小的矩形,然后再排除掉一些矩形区域。只要一个点在限定的区域内,并且不在排除的区域内,则判定成功,否则失败。下图中蓝色区域为限定区域,红色区域为排除区域。

此方法的实现代码如下:

#region rectangleprivate class Rectangle
{public double West;public double North;public double East;public double South;public Rectangle(double latitude1, double longitude1, double latitude2, double longitude2){this.West = Math.Min(longitude1, longitude2);this.North = Math.Max(latitude1, latitude2);this.East = Math.Max(longitude1, longitude2);this.South = Math.Min(latitude1, latitude2);}
}private static Rectangle[] region = new Rectangle[]
{new Rectangle(49.220400, 079.446200, 42.889900, 096.330000),new Rectangle(54.141500, 109.687200, 39.374200, 135.000200),new Rectangle(42.889900, 073.124600, 29.529700, 124.143255),new Rectangle(29.529700, 082.968400, 26.718600, 097.035200),new Rectangle(29.529700, 097.025300, 20.414096, 124.367395),new Rectangle(20.414096, 107.975793, 17.871542, 111.744104),};
private static Rectangle[] exclude = new Rectangle[]
{new Rectangle(25.398623, 119.921265, 21.785006, 122.497559),new Rectangle(22.284000, 101.865200, 20.098800, 106.665000),new Rectangle(21.542200, 106.452500, 20.487800, 108.051000),new Rectangle(55.817500, 109.032300, 50.325700, 119.127000),new Rectangle(55.817500, 127.456800, 49.557400, 137.022700),new Rectangle(44.892200, 131.266200, 42.569200, 137.022700),
};public static bool IsInsideChina(GeoCoordinate pos)
{for (int i = 0; i < region.Length; i++){if (InRectangle(region[i], pos)){for (int j = 0; j < exclude.Length; j++){if (InRectangle(exclude[j], pos)){return false;}}return true;}}return false;
}private static bool InRectangle(Rectangle rect, GeoCoordinate pos)
{return rect.West <= pos.Longitude && rect.East >= pos.Longitude && rect.North >= pos.Latitude && rect.South <= pos.Latitude;
}private void MainPage_Loaded(object sender, RoutedEventArgs e)
{DrawRect(region, new SolidColorBrush(Colors.Blue));DrawRect(exclude, new SolidColorBrush(Colors.Red));myMap.SetView(new GeoCoordinate(39.4916649164344, 108.745576123357), 3.7090892722116307);
}private void DrawRect(Rectangle[] rects, SolidColorBrush brush)
{foreach (var rect in rects){var poly = new MapPolyline();poly.Locations = new LocationCollection();poly.Locations.Add(new GeoCoordinate(rect.North, rect.West));poly.Locations.Add(new GeoCoordinate(rect.North, rect.East));poly.Locations.Add(new GeoCoordinate(rect.South, rect.East));poly.Locations.Add(new GeoCoordinate(rect.South, rect.West));poly.Locations.Add(new GeoCoordinate(rect.North, rect.West));poly.Stroke = brush;myMap.Children.Add(poly);}
}#endregion

这几种方法的比较

  备注
方法1 
通过网络地图接口查询
简单、准确 
信息量大
速度慢 
依赖网络
可用Google、Baidu地图接口
方法2 
根据国家行政边界判定
准确 代码相对较多 
判别计算次数较多 
陆地边界数据误差较大
方法和数据从gMaps中挖掘
方法3 
根据国家行政边界近似判定
快速 
相对准确
边界误差较大 方法和数据从Nokia Maps中挖掘

PS:你可能会问,在方法3中为什么把台湾给排除掉了,这是因为所谓的坐标偏移对台湾地区不适用。对于香港澳门地区,视情况需要单独判定。

判断GPS坐标是否在中国相关推荐

  1. (iOS)判断GPS坐标是否在中国

    由于火星坐标问题,所以需要判断一下经纬度是否在中国. 基本思路是:把整个行政区域划分为几个小的矩形,然后再排除掉一些矩形区域.只要一个点在限定的区域内,并且不在排除的区域内,则判定成功,否则失败.下图 ...

  2. Java判断gps点是否在中国,如何判断一个指定的位置点坐标(GPS上的经纬度)是否落在一个多边形区域内?...

    业务场景举例:快递选择收获区域.车辆电子围栏.运动轨迹路线.地理位置信息检测范围和地图等过滤等等. 比方说地图上有一块区域(抽象成多边形),然后里面每一个位置点(像素点)都有对应的GPS的经纬度坐标值 ...

  3. 实现高德坐标转GPS坐标

    vue 实现高德坐标转GPS坐标 首先介绍一下常见的几种地图的坐标类型: WGS-84:这是一个国际标准,也就是GPS坐标(Google Earth.或者GPS模块采集的都是这个类型). GCJ-02 ...

  4. 高德坐标转GPS坐标 CV版

    WGS-84:是国际标准,GPS坐标(Google Earth使用.或者GPS模块) GCJ-02:中国坐标偏移标准,Google Map.高德.腾讯使用 BD-09:百度坐标偏移标准,Baidu M ...

  5. 百度地图批量转换 GPS坐标转百度地图坐标 问题

    百度地图的官方网址 官方批量转换的demo 花了几天时间了解了一下百度地图,之前是后端的一个小伙伴在负责,他跟我吐槽这是前端的东西,让我来写(之前他们老大交给他了,我也以为是后端的任务(๑′ᴗ‵๑). ...

  6. GPS坐标对应地图坐标偏移问题

    各个地图的坐标和GPS坐标计算方式是不同的,因此有时候需要进行坐标转化再进行定位,坐标转化算法如下: //定义一些常量 var x_PI = 3.14159265358979324 * 3000.0 ...

  7. GPS坐标偏移与转换

    一.坐标类型: 国家规定,中国大陆所有公开地理数据都需要至少用GCJ-02进行加密 WGS84坐标,GPS设备获取的角度坐标,中国地图产品都不允许使用GPS坐标 谷歌地球,googleearth上取到 ...

  8. 坐标计算距离公式 火星坐标系_根据经纬度计算距离的公式、百度坐标转换成GPS坐标(PHP版)...

    //百度坐标转换成GPS坐标 $lnglat = '121.437518,31.224665'; function FromBaiduToGpsXY($lnglat){ // 经度,纬度 $lngla ...

  9. GPS坐标互转:WGS-84(GPS)、GCJ-02(Google地图)、BD-09(百度地图)

    GPS坐标互转:WGS-84(GPS).GCJ-02(Google地图).BD-09(百度地图) WGS-84:是国际标准,GPS坐标(Google Earth使用.或者GPS模块) GCJ-02:中 ...

最新文章

  1. 大学毕业生如何应对“没有工作经验”的难题
  2. 机会:新兴产业+大消费(转载)
  3. 上网速度快的浏览器_2020年浏览器国内市场占有率排行榜
  4. NJUST4316(立体几何投影的面积交)
  5. 计算整数中各位数字之和
  6. SPOJ LCS 后缀自动机
  7. 为SSD加速 实战4KB对齐技巧3/3
  8. python 字符串替换_python字符串替换的2种方法
  9. 在Mac OS X中使用VIM开发STM32(3)
  10. 被神话的Linux, 一文带你看清Linux在多核可扩展性设计上的不足
  11. TurboCAD Mac Pro 12(CAD设计绘图软件)
  12. python 函数注释出现 :expected an indented block
  13. worldpress小工具自定义html,教程 | 用 worldpress 搭建自己的博客 – 小白教程
  14. 浅谈微信小程序和微信公众平台
  15. 沉浸其境,共赴云栖数智硬核美学
  16. wzoi‘‘买苹果‘‘
  17. 2022年高考全国甲卷英语作文保护海洋
  18. 95后国风艺术家联手AI,立志耗时一年,作百米长卷《新西湖繁胜全景图》献礼杭州亚运!...
  19. python爬虫实战小项目
  20. 默纳克调试说明书_默纳克NICE调试说明书修改版

热门文章

  1. virtuoso从电路图导入版图_基于Virtuoso平台的单片射频收发系统PCB电路仿真与版图设计...
  2. 谷歌浏览器皮肤怎么换 修改谷歌浏览器皮肤的方法
  3. bootstrap动态调用select下拉框
  4. 【自然语言处理】 知识融入方法
  5. Matlab调用Cplex的二三事
  6. ORACLE19C中锁的速查表
  7. 当软删除遇到唯一索引
  8. 2021年《故宫日历》青少版亮相,一本可以听的百科日历
  9. 流利阅读Day11 杜克大学道歉
  10. 谁会嫌钱多啊,最适合学生党的Python兼职攻略以及接私活经验!