前提条件是安装了Mapserver for Windows ms4w 软件包中已经包含
1. 设置环境变量
set GDAL_DATA=\ms4w\gdaldata
path=%path%;\ms4w\tools\gdal-ogr;\ms4w\Apache\cgi-bin
设置坐标系,EPSG:90013
在Dos下执行这个命令:
echo ^<900913^> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs ^<^> >> \ms4w\proj\nad\epsg
2. 转换Arcgis 的 SHP 文件(用ogr2ogr实现)
G:\>ogr2ogr
Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update]
[-select field_list] [-where restricted_where]
[-sql <sql statement>]
[-spat xmin ymin xmax ymax] [-preserve_fid] [-fid FID]
[-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def]
[-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]
dst_datasource_name src_datasource_name
[-lco NAME=VALUE] [-nln name] [-nlt type] [layer [layer ...]]
可以用ogrinfo 查看shp 信息
F:\ms4w\apps\webgis\data2>ogrinfo -summary . rivers_region
INFO: Open of `.'
using driver `ESRI Shapefile' successful.
Layer name: rivers_region
Geometry: Polygon
Feature Count: 8370
Extent: (112.980256, 22.538197) - (114.054200, 23.923092)
Layer SRS WKT:
GEOGCS["Lat Long WGS84",
DATUM["WGS84",
SPHEROID["World_Geodetic_System_of_1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]]
序号: String (10.0)
编号: String (10.0)
名称: String (30.0)
类型: String (15.0)
LAT_HI: Real (10.6)
LONG_HI: Real (11.6)
LAT_LOW: Real (10.6)
LONG_LOW: Real (11.6)
现在是LonLat的坐标系,下面我们来转换:
命令格式:
ogr2ogr -t_srs EPSG:900913 <to_file.shp> <from_shape_file.shp>
>ogr2ogr -t_srs EPSG:900913 out_river.shp rivers_region.shp
执行后得到的结果是:
F:\ms4w\apps\webgis\data2>ogrinfo -summary . out_river
INFO: Open of `.'
using driver `ESRI Shapefile' successful.
Layer name: out_river
Geometry: Polygon
Feature Count: 8370
Extent: (12576904.567614, 2576266.091130) - (12696455.466834, 2744039.331212)
Layer SRS WKT:
PROJCS["Google Maps Global Mercator",
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Mercator_2SP"],
PARAMETER["standard_parallel_1",0],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",0],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",1],
EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0
x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs"]]
序号: String (10.0)
编号: String (10.0)
名称: String (30.0)
类型: String (15.0)
LAT_HI: Real (10.6)
LONG_HI: Real (11.6)
LAT_LOW: Real (10.6)
LONG_LOW: Real (11.6)
3.Raster 文件的转换方式。
Raster 文件通过是用GeoTiff的格式,中间包含了地理信息。
命令格式:
gdalwarp -s_srs EPSG:4326 -t_srs EPSG:900913 <from_file.tif> <to_file.tif>
先查看一下图片信息:
F:\ms4w\apps\webgis\data2>gdalinfo fg-900.TIF
Warning 1: TIFFFetchNormalTag:ASCII value for tag "Software" does not end in nul
l byte
Driver: GTiff/GeoTIFF
Files: fg-900.TIF
Size is 2058, 1756
Coordinate System is `'
Origin = (113.204516152319740,23.234743614544065)
Pixel Size = (0.000121859844067,-0.000121859844067)
Metadata:
TIFFTAG_SOFTWARE=
TIFFTAG_XRESOLUTION=300
TIFFTAG_YRESOLUTION=300
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 113.2045162, 23.2347436)
Lower Left ( 113.2045162, 23.0207577)
Upper Right ( 113.4553037, 23.2347436)
Lower Right ( 113.4553037, 23.0207577)
Center ( 113.3299099, 23.1277507)
Band 1 Block=2058x3 Type=Byte, ColorInterp=Palette
Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 128,0,0,255
2: 0,128,0,255
3: 128,128,0,255
4: 0,0,128,255
5: 128,0,128,255
6: 0,128,128,255
执行转换命令:
F:\ms4w\apps\webgis\data2>gdalwarp -s_srs EPSG:4326 -t_srs EPSG:900913 fg-900.TI
F fg-900-n.TIF
Warning 1: TIFFFetchNormalTag:ASCII value for tag "Software" does not end in nul
l byte
Copying color table from fg-900.TIF to new file.
Creating output file that is 1983P x 1840L.
Warning 1: TIFFFetchNormalTag:ASCII value for tag "Software" does not end in nul
l byte
Processing input file fg-900.TIF.
0...10...20...30...40...50...60...70...80...90...100 - done.
执行完成后,我们再检查一下:
F:\ms4w\apps\webgis\data2>gdalinfo fg-900-n.TIF
Driver: GTiff/GeoTIFF
Files: fg-900-n.TIF
Size is 1983, 1840
Coordinate System is:
PROJCS["Google Maps Global Mercator",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.2572235630016,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (12601869.093575137000000,2660431.700036020900000)
Pixel Size = (14.077025171240757,-14.077025171240757)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (12601869.094, 2660431.700) (113d12'16.26"E, 23d22'28.38"N)
Lower Left (12601869.094, 2634529.974) (113d12'16.26"E, 23d 9'34.49"N)
Upper Right (12629783.834, 2660431.700) (113d27'19.00"E, 23d22'28.38"N)
Lower Right (12629783.834, 2634529.974) (113d27'19.00"E, 23d 9'34.49"N)
Center (12615826.464, 2647480.837) (113d19'47.63"E, 23d16'1.59"N)
Band 1 Block=1983x4 Type=Byte, ColorInterp=Palette
Color Table (RGB with 256 entries)
0: 0,0,0,255
1: 128,0,0,255
2: 0,128,0,255
3: 128,128,0,255
4: 0,0,128,255
5: 128,0,128,255
END.
分享到:
相关推荐
我是为了将EPSG:4326坐标转换成EPSG:900913坐标,找到的工具,然后添加了两者之间的转换方式! 使用方式: private static Projection proj = ProjectionFactory .getNamedPROJ4CoordinateSystem("epsg:900913");...
地理2epsg描述提供来自地理坐标的有效 EPSG 代码。 目前仅支持 WGS 84和 UTM。安装 npm install geo2epsg 设置 var Geo2EPSG = require ( 'geo2espg' ) ;var geo2espg = new Geo2EPSG ( ) ;用法十进制度度 分 秒学分...
与同类软件相比,本软件支持多平台,图形数据采集、编辑更简洁,更符合常规操作习惯。 欢迎关注EZSurveyQQ群:523636554 EZSurvey_1_Alpha目前为测试版。请大家测试并反馈测试问题和建议,谢谢! 欢迎加入群:...
`EPSG:4326`与 `EPSG:3857`的坐标系的转换 地图瓦片数量与缩放层级的数量关系 缩放层级和分辨率的数量关系 分辨率与实际距离数量关系 通过经纬度计算瓦片地图的行列号 地图瓦片原点和`EPSG:3857`坐标系原点...
因为openlayer只支持EPSG:4326和EPSG:3857,如果需要用自定义坐标系的话需要转换,proj4js是坐标转换需要的js文件
该系统也称为 EPSG:3857 或 EPSG:900913。 可以在latLonToMeters和metersToLatLon函数中找到纬度/经度和球面墨卡托之间的双向转换:其他一些参考:导出格式导出叠加层会生成一个 gzip 压缩的 tar 存档,其中包含...
使用Python编写的坐标投影转换,需要输入的是源数据的EPSG代号和目标数据的EPSG代号,具体说明文档见文件夹内
EPSG:4326与 EPSG:3857的坐标系的转换 地图瓦片数量与缩放层级的数量关系 缩放层级和分辨率的数量关系 - 不懂运行,下载完可以私聊问,可远程教学 该资源内项目源码是个人的毕设,代码都测试ok,都是运行成功后...
TrasfSianItalia程序与javaproj一起使用,该程序也可以下载,它源于对SIAN上显示的坐标进行转换的需要。 该程序只需要转换epsg代码的坐标并将其显示在Google地图上即可。 转换后的坐标可以保存和打印。 它需要Java 8...
操作系统转换一组JavaScript辅助函数,用于在OSGB36 /英国国家网格(EPSG:27700)和WGS84(EPSG:4326)坐标系之间进行转换; 以及为OS Grid Reference转换提供支持。安装在HTML文档的标题部分中包含以下[removed]...
cli 安装: $ npm install -g reproject使用: $ echo '{"type":"Point","coordinates":[319180, 6399862]}' | reproject --use-epsg-io --from=EPSG:3006 --to=EPSG:4326选项: --from=crs-name是用于转换GeoJSON的...
坐标系转换需要的资源
Paavo postinumerot 和传单 ...Paavo 数据也使用 ogr2ogr 转换为 EPSG:4326 坐标系。 ogr2ogr -f GeoJSON -lco ENCODING=UTF-8 -s_srs http://spatialreference.org/ref/epsg/3067/ -t_srs EPSG:4326 posti
支持转换坐标编码和GeoJSON数据,能在节点环境以及所有现代浏览器(IE8 +)中运行,gzip后仅3kb。 更多信息可以阅读 :police_car_light:注意 在发布,展示,传播数据时,请遵守相关法律规定 (禁止)批准,在测绘...
本范例是将wgs1984坐标 转换为Xian_1980_3_DEGREE_GK_Zone_40(EPSG:2364)
空间坐标系对应EPSG编号,EPSG数据库,EPSG编码数据库、包含国家2000等坐标系SRID.CSV文件, EPSG编码最新数据库,需要手动编写坐标系转换算法必备文件,
将国家地理信息学院标准坐标系(中部起源EPSG:5186)转换为GPS坐标系EPSG:4236 服务说明 输入经度和半径将激活居住在该半径内的动物的栖息地位置标记。 通过单击所需的标记可视化物种,名称和照片 单击以
webgis坐标转换文件
坐标系为EPSG:3875 注意发布服务的时候要将坐标系设置为编码为3875的坐标系,这样在ol.View中直接设置显示的中心点和zoom级别就可以了。
当前支持以下功能: 基于几何误差的图块加载替代和精炼仅支持Box边界体积平铺转换外部磁贴瓷砖类型: 批处理3D模型(b3dm) 点云(PNTS):基本实施目前不支持以下功能: EPSG以外的任何坐标系:3857 区域和球体...