我在《大地经纬度坐标与地心地固坐标的的转换》这篇文章中已经论述了大地坐标系/地理坐标系的概念，简单来说就是由经度、纬度以及高程（BLH）确定的坐标系，它是一种曲面坐标。
然而，在实际使用过程中我们用的最多的还是平面坐标，并且单位最好与常用的长度单位（米）一致。所以就产生了从曲面到平面的转换，这个过程也叫做投影，转换的结果也就是投影平面坐标系。我在《GDAL坐标转换》这篇文章中详细论述了我们国内常用的三种投影平面坐标系：横轴墨卡托投影，高斯-克吕格投影和UTM投影。本质上来说，高斯-克吕格投影和UTM投影其实都是横轴墨卡托投影，横轴墨卡托投影也是用的最为广泛的地图投影方式。
但是在GIS，尤其是WebGIS领域中，横轴墨卡托投影的使用远没有Web墨卡托投影方式用的多。最重要的原因是Web墨卡托投影的转换算法比横轴墨卡托投影要简单很多，符合Web的轻量化的特点。
Web墨卡托投影是横轴墨卡托投影的特化版，要完全搞清楚Web墨卡托投影就必须得先搞清楚横轴墨卡托投影，不过横轴墨卡托投影实在太复杂了，但是我们可以定性地去理解。它的计算过程大概可以这样理解：
在X方向上，为了保证投影到平面后经线和纬线仍然垂直，那么每条纬线都会按照赤道周长展开，也就是\(2*PI*r = 2*20037508.3427892\)。由于原点位于平面中心，那么可以算得X轴的取值范围：[-20037508.3427892,20037508.3427892]。经度与投影后X长是简单的线性关系。
在Y方向上，则需要借助于墨卡托投影公式。为了保证投影的结果是正方形，那么就把Y轴的取值范围也取值成[-20037508.3427892,20037508.3427892]之间。这样做没什么道理，纯粹是为了希望投影的结果是正方形，便于切片。最后，通过墨卡托投影公式进行反算，得到的经纬度范围就是[-85.05112877980659，85.05112877980659]。也就是这种投影方式，大于这个范围是失效的。
参考Cesium的具体实现如下：
#include <iostream>
//#include <eigen3/Eigen/Eigen>
//#include <osgEarth/GeoData>
using namespace std;
const double epsilon = 0.000000000000001;
const double pi = 3.14159265358979323846;
const double d2r = pi / 180;
const double r2d = 180 / pi;
const double a = 6378137.0; //椭球长半轴
const double f_inverse = 298.257223563; //扁率倒数
const double b = a - a / f_inverse;
//const double b = 6356752.314245; //椭球短半轴
const double e = sqrt(a * a - b * b) / a;
//墨卡托范围[-PI, PI]->大地纬度范围[-PI/2, PI/2]
static double mercatorAngleToGeodeticLatitude(double mercatorAngle)
{
return pi / 2.0 - (2.0 * atan(exp(-mercatorAngle)));
//return 2.0 * atan(exp(mercatorAngle)) - pi / 2.0;
}
//Web墨卡托投影所支持的最大纬度（北和南）
static double maximumLatitude = mercatorAngleToGeodeticLatitude(pi);
//大地纬度范围[-PI/2, PI/2]->墨卡托范围[-PI, PI]
static double geodeticLatitudeToMercatorAngle(double latitude)
{
// Clamp the latitude coordinate to the valid Mercator bounds.
if (latitude > maximumLatitude)
{
latitude = maximumLatitude;
}
else if (latitude < -maximumLatitude)
{
latitude = -maximumLatitude;
}
double sinLatitude = sin(latitude);
return 0.5 * log((1.0 + sinLatitude) / (1.0 - sinLatitude));
}
void Blh2Wmc(double &x, double &y, double &z)
{
x = x * d2r * a;
y = geodeticLatitudeToMercatorAngle(y * d2r) * a;
}
void Wmc2Blh(double &x, double &y, double &z)
{
//var oneOverEarthSemimajorAxis = this._oneOverSemimajorAxis;
x = x / a * r2d;
y = mercatorAngleToGeodeticLatitude(y / a) * r2d;
}
int main()
{
double x = 113.6;
double y = 38.8;
double z = 100;
printf("%.10lf\n", maximumLatitude * r2d);
printf("原大地经纬度坐标：%.10lf\t%.10lf\t%.10lf\n", x, y, z);
Blh2Wmc(x, y, z);
printf("Web墨卡托坐标：%.10lf\t%.10lf\t%.10lf\n", x, y, z);
Wmc2Blh(x, y, z);
printf("转回大地经纬度坐标：%.10lf\t%.10lf\t%.10lf\n", x, y, z);
}
最终运行的结果：
通过GlobalMapper中的坐标转换工具对照的结果如下：
两者结果基本一致。
我正在使用Ruby/Mechanize编写一个“自动填写表格”应用程序。它几乎可以工作。我可以使用精彩CharlesWeb代理以查看服务器和我的Firefox浏览器之间的交换。现在我想使用Charles查看服务器和我的应用程序之间的交换。Charles在端口8888上代理。假设服务器位于https://my.host.com。.一件不起作用的事情是:@agent||=Mechanize.newdo|agent|agent.set_proxy("my.host.com",8888)end这会导致Net::HTTP::Persistent::Error:...lib/net/http/pe
Asitcurrentlystands,thisquestionisnotagoodfitforourQ&Aformat.Weexpectanswerstobesupportedbyfacts,references,orexpertise,butthisquestionwilllikelysolicitdebate,arguments,polling,orextendeddiscussion.Ifyoufeelthatthisquestioncanbeimprovedandpossiblyreopened,visitthehelpcenter提供指导。11年前关闭。我是一位精通HTML
对于我的大部分应用程序，我从Geocoder获取城镇/城市的纬度、经度、邮政编码等。我只是投入城市和州，作为返回，我得到了我处于一个有field的场景中。该场所需要一个地址，我正在从另一个来源获取该场所的纬度和经度。使用Geocodergem，我可以通过给地址提供经纬度来获取地址吗？ 最佳答案 在Rails控制台中运行。latitude=40.0397longitude=-76.30144geo_localization="#{latitude},#{longitude}"query=Geocoder.search(geo_loca
我目前正在为一个新网站设计版本化的API。我了解如何为路由命名空间，但我一直坚持在模型中实现版本化方法的最佳方式。下面的代码示例使用的是rails框架，但是事情的原理在大多数web框架之间应该是一致的。目前的路线看起来像这样:MyApp::Application.routes.drawdonamespace:apidonamespace:v1doresources:products,:only=>[:index,:show]endendend和Controller:classApi::V1::ProductsController很明显，我们只是在此处公开Product上可用的属性，如果
我正在尝试使用ruby来使用Sharepoint网络服务。我基本上已经放弃尝试使用NTLM进行身份验证，并暂时将Sharepoint服务器更改为使用基本身份验证。我已成功使用soap4r获得WSDL，但在尝试使用实际Web服务调用时仍然无法进行身份验证。有没有人有过让ruby和Sharepoint对话的经验？ 最佳答案 我是个新手。但经过很多时间并在更多经验编码人员的帮助下，我能够让ruby与Sharepoint2010一起工作。下面的代码需要“ntlm/mechanize”gem。我已经能够使用列表GUID和ListV
2022年10月21日星期五【数据指标】加密货币总市值：$0.95万亿BTC市值占比：38.51%恐慌贪婪指数：23极度恐慌 【今日快讯】1、【政讯】1.1.1、美联储布拉德：市场预期美联储11月会加息75个基点1.1.2、美联储哈克：将维持加息一段时间1.2、美国10年期国债收益率触及4.197%，为2008年6月以来最高1.3、法国数字转型部长：政府将专注于DeFi和Web31.4、巴西ATM机将于11月3日起支持USDT1.5、美众议院副议长将于11月初加入a16zCrypto担任政府事务主管1.6、香港数字资产托管机构FirstDigitalTrust首席执行官：香港仍是安全
WAF可以对网站进行扫描，识别API漏洞。API安全如何设置API安全_Web应用防火墙-阿里云帮助中心API安全如何划分API业务用途？登录认证手机验证码认证数据保存数据查询数据导出数据分享数据更新数据删除数据增加下线注销信息发送信息认证邮件信息发送邮箱验证码认证账号密码认证账号注册API安全支持检测哪些敏感数据？敏感数据级别敏感数据类型非敏感数据（N）不涉及。特级敏感数据（L0）与一级敏感数据（L1）或二级敏感数据（L2）相同。单次响应中一级敏感数据（L1）较多时，升级为特级敏感数据（L0）。单次响应中二级敏感数据（L2）较多时，升级为一级敏感数据（L1）或特级敏感数据（L0）。一级敏感数
我正在使用RubyonRails3，我正在尝试实现API以从Web服务检索帐户信息。也就是说，我想连接到具有Account类的Web服务并从show获取信息。在URIhttp:///accounts/1处路由的操作.此时，在网络服务中accounts_controller.rb我有的文件:classAccountsController@account.to_json}endendend现在我需要一些关于连接到网络服务的建议。在客户端应用程序中，我应该有一个HTTPGET请求，但这是我的问题:连接到发出HTTP请求的Web服务的“最佳”方法是什么？客户端应用程序中的这段代码有效:url=
谁能告诉我实现rubyonrailsweb应用程序30试用期的最佳方法，很像Basecampfrom37signals的方式是吗？目前我有一个用户登录页面，然后用户可以访问显示有关其产品/定价等的当前信息的仪表板。我希望用户能够注册并拥有完整的应用程序功能，然后在30天后过期。谢谢 最佳答案 创建用于x天试用期的Rails应用程序非常容易。您想为您的用户实现30天的试用期，然后执行以下操作:第1步:在application_controller.rb中创建这些方法，例如#application_controller.rbclas
一、扫描原因 （1）寻找到网站后台管理 （2）寻找未授权界面 （3）寻找网站更多隐藏信息 （4）通过使用目录扫描可以让我们发现这个网站存在多少个目录，多少个页面，探索出网站的整体结构。通过目录扫描我们还能扫描敏感文件，后台文件，数据库文件，和信息泄漏文件等等。二、方法1、robots.txt （1）Robots协议(RobotsExclusionProtocol)“网络爬虫排除标准”，网站通过Robots协议告诉搜索引擎哪些页面可以抓取，哪些页面不能抓取。 （2）同时也记录网站所具有基本的目录。