引言

蘭伯特投影簡介參見百科搜索：

蘭伯特投影在氣象數據的處理中，是比較常用的投影坐標系，根據不同區域、范圍進行投影。

proj4是專業的坐標轉換類庫，有各種語言版本的，C++，java，js，python版等，可以很方便的將坐標從一個坐標系轉換到另一個坐標系。

在前端使用的時候，應用場景需要轉換大量的坐標，就會發現使用proj4js存在性能問題，查看了一下proj4js的源代碼，發現類庫每次調用初始化很多不相關的類型，對象等，所以，在基礎上，進行了提取。

轉換代碼及說明

//初始化常用的變量，直接換算成弧度，提升計算性能

varEPSLN?=?(typeofNumber.EPSILON?==='undefined')???1.0e-10?:?Number.EPSILON;

varconv?=?180?/?Math.PI;

varHALF_PI?=?Math.PI?/?2;

varSPI?=?3.14159265359;

varTWO_PI?=?2?*?Math.PI;

vara?=?6378137;

varb?=?6356752.314245179;

vare?=?0.08181919084262157;

varlat1?=?0.52359877559829;

varlat2?=?1.04719755119659;

varlong0?=?1.8029251173101;

varlat0?=?0;

vark0?=?1;

varns;

varf0;

varrh;

//常用的轉換參數，直接提取引用

vartsfnz?=function(eccent,?phi,?sinphi)?{

varcon?=?eccent?*?sinphi;

varcom?=?0.5?*?eccent;

con?=?Math.pow(((1?-?con)?/?(1?+?con)),?com);

return(Math.tan(0.5?*?(HALF_PI?-?phi))?/?con);

};

varsign?=function(x)?{

returnx?

};

varmsfnz?=function(eccent,?sinphi,?cosphi)?{

varcon?=?eccent?*?sinphi;

returncosphi?/?(Math.sqrt(1?-?con?*?con));

};

varadjust_lon?=function(x)?{

return(Math.abs(x)?<=?SPI)???x?:?(x?-?(sign(x)?*?TWO_PI));

};

varphi2z?=function(eccent,?ts)?{

vareccnth?=?0.5?*?eccent;

varcon,?dphi;

varphi?=?HALF_PI?-?2?*?Math.atan(ts);

for(vari?=?0;?i?<=?15;?i++)?{

con?=?eccent?*?Math.sin(phi);

dphi?=?HALF_PI?-?2?*?Math.atan(ts?*?(Math.pow(((1?-?con)?/?(1?+?con)),?eccnth)))?-?phi;

phi?+=?dphi;

if(Math.abs(dphi)?<=?0.0000000001)?{

returnphi;

}

}

return-9999;

};

//根據proj4的坐標系描述字符串，解析其中的參數

functioninit(prjstr)?{

if(prjstr.indexOf("?")?>?-1)?{

var_prjArr?=?prjstr.split("?");

_prjArr.forEach(function(item,?index,?input)?{

if(item.indexOf("lat_0")?>?-1)?{

lat0?=?parseFloat(item.split("=")[1])?/?conv;

}

})

}

varsin1?=?Math.sin(lat1);

varcos1?=?Math.cos(lat1);

varms1?=?msfnz(e,?sin1,?cos1);

varts1?=?tsfnz(e,?lat1,?sin1);

varsin2?=?Math.sin(lat2);

varcos2?=?Math.cos(lat2);

varms2?=?msfnz(e,?sin2,?cos2);

varts2?=?tsfnz(e,?lat2,?sin2);

varts0?=?tsfnz(e,?lat0,?Math.sin(lat0));

if(Math.abs(lat1?-?lat2)?>?EPSLN)?{

ns?=?Math.log(ms1?/?ms2)?/?Math.log(ts1?/?ts2);

}else{

ns?=?sin1;

}

if(isNaN(ns))?{

ns?=?sin1;

}

f0?=?ms1?/?(ns?*?Math.pow(ts1,?ns));

rh?=?a?*?f0?*?Math.pow(ts0,?ns);

}

//經緯度坐標轉換蘭伯特坐標

functionprojCood(lon,?lat)?{

lon?=?lon?/?conv;

lat?=?lat?/?conv;

if(Math.abs(2?*?Math.abs(lat)?-?Math.PI)?<=?EPSLN)?{

lat?=?sign(lat)?*?(HALF_PI?-?2?*?EPSLN);

}

varcon?=?Math.abs(Math.abs(lat)?-?HALF_PI);

varts,?rh1;

if(con?>?EPSLN)?{

ts?=?tsfnz(e,?lat,?Math.sin(lat));

rh1?=?a?*?f0?*?Math.pow(ts,?ns);

}else{

con?=?lat?*?ns;

if(con?<=?0)?{

returnnull;

}

rh1?=?0;

}

vartheta?=?ns?*?adjust_lon(lon?-?long0);

varnlon?=?(rh1?*?Math.sin(theta));

varnlat?=?(rh?-?rh1?*?Math.cos(theta));

return[nlon,?nlat];

}

//蘭伯特坐標轉經緯度坐標

functioninverseProj(x1,?y1)?{

varrh1,?con,?ts;

varlat,?lon;

varx?=?x1?/?k0;

vary?=?(rh?-?y1?/?k0);

if(ns?>?0)?{

rh1?=?Math.sqrt(x?*?x?+?y?*?y);

con?=?1;

}else{

rh1?=?-Math.sqrt(x?*?x?+?y?*?y);

con?=?-1;

}

vartheta?=?0;

if(rh1?!==?0)?{

theta?=?Math.atan2((con?*?x),?(con?*?y));

}

if((rh1?!==?0)?||?(ns?>?0))?{

con?=?1?/?ns;

ts?=?Math.pow((rh1?/?(a?*?f0)),?con);

lat?=?phi2z(e,?ts);

if(lat?===?-9999)?{

returnnull;

}

}else{

lat?=?-HALF_PI;

}

lon?=?adjust_lon(theta?/?ns?+?long0);

return[lon?*?conv,?conv?*?lat];

}

對于在應用中的其他坐標系轉換，同樣能夠在其中進行代碼提取。

總結

以上是生活随笔為你收集整理的proj4经纬度bl转换xy_分享proj4js中经纬度和兰伯特投影的转换代码的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。